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

Divergent Aspergillus flavus corn population is composed of prolific conidium producers: Implications for saprophytic disease cycle

The ascomycete fungus Aspergillus flavus infects and contaminates corn, peanuts, cottonseed, and tree nuts with toxic and carcinogenic aflatoxins. Subdivision between soil and host plant populations suggests that certain A. flavus strains are specialized to infect peanut, cotton, and corn despite having a broad host range. In this study, the ability of strains isolated from corn and/or soil in 11 Louisiana fields to produce conidia (field inoculum and male gamete) and sclerotia (resting bodies and female gamete) was assessed and compared with genotypic single-nucleotide polymorphism (SNP) differences between whole genomes. Corn strains produced upward of 47× more conidia than strains restricted to soil. Conversely, corn strains produced as much as 3000× fewer sclerotia than soil strains. Aspergillus flavus strains, typified by sclerotium diameter (small S-strains, <400 μm; large L-strains, >400 μm), belonged to separate clades. Several strains produced a mixture (M) of S and L sclerotia, and an intermediate number of conidia and sclerotia, compared with typical S-strains (minimal conidia, copious sclerotia) and L-strains (copious conidia, minimal sclerotia). They also belonged to a unique phylogenetic mixed (M) clade. Migration from soil to corn positively correlated with conidium production and negatively correlated with sclerotium production. Genetic differences correlated with differences in conidium and sclerotium production. Opposite skews in female (sclerotia) or male (conidia) gametic production by soil or corn strains, respectively, resulted in reduced effective breeding population sizes when comparing male:female gamete ratio with mating type distribution. Combining both soil and corn populations increased the effective breeding population, presumably due to contribution of male gametes from corn, which fertilize sclerotia on the soil surface. Incongruencies between aflatoxin clusters, strain morphotype designation, and whole genome phylogenies suggest a history of sexual reproduction within this Louisiana population, demonstrating the importance of conidium production, as infectious propagules and as fertilizers of the A. flavus soil population.

Reexamination of Aspergillus cristatus phylogeny in dark tea: Characteristics of the mitochondrial genome

To enhance our understanding of Aspergillus cristatus, an important functional microorganism, the characteristics of its mitochondrial genome were analyzed and compared with related species. The mitochondrial genome of A. cristatus was determined to be 77,649 bp in length, with 15 protein-coding regions. Notably, its length surpassed that of the other species, primarily attributable to the intron length. Gene order exhibited significant variations, with greater conservation observed in the genus Penicillium compared to Aspergillus. Phylogenetic tree analyses indicated that the genera Aspergillus and Penicillium are closely related but monophyletic. Furthermore, the phylogenetic tree constructed based on protein-coding genes effectively distinguished all strains with high branching confidence. This approach provides a robust reflection of the evolutionary relationship between A. cristatus and its related species, offering potential for the development of molecular markers suitable for Aspergillus and Penicillium.

Pangenomics of the death cap mushroom Amanita phalloides, and of Agaricales, reveals dynamic evolution of toxin genes in an invasive range

The poisonous European mushroom Amanita phalloides (the "death cap") is invading California. Whether the death caps' toxic secondary metabolites are evolving as it invades is unknown. We developed a bioinformatic pipeline to identify the MSDIN genes underpinning toxicity and probed 88 death cap genomes from an invasive Californian population and from the European range, discovering a previously unsuspected diversity of MSDINs made up of both core and accessory elements. Each death cap individual possesses a unique suite of MSDINs, and toxin genes are significantly differentiated between Californian and European samples. MSDIN genes are maintained by strong natural selection, and chemical profiling confirms MSDIN genes are expressed and result in distinct phenotypes; our chemical profiling also identified a new MSDIN peptide. Toxin genes are physically clustered within genomes. We contextualize our discoveries by probing for MSDINs in genomes from across the order Agaricales, revealing MSDIN diversity originated in independent gene family expansions among genera. We also report the discovery of an MSDIN in an Amanita outside the "lethal Amanitas" clade. Finally, the identification of an MSDIN gene and its associated processing gene (POPB) in Clavaria fumosa suggest the origin of MSDINs is older than previously suspected. The dynamic evolution of MSDINs underscores their potential to mediate ecological interactions, implicating MSDINs in the ongoing invasion. Our data change the understanding of the evolutionary history of poisonous mushrooms, emphasizing striking parallels to convergently evolved animal toxins. Our pipeline provides a roadmap for exploring secondary metabolites in other basidiomycetes and will enable drug prospecting.

