ABSTRACT: Leaf senescence is an important process of growth and development in plant, and it is a programmed decline controlled by a series of genes. In this study, the biochemical properties and transcriptome at five maturity stages (M1∼M5) of tobacco leaves were analyzed to reveal the dynamic changes in leaf senescence of tobacco. A total of 722, 1,534, 3,723, and 6,933 genes were differentially expressed (DEG) between M1 and M2, M1 and M3, M1 and M4, and M1 and M5, respectively. Significant changes of nitrogen, sugars, and the DEGs related to metabolite accumulation were identified, suggesting the importance of energy metabolism during leaf senescence. Gene Ontology (GO) analysis found that DEGs were enriched in biosynthetic, metabolic, photosynthesis, and redox processes, and especially, the nitrogen metabolic pathways were closely related to the whole leaf senescence process (M1∼M5). All the DEGs were grouped into 12 expression profiles according to their distinct expression patterns based on Short Time-series Expression Miner (STEM) software analysis. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis found that these DEGs were enriched in pathways of carbon metabolism, starch and sucrose metabolism, nitrogen metabolism, and photosynthesis among these expression profiles. A total of 30 core genes were examined by Weight Gene Co-expression Network Analysis (WGCNA), and they appeared to play a crucial role in the regulatory of tobacco senescence. Our results provided valuable information for further functional investigation of leaf senescence in plants.
