Genetic regulatory landscape of gene expression during leaf and stem development in common tobacco

Mou Yin; Zhimeng Zhang; Yuexuan Hou; Yueying Wan; Qiaoling Luo; Chunqiong Wang; Xiaowei Zhang; Haowei Sun; Jie Long; Jieyun Cai; Dan Chen; Shibin Xu; Yingjiao Xiao; Yingyue Li; Qiqing Tang; Chenxi Zhang; Mengdi Zhao; Xinting Yang; Jiamei Chen; Dexin Liu; Fei He; Ke Zhang; Zhijun Tong

doi:10.1016/j.jgg.2026.03.016

Turn off MathJax

Article Contents

Article Navigation > Journal of Genetics and Genomics > 2026 > Accepted Manuscript

PDF( 0 KB)

Genetic regulatory landscape of gene expression during leaf and stem development in common tobacco

doi: 10.1016/j.jgg.2026.03.016

Funds:

This research was funded by the National Natural Science Foundation of China (32460511), Yunnan Daguan Laboratory (YNDG202302QT14) and YNTC (2023530000241027).

Received Date: 2025-12-12
Accepted Date: 2026-03-16
Rev Recd Date: 2026-03-16

Available Online: 2026-03-27

Abstract

Abstract

Common tobacco (Nicotiana tabacum L.) is an important allotetraploid crop, yet the genetic regulation of gene expression during development remains poorly understood. Here, we conduct whole-genome and transcriptome sequencing of 220 diverse tobacco accessions to dissect the genetic architecture and regulatory mechanisms underlying trait variation. We identify 3.24 million SNPs and 231,700 InDels, characterizing population structure and genomic features. Our analysis reveals two genetically distinct subpopulations within the flue-cured tobacco that are differentiated across five large genomic segments on chromosomes 2S, 2T and 12S, which exhibit elevated SNP density, nucleotide diversity and fixation indices. A genetic regulatory atlas of gene expression, constructed using expression quantitative trait locus (eQTL) mapping across four developmental stages, reveals asymmetric gene expression regulation between sub-genomes. Genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) for eight leaf and stem architecture traits identify 43 quantitative trait loci (QTLs) and 153 candidate genes. Furthermore, we jointly analyze TWAS and eQTL summary statistics to construct gene regulatory networks underlying leaf and stem traits, uncovering substantial rewiring of gene regulations across developmental stages. The genomic and transcriptomic datasets generated in this study provide a valuable resource and insights into the genetic architecture of important traits in common tobacco.
- GWAS,
- TWAS,
- eQTL,
- Development,
- Tobacco,
- Allotetraploid

FullText(HTML)

References(67)

References

Aguet, F., Anand, S., Ardlie, K.G., Gabriel, S., Getz, G.A., Graubert, A., Hadley, K., Handsaker, R.E., Huang, K.H., Kashin, S., 2020. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 369, 1318-1330.

Alexander, D.H., Novembre, J., Lange, K., 2009. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655-1664.

Bally, J., Jung, H., Mortimer, C., Naim, F., Philips, J.G., Hellens, R., Bombarely, A., Goodin, M.M., Waterhouse, P.M., 2018. The rise and rise of Nicotiana benthamiana: a plant for all reasons. Annu. Rev. Phytopathol. 56, 405-426.

Bastian, M., Heymann, S., Jacomy, M., 2009. Gephi: an open source software for exploring and manipulating networks, Proceedings of the international AAAI conference on web and social media. 361–362.

Benjamini, Y., Hochberg, Y., 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc., Ser. B. 57, 289-300.

Bolstad, B.M., 2013. preprocessCore: A collection of pre-processing functions. R package, v1. https://bioconductor.org/packages//2.7/bioc/html/preprocessCore.html.

Cai, X., Liu, D., 1989. Identification of a 1B/1R wheat-rye chromosome translocation. Theor. Appl. Genet. 77, 81-83.

Chen, Y., Chen, Y., Shi, C., Huang, Z., Zhang, Y., Li, S., Li, Y., Ye, J., Yu, C., Li, Z., et al. 2018. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. GigaScience 7, gix120.

Chen, Y., Hou, M., Liu, L., Wu, S., Shen, Y., Ishiyama, K., Kobayashi, M., McCarty, D.R., Tan, B.-C., 2014. The maize DWARF1 encodes a gibberellin 3-oxidase and is dual localized to the nucleus and cytosol. Plant Physiol. 166, 2028-2039.