Assembling Quality Genomes of Flax Fungal Pathogens from Oxford Nanopore Technologies Data

Flax (Linum usitatissimum L.) is attacked by numerous devastating fungal pathogens, including Colletotrichum lini, Aureobasidium pullulans, and Fusarium verticillioides (Fusarium moniliforme). The effective control of flax diseases follows the paradigm of extensive molecular research on pathogenicity. However, such studies require quality genome sequences of the studied organisms. This article reports on the approaches to assembling a high-quality fungal genome from the Oxford Nanopore Technologies data. We sequenced the genomes of C. lini, A. pullulans, and F. verticillioides (F. moniliforme) and received different volumes of sequencing data: 1.7 Gb, 3.9 Gb, and 11.1 Gb, respectively. To obtain the optimal genome sequences, we studied the effect of input data quality and genome coverage on assembly statistics and tested the performance of different assembling and polishing software. For C. lini, the most contiguous and complete assembly was obtained by the Flye assembler and the Homopolish polisher. The genome coverage had more effect than data quality on assembly statistics, likely due to the relatively low amount of sequencing data obtained for C. lini. The final assembly was 53.4 Mb long and 96.4% complete (according to the glomerellales_odb10 BUSCO dataset), consisted of 42 contigs, and had an N50 of 4.4 Mb. For A. pullulans and F. verticillioides (F. moniliforme), the best assemblies were produced by Canu–Medaka and Canu–Homopolish, respectively. The final assembly of A. pullulans had a length of 29.5 Mb, 99.4% completeness (dothideomycetes_odb10), an N50 of 2.4 Mb and consisted of 32 contigs. F. verticillioides (F. moniliforme) assembly was 44.1 Mb long, 97.8% complete (hypocreales_odb10), consisted of 54 contigs, and had an N50 of 4.4 Mb. The obtained results can serve as a guideline for assembling a de novo genome of a fungus. In addition, our data can be used in genomic studies of fungal pathogens or plant–pathogen interactions and assist in the management of flax diseases.

Advanced mycotoxin control and decontamination techniques in view of an increased aflatoxin risk in Europe due to climate change

Aflatoxins are toxic secondary metabolites produced by Aspergillus spp. found in staple food and feed commodities worldwide. Aflatoxins are carcinogenic, teratogenic, and mutagenic, and pose a serious threat to the health of both humans and animals. The global economy and trade are significantly affected as well. Various models and datasets related to aflatoxins in maize have been developed and used but have not yet been linked. The prevention of crop loss due to aflatoxin contamination is complex and challenging. Hence, the set-up of advanced decontamination is crucial to cope with the challenge of climate change, growing population, unstable political scenarios, and food security problems also in European countries. After harvest, decontamination methods can be applied during transport, storage, or processing, but their application for aflatoxin reduction is still limited. Therefore, this review aims to investigate the effects of environmental factors on aflatoxin production because of climate change and to critically discuss the present-day and novel decontamination techniques to unravel gaps and limitations to propose them as a tool to tackle an increased aflatoxin risk in Europe.

Asymmetrical lineage introgression and recombination in populations of Aspergillus flavus: Implications for biological control

