Genome-scale mRNA transcriptomic insights into the mature growth dynamic of Cigar tobacco leaves

The harvest time is a key factor for cigar leaves with high quality, which varies greatly depending on the environment. Here, we performed a genome-scale mRNA transcriptomic analysis on the cigar cultivar CX-26 (Nicotiana tabacum L.) to evaluate the relationship between gene expression and growth state. The leaves were harvested with 67 (T1), 71(T2) and 75 (T3)-day growth. A total of 80,502 genes were detected in the CX-26 leaves, of which 64,611 genes were annotated on the reference genome. Principal component analysis showed that T1 and T2 leaves had a high overlapping pattern, while T3 leaves were distinct. Indeed, T1 and T2 leaves had fewer differential expressing genes (DEGs), while T3 leaves had 26,456 DGEs from T2 leaves, supporting the distinct growth of T3 leaves. GO annotations mainly enriched the photosynthesis-related metabolic process, catalytic activity and binding biological processes. KEEG analysis identified the key pathways including photosynthesisantenna proteins, plant hormone/MAPK signaling pathway and plant-pathogen interaction. The maturity regulation and defense response-associated hormones abscisic acid and jasmonate acid were higher in the T3 leaves than that in T1 and T2 leaves confirming the KEEG analysis. Furthermore, several photosynthesis-related enzymes and a transcription factor were highlighted in the gene regulatory network, which might regulate the dynamics of carbohydrate metabolism, lipid metabolism and energy metabolism. In summary, our study provides insight into the growth state of CX-26 cigar leaves.