Genome sequence and comparative analyses of atoxigenic Aspergillus flavus WRRL 1519

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.

Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

Current Approaches for Advancement in Understanding the Molecular Mechanisms of Mycotoxin Biosynthesis

Filamentous fungi are able to synthesise a remarkable range of secondary metabolites, which play various key roles in the interaction between fungi and the rest of the biosphere, determining their ecological fitness. Many of them can have a beneficial activity to be exploited, as well as negative impact on human and animal health, as in the case of mycotoxins contaminating large quantities of food, feed, and agricultural products worldwide and posing serious health and economic risks. The elucidation of the molecular aspects of mycotoxin biosynthesis has been greatly sped up over the past decade due to the advent of next-generation sequencing technologies, which greatly reduced the cost of genome sequencing and related omic analyses. Here, we briefly highlight the recent progress in the use and integration of omic approaches for the study of mycotoxins biosynthesis. Particular attention has been paid to genomics and transcriptomic approaches for the identification and characterisation of biosynthetic gene clusters of mycotoxins and the understanding of the regulatory pathways activated in response to physiological and environmental factors leading to their production. The latest innovations in genome-editing technology have also provided a more powerful tool for the complete explanation of regulatory and biosynthesis pathways. Finally, we address the crucial issue of the interpretation of the combined omics data on the biology of the mycotoxigenic fungi. They are rapidly expanding and require the development of resources for more efficient integration, as well as the completeness and the availability of intertwined data for the research community.

Genome sequence of Aspergillus flavus A7, a marine-derived fungus with antibacterial activity

Due to the specific properties of the marine environment, marine microorganisms have exclusive physicochemical characteristics that are different from those of terrestrial microorganisms, which can produce various secondary metabolites (SMs) with considerable structural diversity and biological activity. In this study, three strains of coepiphytic Aspergillus with potential antibacterial activities, A7 (Aspergillus flavus), B27 (Aspergillus flavipes), and R12 (Aspergillus sydowii), were isolated from the South China Sea. Via the Illumina MiSeq sequencing platform, the genomes of the three strains were sequenced, and genome comparison showed the highest diversity of the biosynthetic gene clusters (BGCs) in A7. In addition, a comparison of physiological and genomic characteristics between A7 and other A. flavus strains demonstrated the superior environmental adaptability of A7, which is apparently consistent with the genetic richness of BGCs. By assigning reads to known BGCs, putative BGCs were allocated in A7 that corresponded to various SMs, including naphthopyrone, pyranonigrin E, and cyclopiazonic acids. Based on gene homology analysis, we surmise that a region is involved in the biosynthesis of ustiloxin-like RiPPs, a less thoroughly studied SM in fungi. Our results provide genetic information for the investigation of marine Aspergillus spp., which may help to elucidate their chemical diversity and adaptive strategies.

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Evolution of Aspergillus oryzae before and after domestication inferred by large-scale comparative genomic analysis

Aspergillus oryzae is an industrially useful species, of which various strains have been identified; however, their genetic relationships remain unclear. A. oryzae was previously thought to be asexual and unable to undergo crossbreeding. However, recent studies revealed the sexual reproduction of Aspergillus flavus, a species closely related to A. oryzae. To investigate potential sexual reproduction in A. oryzae and evolutionary history among A. oryzae and A. flavus strains, we assembled 82 draft genomes of A. oryzae strains used practically. The phylogenetic tree of concatenated genes confirmed that A. oryzae was monophyletic and nested in one of the clades of A. flavus but formed several clades with different genomic structures. Our results suggest that A. oryzae strains have undergone multiple inter-genomic recombination events between A. oryzae ancestors, although sexual recombination among domesticated species did not appear to have occurred during the domestication process, at least in the past few decades. Through inter- and intra-cladal comparative analysis, we found that evolutionary pressure induced by the domestication of A. oryzae appears to selectively cause non-synonymous and gap mutations in genes involved in fermentation characteristics, as well as intra-genomic rearrangements, with the conservation of industrially useful catalytic enzyme-encoding genes.

Aspergillus flavus NRRL 35739, a Poor Biocontrol Agent, May Have Increased Relative Expression of Stress Response Genes

Biocontrol of the mycotoxin aflatoxin utilizes non-aflatoxigenic strains of Aspergillus flavus, which have variable success rates as biocontrol agents. One non-aflatoxigenic strain, NRRL 35739, is a notably poor biocontrol agent. Its growth in artificial cultures and on peanut kernels was found to be slower than that of two aflatoxigenic strains, and NRRL 35739 exhibited less sporulation when grown on peanuts. The non-aflatoxigenic strain did not greatly prevent aflatoxin accumulation. Comparison of the transcriptomes of aflatoxigenic and non-aflatoxigenic A. flavus strains AF36, AF70, NRRL 3357, NRRL 35739, and WRRL 1519 indicated that strain NRRL 35739 had increased relative expression of six heat shock and stress response proteins, with the genes having relative read counts in NRRL 35739 that were 25 to 410 times more than in the other four strains. These preliminary findings tracked with current thought that aflatoxin biocontrol efficacy is related to the ability of a non-aflatoxigenic strain to out-compete aflatoxigenic ones. The slower growth of NRRL 35739 might be due to lower stress tolerance or overexpression of stress response(s). Further study of NRRL 35739 is needed to refine our understanding of the genetic basis of competitiveness among A. flavus strains.