Aspergillus flavus is an agriculturally important fungus that causes ear rot of maize and produces aflatoxins, of which B1 is the most carcinogenic naturally-produced compound. In the US, the management of aflatoxins includes the deployment of biological control agents that comprise two nonaflatoxigenic A. flavus strains, either Afla-Guard (member of lineage IB) or AF36 (lineage IC). We used genotyping-by-sequencing to examine the influence of both biocontrol agents on native populations of A. flavus in cornfields in Texas, North Carolina, Arkansas, and Indiana. This study examined up to 27,529 single-nucleotide polymorphisms (SNPs) in a total of 815 A. flavus isolates, and 353 genome-wide haplotypes sampled before biocontrol application, three months after biocontrol application, and up to three years after initial application. Here, we report that the two distinct A. flavus evolutionary lineages IB and IC differ significantly in their frequency distributions across states. We provide evidence of increased unidirectional gene flow from lineage IB into IC, inferred to be due to the applied Afla-Guard biocontrol strain. Genetic exchange and recombination of biocontrol strains with native strains was detected in as little as three months after biocontrol application and up to one and three years later. There was limited inter-lineage migration in the untreated fields. These findings suggest that biocontrol products that include strains from lineage IB offer the greatest potential for sustained reductions in aflatoxin levels over several years. This knowledge has important implications for developing new biocontrol strategies.

Recombination-aware phylogeographic inference using the structured coalescent with ancestral recombination

Movement of individuals between populations or demes is often restricted, especially between geographically isolated populations. The structured coalescent provides an elegant theoretical framework for describing how movement between populations shapes the genealogical history of sampled individuals and thereby structures genetic variation within and between populations. However, in the presence of recombination an individual may inherit different regions of their genome from different parents, resulting in a mosaic of genealogical histories across the genome, which can be represented by an Ancestral Recombination Graph (ARG). In this case, different genomic regions may have different ancestral histories and so different histories of movement between populations. Recombination therefore poses an additional challenge to phylogeographic methods that aim to reconstruct the movement of individuals from genealogies, although also a potential benefit in that different loci may contain additional information about movement. Here, we introduce the Structured Coalescent with Ancestral Recombination (SCAR) model, which builds on recent approximations to the structured coalescent by incorporating recombination into the ancestry of sampled individuals. The SCAR model allows us to infer how the migration history of sampled individuals varies across the genome from ARGs, and improves estimation of key population genetic parameters such as population sizes, recombination rates and migration rates. Using the SCAR model, we explore the potential and limitations of phylogeographic inference using full ARGs. We then apply the SCAR to lineages of the recombining fungus Aspergillus flavus sampled across the United States to explore patterns of recombination and migration across the genome.

Phylogeographic methods are widely used to reconstruct the historical movement of individuals between different populations. When applied to infectious pathogens, these methods are often used to reconstruct the origin or source of novel pathogen lineages. Most existing phylogeographic methods reconstruct movement based on a single phylogenetic tree, which is assumed to reflect the genetic ancestry of all sampled individuals. However in populations undergoing recombination, genetic material can be exchanged between lineages such that individuals may inherit different regions of their genome from different ancestors. In this case, phylogenetic relationships among individuals can only be captured by a reticulated network rather than any single tree. Ancestral Recombination Graphs (ARGs) provide one way of capturing these reticulate relationships and we develop new models that allow for demographic inference of historical population sizes, recombination rates and migration rates between subpopulations from ARGs. By accounting for recombination, our models not only allow for accurate demographic inference, but can take full advantage of the additional information contained in ARGs about how ancestry varies across genomes to more precisely reconstruct the movement of genetic material between populations.

Post-translational modifications drive secondary metabolite biosynthesis in Aspergillus: a review

Post‐translational modifications (PTMs) are important for protein function and regulate multiple cellular processes and secondary metabolites (SMs) in fungi. Aspergillus species belong to a genus renown for an abundance of bioactive secondary metabolites, many important as toxins, pharmaceuticals and in industrial production. The genes required for secondary metabolites are typically co‐localized in biosynthetic gene clusters (BGCs), which often localize in heterochromatic regions of genome and are ‘turned off’ under laboratory condition. Efforts have been made to ‘turn on’ these BGCs by genetic manipulation of histone modifications, which could convert the heterochromatic structure to euchromatin. Additionally, non‐histone PTMs also play critical roles in the regulation of secondary metabolism. In this review, we collate the known roles of epigenetic and PTMs on Aspergillus SM production. We also summarize the proteomics approaches and bioinformatics tools for PTM identification and prediction and provide future perspectives on the emerging roles of PTM on regulation of SM biosynthesis in Aspergillus and other fungi.

