Population Subdivision and the Frequency of Aflatoxigenic Isolates in Aspergillus flavus in the United States

Microevolution in the pansecondary metabolome of Aspergillus flavus and its potential macroevolutionary implications for filamentous fungi

Fungi produce a wealth of pharmacologically bioactive secondary metabolites (SMs) from biosynthetic gene clusters (BGCs). It is common practice for drug discovery efforts to treat species' secondary metabolomes as being well represented by a single or a small number of representative genomes. However, this approach misses the possibility that intraspecific population dynamics, such as adaptation to environmental conditions or local microbiomes, may harbor novel BGCs that contribute to the overall niche breadth of species. Using 94 isolates of Aspergillus flavus, a cosmopolitan model fungus, sampled from seven states in the United States, we dereplicate 7,821 BGCs into 92 unique BGCs. We find that more than 25% of pangenomic BGCs show population-specific patterns of presence/absence or protein divergence. Population-specific BGCs make up most of the accessory-genome BGCs, suggesting that different ecological forces that maintain accessory genomes may be partially mediated by population-specific differences in secondary metabolism. We use ultra-high-performance high-resolution mass spectrometry to confirm that these genetic differences in BGCs also result in chemotypic differences in SM production in different populations, which could mediate ecological interactions and be acted on by selection. Thus, our results suggest a paradigm shift that previously unrealized population-level reservoirs of SM diversity may be of significant evolutionary, ecological, and pharmacological importance. Last, we find that several population-specific BGCs from A. flavus are present in Aspergillus parasiticus and Aspergillus minisclerotigenes and discuss how the microevolutionary patterns we uncover inform macroevolutionary inferences and help to align fungal secondary metabolism with existing evolutionary theory.

Aflasafe SN01 is the First Biocontrol Product Approved for Aflatoxin Mitigation in Two Nations, Senegal and The Gambia

Aflatoxin contamination is caused by Aspergillus flavus and closely related fungi. In The Gambia, aflatoxin contamination of groundnut and maize, two staple and economically important crops, is common. Groundnut and maize consumers are chronically exposed to aflatoxins, sometimes at alarming levels, and this has severe consequences on their health and productivity. Aflatoxin contamination also impedes commercialization in local and international premium markets. In neighboring Senegal, an aflatoxin biocontrol product containing four atoxigenic isolates of A. flavus, Aflasafe SN01, has been registered and is approved for commercial use in groundnut and maize. We detected that the four genotypes composing Aflasafe SN01 are also native to The Gambia. The biocontrol product was tested during two years in 129 maize and groundnut fields and compared with corresponding untreated fields cropped by smallholder farmers in The Gambia. Treated crops contained up to 100% less aflatoxins than untreated crops. A large portion of the crops could have been commercialized in premium markets due to the low aflatoxin content (in many cases no detectable aflatoxins), both at harvest and after storage. Substantial aflatoxin reductions were also achieved when commercially produced groundnut received treatment. Here we report for the first time the use and effectiveness of an aflatoxin biocontrol product registered for use in two nations. With the current scale-out and -up efforts of Aflasafe SN01, a large number of farmers, consumers, and traders in The Gambia and Senegal will obtain health, income, and trade benefits.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

The Frequency of Sex: Population Genomics Reveals Differences in Recombination and Population Structure of the Aflatoxin-Producing Fungus Aspergillus flavus

