CB-1 and K326 are closely related tobacco cultivars; however, their cold tolerance capacities are different. K326 is much more cold tolerant than CB-1. We studied the transcriptomes and metabolomes of CB-1 and K326 leaf samples treated with cold stress. Totally, we have identified 14,590 differentially expressed genes (DEGs) in CB-1 and 14,605 DEGs in K326; there was also 200 differentially expressed metabolites in CB-1 and 194 in K326. Moreover, there were many overlapping genes (around 50%) that were cold-responsive in both plant cultivars, although there were also many differences in the cold responsive genes between the two cultivars. Importantly, for most of the overlapping cold responsive genes, the extent of the changes in expression were typically much more pronounced in K326 than in CB-1, which may help explain the superior cold tolerance of K326. Similar results were found in the metabolome analysis, particularly with the analysis of primary metabolites, including amino acids, organic acids, and sugars. The large number of specific responsive genes and metabolites highlight the complex regulatory mechanisms associated with cold stress in tobacco. In addition, our work implies that the energy metabolism and hormones may function distinctly between CB-1 and K326. Differences in gene expression and metabolite levels following cold stress treatment seem likely to have contributed to the observed difference in the cold tolerance phenotype of these two tobacco cultivars.

Drought stress is one of the most severe problem limited agricultural productivity worldwide. It has been reported that plants response to drought-stress by sophisticated mechanisms at both transcriptional and post-transcriptional levels. However, the precise molecular mechanisms governing the responses of tobacco leaves to drought stress and water status are not well understood. To identify genes and miRNAs involved in drought-stress responses in tobacco, we performed both mRNA and small RNA sequencing on tobacco leaf samples from the following three treatments: untreated-control (CL), drought stress (DL), and re-watering (WL). In total, we identified 798 differentially expressed genes (DEGs) between the DL and CL (DL vs. CL) treatments and identified 571 DEGs between the WL and DL (WL vs. DL) treatments. Further analysis revealed 443 overlapping DEGs between the DL vs. CL and WL vs. DL comparisons, and, strikingly, all of these genes exhibited opposing expression trends between these two comparisons, strongly suggesting that these overlapping DEGs are somehow involved in the responses of tobacco leaves to drought stress. Functional annotation analysis showed significant up-regulation of genes annotated to be involved in responses to stimulus and stress, (e.g., late embryogenesis abundant proteins and heat-shock proteins) antioxidant defense (e.g., peroxidases and glutathione S-transferases), down regulation of genes related to the cell cycle pathway, and photosynthesis processes. We also found 69 and 56 transcription factors (TFs) among the DEGs in, respectively, the DL vs. CL and the WL vs. DL comparisons. In addition, small RNA sequencing revealed 63 known microRNAs (miRNA) from 32 families and 368 novel miRNA candidates in tobacco. We also found that five known miRNA families (miR398, miR390, miR162, miR166, and miR168) showed differential regulation under drought conditions. Analysis to identify negative correlations between the differentially expressed miRNAs (DEMs) and DEGs revealed 92 mRNA-miRNA interactions between CL and DL plants, and 32 mRNA-miRNA interactions between DL and WL plants. This study provides a global view of the transcriptional and the post-transcriptional responses of tobacco under drought stress and re-watering conditions. Our results establish an empirical foundation that should prove valuable for further investigations into the molecular mechanisms through which tobacco, and plants more generally, respond to drought stress at multiple molecular genetic levels.

Carotenoids are important pigments in plants that play crucial roles in plant growth and in plant responses to environmental stress. Lycopene β cyclase (β-LCY) functions at the branch point of the carotenoid biosynthesis pathway, catalyzing the cyclization of lycopene. Here, a β-LCY gene from Nicotiana tabacum, designated as Ntβ-LCY1, was cloned and functionally characterized. Robust expression of Ntβ-LCY1 was found in leaves, and Ntβ-LCY1 expression was obviously induced by salt, drought, and exogenous abscisic acid treatments. Strong accumulation of carotenoids and expression of carotenoid biosynthesis genes resulted from Ntβ-LCY1 overexpression. Additionally, compared to wild-type plants, transgenic plants with overexpression showed enhanced tolerance to salt and drought stress with higher abscisic acid levels and lower levels of malondialdehyde and reactive oxygen species. Conversely, transgenic RNA interference plants had a clear albino phenotype in leaves, and some plants did not survive beyond the early developmental stages. The suppression of Ntβ-LCY1 expression led to lower expression levels of genes in the carotenoid biosynthesis pathway and to reduced accumulation of carotenoids, chlorophyll, and abscisic acid. These results indicate that Ntβ-LCY1 is not only a likely cyclization enzyme involved in carotenoid accumulation but also confers salt and drought stress tolerance in Nicotiana tabacum.

Salinity is one of the most severe forms of abiotic stress and affects crop yields worldwide. Plants respond to salinity stress via a sophisticated mechanism at the physiological, transcriptional and metabolic levels. However, the molecular regulatory networks involved in salt and alkali tolerance have not yet been elucidated. We developed an RNA-seq technique to perform mRNA and small RNA (sRNA) sequencing of plants under salt (NaCl) and alkali (NaHCO₃) stress in tobacco. Overall, 8064 differentially expressed genes (DEGs) and 33 differentially expressed microRNAs (DE miRNAs) were identified in response to salt and alkali stress. A total of 1578 overlapping DEGs, which exhibit the same expression patterns and are involved in ion channel, aquaporin (AQP) and antioxidant activities, were identified. Furthermore, genes involved in several biological processes, such as “photosynthesis” and “starch and sucrose metabolism,” were specifically enriched under NaHCO₃ treatment. We also identified 15 and 22 miRNAs that were differentially expressed in response to NaCl and NaHCO₃, respectively. Analysis of inverse correlations between miRNAs and target mRNAs revealed 26 mRNA-miRNA interactions under NaCl treatment and 139 mRNA-miRNA interactions under NaHCO₃ treatment. This study provides new insights into the molecular mechanisms underlying the response of tobacco to salinity stress.

Maca (Lepidium meyenii Walp.) is a traditional Andean crop with great potential for various sanitarian and medical functions, which is attracting increased research attention. The majority of previous Maca studies were focused on biochemistry and pharmacodynamics, while the genetic basis of its unique characteristics lagged due to a lack of genome information. The authors perform gas chromatography-mass spectrometry (GC/MS) analysis in the hypocotyls of three Maca ecotypes and identify 79 compounds. Among them, 62 compounds have distinct profiles among Maca ecotypes. To reveal the underlying regulatory mechanism of the chemical composition differences, de novo transcriptome sequencing is performed and the transcription profiles of three Maca ecotypes are comparatively analyzed. Functional analysis indicates several key pathways, including “starch and sucrose metabolism,” “phenylpropanoid biosynthesis,” “phenylalanine metabolism” and “plant-pathogen interaction,” are involved in regulating the chemical compositions of Maca. Combining metabolomics and transcriptomics analysis indicates transcription factors such as MYB and WRKY and mediators such as protein kinase and bifunctional inhibitors might be critical regulators of chemical composition in Maca. The transcriptome reference genome and differentially expressed genes (DEGs) obtained in this study might serve as an initial step to illustrate the genetic differences in nutrient component, secondary metabolites content, medicinal function and stress resistance in Maca.