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

Genomic insights into the evolution of secondary metabolism of Escovopsis and its allies, specialized fungal symbionts of fungus-farming ants

The metabolic intimacy of symbiosis often demands the work of specialists. Natural products and defensive secondary metabolites can drive specificity by ensuring infection and propagation across host generations. But in contrast to bacteria, little is known about the diversity and distribution of natural product biosynthetic pathways among fungi and how they evolve to facilitate symbiosis and adaptation to their host environment. In this study, we define the secondary metabolism of Escovopsis and closely related genera, symbionts in the gardens of fungus-farming ants. We ask how the gain and loss of various biosynthetic pathways correspond to divergent lifestyles. Long-read sequencing allowed us to define the chromosomal features of representative Escovopsis strains, revealing highly reduced genomes composed of seven to eight chromosomes. The genomes are highly syntenic with macrosynteny decreasing with increasing phylogenetic distance, while maintaining a high degree of mesosynteny. An ancestral state reconstruction analysis of biosynthetic pathways revealed that, while many secondary metabolites are shared with non-ant-associated Sordariomycetes, 56 pathways are unique to the symbiotic genera. Reflecting adaptation to diverging ant agricultural systems, we observe that the stepwise acquisition of these pathways mirrors the ecological radiations of attine ants and the dynamic recruitment and replacement of their fungal cultivars. As different clades encode characteristic combinations of biosynthetic gene clusters, these delineating profiles provide important insights into the possible mechanisms underlying specificity between these symbionts and their fungal hosts. Collectively, our findings shed light on the evolutionary dynamic nature of secondary metabolism in Escovopsis and its allies, reflecting adaptation of the symbionts to an ancient agricultural system.IMPORTANCEMicrobial symbionts interact with their hosts and competitors through a remarkable array of secondary metabolites and natural products. Here, we highlight the highly streamlined genomic features of attine-associated fungal symbionts. The genomes of Escovopsis species, as well as species from other symbiont genera, many of which are common with the gardens of fungus-growing ants, are defined by seven chromosomes. Despite a high degree of metabolic conservation, we observe some variation in the symbionts' potential to produce secondary metabolites. As the phylogenetic distribution of the encoding biosynthetic gene clusters coincides with attine transitions in agricultural systems, we highlight the likely role of these metabolites in mediating adaptation by a group of highly specialized symbionts.

Pathogenicity is associated with population structure in a fungal pathogen of humans

Aspergillus flavus is a clinically and agriculturally important saprotrophic fungus responsible for severe human infections and extensive crop losses. We analyzed genomic data from 250 (95 clinical and 155 environmental) A. flavus isolates from 9 countries, including 70 newly sequenced clinical isolates, to examine population and pan-genome structure and their relationship to pathogenicity. We identified five A. flavus populations, including a new population, D, corresponding to distinct clades in the genome-wide phylogeny. Strikingly, > 75% of clinical isolates were from population D. Accessory genes, including genes within biosynthetic gene clusters, were significantly more common in some populations but rare in others. Population D was enriched for genes associated with zinc ion binding, lipid metabolism, and certain types of hydrolase activity. In contrast to the major human pathogen Aspergillus fumigatus, A. flavus pathogenicity in humans is strongly associated with population structure, making it a great system for investigating how population-specific genes contribute to pathogenicity.

Genomic and metabolomic diversity within a familial population of Aspergillus flavus

Aspergillus flavus is an agriculturally significant micro-fungus having potential to contaminate food and feed crops with toxic secondary metabolites such as aflatoxin (AF) and cyclopiazonic acid (CPA). Research has shown A. flavus strains can overcome heterokaryon incompatibility and undergo meiotic recombination as teleomorphs. Although evidence of recombination in the AF gene cluster has been reported, the impacts of recombination on genotype and metabolomic phenotype in a single generation are lacking. In previous studies, we paired an aflatoxigenic MAT1-1 A. flavus strain with a non-aflatoxigenic MAT1-2 A. flavus strain that had been tagged with green fluorescent protein and then 10 F1 progenies (a mix of fluorescent and non-fluorescent) were randomly selected from single-ascospore colonies and broadly examined for evidence of recombination. In this study, we determined four of those 10 F1 progenies were recombinants because they were not vegetatively compatible with either parent or their siblings, and they exhibited other distinctive traits that could only result from meiotic recombination. The other six progenies examined shared genomic identity with the non-aflatoxigenic, fluorescent, and MAT1-2 parent, but were metabolically distinct. This study highlights phenotypic and genomic changes that may occur in a single generation from the outcrossing of sexually compatible strains of A. flavus.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters

BACKGROUND

Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes.

RESULTS

Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites.

CONCLUSION

In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

A Timeline of Biosynthetic Gene Cluster Discovery in Aspergillus fumigatus: From Characterization to Future Perspectives

In 1999, the first biosynthetic gene cluster (BGC), synthesizing the virulence factor DHN melanin, was characterized in Aspergillus fumigatus. Since then, 19 additional BGCs have been linked to specific secondary metabolites (SMs) in this species. Here, we provide a comprehensive timeline of A. fumigatus BGC discovery and find that initial advances centered around the commonly expressed SMs where chemical structure informed rationale identification of the producing BGC (e.g., gliotoxin, fumigaclavine, fumitremorgin, pseurotin A, helvolic acid, fumiquinazoline). Further advances followed the transcriptional profiling of a ΔlaeA mutant, which aided in the identification of endocrocin, fumagillin, hexadehydroastechrome, trypacidin, and fumisoquin BGCs. These SMs and their precursors are the commonly produced metabolites in most A. fumigatus studies. Characterization of other BGC/SM pairs required additional efforts, such as induction treatments, including co-culture with bacteria (fumicycline/neosartoricin, fumigermin) or growth under copper starvation (fumivaline, fumicicolin). Finally, four BGC/SM pairs were discovered via overexpression technologies, including the use of heterologous hosts (fumicycline/neosartoricin, fumihopaside, sphingofungin, and sartorypyrone). Initial analysis of the two most studied A. fumigatus isolates, Af293 and A1160, suggested that both harbored ca. 34-36 BGCs. However, an examination of 264 available genomes of A. fumigatus shows up to 20 additional BGCs, with some strains showing considerable variations in BGC number and composition. These new BGCs present a new frontier in the future of secondary metabolism characterization in this important species.

Insight into the role of niche concept in deciphering the ecological drivers of MPs-associated bacterial communities in mangrove forest

Myriad inherent and variable environmental features are controlling the assembly and succession of bacterial communities colonizing on mangrove microplastics (MPs). However, the mechanisms governing mangrove MPs-associated bacterial responses to environmental changes still remain unknown. Here, we assessed the dissimilarities of MPs-associated bacterial composition, diversity and functionality as well as quantified the niche variations of each taxon on plastispheres along river-mangrove-ocean and mangrove landward-to-seaward gradients in the Beibu Gulf, China, respectively. The bacterial richness and diversity as well as the niche breadth on mangrove sedimentary MPs dramatically decreased from landward to seaward regions. Characterizing the niche variations linked the difference of ecological drivers of MPs-associated bacterial populations and functions between river-mangrove-ocean (microplastic properties) and mangrove landward-to-seaward plastispheres (sediment physicochemical properties) to the trade-offs between selective stress exerted by inherent plastic substrates and microbial competitive stress imposed by environmental conditions. Notably, Rhodococcus erythropolis was predicted to be the generalist species and closely associated to biogeochemical cycles of mangrove plastispheres. Our work provides a reliable pathway for tackling the hidden mechanisms of environmental factors driving MPs-associated microbe from perspectives of niches and highlights the spatial dynamic variations of mangrove MPs-associated bacteria.

The Development of Aspergillus flavus and Biosynthesis of Aflatoxin B1 are Regulated by the Golgi-Localized Mn2+ Transporter Pmr1

