小鼠耳芥MAPKKK基因家族全基因组鉴定及进化与表达

doi:10.16288/j.yczz.22-261

[an error occurred while processing this directive]

Hereditas(Beijing) ›› 2022, Vol. 44 ›› Issue (11): 1044-1055.doi: 10.16288/j.yczz.22-261

• Research Article • Previous Articles Next Articles

Genome-wide identification, phylogenetic and expression of MAPKKK gene family in Arabidopsis pumila

Qianbin Zhu^1,²(), Zhicheng Gan², Xiaocui Li¹, Yingjie Zhang¹, Heming Zhao², Xianzhong Huang²()

1. College of Life Sciences, Shihezi University, Shihezi 832003, China
2. Center for Crop Biotechnology, College of Agriculture, Anhui Science and Technology University, Fengyang 233100, China

Received:2022-08-03 Revised:2022-09-18 Online:2022-11-20 Published:2022-09-30
Contact: Huang Xianzhong E-mail:qianbinzhu@126.com;huangxz@ahstu.edu.cn
Supported by:
the National Natural Science Foundation of China(32270385);the National Natural Science Foundation of China(U1303302);the National Natural Science Foundation of China(31060149);the Talent Introduction Start-up Fund Project of Anhui Science and Technology University(NXYJ202001)

Abstract

Abstract:

Mitogen-activated protein kinase kinase kinases (MAPKKKs) are important components of the MAPK cascade and play crucial roles in development and stress responses. Arabidopsis pumila is an ephemeral Brassicaceae plant growing in Xinjiang desert regions, which possesses salt tolerance. To explore the evolution and function of the MAPKKK gene family in A. pumila, 143 ApMAPKKK genes were identified from A. pumila genome by genome-wide analysis, which were categorized into three subfamilies: ZIK (20), MEKK (36) and RAF (87). There existed 74 and 72 colinear genes between A. thaliana, A. lyrata and A. pumila, respectively, indicating that this gene family expanded obviously in A. pumila genome. Evolutionary analysis revealed that there were 64 duplicated gene pairs with Ka/Ks less than 1, and purifying selection was dominant. RNA-seq data were used to analyze the expression characteristics of ApMAPKKK genes in response to salt stress and in different tissues. The results showed that most ApMAPKKK genes were up-regulated under 250 mmol/L NaCl stress. For example, ApMAPKKK18-1/2 and ApMAPKKK17-1/2 were substantially up-regulated. Tissue expression profiles showed that ApMAPKKK mainly presented six expression patterns. Some duplicated genes were differentially expressed in response to salt stress and in different tissues. These results lay a foundation for further understanding the complex mechanism of MAPKKK gene family transduction pathway in response to abiotic stresses in A. pumila.

Key words: MAPKKK, Arabidopsis pumila, ephemeral plant, salt stress, gene expression

Qianbin Zhu, Zhicheng Gan, Xiaocui Li, Yingjie Zhang, Heming Zhao, Xianzhong Huang. Genome-wide identification, phylogenetic and expression of MAPKKK gene family in Arabidopsis pumila[J]. Hereditas(Beijing), 2022, 44(11): 1044-1055.

Figures/Tables 4

References 53