Dysfunction of FadA-cAMP signalling decreases Aspergillus flavus resistance to antimicrobial natural preservative Perillaldehyde and AFB1 biosynthesis

Aspergillus flavus is an opportunistic fungal pathogen that colonizes agriculture crops with aflatoxin contamination. We found that Perillaldehyde (PAE) effectively inhibited A. flavus viability and aflatoxin production by inducing excess reactive oxygen species (ROS). Transcriptome analysis indicated that the Gα protein FadA was significantly induced by PAE. Functional characterization of FadA showed it is important for asexual development and aflatoxin biosynthesis by regulation of cAMP-PKA signalling. The ΔfadA mutant was more sensitive to PAE, while ΔpdeL and ΔpdeH mutants can tolerate excess PAE compared to wild-type A. flavus. Further RNA-sequence analysis showed that fadA was important for expression of genes involved in oxidation-reduction and cellular metabolism. The flow cytometry and fluorescence microscopy demonstrated that ΔfadA accumulated more concentration of ROS in cells, and the transcriptome data indicated that genes involved in ROS scavenging were downregulated in ΔfadA mutant. We further found that FadA participated in regulating response to extracellular environmental stresses by increasing phosphorylation levels of MAPK Kinase Slt2 and Hog1. Overall, our results indicated that FadA signalling engages in mycotoxin production and A. flavus resistance to antimicrobial PAE, which provide valuable information for controlling this fungus and AF biosynthesis in pre- and postharvest of agricultural crops.

Yeasts from temperate forests

Yeasts are ubiquitous in temperate forests. While this broad habitat is well-defined, the yeasts inhabiting it and their life cycles, niches, and contributions to ecosystem functioning are less understood. Yeasts are present on nearly all sampled substrates in temperate forests worldwide. They associate with soils, macroorganisms, and other habitats, and no doubt contribute to broader ecosystem-wide processes. Researchers have gathered information leading to hypotheses about yeasts' niches and their life cycles based on physiological observations in the laboratory as well as genomic analyses, but the challenge remains to test these hypotheses in the forests themselves. Here we summarize the habitat and global patterns of yeast diversity, give some information on a handful of well-studied temperate forest yeast genera, discuss the various strategies to isolate forest yeasts, and explain temperate forest yeasts' contributions to biotechnology. We close with a summary of the many future directions and outstanding questions facing researchers in temperate forest yeast ecology. Yeasts present an exciting opportunity to better understand the hidden world of microbial ecology in this threatened and global habitat.

Characterization of 260 Isolates of Aspergillus Section Flavi Obtained from Sesame Seeds in Punjab, Pakistan

Sesame Sesamum indicum L. is a major oil-based seed crop that has been widely cultivated and consumed in Pakistan. Unfortunately, sesame is highly prone to Aspergillus fungal growth in the field, and under inappropriate storage conditions can become contaminated with aflatoxins, the most potent carcinogen found in nature. Here, we have isolated a high number of Aspergillus isolates from sesame seeds in fresh and stored conditions obtained from rainfed and irrigated zones of Punjab, Pakistan, and characterized them for aflatoxigenic potentials. Using morphological identification techniques, 260 isolates were grouped as potential Aspergillus section Flavi, with 126 and 134 originating from the rainfed and irrigated zones, respectively. Out of 260 in total, 188 isolates were confirmed to produce aflatoxins. There were no significant differences in potential aflatoxigenic isolates with respect to the rainfed and irrigated zones. However, the number of potential aflatoxigenic isolates was significantly higher (p < 0.05) in stored samples than that of those from fresh sesame seeds in the rainfed and irrigated zone. Whole genome sequencing and comparative analyses of 12 select isolates have revealed that one of the A. flavus isolates, which produced very low aflatoxins (AFP10), has an elevated missense variant rate, numerous high impact mutations, and a 600 base pair deletion in the norB gene. In summary, our study provides insights into aflatoxigenic potential and the associated genetic diversity of indigenous Aspergillus section Flavi isolates and potential management strategies for reducing aflatoxin contamination levels in a major crop consumed in Punjab, Pakistan.