Microorganisms rely on diverse ion transport and trace elements to sustain growth, development, and secondary metabolism. Manganese (Mn2+) is essential for various biological processes and plays a crucial role in the metabolism of human cells, plants, and yeast. In Aspergillus flavus, we confirmed that Pmr1 localized in cis- and medial-Golgi compartments was critical in facilitating Mn2+ transport, fungal growth, development, secondary metabolism, and glycosylation. In comparison to the wild type, the Δpmr1 mutant displayed heightened sensitivity to environmental stress, accompanied by inhibited synthesis of aflatoxin B1, kojic acid, and a substantial reduction in pathogenicity toward peanuts and maize. Interestingly, the addition of exogenous Mn2+ effectively rectified the developmental and secondary metabolic defects in the Δpmr1 mutant. However, Mn2+ supplement failed to restore the growth and development of the Δpmr1Δgdt1 double mutant, which indicated that the Gdt1 compensated for the functional deficiency of pmr1. In addition, our results showed that pmr1 knockout leads to an upregulation of O-glycosyl-N-acetylglucose (O-GlcNAc) and O-GlcNAc transferase (OGT), while Mn2+ supplementation can restore the glycosylation in A. flavus. Collectively, this study indicates that the pmr1 regulates Mn2+ via Golgi and maintains growth and metabolism functions of A. flavus through regulation of the glycosylation.

The KdmB-EcoA-RpdA-SntB (KERS) chromatin regulatory complex controls development, secondary metabolism and pathogenicity in Aspergillus flavus

The filamentous fungus Aspergillus flavus is a plant and human pathogen predominantly found in the soil as spores or sclerotia and is capable of producing various secondary metabolites (SM) such as the carcinogenic mycotoxin aflatoxin. Recently, we have discovered a novel nuclear chromatin binding complex (KERS) that contains the JARID1-type histone demethylase KdmB, a putative cohesion acetyl transferase EcoA, a class I type histone deacetylase RpdA and the PHD ring finger reader protein SntB in the model filamentous fungus Aspergillus nidulans. Here, we show the presence of the KERS complex in A. flavus by immunoprecipitation-coupled mass spectrometry and constructed kdmBΔ and rpdAΔ strains to study their roles in fungal development, SM production and histone post-translational modifications (HPTMs). We found that KdmB and RpdA couple the regulation of SM gene clusters with fungal light-responses and HPTMs. KdmB and RpdA have opposing roles in light-induced asexual conidiation, while both factors are positive regulators of sclerotia development through the nsdC and nsdD pathway. KdmB and RpdA are essential for the productions of aflatoxin (similar to findings for SntB) as well as cyclopiazonic acid, ditryptophenaline and leporin B through controlling the respective SM biosynthetic gene clusters. We further show that both KdmB and RpdA regulate H3K4me3 and H3K9me3 levels, while RpdA also acts on H3K14ac levels in nuclear extracts. Therefore, the chromatin modifiers KdmB and RpdA of the KERS complex are key regulators for fungal development and SM metabolism in A. flavus.

Fungal secondary metabolism is governed by an RNA-binding protein CsdA/RsdA complex

Production of secondary metabolites is controlled by a complicated regulatory network in eukaryotic cells. Several layers of regulators are involved in this process, ranging from pathway-specific regulation, to epigenetic control, to global regulation. Here, we discover that interaction of an RNA-binding protein CsdA with a regulator RsdA coordinates fungal secondary metabolism. Employing a genetic deletion approach and transcriptome analysis as well as metabolomics analysis, we reveal that CsdA and RsdA synergistically regulate fungal secondary metabolism comprehensively. Mechanistically, comprehensive genetic and biochemical studies prove that RsdA and CsdA co-localize in the nucleus and physically interact to achieve their functions. In particular, we demonstrate that CsdA mediates rsdA expression by binding specific motif "GUCGGUAU" of its pre-mRNA at a post-transcriptional level. We thus uncover a mechanism in which RNA-binding protein physically interacts with, and controls the expression level of, the RsdA to coordinate fungal secondary metabolism.

Mining for a new class of fungal natural products: the evolution, diversity, and distribution of isocyanide synthase biosynthetic gene clusters

