- Research article
- Open Access
Growth and metabolism of Beauveria bassiana spores and mycelia
BMC Microbiologyvolume 15, Article number: 267 (2015)
Fungi are ubiquitous in nature and have evolved over time to colonize a wide range of ecosystems including pest control. To date, most research has focused on the hypocrealean genera Beauveria bassiana, which is a typical filamentous fungus with a high potential for insect control. The morphology and components of fungi are important during the spores germination and outgrow to mycelia. However, to the best of our knowledge, there is no report on the morphology and components of B. bassiana spores and mycelia. In the work, the growth and metabolism of Beauveria bassiana spores and mycelia were studied. High performance liquid chromatography-mass spectrometry (HPLC-MS) was employed to study the metabolism of B. bassiana spores and mycelia. Principal component analysis (PCA) based on HPLC-MS was conducted to study the different components of the spores and mycelia of the fungus. Metabolic network was established based on HPLC-MS and KEGG database.
Through Gompertz model based on macroscopic and microscopic techniques, spore elongation length was found to increase exponentially until approximately 23.1 h after cultivation, and then growth became linear. In the metabolic network, the decrease of glyoxylate, pyruvate, fumarate, alanine, succinate, oxaloacetate, dihydrothymine, ribulose, acetylcarnitine, fructose-1, 6-bisphosphate, mycosporin glutamicol, and the increase of betaine, carnitine, ergothioneine, sphingosine, dimethyl guanosine, glycerophospholipids, and in spores indicated that the change of the metabolin can keep spores in inactive conditions, protect spores against harmful effects and survive longer.
Analysis of the metabolic pathway in which these components participate can reveal the metabolic difference between spores and mycelia, which provide the tools for understand and control the process of of spores germination and outgrow to mycelia.
Fungi are ubiquitous in nature and have evolved over time to colonize a wide range of ecosystems including pest control. Over the past century, approximately 1000 fungal species have been reported to kill insects [1–3]. To date, most research has focused on the hypocrealean genera Beauveria bassiana, which is a typical filamentous fungus with a high potential for insect control, because its spores are relatively easy and inexpensive to mass produce for field applications [4–6]. Moreover, it is known to have nontoxic effects on nontarget organisms, including animals and humans . killing insect by fungi is a process of spore germination and outgrow to mycelia , which can be concluded as follows: as soon as a fungal spore is exposed to its favourable conditions, it changes from a dormant state to an actively metabolizing cell, and then a germ tube emerges from the spore. When the germ tube reaches a certain length, the spore is considered to be germinated. In the next steps, germ tube elongation and branching take place until the mycelium forms a colony . The process has been assessed through typical reverse genetic gene-specific studies in terms of cell division , trehalose metabolism , respiration [12, 13], and nucleic and protein synthesis during spore germination . We can conclude that the morphology and components of fungi during spores outgrowth to mycelia are changed [15, 16]. However, to the best of our knowledge, there is no report on the morphology and components of B. bassiana spores and mycelia. The change of morphology and components can reveal the connection between spores and mycelia, and provide a systems-level understanding of the cell.
Despite its importance, only a limited number of methodologies have been developed for morphology and components analysis. This is primarily due to the characteristics of most components that display high polarity, nonvolatility, poor detectability, and overall similar properties . Recently, high performance liquid chromatography − mass spectrometry (HPLC-MS) equipped with electrospray ionization (ESI) detection has been used for components analysis [18–21]. It is a robust, sensitive, and selective technique, and also has become popular for quantitative and qualitative analyses. In the present study, the morphology of B. bassiana spores and mycelia were studied by combining macroscopic and microscopic techniques. And then HPLC-MS coupled with PCA were used to distinguish different metabolites of mycelia and spores. In addition, metabolic pathway was established based on HPLC-MS and KEGG database. Tracking metabolite changes under different conditions not only provides direct information on metabolism but is also complementary to gene expression and proteome analysis [22, 23]. Metabolomics, which can be defined as the measurement of the level of all intracellular metabolites, has become a powerful new tool for gaining insight into cellular function. The aim of the study was to reveal the reason of keep survive longer and infective of spores by compare significant change in metabolites between spores and mycelia. And provide the tools for understand and control the process of spores germination and outgrow to mycelia.
Results and discussion
Spore germination kinetics
The germination of spores takes place when the spores are introduced into a proper environment, which requires proper nutrition and special conditions. The spore germination can be divided into three phases: spore swelling, germ tube emergence and germ tube elongation . In the first phase, spores begin to swell to increase their dormant diameter significantly until a germ tube emerges (second phase). The two phases in early growth are supported by mobilization and utilization of storage compounds in the spores. In the third phase the elongation of the germ tube is observed, which contributes to biosynthesis and extension by uptake and metabolism of nutrients from the medium .
The spore germination kinetics was investigated in the study. The values for hyphal length were measured with the aid of Image-Pro Plus software in a series of images monitoring the growth of B. bassiana spores on PDA at 26 °C, and the duration of the germination phase was estimated. Until the 6th hour of the cultivation, no germ tubes could be spotted, although an increase in the mean diameter of spores due to swelling. (Fig. 1).
Figure 1 showed typical forms of B. bassiana spores and hyphae in their development. Tubes emerged from 8 h to approximately 11 h. About 10 h after cultivation, most of the spores had their tubes emerged. At that moment the spores entered the phase of tube elongation. The hyphae remained unbranched till the 20th hour. After 28 h, the objects became difficult to observe and analyse due to their increasing concentration and aggregation of hyphal elements.
Models were used to describe the behavior of microorganisms . In the study, the modified Gompertz model was employed to detect the relation between spore elongation length and culture time, which can provide useful information about the variance of the germination time for individual spores. Figure 2 showed changes of the most important morphological parameters. It can be seen that the modified Gompertz equation described satisfactorily the length of spore elongation over time for B. bassiana, with coefficients of determination (R2) of 0.9825. The three phases, spore swelling, germ tube emergence and germ tube elongation can be observed distinctly. In the phase of germ tube emergence and in the early hyphal development, the measured length appeared to increase exponentially until approximately 23.1 h after cultivation, when growth became linear. The above results are in agreement with the findings of Trinci [25, 26], who examined the kinetics of hyphal extension of several fungi, and considered the observed transition from exponential to linear growth can be attributed to the weakness of the hyphal tips to incorporate the increasing material that is being supplied or the deficiency of transporting material from distal hyphal regions.
Principal component analysis
PCA based on HPLC-MS data was used to study the metabolomic differences of mycelia and spores extracts from B. bassiana. PCA is an unsupervised pattern recognition method, which means that no prior knowledge concerning groups or tendencies within the data sets is necessary. PCA is usually employed to reduce the dimensionality of the data and extract essential information from large, mixed data sets . A score plot is applied for the grouping of samples by reducing the dimensionality of the data. The complex data was reduced to two principal components PC1 and PC2, which can represent most of the components. The values of the two variables for the observations are called factor scores, two factors scores can be interpreted geometrically as the projections of the observations onto the whole components. The data sets exhibiting similarities are clustered together, and those that are different are placed further apart .
Acid compounds were extracted in negative mode, while alkaline compounds were extracted in positive mode, to study the whole components extracted from spores and mycelia, both positive mode and negative mode must be employed to PCA. As shown in Fig. 3, from both positive mode and negative mode, PCA score plots employed in this study found that the mycelia were clearly separated from the spores of B. bassiana, and the spores or mycelia were grouped correctly together respectively. In the posotive mode, both principle components were significant: PC1 accounts for 35.7 % of the total variance and PC2 accounts for 52.8 %. The first two principal components (PC1 and PC2) explain more than 80 % of the total variance. In the negative mode, PC2 leads to classification of the two groups and accounts for 58.3 % of the total variance, whereas PC1 denotes 29.5 % of the total variance. The high-resolution ESI-MS data for some metabolites are shown in Fig. 4. These findings were in good agreement with HPLC–MS based metabolic profiles (Tables 1 and 2). In Table 1, twenty-eight major components were identified from B. bassiana mycelia by HPLC-MS method, eleven of which were different from that in spores, such as glycerophosphocholine, palmitic acid, linoleic acid, phosphatidylethanolamine. For spores, thirty-six compounds were extracted, seventeen of which were distinguished from mycelia, such as succinic anhydride, dihydrouracil, mannitol, sphinganine (Table 2).
Metabolic pathway analysis
Nutritional limitations or increase in cell density leads B. bassiana to produce dormant, environmentally resistant spores. The complex morphological changes that occur during sporulation are thought to be highly controlled by metabolic networks . However, no comprehensive metabolite profiling approach has been used to demonstrate alterations in large-scale metabolites. Thus, to overall unravel the effects of these metabolites on the metabolic network, HPLC-MS described above was employed to detect the metabolic pathway of some metabolites identified in mycelia and spores.
B. bassiana cells produce spores under conditions of nutrition deprivation. Interestingly, most metabolites (pyruvate, fumarate and ribulose) in the glycolytic, tricarboxylic acid (TCA) cycle and pentose phosphate pathways were markedly decreased in the stage of spores (Fig. 5). In particular, the level of fructose-1, 6-bisphosphate, a key factor in catabolite repression, dropped more rapidly in spores than in mycelia. It was possibly because that the decrease in fructose-1, 6-bisphosphate resulted in slow of metabolism, which kept spores in a dormant, metabolically inactive conditions and help them survive longer.
The pentose phosphate pathway is a biochemical pathway parallel to glycolysis, which can generates NADPH and pentoses (5-carbon sugars). While it does involve in oxidation of glucose, its primary role is anabolic rather than catabolic [30, 31]. The pathway is a major source of reductant for biosynthetic processes such as fatty-acid synthesis and assimilation of inorganic nitrogen , and maintains redox potential necessary to protect against oxidative stress . In the first stage of pentose phosphate pathway, glucose-6-phosphate was converted into ribose 5-phosphate, along with the energy production. And then ribose 5-phosphate was converted to phosphoribosyl pyrophosphate (PRPP). At the step, the metabolism was divided into three different pathways: histidine metabolism, purine metabolism and pyrimidine metabolism (Fig. 5). In the histidine metabolism, ergothioneine was extracted from B. bassiana as a metabolin, which is a naturally occurring amino acid and a thiourea derivative of histidine only produced by fungi and some prokaryotes. Ergothioneine has antioxidant properties in vitro . Under laboratory conditions, it scavenges hydroxyl radicals and hypochlorous acid, inhibits production of oxidants by metal ions, and may participate in metal ion transport and regulation of metalloenzymes [35, 36]. In 2012, Bello identified an ergothioneine biosynthetic gene Egt-1 in Neurospora crass and demonstrated that ergothioneine enhanced conidial survival and protected against peroxide toxicity during spore germination . Thus, the high levels of ergothioneine in the spores play an antioxidant role, and keep active of spores.
Dimethyl guanosine is a product of purine metabolism. Charles and his co-workers have mentioned that dimethyl guanosine as a self-inhibitor serves to regulate spore development in Dictyostelium mucoroides . Thus, the increased dimethyl guanosine in spores might diminish some gene transcription, which would make metabolism slower , and lead to improved shelf life of biological insecticide.
Dihydrothymine is a degradation product of thymine in the pyrimidine metabolism which is catalyzed by dihydropyrimidine dehydrogenase . It associated with valine, leucine and isoleucine metabolism . The decrease of dihydrothymine in spores might mean that the correlative metabolism network became slower than that in mycelia.
A mannitol cycle was first proposed by Hult and Gatenbeck , which comprises four enzymes: mannitol 1-phosphate dehydrogenase (MPDH), mannitol 1-phosphate phosphatase (MPP), NADP+-mannitol 2-dehydrogenase (MDH), and hexokinase (HX), of which MPDH is the main synthetic enzyme and MDH is the main catabolic enzyme. The roles of mannitol vary in different fungi and might act as a scavenger of reactive oxygen species . It was found only in spores (not found in mycelia) of B. bassiana, which implied that mannitol in spores might have an important impact on environment adaptability, germination, and virulence.
Glycerophospholipids, as glycerol-based phospholipids, are the main components of biological membranes, which play an important role in the generation of both extracellular and intracellular signals. As showed in Fig. 5, the increase of glycerophospholipids, esters and carnitine and the decrease of glyoxylate, pyruvate and acetylcarnitine in spores indicated that the hydrolysis of lipids and oxidation of fatty acids were depressed .
Sphingosine is a signaling lipid and may be involved in insect-fungi recognition, which can recognize the receptors on spores and induce the germination . Sphinganine begins with the condensation of serine with palmitoyl-CoA to form 3-ketosphinganine, which is rapidly reduced to sphinganine . Sphinganine is the biosynthetic precursor of sphingosines and sphingolipids, and its emergence in spores may affect fungal signaling and germination.
The TCA is a key metabolic pathway that unifies carbohydrate, fat and protein metabolism. The reactions of the cycle are carried out by eight enzymes that completely oxidize acetyl-CoA into two molecules of carbon dioxide. Through catabolism of sugars, fats and proteins, a two-carbon organic product acetate in the form of acetyl-CoA is produced which enters the citric acid cycle . An increase of carnitine and decrease of acetylcarnitine, succinate, oxaloacetate and fumarate suggested that the TCA was depressed and lipid metabolism was enhanced in the spores. As a result, lipid was accumulated to protect spores against harmful effects of environment or others.
The proposed HPLC-MS methods enable the global determination of charged species, so that they can be used as universal tools for metabolome analysis. Metabolome data, along with macroscopic and microscopic techniques were used to provide important new information of the metabolism and growth of B. bassiana spores and mycelia. Gompertz model based on macroscopic and microscopic techniques was used to detect the relation between spores elongation length and culture time. Spore elongation length was found to increase exponentially until approximately 23 h after cultivation, and then growth became linear. The results of PCA displayed clear differences of the components in mycelia and spores of the fungus. Metabolic pathway of B. bassiana spores and mycelia was established based on HPLC-MS and KEGG database, which revealed the presence of twenty-eight major components in mycelia and thirty-six compounds in spores. In the metabolic network, the decrease of glyoxylate, pyruvate, fumarate, alanine, succinate, oxaloacetate, dihydrothymine, ribulose, acetylcarnitine, fructose-1, 6-bisphosphate, mycosporin glutamicol, and the increase of betaine, carnitine, ergothioneine, sphingosine, dimethyl guanosine, glycerophospholipids, and in spores indicated that the change of the metabolin can keep spores in inactive conditions, protect spores against harmful effects and survive longer. The study provided the tools for understand and control the process of spores germination and outgrow to mycelia.
Microorganisms and culture medium
B. bassiana Bb0062 provided by the Anhui Provincial Key Laboratory of Microbial Control, Anhui Agricultural University, was stored at −20 °C in sterilized cryovials containing 10 % glycerol (in 0.02 % Tween 80 solution). It was cultured on potato dextrose agar (PDA) slants (potato 200 g/L, glucose 20 g/L and agar 20 g/L) at 26 °C for 7 days, and then stored at 4 °C until use.
Preparation of inocula
B. bassiana was cultured on PDA medium mentioned above at 26 °C for 7 days to obtain heavily sporulating cultures. Spores were then suspended in sterile distilled water containing 0.02 % (v/v) Tween 80 by gently scraping the agar surface with a sterile spatula, and then filtered through two layers of gauze to remove any debris (mostly mycelial fragments). The final spore concentration was adjusted to 1 × 106 spores/mL and used as quickly as possible.
Assessment of germination and outgrowth
Portions (100 μL) of the inoculum (B. bassiana), containing approximately 105 spores, were surface plated aseptically on PDA. After inoculation, plates were sealed with parafilm to prevent moisture loss, and stored at 26 °C for 2 days. To monitor the kinetic behavior of the B. bassiana spores for prolonged periods (from germination to mycelium formation), computer morphometry (Leica Microscopy System Ltd, DMI 4000B, Germany) was employed to examine spore germination and mycelia outgrowth. Germination time was defined as the time at which the length of the germ tube was equal to the diameter of the swollen spore. Samples were measured every 2 h until 32 h. Images was analysed by Image-Pro Plus image analysis software version 6.3 (MediaCybernetics Inc., Bethesda, United States) and an auto-focus system.
where t (h) is the time, L indicates length of spores elongation at time t, Lmax represents length of hypha t → +∞, ke signifies the slope of tangent line through the inflection point, c is the time when hyphal extention reach to maximum speed.
Mycelia and spores preparation for HPLC-MS analysis
Mycelia were harvested after 32 h of growth by a scoop, suspended in sterile water, and then filtered through two layers of gauze to remove debris, subsequently frozen in liquid nitrogen to terminate metabolism, and kept at −80 °C. Spores were collected at 7 days post-inoculation by scrapping the colony into 0.01 % Tween 80 solution, and the contents were vortexed and filtered through a 10 μm microfiltration membrane to remove any debris. The filtrate was transferred to a 50 mL centrifuge tube and centrifuged for 10 min at 8000 rpm. The precipitate was resuspended in 2 mL of distilled water and transferred to a 5 mL centrifuge tube, centrifuged for 10 min at 10,000 rpm, and then frozen in liquid nitrogen to terminate metabolism. Finally, the spores were kept at −80 °C. Three independent biological replicates were measured per assay.
Mycelia and spores were lyophilized until a constant weight was attained. Mycelia were then crushed into a fine powder and then kept at 4 °C until extraction. Ten milligrams of sample was extracted with 2 mL of 80 % methyl alcohol (Tedia company, USA, HPLC grade), followed by 1 min of vortexing and subsequent sonication (12-KHz, 8-s exposure followed by a 4-s rest interval) for 1 h. Samples were further kept at 4 °C for 12 h in the dark. After centrifugation at 8000 rpm for 10 min, 1.8 mL of supernatant was collected and dried with a centrifugal concentrator. All samples were stored at −80 °C until analysis. The dried extracts were redissolved ultrasonically in 300 μL of 80 % methanol and filtered through a 0.22 μm polyvinylidene fluoride membrane filter before HPLC − MS analysis.
HPLC − MS conditions
In order to determine retention time and to obtain extract HPLC profiles, HPLC was performed on a Thermo-Fisher UPLC system (Thermo-Fisher, SanJose, CA, USA) coupled with an LTQ XL mass spectrometer. HPLC analyses were conducted on a C18 reversed-phase (RP) column (5 μm, 3 mm × 150 mm, 100A, Luna PFP Phenomenex, Torrance, CA, USA). The parameters were as follows: injection volume, 5 μL; column temperature, 40 °C; flow rate, 0.3 mL/min; and the eluates were monitored by full-length scan from 200 to 600 nm. The mobile phase was (A) 0.1 % formic acid (Anaqua Chemical Supply, USA, HPLC grade) in water and (B) 0.1 % formic acid in acetonitrile (Merck, Germany, HPLC grade), and gradient elution was carried out: 5 % B for 0–3 min, 5–100 % B for 3–50 min, and 100 % B for 50–60 min. The mass spectrometer parameter settings used for the measurement were as follows: ionization mode, for both positive and negative; gas temperature, 350 °C; drying gas, 12 L/min; nebulizer pressure, 45 psi; capillary voltage, 4000 V in positive mode and 3500 V in negative mode; fragmentor voltage, 215 V in positive mode and 170 V in negative mode; skimmer voltage, 60 V; and OCT 1 RF, 250 V. Data acquisition was performed in the m/z range of 50–1100 Da.
Data processing and statistical analysis
All data were processed using the Xcalibur software provided by the manufacturer. After all of the detected peaks were subjected to noise-reduction in both the HPLC and MS domains, the analytical peaks were processed by the software. A list of peak intensities with retention times and m/z data pairs was generated. The intensity of each peak was normalized by the sum of all of the peak intensities. Peaks with signal-to-noise (S/N) ratios lower than 5 were rejected. PCA were performed by SPSS v18.0 (IBM SPSS Statistics, Ontario, Canada) to envisage the different components of B. bassiana spores and mycelia.
The small-molecule inventory (SMI) or metabolome is a pattern of molecules that reflects the cell’s status. The molecular formula calculated by the Xcalibur software was predicted, and based on a general understanding of fungi metabolism pathways by searching web databases (Dictionary of Natural Products, METLIN, PUBCHEM, and CHEMSPIDER). The exact monoisotopic masses of possible metabolites were calculated based on their elemental compositions. Putative biomarkers were verified by its elution order (polarity) and structure characteristics. The ambiguous metabolites were identified by comparison to authentic compounds available or referring to the published literature about fungi, especially entomopathogenic fungi.
High performance liquid chromatography-mass spectrometry
Kyoto Encyclopedia of Genes and Genomes
Mannitol 1-phosphate dehydrogenase
Mannitol 1-phosphate phosphatase
Principal component analysis
Leathers TD, Gupta CS, Alexander NJ. Mycopesticides: status, challenges and potential. J Ind Microbiol Biot. 1993;12:69–75.
Pu SC, Chen MJ, Ma ZY, Xie L, Li ZZ, Bo H. Genotyping isolates of the entomopathogenic fungus Beauveria bassiana sensu lato by multi-locus polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) analysis. Afr J Biotechnol. 2010;9:4290–4.
Xiao GH, Ying SH, Zheng P, Wang ZL, Zhang SW, Xie XQ et al. Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci Rep. 2012;2:483–92.
Tarocco F, Lecuona RE, Couto AS, Arcas JA. Optimization of erythritol and glycerol accumulation in conidia of Beauveria bassiana by solid-state fermentation, using response surface methodology. Appl Microbiol Biotechnol. 2005;68:481–8.
Taylor B, Edgington S, Luke B, Moore D. Yield and germination of the entomopathogenic fungus Beauveria bassiana when grown on different rice preparations. J Stored Prod Res. 2013;53:23–6.
Gouli V, Gouli S, Kim JS. Production of Beauveria bassiana Air conidia by means of optimization of biphasic system technology. Braz Arch Biol Technol. 2014;57:571–7.
Fang WJ, Leng B, Xiao YH, Jin K, Ma JC, Fan YH, et al. Cloning of Beauveria bassiana chitinase gene Bbchit1 and its application to improve fungal strain virulence. Appl Environ Microbiol. 2005;71:363–70.
Zimmermann G. Review on safety of the entomopathogenic fungi Beauveria bassiana and Beauveria brongniartii. Biocontrol Sci Tech. 2007;17:553–96.
Gougouli M, Koutsoumanis KP. Modeling germination of fungal spores at constant and fluctuating temperature conditions. Int J Food Microbiol. 2012;152:153–61.
Bainbridge BW. Macromolecular composition and nuclear division during spore germination in aspergillus nidulans. J Gene Microbiol. 1971;66:319–25.
Fillinger S, Chaveroche MK, Dijck PV, Vries RD, Ruijter G, Thevelein J, et al. Trehalose is required for the acquisition of tolerance to a variety of stresses in the filamentous fungus Aspergillus nidulans. Microbiology. 2001;147:1851–62.
Bradshaw RE, Bird DM, Brown S, Gardiner RE, Hirst P. Cytochrome c is not essential for viability of the fungus Aspergillus nidulans. Mol Genet Genomics. 2001;266:48–55.
Taubitz A, Bauer B, Heesemann J, Ebel F. Role of respiration in the germination process of the pathogenic mold aspergillus fumigatus. Curr Microbiol. 2007;54:354–60.
Osherov N, May G. Conidial germination in aspergillus nidulans requires RAS signaling and protein synthesis. Genetics Soc America. 2000;155:647–56.
Papagianni M. Fungal morphology and metabolite production in submerged mycelial processes. Biotechnol Adv. 2004;22:189–259.
Feofilova EP, Ivashechkinb AA, Alekhinc AI, Sergeeva YE. Fungal Spores: Dormancy, Germination, Chemical Composition, and Role in Biotechnology (Review). Appl Biochem Micro. 2012;48:1–11.
Tomoyoshi S, Yoshiaki O, Yuki U, Hisako N, Masaru T, Takaaki N. Quantitative metabolome analysis using capillary electrophoresis mass spectrometry. J Proteome Res. 2003;2:488–94.
Wagner S, Scholz K, Donegan M, Burton L, Wingate J, Volkel W. Metabonomics and biomarker discovery: LC-MS metabolic profiling and constant neutral loss scanning combined with multivariate data analysis for mercapturic acid analysis. Anal Chem. 2006;78:1296–305.
Lee JH, Kim NS, Han KM, Kim SH, Cho S, Kim WS. Monitoring by LC-MS/MS of 48 compounds of sildenafil, tadalafil, vardenafil and their analogues in illicit health food products in the Korean market advertised as enhancing male sexual performance. Food Addit Contam A. 2013;30:1849–57.
Shi XW, Zhang K, Xue N, Su LF, Ma GX, Qi JL, et al. Differentiation of genuine Inula britannica L. and substitute specimens based on the determination of 15 components using LC–MS/MS and principal components analysis. Food Chem. 2013;141:4019–25.
Napolitano A, Akay S, Mari A, Bedir E, Pizza C, Piacente S. An analytical approach based on ESI-MS, LC–MS and PCA for the quali–quantitative analysis of cycloartane derivatives in Astragalus spp. J Pharmaceut Biomed. 2013;85:46–54.
Raamsdonk LM, Teusink B, Broadhurst D, Zhang N, Hayes A, Walsh MC, et al. A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations. Nat Biotechnol. 2001;19:45–50.
Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, et al. Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science. 2001;292:929–34.
Garcia D, Ramos AJ, Sanchis V, Marín S. Predicting mycotoxins in foods: a review. Food Microbiol. 2009;26:757–69.
Trinci APJ. A kinetic study of the growth of aspergillus nidulans and other fungi. J Gen Microbiol. 1969;57:11–24.
Trinci APJ. Influence of the width of the peripheral growth zone on the radial growth rate of fungal colonies on solid media. J Gen Microbiol. 1971;67:325–44.
Camacho J, Picó J, Ferrer A. Data understanding with PCA: structural and variance information plots. Chemometr Intell Lab. 2010;100:48–56.
Pan ZZ, Gu HW, Talaty N, Chen HW, Shanaiah N, Hainline BE, et al. Principal component analysis of urine metabolites detected by NMR and DESI–MS in patients with inborn errors of metabolism. Anal Bioanal Chem. 2007;387:539–49.
Herrgård MJ, Swainston N, Dobson P, Dunn WB, Arga KY, Arvas M, et al. A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. Nat Biotechnol. 2008;26:1155–60.
Kruger NJ, Schaewen AV. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6:236–46.
Keller MA, Turchyn AV, Ralser M. Non-enzymatic glycolysis and pentose phosphate pathway-like reactions in a plausible Archean ocean. Mol Syst Biol. 2014;10:725–5.
Emes MJ, Neuhaus HE. Metabolism and transport in non-photosynthetic plastids. J Exp Bot. 1997;48:1995–2005.
Juhnke H, Krems B, Kötter P, Entian KD. Mutants that show increased sensitivity to hydrogen peroxide reveal an important role for the pentose phosphate pathway in protection of yeast against oxidative stress. Mol Gen Genet. 1996;252:456–64.
Hand CE, Honek JF. Biological chemistry of naturally occurring thiols of microbial and marine origin. J Nat Prod. 2005;68:293–308.
Misiti F, Castagnola M, Zuppi C, Giardina B, Messana I. Role of ergothioneine on S-nitrosoglutathione catabolism. Biochem J. 2001;356:799–804.
Cheah IK, Halliwell B. Ergothioneine; antioxidant potential, physiological function and role in disease. Biochim Biophys Acta. 1822;2012:784–93.
Bello MH, Barrera-Perez V, Morin D, Epstein L. The Neurospora crassa mutant NcΔEgt-1 identifies an ergothioneine biosynthetic gene and demonstrates that ergothioneine enhances conidial survival and protects against peroxide toxicity during conidial germination. Fungal Genet Biol. 2012;49:160–72.
Bacon CW, Sussman AS, Paul AG. Identification of a self-inhibitor from spores of Dictyostelium discoideum. J Bacteriol. 1973;133:1061–3.
Pluskal T, Hayashi T, Saitoh S, Fujisawa A, Yanagida M. Specific biomarkers for stochastic division patterns and starvation-induced quiescence under limited glucose levels in fission yeast. FEBS J. 2011;278:1299–315.
Nyhan WL. Disorders of purine and pyrimidine metabolism. Mol Genet Metab. 2005;86:25–33.
Furberg S, Jensen LH. The crystal and molecular structure of dihydrothymine. J Am Chem Soc. 1967;90:470–4.
Hult K, Gatenbeck S. Production of NADPH in the mannitol cycle and its relation to polyketide formation in Alternaria alternata. Eur J Biochem. 1978;88:607–12.
Voegele RT, Hahn M, Lohaus G, Link T, Heiser I, Mendgen K. Possible roles for mannitol and mannitol dehydrogenase in the biotrophic plant pathogen uromyces fabae. Plant Physiol. 2005;137:190–8.
Tsitsigiannis DI, Keller NP. Oxylipins as developmental and host–fungal communication signals. Trends Microbiol. 2007;15:109–18.
Jegorov A, Paizs B, Kuzma M, Zabka M, Landa Z, Sulc M, et al. Extraribosomal cyclic tetradepsipeptides beauverolides: profiling and modeling the fragmentation pathways. J Mass Spectrom. 2004;39:949–60.
Takabe K, Paugh SW, Milstien S, Spiegel S. “Inside-Out” signaling of sphingosine-1-phosphate: therapeutic targets. Pharmacol Rev. 2008;60:181–95.
Spiegel S, Merrill AH. Sphingolipid metabolism and cell growth regulation. Faseb J. 1996;10:1388–97.
Vagelos PR, Alberts AW, Martin DB. Studies on the mechanism of activation of acetyl coenzyme a carboxylase by citrate. J Biol Chem. 1963;238:533–40.
Zwietering MH, Jongenburger I, Rombouts FM, Riet KVT. Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990;56:1875–81.
This research was funded by Anhui Provincial Natural Science Foundation (1308085MA07) and key 863 Fund of China (2014AA021704).
The authors declare that they have no competing interests.
Conceived and designed the experiments: HXL. Performed the experiments: HXL, MXG. Measured, photographed and analyzed the experiments: HXL, HL. Analyzed the data: HXL. Contributed reagents/materials/analysis tools: HXL, XSZ, JYQ. Contributed to the writing of the manuscript: HXL, ZMZ. All authors have read and approved the final manuscript.