Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 05 April 2019
Sec. Fungi and Their Interactions
This article is part of the Research Topic The Chemistry, Analysis, and Toxic Effects of Emerging Mycotoxins View all 5 articles

Secondary Metabolite Dereplication and Phylogenetic Analysis Identify Various Emerging Mycotoxins and Reveal the High Intra-Species Diversity in Aspergillus flavus

  • 1Center of Excellence in Mycotoxicology and Public Health, Faculty of Pharmaceutical Sciences, Ghent University, Ghent, Belgium
  • 2Department of Pharmacy, Faculty of Medicine, University of Prishtina, Prishtina, Kosovo
  • 3Southern Regional Research Center, Agricultural Research Service, United States Department of Agriculture, New Orleans, LA, United States
  • 4Department of Analytical Chemistry, Faculty of Chemistry, Regional Campus of International Excellence “Campus Mare-Nostrum”, University of Murcia, Murcia, Spain

Aspergillus flavus is one of the most important mycotoxigenic species from the genus Aspergillus, due to its ability to synthesize the potent hepatocarcinogen, aflatoxin B1. Moreover, this fungus is capable of producing several other toxic metabolites from the class of indole-tetramates, non-ribosomal peptides, and indole-diterpenoids. Populations of A. flavus are characterized by considerable diversity in terms of morphological, functional and genetic features. Although for many years A. flavus was considered an asexual fungus, researchers have shown evidence that at best these fungi can exhibit a predominantly asexual existence. We now know that A. flavus contains functional genes for mating, uncovering sexuality as potential contributor for its diversification. Based on our results, we reconfirm that A. flavus is a predominant producer of B-type aflatoxins. Moreover, this fungus can decisively produce AFM1 and AFM2. We did not observe any clear relationship between mating-type genes and particular class of metabolites, probably other parameters such as sexual/asexual ratio should be investigated. A dynamic secondary metabolism was found also in strains intended to be used as biocontrol agents. In addition we succeeded to provide mass spectrometry fragmentation spectra for the most important classes of A. flavus metabolites, which will serve as identification cards for future studies. Both, metabolic and phylogenetic analysis proved a high intra-species diversity for A. flavus. These findings contribute to our understanding about the diversity of Aspergillus section Flavi species, raising the necessity for polyphasic approaches (morphological, metabolic, genetic, etc.) when dealing with this type of complex group of species.

Introduction

Among the casual agents of food and feed contamination, fungi from the genus Aspergillus play a crucial role. Within this genus, there are species that have the ability to produce a broad assortment of secondary metabolites, including mycotoxins. One very important group of fungi from Aspergillus section Flavi is known to produce the most potent carcinogenic compound of the mycotoxins, aflatoxin. In this connection, Aspergillus flavus is one of the most studied mycotoxigenic fungi (Bennett, 2010). Polyketide-derived mycotoxins like aflatoxin B1 have also been shown to possess additional hepatotoxic and immunosuppressive properties (Turner et al., 2003). Aflatoxin levels exceeding 0.025–15 μg/kg on food or feed commodities are strictly prohibited in the European Union (Commission regulation (EC) No 1881/2006). Moreover, A. flavus has the potential to act as a human pathogen, being the second most frequent cause, after A. fumigatus, of invasive and localized aspergillosis in immunocompromised patients (Hedayati et al., 2007).

Beyond aflatoxins, many other toxic secondary metabolites have been detected that contribute to the overall toxicity of A. flavus, including: other polyketides such as asparasones (Cary et al., 2014) and aflavarins (Cary et al., 2015a); non-ribosomal peptides such as piperazines (Forseth et al., 2013) and ditryptophenaline (Springer et al., 1977); hybrid molecules such as cyclopiazonic acids (Luk et al., 1977; Uka et al., 2017) and leporins (Arroyo-Manzanares et al., 2015; Cary et al., 2015b); and various indole-diterpenes such as aflatrem (Gallagher and Wilson, 1979; Nicholson et al., 2009) and aflavinines (Gallagher et al., 1980). Alpha-cyclopiazonic acid (α-CPA) (Uka et al., 2017) is another important A. flavus mycotoxin, which has been associated with different liver, kidney and gastrointestinal complications in animal health (Cole, 1986; Burdock and Flamm, 2000). Aflatoxins and α-CPA very often co-contaminate food and feed commodities (Cole, 1986; Lee and Hagler, 1991; Urano et al., 1992; Heperkan et al., 2012; Zorzete et al., 2013; Ezekiel et al., 2016). A more complete list of key fungal metabolites identified in A. flavus, based on the literature, is shown in Supplementary Figure S1 and Supplementary Table S1.

In general, A. flavus and related species within the ascomycete group of fungi are considered primarily asexual microorganisms (Debuchy and Turgeon, 2006). However, A. flavus exhibits high genetic diversity, which is reflected by its spatial distribution and by its large number of vegetative compatibility groups (VCGs), but also by diversity of its morphology (e.g., sclerotium size and number, conidium production) and physiology (e.g., secondary metabolite production) (Horn et al., 1996). Fungal individuals belonging to the same VCG are vegetatively compatible, so their hyphae have the ability to anastomose and possibly exchange genetic material via parasexual modes of recombination. Individuals from different VCGs are vegetatively incompatible, so hyphal fusion typically results in cell necrosis (Papa, 1973; Glass and Kaneko, 2003). Parasexuality occurs only among strains of the same VCG, which usually leads to isolation and homogeneity in metabolite patterns and morphology. This is the reason why individuals within a given VCG share similar morphological and physiological features (Grubisha and Cotty, 2009). The high level of diversity exhibited among A. flavus strains, even within localized populations, is difficult to explain by asexuality or parasexuality; therefore, sexuality and subsequent genetic recombination should be considered important for generating observed intra-specific diversity. Hence, research began in the 1990s aiming to investigate a possible sexual stage within the life cycles of filamentous ascomycetes (Metzenberg and Glass, 1990; Coppin et al., 1997). Moreover, it has been proposed that also horizontal gene transfer (HGT) has contributed significantly in the overall chemodiversity of Aspergillus species (Frisvad and Larsen, 2015). While parasexuality mainly occurs between strains of the same species, HGT is more likely to occur between phylogenetically distant species (Rank et al., 2011; Campbell et al., 2012). HGT of a secondary metabolite gene cluster can be a result of species occurring in the same habitat and sharing the same ecological challenges (Frisvad and Larsen, 2015). Specifically, through bioinformatics tools it was shown that A. flavus genome likely harbors more than 500 horizontally transferred genes, 41% of which were found to be engaged in physically linked gene clusters (Nguyen et al., 2015).

In 1998, evidence of cryptic sex was reported for an Australian A. flavus population, which suggested a potential risk for using non-aflatoxigenic strains of A. flavus as biopesticides (Geiser et al., 1998). Nearly a decade after, evidence of recombination was uncovered for a single A. parasiticus population (Carbone et al., 2007). Then it was determined that both A. flavus and A. parasiticus are heterothallic (self-infertile) fungi with individuals possessing either a MAT1-1 or MAT1-2 mating-type (MAT) gene (Ramirez-Prado et al., 2008), followed by additional evidence of recombination for a single population of A. flavus (Moore et al., 2009). MAT1-1 encodes for a conserved alpha (α)-domain protein, whereas MAT1-2 gene encodes for a conserved class of high mobility group (HMG)-domain protein (Turgeon and Yoder, 2000). Gene-disruption studies have shown that MAT genes are required for normal sexual development in heterothallic species (Debuchy and Turgeon, 2006). Moreover, is important to note that sexual reproduction is correlated with a significant increase expression of MAT genes and key genes of a pheromone-response MAP (mitogen-activated protein)-kinase signaling pathway involved in heterothallic outcrossing (Lengeler et al., 2000). Next, the sexual states for A. parasiticus (Horn et al., 2009a) and A. flavus (Horn et al., 2009b) were observed via experimental crosses. Finally, it was reported that the distribution of MAT genes in several global populations of A. flavus and A. parasiticus correlated with the amount of diversity within the population, as well as the distribution of chemotype profiles across a given population. As well, equal distributions (1:1 ratio) of MAT genes in populations of A. flavus and A. parasiticus were found to result in greater likelihood of uncovering evidence for the existence of sexual reproduction (Moore et al., 2013). These reports offered new insights into the biology of these agriculturally significant molds. In contrast to parasexuality, sexual reproduction happens between heterothallic individuals with opposite mating type (MAT1-1 or MAT1-2) that belong to different VCGs, and may differ in their morphological and functional features (Horn et al., 2009a,b). Consequently, it was understood that asexual reproduction fixes metabolic profiles in fungal populations whereas sexual reproduction creates new VCGs with diversity of genetic composition and metabolite chemotypes (Moore et al., 2013). Even inter-specific hybridization via laboratory crosses has been observed, which offers more evidence of how these fungi may be capable of diversifying, morphologically and physiologically, as well as their overall adaptive ability (Olarte et al., 2015). Apart from laboratory experiments, A. flavus sexuality was also demonstrated in experimental field studies (Horn et al., 2013), and potential male/female roles were recently suggested for conidia/sclerotia as well as for each mating type (Horn et al., 2016).

The discovery of sexual states in aflatoxigenic fungi is expected to have a great impact in agriculture, ecology, and food safety. This is because mycotoxigenic species that undergo sexual recombination exhibit increased adaptability and may create more difficulties for their control and monitoring. Specifically, this may interfere with current strategies for developing biopesticides (pre-harvest biological controls) through use of non-aflatoxigenic A. flavus as a replacement for chemical pesticides (Abbas et al., 2011; Probst et al., 2011; Moore et al., 2013). This is highly relevant because its use as biological control is largely based on the paradigm that A. flavus is reproductively stable (i.e., asexual) (Ehrlich and Cotty, 2004). The long-term fate of A. flavus biocontrol strains in fields has not yet been fully addressed (Ehrlich, 2014), and the possibility of inheriting aflatoxin-producing ability for progeny that arise through sexual recombination is a possibility that should not be neglected after long-term treatment (Olarte et al., 2011).

Most investigations reporting correlations between sclerotium size, mating-type and VCGs with mycotoxin production have involved the aflatoxins because of their toxicity (Novas and Cabral, 2002; Amani et al., 2012; Moore et al., 2013). However, although A. flavus has potential to produce many other toxic secondary metabolites, most of the toxic effects of these secondary metabolites have been ascribed to aflatoxins. There is the possibility that some of these toxic metabolites are being produced in higher quantities by non-aflatoxigenic strains, which further might complicate the application of biocontrol strains in pre-harvest decontamination. Hence, the main objective of this study is to evaluate any relationship between secondary metabolite patterns in A. flavus and mating-type on one hand, and sclerotium size on the other hand, including not only aflatoxins but a much broader repertoire of secondary metabolites produced by this fungus. Thus, we aim to get a more complete picture of the mosaic of the A. flavus metabolome as a function of different genetic and morphological features. To address this issue, a secondary metabolite screening program was performed on different biological strains of A. flavus deploying an identification methodology based on high-resolution mass spectrometry (HRMS). Indeed, screening of biological samples through full-scan HRMS can result in an unlimited number of putative metabolites that can be identified. However, the final confirmation of these metabolites is still seen as the major bottleneck in MS-based analytical approaches. MS/MS or fragmentation data allows a stronger absolute metabolite identification. In order to obtain these types of data (both accurate mass measurements and MS/MS) within a single analytical run, an UHPLC hybrid qTOF (quadrupole Time-of-Flight) mass spectrometer was used. In addition, complementary identification tools such as isotope similarity and elution order were applied. Where possible, the retention (RT) time and MS/MS spectra of putative compounds were matched with authentic reference standards, although standards are not available for all of the known or suspected metabolites.

Materials and Methods

Chemicals and Materials

Methanol (MeOH), LC–MS grade, was obtained from Biosolve (Valkenswaard, Netherlands), and HPLC-grade MeOH was obtained from VWR International (Zaventem, Belgium). Ethyl acetate (EtOAc), dichloromethane (DCM) and acetone [dimethyl ketone (DMK)] were purchased from Acros Organics (Geel, Belgium). Sigma-Aldrich (Bornem, Belgium) supplied ammonium formate (HCOONH4). Formic acid (HCOOH, Merck, Darmstadt, Germany) was used. Ultrapure H2O, used during these analyses was produced in-house with a Milli-Q Gradient System (Millipore, Brussels, Belgium). Ultrafree®-MC centrifugal filter units (0.22 μm) from Millipore (Bedford, MA, United States) were used for filtration and reconstitution of fungal extracts. Sigma-Aldrich was our source for fungal growth media, such as agar, corn steep solids, dextrose, peptone, sucrose, yeast extract, dipotassium hydrogen phosphate trihydrate (K2HPO4⋅3H2O), magnesium sulfate heptahydrate (MgSO4⋅7H2O), and iron(II) sulfate heptahydrate (FeSO4⋅7H2O). Other compounds necessary for this study included Triton X-100, potassium chloride (KCl), and sodium nitrate (NaNO3), which were obtained from Merck. Standards of AFB1, AFB2, AFG1, AFG2, AFM1, and ST (sterigmatocystin) were purchased from Oskar Tropitzsch (Marktredwitz, Germany).