The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We amalgamated a pipeline of tools to predict BGCs based on shared promoter motifs and located 3800 ICS BGCs in 3300 genomes, making ICS BGCs the fifth largest class of specialized metabolites compared to canonical classes found by antiSMASH. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1/2 gene cluster family (GCF), which was prior only studied in yeast, is present in ∼30% of all Ascomycetes. The dit variety ICS exhibits greater similarity to bacterial ICS than other fungal ICS, suggesting a potential convergence of the ICS backbone domain. The evolutionary origins of the dit GCF in Ascomycota are ancient and these genes are diversifying in some lineages. Our results create a roadmap for future research into ICS BGCs. We developed a website (https://isocyanides.fungi.wisc.edu/) that facilitates the exploration and downloading of all identified fungal ICS BGCs and GCFs.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Correlative metabologenomics of 110 fungi reveals metabolite-gene cluster pairs

Natural products research increasingly applies -omics technologies to guide molecular discovery. While the combined analysis of genomic and metabolomic datasets has proved valuable for identifying natural products and their biosynthetic gene clusters (BGCs) in bacteria, this integrated approach lacks application to fungi. Because fungi are hyper-diverse and underexplored for new chemistry and bioactivities, we created a linked genomics-metabolomics dataset for 110 Ascomycetes, and optimized both gene cluster family (GCF) networking parameters and correlation-based scoring for pairing fungal natural products with their BGCs. Using a network of 3,007 GCFs (organized from 7,020 BGCs), we examined 25 known natural products originating from 16 known BGCs and observed statistically significant associations between 21 of these compounds and their validated BGCs. Furthermore, the scalable platform identified the BGC for the pestalamides, demystifying its biogenesis, and revealed more than 200 high-scoring natural product-GCF linkages to direct future discovery.

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.

Mining for a New Class of Fungal Natural Products: The Evolution, Diversity, and Distribution of Isocyanide Synthase Biosynthetic Gene Clusters

The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) have notable bioactivities that mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We developed the first genome-mining pipeline to identify ICS BGCs, locating 3,800 ICS BGCs in 3,300 genomes. Genes in these clusters share promoter motifs and are maintained in contiguous groupings by natural selection. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1 / 2 gene cluster family (GCF), which was thought to only exist in yeast, is present in ∼30% of all Ascomycetes, including many filamentous fungi. The evolutionary history of the dit GCF is marked by deep divergences and phylogenetic incompatibilities that raise questions about convergent evolution and suggest selection or horizontal gene transfers have shaped the evolution of this cluster in some yeast and dimorphic fungi. Our results create a roadmap for future research into ICS BGCs. We developed a website ( www.isocyanides.fungi.wisc.edu ) that facilitates the exploration, filtering, and downloading of all identified fungal ICS BGCs and GCFs.

Evolutionary investigations of the biosynthetic diversity in the skin microbiome using lsaBGC