[1]	Zhu JK. Abiotic stress signaling and responses in plants. Cell, 2016, 167(2): 313-324. doi: 10.1016/j.cell.2016.08.029
[2]	Zhang MM, Zhang SQ. Mitogen-activated protein kinase cascades in plant signaling. J Integr Plant Biol, 2022, 64(2): 301-341. doi: 10.1111/jipb.13215
[3]	Ichimura K, Shinozaki K, Tena G, Sheen J, Henry Y, Champion A, Kreis M, Zhang SQ, Hirt H, Wilson C, Heberle-Bors E, Ellis BE, Morris PC, Innes RW, Ecker JR, Scheel D, Klessig DF, Machida Y, Mundy J, Ohashi Y, Walker JC. Mitogen-activated protein kinase cascades in plants: a new nomenclature. Trends Plant Sci, 2002, 7(7): 301-308. pmid: 12119167
[4]	Zhang MM, Su JB, Zhang Y, Xu J, Zhang SQ. Conveying endogenous and exogenous signals: MAPK cascades in plant growth and defense. Curr Opin Plant Biol, 2018, 45(Pt A):1-10. doi: S1369-5266(17)30213-3 pmid: 29753266
[5]	Wang W, Feng BM, Zhou JM, Tang DZ. Plant immune signaling: advancing on two frontiers. J Integr Plant Biol, 2020, 62(1): 2-24. doi: 10.1111/jipb.12898
[6]	Jonak C, Okrész L, Bögre L, Hirt H. Complexity, cross talk and integration of plant MAP kinase signalling. Curr Opin Plant Biol, 2002, 5(5): 415-424. pmid: 12183180
[7]	Rao KP, Richa T, Kumar K, Raghuram B, Sinha AK. In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice. DNA Res, 2010, 17(3): 139-153. doi: 10.1093/dnares/dsq011 pmid: 20395279
[8]	Wu P, Wang WL, Li Y, Hou XL. Divergent evolutionary patterns of the MAPK cascade genes in Brassica rapa and plant phylogenetics. Hortic Res, 2017, 4: 17079. doi: 10.1038/hortres.2017.79
[9]	Rodriguez MCS, Petersen M, Mundy J. Mitogen-activated protein kinase signaling in plants. Annu Rev Plant Biol, 2010, 61: 621-649. doi: 10.1146/annurev-arplant-042809-112252 pmid: 20441529
[10]	Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, Boller T, Ausubel FM, Sheen J. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature, 2002, 415(6875): 977-983. doi: 10.1038/415977a
[11]	Matsuoka D, Yasufuku T, Furuya T, Nanmori T. An abscisic acid inducible Arabidopsis MAPKKK, MAPKKK18 regulates leaf senescence via its kinase activity. Plant Mol Biol, 2015, 87(6): 565-575. doi: 10.1007/s11103-015-0295-0
[12]	Li YY, Cai HX, Liu P, Wang CY, Gao HY, Wu CA, Yan K, Zhang SZ, Huang JG, Zheng CC. Arabidopsis MAPKKK 18 positively regulates drought stress resistance via downstream MAPKK3. Biochem Biophys Res Commun, 2017, 484(2): 292-297. doi: 10.1016/j.bbrc.2017.01.104
[13]	Shao YM, Yu XX, Xu XW, Li Y, Yuan WX, Xu Y, Mao CZ, Zhang SQ, Xu J. The YDA-MKK4/MKK5-MPK3/MPK6 cascade functions downstream of the RGF1-RGI ligand- receptor pair in regulating mitotic activity in root apical meristem. Mol Plant, 2020, 13(11): 1608-1623. doi: 10.1016/j.molp.2020.09.004
[14]	Soyano T, Nishihama R, Morikiyo K, Ishikawa M, Machida Y. NQK1/NtMEK1 is a MAPKK that acts in the NPK1 MAPKKK-mediated MAPK cascade and is required for plant cytokinesis. Genes Dev, 2003, 17(8): 1055-1067. doi: 10.1101/gad.1071103
[15]	Liu YL, Schiff M, Dinesh-Kumar SP. Involvement of MEK1 MAPKK, NTF6 MAPK, WRKY/MYB transcription factors, COI1 and CTR1 in N-mediated resistance to tobacco mosaic virus. Plant J, 2004, 38(5): 800-809. doi: 10.1111/j.1365-313X.2004.02085.x
[16]	Li FJ, Li MY, Wang P, Cox KL, Duan LS, Dever JK, Shan LB, Li ZH, He P. Regulation of cotton (Gossypium hirsutum) drought responses by mitogen-activated protein (MAP) kinase cascade-mediated phosphorylation of GhWRKY59. New Phytol, 2017, 215(4): 1462-1475. doi: 10.1111/nph.14680
[17]	Xu R, Duan PG, Yu HY, Zhou ZK, Zhang BL, Wang RC, Li J, Zhang GZ, Zhuang SS, Lyu J, Li N, Chai TY, Tian ZX, Yao SG, Li YH. Control of grain size and weight by the OsMKKK10-OsMKK4-OsMAPK6 signaling pathway in rice. Mol Plant, 2018, 11(6): 860-873. doi: S1674-2052(18)30133-3 pmid: 29702261
[18]	Guo T, Chen K, Dong NQ, Shi CL, Ye WW, Gao JP, Shan JX, Lin HX. GRAIN SIZE AND NUMBER1 negatively regulates the OsMKKK10-OsMKK4-OsMPK6 cascade to coordinate the trade-off between grain number per panicle and grain size in rice. Plant Cell, 2018, 30(4): 871-888. doi: 10.1105/tpc.17.00959
[19]	Liu ZQ, Mei EY, Tian XJ, He ML, Tang JQ, Xu M, Liu JL, Song L, Li XF, Wang ZY, Guan QJ, Xu QJ, Bu QY. OsMKKK 70 regulates grain size and leaf angle in rice through the OsMKK4-OsMAPK6-OsWRKY53 signaling pathway. J Integr Plant Biol, 2021, 63(12): 2043-2057. doi: 10.1111/jipb.13174
[20]	Wang M, Yue H, Feng KW, Deng PC, Song WN, Nie XJ. Genome-wide identification, phylogeny and expressional profiles of mitogen activated protein kinase kinase kinase (MAPKKK) gene family in bread wheat (Triticum aestivum L.). BMC Genomics, 2016, 17(1): 668. doi: 10.1186/s12864-016-2993-7
[21]	Kong XP, Lv W, Zhang D, Jiang SS, Zhang SZ, Li DQ. Genome-wide identification and analysis of expression profiles of maize mitogen-activated protein kinase kinase kinase. PLoS One, 2013, 8(2): e57714. doi: 10.1371/journal.pone.0057714
[22]	Zhang JB, Wang XP, Wang YC, Chen YH, Luo JW, Li DD, Li XB. Genome-wide identification and functional characterization of cotton (Gossypium hirsutum) MAPKKK gene family in response to drought stress. BMC Plant Biol, 2020, 20(1): 217. doi: 10.1186/s12870-020-02431-2
[23]	Huang XZ, Yang LF, Jin YH, Lin J, Liu F. Generation, annotation, and analysis of a large-scale expressed sequence tag library from Arabidopsis pumila to explore salt-responsive genes. Front Plant Sci, 2017, 8: 955. doi: 10.3389/fpls.2017.00955
[24]	Yang LF, Jin YH, Huang W, Sun Q, Liu F, Huang XZ. Full-length transcriptome sequences of ephemeral plant Arabidopsis pumila provides insight into gene expression dynamics during continuous salt stress. BMC Genomics, 2018, 19(1): 717. doi: 10.1186/s12864-018-5106-y
[25]	Li XC, Kang KC, Huang XZ, Fan YB, Song MM, Huang YJ, Ding JJ. Genome-wide identification, phylogenetic analysis and expression profiling of the MKK gene family in Arabidopsis pumila. Hereditas(Beijing), 2020, 42(4): 403-421.
	李晓翠, 康凯程, 黄先忠, 范永斌, 宋苗苗, 黄韵杰, 丁佳佳. 小拟南芥MKK基因家族全基因组鉴定及进化和表达分析. 遗传, 2020, 42(4): 403-421.
