Transcriptome Analysis of Gene Expression Patterns Potentially Associated with Premature Senescence in Nicotiana tabacum L

Comparative transcriptome analysis reveals genes involved in trichome development and metabolism in tobacco

BACKGROUND

The glandular trichomes of tobacco (Nicotiana tabacum) can efficiently produce secondary metabolites. They act as natural bioreactors, and their natural products function to protect plants against insect-pests and pathogens and are also components of industrial chemicals. To clarify the molecular mechanisms of tobacco glandular trichome development and secondary metabolic regulation, glandular trichomes and glandless trichomes, as well as other different developmental tissues, were used for RNA sequencing and analysis.

RESULTS

By comparing glandless and glandular trichomes with other tissues, we obtained differentially expressed genes. They were obviously enriched in KEGG pathways, such as cutin, suberine, and wax biosynthesis, flavonoid and isoflavonoid biosynthesis, terpenoid biosynthesis, and plant-pathogen interaction. In particular, the expression levels of genes related to the terpenoid, flavonoid, and wax biosynthesis pathway mainly showed down-regulation in glandless trichomes, implying that they lack the capability to synthesize certain exudate compounds. Among the differentially expressed genes, 234 transcription factors were found, including AP2-ERFs, MYBs, bHLHs, WRKYs, Homeoboxes (HD-ZIP), and C2H2-ZFs. These transcription factor and genes that highly expressed in trichomes or specially expressed in GT or GLT. Following the overexpression of R2R3-MYB transcription factor Nitab4.5_0011760g0030.1 in tobacco, an increase in the number of branched glandular trichomes was observed.

CONCLUSIONS

Our data provide comprehensive gene expression information at the transcriptional level and an understanding of the regulatory pathways involved in glandular trichome development and secondary metabolism.

Programmed cell death associated with the formation of schizo-lysigenous aerenchyma in Nelumbo nucifera root

Nelumbo nucifera (N. nucifera) is an important aquatic economic crop with high edible, medicinal, ornamental, and ecological restoration values. Aerenchyma formation in N. nucifera root is an adaptive trait to the aquatic environment in long-term evolution. In this study, light microscopy, electron microscopy, and molecular biology techniques were used to study the process of the aerenchyma development and cytological events in N. nucifera root and the dynamic changes of aerenchyma formation under the treatment of exogenous 21% oxygen, ethylene (ET), and ET synthesis i + nhibitor 1-methylcyclopropene (1-MCP). The results showed that programmed cell death (PCD) occurred during the aerenchyma formation in N. nucifera root. Plasmalemma invagination and vacuole membrane rupture appeared in the formation stage, followed by nuclear deformation, chromatin condensation and marginalization, and terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling (TUNEL) detection was positive at this time. In the expansion stage of the aerenchyma development, cytoplasmic degradation and many vesicles appeared in the cytoplasm, and organelles began to degrade. Then the plasma membrane began to degrade, and the degradation of the cell wall was the last PCD step. After 21% oxygen was continuously filled in the rhizosphere environment of N. nucifera roots, the area of aerenchyma in N. nucifera roots was smaller than that in the control group. Moreover, ET induced the earlier occurrence of aerenchyma in N. nucifera root, but also, the area of aerenchyma became larger than that of the control. On the contrary, 1-MCP inhibited the occurrence of aerenchyma to some extent. Therefore, the formation of aerenchyma in N. nucifera root resulted from PCD, and its formation mode was schizo-lysigenous. A hypoxic environment could induce aerenchyma formation in plants. ET signal was involved in aerenchyma formation in N. nucifera root and had a positive regulatory effect. This study provides relevant data on the formation mechanism of plant aerenchyma and the cytological basis for exploring the regulation mechanism of plant aerenchyma formation.

Ultra-low concentration of chlorine dioxide regulates stress-caused premature leaf senescence in tobacco by modulating auxin, ethylene, and chlorophyll biosynthesis

Exploring novel growth regulators for premature senescence regulation is important for tobacco production. In the present study, chlorine dioxide (ClO2) was explored as a novel plant growth regulator for tobacco growth, particularly its effect on leaf senescence and root development. The results showed that 0.15 μM ClO2 maintained the lushness of detached leaves and whole plants. Also, the leaves of ClO2-treated plants exhibited a chlorophyll content of 58% higher than in CK (control) plants (P < 0.05). Besides, ClO2 treatment increased the biomass of roots and aboveground parts by 54 and 16%, respectively. The ClO2-treated plants also showed enhanced activities of antioxidant enzymes and significantly reduced malondialdehyde contents (P < 0.05). Moreover, ClO2 treatment remarkably alleviated drought-caused premature senescence in the tobacco plants and partly rescued the exogenous ethylene-caused plant dwarfism. The indole-3-acetic acid content in ClO2-treated plants was higher than in non-treated plants (P < 0.05), but ethylene content was significantly lower (P < 0.05). Gene expression analysis showed that ClO2 treatment remarkably suppressed ethylene synthase genes. However, the auxin biosynthesis and transport genes were up-regulated, with NtIAA17 increasing by five folds (P < 0.05). Further, ClO2 remarkably up-regulated the expression of chlorophyll biosynthesis genes, with a >20-fold increase in NtHEMA1 and NtCHLH expressions. These results designate ClO2 as a potential regulator for improving tobacco productivity by retaining higher chlorophyll content and promoting root growth.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

