Skip to main content

Advertisement

Increased production of ganoderic acids by overexpression of homologous farnesyl diphosphate synthase and kinetic modeling of ganoderic acid production in Ganoderma lucidum

Article metrics

  • 417 Accesses

Abstract

Background

Ganoderic acids (GAs), derived from the medicinal mushroom Ganoderma lucidum, possess anticancer and other important pharmacological activities. To improve production of GAs, a homologous farnesyl diphosphate synthase (FPS) gene was overexpressed in G. lucidum. Moreover, the influence of FPS gene overexpression on GA production was investigated by developing the corresponding mathematical models.

Results

The maximum levels of total GAs and individual GAs (GA-T, GA-S, and GA-Me) in the transgenic strain were 2.76 mg/100 mg dry weight (DW), 41 ± 2, 21 ± 5, and 28 ± 1 μg/100 mg DW, respectively, which were increased by 2.28-, 2.27-, 2.62-, and 2.80-folds compared with those in the control. Transcription levels of squalene synthase (SQS) and lanosterol synthase (LS) genes during GA biosynthesis were upregulated by 2.28- and 1.73-folds, respectively, in the transgenic G. lucidum. In addition, the developed unstructured models had a satisfactory fit for the process of GA production in submerged cultures of G. lucidum. Analysis of the kinetic process showed that FPS gene overexpression had a stronger positive impact on GA production compared with its influence on cell growth. Also, FPS gene overexpression led to a higher non-growth-associated-constant β (1.151) over the growth-associated-constant α (0.026) in the developed models.

Conclusions

FPS gene overexpression is an effective strategy to improve the production of GAs in G. lucidum. The developed mathematical models are useful for developing a better GA production process in future large-scale bioreactors.

Background

Mushrooms are rich sources of biologically active compounds and cosmetic ingredients [1, 2]. Ganoderma lucidum, a traditional medicinal mushroom, has been used to treat and prevent various diseases for two millennia in Asian countries [3, 4]. The market size of Ganoderma products is estimated to be worth over US $2.5 billion (Bishop et al. [5]). Ganoderic acids (GAs), a type of tetracyclic triterpenoid, are one of the major ingredients of G. lucidum and have diverse pharmacological activities, including antitumor, anti-metastasis, anti-HIV, antiviral, hepatoprotective, hypocholesterolemic antioxidant, and antiaging effects [6,7,8].

GAs are biosynthesized via the mevalonate/isoprenoid (MVA) pathway from acetyl-coenzyme A to produce lanosterol in Ganoderma species (Fig. 1) [8,9,10]. Further biosynthetic steps to form GAs probably include a series of complex oxidation, reduction, hydroxylation, and acylation reactions [8, 11]. Some genes involved in GA biosynthetic pathways have been characterized in G. lucidum [9, 12,13,14,15]: 3-hydroxy-3-methyglutaryl coenzyme A reductase (HMGR), farnesyl-diphosphate synthase (FPS), squalene synthase (SQS), and lanosterol synthase (LS) have been identified as key enzymes involved in GA biosynthesis [9, 16,17,18,19].

Fig. 1
figure1

The biosynthetic pathway of ganoderic acids

Despite their important biological functions, the low production of GAs is a bottleneck for clinical trials and commercial applications [7, 20]. At present, GAs are mainly obtained from fruiting bodies and mycelia of G. lucidum. Compared to fruiting body cultures, a submerged culture is a promising alternative for the production of triterpenoids, as it is easy to control the product quality and is cost effective [21, 22]. Many attempts, such as by manipulating fermentation conditions [6, 23], developing bioprocessing and elicitor strategies [24,25,26,27,28,29], and metabolic engineering [10, 14, 15, 17, 18, 30, 31], have been conducted to increase the production of GAs by fermentation of mycelia.

Farnesyl diphosphate synthase (FPS) catalyzes the condensation reactions of dimethylallyl diphosphate with two units of isopentenyl pyrophosphate to produce farnesyl pyrophosphate, which is located at the branch point in the MVA pathway for several terpenes, like sesquiterpenes, sterols, and triterpenes [32, 33]. In hairy roots of Artemisia annua, overexpression of the FPS gene from Gossypium arboretum increases the accumulation of the sesquiterpene artemisinin [32]. FPS gene overexpression leads to significant increases in sterol and carotenoids in Nicotiana tabacum [34]. In Panax ginseng adventitious roots, overexpression of the FPS gene enhances the levels of triterpene ginsenosides [35]. The G. lucidum FPS gene has been cloned and characterized, and its transcription level exhibits a positive correlation with the triterpene content during fruiting body development [36]. These results indicate that FPS plays an important role in the control of triterpene biosynthesis. In G. lucidum, enhanced levels of GAs have been achieved by overexpression of some structural genes involved in the upstream biosynthetic pathway [14, 15, 17, 18]. However, the genetic manipulation of the FPS gene, which may regulate GA biosynthesis, has yet to be performed.

