Next Article in Journal
Molecular and Cellular Mechanisms that Induce Arterial Calcification by Indoxyl Sulfate and P-Cresyl Sulfate
Previous Article in Journal
A Polyphasic Approach to Compare the Genomic Profiles of Aflatoxigenic and Non-Aflatoxigenic Isolates of Aspergillus Section Flavi
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gas Chromatography–Mass Spectrometry Profiling of Volatile Compounds Reveals Metabolic Changes in a Non-Aflatoxigenic Aspergillus flavus Induced by 5-Azacytidine

1
School of Life Science, Jiangsu Normal University, Xuzhou 221116, China
2
Key Laboratory of Pathogenic Fungi and Mycotoxins of Fujian Province, Key Laboratory of Biopesticide and Chemical Biology of Education Ministry, and School of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China
*
Authors to whom correspondence should be addressed.
These authors contributed to this work equally.
Submission received: 11 December 2019 / Revised: 8 January 2020 / Accepted: 16 January 2020 / Published: 19 January 2020
(This article belongs to the Section Mycotoxins)

Abstract

:
Aspergillus flavus is one of the most opportunistic pathogens invading many important oilseed crops and foodstuffs with such toxic secondary metabolites as aflatoxin (AF) and Cyclopiazonic acid. We previously used the DNA methylation inhibitor 5-azacytidine to treat with an AF-producing A. flavus A133 strain, and isolated a mutant (NT) of A. flavus, which displayed impaired abilities of AF biosynthesis and fungal development. In this study, gas chromatography–mass spectrometry (GC-MS) analysis was used to reveal the metabolic changes between these two strains. A total of 1181 volatiles were identified in these two strains, among which 490 volatiles were found in these two strains in vitro and 332 volatiles were found in vivo. The NT mutant was found to produce decreasing volatile compounds, among which most of the fatty acid-derived volatiles were significantly downregulated in the NT mutant compared to the A133 strain, which are important precursors for AF biosynthesis. Two antioxidants and most of the amino acids derived volatiles were found significantly upregulated in the NT mutant. Overall, our results reveal the difference of metabolic profiles in two different A. flavus isolates, which may provide valuable information for controlling infections of this fungal pathogen.
Key Contribution: DNA methylation inhibitor induces defects in A. flavus development and AF biosynthesis. The NT mutant produces decreasing volatile compounds. Fatty acid-derived volatiles are decreasing in the NT mutant. Antioxidants and most of the amino acids derived volatiles were upregulated in NT.

1. Introduction

Aspergillus flavus is a soil-born crossover pathogen that could opportunistically invade many important oilseed crops and immunosuppressed patients [1,2]. Importantly, the toxigenic A. flavus produces various harmful secondary metabolites (SMs) when spoiling foodstuffs, among which aflatoxins (AFs) are the most notorious and carcinogenic mycotoxin [3,4]. Foodstuff contaminated with a low dose of AFs that is ingested chronically by animal or human could pose liver cancer, while ingestion of a high dose of AFs contaminating foodstuffs is poisonous and even fatal [5,6,7,8]. Adding to AFs, cyclopiazonic acid (CPA), aflatrem, and kojic acid (KA) were also identified in A. flavus [9,10]. These SMs produced by A. flavus are small bioactive molecules that exhibit harmful for animals and humans [11]. The whole genome of A. flavus has been sequenced [12,13], which reveals 56 secondary metabolite clusters that are regulated by different environmental factors [14]. Although each SM cluster can potentially produce specific metabolites, only the molecular structures of AFs, CPA, aflatrem, and KA are identified in this fungus [15].
SMs are rarely found to be involved in fungal development or reproduction, but they do play critical roles in fungal ecology in response to environmental factors [16]. In Aspergillus, although some SMs exhibit harmful properties (like AFs and CPA), some other SMs like lovastatin also display beneficial properties for humans [9,17]. The SM genes, which are necessary for synthesizing a specific compound, are usually clustered in the genome. In filamentous fungi, most of the SM clusters are located proximately to the telomere or heterochromatin, where they are subject to specific silencing mechanisms [18]. The silencing mechanisms in these chromatin regions are associated with histone acetylation and DNA methylation [19]. In filamentous fungi, the inhibitors of histone deacetylases (HDACs) and DNA methyltransferases (DNMTs) have been utilized to inhibit gene silencing in the genome and activate the expression of SM gene clusters [19,20,21]. The inhibitors of the second class HDACs, such as trichostatin A (TSA), SAHA (vorinostat), and suberoyl bis-hydroxamic acid (SBHA) [19,21,22], have been utilized to successfully identify new SMs in fungi. TSA was the first HDAC inhibitor used to increase SM production in Aspergillus [21]. In A. nidulans, the treatment of 5-methylmellein, an inhibitor of HDAC SirA produced by Didymobotryum rigidum, promoted SM profiles [23].
The DNA methyltransferases inhibitors 5-azacytidine and RG-108 have also been used to identify novel metabolites in fungi [20,24]. Although the Aspergillus species exhibit low concentration of DNA methylation in their genomes, our former study has demonstrated that the use of 5-azacytidine in A. flavus resulted in a dramatic reduction of AF production and concurrent developmental defects [25]. Additionally, a non-aflatoxigenic mutant (NT) of A. flavus induced by 5-azacytidine was also isolated, which displayed impaired abilities of AF biosynthesis and fungal development. Fungal volatiles are important cellular metabolites, which hold a considerable portion of the microbial metabolome [17]. In this study, we performed comparative volatile profiling to analyze the information on intracellular metabolism in an aflatoxigenic A. flavus A133 strain and a non-aflatoxigenic mutant (NT) induced by 5-azacytidine by gas chromatography–mass spectrometry (GC-MS). A total of 1181 volatiles were identified in A133 or NT strains, among which 490 volatiles were found in these two strains in vitro and 332 volatiles were found in vivo.

2. Results

2.1. The NT Mutant Induced by 5-Azacytidine Is Impaired in Conidiation and Mycotoxins Production