Strains and Sample Preparation

The A. flavus strains used in this study are listed in Table 1. Conidia from each strain were inoculated onto 10 cm petri plates containing solid Wickerham medium (≈25 ml per plate), which is comprised of the following: 2.0 g yeast extract, 3.0 g peptone, 5.0 g corn steep solids, 2.0 g dextrose, 30 g sucrose, 2 g NaNO3, 1 g K2HPO4⋅3H2O, 0.5 g MgSO4⋅7H2O, 0.2 g KCl, 0.1 g FeSO4⋅7H2O, 15 g agar per liter (pH 5.5). All cultures were incubated at 28°C in the dark for 7 days. The resulting fungal colonies were cut into small pieces with a scalpel and subsequently transferred to 500 ml screwcap Ehrlenmeyer flasks. Metabolites were extracted with 30 ml MeOH:DCM:EtOAc (10:20:30, v/v/v). The samples were agitated for 60 min on an Agitelec overhead shaker (J. Toulemonde & Cie, Paris, France). A total of 4 ml of extract was transferred to a glass tube and evaporated under a stream of nitrogen gas. The residue was reconstituted with 200 μl MeOH, and centrifuged in an Ultrafree®-MC centrifugal device for 5 min at 14,000 × g.

TABLE 1
www.frontiersin.org

Table 1. List of A. flavus strains used in this study.

HRMS Analysis

The experiments were carried out using a hybrid qTOF MS instrument, the AB SCIEX TripleTOF®4600 (Concord, Ontario, Canada), equipped with DuoSprayTM and coupled to an EksigentekspertTM ultraLC 100-XL system. The DuoSprayTM ion source [consisting of both electrospray ionization (ESI) and atmospheric pressure chemical ionization (APCI) probes] was operated in both, positive and negative ESI modes (ESI+/ESI-). The APCI probe was used for automated mass calibration using the calibrant delivery system (CDS). The CDS injects a calibration solution matching the polarity of ionization, and calibrates the mass axis of the TripleTOF® system in all of the scan functions used (MS and/or MS/MS). The qTOF HRMS method consisted of a full scan TOF survey (dwell time 100 ms, 100–1600 Da) and a maximum number of eight information-dependent acquisition (IDA) MS/MS scans (dwell time 50 ms). The MS parameters were as follows: curtain gas (CUR) 25 psi, nebulizer gas (GS 1) 50 psi, heated gas (GS 2) 60 psi, ion spray voltage (ISVF) 5.5 kV, interface heater temperature (TEM) 500°C, collision energy (CE) 10 V and declustering potential (DP) 70 V. For the IDA MS/MS experiments, a CE of 35 V was applied with a collision energy spread (CES) of 15 V. An EksigentekspertTM ultraLC 100-XL system was used for separation. The column was a ZORBAX RRHD Eclipse Plus C18 (1.8 μm, 2.1 mm × 100 mm) from Agilent Technologies (Diegem, Belgium). The mobile phase consisted of H2O:MeOH (95:5, v/v) containing 0.1% HCOOH and 10 mM HCOONH4 (solvent A) and MeOH:H2O (95:5, v/v) containing 0.1% HCOOH and 10 mM HCOONH4 (solvent B). The gradient elution program for UHPLC-qTOF HRMS analyses was applied as follows: 0–0.5 min: 0% B, 0.5–7 min: 0–99% B, 7–9 min: 99% B, 9–10 min: 99–0% B, 10–14 min: 0% B. The flow rate was 0.4 ml/min. The column temperature was set at 40°C and temperature of the autosampler was 4°C. A volume of 5 μl of each sample was injected. The instrument was controlled by Analyst® TF 1.6 software, while data processing was carried out using PeakView® 2.0 and MasterViewTM 1.0 software (all from AB Sciex).

Metabolite Screening Workflow

Accurate mass qTOF HRMS data were used to identify A. flavus metabolites, putatively by dereplication, using a customized database. Afterward, accurate mass HRMS data were combined with a careful evaluation of fragmentation spectra, to ascertain the presence of previously identified A. flavus metabolites and to establish an unambiguous identification strategy for subsequent screening work. The employed analytical methodology involved an untargeted data acquisition (consisting of full scan TOF HRMS survey and IDA MS/MS scans) and the processing of data using both targeted and untargeted approaches. Where possible, reference standards were used, and in this way compounds were identified in the fungal extracts by comparison of retention time (RT), accurate mass HRMS and HRMS/MS data with the authentic samples. A careful investigation and interpretation of fragmentation data was the basis for the identification of metabolites as described further.

Molecular Studies of 55 Examined A. flavus Strains

DNA was first extracted from each of the 55 A. flavus strains by inoculating their respective spore suspensions into each of 75 ml potato dextrose broth (PDB) in 250 ml Erlenmeyer flasks, which were placed on an orbital shaker (30°C, 135 rpm) for 24 h. Once substantial spore germination occurred, the resulting mycelia were vacuum filtered through a Büchner funnel using Miracloth to retain the mycelial tissue, and a wash was performed with sterile water. Mycelial tissue (200 mg) was collected into a 1.5 mL microcentrifuge tube, and a DNA extraction protocol was followed using the MasterPure Yeast Purification kit (Epicentre Biotechnologies, Madison, WI, United States). Optical density readings were performed on extracted genomic DNA using a Nano Drop (ND-1000) spectrophotometer, and dilutions (if necessary) were made in preparation for PCR amplification of unlinked genomic loci: acetamidase (amdS), beta-tubulin (benA), calmodulin (cmdA), and tryptophan synthase (trpC). The use of sequence variation at multiple, unlinked, conserved genomic loci allows for clone correction and is known as multi-locus sequence typing (MLST), and has been used in studies where haplotype associations account for genetic recombination (Olarte et al., 2011, Olarte et al., 2013; Moore et al., 2013). Another genomic region of interest spanned the area between the aflF and aflU genes from the aflatoxin gene cluster, which in A. flavus reportedly contains deletions associated with loss of G-aflatoxin production as well as with sclerotium morphotype (Ehrlich et al., 2004). For this study, no other aflatoxin cluster genes were examined. The mating-type (MAT) locus was also examined for these strains. Amplification of the aflF/aflU region and the MAT locus were solely diagnostic, and were not used for MLST in this study. The primers for each locus are listed in Supplementary Table S2. We used GoTaq Colorless Master Mix (Promega) to perform the PCR amplifications. Variable cycle conditions for each primer set are also listed in Supplementary Table S2. To test the success of PCR amplification, electrophoresis was performed using 5 μl of PCR product on a 1.5% agarose gel in TAE buffer. DNA amplicons were purified using a QIAquick PCR purification kit (Qiagen, Hilden, Germany) before sequencing. Nucleotide sequences for each locus were aligned, cleaned and trimmed using Sequencher version 5.4.2 (Gene Codes Corporation, Ann Arbor, MI, United States). The alignments for amdS (290 bp), benA (486 bp), cmdA (628 bp), and trpC (235 bp) sequences were then exported as nexus files for concatenation, haplotype associations through MLST, and phylogenetic analysis using SNAP Workbench (Price and Carbone, 2005). Once the nexus files were imported into SNAP Workbench, they were converted into phylip files using SNAP Map and then concatenated using SNAP Combine (Aylor et al., 2006). The concatenated alignment (1639 bp sequences) was collapsed into haplotypes using SNAP Map with the options of recoding indels (insertions/deletions) and exclusion of sites that violated the infinite sites model. Phylogenetic inference was based on heuristic parsimony searches in PAUP 4.0 (Swofford, 1998), which involved a branch-swapping algorithm. An A. parasiticus strain, SU1, was used as the outgroup taxon for our phylogenetic inference. Resulting most parsimonious trees were viewed with branch lengths based on mutational differences. Bootstrap analysis was conducted using PAUP 4 in SNAP Workbench, using default parmeters of 500 replicate iterations. Bootstrap values lower than 70 were considered unsupportive and would not be shown on the inferred tree.

Results

Metabolite Identification

Polyketides

Since aflatoxins remain one of the most studied groups of mycotoxins, with clearly established analytical parameters, their identification was straightforward. The identification of AFB1, AFB2, AFG1, and AFG2 is described in Supplementary Figure S2. Interestingly, we were able to demonstrate the presence of aflatoxin M1 (AFM1) and aflatoxin M2 (AFM2) as native metabolites in some of our fungal cultures of A. flavus (see Supplementary Figure S3 for details), although these compounds are mainly reported as metabolic derivatives of B-type aflatoxins found in different animal and human biological fluids (Prandini et al., 2009). A perfect match between tandem MS spectra of AFM1 and AFM2 found in our data (Supplementary Figure S3) and previous reports was observed (Plattner et al., 1984; Cavaliere et al., 2006; Huang et al., 2010; Díaz et al., 2011). The identified aflatoxins were confirmed by comparison of theirs MS data and RTs with those of authentic standards, highlighting the suitability of our proposed dereplication strategy for analyte detection and identification. In addition, we detected an aflatoxin derivative ([m+H]/z 331.0820, Δ = 2.4 ppm) having the same chemical composition as AFM2, but eluting just after it. The MS fragmentation pattern of this metabolite pointed to aflatoxin B2a or another functional isomer of AFM2. Its MS/MS spectrum clearly shows the loss of the hydroxyl group as a water molecule generating a very prominent ion at m/z 313 (Supplementary Figure S3C). Subsequently, the fragment ion at m/z 313 follows an identical fragmentation pathway as described for AFB1. This fragmentation pattern is perfectly in agreement with previous reports (Inoue et al., 2013; Rushing and Selim, 2017). However, without an authentic standard we could not assign it decisively as aflatoxin B2a.

Four other ultimate and penultimate precursors of aflatoxins could also be detected in fungal cultures of different strains, including: sterigmatocystin (ST), dihydrosterigmatocystin (DHST), O-methylsterigmatocystin (OMST), and dihydro-O-methylsterigmatocystin (DHOMST). All of the aforementioned precursors belong to the xanthone chemical subgroup of aflatoxin intermediates, possessing almost an identical core structure with very few minor chemical differences. These xanthone intermediates in the aflatoxin biosynthetic pathway share comparable and logical fragmentation patterns (Supplementary Figure S4). The fragmentation pattern we observed for these compounds is in strong accordance with previous publications (Bloom et al., 2007; Pfeiffer et al., 2014; Buiarelli et al., 2015). In this line of investigation, another derivative of OMST was putatively identified as being aspertoxin ([m+H]/z 355.0819, Δ = 1.8 ppm), which is actually a bisfuran oxygenated analog of OMST displaying a 16 Da difference between homologous fragment ions, such as 340 (324), 327 (311), 322 (306), 311 (295), and 294 (278) (Supplementary Figure S4E). Other aflatoxin precursors from the xanthone group of metabolites, such as demethylsterigmatocystin (DMST) and dihydro-demethylsterigmatocystin (DHDMST), could not be detected in our observations. Early stage intermediates in aflatoxin biosynthesis such as anthraquinones and other polyketide progenitors were not subject of this study. The extracted ion chromatograms (EICs) of the different A. flavus metabolites detected in this study are shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Extracted ion chromatograms (EICs) and elution order of the Aspergillus flavus metabolites detected in this study. Red arrows indicate the correct retention time.