Mathematical models are important engineering tools, and they are useful for understanding bioprocess behavior, control, and optimization. Unstructured models can describe the fermentation profiles of products such as polysaccharides, gibberellic acid, triterpenoids, and flavonoids with satisfactory accuracy [24, 37,38,39]. Previously, Tang and Zhong described an unstructured mathematical model for GA biosynthesis by G. lucidum in liquid static culture conditions [40]. Feng et al. proposed an unstructured model for triterpenoid production in submerged fermentation of G. lucidum G0119 [24]. However, all those models are about fermentation of wild-type (WT) G. lucidum. Kinetic models of GA production by genetically modified strains remain to be produced. The development of such mathematical models could improve our understanding of the effects of FPS gene overexpression on the kinetics of cell growth, GA production, and sugar consumption.

In this study, the effect of FPS gene overexpression on the level of GAs and the transcription levels of key enzymes (HMGR, SQS, and LS) of the GA biosynthetic pathway were investigated in submerged cultivation of G. lucidum. Moreover, the impact of FPS gene overexpression on the production of GAs by G. lucidum was also analyzed by developing corresponding mathematical models. This study sheds light on GA biosynthesis regulation and will help us to develop efficient production of GAs on a larger scale.

Results

Generation of homologous FPS-overexpressing G. lucidum

The G. lucidum FPS gene was amplified by genomic PCR with the primers FPS-NheI-F and FPS-SmaI-R. To overexpress the FPS gene, the plasmid pJW-EXP-FPS (Additional file 1: Fig. S1) was constructed and introduced into G. lucidum using the reported genetic transformation method [10, 14, 41]. The expression of the FPS gene was driven by the glyceraldehyde-3-phosphate dehydrogenase gene promoter in the pJW-EXP-FPS plasmid. Transformants were screened on selective CYM medium with 2 μg/mL carboxin after five rounds of culture on carboxin-free medium. Those transformants maintaining antibiotic resistance were further confirmed via genomic PCR analysis, followed by sequencing of the obtained PCR products. A 1.11 kb amplification band, representing the fused fragment of the glyceraldehyde-3-phosphate dehydrogenase gene promoter and FPS gene, was detected in the lanes of the FPS transformant and the pJW-EXP-FPS plasmid, whereas no corresponding band was detected in the lane of WT G. lucidum (Additional file 1: Fig. S1). In addition, no differences were observed in morphology between the FPS transformant and WT strains (data not shown). The results of qRT-PCR analysis revealed that FPS gene overexpression induced accumulation of its transcripts in G. lucidum (Fig. 2). The transcript levels of the FPS gene in the transformant under submerged culture conditions were 3.12 and 2.71 times higher than those of the WT strain on days 6 and 9, respectively. Altogether, these results demonstrated that FPS-overexpressing G. lucidum was generated successfully, and that it gives rise to high transcription levels of the FPS gene.

Fig. 2
figure2

Transcriptional levels of GA biosynthetic genes in the WT and the FPS transgenic strain. Expression of samples from the WT strain is defined as 1.0, and expression levels in the FPS strain are displayed as fold increases over the reference sample. d days. * indicates statistical significance (P < 0.05) compared to the WT strain

Increased GA biosynthesis in G. lucidum overexpressing the homologous FPS gene

To determine whether FPS gene overexpression contributes to the accumulation of GAs in G. lucidum, the contents of total GAs and individual GAs in the WT strain and three G. lucidum transformants (namely, FPS1, FPS2, and FPS3) were analyzed in submerged cultivation conditions. Table 1 shows that the total GA levels in FPS1, FPS2, and FPS3 strains were increased to 2.28, 2.15, 2.12 times that of the WT strain (1.21 mg/100 mg DW). GA-Mk, GA-T, GA-S, and GA-Me were detected as major GA components from G. lucidum mycelia [23]. The levels of individual GAs were also measured in different G. lucidum strains. The contents of GA-Mk, GA-T, GA-S, and GA-Me in the FPS-1 transformant were 4 ± 1, 41 ± 2, 21 ± 5, and 28 ± 1 μg/100 mg DW, respectively, which were increased by 2-, 2.23-, 2.74-, and 2.85-fold compared with those of the untransformed WT strain. These results showed that overexpression of the homologous FPS gene increased GA accumulation in G. lucidum.