[26]	Thompson JD, Gibson TJ, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics, 2002, Chapter 2: Unit 2.3.
[27]	Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol, 2013, 30(12): 2725-2729. doi: 10.1093/molbev/mst197 pmid: 24132122
[28]	Wang YP, Tang HB, Debarry JD, Tan X, Li JP, Wang XY, Lee TH, Jin HZ, Marler B, Guo H, Kissinger JC, Paterson AH. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res, 2012, 40(7): e49. doi: 10.1093/nar/gkr1293
[29]	Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant, 2020, 13(8): 1194-1202. doi: S1674-2052(20)30187-8 pmid: 32585190
[30]	Tang K, Dong CJ, Liu JY. Genome-wide comparative analysis of the phospholipase D gene families among allotetraploid cotton and its diploid progenitors. PLoS One, 2016, 11(5): e0156281. doi: 10.1371/journal.pone.0156281
[31]	Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao- Rico S, Librado P, Ramos-Onsins SE, Sánchez-Gracia A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol, 2017, 34(12), 3299-3302. doi: 10.1093/molbev/msx248 pmid: 29029172
[32]	Yadav CB, Bonthala VS, Muthamilarasan M, Pandey G, Khan Y, Prasad M. Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database. DNA Res, 2015, 22(1): 79-90.
[33]	Niu XQ, Luo XY, Kang KC, Huang XZ, Hu NB, Sui YH, Ai H. Genome-wide identification, comparative evolution and expression analysis of PEBP gene family from Capsicum annuum. Acta Hortic Sin, 2021, 48(5), 947-959.
	牛西强, 罗潇云, 康凯程, 黄先忠, 胡能兵, 隋益虎, 艾昊. 辣椒PEBP基因家族的全基因组鉴定、比较进化与组织表达分析. 园艺学报, 2021, 48(5), 947-959.