DNA methylation, an important epigenetic modification, plays a critical role in the regulation of gene transcription involving in development and secondary metabolites in many plants and fungi. Our former study has demonstrated that the treatment of 5-azacytidine, a DNA methyltransferase inhibitor, blocks AF biosynthesis and fungal development in A. flavus [25]. We also obtained a non-aflatoxigenic mutant (NT) of A. flavus by 5-azacytidine treatment, which failed to produce AF. To better understand the changes in this mutant induced by 5-azacytidine, we both examined its developmental phenotype and mycotoxins production compared to the A133 wild-type strain. As shown in Figure 1A, the NT strain produced white fluffy mycelium and failed to produce conidia pigment, while the A133 strain produced dense dark green spores. Further examination of aerial conidiophores by microscope also showed that the NT strain produced fewer conidial heads and could not form normal conidiophores compared to the A133 strain (Figure 1B). The abilities of mycotoxins production of the NT strain in a Czapek liquid medium were also determined. As shown in Figure 1C, the A133 strain produced yellow mycelial pellets and grew with adherence, while the NT strain produced white and bigger mycelial pellets and grew without adherently. We then extracted their secondary metabolites (SMs) from the cultures with chloroform, and the result showed that the extracts of the A133 strain displayed yellow in color, while the NT strain was colorless. Then the SMs extracts were separated using thin-layer chromatography (TLC), which showed that the A133 strain was able to produce aflatoxin B1 (AFB1) and AFB2 both in vivo (volatiles extracted from the mycelium of A. flavus) and in vitro (volatiles emitted from the A. flavus), while non detectable aflatoxins were found in the NT strain in vitro (Figure 1D). Intriguingly, AFB2 but not AFB1 was found in the NT strain in vivo (Figure 1D).

2.2. GC-MS Analysis of Volatiles Extracts

