Gene co-expression network analysis reveals coordinated regulation of three characteristic secondary biosynthetic pathways in tea plant (Camellia sinensis)

Abstract

BackgroundThe leaves of tea plants (Camellia sinensis) are used to produce tea, which is one of the most popular beverages consumed worldwide. The nutritional value and health benefits of tea are mainly related to three abundant characteristic metabolites; catechins, theanine and caffeine. Weighted gene co-expression network analysis (WGCNA) is a powerful system for investigating correlations between genes, identifying modules among highly correlated genes, and relating modules to phenotypic traits based on gene expression profiling. Currently, relatively little is known about the regulatory mechanisms and correlations between these three secondary metabolic pathways at the omics level in tea.ResultsIn this study, levels of the three secondary metabolites in ten different tissues of tea plants were determined, 87,319 high-quality unigenes were assembled, and 55,607 differentially expressed genes (DEGs) were identified by pairwise comparison. The resultant co-expression network included 35 co-expression modules, of which 20 modules were significantly associated with the biosynthesis of catechins, theanine and caffeine. Furthermore, we identified several hub genes related to these three metabolic pathways, and analysed their regulatory relationships using RNA-Seq data. The results showed that these hub genes are regulated by genes involved in all three metabolic pathways, and they regulate the biosynthesis of all three metabolites. It is notable that light was identified as an important regulator for the biosynthesis of catechins.ConclusionOur integrated omics-level WGCNA analysis provides novel insights into the potential regulatory mechanisms of catechins, theanine and caffeine metabolism, and the identified hub genes provide an important reference for further research on the molecular biology of tea plants.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4999-9) contains supplementary material, which is available to authorized users.

Integration of Metabolite Profiling and Transcriptome Analysis Reveals Genes Related to Volatile Terpenoid Metabolism in Finger Citron (C. medica var. sarcodactylis)

“…We used weighted correlation network analysis (WGCNA) to uncover the correlations between terpene content and expression levels of specific genes [44]. This method has been successfully applied in mining potential TFs in grapevine and sugarcane [45,46] and elucidating the coordinated regulation of compounds in tea plants [47]. This strategy reconstructed subnetworks with high connectivities, presenting gene targets for functional characterization.…”

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

“…Another novelty of our study is that we further explored the functional relevance of target genes between tissues and functional divergence of isoforms by building an isoform-based co-expression network. Meanwhile, the reason why we selected WGCNA for coexpression analysis is its diverse functions and maturity of data processing which has been successfully applied in many gene network and functional genomic studies [38,48,49]. Our findings revealed the case that some isoforms from the same gene appear to be divergent in the isoform co-expression network and selectively targeted by miR-NAs due to their difference in sequence structure, suggesting the coordination of both structural and expression alteration during gene regulation.…”

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

“…Another novelty of our study is that we further explored the functional relevance of target genes between tissues and functional divergence of isoforms by building an isoform-based co-expression network. Meanwhile, the reason why we selected WGCNA for co-expression analysis is its diverse functions and maturity of data processing which has been successfully applied in many gene network and functional genomic studies [37,47,48]. Our findings revealed the case that some isoforms from the same gene appear to be divergent in the isoform coexpression network and selectively targeted by miRNAs due to their difference of sequence structure, suggesting the coordination of both structural and expression alteration during gene regulation.…”