Full-Length Transcriptome Sequencing and Modular Organization Analysis of the Naringin/Neoeriocitrin-Related Gene Expression Pattern in Drynaria roosii

• Published in: Plant and cell physiology Vol. 59; no. 7; pp. 1398 - 1414
• Main Authors: Sun, Mei-Yu, Li, Jing-Yi, Li, Dong, Huang, Feng-Jie, Wang, Di, Li, Hui, Xing, Quan, Zhu, Hui-Bin, Shi, Lei
• Format: Journal Article
• Language: English
• Published: Japan Oxford University Press 01.07.2018
• Subjects: Chromatography, High Pressure Liquid; Disaccharides - genetics; Disaccharides - metabolism; Flavanones - genetics; Flavanones - metabolism; Gene Expression Regulation, Plant; Gene Regulatory Networks; Malonyl Coenzyme A - genetics; Malonyl Coenzyme A - metabolism; Phenylalanine - genetics; Phenylalanine - metabolism; Plant Proteins - genetics; Polypodiaceae - genetics; Polypodiaceae - metabolism; Sequence Analysis, RNA - methods; Tandem Mass Spectrometry; Transcription Factors - genetics; Transcriptome; Naringin/neoeriocitrin; Fern; Modular organization; SMRT; WGCNA
• ISSN: 0032-0781; 1471-9053
• DOI: 10.1093/pcp/pcy072

Abstract Drynaria roosii (Nakaike) is a traditional Chinese medicinal fern, known as ‘GuSuiBu’. The effective components, naringin and neoeriocitrin, share a highly similar chemical structure and medicinal function. Our HPLC–tandem mass spectrometry (MS/MS) results showed that the accumulation of naringin/neoeriocitrin depended on specific tissues or ages. However, little was known about the expression patterns of naringin/neoeriocitrin-related genes involved in their regulatory pathways. Due to a lack of basic genetic information, we applied a combination of single molecule real-time (SMRT) sequencing and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of D. roosii. According to the SGS data, the differentially expressed gene (DEG)-based heat map analysis revealed that naringin/neoeriocitrin-related gene expression exhibited obvious tissue- and time-specific transcriptomic differences. Using the systems biology method of modular organization analysis, we clustered 16,472 DEGs into 17 gene modules and studied the relationships between modules and tissue/time point samples, as well as modules and naringin/neoeriocitrin contents. We found that naringin/neoeriocitrin-related DEGs distributed in nine distinct modules, and DEGs in these modules showed significantly different patterns of transcript abundance to be linked to specific tissues or ages. Moreover, weighted gene co-expression network analysis (WGCNA) results further identified that PAL, 4CL and C4H, and C3H and HCT acted as the major hub genes involved in naringin and neoeriocitrin synthesis, respectively, and exhibited high co-expression with MYB- and basic helix–leucine-helix (bHLH)-regulated genes. In this work, modular organization and co-expression networks elucidated the tissue and time specificity of the gene expression pattern, as well as hub genes associated with naringin/neoeriocitrin synthesis in D. roosii. Simultaneously, the comprehensive transcriptome data set provided important genetic information for further research on D. roosii.