- Research article
- Open Access
Combining transcriptomics and metabolomics to reveal the underlying molecular mechanism of ergosterol biosynthesis during the fruiting process of Flammulina velutipes
BMC Genomics volume 20, Article number: 999 (2019)
Flammulina velutipes has been recognized as a useful basidiomycete with nutritional and medicinal values. Ergosterol, one of the main sterols of F. velutipes is an important precursor of novel anticancer and anti-HIV drugs. Therefore, many studies have focused on the biosynthesis of ergosterol and have attempted to upregulate its content in multiple organisms. Great progress has been made in understanding the regulation of ergosterol biosynthesis in Saccharomyces cerevisiae. However, this molecular mechanism in F. velutipes remains largely uncharacterized.
In this study, nine cDNA libraries, prepared from mycelia, young fruiting bodies and mature fruiting bodies of F. velutipes (three replicate sets for each stage), were sequenced using the Illumina HiSeq™ 4000 platform, resulting in at least 6.63 Gb of clean reads from each library. We studied the changes in genes and metabolites in the ergosterol biosynthesis pathway of F. velutipes during the development of fruiting bodies. A total of 13 genes (6 upregulated and 7 downregulated) were differentially expressed during the development from mycelia to young fruiting bodies (T1), while only 1 gene (1 downregulated) was differentially expressed during the development from young fruiting bodies to mature fruiting bodies (T2). A total of 7 metabolites (3 increased and 4 reduced) were found to have changed in content during T1, and 4 metabolites (4 increased) were found to be different during T2. A conjoint analysis of the genome-wide connection network revealed that the metabolites that were more likely to be regulated were primarily in the post-squalene pathway.
This study provides useful information for understanding the regulation of ergosterol biosynthesis and the regulatory relationship between metabolites and genes in the ergosterol biosynthesis pathway during the development of fruiting bodies in F. velutipes.
Edible fungi are the sixth largest crop in China with a total output of 33 million tons in 2015 . Flammulina velutipes (F. velutipes) has been recognized as a model industrial basidiomycete; it is one of the most commonly used edible fungi, serving as an excellent source of vitamins, amino acids, polysaccharides, fibre, terpenoids, phenolic acids, steroids, fatty acids and other metabolites, and is widely cultivated worldwide [2,3,4,5]. Compounds with pharmaceutical value can be isolated from the fruiting bodies or mycelia of F. velutipes, including anti-inflammatory and immunomodulatory proteins , antitumour, antioxidant and acetylcholinesterase inhibitory polysaccharides, antitumour agglutinins and immunomodulatory compounds , antimicrobial terpenoids , and antitumour and antioxidant sterols [9, 10]. The active antitumour sterols include ergosterol, 22,23-dihydroergosterol, ergosta-5,8,22-trien-3-ol and ergo-8(14)-ene-3-ol [11, 12]. The chemical composition of sterols is mainly ergosterol (54.8%) and 22,23-dihydroergosterol (27.9%) . GC-MS or HPLC studies of saponification extraction have revealed that the ergosterol content in F. velutipes was 35.5 mg/100 g in wet weight or 68.0 mg/100 g in dry weight [13, 14].
Ergosterol (C28H43OH) is a typical fugal sterol and an important constituent of various membrane structures of fungal cells, and it contributes to multiple physiological functions in cells, such as cell viability, membrane permeability, membrane fluidity, membrane integrity and intracellular transport. Therefore, when ergosterol is lacking, abnormal cell membrane function and even cell rupture may occur . In recent years, a variety of fungicides collectively known as sterol biosynthesis inhibitors (SBIs) have been successfully developed to target certain enzymes or end products of the ergosterol biosynthesis pathway and have been widely used in medicine and agricultural production . More importantly, ergosterol and some of its biosynthetic intermediates have great economic value. In the pharmaceutical industry, ergosterol is an important precursor of vitamin D2, progesterone, hydrocortisone, and brassinolide, and the products of almost all steps of its biosynthesis are drug precursors [17, 18].
Ergosterol and its derivatives are obtained mainly by chemical synthesis, genetic engineering and metabolic engineering . Because of the various steps, long route, low efficiency and high cost involved, chemical synthesis of ergosterol and its derivatives is not the preferable way to obtain these compounds. One of the main approaches for producing ergosterol and its derivatives includes metabolic engineering of yeast, but because the content of ergosterol in cells is low, this production method is not efficient . The biosynthesis of ergosterol is an extremely complicated process. Transcriptional regulation of the expression of related genes is one of the main means of adjusting ergosterol biosynthesis, and feedback regulation can play an important role in ergosterol production [21, 22]. Sterol regulatory element-binding proteins (SREBPs) are transcription factors that bind to the sterol regulatory element DNA sequence. Therefore, the manipulation of biosynthesis genes by genetic engineering may be an effective way to modulate sterol biosynthesis and intracellular sterol components. Although progress has been made in the metabolic and genetic engineering of synthetic pathways in Saccharomyces cerevisiae (S. cerevisiae), the roles of ergosterol biosynthesis genes in fruiting body growth and associated metabolic changes remain a mystery.
The ergosterol biosynthetic pathway can be divided into two parts: the mevalonate pathway and the post-squalene pathway (Fig. 1). Part 1 includes nine steps (Fig. 1a) in the synthesis of farnesyl pyrophosphate from acetyl-CoA. The first step produces acetoacetyl-CoA from two acetyl-CoA molecules whose formation was previously catalysed by acetoacetyl-CoA thiolase (ERG10) . Then, ERG13, HMG, ERG12, ERG8, ERG19, IDI1 and ERG20 successively catalyse eight reactions to produce farnesyl pyrophosphate from acetoacetyl-CoA. The enzymes in the mevalonate pathway are essential genes that are conserved in eukaryotes [24, 25]. Part 2 comprises 14 steps in the production of ergosterol from farnesyl pyrophosphate (Fig. 1b). The first step forms squalene from farnesyl pyrophosphate, and the squalene is then converted into lanosterol by squalene cyclization. Ergosterol is derived from lanosterol through steps regulated or catalysed by ERG7, ERG11, ERG24, ERG25, ERG26, ERG27, ERG6, ERG2, ERG3, ERG5 and ERG4 . As ergosterol biosynthesis is regulated by both biosynthesis regulatory genes and environmental factors, genetic engineering and the optimization of culture conditions are the two main methods for increasing ergosterol productivity. For example, oxidative-fermentative growth combined with ethanol stimulation can increase ergosterol productivity . Thus, the regulation of ergosterol biosynthesis is a complex process involving multiple factors.
Multi-omics has become a common biological approach for systematic genome analyses [28, 29]. In this study, we studied the first transcriptome and metabolome of F. velutipes samples from three developmental stages: the mycelia stage (FrI), the young fruiting bodies stage (FrII) and the mature fruiting bodies stage (FrIII). The transcriptome technique was used to identify changes in the expression of genes involved in the ergosterol biosynthesis pathway during fruiting body development. Thereafter, the metabolites in this pathway were completely scanned by nontargeted metabolomic techniques. We explored the regulatory relationship between genes and ergosterol biosynthesis during fruiting body development. The results had vital significance for understanding the metabolic pathway of ergosterol biosynthesis in F. velutipes.
The analysis of RNA-Seq data
In this study, nine libraries (FrI_1, FrI_2, FrI_3, FrII_1, FrII_2, FrII_3, FrIII_1, FrIII_2 and FrIII_3) from F. velutipes at three different developmental stages were prepared and sequenced using the Illumina HiSeq™ 4000 platform (Fig. 2). An overview of sequencing is given in Additional file 1: Table S1. After data filtering, approximately 60.29 Gb of clean reads was obtained, and at least 6.63 Gb of clean reads was generated for every library. The Q30 of each sample was approximately 92%, suggesting that the sequence data were accurate. These results demonstrated that the transcriptional profiling datasets presented satisfactory reliability for further analysis. After data filtering, the clean reads were aligned to the reference genome, and the statistical results are shown in Table 1. The ratio of mapped reads to the reference genome was approximately 82.0%.
Functional annotation and pathway enrichment of differentially expressed genes (DEGs)
A total of 4907 (2798 downregulated and 2109 upregulated) DEGs were identified during the first developmental transition (T1), and 1383 (551 downregulated and 832 upregulated) DEGs were identified during the second developmental transition (T2) (Additional file 1: Figure S1). COG assignments were used to predict and classify the possible functions of the unique sequences, and describe gene evolution processes. In this study, COG annotation functions and the COG-annotated putative proteins were classified into 24 functional groups. As shown in Fig. 3, 26.83% (2151) of the DEGs did not have COG or belonged to the category with unknown function. In total, 7.96% of the DEGs were annotated with post-translational modification, protein turnover, and chaperones; 7.51% were annotated with carbohydrate transport and metabolism; 7.27% were annotated with signal transduction mechanisms and 5.14% were annotated with secondary metabolite biosynthesis, transport and catabolism.
The GO functional annotation and classification of DEGs during T1 and T2 were assigned 45 and 42 significant shared terms, respectively, which are displayed in Additional file 1: Figure S2 and S3. The results showed that the metabolic process, cellular process, single-organism process, localization, biological regulation, cellular component organization or biogenesis and regulation of biological process terms were significantly shared GO terms in the biological process category. Membrane, cell, cell part, organelle, membrane part, macromolecular complex and organelle part were the most shared terms in the cellular component category. Catalytic activity, binding, transporter activity and structural molecule activity were markedly shared terms in the molecular function category.
KEGG pathway analysis revealed that diverse pathways were represented in the transcriptome dataset, with 5444 DEGs assigned to 121 pathways. From the bubble map of the DEG pathway enrichment analysis (only the top 20 metabolic pathways are shown) (Fig. 4), we found that two metabolic pathways related to ergosterol biosynthesis with significant enrichment were the terpenoid backbone biosynthesis pathway (ko00900) and the steroid biosynthesis pathway (ko00100).
DEGs related to ergosterol biosynthesis
We analysed the essential genes involved in the terpenoid backbone biosynthesis pathway and the steroid biosynthesis pathway in F. velutipes. The results, revealed in Table 2 show that 13 genes (6 upregulated and 7 downregulated) were differentially expressed during T1. In addition, only 1 gene (1 downregulated) in the two metabolic pathways was differentially expressed during T2. It was found that the DEGs involved in the ergosterol biosynthesis process were concentrated mainly in T1. To validate the reliability of the transcriptome data, the sequences of 12 DEGs were analysed with RT-qPCR primers. The results of the RT-qPCR analysis exhibited close similarity to the RNA-Seq results, as shown in Additional file 1: Figure S4.
Metabolic differences among the three different developmental stages of F. velutipes
RNA-Seq analysis results indicated significant differences in metabolism during the development of F. velutipes; therefore, we investigated the changes in metabolic constituents over the three developmental stages. In this study, we used 18 samples (three stages × 6 biological replicates) to observe differences in metabolic constituents among the three developmental stages of F. velutipes. The metabolome used the VIP values of the first two principal components of the multivariate PLS-DA model and a combined univariate analysis of fold change and p-value to screen for differentially expressed metabolites. The screening conditions are as follows: 1) VIP ≥ 1; 2) fold change ≥1.2 or ≤ 0.83 and 3) p-value < 0.05. These three factors were taken into account to obtain a common ion. Metabolic pathway analysis was based on the KEGG database. To compare the metabolic constituents in the three developmental stages, datasets obtained from UPLC-TOF-MS in the ESI+ (ESI−) mode were subjected to PCA. The results showed different metabolic profiles among the three groups (Fig. 5). Indeed, the first principal component (PC2) in ESI+ mode (15.45% of the total variables) and PC1 in ESI− (42.56%) were clearly separated between the FrI and FrII groups. The differences between the FrII and FrIII groups resulted from PC2 (15.45% variables) in ESI+ mode and PC2 (16.69%) in ESI− mode. A total of 1742 (2154) and 751 (944) mass ions were selected between the FrI and FrII groups and between the FrII and FrIII groups in the ESI+ (ESI−) mode, respectively (Additional file 1: Table S2).
Different accumulation of sterol derivatives at three developmental stages of F. velutipes
To understand the metabolic changes in ergosterol biosynthesis, we compared the metabolic profiles of F. velutipes at different developmental stages (Tables 3 and 4). In this study, we identified 17 metabolites involved in ergosterol biosynthesis, namely, mevalonate, mevolonate-5-phosphate, isopentenyl pyrophosphate, dimethylallyl pyrophosphate, farnesyl pyrophosphate, squalene-2-3-epoxide, lanosterol, 4,4-dimethyl-cholesta-8,14,24-trienol, 14-demethyl lanosterol, 4-methylzymosterol-carboxylate, 3-keto-4-methylzymosterol, 4-methylzymosterol, fecosterol, episterol, ergosta-5,7,24(28)-trienol, ergosta-5,7,22,24(28)-tetraenol and ergosterol, which are listed in Additional file 1: Table S3 and S4. Seven of the 17 metabolites exhibited significantly different expression levels during T1 (Table 3). Among these metabolites, the expression levels of three metabolites (isopentenyl pyrophosphate, dimethylallyl pyrophosphate and 4-methylzymosterol) were significantly increased, and the expression levels of 4 metabolites (ergosta-5,7,22,24(28)-tetraenol, 4,4-dimethyl-cholesta-8,14,24-trienol, 4-methylzymosterol-carboxylate and squalene-2-3-epoxide) were significantly decreased. The UPLC-MS profile of the change in metabolites in T2 is listed in Table 4. A total of 4 metabolites (3-keto-4-methyzymosterol, 4-methylzymosterol, episterol and ergosterol) showed significantly different concentrations. The results reveal that the expression levels of these metabolites significantly varied among the different developmental stages. To assess metabolomic performance, we measured the end product ergosterol. The results are shown in Additional file 1: Figure S5. The m/z values and retention times of the metabolomic results are consistent with the validation measurements, indicating that the metabolomic results are reliable.
Correlation analysis between transcripts and sterol derivatives reveals the regulatory network of ergosterol biosynthesis in F. velutipes
Systems biology approaches have recently emerged as highly powerful tools for discovering links between regulated genes and metabolites . To unveil the underlying regulatory mechanism in sterol derivative metabolism during the development of F. velutipes, we performed correlation analyses of the metabolites related to ergosterol biosynthesis and the transcripts at three developmental stages of F. velutipes. We compared the profiles of metabolites and gene expression at three different developmental stages of F. velutipes using Pearson’s correlation coefficient (Additional file 1: Excel S1 and S2). The regulatory network analysis that helps us to understand the correlations between the metabolites and genes is shown in Fig. 6. The results indicated that metabolites such as ergosta-5, 7, 22, 24-butenol, lanosterol, decanoate, squalene-2, 3-epoxide and 14-desmethyllenol were more likely to be regulated in the post-squalene pathway. These results could provide insight into the relationship between the genetic control of metabolite levels and metabolic impact on gene expression.
In this study, a high-quality database of the F. velutipes transcriptome was generated based on NGS technology to illustrate the gene expression reprogramming of F. velutipes at different developmental stages. RT-qPCR was used to check the reliability of the transcriptomic results for F. velutipes. The sterol profiles of F. velutipes from three development stages were generated via a UPLC-Q-TOF-MS approach. We studied the changes in the expression levels of genes and metabolic content during fruiting body development and investigated regulatory networks in the fruiting process using correlation analysis. In this study, to explore the regulation of ergosterol biosynthesis, a correlation analysis was performed on metabolites and genes at three developmental stages of F. velutipes.
In this study, through metabonomic analysis, we found that metabolite profiles were significantly different and that the contents of ergosterol biosynthesis-related metabolites significantly changed among the three developmental stages (Tables 3 and 4). These results indicated that some of the metabolites (isopentenyl pyrophosphate and dimethylallyl pyrophosphate) present in F. velutipes accumulated in young fruiting bodies, while others (3-keto-4-methylzymosterol, episterol, 4-methylzymosterol and ergosterol) accumulated in mature fruiting bodies. In early fruiting body development, the accumulation of metabolites greatly contributes to the acquisition of fruiting traits . In most cases, fruiting body development and metabolism are clearly interconnected and undergo major transitions that coincide with successive phases of fruiting body development [32, 33]. In addition, in this experiment, the culture media of the mycelia stage and fruiting bodies stage were two different media of PDA and sawdust, respectively. Park et al. found that the complexity of the respective culture media indicates a possible correlation between complexity and the number of expressed genes and metabolites (F. velutipes, PDA, MCM and sawdust) . However, there is currently no clear explanation for the exceptional expression levels in F. velutipes.
The ERG10 gene encodes an acetoacetyl-CoA thiolase that catalyses the formation of acetoacetyl-CoA from two acetyl-CoA molecules. When the levels of some sterols in the cell are low, the ERG10 gene is expressed at a higher level and then regulates the activation pathway . In this study, the expression level of ERG10 was downregulated during T1, and the results indicated that the gene may be subject to feedback regulation by sterols. In previous studies, ERG1 was identified as the key regulator of post-squalene biosynthesis in S. cerevisiae and Trichoderma harzianum [35, 36]. For example, the overexpression of ERG1 could significantly increase ergosterol biosynthesis . In S. cerevisiae, the deletion of ERG26 is lethal and disrupts the synthesis of ergosterol [38, 39]. These results indicated that ERG26 is essential for cell growth and impacts the synthesis of ergosterol. The various enzymes in the ergosterol biosynthesis pathway cooperate to tightly regulate the ergosterol content. Moreover, the genetic engineering of F. velutipes has been very successful [40, 41]. Genetic modification of the ergosterol pathway can be used for the production of sterols. Therefore, the study of ergosterol biosynthesis provides not only new ideas for enhancing ergosterol production but also findings applicable to the production of other economically interesting steroid molecules.
Effective genetic engineering approaches for efficient ergosterol production from the mycelia or fruiting bodies of a fungus cannot be devised until the metabolic pathway and regulation mechanism are well understood. Although the biosynthesis pathway of ergosterol in S. cerevisiae has been well characterized, few efforts have been made to examine ergosterol biosynthesis in F. velutipes . The results in this paper could contribute to the improvement of the production of ergosterol and its derivatives. As shown in Fig. 6, a combined analysis of the differentially produced metabolites and genes was performed with the aim of identifying regulatory relationships. This could be a useful method for comparing the correlations of metabolites or genes between different groups . A total of 13 CYPs were identified in Additional file 1: Excel S2, and the expression profiles of these genes were highly correlated with that of sterols (correlation coefficient R > 0.95). Furthermore, the majority of their expression profiles are positively correlated with sterol content profiles during development, suggesting that some of these 13 CYPs might be involved in ergosterol biosynthesis. The exact roles of these CYPs will be investigated further.
The efficiency of ergosterol biosynthesis is determined by rate-limiting enzymes, and more crucially by the optimal coordination of all enzymes . The transcription factor UPC2 has been reported to upregulate target genes involved in the biosynthesis of sterols by activating the sterol response elements in their promoter regions [43, 44]. The core motifs of the sterol-response elements have been identified in nine responsive Candida albicans ERGs (ERG1, ERG2, ERG5, ERG6, ERG10, ERG11, ERG24, ERG26 and ERG27). Our results revealed that the expression levels of most of these genes differed, which was consistent with the regulatory effect of UP2C on these genes. Most of these genes were found to be related to the post-squalene pathway, which is promising for the improvement of sterol biosynthesis in F. velutipes [44, 45]. The results are consistent with the results of the gene and sterol joint analysis network. However, the overexpression of specific enzymes could result in an imbalance in cellular sterol homeostasis in sterol intermediate accumulation, significantly reducing the total cell biomass yield and repressing end product formation. Therefore, it is necessary to avoid the overexpression of specific enzymes and to control the accumulation of some cytotoxic intermediates . In summary, a balance between the precursor supply and the catalytic activities of enzymes should be achieved for the optimal production of ergosterol.
The in vivo biosynthesis of ergosterol is a complex metabolic process involving a variety of enzymes (at least 20), reactions and genes. Due to the different positions of each reaction in the metabolic pathway, the role of metabolic regulation is different, and therefore the effects of high expression of related genes are not the same. Some genes function well in ergosterol synthesis when properly expressed, while the overexpression of these genes may inhibit ergosterol biosynthesis. These results will shed light on the molecular mechanisms responsible for ergosterol biosynthesis in F. velutipes.
This study explored the regulatory relationship between genes and ergosterol biosynthesis by jointly analysing transcriptional and metabolic changes at three developmental stages of F. velutipes. Correlation analysis highlighted regulatory genes closely associated with individual metabolites or much larger networks of genes and metabolites, thereby suggesting that a strategy based on the combined analysis of different developmental stages can be very helpful for pinpointing candidate regulatory genes linked to compositional changes and fruiting body development in F. velutipes. Researchers could utilize this dataset in genetic approaches to clarify the mechanism of ergosterol regulation. In summary, this study will be instrumental for further research on the biosynthesis of sterol-related metabolites in F. velutipes.
Strains and growth conditions
This article used F. velutipes samples from three different developmental stages, the mycelia stage (FrI), the young fruiting bodies stage (FrII) and the mature fruiting bodies stage (FrIII), as experimental materials [5, 29, 47]. T1 denotes the stage transition from mycelia to young fruiting bodies and T2 denotes the transition from young fruiting bodies to mature fruiting bodies. Mycelia, young fruiting bodies and mature fruiting bodies were acquired from Shaanxi Zhongxing Gaoke Biological Technology Co., Ltd. in Xianyang City, Shaanxi Province, China. F. velutipes mycelia were grown on potato dextrose agar (PDA) (20% potato extract, 2% glucose, 2% agar, pH 7.0) plates at 23–25 °C in darkness . For the production of fruiting bodies, the mycelia were inoculated into 570 g of poplar sawdust, 120 g of rice bran and 65% water in a 1000 ml bottle. The cultures were incubated at 25 °C in the dark for 25~40 days and then transferred to conditions that induced fruiting. The resulting F. velutipes samples at three developmental stages were collected under sterile conditions, immediately frozen and stored in liquid nitrogen at − 80 °C before the metabolic content determination and RNA isolation.
RNA extraction, library preparation and sequencing
Total RNA was extracted with the RNAprep Pure Plant Kit (Tiangen Biotech Co., Ltd., Beijing, China). The total RNA concentration, RIN value, 28S/18S and fragment size were determined using an Agilent 2100 Bioanalyser (Agilent Technologies Co. Ltd., Santa Clara, CA, USA) with the Agilent RNA 6000 Nano Kit. The purity of the samples was measured using a NanoDrop™ ultraviolet spectrophotometer. After the isolation and fragmentation of total RNA, eukaryotic mRNA was enriched by using Oligo (dT) coupled to magnetic beads. In this work, nine cDNA libraries produced from the three stages of F. velutipes (three replicates of each stage) were sequenced by using the Illumina HiSeq™ 4000 platform.
Gene prediction and functional annotation of the F. velutipes reference genome
The reference genome of F. velutipes KACC42780 was downloaded from the NCBI under accession number PRJNA191921. GeneMark-ES was used to predict genes from the genome sequences . The CDS sequences of each chromosome were extracted using in-house python scripts. The CDS sequences were translated using the Biopython package using the standard codon table . The functions of those genes were predicted using Eggnog .
Transcriptome preprocessing and differential gene expression analysis
Reads were first processed using in-house Perl scripts and filtered to remove low-quality sequences, contaminated adapters and poly-N sequences. The raw reads were mapped to the reference genome  using Hisat2 software . HT-Seq was used to calculate the counts of each gene . The Transcripts Per Kilobase Million (TPM) was calculated using the edgeR package in R programming software. Additionally, the differentially expressed genes (DEGs) were identified using the EdgeR package . A gene with log2(fold change) ≥ 1 and a false discovery rate (FDR) < 0.01 were identified as a DEG.
Functional annotation and enrichment analysis of DEGs
The protein sequences of DEGs were used to perform functional prediction and classification using Eggnog . Functional annotation by gene ontology terms (GO) was downloaded from the UniProt (http://www.uniprot.org/uniprot) database. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and clusters of orthologous groups (COG) annotation was performed using Blastall software against the KEGG database and the COG database, respectively, with a cut-off E value of 1e− 5. For functional enrichment analysis, GO terms of DEGs were compared to the genome background, and the corrected p-value less than 0.05 was set to judge the significantly enriched GO terms. The pathway enrichment analysis was performed similarly using the KEGG database, and a p-value less than 0.05 was the threshold.
To validate internal control genes for expression analysis, the ABI StepOnePlusTM RT-PCR System (Applied Biosystems, USA) was used to perform RT-qPCR and analyse the DEGs related to ergosterol biosynthesis. Amplification primers were designed using GenScript Real-time PCR (TaqMan) Primer Design, and the primer sequences are provided in Additional file 1: Table S5. cDNA and RT-qPCR were acquired using the PrimeScript™ RT Reagent Kit and SYBR Green™ Premix Ex Taq™ П (Takara Biotech Co., Ltd., Japan). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin were used as internal controls . Each reaction was performed in a total reaction mixture volume of 20 μl containing 2 μl of first-strand cDNA as the template. The amplification program was as follows: 3 min at 95 °C followed by 40 cycles of 10 s at 95 °C and 30 s at 60 °C. Three independent technical replicates and three biological replicates for each sample were run to measure and assess the performance of RT-qPCR. The expression levels of DEGs were determined using the 2−△△Ct method.
Metabolite profiling using ultraperformance LC quadrupole time-of-flight tandem MS (UPLC-Q-TOF-MS)
In this study, 18 samples produced from F. velutipes at three developmental stages (six biological replicates of each stage) were tested by using UPLC-MS/MS. All samples stored at − 80 °C were placed in a freezer at − 20 °C for 30 min and then thawed at 4 °C. A 25 mg tissue sample was weighed, and placed in an EP tube, 800 μl of a cooled solution of methanol/water (1:1) was added to each EP tube, and two small steel balls were frozen and placed in each EP tube. The sample was then placed in a tissue lyser and the parameters were set to 35 Hz for 4 min. After grinding, the ball was removed and the tube was placed in a − 20 °C refrigerator for 2 h. The sample was centrifuged at 30,000×g for 15 min at 4 °C. The EP tube was carefully removed from the centrifuge, and 550 μl of each sample was transferred to a new EP tube. The sample was placed in the rack of a centrifuge, and an image was obtained according to the order of the samples on the task sheet. The samples were subsequently analysed by LC-MS, and the remaining original samples were provided to the sample manager for storage.
All samples were produced and processed following protocols for LC-MS data analysis. First, chromatographic isolation was performed by using a UPLC system (Waters, UK). Default settings were adopted for the column type, column temperature, mobile phase flow rate, mobile phase solvent ratio and other program settings . The injection volume was 10 μl. The Q-TOF was run with the capillary tube and sampling voltages set at 3 kV and 40 V, respectively, in positive-ion mode and at 1 kV and 40 V in negative-ion mode. MS data were generated by Xevo G2 XS QTOF with a TOF mass ranging from 50 to 1200 Da and were scanned at 0.2 s. The MS/MS analysis was conducted to select, separate and detect precursor ions using 20 and 40 eV in different steps with a scan time of 0.2 s. When processing the data, the LE signal was received at 3 s intervals to measure the mass accuracy. A control sample was picked from every 10 samples for evaluation of the data acquisition performance of the LC-MS.
Integrative analysis of metabolome and transcriptome data
To unveil the regulatory mechanism of ergosterol biosynthesis, we measured the associations between metabolome and transcriptome data by calculating the cosine coefficients between fold changes of the metabolite-gene pairs. For this purpose, the mean of all biological replicates from each stage in the metabolome data and the mean value of expression of each transcript in the transcriptome data were calculated. The fold changes in developmental transitions were calculated for both the metabolome and transcriptome data. Finally, the correlation matrix was imported into Cytoscape software to visualize the gene-metabolite network .
We used F statistics and the James-Stein estimator to calculate the error variance and selected genes (fold change ≥2 and FDR < 0.01) with significant differential expression. Principal component analysis (PCA) of data from F. velutipes (three stages×n biological replicates) was performed to observe differences in metabolic composition and expressed genes from the three developmental stages. The gene expression data and metabolite data were standardized as Z scores. Pearson’s correlation coefficient was calculated both among the gene expression data and between the gene expression and metabolite accumulation data. For correlations with p < 0.01, correlation matrices were visualized using heat maps generated with MultiExperiment Viewer software version 4.0 . The p value was calculated using an independent two-sample t-test.
Availability of data and materials
The datasets supporting the conclusions of this article are included in the article and its additional file. Raw reads data has been deposited in NCBI under accession number PRJNA592256.
Clusters of orthologous groups
Differentially expressed genes
False discovery rate
- Flammulina velutipes :
Young fruiting bodies stage
Mature fruiting bodies stage
F. velutipes sterols
Kyoto Encyclopedia of Genes and Genomes
Next generation sequencing
- Q30 percentage:
Percentage of bases with sequencing error rate lower than 1‰
Real-time quantitative polymerase chain reaction
Sterol biosynthesis inhibitors
Ultra-performance liquid chromatography
Li HP, Yang WJ, Qu SX, Pei F, Luo X, Mariga AM, Ma L. Variation of volatile terpenes in the edible fungi mycelia Flammulina velutipes and communications in fungus-mite interactions. Food Res Int. 2018;103:150–5.
Cai H, Liu X, Chen Z, Liao S, Zou Y. Isolation, purification and identification of nine chemical compounds from Flammulina velutipes fruiting bodies. Food Chem. 2013;141(3):2873–9.
Liu JY, Chang MC, Meng JL, Feng CP, Zhao H, Zhang ML. Comparative proteome reveals metabolic changes during the fruiting process in Flammulina velutipes. J Agric Food Chem. 2017;65(24):5091–100.
Smiderle FR, Carbonero ER, Sassaki GL, Gorin PAJ, Iacomini M. Characterization of a heterogalactan: some nutritional values of the edible mushroom Flammulina velutipes. Food Chem. 2008;108(1):329–33.
Park YJ, Baek JH, Lee S, Kim C, Rhee H, Kim H, et al. Whole genome and global gene expression analyses of the model mushroom Flammulina velutipes reveal a high capacity for lignocellulose degradation. PLoS One. 2014;9(4):e93560.
Chang YC, Hsiao YM, Wu MF, Ou CC, Lin YW, Lue KH, Ko JL. Interruption of lung cancer cell migration and proliferation by fungal immunomodulatory protein FIP-fve from Flammulina velutipes. J Agric Food Chem. 2013;61(49):12044–52.
El Enshasy HA, Hatti-Kaul R. Mushroom immunomodulators: unique molecules with unlimited applications. Trends Biotechnol. 2013;31(12):668–77.
Xiao H, Zhong JJ. Production of useful terpenoids by higher-fungus cell factory and synthetic biology approaches. Trends Biotechnol. 2016;34(3):242–55.
Shao S, Hernandez M, Kramer JK, Rinker DL, Tsao R. Ergosterol profiles, fatty acid composition, and antioxidant activities of button mushrooms as affected by tissue part and developmental stage. J Agric Food Chem. 2010;58(22):11616–25.
Yi C, Zhong H, Tong S, Cao X, Firempong CK, Liu H, et al. Enhanced oral bioavailability of a sterol-loaded microemulsion formulation of Flammulina velutipes, a potential antitumor drug. Int J Nanomedicine. 2012;7:5067–78.
Tong S, Zhong H, Yi C, Cao X, Firempong CK, Zheng Q, Feng Y, Yu J, Xu X. Simultaneous HPLC determination of ergosterol and 22,23-dihydroergosterol in Flammulina velutipes sterol-loaded microemulsion. Biomed Chromatogr. 2014;28(2):247–54.
Yi C, Sun C, Tong S, Cao X, Feng Y, Firempong CK, Jiang X, Xu X, Yu J. Cytotoxic effect of novel Flammulina velutipes sterols and its oral bioavailability via mixed micellar nanoformulation. Int J Pharm. 2013;448(1):44–50.
Phillips KM, Ruggio DM, Horst RL, Minor B, Simon RR, Feeney MJ, Byrdwell WC, Haytowitz DB. Vitamin D and sterol composition of 10 types of mushrooms from retail suppliers in the United States. J Agric Food Chem. 2011;59(14):7841–53.
Jasinghe VJ, Perera CO. Distribution of ergosterol in different tissues of mushrooms and its effect on the conversion of ergosterol to vitamin D2 by UV irradiation. Food Chem. 2005;92(3):541–6.
Ma L, Chen H, Dong P, Lu X. Anti-inflammatory and anticancer activities of extracts and compounds from the mushroom Inonotus obliquus. Food Chem. 2013;139(1–4):503–8.
Palmie-Peixoto IV, Rocha MR, Urbina JA, de Souza W, Einicker-Lamas M, Motta MC. Effects of sterol biosynthesis inhibitors on endosymbiont-bearing trypanosomatids. FEMS Microbiol Lett. 2006;255(1):33–42.
Rhee YH, Jeong SJ, Lee HJ, Lee HJ, Koh W, Jung JH, Kim SH, Sung-Hoon K. Inhibition of STAT3 signaling and induction of SHP1 mediate antiangiogenic and antitumor activities of ergosterol peroxide in U266 multiple myeloma cells. BMC Cancer. 2012;12:28.
Pluchino LA, Liu AK, Wang HC. Reactive oxygen species-mediated breast cell carcinogenesis enhanced by multiple carcinogens and intervened by dietary ergosterol and mimosine. Free Radic Biol Med. 2015;80:12–26.
Ma BX, Ke X, Tang XL, Zheng RC, Zheng YG. Rate-limiting steps in the Saccharomyces cerevisiae ergosterol pathway: towards improved ergosta-5,7-dien-3beta-ol accumulation by metabolic engineering. World J Microbiol Biotechnol. 2018;34(4):55.
Souza CM, Schwabe TM, Pichler H, Ploier B, Leitner E, Guan XL, Wenk MR, Riezman I, Riezman H. A stable yeast strain efficiently producing cholesterol instead of ergosterol is functional for tryptophan uptake, but not weak organic acid resistance. Metab Eng. 2011;13(5):555–69.
Zhang K, Tong M, Gao K, Di Y, Wang P, Zhang C, Wu X, Zheng D. Genomic reconstruction to improve bioethanol and ergosterol production of industrial yeast Saccharomyces cerevisiae. J Ind Microbiol Biotechnol. 2015;42(2):207–18.
Yuan J, Ching CB. Dynamic control of ERG9 expression for improved amorpha-4,11-diene production in Saccharomyces cerevisiae. Microb Cell Factories. 2015;14:38.
Hu Z, He B, Ma L, Sun Y, Niu Y, Zeng B. Recent advances in Ergosterol biosynthesis and regulation mechanisms in Saccharomyces cerevisiae. Indian J Microbiol. 2017;57(3):270–7.
Layer JV, Barnes BM, Yamasaki Y, Barbuch R, Li L, Taramino S, Balliano G, Bard M. Characterization of a mutation that results in independence of oxidosqualene cyclase (Erg7) activity from the downstream 3-ketoreductase (Erg27) in the yeast ergosterol biosynthetic pathway. Biochim Biophys Acta. 2013;1831(2):361–9.
Long N, Xu X, Zeng Q, Sang H, Lu L. Erg4A and Erg4B Are Required for Conidiation and Azole Resistance via Regulation of Ergosterol Biosynthesis in Aspergillus fumigatus. Appl Environ Microbiol. 2017;83(4):e02924-16. https://doi.org/10.1128/AEM.02924-16.
Wriessnegger T, Pichler H. Yeast metabolic engineering--targeting sterol metabolism and terpenoid formation. Prog Lipid Res. 2013;52(3):277–93.
Nahlik J, Hrncirik P, Mares J, Rychtera M, Kent CA. Towards the design of an optimal strategy for the production of ergosterol from Saccharomyces cerevisiae yeasts. Biotechnol Prog. 2017;33(3):838–48.
Szymanski J, Brotman Y, Willmitzer L, Cuadros-Inostroza A. Linking gene expression and membrane lipid composition of Arabidopsis. Plant Cell. 2014;26(3):915–28.
Chen S, Xu J, Liu C, Zhu Y, Nelson DR, Zhou S, et al. Genome sequence of the model medicinal mushroom Ganoderma lucidum. Nat Commun. 2012;3:913.
Hsu HH, Araki M, Mochizuki M, Hori Y, Murata M, Kahar P, Yoshida T, Hasunuma T, Kondo A. A systematic approach to time-series metabolite profiling and RNA-seq analysis of Chinese hamster ovary cell culture. Sci Rep. 2017;7:43518.
Guillon F, Philippe S, Bouchet B, Devaux MF, Frasse P, Jones B, Bouzayen M, Lahaye M. Down-regulation of an Auxin response factor in the tomato induces modification of fine pectin structure and tissue architecture. J Exp Bot. 2008;59(2):273–88.
Carrari F, Baxter C, Usadel B, Urbanczyk-Wochniak E, Zanor MI, Nunes-Nesi A, et al. Integrated analysis of metabolite and transcript levels reveals the metabolic shifts that underlie tomato fruit development and highlight regulatory aspects of metabolic network behavior. Plant Physiol. 2006;142(4):1380–96.
Carrari F, Fernie AR. Metabolic regulation underlying tomato fruit development. J Exp Bot. 2006;57(9):1883–97.
Dimster-Denk D, Rine J. Transcriptional regulation of a sterol-biosynthetic enzyme by sterol levels in Saccharomyces cerevisiae. Mol Cell Biol. 1996;16(8):3981–9.
Veen M, Stahl U, Lang C. Combined overexpression of genes of the ergosterol biosynthetic pathway leads to accumulation of sterols in Saccharomyces cerevisiae. FEMS Yeast Res. 2003;4(1):87–95.
Cardoza RE, Vizcaino JA, Hermosa MR, Sousa S, Gonzalez FJ, Llobell A, Monte E, Gutierrez S. Cloning and characterization of the erg1 gene of Trichoderma harzianum: effect of the erg1 silencing on ergosterol biosynthesis and resistance to terbinafine. Fungal Genet Biol. 2006;43(3):164–78.
Ploier B, Korber M, Schmidt C, Koch B, Leitner E, Daum G. Regulatory link between steryl ester formation and hydrolysis in the yeast Saccharomyces cerevisiae. Biochim Biophys Acta. 2015;1851(7):977–86.
Gachotte D, Barbuch R, Gaylor J, Nickel E, Bard M. Characterization of the Saccharomyces cerevisiae ERG26 gene encoding the C-3 sterol dehydrogenase (C-4 decarboxylase) involved in sterol biosynthesis. Proc Natl Acad Sci U S A. 1998;95(23):13794–9.
Swain E, Baudry K, Stukey J, McDonough V, Germann M, Nickels JT Jr. Sterol-dependent regulation of sphingolipid metabolism in Saccharomyces cerevisiae. J Biol Chem. 2002;277(29):26177–84.
Hatoh K, Izumitsu K, Morita A, Shimizu K, Ohta A, Kawai M, et al. Transformation of the mushroom species Hypsizigus marmoreus, Flammulina velutipes, and Grifola frondosa by an agrobacterium-mediated method using a universal transformation plasmid. Mycoscience. 2013;54(1):8–12.
Cho JH, Lee SE, Chang WB, Cha JS. Agrobacterium-mediated transformation of the winter mushroom, Flammulina velutipes. Mycobiology. 2006;34(2):104–7.
Ghodasara A, Voigt CA. Balancing gene expression without library construction via a reusable sRNA pool. Nucleic Acids Res. 2017;45(13):8116–27.
Yang H, Tong J, Lee CW, Ha S, Eom SH, Im YJ. Structural mechanism of ergosterol regulation by fungal sterol transcription factor Upc2. Nat Commun. 2015;6:6129.
Dunkel N, Liu TT, Barker KS, Homayouni R, Morschhauser J, Rogers PD. A gain-of-function mutation in the transcription factor Upc2p causes upregulation of ergosterol biosynthesis genes and increased fluconazole resistance in a clinical Candida albicans isolate. Eukaryot Cell. 2008;7(7):1180–90.
MacPherson S, Akache B, Weber S, De Deken X, Raymond M, Turcotte B. Candida albicans zinc cluster protein Upc2p confers resistance to antifungal drugs and is an activator of ergosterol biosynthetic genes. Antimicrob Agents Chemother. 2005;49(5):1745–52.
Zhou P, Xie W, Li A, Wang F, Yao Z, Bian Q, Zhu Y, Yu H, Ye L. Alleviation of metabolic bottleneck by combinatorial engineering enhanced astaxanthin synthesis in Saccharomyces cerevisiae. Enzym Microb Technol. 2017;100:28–36.
Yin Y, Yu G, Chen Y, Jiang S, Wang M, Jin Y, Lan X, Liang Y, Sun H. Genome-wide transcriptome and proteome analysis on different developmental stages of Cordyceps militaris. PLoS One. 2012;7(12):e51853.
Liu JY, Men JL, Chang MC, Feng CP, Yuan LG. iTRAQ-based quantitative proteome revealed metabolic changes of Flammulina velutipes mycelia in response to cold stress. J Proteome. 2017;156:75–84.
Lomsadze A, Ter-Hovhannisyan V, Chernoff YO, Borodovsky M. Gene identification in novel eukaryotic genomes by self-training algorithm. Nucleic Acids Res. 2005;33(20):6494–506.
Cock PJ, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25(11):1422–3.
Jensen LJ, Julien P, Kuhn M, von Mering C, Muller J, Doerks T, Bork P. eggNOG: automated construction and annotation of orthologous groups of genes. Nucleic Acids Res. 2008;36(Database issue):D250–4.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Yun YH, Koo JS, Kim SH, Kong WS. Cloning and expression analysis of phenylalanine ammonia-lyase gene in the mycelium and fruit body of the edible mushroom Flammulina velutipes. Mycobiology. 2015;43(3):327–32.
Huang JF, Shen ZY, Mao QL, Zhang XM, Zhang B, Wu JS, Liu ZQ, Zheng YG. Systematic analysis of bottlenecks in a multibranched and multilevel regulated pathway: the molecular fundamentals of l-methionine biosynthesis in Escherichia coli. ACS Synth Biol. 2018;7(11):2577–89.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34(2):374–8.
We are grateful to Professor Shiheng Tao for his advice and technical support. We thank Dr. Hao Wang, Dr. Chao Huang and Dr. Zefeng Wu for sharing their expertise in bioinformatic analysis, and we are grateful to Dr. Min Jia for modifying the image of the article. We would like to thank Professor Hu Tao, Assistant Professor Hongbo Guo and Yun Jiang for helping with the final versions of the manuscript. We greatly appreciate Professor Min Duan and Dr. Ningjuan Fan from the Biochemistry and Molecular Platform in College of Life Sciences, Northwest A&F University for help with providing experimental instruments and instructions.
This study was financially supported by the Science and Technology Coordination Innovation Project of Shaanxi Province (Grant No. 2015KTTSNY03–07) and the National Natural Science Foundation of China (Grant No. 31670301 and 31670295). JED designed, interpreted the data and assisted in revising the manuscript. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
This study did not require ethics approval and consent as wild, no endangered or protected plant species were involved.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The statistics of RNA-Seq data at different developmental stages of F. velutipes. Figure S1. The DEG statistics at different stages of F. velutipes. Figure S2. GO functional annotation and classification of DEGs during T1 of F. velutipes. Figure S3. GO functional annotation and classification of DEGs during T2 of F. velutipes. Figure S4. Differential expression profiles of nine DEGs at different developmental stages of F. velutipes. GAPDH and β-actin were used as internal controls. There were three biological replicates and three technical replicates for each gene. Table S2. The results of differential ions identification at different developmental stages of F. velutipes. Table S3. Identified metabolites in the ergosterol biosynthesis during T1 of F. velutipes. Table S4. Identified metabolites in the ergosterol biosynthesis during T2 of F. velutipes. Figure S5. Verification results of the end product ergosterol in metabolomics. (A) LC-MS spectrum of the ergosterol standard. (B) LC-MS spectrum of ergosterol from F. velutipes. Their corresponding m/z and retention time (RT) are indicated in the figure. Excel S1. The profiles of metabolites and gene expression at three different developmental stages of F. velutipes using Pearson’s correlation coefficient. Excel S2. Network analysis results of metabolites and genome-wide genes at three different developmental stages of F. velutipes by Cytoscape. Table S5. Primer sequences of DEGs related to ergosterol biosynthesis in F. velutipes.
About this article
Cite this article
Wang, R., Ma, P., Li, C. et al. Combining transcriptomics and metabolomics to reveal the underlying molecular mechanism of ergosterol biosynthesis during the fruiting process of Flammulina velutipes. BMC Genomics 20, 999 (2019). https://doi.org/10.1186/s12864-019-6370-1
- Flammulina velutipes
- Combined analysis
- Ergosterol biosynthesis
- Fruiting process