Table 1 Content of total GAs, GA-Mk, GA-T, GA-S and GA-Me in different cell lines

To investigate the regulatory role of the FPS gene in the GA biosynthetic pathway, the mRNA accumulations of HMGR, SQS, and LS genes were examined in the FPS transformant and WT strains by qRT-PCR. Samples were measured on days 6 and 9 in accordance with our preliminary experiments. Figure 2 shows that overexpression of the FPS gene significantly improved the transcription levels of SQS and LS genes. The maximum transcription levels of SQS and LS genes in the FPS transformant were 2.28- and 1.73- fold higher than those in the WT strain. However, no differences were detected in transcription levels of the HMGR gene between the FPS transformant and control strain. From the above information, the transformant FPS1 and the WT strain were selected to further study the kinetic profiles of cell growth, residual sugar, and total GA production.

Mathematical model development for production kinetic profiles in G. lucidum

The time-course of cell growth, sugar utilization, and GA production by G. lucidum under submerged culture conditions is shown in Fig. 3a, b. A maximum biomass of 9.30 g/L was obtained after 12 days of fermentation when the residual sugar concentration was decreased to 9.25 g/L. The total GA production increased slowly with cell growth and reached its maximum of 0.112 g/L at day 12.

Fig. 3
figure3

Simulation of GA fermentation process for biomass, residual sugar, and GA production in submerged cultivation of the WT (a, b) and the G. lucidum FPS-overexpressing strain (c, d)

These experimental data were used to determine the kinetic parameters and develop mathematical models to describe cell growth, GA production, and sugar consumption. After fitting the data for cell growth (DW, g/L), GA production (g/L), and sugar consumption (g/L) to Eqs. (1)–(3) using 1st Opt statistical analysis software, the model parameters were calculated to describe the responses as a function of time. The experimental data and model simulation of G. lucidum culture are shown in Fig. 3a, b. The calculated parameters and the correlation coefficients (R2) are summarized in Table 2. The maximum specific growth rate (μmax) and maximum biomass concentration (Xmax) were 0.45 (day−1) and 9.41 (g/L), respectively. The other parameters derived from the model were as follows: α = 0.015 (g GAs g−1 DW), β = 0 (g GAs/g−1 DW day−1), YX/S = 0.63 (g DW g−1), Ms = 0.12 (g sugar g−1 DW day−1). As evident from the results shown in Table 2, the correlation coefficient (R2) for the kinetic model was 0.99, which showed that the kinetic model gave a satisfactory fit to the experimental data. Thus, the developed mathematical models accurately described the process of GA production in submerged culture conditions and were suitable for predicting experimental results.

Table 2 Estimated parameter values of kinetic model equations

The effect of overexpression of homologous FPS on GA production through the mathematical models

To study the effect of FPS gene overexpression on GA biosynthesis, the production of total GAs was assessed in submerged cultivation of the FPS transformant. A comparison between the experimental data and simulated values was made to test the mathematical models (Fig. 3c, d). All simulated values showed a high nonlinear coefficient with the observed values (R2 = 0.99). Fitness of the experimental data with respect to F values in the models was statistically satisfactory (Fig. 3c, d and Table 2). Therefore, the data fitted the mathematical models well according to the statistical results (Table 2). These results suggested that these mathematical models can describe GA production by the FPS transformant in submerged culture conditions.

Figure 3c indicates that the FPS transformant grows at the same rate as the WT strain [18]. The maximum biomass value (Xmax) of the FPS transformant was 9.72 g/L, showing a slight increase compared with that of the WT strain (9.41 g/L). Table 2 shows that the values of μmax (day−1) and Vmax (g day−1) from the FPS transformant only had 8% and 12% increases in comparison to those from the WT strain, respectively. These results suggest that FPS gene overexpression does not significantly affect cell growth during GA production. Analysis of sugar utilization indicated that sugar consumption patterns showed a continuous decline in both the FPS transformant and WT strains during fermentation, which is consistent with the cell growth profiles. The Ms and YX/S in the FPS transformant were 0.13 (g sugar g−1 DW day−1) and 0.70 (g DW g−1 sugar), respectively, which are comparable to those obtained in the WT strain. These results indicated that overexpression of the FPS gene had no significant impact on sugar consumption. The final maximum production of GAs reached 0.209 g/L in the FPS transformant, which was higher than that in the control (0.171 g/L). FPS overexpression resulted in an improvement in productivity, with a value of 0.018 g day−1, twofold higher than that in the control (0.009 g day−1), which improved the process efficiency. The growth-associated parameter α and non-growth-associated parameter β for the FPS transformant were 0.026 (g GAs g−1 CDW) and 1.151 (g GAs g−1 CDW day−1) in submerged culture conditions, respectively. These results suggested that FPS overexpression has a positive impact on GA formation in terms of production and productivity.