Another important group of A. flavus polyketides are bicoumarins, belonging to the C3-C8′-linked chemical class of fungal bicoumarins recognized for their anti-insectan properties. Moreover, these molecules possess a very intriguing pathway of biosynthesis involving a dimerization mechanism of two coumaryl subunits via oxidative phenol-coupling reactions. These reactions are fundamentally controlled by dirigent P450 proteins that enable regio- and stereo-selective C-C-cross coupling of coumaryl monomers (Hüttel and Müller, 2007; Gil Girol et al., 2012). Six bicoumaryl analogs and four monomers are reportedly produced via the aflavarin gene cluster (cluster #39) in A. flavus (Cary et al., 2015a) (Supplementary Figure S1 and Supplementary Table S1). Those six bicoumarins (i.e., aflavarin and derived compounds) could also be detected in the present study (Supplementary Figures S5, S6); the MS data were as previously reported (Cary et al., 2015a) thereby unambiguously confirming the identity of these compounds.

A group of anthraquinones, different from those produced in the aflatoxin biosynthetic cascade, is reportedly synthesized via a polyketide gene cluster (cluster #27) in A. flavus (Cary et al., 2014). We detected and confirmed the identity of four of these compounds, i.e., asparasone A ([m-H]/z 357.0617, Δ = 3.5 ppm), oxyasparasone A ([m+H]/z 375.0723, Δ = 3.5 ppm), dehydroxyasparasone A ([m+H]/z 341.0669, Δ = 3.8 ppm) and deacetylasparasone A ([m-H]/z 315.0511, Δ = 2.4 ppm) (Supplementary Figure S7). This identification was also supported by findings from previous reports (Malysheva et al., 2014).

As described in Supplementary Figure S8, a series of isocoumarin metabolites were also identified in our investigations, including four derivatives of asperentin (6-O-methylasperentin, 8-O-methylasperentin, 4′-hydro- xyasperentin and 5′-hydroxyasperentin) and 6,8-dimethylcitreoisocoumarin. Asperentin and its derivatives chemically are structural analogs of isocoumarin and possess a 6′-methyltetrahydropyran-2′-methylene-2′-yl moiety attached to the C8 position of the isocoumaryl core system. Each of the pairs 6-O-methylasperentin/8-O-methylasperentin and 4′-hydroxyasperentin/5′-hydroxyasperentin represents functional isomers eluting at very close RT points (Figure 1) and possessing almost identical MS/MS spectra (Supplementary Figure S8).

Kojic acid is a small yet well-known metabolite produced by various Aspergillus species with wide application in cosmetics. The presence of this molecule (whose diagnostic fragment ions have been previously reported by Varga et al. (2013) and that of its acetylated derivative was confirmed in some of the examined A. flavus cultures as described in Supplementary Figure S9. Moreover, in ESI negative mode we were able to detect orsellinic acid in several A. flavus isolates ([m-H]/z 167.0354, Δ = 2.4 ppm). The MS/MS spectrum shown in Supplementary Figure S9C was in full accordance with spectra from a previously reported study of A. nidulans (Sanchez et al., 2010). To the best of our knowledge this is the first report of orsellinic acid in A. flavus.

Non-ribosomal Peptides

Within the heterogeneous class of non-ribosomal peptide metabolites produced by A. flavus, the aspergillic acid group comprises a number of closely-related pyrazinones with pronounced antibiotic activity (Dutcher, 1947, 1958). Moreover, this group of metabolites is recognized as hydroxamate siderophores of A. flavus due to the iron-chelating activity exerted by the cyclic hydroxamic acid motif in their chemical scaffolds (Assante et al., 1981; Zhu et al., 2011). The presence of these comppunds (i.e., aspergillic acid, neoaspergillic acid, hydroxyaspergillic acid, and neohydroxyaspergillic acid) in the analyzed A. flavus cultures was demonstrated as described in Supplementary Figure S10. The obtained spectra are in accordance with data from literature (Assante et al., 1981; Perry et al., 1984). RT assignment for isomeric compounds was supported by literature, where it has been shown that neo-forms of this class of pyrazinones elute earlier compared to their corresponding analogs. None of the metal-chelating complexes characteristic for this class of compounds, such as ferriaspergillin or aluminiumaspergillin, could be detected in our observations. Two homologous NRPS-like gene clusters in A. flavus (lna and lnb) have been shown to be responsible for the convergent biosynthesis of a pair of diastereomeric piperazines known as the cis- and trans-2,5-di-(p-hydroxybenzyl) piperazines. These compounds were detected in the investigated strains as described in Supplementary Figure S11A. We also deteted a compound (see Supplementary Figure S11B) that was previously demonstrated to be linked to the same biosynthetic machinery as the above piperazines (Forseth et al., 2013).

Ditryptophenaline, a homodimeric diketopiperazine alkaloid, has been identified as an important secondary metabolite in several A. flavus strains (Springer et al., 1977). This class of metabolites contains two hexahydropyrroloindole substructures, connected at vicinal quaternary stereocenters (C3-C3′ linkage), to form the dimeric molecules. Due to their astonishing chemical architecture, this and other related dimeric natural products exhibit various biological activities. In this context, it has been shown that ditryptophenaline inhibits substance P receptors, implying a potential analgesic or anti-inflammatory application. Moreover, it was recently demonstrated that a single cytochrome P450, DtpC, located within an A. flavus NRPS gene cluster, is responsible not only for pyrroloindole ring formation but also for catalyzing the C3-C3′ dimerization of N-methylphenylalanyl-tryptophanyl diketopiperazine monomers into a homodimeric final product (Saruwatari et al., 2014). The identification of ditryptophenaline in our study was quite straightforward due to its unique mode of fragmentation. The observed fragmentation pattern is described in Supplementary Figure S11C and is in accordance with previous studies (Springer et al., 1977; Barrow and Sedlock, 1994).

Hybrid Metabolites (PK-NRPs)

Chemical entities sculpted by mixed PKS/NRPS enzymatic systems have always been characterized as variable in terms of structural scaffolds and bioactive features. Genome sequencing has demonstrated that A. flavus contains two hybrid PKS/NRPS gene clusters (clusters #55 and #23). Cluster #55 is responsible for the synthesis of α-CPA (Seshime et al., 2009), while cluster #23 has been shown to be responsible for the production of leporins, a group of 2-pyridone metabolites (Cary et al., 2015b). The CPA class of alkaloids represents a heterogeneous group of indole-tetramate mycotoxins with pronounced implications in food safety and was indeed investigated in the present study. The identification of five of these CPA-type alkaloids (i.e., α-CPA, β-CPA, α-CPA imine, 2-oxoCPA, and speradine A) is described in Supplementary Figure S12. Screening of other CPA-type mycotoxins in different A. flavus strains, as well as their fragmentation behaviors, has been previously discussed in detail (Uka et al., 2017). With regard to the biosynthetic products of cluster 23, leporin B (a cyclic hydroxamic acid) was shown to be especially important due to its pronounced iron-chelating activity. MS analysis (Supplementary Figure S13) confirmed the production of leporin A ([m+H]/z 366.2078, Δ = 3,8 ppm), leporin B ([m+H]/z352.1905, Δ < 0,1 ppm), leporin C ([m+H]/z 336.1968, Δ = 2,9 ppm), 8-demethyl-leporin C ([m+H]/z 322.1814, Δ = 4,0 ppm), iron-trioxoleporin B ([m+H]/z 1107.4674, Δ = -1.5 ppm) as well as a hydroxylated derivative of leporin B ([m+H]/z 370.2012, Δ < 0,1 ppm) by the different A. flavus cultures investigated. The MS data were in accordance with previous reports (Arroyo-Manzanares et al., 2015; Cary et al., 2015b).

Indole-Diterpenoids

For many years it has been known that A. flavus produces metabolites with potent tremorgenic activity. Through subsequent investigations, it was determined that this class of molecules share a basic indole-diterpene ring system generated by the fusion of a polycyclic diterpene, derived from a geranylgeranyldiphosphate (ggpp), and an indole moiety with tryptophan origin. These intriguing molecules are synthesized by terpene/terpenoid cyclases, very complex groups of enzymes which contribute crucially to the tremendous structural diversity of this large family of natural products. In the A. flavus secondary metabolome, two different chemical groups of indole-diterpenes have been described, aflavinines and aflatrem, along with their derivatives. Aflavinines possess a very interesting and sterically-congested chemical structure, in which the indole moiety and the diterpenoid ring system are coupled via a single C-C chemical bond. We identified aflavinine ([m+H]/z 406.3107, Δ = 0,74 ppm), 14-hydroxyaflavinine ([m+H]/z 422.3053, Δ < 0,1 ppm) and its isomers, as well as 14,25-dihydroxyaflavinine ([m+H]/z 438.3017, Δ = 3.2 ppm) in cultures of the A. flavus strains analyzed in this study. A detailed description of the confirmation of the identity of those compounds is given in Supplementary Figure S14. Regarding 14-hydroxyaflavinine, interestingly, there are four metabolites eluting at different RTs (RT 5.7, 5.9, 6.3, and 6.5 min), all having the same exact mass and sharing the fragmentation pattern of monohydroxyaflavinine. This is most likely due to the position of the double bond in ring C, as well as the epimerization of the hydroxyl group at C14, resulting in four isomers of monohydroxyaflavinine: 14-hydroxyaflavinine, 14-epi-hydroxyaflavinine, 14-hydroxyisoaflavinine and 14-epi-hydroxyisoaflavinine. The elution order of these metabolites cannot be determined unambiguously; however, reports from the literature indicate that iso-forms elute earlier, as assigned in Figure 1. A similar mode of fragmentation has been encountered for 14,25-dihydroxyaflavinine with three consecutive losses of hydroxyl functionalities, followed by the chronological fragmentation of the C20-terpenoid moiety (Supplementary Figure S14C).

Aflatrem is another A. flavus metabolite possessing a prenylated heptacyclic indole-diterpenoid skeleton with very potent tremorgenic toxicity. Moreover, it was demonstrated that aflatrem is synthesized through a very complex and convergent biosynthetic pathway, governed by two different gene clusters located on separate chromosomes of A. flavus (Nicholson et al., 2009). In our investigation of the A. flavus strains, a compound with [m+H]/z 502.2952, cooresponding to aflatrem (Δ < 0.1 ppm) could be conspicuously detected in significant abundance. Accurate mass data, combined with the fragmentation behavior, decisively confirmed the presence of aflatrem as described in Supplementary Figure S15. Apart from aflatrem, β-aflatrem is another structural isomer described in the literature, which was also detected in our study. The only structural difference between these two isomers is the position of the dimethylallyl moiety in the indole core system, thus exhibiting the same fragmentation pattern as described for aflatrem. Based on previous reports, aflatrem elutes earlier in reverse-phase chromatography (Nicholson et al., 2009). Similarly, hydroxyaflatrem ([m+H]/z 518.2917, Δ = 3.0 ppm), an oxygenated analog of aflatrem, could be detected in several strains of A. flavus. Assignment of the position of this hydroxyl group in the aflatrem prototype structure is not so straightforward, due to the number of possible options. Nevertheless, MS/MS data (Supplementary Figure S15B) indicated that this functionality most likely is attached in the indole moiety of the core structure. The identity of other aflatrem derivatives such as paspaline, paspalinine, or paspalicine could not be decisively confirmed in our experiments.

Distribution of Metabolites as a Function of Genetic and Morphological Features

Occurrence of A. flavus secondary metabolites as function of mating-type genes in different strains is determined as shown in Figure 2 and Supplementary Tables S3S5. The heat map data visualization on a targeted list of metabolites shows the clustering of the 55 A. flavus isolates into four groups (clades). From right to left of the heat map, the first group seems to produce bicoumarins, anthraquinones, NRPs, PK-NRPs, and indole-diterpenoids but less often aflatoxins. The second group, seems to lack the ability in producing bicoumarins in addition to aflatoxin absence. The third group produces almost all classes of secondary metabolites. Whereas in the 4th group a higher frequency for NRPs and PK-NRPs is observed but less polyketides and indole-diterpenoids. Based on this limited list of known metabolites, there was no clear correlation/clustering between mating-type genes and secondary metabolite production. On the other hand distribution of metabolites as function of sclerotium size showed a better clustering pattern (Figure 3).

FIGURE 2
www.frontiersin.org

Figure 2. Heat map representing metabolite distribution in 55 A. flavus isolates. On the top of the figure, a dendogram is depicted, clustering A. flavus isolates as function of mating-type genes (MAT1-1 and MAT1-2) and metabolite occurrence. This figure was generated using the online metabolomics platform MetaboAnalyst 3.0 after the data were treated accordingly using sum normalization, log transformations and Pareto scaling. The data table was prepared using area under the curve (AUC) for each chromatographic peak of respective metabolites. Anq, anthraquinones; PKs, polyketides; NRPs, non-ribosomal peptides; IND-DITs, indole-diterpenoids; ST, sterigmatocystin; OMST, O-methylsterigmatocystin; DHST, dihydro-ST; DHOMST, dihydro-OMST; Asparas., asparasone A; Asper. OCH3, asperentin methyl-ether; 6,8-DMCIC, 6,8-dimethylcitreoisocoumarin; Asp. acid, aspergillic acid; Neoasp. acid, neoaspergillic acid; α-CPA, alfa-cyclopiazonic acid; Lep. B, leporin B.

FIGURE 3
www.frontiersin.org

Figure 3. Heat map representing metabolite distribution in 27 A. flavus isolates with assigned sclerotium morphotype. On the top of the figure, a dendogram is depicted, clustering A. flavus isolates as function of sclerotium morphotype: L (large sclerotia > 400 μm) and S (small sclerotia < 400 μm). This figure was generated using the online metabolomics platform MetaboAnalyst 3.0 after the data were treated accordingly using sum normalization, log transformations and Pareto scaling. The data table was prepared using AUC foreach chromatographic peak of respective metabolites. Anq, anthraquinones; PKs, polyketides; NRPs, non-ribosomal peptides; IND-DITs, indole-diterpenoids; ST, sterigmatocystin; OMST, O-methylsterigmatocystin; DHST, dihydro-ST; DHOMST, dihydro-OMST; Asparas., asparasone A; Asper. OCH3, asperentin methyl-ether; 6,8-DMCIC, 6,8-dimethylcitreoisocoumarin; Asp. acid, aspergillic acid; Neoasp. acid, neoaspergillic acid; alfa-CPA, alfa-cyclopiazonic acid; Lep. B, leporin B.

Molecular Studies of 55 Examined A. flavus Strains

Some of the examined Aspergillus strains had been previously accessioned for one or more of the four conserved loci of interest. A list of all GenBank accession numbers is shown in Supplementary Table S6. Of the 55 A. flavus strains examined, four strains failed to properly amplify and/or sequence for one of the four loci. Despite repeated attempts, strain 150 failed to result in a useable sequence for the trpC locus, while strains 1568, 1574, and 1575 failed to result in useable sequences for the cmdA locus. The sequences obtained for one strain, 1118, across all four loci appeared highly polymorphic in relation to those of the other strains. BLAST query of each locus for 1118 revealed it to be an A. nomius strain. This can only be explained by an inaccurate morphological identification, so 1118’s species identification has been changed in the SRRC’s fungal database to accurately reflect our genomic identifications. Results of BLAST queries for each strain’s locus sequences are shown in Supplementary Table S7. Most of the strains were found to be A. flavus; however, nucleotide sequences for some strains appeared to have variable species designations across the four loci (1576, 1591, 2000, and 2033). Amplification of their aflF/aflU genomic regions revealed deletions that show 1576 to be an A. flavus S-strain (1500 bp deletion) and 2000 to be an A. flavus L-strain (1000 bp deletion). This indicates that strains 1576 and 2000 should not produce G aflatoxins, yet they share species identity for the benA locus with A. minisclerotigenes and A. parvisclerotigenus, respectively, which reportedly produce B and G aflatoxins (Frisvad et al., 2005; Pildain et al., 2008). It is possible for A. parvsclerotigenus to lack production of G aflatoxins (Frisvad et al., 2005). Strains 1591 and 2033 do produce G aflatoxins, therefore their respective species accessions for the benA locus (A. parvisclerotigenus) and amdS locus (A. minisclerotigenes) could be accurate. The amplification of intact aflF/aflU regions (i.e., no deletions present) in 1591 and 2033 further corroborated their production of G aflatoxins.

One of the consequences of a strain not being represented for one or more loci in a concatenated sequence is exclusion from population analyses, therefore, the four strains mentioned above (150, 1568, 1574, 1575) were not included in our multi-locus inferences. We also left out the newly-identified A. nomius strain, 1118. The haplotype designations from our revised sample of 51 Aspergillus strains (including SU1) can be seen in Table 2. A total of 43 haplotypes (i.e., individuals) resulted from collapsing 51 concatenated sequences. Some haplotypes include multiple isolates that share sequence homology across all four unlinked genomic loci. A potential bottleneck with population analysis results when high complexity of sequence data proves too much for certain analysis software. We experienced this during phylogenetic analysis using the stand-alone version of SNAP Workbench. Bootstrap analysis would not complete. To addres this issue we used an online version of the software that allows runs to be placed on a cluster of networks, which accelerates analysis. We were able to complete bootstrap analysis (based on 1000 runs), and then view/export our tree using RAxML 8 (Stamatakis, 2014), accessible through the CIPRES RESTful API (Miller et al., 2015), implemented in the Mobyle SNAP Workbench v1.55 (Monacell and Carbone, 2014); the tree was visualized using an online portal software known as T-BAS v. 2.1 (Carbone et al., 2017). Figure 4 shows inferred phylogenetic associations examined in this study. No clear cladal distinctions could be made on the basis of mating type, sclerotium morphotype, aflF/aflU-deletion type, or metabolic profile. Clade I includes the A. parasiticus Type strain, SU1. It was inferred to share a most recent common ancestor with strain 1591, which is a B+G aflatoxing-producing S-strain that had been identified as both A. flavus, based on BLAST query of its amdS, cmdA and trpC loci, and as A. parvisclerotigenus based on its benA locus. Bootstrap support was moderate, having a value of 70. Strain 1591 was the only one in this entire study found to produce leporin A, a compound that has only been described in A. flavus and A. leporis (Frisvad et al., 2019). Strain 2033 shares a common ancestor with that of SU1 and 1591, having a boostrap value of 73. Although there is much metabolic overlap across the three strains, including the production of M aflatoxin, 2033 fails to produce any of the examined bicoumarins or non-ribosomal peptides. Strain 2033 is a B+G aflatoxin producing S-strain that had been identified as both A. flavus, based on BLAST query of its benA, cmdA and trpC loci, and as A. minisclerotigenes based on its cmdA locus. The remaining strains in Clade I are mostly non-aflatoxigenic and share the S-type deletion in their aflF/aflU region, supporting their lack of G-aflatoxin production, but not all exhibit the S-morphotype for sclerotium production. Strain 1547 produces L-type sclerotia, while strains 283 and 2111 failed to produce any sclerotia. Strain 2111 is a high producer of B and M1 aflatoxins. Three of the strains in this clade (1544, 1557, 1571) share sequence identity across all four conserved loci, are all of the same sclerotium morphotype, have the same S-type deletion in aflF/aflU, and are all non-aflatoxigenic. One observed difference among them relates to mating type, since strain 1544 is MAT1-2 while the others are MAT1-1. Other observed differences are metabolic, such as our observation of 1544 producing only two isocoumarin compounds and six PK-NRP hybrid molecules, or that only strain 1557 produced any of the examined bicoumarin and indole-diterpene compounds, or that only 1571 produced orsellinic acid.

TABLE 2
www.frontiersin.org

Table 2. Haplotype identities for 51 Aspergillus strains examined in this study.

FIGURE 4
www.frontiersin.org

Figure 4. Phylogenetic tree inferred for the A. flavus isolates examined in this study. Strain SU1 is the outgroup taxon, A. parasiticus. Phylogenetic inference is based on the concantenation of four unlinked genomic loci (amdS, benA, cmdA, and trpC).

Clade II encompasses most of the examined isolates (n = 35), representing 29 haplotypes, and includes the two commercially-available biocontrol strains (AF36 = 1533 and NRRL 21882 = 1534). Although this clade has a nearly equal distribution of aflF/aflU-deletion types (800 bp/L-type = 18; 1500 bp/S-type = 17), no partitioning was observed based on this phenotype. Similarly, no branch partitioning could be observed that related to mating type, sclerotium morphotype or aflatoxin chemotype. Interestingly, one strain (2114) has the L-type deletion in aflF/aflU, but was observed to produce G aflatoxins. This strain produced no detectable bicoumarin compounds, nor did it produce aflatrems, NRPs or PK-NRPs. Clades I and II share a most recent common ancestor, but there was a lack of bootstrap support for this. There was, however, bootstrap support for the convergence of Clades I and II with Clade III, having a value of 72.

Clade III includes only four strains (295, 1357, 1541, 1552). These strains mainly share one phenotype, which is the S-type deletion in their aflF/aflU regions. With regard to aflatoxin producing ability, 1357 and 1541 are non-aflatoxigenic, while 295 and 1552 produce a small amount of B1 aflatoxin. In comparing the two non-aflatoxigenice strains, 1357 lacks production of all but two of the examined polyketides (Supplementary Table S3), and 1541 produces only one polyketide (orsellinic acid). Strain 1541 produces greater numbers and quantities of aspergillic acid compounds (Supplementary Table S4), but no flavocol or ditryptophenaline (all non-ribosomal peptides). Nor does it produce hybrid molecules beyond α-CPA, leporins B and C, and iron trioxoleporin B. Neither strain produces indole-diterpenes (Supplementary Table S5). Aflatoxin B1 producers, 295 and 1552, share similar polyketide production profiles, except for the production of asparsone 317 and orsellinic acid (only in 295) and 6,8-dimethylcitreo isocourmarin (only in 1552). Strain 1552 only produces the non-ribosomal peptide ditryptophenaline, and it lacks production of β-CPA among the hybrid molecules, while 295 lacks production of leporins B and C (Supplementary Table S4). Among the indole-diterpenes (Supplementary Table S5), both strains produce aflavinines, but only 295 produces aflatrems. There is also no clustering based on mating type, as only 1357 is a MAT1-2 strain while the other three are MAT1-1.

Clade IV includes two strains (141, 1626) that appear basal to the remaining sample population. These strains share the same mating type, aflF/aflU deletion type (S), and both fail to produce sclerotia. Among the polyketides, strain 141 was observed to produce B and M aflatoxins, while 1626 is non-aflatoxigenic (Supplementary Table S3). Neither strain produces bicoumarin compounds, and only 141 produces the asparasone compounds and orsellinic acid. Both strains produce the isocoumarin compound known as asperentin-6(8)-methyl ether, but only 1626 produces 6,8-dimethylcitreo isocourmarin. Both strains produce only a single non-ribosomal peptide (ditryptophenaline), but 1626 produces a greater number of hybrid molecule compounds (Supplementary Table S4). Strain 1626 produces no indole-diterpenes, while 141 produces all but hydroxyaflatrem (Supplementary Table S5).

Discussion

Occurrence of various metabolites in different A. flavus individuals was unambiguously determined using a dereplication approach, based on HRMS data combined with tandem MS spectra thereby allowing a confident identification of several compounds. Dereplication is an especially important and useful tool for natural product chemists to avoid re-isolation of already known metabolites, thus saving time and resources. Dereplication analytical methodologies are most often implemented through hybrid HRMS mass spectrometers, which are automated approaches known to be fast, robust and have a strong identification potential (Wolfender et al., 2000; Fredenhagen et al., 2005; Nielsen et al., 2011; El-Elimat et al., 2013). Less often, other analytical techniques such as UV spectroscopy and NMR have been used in dereplication strategies (Bobzin et al., 2000; Wolfender et al., 2003). In the present study we constructed our in-house database with specific A. flavus secondary metabolites, gathering from literature all possible parameters such as chemical formulas, MS/MS spectra and elution order, which then were integrated in MS-software packages to proceed with the screening work. As expected, we observed a scattered and heterogeneous occurrence of compounds in different strains of A. flavus, which correlates with its high genetic diversity. This was especially true for the polyketide class of metabolites (Supplementary Table S3). Publications in the literature indicate great variability between different strains of A. flavus in their ability to produce aflatoxins (Hesseltine et al., 1968; Geiser et al., 2000; Ehrlich et al., 2007). In general, aflatoxigenic A. flavus isolates are predominantly known as B-type aflatoxin producers that are incapable of producing G-type aflatoxins (Taber and Schroeder, 1967; Ehrlich et al., 2007). In this study of 55 A. flavus strains, 28 were aflatoxin producers. Only four of those 28 isolates produced all (B-, G-, and M-types) of the aflatoxins; 13 isolates produced only B- and M-type aflatoxins; in one strain we detected production of only AFB1 and AFB2a; 10 strains produced B-type aflatoxins alone; and 27 strains produced no detectable aflatoxins. We did not observe a single strain that could produce G- or M-type aflatoxins without also producing B-type aflatoxins. Our detection of aflatoxins M1 and M2 in fungal cultures is of particular interest because M-type aflatoxins were originally found, and are most often reported being sampled, in milk and other biological fluids (Holzapfel et al., 1966; Masri et al., 1967; Hesseltine et al., 1968; Stubblefield et al., 1970; Diener et al., 1987). These findings indicate that occurrence of M-type aflatoxins in different food matrices does not necessarily require animal or human involvement for their production, because we observed these hydroxylated aflatoxins as endogenous metabolites of in situ A. flavus. The ability of Aspergillus spp. to produce M-aflatoxins has already been previously demonstrated (Yabe et al., 2012). The authors induced AFM1 and AFM2 production in A. parasiticus by feeding experiments. Feeding of A. parasiticus with aspertoxin (12c-hydroxyOMST) resulted in the production among others of these two metabolites by the fungus. Furthermore, cell-free experiments using the microsomal fraction of A. parasiticus and aspertoxin also led to the production of AFM1 and AFM2, indicating that aspertoxin is a precursor of these metabolites. Interestingly, incubation of either the microsomal fraction or OrdA enzyme-expressing yeast with OMST or DHOMST resulted in the production of aspertoxin together with the M-type aflatoxins, indicating that the OrdA enzyme catalyzes both 12c-hydroxylation reaction from OMST/DHOMST to aspertoxin/dihydroaspertoxin and the subsequent reaction from aspertoxin/ dihydroaspertoxin to AFM1/AFM2. Production M-aflatoxins by the investigated A. flavus strains should theorically be possible through mediation of OrdA enzyme; our results decisively confirming this, highlighting the power and the usefulness of mass spectrometry untargeted analysis in providing a global and better insight into the different metabolites that could be present in a sample. Several survey studies have reported presence of AFM1 in different food matrices (Chala et al., 2014; Sulyok et al., 2015; Abdallah et al., 2017). Furthermore, based on our findings it seems that A. flavus more frequently produces M-type aflatoxins compared to aflatoxin G1 and G2. Actually, the ability/inability of A. flavus isolates to produce G-type aflatoxins may relate to a deletion in the encoding gene for a cytochrome P450 enzyme (aflU) within the pathway of the aflatoxin biosynthetic cluster (Ehrlich et al., 2004). More precisely, these deletions have been reported in the aflF/aflU region of the aflatoxin gene cluster. Three out of the four G-producing strains examined in this study had no deletion in their aflF/aflU regions, which justifies their G-producing ability. Our finding of G-aflatoxin production by strain 2114, despite possessing a 1000 bp deletion in its aflF/aflU region, should be explored further. Perhaps there is another gene present that allows for the production of G aflatoxins when the aflU gene is deleted.

Aflatoxin B- and G-producing isolates from Thailand, which morphologically resemble A. flavus, were definitively assigned as a novel clade of A. nomius (Ehrlich et al., 2007). Perhaps this morphological similarity is why strain 1118 was misidentified as A. flavus. It is becoming more commonly accepted that A. flavus S-type isolates, that are also B+G aflatoxin producers, belong to one of the three phylogenetically related species: A. minisclerotigenes (Pildain et al., 2008), A. parvisclerotigenus (Frisvad et al., 2005, 2019), and A. korhogoensis (Carvajal-Campos et al., 2017). Several recent studies have characterized “novel” species that produce B and G aflatoxins such as A. mottae, A. transmontanensis, A. sergii, A. pseudonomius, and A. novoparasiticus (Varga et al., 2011; Gonçalves et al., 2012; Soares et al., 2012). Moore and co-workers are seeking to determine whether or not these novel species, which exhibit striking morphological and metabolic similarities to A. flavus and A. parasiticus, are actually hybrid species through genomic sequencing of the Type strains for all aflatoxigenic species (Moore et al., 2015a,b, 2016). We may find that aflatoxigenic fungi from Section Flavi form a species complex whose constituent “species” are capable of inter-specific recombination, and that recombination between closely-related Aspergillus species contributes to the genus’ genetic and metabolic diversity. As well, evolutionary relationships for their respective metabolomes can be further investigated.

Regarding other polyketide metabolites, we found a consistent presence of bicoumarins and anthraquinones. Among non-ribosomal peptides, ditryptophenaline was the most prevalent metabolite, occurring in 47 of our examined strains, followed by the aspergillic acid group of mycotoxins (Supplementary Table S4). Regarding hybrid polyketide-aminoacid molecules, both the CPA and 2-pyridone groups of metabolites were highly present (Supplementary Table S4). Detailed occurrence of CPA-type alkaloids in this set of A. flavus isolates was reported previously (Uka et al., 2017). Leporin B, a cyclic hydroxamic acid iron chelator, was the most frequent metabolite observed, from within the 2-pyridone class of metabolites. However, we determined the presence of indole-diterpenoids in 42 of the examined strains (Supplementary Table S5). Among these 42, 38 individuals were capable of producing both the aflavinine- and aflatrem-type of indole-diterpenoids.

All A. flavus soil populations have the potential to include strains representing two sclerotial morphotypes: the L-strain isolates with average sclerotial size greater than 400 μm, and the S-strain isolates with sclerotial size less than 400 μm (Cotty, 1997). S-strain isolates reportedly produce higher levels of aflatoxins, more abundant sclerotia and fewer conidia (Probst et al., 2012). In this study we were unable to definitely ascertain sclerotium production for all examined strains; nevertheless, based on the limited number of strains the dendogram clustering analysis show a quite satisfactory grouping of individual isolates by function of sclerotium morphotype (Figure 3). This indicates that sclerotium morphotype is likely an important morphological feature influencing production of secondary metabolites. The aflF/aflU region of the aflatoxin cluster has been associated with sclerotium morphotype (Ehrlich et al., 2004), but we found at least one contradiction to this association in our examination of these A. flavus strains.

Although mycotoxin contamination may occur during transport or storage of food products, pre-harvest contamination from fungal infections in plant hosts is another source of mycotoxins. In order to reduce aflatoxin contamination during pre-harvest periods, a number of mitigation strategies like fungicide application, development of resistant plants and biological control have been explored. Biocontrol approaches, involving the introduction of non-aflatoxigenic A. flavus isolates into agricultural fields to displace aflatoxin-producing strains, have been tested extensively in the last two decades (Dorner et al., 2003; Abbas et al., 2011; Atehnkeng et al., 2014). Apart from aflatoxins, the inability of A. flavus biocontrol agents to produce other relevant mycotoxins, among which are α-CPA, ditryptophenaline, aflatrem and aflavinines, is highly desirable due to the unknown long-term and cumulative toxicological effects of these metabolites. Two commercially available biocontrol agents, AF36 (SRRC strain 1533) and the component strain in Afla-Guard®, NRRL 21882 (SRRC strain 1534), were included in this study. AF36 was originally isolated from a cotton field in Arizona and is approved for application on various commodities in Arizona and California. Loss of aflatoxigenicity in AF36 is the result of a nonsense mutation in its aflC, a pathway gene involved in aflatoxin biosynthesis (Ehrlich and Cotty, 2004). AF36 has an otherwise complete aflatoxin cluster. NRRL 21882 lacks the entire aflatoxin gene cluster (Chang et al., 2005). Apart from CPA-type mycotoxins, most of which were recently investigated (Uka et al., 2017), we could detect many more metabolites as part of the secondary metabolome of these aforementioned biological controls. Secondary metabolites like ditryptophenaline, aflatrem-type tremorgens, the aspergillic acid group of mycotoxins and 2-pyridones were detected in both biocontrol agents. However, aflavarin and aflavinine-type tremorgenic compounds were exclusively detected in 1534 (Supplementary Tables S3S5). Our findings suggest that these biocontrol agents are very active in producing secondary metabolites with different toxicological profiles, an issue which should not be ignored when considering their application to agricultural fields.

As was stated above, populations of A. flavus are characterized by considerable diversity in terms of morphological, functional, and genetic features, which is clearly supported by the large number of VCGs and different mycotoxin chemotypes. The large number of VCGs for A. flavus should significantly limit the frequencies of hyphal anastomosis, except when the necessity to undergo sexual out-crossing (i.e., to circumvent unfavorable conditions) relaxes those incompatibility barriers and contributes greatly to high genetic diversity among individual isolates. It has been reported that species populations such as A. tamarii or A. parasiticus possess lower VCG diversity, and consequently experience a higher frequency of hyphal contact (Horn and Greene, 1995). The high variation between A. flavus strains is not easily justified in the absence of sexuality. Although for many years A. flavus was considered an asexual fungus, researchers have shown evidence that at best these fungi can exhibit a predominantly asexual existence. We now know that A. flavus contains functional genes for mating. The expression of these genes has been demonstrated at the mRNA and protein levels. Following the same logic, it would be of great interest to decipher the expression of these genes in relation to secondary metabolism, aiming to establish any type of relationship between a single or a group of secondary metabolites with one or the other MAT gene. In addition, elucidation of these correlations between MAT genes and secondary metabolism is of great importance since the general indications that sexual reproduction is occurring between strains with differing capacities for mycotoxin production. Based on our analyses, inferred from a limited list of known metabolites, we could not observe any particular relationship between any of the mating-type genes and secondary metabolite production (Figure 2). This could be also due to the fact that other possible unknown metabolites not included in our list are responsible for making the real difference. Or maybe these results are just reflecting the reality that individual mating-type genes do not influence expression of any specific secondary metabolite gene cluster. In this context, it could be that comparing the secondary metabolomes of individual strains possessing opposite mating-type genes is not the right strategy to judge about the real contribution of sexuality in secondary metabolite chemical diversity. Hence, other biological parameters in A. flavus populations like sexual/asexual ratio could be a more relevant feature to track the correlations between sex and chemical diversity. An example could be to compare the secondary metabolome of a field population of A. flavus isolates with an almost equal distribution of MAT1-1 and MAT1-2 genes (sexual recombination) with another population which predominantly undergoes asexual recombination (unequal occurrence of MAT1-1 and MAT1-2 genes). In this way, we would be able to evaluate more precisely the real impact of sexuality in secondary metabolite and mycotoxin production.

Phylogenetic analysis did not offer much resolution for the examined strains. Lack of strong boostrap support and partitioning of haplotypes based on a particular phenotype, underscore the complexity of A. flavus strains with regard to morphology, genetics and production of metabolites. Many of the strains examined are from different geographic locations and/or different hosts. Perhaps the best inferences are made from sample populations that are considered panmictic (i.e., from the same location and host type), which might add more power to the analyses. Our findings support the importance of incorporating a holistic approach to strain/species identification. If we had not undertaken molecular research, we would have assumed all of our strains were A. flavus. What fungal researchers have to be careful of is assuming that just because two strains produce different metabolites it does not mean they are different species. The same should be said for assuming that two strains that are similar in appearance are the same species. With the new found concepts of intra- and inter-specific hybridization for Aspergillus fungi, the possibilities are endless. Therefore, species identifications should encompass morphological, metabolic and genetic comparisons. Furthermore, genetic associations should include multiple loci (both individual and combined).

Author Contributions

JM, GM, and VU, conceived and designed the experiments. JM and GM, generated the data. VU, GM, NA-M, and JM, analyzed the data. VU, GM, DN, SS, and JM, wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00667/full#supplementary-material

References

Abbas, H. K., Zablotowicz, R. M., Horn, B. W., Phillips, N. A., Johnson, B. J., Jin, X., et al. (2011). Comparison of major biocontrol strains of non-aflatoxigenic Aspergillus flavus for the reduction of aflatoxins and cyclopiazonic acid in maize. Food Addit. Contam. Part A Chem. Anal. Control Expo Risk Assess. 28, 198–208. doi: 10.1080/19440049.2010.544680

PubMed Abstract | CrossRef Full Text | Google Scholar

Abdallah, M. F., Girgin, G., Baydar, T., Krska, R., and Sulyok, M. (2017). Occurrence of multiple mycotoxins and other fungal metabolites in animal feed and maize samples from Egypt using LC-MS/MS. J. Sci. Food Agric. 97, 4419–4428. doi: 10.1002/jsfa.8293

PubMed Abstract | CrossRef Full Text | Google Scholar

Amani, S., Shams-Ghahfarokhi, M., Banasaz, M., and Razzaghi-Abyaneh, M. (2012). Mycotoxin-producing ability and chemotype diversity of Aspergillus section Flavi from soils of peanut-growing regions in Iran. Indian J. Microbiol. 52, 551–556. doi: 10.1007/s12088-012-0275-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Arroyo-Manzanares, N., Diana Di Mavungu, J., Uka, V., Malysheva, S. V., Cary, J. W., Ehrlich, K. C., et al. (2015). Use of UHPLC high-resolution Orbitrap mass spectrometry to investigate the genes involved in the production of secondary metabolites in Aspergillus flavus. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 32, 1656–1673. doi: 10.1080/19440049.2015.1071499

PubMed Abstract | CrossRef Full Text | Google Scholar

Assante, G., Camarda, L., Locci, R., Merlini, L., Nasini, G., and Papadopoulos, E. (1981). Isolation and structure of red pigments from Aspergillus flavus and related species, grown on a differential medium. J. Agric. Food Chem. 29, 785–787. doi: 10.1021/jf00106a023

CrossRef Full Text | Google Scholar

Atehnkeng, J., Ojiambob, P. S., Cotty, P. J., and Bandyopadhyay, R. (2014). Field efficacy of a mixture of atoxigenic Aspergillus flavus Link: FR vegetative compatibility groups in preventing aflatoxin contamination in maize (Zea mays L.). Biol. Control 72, 62–70. doi: 10.1016/j.biocontrol.2014.02.009

CrossRef Full Text | Google Scholar

Aylor, D. L., Price, E. W., and Carbone, I. (2006). SNAP: combine and map modules for multilocus population genetic analysis. Bioinformatics 22, 1399–1401. doi: 10.1093/bioinformatics/btl136

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrow, C. J., and Sedlock, D. M. (1994). 1’-(2-Phenyl-ethylene)-ditryptophenaline, a New dimeric diketopiperazine from Aspergillus flavus. J. Nat. Prod. 57, 1239–1244. doi: 10.1021/np50111a008

CrossRef Full Text | Google Scholar

Basappa, S. C., Sreenivasamurthy, V., and Parpia, H. A. (1970). Aflatoxin and kojic acid production by resting cells of Aspergillus flavus Link. J. Gen. Microbiol. 61, 81–86 doi: 10.1099/00221287-61-1-81

PubMed Abstract | CrossRef Full Text | Google Scholar

Bennett, J. W. (2010). “An overview of the genus Aspergillus,” in Aspergillus: Molecular Biology and Genomics, eds M. Machida and K. Gomi (Norfolk, VA: Caister Academic Press), 1–17.

Google Scholar

Bloom, E., Bal, K., Nyman, E., Must, A., and Larsson, L. (2007). Mass spectrometry-based strategy for direct detection and quantification of some mycotoxins produced by Stachybotrys and Aspergillus spp. in indoor environments. Appl. Environ. Microbiol. 73, 4211–4217. doi: 10.1128/aem.00343-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobzin, S. C., Yang, S., and Kasten, T. P. (2000). LC-NMR: a new tool to expedite the dereplication and identification of natural products. J. Ind. Microbiol. Biotechnol. 25, 342–345. doi: 10.1038/sj.jim.7000057

PubMed Abstract | CrossRef Full Text | Google Scholar

Bűchi, G., Kobbe, B., Demain, A. L., and Clardy, J. (1977). The structure of ditryptophenaline - a new metabolite of Aspergillus flavus. Tetrahedron Lett. 18, 2403–2406. doi: 10.1016/s0040-4039

CrossRef Full Text | Google Scholar

Buiarelli, F., Di Filippo, P., Riccardi, C., Pomata, D., Rumolo, E., Giannetti, L., et al. (2015). Analytical method for the determination of mycotoxins in indoor/outdoor airborne particulate matter by HPLC-MS-MS. Int. J. Environ. Anal. Chem. 95, 713–729. doi: 10.1080/03067319.2015.1055468

CrossRef Full Text | Google Scholar

Burdock, G. A., and Flamm, W. G. (2000). Safety assessment of the mycotoxin cyclopiazonic acid. Int. J. Toxicol. 19, 195–218. doi: 10.1080/10915810050074964

CrossRef Full Text | Google Scholar

Bush, M. T., Touster, O., and Brockman, J. E. (1951). The production of β-nitropropionic acid by a strain of Aspergillus flavus. J. Biol. Chem. 188, 685–693.

Google Scholar

Campbell, M. A., Rokas, A., and Slot, J. C. (2012). Horizontal transfer and death of a fungal secondary metabolic gene cluster. Genome Biol. Evol. 4, 289–293. doi: 10.1093/gbe/evs011

PubMed Abstract | CrossRef Full Text | Google Scholar

Carbone, I., Jakobek, J. L., Ramirez-Prado, J. H., and Horn, B. W. (2007). Recombination, balancing selection and adaptive evolution in the aflatoxin gene cluster of Aspergillus parasiticus. Mol. Ecol. 16, 4401–4417. doi: 10.1111/j.1365-294x.2007.03464.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Carbone, I., White, J. B., Miadlikowska, J., Arnold, A. E., Miller, M. A., Kauff, F., et al. (2017). T-BAS: tree-Based Alignment Selector toolkit for phylogenetic-based placement, alignment downloads, and metadata visualization: an example with the Pezizomycotina tree of life. Bioinformatics 33, 1160–1168. doi: 10.1093/bioinformatics/btw808

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvajal-Campos, A., Manizan, A. L., Tadrist, S., Akaki, D. K., Koffi-Nevry, R., Moore, G. G., et al. (2017). Aspergillus korhogoensis, a novel aflatoxin producing species from the Côte d’Ivoire. Toxins 9:E353. doi: 10.3390/toxins9110353

PubMed Abstract | CrossRef Full Text | Google Scholar

Cary, J. W., Han, Z., Yin, Y., Lohmar, J. M., Shantappa, S., Harris-Coward, P. Y., et al. (2015a). Transcriptome analysis of Aspergillus flavus reveals veA-dependent regulation of secondary metabolite gene clusters, including the novel aflavarin cluster. Eukaryot. Cell 14, 983–997.

PubMed Abstract | Google Scholar

Cary, J. W., Uka, V., Han, Z., Buyst, D., Harris-Coward, P. Y., Ehrlich, K. C., et al. (2015b). An Aspergillus flavus secondary metabolic gene cluster containing a hybrid PKS–NRPS is necessary for synthesis of the 2-pyridones, leporins. Fungal Genet. Biol. 81, 88–97.

PubMed Abstract | Google Scholar

Cary, J. W., Harris-Coward, P. Y., Ehrlich, K. C., Diana, Di Mavungu, J., Malysheva, S. V., et al. (2014). Functional characterization of a veA-dependent polyketide synthase gene in Aspergillus flavus necessary for the synthesis of asparasone, a sclerotium-specific pigment. Fungal Genet. Biol. 64, 25–35. doi: 10.1016/j.fgb.2014.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavaliere, C., Foglia, P., Pastorini, E., Samperi, R., and Laganà, A. (2006). Liquid chromatography/tandem mass spectrometric confirmatory method for determining aflatoxin M1 in cow milk: comparison between electrospray and atmospheric pressure photoionization sources. J. Chromatogr. A 1101, 69–78. doi: 10.1016/j.chroma.2005.09.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Chala, A., Taye, W., Ayalew, A., Krska, R., Sulyok, M., and Logrieco, A. (2014). Multimycotoxin analysis of sorghum (Sorghum bicolor L. Moench) and finger millet (Eleusine coracana L. Garten) from Ethiopia. Food Control 45, 29–35. doi: 10.1016/j.foodcont.2014.04.018

CrossRef Full Text | Google Scholar

Chang, P. K., Horn, B. W., and Dorner, J. W. (2005). Sequence breakpoints in the aflatoxin biosynthesis gene cluster and flanking regions in nonaflatoxigenic Aspergillus flavus isolates. Fungal Genet. Biol. 42, 914–923. doi: 10.1016/j.fgb.2005.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, R. J. (1986). Etiology of Turkey “X” disease in retrospect: a case for the involvement of cyclopiazonic acid. Mycotoxin Res. 2, 3–7. doi: 10.1007/bf03191956

PubMed Abstract | CrossRef Full Text | Google Scholar

Commission regulation (EC) No 1881/2006 (2006). Setting Maximum Levels for Certain Contaminants in Foodstuffs. Available at: https://ec.europa.eu/food/safety/chemical_safety/contaminants/catalogue/aflatoxins_en (accessed December 19, 2006).

Google Scholar

Coppin, E., Debuchy, R., Arnaise, S., and Picard, M. (1997). Mating types and sexual development in filamentous ascomycetes. Microbiol. Mol. Biol. Rev. 61, 411–428.

Google Scholar

Cotty, P. J. (1997). Aflatoxin-producing potential of communities of Aspergillus section Flavi from cotton producing areas in the United States. Mycol. Res. 101, 698–704. doi: 10.1017/s0953756296003139

CrossRef Full Text | Google Scholar

Debuchy, R., and Turgeon, B. (2006). “Mating-type structure, evolution, and function in Euascomycetes,” in The Mycota, Growth, Differentiation, and Sexuality, Vol. 1, eds U. Kües and R. Fischer (Berlin: Springer), 293–323.

Google Scholar

Díaz, R., Ibáñez, M.,Sancho, J. V., and Hernandez, F. (2011). Building an empirical mass spectra library for screening of organic pollutants by ultra-high-pressure liquid chromatography/hybrid quadrupole time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 25, 355–369. doi: 10.1002/rcm.4860

PubMed Abstract | CrossRef Full Text | Google Scholar

Diener, U. L., Cole, R. J., Sanders, T. H., Payne, G. A., Lee, L. S., and Klich, M. A. (1987). Epidemiology of aflatoxin formation by Aspergillus flavus. Annu. Rev. Phytopathol. 25, 249–270. doi: 10.1146/annurev.phyto.25.1.249

CrossRef Full Text | Google Scholar

Dorner, J. W., Cole, R., Connick, W. J., Daigle, D. J., Mcguire, M. R., and Shasha, B. S. (2003). Evaluation of biological control formulations to reduce aflatoxin contamination in peanuts. Biol. Control 26, 318–324. doi: 10.1016/s1049-9644%2802%2900139-1

CrossRef Full Text | Google Scholar

Dunn, G., Newbold, G. T., and Spring, F. S. (1949). Synthesis of flavacol, a metabolic product of Aspergillus flavus. J. Chem. Soc. 2586–2587.

Google Scholar

Dutcher, J. D. (1947). Aspergillic acid: an antibiotic substance produced by Aspergillus flavus I. General properties; formation of desoxy aspergillic acid; structural conclusions. J. Biol. Chem. 171, 321–339.

Google Scholar

Dutcher, J. D. (1958). Aspergillic acid: an antibiotic substance produced by Aspergillus flavus III. The structure of hydroxy aspergillic acid. J. Biol. Chem. 232, 785–796.

Google Scholar

Ehrlich, K. C. (2014). Non-aflatoxigenic Aspergillus flavus to prevent aflatoxin contamination in crops: advantages and limitations. Front. Microbiol. 5:50. doi: 10.3389/fmicb.2014.00050

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrlich, K. C., Chang, P. M., Yu, J., and Cotty, P. J. (2004). Aflatoxin biosynthesis cluster gene cypA is required for G aflatoxin formation. Appl. Environ. Microbiol. 70, 6518–6524. doi: 10.1128/aem.70.11.6518-6524.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrlich, K. C., and Cotty, P. J. (2004). An isolate of Aspergillus flavus used to reduce aflatoxin contamination in cotton seed has a defective polyketide synthase gene. Appl. Microbiol. Biotechnol. 65, 473–478. doi: 10.1007/s00253-004-1670-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrlich, K. C., Kobbeman, K., Montalbano, B. G., and Cotty, P. J. (2007). Aflatoxin-producing Aspergillus species from Thailand. Int. J. Food Microbiol. 114, 153–159. doi: 10.1016/j.ijfoodmicro.2006.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Elimat, T., Figueroa, M., Ehrmann, B. M., Cech, N. B., Pearce, C. J., and Oberlies, N. H. (2013). High-resolution MS, MS/MS, and UV database of fungal secondary metabolites as a dereplication protocol for bioactive natural products. J. Nat. Prod. 76, 1709–1716. doi: 10.1021/np4004307

PubMed Abstract | CrossRef Full Text | Google Scholar

Ezekiel, C. N., Sulyok, M., Somorin, Y., Odutayo, F. I., Nwabekee, S. U., Balogun, A. T., et al. (2016). Mould and mycotoxin exposure assessment of melon and bush mango seeds, two common soup thickeners consumed in Nigeria. Int. J. Food Microbiol. 237, 83–91. doi: 10.1016/j.ijfoodmicro.2016.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Forseth, R. R., Amaike, S., Schwenk, D., Affeldt, K. J., Hoffmeister, D., Schroeder, F. C., et al. (2013). Homologous NRPS-like gene clusters mediate redundant small-molecule biosynthesis in Aspergillus flavus. Angew. Chem. Int. Ed. 52, 1590–1594. doi: 10.1002/anie.201207456

PubMed Abstract | CrossRef Full Text | Google Scholar

Fredenhagen, A., Derrien, C., and Gassmann, E. (2005). An MS/MS library on an ion-trap instrument for efficient dereplication of natural products. Different fragmentation patterns for [M+ H]+ and [M+ Na]+ ions. J. Nat. Prod. 68, 385–291.

Google Scholar

Frisvad, J. C., Hubka, V., Ezekiel, C. N., Hong, S. B., Nováková, A.,Chen, A. J., et al. (2019). Taxonomy of Aspergillus section Flavi and their production of aflatoxins, ochratoxins and other mycotoxins. Stud. Mycol. 93, 1–63. doi: 10.1016/j.simyco.2018.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Frisvad, J. C., and Larsen, T. O. (2015). Chemodiversity in the genus Aspergillus. Appl. Microbiol. Biotechnol. 99, 7859–7877. doi: 10.1007/s00253-015-6839-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Frisvad, J. C., Skouboe, P., and Samson, R. A. (2005). Taxonomic comparison of three different groups of aflatoxin producers and a new efficient producer of aflatoxin B1, sterigmatocystin and 3-O-methylsterigmatocystin, Aspergillus rambellii sp. nov. Syst. Appl. Microbiol. 28, 442–453. doi: 10.1016/j.syapm.2005.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Gallagher, R. T., McCabe, T., Hirotsu, K., Clardy, J., Nicholson, J., and Wilson, B. J. (1980). Aflavinine, a novel indole-mevalonate metabolite from tremorgen-producing Aspergillus flavus species. Tetrahedron Lett. 21, 243–246. doi: 10.1016/s0040-403928002971179-8

CrossRef Full Text | Google Scholar

Gallagher, R. T., and Wilson, B. J. (1979). Aflatrem, the tremorgenic mycotoxin from Aspergillus flavus. Mycopathologia 66, 183–185. doi: 10.1007/bf00683969

PubMed Abstract | CrossRef Full Text | Google Scholar

Geiser, D. M., Dorner, J. W., Horn, B. W., and Taylor, J. W. (2000). The phylogenetics of mycotoxin and sclerotium production in Aspergillus flavus and Aspergillus oryzae. Fungal Genet. Biol. 31, 169–179. doi: 10.1006/fgbi.2000.1215

PubMed Abstract | CrossRef Full Text | Google Scholar

Geiser, D. M., Pitt, J. I., and Taylor, J. W. (1998). Cryptic speciation and recombination in the aflatoxin-producing fungus Aspergillus flavus. Proc. Natl. Acad. Sci. U.S.A. 95, 388–393. doi: 10.1073/pnas.95.1.388

PubMed Abstract | CrossRef Full Text | Google Scholar

Gil Girol, C., Fisch, K. M., Heinekamp, T., Günther, S., Hüttel, W., Piel, J., et al. (2012). Regio-and stereoselective oxidative phenol coupling in Aspergillus niger. Angew. Chem. Int. Ed. 51, 9788–9791. doi: 10.1002/anie.201203603

PubMed Abstract | CrossRef Full Text | Google Scholar

Glass, N. L., and Kaneko, I. (2003). Fatal attraction: nonself recognition and heterokaryon incompatibility in filamentous fungi. Eukaryot. Cell 2, 1–8. doi: 10.1128/ec.2.1.1-8.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonçalves, S. S., Stchigel, A. M., Cano, J. F., Godoy-Martinez, P. C., Colombo, A. L., and Guarro, J. (2012). Aspergillus novoparasiticus: a new clinical species of the section Flavi. Med. Mycol. 50, 152–160. doi: 10.3109/13693786.2011.593564

PubMed Abstract | CrossRef Full Text | Google Scholar

Grove, J. F. (1972a). New metabolic products of Aspergillus flavus. Part I. Asperentin, its methyl ethers, and 5′-hydroxyasperentin. J. Chem. Soc. Perkin 1, 2400–2406.

Google Scholar

Grove, J. F. (1972b). New metabolic products of Aspergillus flavus. Part II. Asperflavin, anhydroasperflavin, and 5,7-dihydroxy-4-methylphthalide. J. Chem. Soc. Perkin 19, 2406–2411.

Google Scholar

Grove, J. F. (1973). New metabolic products of Aspergillus flavus. Part IV. 4′-Hydroxyasperentin and 5′-hydroxyasperentin 8-methyl ether. J. Chem. Soc. Perkin 22, 2704–2706. doi: 10.1039/p19730002704

CrossRef Full Text | Google Scholar

Grubisha, L. C., and Cotty, P. J. (2009). Genetic isolation among sympatric vegetative compatibility groups of the aflatoxin-producing fungus Aspergillus flavus. Mol. Ecol. 19, 269–280.

PubMed Abstract | Google Scholar

Hedayati, M. T., Pasqualotto, A. C., Warn, P. A., Bowyer, P., and Denning, D. W. (2007). Aspergillus flavus: human pathogen, allergen and mycotoxin producer. Microbiology 153, 1677–1692. doi: 10.1099/mic.0.2007/007641-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Heperkan, D., Somuncuoglu, S., Karbancioglu-Güler, F., and Mecik, N. (2012). Natural contamination of cyclopiazonic acid in dried figs and co-occurrence of aflatoxin. Food Control 23, 82–86. doi: 10.1016/j.foodcont.2011.06.015

CrossRef Full Text | Google Scholar

Hesseltine, C. W., Shotwell, O. L., Ellis, J. J., and Stubblefield, R. D. (1966). Aflatoxin formation by Aspergillus flavus. Bacteriol. Rev. 30, 795–805.

Google Scholar

Hesseltine, C. W., Shotwell, O. L., Smith, M., Ellis, J. J., Vandegraft, E., and Shannon, G. (1968). “Production of various aflatoxins by strains of the Aspergillus flavus series,” in Proceedings of the 1st U.S.-Japan Conference on Toxic Micro-Organisms; 1968 October 7-10, (Honolulu, HI: U.S. Department of the Interior), 202–210.

Google Scholar

Hillis, D. M., and Bull, J. J. (1993). An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst. Biol. 42, 182–192. doi: 10.1093/sysbio/42.2.182

CrossRef Full Text | Google Scholar

Holzapfel, C. W., Steyn, P. S., and Purchase, I. F. (1966). Isolation and structure of aflatoxins M1 and M2. Tetrahedron Lett. 7, 2799–2803. doi: 10.1016/s0040-4039

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, B. W., Gell, R. M., Singh, R., Sorensen, R. B., and Carbone, I. (2016). Sexual reproduction in Aspergillus flavus sclerotia: acquisition of novel alleles from soil populations and uniparental mitochondrial inheritance. PLoS One 11:e0146169. doi: 10.1371/journal.pone.0146169

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, B. W., and Greene, R. L. (1995). Vegetative compatibility within populations of Aspergillus flavus, A. parasiticus, and A. tamarii from a peanut field. Mycologia 87, 324–332. doi: 10.1080/00275514.1995.12026537

CrossRef Full Text | Google Scholar

Horn, B. W., Greene, R. L., Sobolev, V. S., Dorner, J. W., Powell, J. H., and Layton, R. C. (1996). Association of morphology and mycotoxin production with vegetative compatibility groups in Aspergillus flavus, A. parasiticus, and A. tamarii. Mycologia 88, 574–587. doi: 10.1080/00275514.1996.12026688

CrossRef Full Text | Google Scholar

Horn, B. W., Moore, G. G., and Carbone, I. (2009a). Sexual reproduction in Aspergillus flavus. Mycologia 101, 423–429.

Google Scholar

Horn, B. W., Ramirez-Prado, J. H., and Carbone, I. (2009b). Sexual reproduction and recombination in the aflatoxin-producing fungus Aspergillus parasiticus. Fungal Genet. Biol. 46, 169–175.

PubMed Abstract | Google Scholar

Horn, B. W., Sorensen, R., Lamb, M., Sobolev, V., Olarte, R., Worthington, C., et al. (2013). Sexual reproduction in Aspergillus flavus sclerotia naturally produced in corn. Phytopathology 104, 75–85.

PubMed Abstract | Google Scholar

Huang, B., Han, Z., Cai, Z., Wu, Y., and Ren, Y. (2010). Simultaneous determination of aflatoxins B1, B2, G1, G2, M1 and M2 in peanuts and their derivative products by ultra-high-performance liquid chromatography–tandem mass spectrometry. Anal. Chim. Acta 662, 62–68. doi: 10.1016/j.aca.2010.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Hüttel, W., and Müller, M. (2007). Regio-and stereoselective intermolecular oxidative phenol coupling in kotanin biosynthesis by Aspergillus niger. ChemBioChem 8, 521–529. doi: 10.1002/cbic.200600434

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoue, T., Nagatomi, Y., Uyama, A., and Mochizuki, N. (2013). Degradation of aflatoxin B1 during the fermentation of alcoholic beverages. Toxins 5, 1219–1229. doi: 10.3390/toxins5071219

PubMed Abstract | CrossRef Full Text | Google Scholar

Klausmeyer, P., McCloud, T. G., Tucker, K. D., Cardellina, J. H., and Shoemaker, R. H. (2005). Aspirochlorine class compounds from Aspergillus flavus inhibit azole-resistant Candida albicans. J. Nat. Prod. 68, 1300–1302.

PubMed Abstract | Google Scholar

Lee, Y. J., and Hagler, W. M. J. (1991). Aflatoxin and cyclopiazonic acid production by Aspergillus flavus isolated from contaminated maize. J. Food Sci. 56, 871–872. doi: 10.1111/j.1365-2621.1991.tb05406.x

CrossRef Full Text | Google Scholar

Lengeler, K. B., Davidson, R. C., D’souza, C., Harashima, T., Shen, W. C., Wang, P., et al. (2000). Signal transduction cascades regulating fungal development and virulence. Microbiol. Mol. Biol. Rev. 64, 746–785. doi: 10.1128/mmbr.64.4.746-785.2000

CrossRef Full Text | Google Scholar

Luk, K. C., Kobbe, B., and Townsend, J. M. (1977). Production of cyclopiazonic acid by Aspergillus flavus Link. Appl. Environ. Microbiol. 33, 211–212.

PubMed Abstract | Google Scholar

Ma, Y. M., Ma, C. C., Li, T., and Wang, J. (2016). A new furan derivative from an endophytic Aspergillus flavus of Cephalotaxus fortunei. Nat. Prod. Res. 30, 79–84.

PubMed Abstract | Google Scholar

Malysheva, S. V., Arroyo-Manzanares, N., Cary, J. W., Ehrlich, K. C., Vanden Bussche, J., Vanhaecke, L., et al. (2014). Identification of novel metabolites from Aspergillus flavus by high resolution and multiple stage mass spectrometry. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 31, 111–120. doi: 10.1080/19440049.2013.859743

PubMed Abstract | CrossRef Full Text | Google Scholar

Masri, M. S., Lundin, R. E., Page, J. R., and Garcia, V. C. (1967). Crystalline aflatoxin M1 from urine and milk. Nature 215, 753–755. doi: 10.1038/215753a0

CrossRef Full Text | Google Scholar

Metzenberg, R. L., and Glass, N. L. (1990). Mating type and mating strategies in Neurospora. Bioessays 12, 53–59. doi: 10.1002/bies.950120202

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, M. A., Schwartz, T., Pickett, B. E., He Sherry, Klem, E. B., Scheuermann, R. H., et al. (2015). A RESTful API for access to phylogenetic tools via the CIPRES science gateway. Evol. Bioinform. 11, 43–48.

PubMed Abstract | Google Scholar

Monacell, J. T., and Carbone, I. (2014). Mobyle SNAP Workbench: a web-based analysis portal for population genetics and evolutionary genomics. Bioinformatics 30, 1488–1490. doi: 10.1093/bioinformatics/btu055

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, G. G., Elliott, J. L., Singh, R., Dorner, J. W., Stone, E. A., Chulze, S. N., et al. (2013). Sexuality generates diversity in the aflatoxin gene cluster: evidence on a global scale. PLoS Pathog. 9:e1003574. doi: 10.1371/journal.ppat.1003574

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, G. G., Mack, B. M., and Beltz, S. B. (2015a). Draft genome sequences of two closely related aflatoxigenic Aspergillus species obtained from the Ivory Coast. Genome Biol. Evol. 8, 729–732.

PubMed Abstract | Google Scholar

Moore, G. G., Mack, B. M., and Beltz, S. B. (2015b). Genomic sequence of the aflatoxigenic filamentous fungus Aspergillus nomius. BMC Genomics 16:551. doi: 10.1186/s12864-015-1719-6

PubMed Abstract | Google Scholar

Moore, G. G., Mack, B. M., Beltz, S. B., and Gilbert, M. K. (2016). Draft genome sequence of an aflatoxigenic Aspergillus species. A. bombycis. Genome Biol. Evol. 8, 3297–3300. doi: 10.1093/gbe/evw238

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, G. G., Singh, R., Horn, B. W., and Carbone, I. (2009). Recombination and lineage-specific gene loss in the aflatoxin gene cluster of Aspergillus flavus. Mol. Ecol. 18, 4870–4887. doi: 10.1111/j.1365-294x.2009.04414.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakazato, M., Morozumi, S., Saito, K., Fujinuma, K., Nishima, T., and Kasai, N. (1990). Interconversion of aflatoxin B1 and aflatoxicol by several fungi. Appl. Environ. Microbiol. 56, 1465–1470.

PubMed Abstract | Google Scholar

Nguyen, M., Ekstrom, A., Li, X., and Yin, Y. (2015). HGT-Finder: a new tool for horizontal gene transfer finding and application to Aspergillus genomes. Toxins 7, 4035–4053. doi: 10.3390/toxins7104035

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicholson, M. J., Koulman, A., Monahan, B. J., Pritchard, B. L., Payne, G. A., and Scott, B. (2009). Identification of two aflatrem biosynthesis gene loci in Aspergillus flavus and metabolic engineering of Penicillium paxilli to elucidate their function. Appl. Environ. Microbiol. 75, 7469–7481. doi: 10.1128/aem.02146-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, K. F., Månsson, M., Rank, C., Frisvad, J. C., and Larsen, T. O. (2011). Dereplication of microbial natural products by LC-DAD-TOFMS. J. Nat. Prod. 74, 2338–2348. doi: 10.1021/np200254t

PubMed Abstract | CrossRef Full Text | Google Scholar

Novas, V. M., and Cabral, D. (2002). Association of mycotoxin and sclerotia production with compatibility groups in Aspergillus flavus from peanut in Argentina. Plant Dis. 86, 215–219. doi: 10.1094/pdis.2002.86.3.215

PubMed Abstract | CrossRef Full Text | Google Scholar

Nozawa, K., Sekita, S., Harada, M., Udagawa, S. I., and Kawai, K. I. (1989). Isolation and structures of two new indoloditerpenes related to aflavinine from a microsclerotium-producing strain of Aspergillus flavus. Chem. Pharm. Bull. 37, 626–630. doi: 10.1248/cpb.37.626

CrossRef Full Text | Google Scholar

Olarte, R. A., Horn, B. W., Dorner, J. W., Monacell, J. T., Singh, R., Stone, E. A., et al. (2011). Effect of sexual recombination on population diversity in aflatoxin production by Aspergillus flavus and evidence for cryptic heterokaryosis. Mol. Ecol. 21, 1453–1476.

PubMed Abstract | Google Scholar

Olarte, R. A., Worthington, C. J., Horn, B. W., Moore, G. G., Singh, R., Monacell, J. T., et al. (2015). Enhanced diversity and aflatoxigenicity in interspecific hybrids of Aspergillus flavus and Aspergillus parasiticus. Mol. Ecol. 24, 1889–1909.

PubMed Abstract | Google Scholar

Papa, K. E. (1973). The parasexual cycle in Aspergillus flavus. Mycologia 65, 1201–1205. doi: 10.1080/00275514.1973.12019538

CrossRef Full Text | Google Scholar

Perry, M. J., Makins, J. F., Adlard, M. W., and Holt, G. (1984). Aspergillic acids produced by mixed cultures of Aspergillus flavus and Aspergillus nidulans. Microbiology 130, 319–323. doi: 10.1099/00221287-130-2-319

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfeiffer, E., Fleck, S. C., and Metzler, M. (2014). Catechol formation: a novel pathway in the metabolism of sterigmatocystin and 11-methoxysterigmatocystin. Chem. Res. Toxicol. 27, 2093–2099. doi: 10.1021/tx500308k

PubMed Abstract | CrossRef Full Text | Google Scholar

Pildain, M. B., Frisvad, J. C., Vaamonde, G., Cabral, D., Varga, J., and Samson, R. A. (2008). Two novel aflatoxin-producing Aspergillus species from Argentinean peanuts. Int. J. Syst. Evol. Microbiol. 58, 725–735. doi: 10.1099/ijs.0.65123-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Plattner, R. D., Bennett, G. A., and Stubblefield, R. D. (1984). Identification of aflatoxins by quadrupole mass spectrometry/mass spectrometry. J. Assoc. Off. Anal. Chem. 67, 734–738.

Google Scholar

Prandini, A., Transini, G., Sigolo, S., Filippi, L., Laporta, M., and Piva, G. (2009). On the occurrence of aflatoxin M1 in milk and dairy products. Food Chem. Toxicol. 47, 984–991. doi: 10.1016/j.fct.2007.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, E. W., and Carbone, I. (2005). SNAP: workbench management tool for evolutionary population genetic analysis. Bioinformatics 21, 402–404. doi: 10.1093/bioinformatics/bti003

PubMed Abstract | CrossRef Full Text | Google Scholar

Probst, C., Bandyopadhyay, R., Price, L. E., and Cotty, P. J. (2011). Identification of atoxigenic Aspergillus flavus isolates to reduce aflatoxin contamination of maize in Kenya. Plant Dis. 95, 212–218. doi: 10.1094/pdis-06-10-0438

PubMed Abstract | CrossRef Full Text | Google Scholar

Probst, C., Callicott, K. A., and Cotty, P. J. (2012). Deadly strains of Kenyan Aspergillus are distinct from other aflatoxin producers. Eur. J. Plant Pathol. 132, 419–429. doi: 10.1007/s10658-011-9887-y

CrossRef Full Text | Google Scholar

Ramirez-Prado, J. H., Moore, G. G., Horn, B. W., and Carbone, I. (2008). Characterization and population analysis of the mating-type genes in Aspergillus flavus and Aspergillus parasiticus. Fungal Genet. Biol. 45, 1292–1299. doi: 10.1016/j.fgb.2008.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Rank, C., Nielsen, K. F., Larsen, T. O., Varga, J., Samson, R. A., and Frisvad, J. C. (2011). Distribution of sterigmatocystin in filamentous fungi. Fungal Biol. 115, 406–420. doi: 10.1016/j.funbio.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodricks, J. V., Lustig, E., Campbell, A. D., Stoloff, L., and Henery-Logan, K. R. (1968). Aspertoxin, a hydroxy derivative of O-methylsterigmatocystin from aflatoxin-producing cultures of Aspergillus flavus. Tetrahedron Lett. 9, 2975–2978. doi: 10.1016/s0040-4039

PubMed Abstract | CrossRef Full Text | Google Scholar

Rushing, B. R., and Selim, M. I. (2017). Structure and oxidation of pyrrole adducts formed between aflatoxin B2a and biological amines. Chem. Res. Toxicol. 30, 1275–1285. doi: 10.1021/acs.chemrestox.7b00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez, J. F., Chiang, Y. M., Szewczyk, E., Davidson, A. D., Ahuja, M., Oakley, C. E., et al. (2010). Molecular genetic analysis of the orsellinic acid/F9775 gene cluster of Aspergillus nidulans wz. Mol. Biosyst. 6, 587–593. doi: 10.1039/b904541d

PubMed Abstract | CrossRef Full Text | Google Scholar

Saruwatari, T., Yagishita, F., Mino, T., Noguchi, H., Hotta, K., and Watanabe, K. (2014). Cytochrome P450 as dimerization catalyst in diketopiperazine alkaloid biosynthesis. ChemBioChem 15, 656–659. doi: 10.1002/cbic.201300751

PubMed Abstract | CrossRef Full Text | Google Scholar

Seshime, Y., Juvvadi, P. R., Tokuoka, M., Koyama, Y., Kitamoto, K., Ebizuka, Y., et al. (2009). Functional expression of the Aspergillus flavus PKS-NRPS hybrid CpaA involved in the biosynthesis of cyclopiazonic acid. Bioorg. Med. Chem. Lett. 19, 3288–3292. doi: 10.1016/j.bmcl.2009.04.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiomi, K., Hatae, K., Yamaguchi, Y., Masuma, R., Tomoda, H., Kobayashi, S., et al. (2002). New antibiotics miyakamides produced by a fungus. J. Antibiot. 55, 952–961. doi: 10.7164/antibiotics.55.952

PubMed Abstract | CrossRef Full Text | Google Scholar

Soares, C., Rodrigues, P., Peterson, S. W., Lima, N., and Venâncio, A. (2012). Three new species of Aspergillus section Flavi isolated from almonds and maize in Portugal. Mycologia 104, 682–697. doi: 10.3852/11-088

PubMed Abstract | CrossRef Full Text | Google Scholar

Springer, J. P., Bűchi, G., Kobbe, B., Demain, A. L., and Clardy, J. (1977). The structure of ditryptophenaline- a new metabolite of Aspergillus flavus. Tetrahedron Lett. 18, 2403–2406. doi: 10.1016/s0040-4039

CrossRef Full Text | Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Stubblefield, R. D., Shannon, G. M., and Shotwell, O. L. (1970). Aflatoxins M1 and M2: preparation and purification. J. Am. Oil Chem. Soc. 47, 389–390. doi: 10.1007/bf02632471

PubMed Abstract | CrossRef Full Text | Google Scholar

Sulyok, M., Beed, F., Boni, S., Abass, A., Mukunzi, A., and Krska, R. (2015). Quantitation of multiple mycotoxins and cyanogenic glucosides in cassava samples from Tanzania and Rwanda by an LC-MS/MS-based multi-toxin method. Food Addit. Contam. Part A. Chem. Anal. Control Expo. Risk Assess. 32, 488–502. doi: 10.1080/19440049.2014.975752

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, K., Li, Y., Guo, L., Wang, Y., Liu, P., and Zhu, W. (2014). Indole-diterpenoids and isocoumarin from the fungus, Aspergillus flavus, isolated from the prawn, Penaeus vannamei. Mar. Drugs 12, 3970–3981. doi: 10.3390/md12073970

PubMed Abstract | CrossRef Full Text | Google Scholar

Swofford, D. L. (1998). PAUP Phylogenetic Analysis Using Parsimony (and other Methods). Sunderland, MA: Sinauer Associates.

Google Scholar

Taber, R. A., and Schroeder, H. W. (1967). Aflatoxin-producing potential of isolates of the Aspergillus flavus-oryzae group from peanuts (Arachis hypogaea). Appl. Microbiol. 15, 140–144.

PubMed Abstract | Google Scholar

Tepaske, M., Gloer, J. B., Wicklow, D. T., and Dowd, P. F. (1990). Aflavazole: a new antiinsectan carbazole metabolite from the sclerotia of Aspergillus flavus. J. Org. Chem. 55, 5299–5301. doi: 10.1021/jo00305a030

CrossRef Full Text | Google Scholar

Tokuoka, M., Koyama, Y., Kitamoto, K., and Ebizuka, Y. (2009). Functional expression of the Aspergillus flavus PKS–NRPS hybrid CpaA involved in the biosynthesis of cyclopiazonic acid. Bioorg. Med. Chem. Lett. 19, 3288–3292. doi: 10.1016/j.bmcl.2009.04.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Turgeon, B. G., and Yoder, O. C. (2000). Proposed nomenclature for mating type genes of filamentous ascomycetes. Fungal Genet. Biol. 31, 1–5. doi: 10.1006/fgbi.2000.1227

PubMed Abstract | CrossRef Full Text | Google Scholar

Turner, P. C., Moore, S. E., Hall, A. J., Prentice, A. M., and Wild, C. P. (2003). Modification of immune function through exposure to dietary aflatoxin in Gambian children. Environ. Health Perspect. 111, 217–220. doi: 10.1289/ehp.5753

PubMed Abstract | CrossRef Full Text | Google Scholar

Uka, V., Moore, G. G., Arroyo-Manzanares, N., Nebija, D., De Saeger, S., Diana, et al. (2017). Unravelling the diversity of CPA family of mycotoxins in Aspergillus flavus by UHPLC Triple-TOF HRMS. Toxins 9, 1–21.

Google Scholar

Umemura, M., Nagano, N., Koike, H., Kawano, J., Ishii, T., Miyamura, Y., et al. (2014). Characterization of the biosynthetic gene cluster for the ribosomally synthesized cyclic peptide ustiloxin B in Aspergillus flavus. Fungal Genet. Biol. 68, 23–30. doi: 10.1016/j.fgb.2014.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Urano, T., Trucksess, M. W., Beaver, R. W., Wilson, D. M., Dorner, J. W., and Dowell, F. E. (1992). Co-occurrence of cyclopiazonic acid and aflatoxins in corn and peanuts. J. Off. Anal. Chem. Int. 75, 838–841.

PubMed Abstract | Google Scholar

Varga, E., Glauner, T., Berthiller, F., Krska, R., Schuhmacher, R., and Sulyok, M. (2013). Development and validation of a (semi-) quantitative UHPLC-MS/MS method for the determination of 191 mycotoxins and other fungal metabolites in almonds, hazelnuts, peanuts and pistachios. Anal. Bioanal. Chem. 405, 5087–5104. doi: 10.1007/s00216-013-6831-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Varga, J., Frisvad, J. C., and Samson, R. A. (2011). Two new aflatoxin producing species, and an overview of Aspergillus section Flavi. Stud. Mycol. 69, 57–80. doi: 10.3114/sim.2011.69.05

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfender, J. L., Ndjoko, K., and Hostettmann, K. (2003). Liquid chromatography with ultraviolet absorbance–mass spectrometric detection and with nuclear magnetic resonance spectrometry: a powerful combination for the on-line structural investigation of plant metabolites. J. Chromatogr. A 1000, 437–455. doi: 10.1016/s0021-9673

CrossRef Full Text | Google Scholar

Wolfender, J. L., Waridel, P., Ndjoko, K., Hobby, K. R., Major, H. J., and Hostettmann, K. (2000). Evaluation of Q-TOF-MS/MS and multiple stage IT-MSn for the dereplication of flavonoids and related compounds in crude plant extracts. Analusis 28, 895–906. doi: 10.1051/analusis%3A2000280895

CrossRef Full Text | Google Scholar

Yabe, K., Chihaya, N., Hatabayashi, H., Kito, M., Hoshino, S., Zeng, H., et al. (2012). Production of M-/GM-group aflatoxins catalyzed by the OrdA enzyme in aflatoxin biosynthesis. Fungal Genet. Biol. 49, 744–754. doi: 10.1016/j.fgb.2012.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J., Chang, P. K., Ehrlich, K. C., Cary, J. W., Bhatnagar, D., Cleveland, T. E., et al. (2004). Clustered pathway genes in aflatoxin biosynthesis. Appl. Environ. Microbiol. 70, 1253–1262. doi: 10.1128/aem.70.3.1253-1262.2004

CrossRef Full Text | Google Scholar

Zhu, F., Wu, J., Chen, G., Lu, W., and Pan, J. (2011). Biosynthesis, characterization and biological evaluation of Fe (III) and Cu (II) complexes of neoaspergillic acid, a hydroxamate siderophore produced by co-cultures of two marine-derived mangrove epiphytic fungi. Nat. Prod. Commun. 6, 1137–1140.

Google Scholar

Zorzete, P., Baquiao, A. C., Atayde, D. D., Reis, T. A., Goncalez, E., and Correa, B. (2013). Mycobiota, aflatoxins and cyclopiazonic acid in stored peanut cultivars. Food Res. Int. 52, 380–386. doi: 10.1016/j.foodres.2013.03.029

CrossRef Full Text | Google Scholar

Keywords: Aspergillus, sexuality, chemical diversity, genetic diversity, mass spectrometry

Citation: Uka V, Moore GG, Arroyo-Manzanares N, Nebija D, De Saeger S and Diana Di Mavungu J (2019) Secondary Metabolite Dereplication and Phylogenetic Analysis Identify Various Emerging Mycotoxins and Reveal the High Intra-Species Diversity in Aspergillus flavus. Front. Microbiol. 10:667. doi: 10.3389/fmicb.2019.00667

Received: 06 November 2018; Accepted: 18 March 2019;
Published: 05 April 2019.

Edited by:

Isabelle P. Oswald, Institut National de la Recherche Agronomique de Toulouse, France

Reviewed by:

Jéssica Gil-Serna, Complutense University of Madrid, Spain
Masoomeh Shams-Ghahfarokhi, Tarbiat Modares University, Iran

Copyright © 2019 Uka, Moore, Arroyo-Manzanares, Nebija, De Saeger and Diana Di Mavungu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: José Diana Di Mavungu, jose.dianadimavungu@ugent.be

For the purposes of this publication, all references, including in the bibliography, to “Kosovo” should be understood/read as “Kosovo [in accordance with Security Council Resolution 1244 (1999)]”

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.