Complete Mitochondrial Genome Sequences of Nine Aspergillus flavus Strains

Nuclear genome sequences incompletely characterize the genomic content and thus the genetic diversity of fungal species. Here, we present the complete mitochondrial genome sequences of nine Aspergillus flavus strains, providing useful information for inter- and intraspecific analyses.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Genetic fingerprinting and aflatoxin production of Aspergillus section Flavi associated with groundnut in eastern Ethiopia

BACKGROUND

Aspergillus species cause aflatoxin contamination in groundnut kernels, being a health threat in agricultural products and leading to commodity rejection by domestic and international markets. Presence of Aspergillus flavus and A. parasiticus colonizing groundnut in eastern Ethiopia, as well as presence of aflatoxins have been reported, though in this region, no genetic studies have been done of these species in relation to their aflatoxin production.

RESULTS

In this study, 145 Aspergillus isolates obtained from groundnut kernels in eastern Ethiopia were genetically fingerprinted using 23 Insertion/Deletion (InDel) markers within the aflatoxin-biosynthesis gene cluster (ABC), identifying 133 ABC genotypes. Eighty-four isolates were analyzed by Ultra-Performance Liquid Chromatography (UPLC) for in vitro aflatoxin production. Analysis of genetic distances based on the approximately 85 kb-ABC by Neighbor Joining (NJ), 3D-Principal Coordinate Analysis (3D-PCoA), and Structure software, clustered the isolates into three main groups as a gradient in their aflatoxin production. Group I, contained 98% A. flavus, including L- and non-producers of sclerotia (NPS), producers of B₁ and B₂ aflatoxins, and most of them collected from the lowland-dry Babile area. Group II was a genetic admixture population of A. flavus (NPS) and A. flavus S morphotype, both low producers of aflatoxins. Group III was primarily represented by A. parasiticus and A. flavus S morphotype isolates both producers of B₁, B₂ and G₁, G₂ aflatoxins, and originated from the regions of Darolabu and Gursum. The highest in vitro producer of aflatoxin B₁ was A. flavus NPS N1436 (77.98 μg/mL), and the highest producer of aflatoxin G₁ was A. parasiticus N1348 (50.33 μg/mL), these isolates were from Gursum and Darolabu, respectively.

CONCLUSIONS

To the best of our knowledge, this is the first study that combined the use of InDel fingerprinting of the ABC and corresponding aflatoxin production capability to describe the genetic diversity of Aspergillus isolates from groundnut in eastern Ethiopia. Three InDel markers, AFLC04, AFLC08 and AFLC19, accounted for the main assignment of individuals to the three Groups; their loci corresponded to aflC (pksA), hypC, and aflW (moxY) genes, respectively. Despite InDels within the ABC being often associated to loss of aflatoxin production, the vast InDel polymorphism observed in the Aspergillus isolates did not completely impaired their aflatoxin production in vitro.

Updates on the Functions and Molecular Mechanisms of the Genes Involved in Aspergillus flavus Development and Biosynthesis of Aflatoxins

Aspergillus flavus (A. flavus) is a ubiquitous and opportunistic fungal pathogen that causes invasive and non-invasive aspergillosis in humans and animals. This fungus is also capable of infecting a large number of agriculture crops (e.g., peanuts, maze, cotton seeds, rice, etc.), causing economic losses and posing serious food-safety concerns when these crops are contaminated with aflatoxins, the most potent naturally occurring carcinogens. In particular, A. flavus and aflatoxins are intensely studied, and they continue to receive considerable attention due to their detrimental effects on humans, animals, and crops. Although several studies have been published focusing on the biosynthesis of the aforementioned secondary metabolites, some of the molecular mechanisms (e.g., posttranslational modifications, transcription factors, transcriptome, proteomics, metabolomics and transcriptome, etc.) involved in the fungal development and aflatoxin biosynthesis in A. flavus are still not fully understood. In this study, a review of the recently published studies on the function of the genes and the molecular mechanisms involved in development of A. flavus and the production of its secondary metabolites is presented. It is hoped that the information provided in this review will help readers to develop effective strategies to reduce A. flavus infection and aflatoxin production.