[34]	Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 (-Delta Delta C(T)) method. Methods, 2001, 25(4): 402- 408. doi: 10.1006/meth.2001.1262 pmid: 11846609
[35]	Altenhoff AM, Studer RA, Robinson-Rechavi M, Dessimoz C. Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs. PLoS Comput Biol, 2012, 8(5): e1002514. doi: 10.1371/journal.pcbi.1002514
[36]	Mattick JS. Introns: evolution and function. Curr Opin Genet Dev, 1994, 4(6): 823-831. pmid: 7888751
[37]	Ye JQ, Yang H, Shi HT, Wei YX, Tie WW, Ding ZH, Yan Y, Luo Y, Xia ZQ, Wang WQ, Peng M, Li KM, Zhang H, Hu W. The MAPKKK gene family in cassava: genome-wide identification and expression analysis against drought stress. Sci Rep, 2017, 7(1): 14939. doi: 10.1038/s41598-017-13988-8 pmid: 29097722
[38]	Zhang TC, Qiao Q, Novikova PY, Wang Q, Yue JP, Guan YL, Ming SP, Liu TM, De J, Liu YX, Al-Shehbaz IA, Sun H, Van Montagu M, Huang JL, Van de Peer Y, Qiong L. Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude. Proc Natl Acad Sci USA, 2019, 116(14): 7137-7146. doi: 10.1073/pnas.1817580116
[39]	Mao ZM, An ZX, Zhou GL, Yang CY, Han YL, Li XY, Zhang YF. Flora of Xinjiang (vol. 2, part 2). Urumqi: Xinjiang Science and Technology Health Press, 1995, 145-146.
	毛祖美, 安争夕, 周桂玲, 杨昌友, 韩英兰, 李学禹, 张彦福. 《新疆植物志》(第二卷第二分册). 乌鲁木齐: 新疆科技卫生出版社, 1995, 145-146.
