Background: Nicotiana tabacum is an important economic crop. Topping, a common agricultural practice employed with flue-cured tobacco, is designed to increase leaf nicotine contents by increasing nicotine biosynthesis in roots. Many genes are found to be differentially expressed in response to topping, particularly genes involved in nicotine biosynthesis, but comprehensive analyses of early transcriptional responses induced by topping are not yet available. To develop a detailed understanding of the mechanisms regulating nicotine biosynthesis after topping, we have sequenced the transcriptomes of Nicotiana tabacum roots at seven time points following topping. Results: Differential expression analysis revealed that 4830 genes responded to topping across all time points. Amongst these, nine gene families involved in nicotine biosynthesis and two gene families involved in nicotine transport showed significant changes during the immediate 24 h period following topping. No obvious preference to the parental species was detected in the differentially expressed genes (DEGs). Significant changes in transcript levels of nine genes involved in nicotine biosynthesis and phytohormone signal transduction were validated by qRT-PCR assays. 549 genes encoding transcription factors (TFs), found to exhibit significant changes in gene expression after topping, formed 15 clusters based on similarities of their transcript level time-course profiles. 336 DEGs involved in phytohormone signal transduction, including genes functionally related to the phytohormones jasmonic acid, abscisic acid, auxin, ethylene, and gibberellin, were identified at the earliest time point after topping. Conclusions: Our research provides the first detailed analysis of the early transcriptional responses to topping in N. tabacum, and identifies excellent candidates for further detailed studies concerning the regulation of nicotine biosynthesis in tobacco roots.

Tobacco alkaloid metabolism is regulated by various transcription factors (TFs). Here, we have characterized a non-NIC2 locus gene, Ethylene Response Factor 91 (ERF91), function in regulation of alkaloid accumulation in tobacco. NtERF91 was preferentially expressed in roots and induced by jasmonic acid. Additionally, NtERF91 was able to in vitro bind to the NtPMT2 and NtQPT2 promoters via directly targeting the GCC-box elements and transactivate NtQPT2 gene expression. Ectopic overexpression of NtERF91 not only increased the expression of most nicotine biosynthetic genes, but also altered alkaloid accumulation profile, resulting in dramatically anatabine accumulation. We conclude that NtERF91 plays an overlapped but distinct role in regulating tobacco alkaloid accumulations.

Stomatal movement could be regulated by ABA signaling through synthesis of reactive oxygen species (ROS) in guard cells. By contrast, ethylene triggers biosynthesis of antioxidant flavonol to suppress ROS accumulation and prevent ABA-induced stomatal closure, but the underlying mechanism remains largely unknown. We isolated and characterized tobacco R2R3-MYB NtMYB184, belonging to the flavonol-specific SG7 subgroup. The RNAi suppression (KD line) and knockout (myb184) of NtMYB184 caused a down-regulation of flavonol biosynthetic genes and decreased the concentration of flavonol in tobacco leaves. Yeast one-hybrid assay, transactivation assay, EMSA and ChIP-qPCR demonstrated that NtMYB184 specifically binds to the promoters of flavonol biosynthetic genes via the MYBPLANT motifs. In tobacco, NtMYB184 regulates flavonol biosynthesis in guard cells to modulate ROS homeostasis and stomatal aperture. ABA-induced ROS production was accompanied by the suppression of NtMYB184 and flavonol biosynthesis, which may accelerate ABA-induced stomatal closure. Furthermore, ethylene stimulated the NtMYB184 expression and flavonol biosynthesis to suppress ROS accumulation and curb ABA-induced stomatal closure. In myb184, however, neither the flavonol and ROS concentrations nor the stomatal aperture varied between the ABA and ABA+ethylene treatments, indicating that NtMYB184 was indispensable for the antagonism between ethylene and ABA via regulating flavonol and ROS concentrations in guard cells.

Introduction: Heated tobacco (Nicotiana tabacum L.) products are heating tobacco plug at a temperature of 350°C and produce different emissions in aerosol and sensory perceptions of tobacco leaf compared with combustible tobacco. Previous study assessed different tobacco varieties in heated tobacco for sensory quality and analyzed the links between sensory scores of the final products and certain chemical classes in tobacco leaf. However, contribution of individual metabolites to sensory quality of heated tobacco remains largely open for investigation.Methods: In present study, five tobacco varieties were evaluated as heated tobacco for sensory quality by an expert panel and the volatile and non-volatile metabolites were analyzed by non-targeted metabolomics profiling.Results: The five tobacco varieties had distinct sensory qualities and can be classified into higher and lower sensory rating classes. Principle component analysis and hierarchical cluster analysis showed that leaf volatile and non-volatile metabolome annotated were grouped and clustered by sensory ratings of heated tobacco. Orthogonal projections to latent structures discriminant analysis followed by variable importance in projection and fold-change analysis revealed 13 volatiles and 345 non-volatiles able to discriminate the tobacco varieties with higher and lower sensory ratings. Some compounds such as β-damascenone, scopoletin, chlorogenic acids, neochlorogenic acids, and flavonol glycosyl derivatives had strong contribution to the prediction of sensory quality of heated tobacco. Several lyso-phosphatidylcholine and lyso-phosphatidylethanolamine lipid species, and reducing and non-reducing sugar molecules were also positively related to sensory quality.Discussion: Taken together, these discriminating volatile and non-volatile metabolites support the role of leaf metabolites in affecting the sensory quality of heated tobacco and provide new information on the types of leaf metabolites that can be used to predict applicability of tobacco varieties for heated tobacco products.

Advances in genomics technologies are making it increasingly feasible to characterize breeding lines that carry traits of agronomic interest. Tobacco germplasm lines that carry loci designated VAM and va have been extensively investigated due to their association with potyvirus resistance (both VAM and va) and defects in leaf surface compounds originating from glandular trichomes (VAM only). Molecular studies and classical genetic analyses are consistent with the model that VAM and va represent deletion mutations in the same chromosomal region. In this study, we used RNA-seq analysis, together with emerging tobacco reference genome sequence information to characterize the genomic regions deleted in tobacco lines containing VAM and va. Tobacco genotypes TI 1406 (VAM), K326-va and K326 (wild type) were analyzed using RNA-seq to generate a list of genes differentially expressed in TI 1406 and K326-va, versus the K326 control. Candidate genes were localized onto tobacco genome scaffolds and validated as being absent in only VAM, or missing in both VAM and va, through PCR analysis. These results enabled the construction of a map that predicted the relative extent of the VAM and va mutations on the distal end of chromosome 21. The RNA-seq analyses lead to the discovery that members of the cembratrienol synthase gene family are deleted in TI 1406. Transformation of TI 1406 with a cembratrienol synthase cDNA, however, did not recover the leaf chemistry phenotype. Common to both TI 1406 and K326-va was the absence of a gene encoding a specific isoform of a eukaryotic translation initiation factor (eiF4E1.S). Transformation experiments showed that ectopic expression of eiF4E1.S is sufficient to restore potyvirus susceptibility in plants possessing either the va or VAM mutant loci. We have demonstrated the feasibility of using RNA-seq and emerging whole genome sequence resources in tobacco to characterize the VAM and va deletion mutants. These results lead to the discovery of genes underlying some of the phenotypic traits associated with these historically important loci. Additionally, initial size estimations were made for the deleted regions, and dominant markers were developed that are very close to one of the deletion junctions that defines va.