Chromosome assembled and annotated genome sequence of Aspergillus flavus NRRL 3357

Aspergillus flavus is an opportunistic pathogen of crops, including peanuts and maize, and is the second leading cause of aspergillosis in immunocompromised patients. A. flavus is also a major producer of the mycotoxin, aflatoxin, a potent carcinogen, which results in significant crop losses annually. The A. flavus isolate NRRL 3357 was originally isolated from peanut and has been used as a model organism for understanding the regulation and production of secondary metabolites, such as aflatoxin. A draft genome of NRRL 3357 was previously constructed, enabling the development of molecular tools and for understanding population biology of this particular species. Here, we describe an updated, near complete, telomere-to-telomere assembly and re-annotation of the eight chromosomes of A. flavus NRRL 3357 genome, accomplished via long-read PacBio and Oxford Nanopore technologies combined with Illumina short-read sequencing. A total of 13,715 protein-coding genes were predicted. Using RNA-seq data, a significant improvement was achieved in predicted 5’ and 3’ untranslated regions, which were incorporated into the new gene models.

Transcription Factor Repurposing Offers Insights into Evolution of Biosynthetic Gene Cluster Regulation

The fungal kingdom has provided advances in our ability to identify biosynthetic gene clusters (BGCs) and to examine how gene composition of BGCs evolves across species and genera. However, little is known about the evolution of specific BGC regulators that mediate how BGCs produce secondary metabolites (SMs). A bioinformatics search for conservation of the Aspergillus fumigatus xanthocillin BGC revealed an evolutionary trail of xan-like BGCs across Eurotiales species. Although the critical regulatory and enzymatic genes were conserved in Penicillium expansum, overexpression (OE) of the conserved xan BGC transcription factor (TF) gene, PexanC, failed to activate the putative xan BGC transcription or xanthocillin production in P. expansum, in contrast to the role of AfXanC in A. fumigatus. Surprisingly, OE::PexanC was instead found to promote citrinin synthesis in P. expansum via trans induction of the cit pathway-specific TF, ctnA, as determined by cit BGC expression and chemical profiling of ctnA deletion and OE::PexanC single and double mutants. OE::AfxanC results in significant increases of xan gene expression and metabolite synthesis in A. fumigatus but had no effect on either xanthocillin or citrinin production in P. expansum. Bioinformatics and promoter mutation analysis led to the identification of an AfXanC binding site, 5'-AGTCAGCA-3', in promoter regions of the A. fumigatus xan BGC genes. This motif was not in the ctnA promoter, suggesting a different binding site of PeXanC. A compilation of a bioinformatics examination of XanC orthologs and the presence/absence of the 5'-AGTCAGCA-3' binding motif in xan BGCs in multiple Aspergillus and Penicillium spp. supports an evolutionary divergence of XanC regulatory targets that we speculate reflects an exaptation event in the Eurotiales. IMPORTANCE Fungal secondary metabolites (SMs) are an important source of pharmaceuticals on one hand and toxins on the other. Efforts to identify the biosynthetic gene clusters (BGCs) that synthesize SMs have yielded significant insights into how variation in the genes that compose BGCs may impact subsequent metabolite production within and between species. However, the role of regulatory genes in BGC activation is less well understood. Our finding that the bZIP transcription factor XanC, located in the xanthocillin BGC of both Aspergillus fumigatus and Penicillium expansum, has functionally diverged to regulate different BGCs in these two species emphasizes that the diversification of BGC regulatory elements may sometimes occur through exaptation, which is the co-option of a gene that evolved for one function to a novel function. Furthermore, this work suggests that the loss/gain of transcription factor binding site targets may be an important mediator in the evolution of secondary-metabolism regulatory elements.

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.

Chromosome-Level Aspergillus flavus Strain CA14 Genome Assembly

We report here a chromosome-level genome assembly of the aflatoxigenic fungus Aspergillus flavus strain CA14. This strain is the basis for numerous studies in fungal physiology and secondary metabolism. This full-length assembly will aid in subsequent genomics research.