Discussion

Although the FPS gene has been characterized from the basidiomycete G. lucidum, there have been no reports on the manipulation of GA biosynthesis based on the expression of the FPS gene. In this study, the homologous FPS gene was chosen for overexpression in order to increase production of GAs in G. lucidum. Moreover, an unstructured mathematical model was developed to investigate the effect of FPS gene overexpression on the production of total GAs by G. lucidum.

The FPS gene can be successfully introduced into G. lucidum as confirmed by genomic PCR and qRT-PCR analysis. FPS gene overexpression resulted in improvement of total and individual GA contents. The obtained contents of GAs in the transgenic strain were comparable to those in G. lingzhi overexpressing the squalene epoxidase gene [14]. The results indicated that FPS plays a regulatory role in GA biosynthesis. This report is in agreement with previous findings in P. ginseng and N. tabacum. Overexpression of the Centella asiatica FPS gene improved the contents of seven total ginsenosides in P. ginseng [35]. The accumulation of tetraterpene carotenoids and sterols has also been increased by heterologous expression of the FPS gene from Saccharomyces cerevisiae in tobacco [34]. Overexpression of the FPS gene may increase the metabolic flux of intermediates towards the triterpene biosynthetic pathway, thus leading to the increased accumulation of GAs in G. lucidum. The accumulation of GAs coincided with the increased transcription levels of the FPS gene, suggesting that the increased transcription level may, at least in part, be responsible for the elevation in GAs in the transgenic strain. In previous studies, a positive correlation between the transcription level of the FPS gene and the content of total triterpenoids during development was reported in G. lucidum and Poria cocos [36, 42]. Overexpression of the FPS gene induced the upregulation of SQS and LS genes rather than the HMGR gene in the transgenic strain, which is the same as previous reports showing that the manipulation of downstream genes in triterpene biosynthesis does not significantly affect the transcription levels of upstream genes [15, 18, 43].

A correlation with high significance (R2) was obtained between experimental results (for biomass, GA production, and sugar consumption) and model parameters with significant statistical fitness (F test) values. Calculated and observed time-courses exhibited good agreement in both the WT and FPS transgenic strains in this work. The developed mathematical models fit the experimental data well, and these models provide new information on the effects of FPS gene overexpression on GA production by G. lucidum with respect to the kinetic process. This information will be useful for optimizing and scaling-up of GA production in future bioreactors. Analysis of GA production showed that FPS gene overexpression had a stronger positive impact on GA production than on cell growth and consumption of sugar. Overexpression of the FPS gene led to the enhancement of GA accumulation and productivity. The growth coefficient YX/S (0.63 and 0.70 for the WT and FPS strains, respectively) obtained was comparable to that (0.78 g g−1) reported by Tang and Zhong [40]. The maximum specific growth rate of 0.448–0.483 day−1 in liquid submerged culture of G. lucidum was higher than the μmax of 0.23 day−1 in a previous study that modeled the growth profiles in liquid static culture conditions [40]. The longer fermentation time in the liquid static culture model may be responsible for the difference in the specific growth rate. The obtained higher α value over β (α = 0.015 g GAs g−1 DW and β = 0 g GAs g−1 DW day−1) in the WT strain indicated that a growth-associated factor plays an important role in GA production. These results are consistent with previous reports in liquid-submerged cultivation of G. lucidum G0119 and in liquid static culture of G. lucidum [24, 40]. However, the much higher value of β over α (α = 0.026 g GAs g−1 DW and β = 1.151 g GAs g−1 DW day−1) in the transgenic strain suggested that GA production is more related to the obtained biomass than to the rate of cell growth. This may be due to the higher content of total GAs in G. lucidum overexpressing the FPS gene.

Conclusions

Overexpression of the homologous FPS gene upregulated the transcription levels of SQS and LS genes, and increased production of GAs in a submerged G. lucidum culture. Moreover, empirical mathematical models were developed to explore the effect of FPS gene overexpression on GA production. The results indicated that FPS gene overexpression is an effective strategy to improve the production of GAs in G. lucidum. The developed mathematical models may be useful for developing a better GA production process in future large-scale bioreactors.

Methods

Strains and culture conditions

Ganoderma lucidum CGMCC 5.616 culture was purchased from the China General Microbiological Culture Collection Center. Escherichia coli strain JM109 is routinely grown in our laboratory. G. lucidum was cultured as described by Xu et al. [23]. Protoplasts of G. lucidum were regenerated in CYM media that contained 20 g/L glucose, 10 g/L maltose, 0.6 M mannitol, 2 g/L yeast extract, 2 g/L tryptone, 0.5 g/L MgSO4, 4.6 g/L KH2PO4, and 10 g/L agar.