Bacterial secondary metabolites, synthesized by enzymes encoded in biosynthetic gene clusters (BGCs), can underlie microbiome homeostasis and serve as commercialized products, which have historically been mined from a select group of taxa. While evolutionary approaches have proven beneficial for prioritizing BGCs for experimental characterization efforts to uncover new natural products, dedicated bioinformatics tools designed for comparative and evolutionary analysis of BGCs within focal taxa are limited. We thus developed lineage specific analysis of BGCs (lsaBGC; https://github.com/Kalan-Lab/lsaBGC) to aid exploration of microdiversity and evolutionary trends across homologous groupings of BGCs, gene cluster families (GCFs), in any bacterial taxa of interest. lsaBGC enables rapid and direct identification of GCFs in genomes, calculates evolutionary statistics and conservation for BGC genes, and builds a framework to allow for base resolution mining of novel variants through metagenomic exploration. Through application of the suite to four genera commonly found in skin microbiomes, we uncover new insights into the evolution and diversity of their BGCs. We show that the BGC of the virulence-associated carotenoid staphyloxanthin in Staphylococcus aureus is ubiquitous across the genus Staphylococcus . While one GCF encoding the biosynthesis of staphyloxanthin showcases evidence for plasmid-mediated horizontal gene transfer (HGT) between species, another GCF appears to be transmitted vertically amongst a sub-clade of skin-associated Staphylococcus . Further, the latter GCF, which is well conserved in S. aureus , has been lost in most Staphylococcus epidermidis , which is the most common Staphylococcus species on human skin and is also regarded as a commensal. We also identify thousands of novel single-nucleotide variants (SNVs) within BGCs from the Corynebacterium tuberculostearicum sp. complex, a narrow, multi-species clade that features the most prevalent Corynebacterium in healthy skin microbiomes. Although novel SNVs were approximately 10 times as likely to correspond to synonymous changes when located in the top five percentile of conserved sites, lsaBGC identified SNVs that defied this trend and are predicted to underlie amino acid changes within functionally key enzymatic domains. Ultimately, beyond supporting evolutionary investigations of BGCs, lsaBGC also provides important functionalities to aid efforts for the discovery or directed modification of natural products.

Toxicity sharing model of earthworm intestinal microbiome reveals shared functional genes are more powerful than species in resisting pesticide stress

Earthworm intestinal bacteria and indigenous soil bacteria work closely during various biochemical processes and play a crucial role in maintaining the internal stability of the soil environment. However, the response mechanism of these bacterial communities to external pesticide disturbance is unknown. In this study, soil and earthworm gut contents were metagenomically sequenced after exposure to various concentrations of nitrochlorobenzene (0-1026.7 mg kg^-1). A high degree of similarity was found between the microbial community composition and abundance in the worm gut and soil, both of which decreased significantly (P < 0.05) under elevated pesticide stress. The toxicity sharing model (TSM) showed that the toxicity sharing capacity was 97.4-125.7 % and 100.4-130.2 % for E_genes (genes in the worm gut) and E_met(degradation genes in the worm gut) in the earthworm intestinal microbiome, respectively. This indicated that the earthworm intestinal microbiome assisted in relieving the pesticide toxicity of the indigenous soil microbiome. This study showed that the TSM could quantitatively describe the toxic effect of pesticides on the earthworm intestinal microbiome. It provides a new analytical model for investigating the ecological alliance between earthworm intestinal microbiome and indigenous soil microbiome under pesticide stress while contributing a more profound understanding of the potential to use earthworms to mitigate pesticide pollution in soils and develop earthworm-based soil remediation techniques.

Lysine Succinylation of VBS Contributes to Sclerotia Development and Aflatoxin Biosynthesis in Aspergillus flavus

Aspergillus flavus is a common saprophytic and pathogenic fungus, and its secondary metabolic pathways are one of the most highly characterized owing to its aflatoxin (AF) metabolite affecting global economic crops and human health. Different natural environments can cause significant variations in AF synthesis. Succinylation was recently identified as one of the most critical regulatory post-translational modifications affecting metabolic pathways. It is primarily reported in human cells and bacteria with few studies on fungi. Proteomic quantification of lysine succinylation (Ksuc) exploring its potential involvement in secondary metabolism regulation (including AF production) has not been performed under natural conditions in A. flavus. In this study, a quantification method was performed based on tandem mass tag labeling and antibody-based affinity enrichment of succinylated peptides via high accuracy nano-liquid chromatography with tandem mass spectrometry to explore the succinylation mechanism affecting the pathogenicity of naturally isolated A. flavus strains with varying toxin production. Altogether, 1240 Ksuc sites in 768 proteins were identified with 1103 sites in 685 proteins quantified. Comparing succinylated protein levels between high and low AF-producing A. flavus strains, bioinformatics analysis indicated that most succinylated proteins located in the AF biosynthetic pathway were downregulated, which directly affected AF synthesis. Versicolorin B synthase is a key catalytic enzyme for heterochrome B synthesis during AF synthesis. Site-directed mutagenesis and biochemical studies revealed that versicolorin B synthase succinylation is an important regulatory mechanism affecting sclerotia development and AF biosynthesis in A. flavus. In summary, our quantitative study of the lysine succinylome in high/low AF-producing strains revealed the role of Ksuc in regulating AF biosynthesis. We revealed novel insights into the metabolism of AF biosynthesis using naturally isolated A. flavus strains and identified a rich source of metabolism-related enzymes regulated by succinylation.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi

Introduction

Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria.

Methods

The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation.

Results

In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations.

Discussion

We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Genomic and Phenotypic Trait Variation of the Opportunistic Human Pathogen Aspergillus flavus and Its Close Relatives

Fungal diseases affect millions of humans annually, yet fungal pathogens remain understudied. The mold Aspergillus flavus can cause both aspergillosis and fungal keratitis infections, but closely related species are not considered clinically relevant. To study the evolution of A. flavus pathogenicity, we examined genomic and phenotypic traits of two strains of A. flavus and three closely related species, Aspergillus arachidicola (two strains), Aspergillus parasiticus (two strains), and Aspergillus nomiae (one strain). We identified >3,000 orthologous proteins unique to A. flavus, including seven biosynthetic gene clusters present in A. flavus strains and absent in the three nonpathogens. We characterized secondary metabolite production for all seven strains under two clinically relevant conditions, temperature and salt concentration. Temperature impacted metabolite production in all species, whereas salinity did not affect production of any species. Strains of the same species produced different metabolites. Growth under stress conditions revealed additional heterogeneity within species. Using the invertebrate fungal disease model Galleria mellonella, we found virulence of strains of the same species varied widely; A. flavus strains were not more virulent than strains of the nonpathogens. In a murine model of fungal keratitis, we observed significantly lower disease severity and corneal thickness for A. arachidicola compared to other species at 48 h postinfection, but not at 72 h. Our work identifies variations in key phenotypic, chemical, and genomic attributes between A. flavus and its nonpathogenic relatives and reveals extensive strain heterogeneity in virulence that does not correspond to the currently established clinical relevance of these species. IMPORTANCE Aspergillus flavus is a filamentous fungus that causes opportunistic human infections, such as aspergillosis and fungal keratitis, but its close relatives are considered nonpathogenic. To begin understanding how this difference in pathogenicity evolved, we characterized variation in infection-relevant genomic, chemical, and phenotypic traits between strains of A. flavus and its relatives. We found extensive variation (or strain heterogeneity) within the pathogenic A. flavus as well as within its close relatives, suggesting that strain-level differences may play a major role in the ability of these fungi to cause disease. Surprisingly, we also found that the virulence of strains from species not considered to be pathogens was similar to that of A. flavus in both invertebrate and murine models of disease. These results contrast with previous studies on Aspergillus fumigatus, another major pathogen in the genus, for which significant differences in infection-relevant chemical and phenotypic traits are observed between closely related pathogenic and nonpathogenic species.

Why Do Plant-Pathogenic Fungi Produce Mycotoxins? Potential Roles for Mycotoxins in the Plant Ecosystem

For many plant-pathogenic or endophytic fungi, production of mycotoxins, which are toxic to humans, may present a fitness gain. However, associations between mycotoxin production and plant pathogenicity or virulence is inconsistent and difficult due to the complexity of these host-pathogen interactions and the influences of environmental and insect factors. Aflatoxin receives a lot of attention due to its potent toxicity and carcinogenicity but the connection between aflatoxin production and pathogenicity is complicated by the pathogenic ability and prevalence of nonaflatoxigenic isolates in crops. Other toxins directly aid fungi in planta, trichothecenes are important virulence factors, and ergot alkaloids limit herbivory and fungal consumption due to insect toxicity. We review a panel discussion at the American Phytopathological Society's Plant Health 2021 conference, which gathered diverse experts representing different research sectors, career stages, ethnicities, and genders to discuss the diverse roles of mycotoxins in the lifestyles of filamentous fungi of the families Clavicipitaceae, Trichocomaceae (Eurotiales), and Nectriaceae (Hypocreales).

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.

Ecological generalism drives hyperdiversity of secondary metabolite gene clusters in xylarialean endophytes

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Comprehensive Guide to Extracting and Expressing Fungal Secondary Metabolites with Aspergillus fumigatus as a Case Study

Fungal secondary metabolites (SMs) have captured the interest of natural products researchers in academia and industry for decades. In recent years, the high rediscovery rate of previously characterized metabolites is making it increasingly difficult to uncover novel compounds. Additionally, the vast majority of fungal SMs reside in genetically intractable fungi or are silent under normal laboratory conditions in genetically tractable fungi. The fungal natural products community has broadly overcome these barriers by altering the physical growth conditions of the fungus and heterologous/homologous expression of biosynthetic gene cluster regulators or proteins. The protocols described here summarize vital methodologies needed when researching SM production in fungi. We also summarize the growth conditions, genetic backgrounds, and extraction protocols for every published SM in Aspergillus fumigatus, enabling readers to easily replicate the production of previously characterized SMs. Readers will also be equipped with the tools for developing their own strategy for expressing and extracting SMs from their given fungus or a suitable heterologous model system. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Making glycerol stocks from spore suspensions Alternate Protocol 1: Creating glycerol stocks from non-sporulating filamentous fungi Basic Protocol 2: Activating spore-suspension glycerol stocks Basic Protocol 3: Extracting secondary metabolites from Aspergillus spp grown on solid medium Alternate Protocol 2: Extracting secondary metabolites from Aspergillus spp using ethyl acetate Alternate Protocol 3: High-volume metabolite extraction using ethyl acetate Alternate Protocol 4: Extracting secondary metabolites from Aspergillus spp in liquid medium Support Protocol: Creating an overlay culture Basic Protocol 4: Extracting DNA from filamentous fungi Basic Protocol 5: Creating a DNA construct with double-joint PCR Alternate Protocol 5: Creating a DNA construct with yeast recombineering Basic Protocol 6: Transformation of Aspergillus spp Basic Protocol 7: Co-culturing fungi and bacteria for extraction of secondary metabolites.

The confluence of big data and evolutionary genome mining for the discovery of natural products

This review covers literature between 2003-2021The development and application of genome mining tools has given rise to ever-growing genetic and chemical databases and propelled natural products research into the modern age of Big Data. Likewise, an explosion of evolutionary studies has unveiled genetic patterns of natural products biosynthesis and function that support Darwin's theory of natural selection and other theories of adaptation and diversification. In this review, we aim to highlight how Big Data and evolutionary thinking converge in the study of natural products, and how this has led to an emerging sub-discipline of evolutionary genome mining of natural products. First, we outline general principles to best utilize Big Data in natural products research, addressing key considerations needed to provide evolutionary context. We then highlight successful examples where Big Data and evolutionary analyses have been combined to provide bioinformatic resources and tools for the discovery of novel natural products and their biosynthetic enzymes. Rather than an exhaustive list of evolution-driven discoveries, we highlight examples where Big Data and evolutionary thinking have been embraced for the evolutionary genome mining of natural products. After reviewing the nascent history of this sub-discipline, we discuss the challenges and opportunities of genomic and metabolomic tools with evolutionary foundations and/or implications and provide a future outlook for this emerging and exciting field of natural product research.

Developing fungal heterologous expression platforms to explore and improve the production of natural products from fungal biodiversity

Natural products from fungi represent an important source of biologically active metabolites notably for therapeutic agent development. Genome sequencing revealed that the number of biosynthetic gene clusters (BGCs) in fungi is much larger than expected. Unfortunately, most of them are silent or barely expressed under laboratory culture conditions. Moreover, many fungi in nature are uncultivable or cannot be genetically manipulated, restricting the extraction and identification of bioactive metabolites from these species. Rapid exploration of the tremendous number of cryptic fungal BGCs necessitates the development of heterologous expression platforms, which will facilitate the efficient production of natural products in fungal cell factories. Host selection, BGC assembly methods, promoters used for heterologous gene expression, metabolic engineering strategies and compartmentalization of biosynthetic pathways are key aspects for consideration to develop such a microbial platform. In the present review, we summarize current progress on the above challenges to promote research effort in the relevant fields.

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.

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.