[40]	Teige M, Scheikl E, Eulgem T, Dóczi R, Ichimura K, Shinozaki K, Dangl JL, Hirt H. The MKK2 pathway mediates cold and salt stress signaling in Arabidopsis. Mol Cell, 2004, 15(1): 141-152. pmid: 15225555
[41]	Julkowska MM, Testerink C. Tuning plant signaling and growth to survive salt. Trends Plant Sci, 2015, 20(9): 586-594. doi: 10.1016/j.tplants.2015.06.008 pmid: 26205171
[42]	Hong CY, Chao YY, Yang MY, Cheng SY, Cho SC, Kao CH. NaCl-induced expression of glutathione reductase in roots of rice (Oryza sativa L.) seedlings is mediated through hydrogen peroxide but not abscisic acid. Plant Soil, 2009, 320: 103-115. doi: 10.1007/s11104-008-9874-z
[43]	Bienert GP, Schjoerring JK, Jahn TP. Membrane transport of hydrogen peroxide. Biochim Biophys Acta, 2006, 1758(8): 994-1003. pmid: 16566894
[44]	Choi SW, Lee SB, Na YJ, Jeung SG, Kim SY. Arabidopsis MAP3K16 and other salt-inducible MAP3Ks regulate ABA response redundantly. Mol Cells, 2017, 40(3): 230-242.
[45]	Matsuoka D, Soga K, Yasufuku T, Nanmori T.Control of plant growth and development by overexpressing MAP3K17, an ABA-inducible MAP3K, in Arabidopsis. Plant Biotechnol (Tokyo), 2018, 35(2): 171-176. doi: 10.5511/plantbiotechnology.18.0412a pmid: 31819720
[46]	Virk N, Li DY, Tian LM, Huang L, Hong YB, Li XH, Zhang YF, Liu B, Zhang HJ, Song FM. Arabidopsis Raf-Like mitogen-activated protein kinase kinase kinase gene Raf43 is required for tolerance to multiple abiotic stresses. PLoS One, 2015, 10(7): e0133975. doi: 10.1371/journal.pone.0133975
[47]	Kim JM, Woo DH, Kim SH, Lee SY, Park HY, Seok HY, Chung WS, Moon YH. Arabidopsis MKKK 20 is involved in osmotic stress response via regulation of MPK6 activity. Plant Cell Rep, 2012, 31(1): 217-224. doi: 10.1007/s00299-011-1157-0
[48]	Schmidt R, Mieulet D, Hubberten HM, Obata T, Hoefgen R, Fernie AR, Fisahn J, San Segundo B, Guiderdoni E, Schippers JHM, Mueller-Roeber B. Salt-responsive ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell, 2013, 25(6): 2115-2131. doi: 10.1105/tpc.113.113068
[49]	Kondrashov FA. Gene duplication as a mechanism of genomic adaptation to a changing environment. Proc Biol Sci, 2012, 279(1749): 5048-5057.
[50]	Zhao CZ, Wang PC, Si T, Hsu CC, Wang L, Zayed O, Yu ZP, Zhu YF, Dong J, Tao WA, Zhu JK. MAP kinase cascades regulate the cold response by modulating ICE1 protein stability. Dev Cell, 2017, 43(5): 618-629.e5. doi: S1534-5807(17)30783-9 pmid: 29056551
[51]	Qiu JL, Zhou L, Yun BW, Nielsen HB, Fiil BK, Petersen K, Mackinlay J, Loake GJ, Mundy J, Morris PC.Arabidopsis mitogen-activated protein kinase kinases MKK1 and MKK2 have overlapping functions in defense signaling mediated by MEKK1, MPK4, and MKS1. Plant Physiol, 2008, 148(1): 212-222. doi: 10.1104/pp.108.120006
[52]	Kong Q, Qu N, Gao MH, Zhang ZB, Ding XJ, Yang F, Li YZ, Dong OX, Chen S, Li X, Zhang YL. The MEKK1-MKK1/MKK2-MPK4 kinase cascade negatively regulates immunity mediated by a mitogen-activated protein kinase kinase kinase in Arabidopsis. Plant Cell, 2012, 24(5): 2225-2236. doi: 10.1105/tpc.112.097253
[53]	Xu HY, Zhang C, Li ZC, Wang ZR, Jiang XX, Shi YF, Tian SN, Braun E, Mei Y, Qiu WL, Li S, Wang B, Xu J, Navarre D, Ren DT, Cheng NH, Nakata PA, Graham MA, Whitham SA, Liu JZ. The MAPK kinase kinase GmMEKK1 regulates cell death and defense responses. Plant Physiol, 2018, 178(2): 907-922. doi: 10.1104/pp.18.00903

Metrics

Comments

www.chinagene.cn
备案号：京ICP备09063187号

Genome-wide identification, phylogenetic and expression of MAPKKK gene family in Arabidopsis pumila

RichHTML

PDF (PC)

Knowledge

Abstract

Cite this article

share this article

Figures/Tables 4

References 53

Related Articles 15

Recommended Articles

Metrics

Comments