Choe, S., Dilkes, B.P., Gregory, B.D., Ross, A.S., Yuan, H., Noguchi, T., Fujioka, S., Takatsuto, S., Tanaka, A., Yoshida, S., et al. 1999. The Arabidopsis dwarf1 mutant is defective in the conversion of 24-methylenecholesterol to campesterol in brassinosteroid biosynthesis. Plant Physiol. 119, 897-907.

Cingolani, P., Platts, A., Wang L.L., Coon, M., Nguyen, T., Wang, L., Land, S.J., Lu, X., Ruden, D.M., 2012. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly 6, 80-92.

Comai, L., 2005. The advantages and disadvantages of being polyploid. Nat. Rev. Genet. 6, 836-846.

Cubas, P., Lauter, N., Doebley, J., Coen, E., 1999. The TCP domain: a motif found in proteins regulating plant growth and development. Plant J. 18, 215-222.

Danecek, P., Auton, A., Abecasis, G., Albers, C.A., Banks, E., DePristo, M.A., Handsaker, R.E., Lunter, G., Marth, G.T., Sherry, S.T., et al. Genomes Project Analysis, G., 2011. The variant call format and VCFtools. Bioinformatics 27, 2156-2158.

Debernardi, J.M., Mecchia, M.A., Vercruyssen, L., Smaczniak, C., Kaufmann, K., Inze, D., Rodriguez, R.E., Palatnik, J.F., 2014. Post-transcriptional control of GRF transcription factors by microRNA miR396 and GIF co-activator affects leaf size and longevity. Plant J. 79, 413-426.

Debernardi, J.M., Tricoli, D.M., Ercoli, M.F., Hayta, S., Ronald, P., Palatnik, J.F., Dubcovsky, J., 2020. A GRF-GIF chimeric protein improves the regeneration efficiency of transgenic plants. Nat. Biotechnol. 38, 1274-1279.

De Clercq, I., Van de Velde, J., Luo, X., Liu, L., Storme, V., Van Bel, M., Pottie, R., Vaneechoutte, D., Van Breusegem, F., Vandepoele, K., 2021. Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators. Nat. Plants 7, 500-513.

Efroni, I., Blum, E., Goldshmidt, A., Eshed, Y., 2008. A protracted and dynamic maturation schedule underlies Arabidopsis leaf development. Plant Cell 20, 2293-2306.

Emms, D.M., Kelly, S., 2019. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20, 238.

Gao, F., Dubos, C., 2024. The arabidopsis bHLH transcription factor family. Trends Plant Sci. 29, 668-680.

Gusev, A., Ko, A., Shi, H., Bhatia, G., Chung, W., Penninx, B.W.J.H., Jansen, R., de Geus, E.J.C., Boomsma, D.I., Wright, F.A., et al. 2016. Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245-252.

He, F., Wang, W., Rutter, W.B., Jordan, K.W., Ren, J., Taagen, E., DeWitt, N., Sehgal, D., Sukumaran, S., Dreisigacker, S., et al. 2022. Genomic variants affecting homoeologous gene expression dosage contribute to agronomic trait variation in allopolyploid wheat. Nat. Commun. 13, 826.

Johnson, C. S., 1999. Tobacco: Production, Chemistry and Technology. Davis D.L., Nielsen. M.T., Oxford.

Kakimoto, T., 1996. CKI1, a histidine kinase homolog implicated in cytokinin signal transduction. Science 274, 982-985.

Kirik, V., Baumlein, H., 1996. A novel leaf-specific myb-related protein with a single binding repeat. Gene 183, 109-113.

Klahre, U., Noguchi, T., Fujioka, S., Takatsuto, S., Yokota, T., Nomura, T., Yoshida, S., Chua, N.H., 1998. The Arabidopsis DIMINUTO/DWARF1 gene encodes a protein involved in steroid synthesis. Plant Cell 10, 1677-1690.

Leitch, I.J., Hanson, L., Lim, K.Y., Kovarik, A., Chase, M.W., Clarkson, J.J., Leitch, A.R., 2008. The ups and downs of genome size evolution in polyploid species of Nicotiana (Solanaceae). Ann. Bot. 101, 805-814.

Li, H., Durbin, R., 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760.

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., Durbin, R., 1000 Genome Project Data Processing, S., 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078-2079.

Li, Y., Ge, X., Peng, F., Li, W., Li, J.J., 2022. Exaggerated false positives by popular differential expression methods when analyzing human population samples. Genome Biol. 23, 79.

Liebsch, D., Palatnik, J.F., 2020. MicroRNA miR396, GRF transcription factors and GIF co-regulators: a conserved plant growth regulatory module with potential for breeding and biotechnology. Curr. Opin. Plant Biol. 53, 31-42.