Vector construction and transformation

The FPS gene of G. lucidum (GenBank accession number: EU399544) was acquired by genomic PCR using primers FPS-NheI-F (5′-GCTAGCATGGCCGATGCAAAGGCTCAG-3′) and FPS-SmaI-R (5′-CCCGGGTCACTTCTGCCGCTTGTAGATCTTGTC-3′). The G. lucidum FPS gene was double-digested with NheI and SmaI, and then ligated into the pJW-EXP vector [41], forming pJW-EXP-FPS. G. lucidum protoplasts were transformed with 10 μg of pJW-EXP-FPS by polyethylene glycol-mediated genetic transformation as described previously [10, 41, 44]. Transformants were plated on CYM regeneration plates containing 2 μg carboxin/mL to isolate resistant clones. The genomic DNA of G. lucidum was extracted using the cetyltrimethylammonium bromide (CTAB) method. Stable transformants were confirmed by genomic PCR analysis using the primer pair GPD-FPS-F (5′-CGTAGCAATGCCAGGAAA-3′) and GPD-FPS-R (5′-TTGAGGATCTCGACGGAGT-3′) and sequencing.

Gene expression analysis by quantitative real-time PCR (qRT-PCR)

Total RNAs of G. lucidum were extracted using TriZol Reagent (Invitrogen, Carlsbad, CA), treated with DNAase I (Fermentas, Canada), and then reverse-transcribed to cDNA using a Superscript RNAase H-First-strand synthesis kit (Invitrogen). The transcriptional levels of HMGR, FPS, SQS, and LS genes were determined by quantitative real-time PCR. Primers and PCR protocols used in the quantitative real-time PCR were described previously [17, 30]. The transcription levels of HMGR, FPS, SQS, and LS genes were normalized relative to the transcription level of the 18S rRNA gene. The transcription levels of samples from the wild-type (WT) strain were defined as 1.0, and the transcription levels in the transgenic strain are presented as fold changes with respect to these reference levels. Post-quantitative real-time PCR calculations were performed according to the 2−ΔΔCt method.

Analysis of biomass, total GAs and individual GAs

Ganoderma lucidum mycelia were collected during cultivation and were centrifuged at 10,000×g for 10 min, washed three times with water, and dried at 45 °C. Biomass was determined by the gravimetric method. The total GAs and individual GAs of mycelia were extracted and measured using methods described in our previous work [17, 23].

Statistical analysis

