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.

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.