Liu, L., Shang-Guan, K., Zhang, B., Liu, X., Yan, M., Zhang, L., Shi, Y., Zhang, M., Qian, Q., Li, J., et al. 2013. Brittle Culm1, a COBRA-like protein, functions in cellulose assembly through binding cellulose microfibrils. PLoS Genet. 9, e1003704.

Liu, Y., Yuan, G., Si, H., Sun, Y., Jiang, Z., Liu, D., Jiang, C., Pan, X., Yang, J., Luo, Z., et al. 2022. Identification of QTLs associated with agronomic traits in tobacco via a biparental population and an eight-way magic population. Front. Plant Sci. 13, 878267.

McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., Garimella, K., Altshuler, D., Gabriel, S., Daly, M., et al. 2010. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297-1303.

Murray, M.G., Thompson, W.F., 1980. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321-4326.

Pai, A.A., Pritchard, J.K., Gilad, Y., 2015. The genetic and mechanistic basis for variation in gene regulation. PLoS Genet. 11, e1004857.

Patro, R., Duggal, G., Love, M.I., Irizarry, R.A., Kingsford, C., 2017. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417-419.

Peng, J., Richards, D.E., Hartley, N.M., Murphy, G.P., Devos, K.M., Flintham, J.E., Beales, J., Fish, L.J., Worland, A.J., Pelica, F., et al. 1999. ‘Green revolution’ genes encode mutant gibberellin response modulators. Nature 400, 256-261.

Perreau, F., Frey, A., Effroy-Cuzzi, D., Savane, P., Berger, A., Gissot, L., Marion-Poll, A., 2020. ABSCISIC ACID-DEFICIENT4 has an essential function in both cis-violaxanthin and cis-neoxanthin synthesis. Plant Physiol. 184, 1303-1316.

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M.A.R., Bender, D., Maller, J., Sklar, P., de Bakker, P.I.W., Daly, M.J., Sham, P.C., 2007. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575.

Roudier, F., Fernandez, A.G., Fujita, M., Himmelspach, R., Borner, G.H., Schindelman, G., Song, S., Baskin, T.I., Dupree, P., Wasteneys, G.O., et al. 2005. COBRA, an Arabidopsis extracellular glycosyl-phosphatidyl inositol-anchored protein, specifically controls highly anisotropic expansion through its involvement in cellulose microfibril orientation. Plant Cell 17, 1749-1763.

Sasaki, A., Ashikari, M., Ueguchi-Tanaka, M., Itoh, H., Nishimura, A., Swapan, D., Ishiyama, K., Saito, T., Kobayashi, M., Khush, G.S., et al. 2002. A mutant gibberellin-synthesis gene in rice. Nature 416, 701-702.

Sayols, S., 2023. rrvgo: a Bioconductor package for interpreting lists of Gene Ontology terms. MicroPubl. Biol. 10.17912/micropub.biology.000811.

Schomburg, F.M., Bizzell, C.M., Lee, D.J., Zeevaart, J.A., Amasino, R.M., 2003. Overexpression of a novel class of gibberellin 2-oxidases decreases gibberellin levels and creates dwarf plants. Plant Cell 15, 151-163.

Schommer, C., Debernardi, J.M., Bresso, E.G., Rodriguez, R.E., Palatnik, J.F., 2014. Repression of cell proliferation by miR319-regulated TCP4. Mol. Plant. 7, 1533-1544.

Shabalin, A.A., 2012. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353-1358.

Sierro, N., Battey, J.N., Ouadi, S., Bakaher, N., Bovet, L., Willig, A., Goepfert, S., Peitsch, M.C., Ivanov, N.V., 2014. The tobacco genome sequence and its comparison with those of tomato and potato. Nat. Commun. 5, 3833.

Sierro, N., Battey, J.N., Ouadi, S., Bovet, L., Goepfert, S., Bakaher, N., Peitsch, M.C., Ivanov, N.V., 2013. Reference genomes and transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis. Genome Biol. 14, R60.

Song, Y., Wang, Y., Liu, Y., Li, H., Ding, J., Wu, X., Li, Y., Jiao, F., Yang, L., 2025. Whole genome re-sequencing in 437 tobacco germplasms identifies plant height candidate genes. Sci. Rep. 15, 4734.

Sreeramulu, S., Mostizky, Y., Sunitha, S., Shani, E., Nahum, H., Salomon, D., Hayun, L.B., Gruetter, C., Rauh, D., Ori, N., Sessa, G., 2013. BSKs are partially redundant positive regulators of brassinosteroid signaling in Arabidopsis. Plant J. 74, 905-919.