Data are the averages of three independent sample measurements. The results were analyzed for statistical significance by Student’s t test. The differences were considered significant at P < 0.05 in a two-tailed analysis. All data are represented as mean ± standard deviation. The 1st Opt statistical analysis software (http://www.7d-soft.com) was used to fit the developed models and estimate the parameters based on universal global optimization (UGO) algorithms. 1st Opt can be used to preform nonlinear regression, curve fitting, and for the optimization of nonlinear model parameters [45].

Mathematical model development

Mycelia growth

For mycelia growth, a logistic model was used to describe the growth of G. lucidum:

$$\frac{{{\text{d}}X}}{{{\text{d}}t}} = \mu_{\text{max} } \times \left( {1 - \frac{X}{{X_{\text{max} } }}} \right)$$
(1)

where \(\frac{{{\text{d}}X}}{{{\text{d}}t}}\) is the rate of cell growth, X is the cell dry weight (g DW L−1), μmax is the maximum specific growth rate (day−1), and Xmax is the maximum cell dry weight (g DW L−1).

Production of total GAs

The kinetics of the formation of total GAs were simulated by the Luedeking–Piret model:

$$\frac{{{\text{d}}P}}{{{\text{d}}t}} = \alpha \frac{{{\text{d}}X}}{{{\text{d}}t}} + \beta X,$$
(2)

where \(\frac{{{\text{d}}P}}{{{\text{d}}t}}\) is the GA formation rate, P is the production of total GAs (mg L−1), α is the constant for growth-associated kinetics of GA production (g GAs g−1 DW), and β is the non-growth-associated-constant (g GAs/g−1 DW day−1).

Sugar consumption

The sugar consumption equation is a Luedeking–Piret-like equation in which the amount of sugar used for GA production is assumed to be negligible:

$$- \frac{{{\text{d}}S}}{{{\text{d}}t}} = \frac{1}{{Y_{X/S} }}\frac{{{\text{d}}X}}{{{\text{d}}t}} + M_{\text{s}} X ,$$
(3)

where \(- \frac{{{\text{d}}S}}{{{\text{d}}t}}\) is the consumption rate of sugar, S is the residual sugar concentration (g L−1), YX/S (g DW g−1 sugar) is the maximum cell growth coefficient, and Ms is the maintenance constant (g sugar g−1 DW day−1).

Availability of data and materials

Identification and characterization of the FPS gene-overexpressing strain.

Abbreviations

GA:

ganoderic acid

FPS:

farnesyl-diphosphate synthase

WT:

wild-type

DW:

dry cell weight

HMGR:

3-hydroxy-3-methylglutaryl coenzyme A

SQS:

squalene synthase

LS:

lanosterol synthase

CYM:

complete yeast medium

gpd :

glyceraldehyde-3-phosphate dehydrogenase gene

UGO:

universal global optimization

References

  1. 1.

    Chaturvedi VK, Agarwal S, Gupta KK, Ramteke P, Singh MP. Medicinal mushroom: boon for therapeutic applications. 3 Biotech. 2018;8:334.

  2. 2.

    Taofiq O, Gonzalez-Paramas AM, Martins A, Barreiro MF, Ferreira ICFR. Mushrooms extracts and compounds in cosmetics, cosmeceuticals and nutricosmetics—a review. Ind Crops Prod. 2016;90:38–48.

  3. 3.

    Hsu KD, Cheng KC. From nutraceutical to clinical trial: frontiers in Ganoderma development. Appl Microbiol Biotechnol. 2018;102:9037–51.

  4. 4.

    Paterson RR. Ganoderma-a therapeutic fungal biofactory. Phytochemistry. 2006;67:1985–2001.

  5. 5.

    Bishop KS, Kao CHJ, Xu Y, Glucina MP, Paterson RRM, Ferguson LR. From 2000 years of Ganoderma lucidum to recent developments in nutraceuticals. Phytochemistry. 2015;114:56–65.

  6. 6.

    Upadhyay M, Shrivastava B, Jain A, Kidwai M, Kumar S, Gomes J, Goswami DG, Panda AK, Kuhad RC. Production of ganoderic acid by Ganoderma lucidum RCKB-2010 and its therapeutic potential. Ann Microbiol. 2014;64:839–46.

  7. 7.

    Xiao H, Zhong JJ. Production of useful terpenoids by higher-fungus cell factory and synthetic biology approaches. Trends Biotechnol. 2016;34:242–55.

  8. 8.

    Xu JW, Zhao W, Zhong JJ. Biotechnological production and application of ganoderic acids. Appl Microbiol Biotechnol. 2010;87:457–66.

  9. 9.

    Shi L, Ren A, Mu DS, Zhao MW. Current progress in the study on biosynthesis and regulation of ganoderic acids. Appl Microbiol Biotechnol. 2010;88:1243–51.

  10. 10.

    Xu JW, Zhong JJ. Genetic engineering of Ganoderma lucidum for the efficient production of ganoderic acids. Bioengineered. 2015;6:357–60.

  11. 11.

    Chen SL, Xu J, Liu C, Zhu YJ, Nelson DR, Zhou SG, Li CF, Wang LZ, Guo X, Sun YZ, Luo HM, Li Y, Song JY, Henrissat B, Levasseur A, Qian J, Li JQ, Luo X, Shi LC, He L, Xiang L, Xu XL, Niu YY, Li QS, Han MV, Yan HX, Zhang J, Chen HM, Lv AP, Wang Z, Liu MZ, Schwartz DC, Sun C. Genome sequence of the model medicinal mushroom Ganoderma lucidum. Nat Commun. 2012;3:913.

  12. 12.

    Wang WF, Xiao H, Zhong JJ. Biosynthesis of a ganoderic acid in Saccharomyces cerevisiae by expressing a cytochrome P450 gene from Ganoderma lucidum. Biotechnol Bioeng. 2018;115:1842–54.

  13. 13.

    Yang CS, Li WC, Li C, Zhou ZH, Xiao YL, Yang X. Metabolism of ganoderic acids by a Ganoderma lucidum cytochrome P450 and the 3-keto sterol reductase ERG27 from yeast. Phytochemistry. 2018;155:83–92.

  14. 14.

    Zhang DH, Jiang LX, Li N, Yu XY, Zhao P, Li T, Xu JW. Overexpression of the squalene epoxidase gene alone and in combination with the 3-hydroxy-3-methylglutaryl coenzyme A gene increases ganoderic acid production in Ganoderma lingzhi. J Agr Food Chem. 2017;65:4683–90.

  15. 15.

    Zhou JS, Ji SL, Ren MF, He YL, Jing XR, Xu JW. Enhanced accumulation of individual ganoderic acids in a submerged culture of Ganoderma lucidum by the overexpression of squalene synthase gene. Biochem Eng J. 2014;90:178–83.

  16. 16.

    Hu YR, Ahmed S, Li JW, Luo BB, Gao ZY, Zhang QY, Li XH, Hu XB. Improved ganoderic acids production in Ganoderma lucidum by wood decaying components. Sci Rep. 2017;7:46623.

  17. 17.

    Xu JW, Xu YN, Zhong JJ. Enhancement of ganoderic acid accumulation by overexpression of an N-terminally truncated 3-hydroxy-3-methylglutaryl coenzyme A reductase gene in the basidiomycete Ganoderma lucidum. Appl Environ Microbiol. 2012;78:7968–76.

  18. 18.

    Zhang DH, Li N, Yu XY, Zhao P, Li T, Xu JW. Overexpression of the homologous lanosterol synthase gene in ganoderic acid biosynthesis in Ganoderma lingzhi. Phytochemistry. 2017;134:46–53.

  19. 19.

    Zhao W, Xu JW, Zhong JJ. Enhanced production of ganoderic acids in static liquid culture of Ganoderma lucidum under nitrogen-limiting conditions. Bioresour Technol. 2011;102:8185–90.

  20. 20.

    Xu YN, Zhong JJ. Impacts of calcium signal transduction on the fermentation production of antitumor ganoderic acids by medicinal mushroom Ganoderma lucidum. Biotechnol Adv. 2012;30:1301–8.

  21. 21.

    Wagner R, Mitchell DA, Sassaki GL, De Almeida Amazonas MAL, Berovic M. Current techniques for the cultivation of Ganoderma lucidum for the production of biomass, ganoderic acid and polysaccharides. Food Technol Biotechnol. 2003;41:371–82.

  22. 22.

    Zhong JJ, Xiao JH. Secondary metabolites from higher fungi: discovery, bioactivity, and bioproduction. Adv Biochem Eng Biotechnol. 2009;113:79–150.

  23. 23.

    Xu JW, Xu YN, Zhong JJ. Production of individual ganoderic acids and expression of biosynthetic genes in liquid static and shaking cultures of Ganoderma lucidum. Appl Microbiol Biotechnol. 2010;85:941–8.

  24. 24.

    Feng J, Zhang JS, Jia W, Yang Y, Liu F, Lin CC. An unstructured kinetic model for the improvement of triterpenes production by Ganoderma lucidum G0119 based on nitrogen source effect. Biotechnol Bioproc Eng. 2014;19:727–32.

  25. 25.

    Ren A, Qin L, Shi L, Dong X, Mu DS, Li YX, Zhao MW. Methyl jasmonate induces ganoderic acid biosynthesis in the basidiomycetous fungus Ganoderma lucidum. Bioresour Technol. 2010;101:6785–90.

  26. 26.

    You BJ, Tien N, Lee MH, Bao BY, Wu YS, Hu TC, Lee HZ. Induction of apoptosis and ganoderic acid biosynthesis by cAMP signaling in Ganoderma lucidum. Sci Rep. 2017;7:318.

  27. 27.

    Wang XL, Yang HL, Liu GQ. Enhanced triterpene acid production by Ganoderma lucidum using a feeding stimulus integrated with a two-stage pH-control strategy. J Chem Technol Biotechnol. 2016;91:2974–82.

  28. 28.

    Zhang W, Tang YJ. A novel three-stage light irradiation strategy in the submerged fermentation of medicinal mushroom Ganoderma lucidum for the efficient production of ganoderic acid and Ganoderma polysaccharids. Biotechnol Prog. 2008;24:1249–61.

  29. 29.

    Zhang JM, Zhong JJ, Geng AL. Improvement of ganoderic acid production by fermentation of Ganoderma lucidum with cellulase as an elicitor. Process Biochem. 2014;249:1580–6.

  30. 30.

    Li HJ, He YL, Zhang DH, Yue TH, Jiang LX, Li N, Xu JW. Enhancement of ganoderic acid production by constitutively expressing Vitreoscilla hemoglobin gene in Ganoderma lucidum. J Biotechnol. 2016;227:35–40.

  31. 31.

    Xu JW, Yue TH, Yu X, Zhao P, Li T, Li N. Enhanced production of individual ganoderic acids by integrating Vitreoscilla haemoglobin expression and calcium ion induction in liquid static cultures of Ganoderma lingzhi. Microb Biotechnol. 2019. https://doi.org/10.1111/1751-7915.1338.

  32. 32.

    Chen DH, Liu CJ, Ye HC, Li GF, Liu BY, Meng YL, Chen XY. Ri-mediated transformation of Artemisia annua with a recombinant farnesyl diphosphate synthase gene for artemisinin production. Plant Cell Tiss Org. 1999;57:157–62.

  33. 33.

    Dhar MK, Koul A, Kaul S. Farnesyl pyrophosphate synthase: a key enzyme in isoprenoid biosynthetic pathway and potential molecular target for drug development. New Biotechnol. 2013;30:114–23.

  34. 34.

    Daudonnet S, Karst F, Tourte Y. Expression of the farnesyldiphosphate synthase gene of Saccharomyces cerevisiae in tobacco. Mol Breed. 1997;3:137–45.

  35. 35.

    Park HW, Kim OT, Hyun DY, Kim YB, Kim JU, Kim YC, Bang KH, Cha SW, Choi JE. Overexpression of farnesyl diphosphate dynthase by introducing CaFPS gene in Panax ginseng C. A. Mey. Korean J Med Crop Sci. 2013;21:32–8.

  36. 36.

    Ding YX, Yang XO, Shang CH, Ren A, Shi L, Li Y, Zhao MW. Molecular cloning, characterization, and differential expression of a farnesyl diphosphate synthase gene from the basidiomycetous fungus Ganoderma lucidum. Biosci Biotechnol Biochem. 2008;72:1571–9.

  37. 37.

    Ma XK, Peterson EC, Rung TT, Duan XY. The influence of naphthaleneacetic acid (NAA) and coumarin on flavonoid production by fungus Phellinus sp.: modeling of production kinetic profiles. Appl Microbiol Biotechnol. 2015;99:9417–26.

  38. 38.

    Shukla R, Chand S, Srivastava AK. Batch kinetics and modeling of gibberellic acid production by Gibberella fujikuroi. Enzyme Microbio Technol. 2005;36:492–7.

  39. 39.

    Xu DQ, Fu ML, Chen QH, Liu J. Effect of submerged culture conditions on exopolysaccharides production by Armillaria luteo-virens Sacc QH and kinetic modeling. Bioprocess Biosyst Eng. 2011;34:103–11.

  40. 40.

    Tang YJ, Zhong JJ. Modeling the kinetics of cell growth and ganoderic acid production in liquid static cultures of the medicinal mushroom Ganoderma lucidum. Biochem Eng J. 2004;21:259–64.

  41. 41.

    Yu X, Ji SL, He YL, Ren MF, Xu JW. Development of an expression plasmid and its use in genetic manipulation of Lingzhi or Reishi medicinal mushroom, Ganoderma lucidum (higher basidiomycetes). Int J Med Mushrooms. 2014;16:161–8.

  42. 42.

    Wang JR, Li YY, Liu D. Cloning and characterization of farnesyl diphosphate synthase gene involved in triterpenoids biosynthesis from Poria cocos. Int J Mol Sci. 2014;15:2218–22202.

  43. 43.

    Han JY, In JG, Kwon YS, Choi YE. Regulation of ginsenoside and phytosterol biosynthesis by RNA interferences of squalene epoxidase gene in Panax ginseng. Phytochemistry. 2010;71:36–46.

  44. 44.

    Xu JW, Ji SL, Li HJ, Ren MF, Duan YQ, Dang LZ, Mo MH. Increased polysaccharide production and biosynthetic gene expressions in a submerged culture of Ganoderma lucidum by the overexpression of the homologous alpha-phosphoglucomutase gene. Bioproc Biosys Eng. 2015;38:399–405.

  45. 45.

    Zhang HJ, Liu JT, Gao YJ, Wang YT. Effects of particle size on lignite reverse flotation kinetics in the presence of sodium chloride. Powder Technol. 2013;246:658–63.

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (Nos. 21566016, 81860668 and 11561071) and Yunnan Applied Basic Research Project (2018FB065).

Author information

JWX and YF designed the study. YF, NL and DHZ conducted the experiments. JWX, YF and NL analyzed the data. YF, NL and JWX wrote the manuscript. All authors read and approved the final manuscript.

Correspondence to Jun-Wei Xu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

12934_2019_1164_MOESM1_ESM.docx

Additional file 1: Figure S1. (A) The expression vector pJW-EXP-FPS. (B) Identification and characterization of the FPF gene-overexpressing strain. Amplification pattern of genomic PCR obtained with primers for the fusion of the gpd promoter and FPS gene fragment from different strains.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Ganoderma lucidum
  • Ganoderic acid
  • Farnesyl diphosphate synthase
  • Genetic engineering
  • Modeling
  • Production kinetics