[1]	Shunze Wang, Feng Jiang, Dongli Zhu, Tie-Lin Yang, Yan Guo. Application of Hi-C technology in three-dimensional genomics research and disease pathogenesis analysis [J]. Hereditas(Beijing), 2023, 45(4): 279-294.
[2]	Siyuan Xu, Jia Shou, Qiang Wu. Additional evidence of HS5-1 enhancer eRNA PEARL for protocadherin alpha gene regulation [J]. Hereditas(Beijing), 2022, 44(8): 695-764.
[3]	Yan Zhao, Chenxin Wang, Tianming Yang, Chunshuang Li, Lihong Zhang, Dongni Du, Ruoxi Wang, Jing Wang, Min Wei, Xueqing Ba. Linking oxidative DNA lesion 8-OxoG to tumor development and progression [J]. Hereditas(Beijing), 2022, 44(6): 466-477.
[4]	Yongguang Li, Yuhuan Jin, Li Guo, Hao Ai, Ruining Li, Xianzhong Huang. Genome-wide identification and expression analysis of the PEBP genes in Arabidopsis pumila [J]. Hereditas(Beijing), 2022, 44(1): 80-91.
[5]	Cong Zhou, Qiangwei Zhou, Sheng Cheng, Guoliang Li. Research progress of CTCF in mediating 3D genome formation and regulating gene expression [J]. Hereditas(Beijing), 2021, 43(9): 816-821.
[6]	Haidong Xu, Bolin Ning, Fang Mu, Hui Li, Ning Wang. Advances of functional consequences and regulation mechanisms of alternative cleavage and polyadenylation [J]. Hereditas(Beijing), 2021, 43(1): 4-15.
[7]	Xiaocui Li, Kaicheng Kang, Xianzhong Huang, Yongbin Fan, Miaomiao Song, Yunjie Huang, Jiajia Ding. Genome-wide identification, phylogenetic analysis and expression profiling of the MKK gene family in Arabidopsis pumila [J]. Hereditas(Beijing), 2020, 42(4): 403-421.
[8]	Taotao Wang, Yong Yang, Wei Wei, Chentao Lin, Liuyin Ma. Identification and expression analyses of the NAC transcription factor family in Spartina alterniflora [J]. Hereditas(Beijing), 2020, 42(2): 194-211.
[9]	Huiyou Chen, Jianmin Zhang, Baisen Li, Yonglin Deng, Gongwei Zhang. Progress on meiotic gene expression and epigenetic regulation of male sterility in Dzo cattle [J]. Hereditas(Beijing), 2020, 42(11): 1081-1092.
[10]	Xiaomeng Gao, Zhihua Zhang. Three-dimensional structure and function of chromatin regulated by “liquid-liquid phase separation” of biological macromolecules. [J]. Hereditas(Beijing), 2020, 42(1): 45-56.
[11]	Qichao Yu,Bin Song,Xuanxuan Zou,Ling Wang,Dequan Liu,Bo Li,Kun Ma. Analysis of normal tissues adjacent to the tumour-specific expressed genes in breast cancer [J]. Hereditas(Beijing), 2019, 41(7): 625-633.
[12]	Tianpei Shi,Li Zhang. Application of whole transcriptomics in animal husbandry [J]. Hereditas(Beijing), 2019, 41(3): 193-205.
[13]	Qingqian Ding,Xiaoting Wang,Liqin Hu,Xin Qi,Linhao Ge,Weiya XU,Zhaoshi Xu,Yongbin Zhou,Guanqing Jia,Xianmin Diao,Donghong Min,Youzhi Ma,Ming Chen. MYB-like transcription factor SiMYB42 from foxtail millet (Setaria italica L.) enhances Arabidopsis tolerance to low-nitrogen stress [J]. Hereditas(Beijing), 2018, 40(4): 327-338.
[14]	Yingxia Li, Tingting Zhang, Lei Ma. Structural characteristics of natural chimeric genes and their implications for gene design [J]. Hereditas(Beijing), 2018, 40(2): 135-144.
[15]	Zongchang Xu,Yingzhen Kong. Genome-wide identification, subcellular localization and gene expression analysis of the members of CESA gene family in common tobacco (Nicotiana tabacum L.) [J]. Hereditas(Beijing), 2017, 39(6): 512-524.