Integrative Analyses of Biochemical Properties and Transcriptome Reveal the Dynamic Changes in Leaf Senescence of Tobacco (Nicotiana tabacum L.)

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.

Contrasting transcriptional responses to Fusarium virguliforme colonization in symptomatic and asymptomatic hosts

The broad host range of Fusarium virguliforme represents a unique comparative system to identify and define differentially induced responses between an asymptomatic monocot host, maize (Zea mays), and a symptomatic eudicot host, soybean (Glycine max). Using a temporal, comparative transcriptome-based approach, we observed that early gene expression profiles of root tissue from infected maize suggest that pathogen tolerance coincides with the rapid induction of senescence dampening transcriptional regulators, including ANACs (Arabidopsis thaliana NAM/ATAF/CUC protein) and Ethylene-Responsive Factors. In contrast, the expression of senescence-associated processes in soybean was coincident with the appearance of disease symptom development, suggesting pathogen-induced senescence as a key pathway driving pathogen susceptibility in soybean. Based on the analyses described herein, we posit that root senescence is a primary contributing factor underlying colonization and disease progression in symptomatic versus asymptomatic host-fungal interactions. This process also supports the lifestyle and virulence of F. virguliforme during biotrophy to necrotrophy transitions. Further support for this hypothesis lies in comprehensive co-expression and comparative transcriptome analyses, and in total, supports the emerging concept of necrotrophy-activated senescence. We propose that F. virguliforme conditions an environment within symptomatic hosts, which favors susceptibility through transcriptomic reprogramming, and as described herein, the induction of pathways associated with senescence during the necrotrophic stage of fungal development.

Transcriptome divergence between developmental senescence and premature senescence in Nicotiana tabacum L

Senescence is a degenerative process triggered by intricate and coordinated regulatory networks, and the mechanisms of age-dependent senescence and stress-induced premature senescence still remain largely elusive. Thus we selected leaf samples of developmental senescence (DS) and premature senescence (PS) to reveal the regulatory divergence. Senescent leaves were confirmed by yellowing symptom and physiological measurement. A total of 1171 and 309 genes (DEGs) were significantly expressed respectively in the whole process of DS and PS. Up-regulated DEGs in PS were mostly related to ion transport, while the down-regulated DEGs were mainly associated with oxidoreductase activity and sesquiterpenoid and triterpenoid biosynthesis. In DS, photosynthesis, precursor metabolites and energy, protein processing in endoplasmic reticulum, flavonoid biosynthesis were notable. Moreover, we found the vital pathways shared by DS and PS, of which the DEGs were analyzed further via protein-protein interaction (PPI) network analysis to explore the alteration responding to two types of senescence. In addition, plant hormone transduction pathway was mapped by related DEGs, suggesting that ABA and ethylene signaling played pivotal roles in formulating the distinction of DS and PS. Finally, we conducted a model containing oxidative stress and ABA signaling as two hub points, which highlighted the major difference and predicted the possible mechanism under DS and PS. This work gained new insight into molecular divergence of developmental senescence and premature senescence and would provide reference on potential mechanism initiating and motivating senescence for further study.

iTRAQ-based comparative proteomic analysis reveals high temperature accelerated leaf senescence of tobacco (Nicotiana tabacum L.) during flue-curing

Tobacco (Nicotiana tabacum) is extensively cultivated all over the world for its economic value. During curing and storage, senescence occurs, which is associated with physiological and biochemical changes in postharvest plant organs. However, the molecular mechanisms involved in accelerated senescence due to high temperatures in tobacco leaves during curing need further elaboration. We studied molecular mechanisms of senescence in tobacco leaves exposed to high temperature during curing (Fresh, 38 °C and 42 °C), revealed by isobaric tags for relative and absolute quantification (iTRAQ) for the proteomic profiles of cultivar Bi'na1. In total, 8903 proteins were identified, and 2034 (1150 up-regulated and 1074 down-regulated) differentially abundant proteins (DAPs) were obtained from tobacco leaf samples. These DAPs were mainly involved in posttranslational modification, protein turnover, energy production and conversion. Sugar- and energy-related metabolic biological processes and pathways might be critical regulators of tobacco leaves exposed to high temperature during senescence. High-temperature stress accelerated tobacco leaf senescence mainly by down-regulating photosynthesis-related pathways and degrading cellular constituents to maintain cell viability and nutrient recycling. Our findings provide a valuable inventory of novel proteins involved in senescence physiology and elucidate the protein regulatory network in postharvest organs exposed to high temperatures during flue-curing.

Integrating transcriptome and metabolome reveals molecular networks involved in genetic and environmental variation in tobacco