We report here a chromosome-level genome assembly of the aflatoxigenic fungus Aspergillus flavus strain CA14. This strain is the basis for numerous studies in fungal physiology and secondary metabolism. This full-length assembly will aid in subsequent genomics research.

Ssu72 Regulates Fungal Development, Aflatoxin Biosynthesis and Pathogenicity in Aspergillus flavus

The RNA polymerase II (Pol II) transcription process is coordinated by the reversible phosphorylation of its largest subunit-carboxy terminal domain (CTD). Ssu72 is identified as a CTD phosphatase with specificity for phosphorylation of Ser5 and Ser7 and plays critical roles in regulation of transcription cycle in eukaryotes. However, the biofunction of Ssu72 is still unknown in Aspergillus flavus, which is a plant pathogenic fungus and produces one of the most toxic mycotoxins-aflatoxin. Here, we identified a putative phosphatase Ssu72 and investigated the function of Ssu72 in A. flavus. Deletion of ssu72 resulted in severe defects in vegetative growth, conidiation and sclerotia formation. Additionally, we found that phosphatase Ssu72 positively regulates aflatoxin production through regulating expression of aflatoxin biosynthesis cluster genes. Notably, seeds infection assays indicated that phosphatase Ssu72 is crucial for pathogenicity of A. flavus. Furthermore, the Δssu72 mutant exhibited more sensitivity to osmotic and oxidative stresses. Taken together, our study suggests that the putative phosphatase Ssu72 is involved in fungal development, aflatoxin production and pathogenicity in A. flavus, and may provide a novel strategy to prevent the contamination of this pathogenic fungus.

Chemical repertoire and biosynthetic machinery of the Aspergillus flavus secondary metabolome: A review

Filamentous fungi represent a rich source of extrolites, including secondary metabolites (SMs) comprising a great variety of astonishing structures and interesting bioactivities. State-of-the-art techniques in genome mining, genetic manipulation, and secondary metabolomics have enabled the scientific community to better elucidate and more deeply appreciate the genetic and biosynthetic chemical arsenal of these microorganisms. Aspergillus flavus is best known as a contaminant of food and feed commodities and a producer of the carcinogenic family of SMs, aflatoxins. This fungus produces many SMs including polyketides, ribosomal and nonribosomal peptides, terpenoids, and other hybrid molecules. This review will discuss the chemical diversity, biosynthetic pathways, and biological/ecological role of A. flavus SMs, as well as their significance concerning food safety and security.

Transcriptome Sequencing Revealed an Inhibitory Mechanism of Aspergillus flavus Asexual Development and Aflatoxin Metabolism by Soy-Fermenting Non-Aflatoxigenic Aspergillus

Aflatoxins (AFs) have always been regarded as the most effective carcinogens, posing a great threat to agriculture, food safety, and human health. Aspergillus flavus is the major producer of aflatoxin contamination in crops. The prevention and control of A. flavus and aflatoxin continues to be a global problem. In this study, we demonstrated that the cell-free culture filtrate of Aspergillus oryzae and a non-aflatoxigenic A. flavus can effectively inhibit the production of AFB1 and the growth and reproduction of A. flavus, indicating that both of the non-aflatoxigenic Aspergillus strains secrete inhibitory compounds. Further transcriptome sequencing was performed to analyze the inhibitory mechanism of A. flavus treated with fermenting cultures, and the results revealed that genes involved in the AF biosynthesis pathway and other biosynthetic gene clusters were significantly downregulated, which might be caused by the reduced expression of specific regulators, such as AflS, FarB, and MtfA. The WGCNA results further revealed that genes involved in the TCA cycle and glycolysis were potentially involved in aflatoxin biosynthesis. Our comparative transcriptomics also revealed that two conidia transcriptional factors, brlA and abaA, were found to be significantly downregulated, which might lead to the downregulation of conidiation-specific genes, such as the conidial hydrophobins genes rodA and rodB. In summary, our research provides new insights for the molecular mechanism of controlling AF synthesis to control the proliferation of A. flavus and AF pollution.

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.