The apparent rarity of sex in many fungal species has raised questions about how much sex is needed to purge deleterious mutations and how differences in frequency of sex impact fungal evolution. We sought to determine how differences in the extent of recombination between populations of Aspergillus flavus impact the evolution of genes associated with the synthesis of aflatoxin, a notoriously potent carcinogen. We sequenced the genomes of, and quantified aflatoxin production in, 94 isolates of A. flavus sampled from seven states in eastern and central latitudinal transects of the United States. The overall population is subdivided into three genetically differentiated populations (A, B, and C) that differ greatly in their extent of recombination, diversity, and aflatoxin-producing ability. Estimates of the number of recombination events and linkage disequilibrium decay suggest relatively frequent sex only in population A. Population B is sympatric with population A but produces significantly less aflatoxin and is the only population where the inability of nonaflatoxigenic isolates to produce aflatoxin was explained by multiple gene deletions. Population expansion evident in population B suggests a recent introduction or range expansion. Population C is largely nonaflatoxigenic and restricted mainly to northern sampling locations through restricted migration and/or selection. Despite differences in the number and type of mutations in the aflatoxin gene cluster, codon optimization and site frequency differences in synonymous and nonsynonymous mutations suggest that low levels of recombination in some A. flavus populations are sufficient to purge deleterious mutations.IMPORTANCE Differences in the relative frequencies of sexual and asexual reproduction have profound implications for the accumulation of deleterious mutations (Muller's ratchet), but little is known about how these differences impact the evolution of ecologically important phenotypes. Aspergillus flavus is the main producer of aflatoxin, a notoriously potent carcinogen that often contaminates food. We investigated if differences in the levels of production of aflatoxin by A. flavus could be explained by the accumulation of deleterious mutations due to a lack of recombination. Despite differences in the extent of recombination, variation in aflatoxin production is better explained by the demography and history of specific populations and may suggest important differences in the ecological roles of aflatoxin among populations. Furthermore, the association of aflatoxin production and populations provides a means of predicting the risk of aflatoxin contamination by determining the frequencies of isolates from low- and high-production populations.

A Liquid Chromatographic Method for Rapid and Sensitive Analysis of Aflatoxins in Laboratory Fungal Cultures

Culture methods supplemented with high-performance liquid chromatography (HPLC) technique provide a rapid and simple tool for detecting levels of aflatoxins (AFs) produced by fungi. This study presents a robust method for simultaneous quantification of aflatoxin (AF) B1, B2, G1, and G2 levels in several fungal cultivation states: submerged shake culture, liquid slant culture, and solid-state culture. The recovery of the method was evaluated by spiking a mixture of AFs at several concentrations to the test medium. The applicability of the method was evaluated by using aflatoxigenic and non-aflatoxigenic Aspergilli. A HPLC coupled with the diode array (DAD) and fluorescence (FLD) detectors was used to determine the presence and amounts of AFs. Both detectors showed high sensitivity in detecting spiked AFs or AFs produced in situ by toxigenic fungi. Our methods showed 76%–88% recovery from medium spiked with 2.5, 10, 50, 100, and 500 ng/mL AFs. The limit of quantification (LOQ) for AFs were 2.5 to 5.0 ng/mL with DAD and 0.025 to 2.5 ng/mL with FLD. In this work, we described in detail a protocol, which can be considered the foremost and only verified method, to extract, detect, and quantify AFs employing both aflatoxigenic and non-toxigenic Aspergilli.

Biocontrol Strains Differentially Shift the Genetic Structure of Indigenous Soil Populations of Aspergillus flavus

Biocontrol using non-aflatoxigenic strains of Aspergillus flavus has the greatest potential to mitigate aflatoxin contamination in agricultural produce. However, factors that influence the efficacy of biocontrol agents in reducing aflatoxin accumulation under field conditions are not well-understood. Shifts in the genetic structure of indigenous soil populations of A. flavus following application of biocontrol products Afla-Guard and AF36 were investigated to determine how these changes can influence the efficacy of biocontrol strains in reducing aflatoxin contamination. Soil samples were collected from maize fields in Alabama, Georgia, and North Carolina in 2012 and 2013 to determine changes in the population genetic structure of A. flavus in the soil following application of the biocontrol strains. A. flavus L was the most dominant species of Aspergillus section Flavi with a frequency ranging from 61 to 100%, followed by Aspergillus parasiticus that had a frequency of <35%. The frequency of A. flavus L increased, while that of A. parasiticus decreased after application of biocontrol strains. A total of 112 multilocus haplotypes (MLHs) were inferred from 1,282 isolates of A. flavus L using multilocus sequence typing of the trpC, mfs, and AF17 loci. A. flavus individuals belonging to the Afla-Guard MLH in the IB lineage were the most dominant before and after application of biocontrol strains, while individuals of the AF36 MLH in the IC lineage were either recovered in very low frequencies or not recovered at harvest. There were no significant (P > 0.05) differences in the frequency of individuals with MAT1-1 and MAT1-2 for clone-corrected MLH data, an indication of a recombining population resulting from sexual reproduction. Population mean mutation rates were not different across temporal and spatial scales indicating that mutation alone is not a driving force in observed multilocus sequence diversity. Clustering based on principal component analysis identified two distinct evolutionary lineages (IB and IC) across all three states. Additionally, patristic distance analysis revealed phylogenetic incongruency among single locus phylogenies which suggests ongoing genetic exchange and recombination. Levels of aflatoxin accumulation were very low except in North Carolina in 2012, where aflatoxin levels were significantly (P < 0.05) lower in grain from treated compared to untreated plots. Phylogenetic analysis showed that Afla-Guard was more effective than AF36 in shifting the indigenous soil populations of A. flavus toward the non-toxigenic or low aflatoxin producing IB lineage. These results suggest that Afla-Guard, which matches the genetic and ecological structure of indigenous soil populations of A. flavus in Alabama, Georgia, and North Carolina, is likely to be more effective in reducing aflatoxin accumulation and will also persist longer in the soil than AF36 in the southeastern United States.

Fitness Cost of Aflatoxin Production in Aspergillus flavus When Competing with Soil Microbes Could Maintain Balancing Selection

Selective forces that maintain the polymorphism for aflatoxigenic and nonaflatoxigenic individuals of Aspergillus flavus are largely unknown. As soils are widely considered the natural habitat of A. flavus, we hypothesized that aflatoxin production would confer a fitness advantage in the soil environment. To test this hypothesis, we used A. flavus DNA quantified by quantitative PCR (qPCR) as a proxy for fitness of aflatoxigenic and nonaflatoxigenic field isolates grown in soil microcosms. Contrary to predictions, aflatoxigenic isolates had significantly lower fitness than did nonaflatoxigenic isolates in natural soils across three temperatures (25, 37, and 42°C). The addition of aflatoxin to soils (500 ng/g) had no effect on the growth of A. flavus Amplicon sequencing showed that neither the aflatoxin-producing ability of the fungus nor the addition of aflatoxin had a significant effect on the composition of fungal or bacterial communities in soil. We argue that the fitness disadvantage of aflatoxigenic isolates is most likely explained by the metabolic cost of producing aflatoxin. Coupled with a previous report of a selective advantage of aflatoxin production in the presence of some insects, our findings give an ecological explanation for balancing selection resulting in persistent polymorphisms in aflatoxin production.IMPORTANCE Aflatoxin, produced by the fungus Aspergillus flavus, is an extremely potent hepatotoxin that causes acute toxicosis and cancer, and it incurs hundreds of millions of dollars annually in agricultural losses. Despite the importance of this toxin to humans, it has remained unclear what the fungus gains by producing aflatoxin. In fact, not all strains of A. flavus produce aflatoxin. Previous work has shown an advantage to producing aflatoxin in the presence of some insects. Our current work demonstrates the first evidence of a disadvantage to A. flavus in producing aflatoxin when competing with soil microbes. Together, these opposing evolutionary forces could explain the persistence of both aflatoxigenic and nonaflatoxigenic strains through evolutionary time.