Fungal volatiles are important cellular metabolites which hold a considerable portion of the microbial metabolome. To analyze the diversity of volatile profiles within the A133 (WT) and NT strains, GC-MS analysis was performed to detect the volatiles emitted from these two strains (in vitro) or extracted from their mycelium (in vivo). A representative total ion chromatogram (TIC) was shown in Figure 2A (in vitro) and Figure 2B (in vivo), and indicates the different volatile profiles among the A133 and NT strains. The resulted GC-MS profiles were subjected to XCMS online program (https://xcmsonline.scripps.edu/) for further analysis. As shown in Figure 3A,B, the results of principal component analysis (PCA) demonstrated that NT mutants clustered separately with A133 wild-type strain, which implies that they have specific metabolic profiles under aflatoxin-producing condition in comparison to wild type (PCA 1, 90%), as well as to each other (PCA 2, 3%). Metabolite with p value ≤ 0.01 and fold change ≥ 1.5 in each comparison was selected as a changed metabolite (Figure 3C,D).

2.3. The A. flavus Non-Aflatoxigenic Mutant Produces Decreasing Volatile Compounds

The NT mutant was compared to the wild type both in vitro and in vivo within the XCMS program and visualized using volcano plots according to log10 of p value (y-axis) and log2 of fold change (x-axis) (Figure 4A,B). Here we identified a total of 849 and 691 volatiles in A133 or NT strains in vitro and in vivo, respectively (Figure 4C), among which 490 specific volatiles were identified in these two strains in vitro and 332 specific volatiles in vivo. Intriguingly, the number of volatile compounds observed in the NT strain (531/535) was reduced compared to the wild-type strain (567/559) in vitro and in vivo (Figure 4A,B).

2.4. Screening of the Changed Volatile Compounds in NT Mutant

The identified metabolites were clustered and visualized with a heat map, and the results showed that, among the 849 volatile metabolites emitted extracellularly from these two strains (in vitro), 376 compounds showed increasing expression in the A133 strain compared to the NT mutant, while 473 compounds were significantly up-regulated in the NT mutants compared to the A133 strain (Figure S1A in Supplementary Materials). A total of 691 volatile metabolites were identified from the mycelium of these two strains (in vivo), among which 450 compounds were significantly up-regulated in the A133 strain compared to the NT mutant, while 241 compounds were shown to be over-expressed in the NT mutants compared to the A133 strain (Figure S1B). The top 30 volatile metabolites produced by the A133 strain and NT strain in vitro, and the top 20 volatile metabolites in vivo, differing according to t-test were also clustered and visualized with a heat map (Figure 5A,B). A total of 21 compounds were up-regulated in the A133 strain among the top 30 metabolites, and 9 compounds were found increasing expression in NT mutants (Figure 5A). These top differential expressed metabolites were classified, which included fatty, amino, organic acid, as well as the saccharides. The fatty acid-derived volatiles and most of the phenylpropanoid derivatives were significantly up-regulated in the A133 strain compared to the NT strain (Table 1 and Table 2). The A133 strain accumulated extremely higher levels of n-Hexadecanoic acid (more than 2846 fold levels vs. in the NT mutant), 9,12-Octadecadienoicacid(Z,Z)-(more than 1980 fold levels vs. in the NT mutant), Ergosterol (more than 106 fold levels vs. in the NT mutant) and Flopropione (more than 140 fold levels vs. in the NT mutant) extracellularly. Among these top metabolites, M69T60 (2,6,10,14,18,22-Tetracosahexaene,2,6,10,15, 19,23-hexamethyl-,(all-E)-) was the only triterpene identified in the A133 strain, which was detected with more than 82 fold levels in A133 vs. in the NT mutant (Table 1). Most of the differential expressed amino acids derived volatiles were detected highly expressed in the NT mutant both in vitro and in vivo with more than 66 fold levels vs in the A133 strain (Table 1 and Table 2). Among the top differential expressed metabolites extracted from the fungal mycelium (in vivo), accumulation of three phenylpropanoid derivatives (4 in total) and one saccharide derived volatile (D-Glucoside), the only saccharide derived volatile identified in these two strain in vivo, were found in the NT mutant compared to the A133 strain (Table 2). In the NT mutant, one of the phenylpropanoid derivatives, M353T55 (Phenol,4,4′-methylenebis [2,6-bis (1,1-dimethylethyl)]-) was detected with more than 9 fold levels vs. the A133 strain, which is considered as an antioxidant.
The volatile metabolites were also analyzed with the Agilent ChemStation software, which showed that a total number of 84 compounds emitted extracellularly were identified from these two A. flavus strains (in vitro), and a total amount of 54 compounds were identified in vivo from the mycelium of these two strains (Tables S1 and S2). The A133 strain and NT mutant shared 18 and 28 same compounds produced in vitro and in vivo, respectively. Among these identified volatile metabolites, most of the fatty acid-derived volatiles were up-regulated in the A133 strain produced in vitro (Table S1). While most of the amino acid-derived volatiles, intriguingly, were detected increasingly in the NT mutant (Tables S1 and S2). The properties of the metabolites were identified both in the websites of PubChem database (https://pubchem.ncbi.nlm.nih.gov) and Guidechem (https://www.guidechem.com/). Here, we found that two metabolites, Phenol,2,2′-methylenebis[6-(1,1-dimethylethyl)-4-methyl-] (CAS: 000119-47-1) and Phenol,4,4′-methylenebis[2,6-bis(1,1-dimethylethyl)-] (CAS: 000118-82-1), both of which were considered as antioxidants, were found significantly upregulated in the NT mutant (Table S2), which is consistent with the result analyzed by the XCMS program.

3. Discussion

A. flavus is a notorious crossover pathogen, disseminating in agricultural soils, that could not only spoil crop seeds and foodstuffs but also be a pathogen of immunosuppressed patients [2]. The contamination of aflatoxins (AFs) and other mycotoxins produced by this fungus on foodstuff or crop seeds has been the cause of significant concern within the food industry and farmers [26,27]. This fungus is hard to eliminate from the environment as the formation of resistant overwintering structures, sclerotia. Our former work indicated that treatment with the DNA methylation inhibitor 5-azacytidine (5-AC) inhibited the reproduction of A. flavus and blocked AFs biosynthesis [25]. In our previous study, the 5-AC was further used to successfully isolate an A. flavus mutant [25], which displayed impaired abilities of AF biosynthesis and fungal development. Here, we performed a GC-MS analysis to reveal the metabolic changes in this mutant. A total number of 1181 volatiles were identified in the NT mutant and wild-type strain, among which 490 emitted compounds were identified in these two strains, and 332 intracellular compounds were found in these two strains intracellularly. Intriguingly, the decreasing volatile compounds were observed in the NT mutant compared to the wild-type strain, which could result in such physiological changes as reducing AF production and asporulate phenotype in the NT mutant.
In filamentous fungi, DNA methyltransferases inhibitors like 5-azacytidine have been used to identify novel metabolites re-expressing the methylation-silenced SM genes [20,24]. Although the existence of DNA methylation in Aspergillus species is controversial, we have previously demonstrated that treatment with 5-AC, that could inhibit the activity of DNA methyltransferase, resulted in inhibition of AF biosynthesis and developmental defects in A. flavus [25], which is consistent with the phenotypes defect of the disruption function of methyltransferase DmtA in A. flavus [28]. Intriguingly, here, we found that NT mutant produced decreasing volatile compounds compared to the wild-type strain. Our former study showed that the expression level of veA, a core component of the velvet complex, VeA/LaeA/VelB, was significantly downregulated in the NT mutant at a late growing stage. VeA is important for the regulation of fungal development and SMs with LaeA in Aspergillus, which is necessary for the activation of LaeA [29]. In Aspergillus, LaeA, as a global regulator of many secondary metabolisms also known as a methyltransferase, is involved in chromatin remodeling at the site of SMs gene clusters [9]. It is reasonable to see that, as the transcriptional expression level of veA was significantly decreased in the NT mutant, which is critical for the activation of the global regulator of SMs, LaeA, the volatile compounds produced by NT mutant are markedly reduced compared to the wild type.
Aflatoxin is the most important SMs produced by A. flavus which has caused severe food safety issues. Here, we found that the NT mutant was a defect in aflatoxins biosynthesis. Former studies have demonstrated that specific nutrients, like vitamins, amino acids and fatty acids are necessary for AFs formation [30]. In the AF biosynthesis pathway, it starts from acetyl-CoA, which is synthesized via oxidative decarboxylation of pyruvate, β-oxidation of fatty acids, or from ketogenic amino acids in primary metabolism [31]. In this study, we found that most of the fatty acid-derived volatiles and some of the amino acid-derived volatiles were significantly decreased in the NT mutant compared to the wild-type strain, which could explain the decreasing level of AFs in the mutant. Intriguingly, here, we found that AFB2 was detected in the hyphae of the NT strain, indicating that the secretion of AF might be impaired in the mutant. As previously reported, A. flavus preferentially invades oil-rich crops, which, as well as the Aspergilli, contain high levels of the unsaturated fatty acids, that are substrates for oxygenases [9]. These fatty acids, on one hand, are important synthesis intermediates for AF, they affect the balance between sclerotia formation and conidiation in A. flavus as well. The oxidative fatty acid, like oxylipins, encoded by ppo genes, are the signaling molecules that are involved in the regulation of developmental and SMs in A. nidulans [32]. Here, we found that the NT mutant failed to produce conidia, which could potentially be affected by the decreasing levels of fatty acid-derived compounds in the NT mutant.
Oxidative stress has been shown to be related to AFs biosynthesis and fungal development [33,34]. In Aspergillus or other filamentous fungi, the production of reactive oxygen species (ROS) and conidial germination/conidiogenesis are intimately correlated events. The mutant, which is defect in conidiation, was found to accumulate high levels of two antioxidants, including Phenol,2,2′-methylenebis[6-(1,1-dimethylethyl)-4-methyl-] and Phenol,4,4′-methylenebis[2,6-bis(1,1-dimethylethyl)-]. Our previous study demonstrated that the mutant induced by 5-AC, which showed a significant reduction in ROS accumulation, was sensitive to a strong oxidant, hydrogen peroxide. In A. nidulans, conidial pigmentation, mainly composed of naphthopyrones, is important for eliminating ROS for the protection of conidia against oxidative damage [35]. These data indicate that conidiogenesis might have a close relationship with oxidative stress.

4. Conclusions

In conclusion, we demonstrated that the NT mutant induced by the DNA methylation inhibitor, 5-azacytidine, was a severe defect in AF biosynthesis and fungal development. Additionally, we revealed that the NT mutant produced less volatile compounds, among which were the fatty acid-derived volatiles, which are important precursors for AF biosynthesis. This study may provide valuable information for controlling the contamination of mycotoxins produced by this pathogen.

5. Materials and Methods

5.1. Strains and Culture Conditions

Aspergillus flavus 3.4409 (A133) strain, which produces aflatoxin B1, B2, G1, and G2, was bought from the Culture Collection Center of the Chinese Academy of Sciences (Beijing, China). The NT strain is a non-aflatoxigenic mutant induced by 5-azacytidine [25]. For colony analysis, the A. flavus strains were cultured on PDA agar (BD DifcoTM, Franklin Lakes, NJ, USA) at 29 °C. The Czapek liquid medium was used for the extraction of aflatoxins (AFs) and volatile metabolites. To detect the conidiophore formation, after grown on PDA medium for two days, A. flavus colony was observed under a light microscope.

5.2. Aflatoxins Analysis

For aflatoxins and volatile metabolites extraction, a total amount of 107 A. flavus conidia was cultured into 50 mL of Czapek liquid medium at 29 °C for 5 days. After a five-day inoculation, the mycelium was harvested for the in vivo analysis of aflatoxins and volatile metabolites. Chloroform was used to extract AFs production according to the previously described method [28]. Thin Layer Chromatography (TLC) was performed to analyze AFs production.

5.3. GC-MS Volatile Analysis

An Agilent Model 7890A gas chromatograph (Agilent Technologies, Palo Alto, CA, USA) coupled to an Agilent 5975C Mass Selective Detector was used to detect the Aspergillus volatile metabolites. After being held at 160 °C for 3 min, the initial column temperature ramp at 10 °C min−1 to 225 °C for 3.5 min, then ramp 3 °C min−1 to the final temperature of 280 °C. The total run time was 70 min. The transfer line between the GC and MS systems was held at 290 °C. The ion source temperature was 230 °C.

5.4. Data Analysis

The resulted masses data were further analyzed by XCMS online program (https://xcmsonline.scripps.edu/) according to the published protocol [36]. The chemical structure or the chemical structural formula was analyzed in the PubChem database (https://pubchem.ncbi.nlm.nih.gov). Clustering and heat maps of the identified volatile metabolites after being analyzed with XCMS were drawn with the heatmap function in R (www.r-project.org). The top 30 volatile metabolites produced by A. flavus in vitro and the top 20 volatile metabolites produced by A. flavus in vivo that differed according to the t-test were further selected for analysis of heat map.

5.5. Statistical Analysis

In this study, the GraphPad Prism 6 was used for the analysis of statistics and significance. Student’s t-test was used for comparison of two different groups. One-way ANOVA multiple comparisons test was performed for significance analysis of multiple comparisons.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/1/57/s1, Figure S1: Clustering and heat map visualization of the identified volatile metabolites in WT and NT strains both in vitro and in vivo. (A) Clustering and heat map visualization of the identified volatile metabolites in WT (A133 strain) and NT strains in vitro. (B) Clustering and heat map visualization of the identified volatile metabolites in WT and NT strains in vivo, Table S1: Volatile metabolites assessed via GC-MS analysis in WT (A133 strain) and NT strains in vitro were analyzed with the Agilent ChemStation software, Table S2: Volatile metabolites assessed via GC-MS analysis in WT and NT strains in vivo were analyzed with the Agilent ChemStation software.

Author Contributions

K.Y., J.T., and S.W. designed the experiments and wrote the manuscript. F.S., K.Y., and Q.G. performed all the experiments. X.W., X.G., X.H., W.Z., and H.L. performed a few experiments and data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the grants of the National Natural Science Foundation of China (31900036, 31972171, 31671944), the Program of Natural Science Foundation of the Jiangsu Higher Education Institutions of China (19KJB180016), the Program of Natural Science Foundation of Jiangsu Normal University (18XLRX029), Natural Science Foundation of Jiangsu Province (BK20190994), Postgraduate Research & Practice Innovation Program of Jiangsu Province, Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Acknowledgments

The authors would like to thank Jiebo Chen for GCMS detection and Bo Yuan for XCMS guidance.

Conflicts of Interest

The authors report no conflicts of interest. The authors alone are responsible for the content and writing of the paper.

References

  1. Yang, K.; Shadkchan, Y.; Tannous, J.; Landero Figueroa, J.A.; Wiemann, P.; Osherov, N.; Wang, S.; Keller, N.P. Contribution of atpase copper transporters in animal but not plant virulence of the crossover pathogen aspergillus flavus. Virulence 2018, 9, 1273–1286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Chakrabarti, A.; Singh, R. The emerging epidemiology of mould infections in developing countries. Curr. Opin. Infect. Dis. 2011, 24, 521–526. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, K.; Liu, Y.; Wang, S.; Wu, L.; Xie, R.; Lan, H.; Fasoyin, O.E.; Wang, Y.; Wang, S. Cyclase-associated protein cap with multiple domains contributes to mycotoxin biosynthesis and fungal virulence in aspergillus flavus. J. Agric. Food Chem. 2019, 67, 4200–4213. [Google Scholar] [CrossRef]
  4. Rushing, B.R.; Selim, M.I. Aflatoxin b1: A review on metabolism, toxicity, occurrence in food, occupational exposure, and detoxification methods. Food Chem. Toxicol. 2019, 124, 81–100. [Google Scholar] [CrossRef] [PubMed]
  5. Hedayati, M.T.; Pasqualotto, A.C.; Warn, P.A.; Bowyer, P.; Denning, D.W. Aspergillus flavus: Human pathogen, allergen and mycotoxin producer. Microbiology 2007, 153, 1677–1692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ruadrew, S.; Craft, J.; Aidoo, K. Occurrence of toxigenic aspergillus spp. And aflatoxins in selected food commodities of asian origin sourced in the west of scotland. Food Chem. Toxicol. 2013, 55, 653–658. [Google Scholar] [CrossRef]
  7. Tian, J.; Gan, Y.; Pan, C.; Zhang, M.; Wang, X.; Tang, X.; Peng, X. Nerol-induced apoptosis associated with the generation of ros and Ca2+ overload in saprotrophic fungus aspergillus flavus. Appl. Microbiol. Biotechnol. 2018, 102, 6659–6672. [Google Scholar] [CrossRef]
  8. Tian, J.; Wang, Y.; Lu, Z.; Sun, C.; Zhang, M.; Zhu, A.; Peng, X. Perillaldehyde, a promising antifungal agent used in food preservation, triggers apoptosis through a metacaspase-dependent pathway in aspergillus flavus. J. Agric. Food Chem. 2016, 64, 7404–7413. [Google Scholar] [CrossRef]
  9. Amaike, S.; Keller, N.P. Aspergillus flavus. Annu. Rev. Phytopathol. 2011, 49, 107–133. [Google Scholar] [CrossRef]
  10. Chalivendra, S.C.; DeRobertis, C.; Chang, P.K.; Damann, K.E. Cyclopiazonic acid is a pathogenicity factor for aspergillus flavus and a promising target for screening germplasm for ear rot resistance. Mol. Plant-Microbe Interact. MPMI 2017, 30, 361–373. [Google Scholar] [CrossRef] [Green Version]
  11. Chang, P.K.; Horn, B.W.; Dorner, J.W. Clustered genes involved in cyclopiazonic acid production are next to the aflatoxin biosynthesis gene cluster in aspergillus flavus. Fungal Genet. Biol. 2009, 46, 176–182. [Google Scholar] [CrossRef]
  12. Gilbert, M.K.; Mack, B.M.; Moore, G.G.; Downey, D.L.; Lebar, M.D.; Joardar, V.; Losada, L.; Yu, J.; Nierman, W.C.; Bhatnagar, D. Whole genome comparison of aspergillus flavus l-morphotype strain nrrl 3357 (type) and s-morphotype strain af70. PLoS ONE 2018, 13, e0199169. [Google Scholar] [CrossRef] [PubMed]
  13. Payne, G.A.; Nierman, W.C.; Wortman, J.R.; Pritchard, B.L.; Brown, D.; Dean, R.A.; Bhatnagar, D.; Cleveland, T.E.; Machida, M.; Yu, J. Whole genome comparison ofaspergillus flavusanda. Oryzae. Med. Mycol. 2006, 44, 9–11. [Google Scholar] [CrossRef] [PubMed]
  14. Amare, M.G.; Keller, N.P. Molecular mechanisms of aspergillus flavus secondary metabolism and development. Fungal Genet. Biol. 2014, 66, 11–18. [Google Scholar] [CrossRef] [PubMed]
  15. Georgianna, D.R.; Fedorova, N.D.; Burroughs, J.L.; Dolezal, A.L.; Bok, J.W.; Horowitz-Brown, S.; Woloshuk, C.P.; Yu, J.; Keller, N.P.; Payne, G.A. Beyond aflatoxin: Four distinct expression patterns and functional roles associated with aspergillus flavus secondary metabolism gene clusters. Mol. Plant Pathol. 2010, 11, 213–226. [Google Scholar] [CrossRef] [PubMed]
  16. Khaldi, N.; Collemare, J.; Lebrun, M.H.; Wolfe, K.H. Evidence for horizontal transfer of a secondary metabolite gene cluster between fungi. Genome Biol. 2008, 9, R18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Roze, L.V.; Chanda, A.; Laivenieks, M.; Beaudry, R.M.; Artymovich, K.A.; Koptina, A.V.; Awad, D.W.; Valeeva, D.; Jones, A.D.; Linz, J.E. Volatile profiling reveals intracellular metabolic changes in aspergillus parasiticus: Vea regulates branched chain amino acid and ethanol metabolism. BMC Biochem. 2010, 11, 33. [Google Scholar] [CrossRef] [Green Version]
  18. Osbourn, A. Secondary metabolic gene clusters: Evolutionary toolkits for chemical innovation. Trends Genet. 2010, 26, 449–457. [Google Scholar] [CrossRef]
  19. Pfannenstiel, B.T.; Keller, N.P. On top of biosynthetic gene clusters: How epigenetic machinery influences secondary metabolism in fungi. Biotechnol. Adv. 2019, 37, 107345. [Google Scholar] [CrossRef]
  20. Sun, J.; Awakawa, T.; Noguchi, H.; Abe, I. Induced production of mycotoxins in an endophytic fungus from the medicinal plant Datura stramonium L. Bioorganic Med. Chem. Lett. 2012, 22, 6397–6400. [Google Scholar] [CrossRef] [PubMed]
  21. Shwab, E.K.; Bok, J.W.; Tribus, M.; Galehr, J.; Graessle, S.; Keller, N.P. Histone deacetylase activity regulates chemical diversity in aspergillus. Eukaryot. Cell 2007, 6, 1656–1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Sun, Y.; Sun, Y.; Yue, S.; Wang, Y.; Lu, F. Histone deacetylase inhibitors in cancer therapy. Curr. Top. Med. Chem. 2018, 18, 2420–2428. [Google Scholar] [CrossRef] [PubMed]
  23. Shigemoto, R.; Matsumoto, T.; Masuo, S.; Takaya, N. 5-methylmellein is a novel inhibitor of fungal sirtuin and modulates fungal secondary metabolite production. J. Gen. Appl. Microbiol. 2018, 64, 240–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Savickiene, J.; Treigyte, G.; Jazdauskaite, A.; Borutinskaite, V.V.; Navakauskiene, R. DNA methyltransferase inhibitor rg108 and histone deacetylase inhibitors cooperate to enhance nb4 cell differentiation and e-cadherin re-expression by chromatin remodelling. Cell Biol. Int. 2012, 36, 1067–1078. [Google Scholar] [CrossRef]
  25. Yang, K.; Zhuang, Z.; Zhang, F.; Song, F.; Zhong, H.; Ran, F.; Yu, S.; Xu, G.; Lan, F.; Wang, S. Inhibition of aflatoxin metabolism and growth of aspergillus flavus in liquid culture by a DNA methylation inhibitor. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 2015, 32, 554–563. [Google Scholar] [CrossRef]
  26. Carvajal-Moreno, M.; Vargas-Ortiz, M.; Hernandez-Camarillo, E.; Ruiz-Velasco, S.; Rojo-Callejas, F. Presence of unreported carcinogens, aflatoxins and their hydroxylated metabolites, in industrialized oaxaca cheese from mexico city. Food Chem. Toxicol. 2019, 124, 128–138. [Google Scholar] [CrossRef]
  27. Qu, S.; Yang, K.; Chen, L.; Liu, M.; Geng, Q.; He, X.; Li, Y.; Liu, Y.; Tian, J. Cinnamaldehyde, a promising natural preservative against aspergillus flavus. Front. Microbiol. 2019, 10, 2895. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, K.; Liang, L.; Ran, F.; Liu, Y.; Li, Z.; Lan, H.; Gao, P.; Zhuang, Z.; Zhang, F.; Nie, X.; et al. The dmta methyltransferase contributes to aspergillus flavus conidiation, sclerotial production, aflatoxin biosynthesis and virulence. Sci. Rep. 2016, 6, 23259. [Google Scholar] [CrossRef] [Green Version]
  29. Chang, P.K.; Scharfenstein, L.L.; Li, P.; Ehrlich, K.C. Aspergillus flavus velb acts distinctly from vea in conidiation and may coordinate with flug to modulate sclerotial production. Fungal Genet. Biol. 2013, 58–59, 71–79. [Google Scholar] [CrossRef]
  30. Han, X.; Qiu, M.; Wang, B.; Yin, W.B.; Nie, X.; Qin, Q.; Ren, S.; Yang, K.; Zhang, F.; Zhuang, Z.; et al. Functional analysis of the nitrogen metabolite repression regulator gene nmra in aspergillus flavus. Front. Microbiol. 2016, 7, 1794. [Google Scholar] [CrossRef] [Green Version]
  31. Yabe, K.; Nakajima, H. Enzyme reactions and genes in aflatoxin biosynthesis. Appl. Microbiol. Biotechnol. 2004, 64, 745–755. [Google Scholar] [CrossRef] [PubMed]
  32. Tsitsigiannis, D.I.; Keller, N.P. Oxylipins act as determinants of natural product biosynthesis and seed colonization in aspergillus nidulans. Mol. Microbiol. 2006, 59, 882–892. [Google Scholar] [CrossRef] [PubMed]
  33. Reverberi, M.; Ricelli, A.; Zjalic, S.; Fabbri, A.A.; Fanelli, C. Natural functions of mycotoxins and control of their biosynthesis in fungi. Appl. Microbiol. Biotechnol. 2010, 87, 899–911. [Google Scholar] [CrossRef] [PubMed]
  34. Narasaiah, K.V.; Sashidhar, R.B.; Subramanyam, C. Biochemical analysis of oxidative stress in the production of aflatoxin and its precursor intermediates. Mycopathologia 2006, 162, 179–189. [Google Scholar] [CrossRef] [PubMed]
  35. Jahn, B.; Boukhallouk, F.; Lotz, J.; Langfelder, K.; Wanner, G.; Brakhage, A.A. Interaction of human phagocytes with pigmentless aspergillus conidia. Infect. Immun. 2000, 68, 3736–3739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Gowda, H.; Ivanisevic, J.; Johnson, C.H.; Kurczy, M.E.; Benton, H.P.; Rinehart, D.; Nguyen, T.; Ray, J.; Kuehl, J.; Arevalo, B.; et al. Interactive xcms online: Simplifying advanced metabolomic data processing and subsequent statistical analyses. Anal. Chem. 2014, 86, 6931–6939. [Google Scholar] [CrossRef]
Figure 1. Colony morphology and aflatoxins analysis of A. flavus A133 and NT strains. (A) Colony morphology of A. flavus A133 and NT strains grown on PDA medium for three days. (B) Conidiophore formation was observed under a light microscope after 24 h incubation on PDA agar medium. (C) Colony morphology of A. flavus A133 and NT strains grown in Czapek liquid medium. The bottom panel was the liquid extracts of A. flavus A133 and NT strains which displayed a different color. (D) TLC detection of aflatoxins extracted from A. flavus A133 and NT strains in vivo and in vitro.
Figure 1. Colony morphology and aflatoxins analysis of A. flavus A133 and NT strains. (A) Colony morphology of A. flavus A133 and NT strains grown on PDA medium for three days. (B) Conidiophore formation was observed under a light microscope after 24 h incubation on PDA agar medium. (C) Colony morphology of A. flavus A133 and NT strains grown in Czapek liquid medium. The bottom panel was the liquid extracts of A. flavus A133 and NT strains which displayed a different color. (D) TLC detection of aflatoxins extracted from A. flavus A133 and NT strains in vivo and in vitro.
Toxins 12 00057 g001
Figure 2. GC-MS spectra of volatiles extracted from A. flavus A133 and NT strains. (A) GC-MS spectra of volatiles emitted from the A. flavus A133 and NT strains (in vitro). (B) GC-MS spectra of volatiles extracted from the mycelium of A. flavus A133 and NT strains (in vivo).
Figure 2. GC-MS spectra of volatiles extracted from A. flavus A133 and NT strains. (A) GC-MS spectra of volatiles emitted from the A. flavus A133 and NT strains (in vitro). (B) GC-MS spectra of volatiles extracted from the mycelium of A. flavus A133 and NT strains (in vivo).
Toxins 12 00057 g002
Figure 3. Principal components analysis (PCA) and cloud plot of the identified metabolic compounds by XCMS online program (https://xcmsonline.scripps.edu/). (A) Principal components analysis (PCA) of the identified in vitro volatile metabolites in WT and NT strains by XCMS online program. The PCA is calculated using the feature intensities from all samples. The colors (red/green) are assigned based on the sample class. (B) PCA analysis of the identified in vivo volatile metabolites in WT and NT strains by the XCMS online program. (C) Cloud plot of the identified in vitro volatile metabolites in WT and NT strains. (D) Cloud plot of the identified in vivo volatile metabolites in WT and NT strains. Only features that are dysregulated (p-value ≤ 0.01, fold change ≥ 1.5) are displayed. Upregulated features are shown in green, downregulated features in red. The size of each bubble corresponds to the log fold change of the feature. The shade of the bubbles corresponds to the magnitude of the p-value (the darker the color, the smaller the p-value).
Figure 3. Principal components analysis (PCA) and cloud plot of the identified metabolic compounds by XCMS online program (https://xcmsonline.scripps.edu/). (A) Principal components analysis (PCA) of the identified in vitro volatile metabolites in WT and NT strains by XCMS online program. The PCA is calculated using the feature intensities from all samples. The colors (red/green) are assigned based on the sample class. (B) PCA analysis of the identified in vivo volatile metabolites in WT and NT strains by the XCMS online program. (C) Cloud plot of the identified in vitro volatile metabolites in WT and NT strains. (D) Cloud plot of the identified in vivo volatile metabolites in WT and NT strains. Only features that are dysregulated (p-value ≤ 0.01, fold change ≥ 1.5) are displayed. Upregulated features are shown in green, downregulated features in red. The size of each bubble corresponds to the log fold change of the feature. The shade of the bubbles corresponds to the magnitude of the p-value (the darker the color, the smaller the p-value).
Toxins 12 00057 g003
Figure 4. Differential volatile metabolites in the A133 (WT) strain compared to the NT strain. (A) Volcano plots of the differential in vitro volatile metabolites in the A133 strain compared to the NT strain. The gray spots inside the blue dash lines indicate the non-significant volatile metabolites. Volatile metabolites with p-value > 0.05 (−log10 (p-value) < 1.3) or log2 (Fold Change) < 1 are considered as non-significantly different volatile metabolites. (B) Volcano plots of the differential in vivo volatile metabolites in the A133 strain compared to the NT strain. The dashed line indicates p = 0.05. The gray spots inside the blue dash lines indicate the non-significant volatile metabolites. (C) The total specific in vitro volatile metabolites produced by A133 and NT strains compared to the in vivo volatile metabolites produced by these two strains.
Figure 4. Differential volatile metabolites in the A133 (WT) strain compared to the NT strain. (A) Volcano plots of the differential in vitro volatile metabolites in the A133 strain compared to the NT strain. The gray spots inside the blue dash lines indicate the non-significant volatile metabolites. Volatile metabolites with p-value > 0.05 (−log10 (p-value) < 1.3) or log2 (Fold Change) < 1 are considered as non-significantly different volatile metabolites. (B) Volcano plots of the differential in vivo volatile metabolites in the A133 strain compared to the NT strain. The dashed line indicates p = 0.05. The gray spots inside the blue dash lines indicate the non-significant volatile metabolites. (C) The total specific in vitro volatile metabolites produced by A133 and NT strains compared to the in vivo volatile metabolites produced by these two strains.
Toxins 12 00057 g004
Figure 5. Clustering and heat map visualization of volatile metabolites in the A133 (WT) strain and NT strain in vitro and vivo. (A) Clustering/heat map visualization of the top 30 volatile metabolites produced by the A133 strain and NT strain in vitro differing according to t-test. (B) Clustering/heat map visualization of the top 20 volatile metabolites produced by the A133 strain and NT strain in vivo differing according to t-test.
Figure 5. Clustering and heat map visualization of volatile metabolites in the A133 (WT) strain and NT strain in vitro and vivo. (A) Clustering/heat map visualization of the top 30 volatile metabolites produced by the A133 strain and NT strain in vitro differing according to t-test. (B) Clustering/heat map visualization of the top 20 volatile metabolites produced by the A133 strain and NT strain in vivo differing according to t-test.
Toxins 12 00057 g005
Table 1. The top 30 volatile metabolites produced by the A133 strain and NT strain in vitro with the XCMS program. Results are derived from 3 replicates. The metabolites produced by NT mutants were compared to the A133 strain. Log2fold = log2((NT intensity)/(WT intensity)).
Table 1. The top 30 volatile metabolites produced by the A133 strain and NT strain in vitro with the XCMS program. Results are derived from 3 replicates. The metabolites produced by NT mutants were compared to the A133 strain. Log2fold = log2((NT intensity)/(WT intensity)).
NameVolatile CompoundRetention Time (min)Quantification Ion (m/z)Intensitylog2fold
A133NT
Fatty acid-derived volatiles
M256T43n-Hexadecanoic acid42.53256173,803 ± 589942 ± 22−11.48
M284T43Hexadecanoicacid,ethylester43.3128410,182 ± 273162 ± 57−5.97
M59T44Normeperidinic acid44.3320547,057 ± 1980711 ± 114−6.05
M280T469,12-Octadecadienoicacid(Z,Z)-46.47280168,927 ± 10,05748 ± 48−10.95
M264T47_26-Octadecenoic acid46.6026449,278 ± 16970 ± 0−9.77
M308T479,12-Octadecadienoicacid,ethylester47.0830814,919 ± 3360 ± 0−8.05
M284T47n-Propyl9-octadecenoate47.1428418,212 ± 16210 ± 0−8.34
M312T48Heptadecanoicacid,15-methyl-,ethylester47.873126804 ± 33030 ± 15−7.12
M253T606,9-hexadecadienoic acid60.2825342,150 ± 672291 ± 50−7.18
M143T66Ergosterol66.43143407,309 ± 17,0883840 ± 72−6.73
Phenylpropanoid derivatives/benzenoids
M159T29N-Methyl-N-methoxy-5,6,7,8-tetrahydro-1-naphtamide29.3315973,675 ± 29201314 ± 197−5.81
M207T30Phenol,2,4-bis(1,1-dimethylethyl)-30.16207636 ± 37237,011 ± 43545.82
M201T344′-Azidobenzo[1′,2′-b]-1,4-diazabicyclo[2.2.2]octene33.9120110,246 ± 44664 ± 64−6.66
M254T39Phenol,3,5-dimethoxy-38.5125485 ± 153599 ± 1015.4
M182T39Flopropione38.63182340,108 ± 32,3662423 ± 111−7.13
M59T39Phthalicacid,isobutylnonylester39.405932,073 ± 2209818 ± 191−5.29
M261T49Tiaprofenic acid49.0126037 ± 3710,811 ± 6277.18
M371T542,4-Difluorophenol53.6437118,882 ± 109860 ± 60-7.6
M353T55Phenol,4,4′-methylenebis[2,6-bis(1,1-dimethylethyl)]-55.184241710 ± 25016,636 ± 5813.28
M130T571H-Indole,4-methyl-57.4513087,806 ± 14,2511825 ± 109−5.59
M251T62Propiophenone,2′-(trimethylsiloxy)-62.0022240,280 ± 13501112 ± 147−5.18
Triterpene
M69T602,6,10,14,18,22-Tetracosahexaene,2,6,10,15,19,23-hexamethyl-,(all-E)-59.85410694,833 ± 19,8218465 ± 313−6.36
Amino acids derived volatiles
M182T34L-Tyrosine33.6618268 ± 347028 ± 3156.34
M240T366-Methyl-3,5-heptadien-2-one35.8112415 ± 154429 ± 2006.4
M269T411-Hexadecanamine,N,N-dimethyl-41.182690 ± 039,265 ± 1459.44
M226T46N-Methyl-beta-carboline-3-carboxamide46.12250 ± 04410 ± 7946.29
M265T479-Octadecenamide,(Z)-47.2526518,661 ± 2130 ± 0−8.37
M179T512-phenyl-N-(5-propan-2-yl-1,3-thiazol-2-yl)acetamide50.522601332 ± 17389,236 ± 3706.07
M120T544-Morpholineethanol54.4013157,950 ± 3565562 ± 108−6.69
M284T551-(2-Thiazolylazo)-2-naphthol54.7425533,915 ± 1238170 ± 25−7.64
Table 2. The top 20 volatile metabolites produced by the A133 strain and NT strain in vivo with the XCMS program. Results are derived from 3 replicates. The metabolites produced by the NT mutants were compared to the A133 strain. Log2fold = log2((NT intensity)/(WT intensity)).
Table 2. The top 20 volatile metabolites produced by the A133 strain and NT strain in vivo with the XCMS program. Results are derived from 3 replicates. The metabolites produced by the NT mutants were compared to the A133 strain. Log2fold = log2((NT intensity)/(WT intensity)).
NameVolatile CompoundRetention Time (min)Quantification Ion (m/z)Intensitylog2fold
A133NT
Fatty acid-derived volatiles
M119T27Trimethylammonioacetate27.35117.1532,031 ± 2461471 ± 71−4.44
M100T54Hexadecanoicacid,2-hydroxy-1-(hydroxymethyl)ethylester54.40330.536,178 ± 7533406 ± 1018−3.41
M284T551,2-Benzenedicarboxylicacid,mono(2-ethylhexyl)ester54.7454.7453,296 ± 2260378 ± 61−7.14
M158T66Ergosterol66.4414316,097 ± 6460897 ± 897−3.47
Phenylpropanoid derivatives/benzenoids
M355T182,6-Dichloroindophenol17.59268.133,856 ± 3102579 ± 469−3.71
M167T35_1Benzo[b]tetrahydrofuran-3-one,5,6-dihydroxy-34.79166.0271110 ± 35111,478 ± 74176.65
M166T35_12,6-toluenediamine35.23122.174043 ± 615211,094 ± 56065.71
M209T40Isoelemicin40.11208.2513,703 ± 254353,432 ± 42584.69
Saccharide derived volatiles
M195T41D-Glucoside41.2194.1834,031 ± 1851,301,860 ± 92105.26
Amino acids derived volatiles
M71T30Caulophylline29.66204.27112,547 ± 86812,794 ± 455−3.14
M181T36Theobromine36.27180.1611,780 ± 637382 ± 66−4.95
M206T38N-Acetyl-D-quinovosamine37.95205.21170 ± 647878 ± 6345.54
M195T39N,N′-Diacetyl-2-nitro-p-phenylenediamine38.6237.21574,864 ± 10,2075272 ± 227−6.77
M164T392,2′-(2-Hydroxy-2-nitrosohydrazinylidene)bis-ethanamine38.61163.1754,132 ± 22061027 ± 49−5.72
M224T414-(Methylnitrosamino)-1-(1-oxido-3-pyridinyl)-1-butanone40.90223.23293 ± 387571 ± 1964.69
M196T419-Aminoacridine41.19194.234102 ± 160145,064 ± 24685.14
M168T43_22,8-Dihydroxyadenine42.71167.12455 ± 17411,735 ± 8054.69
M194T44Phenylacetylglycine43.50194.11457 ± 10318,559 ± 3713.67
M238T44N-Benzylphthalimide44.31237.25469 ± 20132,530 ± 17026.12
M241T55Tetramisole hydrochloride54.74240.759247 ± 476243 ± 101−5.25

Share and Cite

MDPI and ACS Style

Song, F.; Geng, Q.; Wang, X.; Gao, X.; He, X.; Zhao, W.; Lan, H.; Tian, J.; Yang, K.; Wang, S. Gas Chromatography–Mass Spectrometry Profiling of Volatile Compounds Reveals Metabolic Changes in a Non-Aflatoxigenic Aspergillus flavus Induced by 5-Azacytidine. Toxins 2020, 12, 57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010057

AMA Style

Song F, Geng Q, Wang X, Gao X, He X, Zhao W, Lan H, Tian J, Yang K, Wang S. Gas Chromatography–Mass Spectrometry Profiling of Volatile Compounds Reveals Metabolic Changes in a Non-Aflatoxigenic Aspergillus flavus Induced by 5-Azacytidine. Toxins. 2020; 12(1):57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010057

Chicago/Turabian Style

Song, Fengqin, Qingru Geng, Xuewei Wang, Xiaoqing Gao, Xiaona He, Wei Zhao, Huahui Lan, Jun Tian, Kunlong Yang, and Shihua Wang. 2020. "Gas Chromatography–Mass Spectrometry Profiling of Volatile Compounds Reveals Metabolic Changes in a Non-Aflatoxigenic Aspergillus flavus Induced by 5-Azacytidine" Toxins 12, no. 1: 57. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12010057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop