Comprehensive transcriptome analysis provides insights into metabolic and gene regulatory networks in trichomes of Nicotiana tabacum.

KEY MESSAGE: Comprehensive transcriptome analysis suggested that the primary metabolism is modulated to augment the supply of substrates towards secondary metabolism operating in the glandular trichomes of Nicotiana tabacum. The comparative gene expression and co-expression network analysis revealed that certain members of transcription factor genes belonging to the MYB, HD-ZIP, ERF, TCP, SRS, WRKY and DOF families may be involved in the regulation of metabolism and/other aspects in the glandular trichomes of N. tabacum The glandular trichomes of Nicotiana tabacum are highly productive in terms of secondary metabolites and therefore have been projected to be used as a prognostic platform for metabolic engineering of valuable natural products. For obvious reasons, detailed studies pertaining to the metabolic and gene regulatory networks operating in the glandular trichomes of N. tabacum are of pivotal significance to be undertaken. We have carried out next-generation sequencing of glandular trichomes of N. tabcaum and investigated differential gene expression among different tissues, including trichome-free leaves. We identified a total of 37,269 and 37,371 genes, expressing in trichome free leaf and glandular trichomes, respectively, at a cutoff of FPKM ≥ 1. The analysis revealed that different pathways involved with the primary metabolism are modulated in glandular trichomes of N. tabacum, providing a plausible explanation for the enhanced biosynthesis of secondary metabolism in the glandular trichomes. Further, comparative gene expression analysis revealed several genes, which display preferential expression in the glandular trichomes and thereby seem to be potential candidate genes for future studies in connection to the discovery of novel trichome specific promoters. The present study also led to the comprehensive identification of 1750 transcription factor genes expressing at a cutoff of FPKM ≥ 1 in the glandular trichomes of N. tabacum. The clustering and co-expression analysis suggested that transcription factor genes belonging to HD-ZIP, ERF, WRKY, MYB, TCP, SRS and DOF families may be the major players in the regulation of gene expression in the glandular trichomes of N. tabacum. To the best of our knowledge, the present work is the first effort towards detailed identification of genes, especially regulatory genes expressing in the glandular trichomes of N. tabacum. The data resource and the empirical findings from present work in all probability must, therefore, provide a reference and background context for future work aiming at deciphering molecular mechanism of regulation of secondary metabolism and gene expression in the glandular trichomes of N. tabacum.