Stegle, O., Parts, L., Piipari, M., Winn, J., Durbin, R., 2012. Using probabilistic estimation of expression residuals (PEER) to obtain increased power and interpretability of gene expression analyses. Nat. Protoc. 7, 500-507.

Takahashi, T., Gasch, A., Nishizawa, N., Chua, N.H., 1995. The DIMINUTO gene of Arabidopsis is involved in regulating cell elongation. Genes Dev. 9, 97-107.

Tang, H., Krishnakumar, V., Zeng, X., Xu, Z., Taranto, A., Lomas, J.S., Zhang, Y., Huang, Y., Wang, Y., Yim, W.C., et al. 2024. JCVI: A versatile toolkit for comparative genomics analysis. iMeta 3, e211.

Tang, W., Kim, T.W., Oses-Prieto, J.A., Sun, Y., Deng, Z., Zhu, S., Wang, R., Burlingame, A.L., Wang, Z.Y., 2008. BSKs mediate signal transduction from the receptor kinase BRI1 in Arabidopsis. Science 321, 557-560.

Ueguchi-Tanaka, M., Ashikari, M., Nakajima, M., Itoh, H., Katoh, E., Kobayashi, M., Chow, T.-y., Hsing, Y.-i.C., Kitano, H., Yamaguchi, I., et al. 2005. GIBBERELLIN INSENSITIVE DWARF1 encodes a soluble receptor for gibberellin. Nature 437, 693-698.

Van de Peer, Y., Mizrachi, E., Marchal, K., 2017. The evolutionary significance of polyploidy. Nat. Rev. Genet. 18, 411-424.

Vosa, U., Claringbould, A., Westra, H.-J., Bonder, M.J., Deelen, P., Zeng, B., Kirsten, H., Saha, A., Kreuzhuber, R., Yazar, S., et al., 2021. Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression. Nat. Genet. 53, 1300-1310.

Vukasinovic, N., Nolan, T.M., Russinova, E., 2025. Unlocking the potential of brassinosteroids: A path to precision plant engineering. Science 390, eadu9798.

Wang, J., Zhang, Q., Tung, J., Zhang, X., Liu, D., Deng, Y., Tian, Z., Chen, H., Wang, T., Yin, W., Li, B., Lai, Z., Dinesh-Kumar, S.P., Baker, B., Li, F., 2024. High-quality assembled and annotated genomes of Nicotiana tabacum and Nicotiana benthamiana reveal chromosome evolution and changes in defense arsenals. Mol. Plant 17, 423-437.

Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., Feng, T., Zhou, L., Tang, W., Zhan, L., Fu, X., Liu, S., Bo, X., Yu, G., 2021. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation 2, 100141.

Yang, J., Benyamin, B., McEvoy, B.P., Gordon, S., Henders, A.K., Nyholt, D.R., Madden, P.A., Heath, A.C., Martin, N.G., Montgomery, G.W., et al. 2010. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42, 565-569.

Yang, J., Lee, S.H., Goddard, M.E., Visscher, P.M., 2011a. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76-82.

Yang, J., Manolio, T.A., Pasquale, L.R., Boerwinkle, E., Caporaso, N., Cunningham, J.M., de Andrade, M., Feenstra, B., Feingold, E., Hayes, M.G., et al. 2011b. Genome partitioning of genetic variation for complex traits using common SNPs. Nat. Genet. 43, 519-525.

You, J., Liu, Z., Qi, Z., Ma, Y., Sun, M., Su, L., Niu, H., Peng, Y., Luo, X., Zhu, M., et al. 2023. Regulatory controls of duplicated gene expression during fiber development in allotetraploid cotton. Nat. Genet. 55, 1987-1997.

Youn, J.H., Kim, T.W., Joo, S.H., Son, S.H., Roh, J., Kim, S., Kim, T.W., Kim, S.K., 2018. Function and molecular regulation of DWARF1 as a C-24 reductase in brassinosteroid biosynthesis in Arabidopsis. J. Exp. Bot. 69, 1873-1886.

Zan, Y., Chen, S., Ren, M., Liu, G., Liu, Y., Han, Y., Dong, Y., Zhang, Y., Si, H., Liu, Z., et al. 2025. The genome and GeneBank genomics of allotetraploid Nicotiana tabacum provide insights into genome evolution and complex trait regulation. Nat. Genet. 57, 986-996.

Zhang, C., Dong, S.-S., Xu, J.-Y., He, W.-M., Yang, T.-L., 2019. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 35, 1786-1788.

Relative Articles

Supplements(0)

Cited By

Proportional views