Tobacco (Nicotiana tabacum) is one of the most widely cultivated commercial non-food crops with significant social and economic impacts. Here we profiled transcriptome and metabolome from 54 tobacco samples (2-3 replicates; n = 151 in total) collected from three varieties (i.e. genetic factor), three locations (i.e. environmental factor), and six developmental stages (i.e. developmental process). We identified 3,405 differentially expressed (DE) genes (DEGs) and 371 DE metabolites, respectively. We used quantitative real-time PCR to validate 20 DEGs, and confirmed 18/20 (90%) DEGs between three locations and 16/20 (80%) with the same trend across developmental stages. We then constructed nine co-expression gene modules and four co-expression metabolite modules , and defined seven de novo regulatory networks, including nicotine- and carotenoid-related regulatory networks. A novel two-way Pearson correlation approach was further proposed to integrate co-expression gene and metabolite modules to identify joint gene-metabolite relations. Finally, we further integrated DE and network results to prioritize genes by its functional importance and identified a top-ranked novel gene, LOC107773232, as a potential regulator involved in the carotenoid metabolism pathway. Thus, the results and systems-biology approaches provide a new avenue to understand the molecular mechanisms underlying complex genetic and environmental perturbations in tobacco.

Metatranscriptomic Analysis of Multiple Environmental Stresses Identifies RAP2.4 Gene Associated with Arabidopsis Immunity to Botrytis cinerea

In this study, we aimed to identify common genetic components during stress response responsible for crosstalk among stresses, and to determine the role of differentially expressed genes in Arabidopsis-Botrytis cinerea interaction. Of 1,554 B. cinerea up-regulated genes, 24%, 1.4% and 14% were induced by biotic, abiotic and hormonal treatments, respectively. About 18%, 2.5% and 22% of B. cinerea down-regulated genes were also repressed by the same stress groups. Our transcriptomic analysis indicates that plant responses to all tested stresses can be mediated by commonly regulated genes; and protein-protein interaction network confirms the cross-interaction between proteins regulated by these genes. Upon challenges to individual or multiple stress(es), accumulation of signaling molecules (e.g. hormones) plays a major role in the activation of downstream defense responses. In silico gene analyses enabled us to assess the involvement of RAP2.4 (related to AP2.4) in plant immunity. Arabidopsis RAP2.4 was repressed by B. cinerea, and its mutants enhanced resistance to the same pathogen. To the best of our knowledge, this is the first report demonstrating the role of RAP2.4 in plant defense against B. cinerea. This research can provide a basis for breeding programs to increase tolerance and improve yield performance in crops.

Characterization and Mapping of a Novel Premature Leaf Senescence Mutant in Common Tobacco (Nicotiana tabacum L.)

As the last stage of plant development, leaf senescence has a great impact on plant’s life cycle. Genetic manipulation of leaf senescence has been used as an efficient approach in improving the yield and quality of crop plants. Here we describe an ethyl methane sulfonate (EMS) mutagenesis induced premature leaf senescence mutant yellow leaf 1 (yl1) in common tobacco (Nicotiana tabacum L.). The yl1 plants displayed early leaf yellowing. Physiological parameters and marker genes expression indicated that the yl1 phenotype was caused by premature leaf senescence. Genetic analyses indicated that the yl1 phenotype was controlled by a single recessive gene that was subsequently mapped to a specific interval of tobacco linkage group 11 using simple sequence repeat (SSR) markers. Exogenous plant hormone treatments of leaves showed that the yl1 mutant was more sensitive to ethylene and jasmonic acid than the wild type. No similar tobacco premature leaf senescence mutants have been reported. This study laid a foundation for finding the gene controlling the mutation phenotype and revealing the molecular regulation mechanism of tobacco leaf senescence in the next stage.

Sedoheptulose-1,7-Bisphosphatase is Involved in Methyl Jasmonate- and Dark-Induced Leaf Senescence in Tomato Plants

Leaf senescence represents the final stage of leaf development and is regulated by diverse internal and environmental factors. Jasmonates (JAs) have been demonstrated to induce leaf senescence in several species; however, the mechanisms of JA-induced leaf senescence remain largely unknown in tomato plants (Solanum lycopersicum). In the present study, we tested the hypothesis that sedoheptulose-1,7-bisphosphatase (SBPase), an enzyme functioning in the photosynthetic carbon fixation in the Calvin⁻Benson cycle, was involved in methyl jasmonate (MeJA)- and dark-induced leaf senescence in tomato plants. We found that MeJA and dark induced senescence in detached tomato leaves and concomitantly downregulated the expression of SlSBPASE and reduced SBPase activity. Furthermore, CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9)-mediated mutagenesis of SlSBPASE led to senescence-associated characteristics in slsbpase mutant plants, including loss of chlorophyll, repressed photosynthesis, increased membrane ion leakage, and enhanced transcript abundance of senescence-associated genes. Collectively, our data suggest that repression of SBPase by MeJA and dark treatment plays a role in JA- and dark-induced leaf senescence.