The evolutionary imprint of domestication on genome variation and function of the filamentous fungus Aspergillus oryzae

Genome diversity, population structure and MALDI-TOF MS profiling of Aspergillus oryzae/flavus strains from fermentation and wild environments

Various strains of Aspergillus oryzae, regarded as a domesticated variant of aflatoxigenic Aspergillus flavus, are utilized in the soybean fermentation industry of Korea. This study compared A. oryzae/flavus strains isolated from various environments in Korea including industrial settings, Meju (brick of dried fermented soybeans), and wild conditions with globally reported strains using genomic analysis to determine their taxonomic positions and risk of mycotoxicity. Using population genomics, five distinct groups (A to E) were identified, with all aflatoxigenic Korean strains in Group C and non-aflatoxigenic Korean strains in Groups A, B, and E. Korean strains from Meju and wild conditions are distributed across Groups A and B, and most of the Korean industrial strains form a sub-cluster with Japanese industrial strains in Group A. Comparing secondary metabolite gene cluster mutation pattern, three gene clusters (Aflatoxin, Cyclopiazonic acid and Ditryptophenaline) were revealed as group specific ones. In aflatoxin and cyclopiazonic acid clusters, most of the Group C strains had intact regions compared to strains in other groups. Since most of the Group C strains produce aflatoxin and have intact Aflatoxin and Cyclopiazonic acid gene clusters, we considered that this group represent A. flavus. Profiling using MALDI-TOF MS analysis also distinguished Group C from Groups A, B and E by specific three proteomic peaks. Among the three peaks, those around 12,700 to 12,900 m/z (Da) are expected to correspond to AflF (nor B), an enzyme involved in Aflatoxin metabolism. These results showed taxonomic positions of Korean strains of A. oryzae/flavus from various environments and also showed possibility to differentiate between A. oryzae and A. flavus with genome and MALDI-TOF MS analysis.

Adaptive evolution of sesquiterpene deoxyphomenone in mycoparasitism by Hansfordia pulvinata associated with horizontal gene transfer from Aspergillus species

Leaf mold caused by the ascomycete fungus Cladosporium fulvum is a devastating disease of tomato plants. The mycoparasitic fungus Hansfordia pulvinata is an effective biocontrol agent that parasitizes C. fulvum hyphae on leaves and secretes 13-deoxyphomenone, an eremophilane-type sesquiterpene, which was also identified as a sporulation-inducing factor in Aspergillus oryzae. Here, we identified deoxyphomenone biosynthesis (DPH) gene clusters conserved in both H. pulvinata and Aspergillus section Flavi, including A. oryzae and A. flavus. Functional disruption of DPH1 orthologous genes encoding sesquiterpene cyclase in H. pulvinata, A. oryzae, and its close relative A. flavus revealed that deoxyphomenone in H. pulvinata had exogenic antifungal activity against C. fulvum and controlled endogenic sporulation in Aspergillus species. Complete DPH clusters, highly similar to those in H. pulvinata, were exclusive to Aspergillus section Flavi, while species in other Aspergillus sections contained fragmented DPH clusters. A comparative genomics analysis revealed that these DPH gene clusters share a common origin and are horizontally transferred from an ancestor of Aspergillus to H. pulvinata. Our results suggest that after horizontal transfer, H. pulvinata maintained the DPH cluster as the inhibitory effect of deoxyphomenone on spore germination and mycelial growth contributed to its mycoparasitism on the host fungus C. fulvum.

IMPORTANCE

Tomato leaf mold disease caused by C. fulvum poses a significant economic threat to tomato production globally. Breeders have developed tomato cultivars with Cf resistance genes. C. fulvum frequently evolves new races that overcome these genetic defenses, complicating control efforts. Additionally, the pathogen has developed resistance to chemical fungicides, prompting the need for sustainable alternatives like biocontrol agents. The mycoparasitic fungus H. pulvinata is crucial as an effective agent against C. fulvum. Clarifying the mechanism of mycoparasitism is significant, as it enhances its application as a biocontrol agent against plant pathogens. This study revealed how H. pulvinata produces deoxyphomenone, an antifungal compound, through horizontal gene transfer from Aspergillus species. It is hypothesized that mycoparasitism could be one of the mechanisms that facilitated horizontal gene transfer between fungi. These insights facilitate the development of eco-friendly, sustainable agricultural practices by reducing dependence on chemical fungicides and promoting natural pathogen control methods.

Population Genomics of Aspergillus sojae is Shaped by the Food Environment

Traditional fermented foods often contain specialized microorganisms adapted to their unique environments. For example, the filamentous mold Aspergillus oryzae, used in saké fermentation, has evolved to thrive in starch-rich conditions compared to its wild ancestor, Aspergillus flavus. Similarly, Aspergillus sojae, used in soybean-based fermentations like miso and shochu, is hypothesized to have been domesticated from Aspergillus parasiticus. Here, we examined the effects of long-term A. sojae use in soybean fermentation on population structure, genome variation, and phenotypic traits. We analyzed 17 A. sojae and 24 A. parasiticus genomes (23 of which were sequenced for this study), alongside phenotypic traits of 9 isolates. Aspergillus sojae formed a distinct, low-diversity population, suggesting a recent clonal expansion. Interestingly, a population of A. parasiticus was more closely related to A. sojae than other A. parasiticus populations. Genome comparisons revealed loss-of-function mutations in A. sojae, notably in biosynthetic gene clusters encoding secondary metabolites, including the aflatoxin cluster. Interestingly though, A. sojae harbored a partial duplication of a siderophore biosynthetic cluster. Phenotypic assays showed A. sojae lacked aflatoxin production, while it was variable in A. parasiticus isolates. Additionally, certain A. sojae strains exhibited larger colony diameters under miso-like salt conditions. These findings support the hypothesis that A. parasiticus is the progenitor of A. sojae and that domestication significantly reduced genetic diversity. Future research should explore how wild and food-associated strains influence sensory attributes and microbial community dynamics in fermented soy products.

A microscopy-based image analysis pipeline for the quantification of germination of filamentous fungi

Germination is the fundamental process whereby fungi transition from the dormant and stress resistant spores into actively replicating cells such as hyphae. Germination is essential for fungal colonization of new environments and pathogenesis, yet this differentiation process remains relatively poorly understood. For filamentous fungi, the study of germination has been limited by the lack of high-throughput, temporal, low cost, and easy-to-use methods of quantifying germination. To this end we have developed an image analysis pipeline to automate the quantification of germination from microscopy images. We have optimized this tool for the fungal pathogen Aspergillus fumigatus and demonstrated its potential applications by evaluating different strains, germination inhibitors, and auxotrophic and antifungal resistant mutants. Finally, we have expanded this tool to a variety of filamentous fungi and developed an easy-to-use web app for the fungal research community.

Species Identification and Orthologous Allergen Prediction and Expression in the Genus Aspergillus

The genus Aspergillus comprises a diverse group of fungi that can cause a range of health issues, including systemic infections and allergic reactions. In this regard, A. fumigatus has been recognized as the most prevalent allergen-producing species. This genus taxonomic classification has been subject to frequent updates, which has generated considerable difficulties for its classification when traditional identification methodologies are employed. To demonstrate the feasibility of this approach, we sequenced the whole genomes of 81 Aspergillus isolates and evaluated a WGS-based pipeline for precise species identification. This pipeline employed two methodologies: (i) BLASTn web using four barcode genes and (ii) species tree inference by OrthoFinder. Furthermore, we conducted a prediction of allergenic capacity based on a homology analysis across all the isolated species and confirmed by RT-qPCR the expression of three orthologous allergens (Asp f 1, Asp f 3 and Asp f 22) in fifteen different Aspergillus species. The species-level identification rate with the barcoding and the species tree were calculated at 64.2% and 100%, respectively. The results demonstrated that A. fumigatus, A. flavus and A. niger were the most prevalent species. The species A. hortae, A. uvarum, A. spinulosporus, A. sydowii, A. westerdijkiae, A. amoenus and A. rhizopodus identified in this study represent the inaugural report of their presence in our region. The results of the homology analysis indicated the presence of orthologous allergens in a wide range of non-fumigatus species. This study presents a novel approach based on WGS that enables the classification of new species within the genus Aspergillus and reports the genomic sequences of a great diversity of species isolated in our geographic area that had never been reported before. Additionally, this approach enables the prediction of allergens in species other than A. fumigatus and demonstrates their genetic expression, thereby contributing to the understanding of the allergenic potential of different species within this fungal genus.

Different metabolite profiles across Penicillium roqueforti populations associated with ecological niche specialisation and domestication

Fungi are known to produce many chemically diversified metabolites, yet their ecological roles are not always fully understood. The blue cheese fungus Penicillium roqueforti thrives in different ecological niches and is known to produce a wide range of metabolites, including mycotoxins. Three P. roqueforti populations have been domesticated for cheese production and two populations thrive in other anthropized environments, i.e., food, lumber and silage. In this study, we looked for differences in targeted and untargeted metabolite production profiles between populations using HPLC-HR-Q-TOF and UHPLC-Q-TOF-HR-MS/MS. The non-cheese populations produced several fatty acids and different terpenoids, lacking in cheese strains. The Termignon cheese population displayed intermediate metabolite profiles between cheese and non-cheese populations, as previously shown for other traits. The non-Roquefort cheese population with the strongest domestication syndrome, produced the lowest quantities of measured metabolites, including mycophenolic acid (MPA), andrastin A and PR toxin. Its inability to produce MPA was due to a deletion in the mpaC gene, while a premature stop codon in ORF 11 of the PR toxin gene cluster explained PR toxin absence and the accumulation of its intermediates, i.e., eremofortins A and B. In the Roquefort population, we detected no PR toxin nor eremofortins A or B, but found no indel or frameshift mutation, suggesting downregulation. The hypotoxigenic trait of domesticated cheese populations can be hypothesized to be linked to the loss of this ability through trait degeneration and/or the selection of low toxin producers. It may also be due to the fact that populations from other anthropized environments maintained high metabolite diversity as the bioactivities of these compounds are likely important in these ecological niches.

Isolation and characterization of koji mold (Aspergillus oryzae) from nature in Niigata

Koji mold (Aspergillus oryzae) is a key microorganism in brewing and fermentation in Japan. We isolated koji molds from the environment in Niigata Prefecture. Eighty-one environmental samples were placed on isolation medium made from steamed rice with wood ash and 36 Aspergillus section Flavi-like strains were obtained. Of those, 26 strains did not produce aflatoxin. We studied their morphology, sequence of ITS region, calmodulin gene, aflatoxin biosynthetic homologous gene cluster and α-amylase gene and fermentation-related enzyme activities. Furthermore, DNA-seq analysis of 14 strains from 26 non-aflatoxin producing strains were conducted and compared the three mycotoxin biosynthetic gene clusters (aflatoxin, cyclopiazonic acid, and aflatrem) and fermentation-related genes against those of reference strain A. oryzae RIB40. In some strains, gene sequences confirmed the absence of mycotoxin production, but differences in fermentation-related enzyme activities could not be explained well by amino acid substitutions. We classified the 26 isolates into 6 morphology types based on the appearance of colonies and mating types, and it was found that strains of the same morphology type had similar enzymatic profiles and gene sequences. Our results show that koji molds with various properties occur in the environment, and it will expand the possibilities of koji mold in industrial use.

Genomic Localization Bias of Secondary Metabolite Gene Clusters and Association with Histone Modifications in Aspergillus

Fungi are well-known producers of bioactive secondary metabolites (SMs), which have been exploited for decades by humankind for various medical applications like therapeutics and antibiotics. SMs are synthesized by biosynthetic gene clusters (BGCs)-physically co-localized and co-regulated genes. Because BGCs are often regulated by histone post-translational modifications (PTMs), it was suggested that their chromosomal location is important for their expression. Studies in a few fungal species indicated an enrichment of BGCs in sub-telomeric regions; however, there is no evidence that BGCs with distinct genomic localization are regulated by different histone PTMs. Here, we used 174 Aspergillus species covering 22 sections to determine the correlation between BGC genomic localization, gene expression, and histone PTMs. We found a high abundance and diversity of SM backbone genes across the Aspergillus genus, with notable unique genes within sections. Being unique or conserved in many species, BGCs showed a strong bias for being localized in low-synteny regions, regardless of their position in chromosomes. Using chromosome-level assemblies, we also confirmed a significantly biased localization in sub-telomeric regions. Notably, SM backbone genes in sub-telomeric regions and about half of those in low-synteny regions exhibit higher gene expression variability, likely due to the similar higher variability in H3K4me3 and H3K36me3 histone PTMs; while variations in histone H3 acetylation and H3K9me3 are not correlated to genomic localization and expression variation, as analyzed in two Aspergillus species. Expression variability across four Aspergillus species further supports that BGCs tend to be located in low-synteny regions and that regulation of expression in those regions likely involves different histone PTMs than the most commonly studied modifications.

Genomic and phenotypic imprints of microbial domestication on cheese starter cultures

Domestication - the artificial selection of wild species to obtain variants with traits of human interest - was integral to the rise of complex societies. The oversupply of food was probably associated with the formalization of food preservation strategies through microbial fermentation. While considerable literature exists on the antiquity of fermented food, only few eukaryotic microbes have been studied so far for signs of domestication, less is known for bacteria. Here, we tested if cheese starter cultures harbour typical hallmarks of domestication by characterising over 100 community samples and over 100 individual strains isolated from historical and modern traditional Swiss cheese starter cultures. We find that cheese starter cultures have low genetic diversity both at the species and strain-level and maintained stable phenotypic traits. Molecular clock dating further suggests that the evolutionary origin of the bacteria approximately coincided with the first archaeological records of cheese making. Finally, we find evidence for ongoing genome decay and pseudogenization via transposon insertion related to a reduction of their niche breadth. Future work documenting the prevalence of these hallmarks across diverse fermented food systems and geographic regions will be key to unveiling the joint history of humanity with fermented food microbes.

Neurospora intermedia from a traditional fermented food enables waste-to-food conversion

Fungal fermentation of food and agricultural by-products holds promise for improving food sustainability and security. However, the molecular basis of fungal waste-to-food upcycling remains poorly understood. Here we use a multi-omics approach to characterize oncom, a fermented food traditionally produced from soymilk by-products in Java, Indonesia. Metagenomic sequencing of samples from small-scale producers in Western Java indicated that the fungus Neurospora intermedia dominates oncom. Further transcriptomic, metabolomic and phylogenomic analysis revealed that oncom-derived N. intermedia utilizes pectin and cellulose degradation during fermentation and belongs to a genetically distinct subpopulation associated with human-generated by-products. Finally, we found that N. intermedia grew on diverse by-products such as fruit and vegetable pomace and plant-based milk waste, did not encode mycotoxins, and could create foods that were positively perceived by consumers outside Indonesia. These results showcase the traditional significance and future potential of fungal fermentation for creating delicious and nutritious foods from readily available by-products.

GastronOmics: Edibility and safety of mycelium of the oyster mushroom Pleurotus ostreatus

Food production is one of the most environmentally damaging human activities. In the face of climate change, it is essential to rethink our dietary habits and explore potential alternative foods catering both towards human and planetary needs. Fungal mycelium might be an attractive alternative protein source due to its rapid growth on sustainable substrates as well as promising nutritional and organoleptic properties. The natural biodiversity of filamentous fungi is vast and represents an untapped reservoir for food innovation. However, fungi are known to produce bioactive compounds that may affect human health, both positively and negatively. To narrow the search for safe and culinarily attractive fungal species, mycelia of edible fruiting-body forming fungi provide a promising starting point. Here, we explore whether the culinary attractiveness and safety of the commonly eaten mushroom, Pleurotus ostreatus, can also be translated to its mycelium. Whole-genome sequencing and pan-genome analysis revealed a high degree of genetic variability within the genus Pleurotus, suggesting that gastronomic traits as well as food safety may differ between strains. A representative strain, P. ostreatus M2191, was further analyzed for the food safety, nutritional properties and culinary applicability of its mycelium. No regulated mycotoxins were detected in either the fruiting body nor the mycelium. Yet, P. ostreatus is known to produce four peptide toxins, Ostreatin, Ostreolysin and Pleurotoysin A/B. These were found to be lower in the mycelium compared to fruiting bodies, which are already considered safe for consumption. Instead, a number of secondary metabolites with potential health benefits were detected in the fungal mycelium. In silico analysis of the proteome suggested low allergenicity. In addition, the fruiting body and the mycelium showed similar nutritional value, which was dependent on the growth substrate. To highlight the culinary potential of mycelium, we created a dish served at the two-star restaurant the Alchemist in Copenhagen, Denmark. Sensory analysis of the mycelium dish by an untrained consumer panel indicated consumer liking and openness to fungal mycelia. Based on sustainability, safety, culinary potential, and consumer acceptance, our findings suggest that P. ostreatus mycelium has great potential for use as a novel food source.

Comparative Genome Analysis of Japanese Field-Isolated Aspergillus for Aflatoxin Productivity and Non-Productivity

Aspergillus flavus produces aflatoxin, a carcinogenic fungal toxin that poses a threat to the agricultural and food industries. There is a concern that the distribution of aflatoxin-producing A. flavus is expanding in Japan due to climate change, and it is necessary to understand what types of strains inhabit. In this study, we sequenced the genomes of four Aspergillus strains isolated from agricultural fields in the Ibaraki prefecture of Japan and identified their genetic variants. Phylogenetic analysis based on single-nucleotide variants revealed that the two aflatoxin-producing strains were closely related to A. flavus NRRL3357, whereas the two non-producing strains were closely related to the RIB40 strain of Aspergillus oryzae, a fungus widely used in the Japanese fermentation industry. A detailed analysis of the variants in the aflatoxin biosynthetic gene cluster showed that the two aflatoxin-producing strains belonged to different morphotype lineages. RT-qPCR results indicated that the expression of aflatoxin biosynthetic genes was consistent with aflatoxin production in the two aflatoxin-producing strains, whereas the two non-producing strains expressed most of the aflatoxin biosynthetic genes, unlike common knowledge in A. oryzae, suggesting that the lack of aflatoxin production was attributed to genes outside of the aflatoxin biosynthetic gene cluster in these strains.

Starter Culture Development and Innovation for Novel Fermented Foods

Interest in fermented foods is increasing because fermented foods are promising solutions for more secure food systems with an increased proportion of minimally processed plant foods and a smaller environmental footprint. These developments also pertain to novel fermented food for which no traditional template exists, raising the question of how to develop starter cultures for such fermentations. This review establishes a framework that integrates traditional and scientific knowledge systems for the selection of suitable cultures. Safety considerations, the use of organisms in traditional food fermentations, and the link of phylogeny to metabolic properties provide criteria for culture selection. Such approaches can also select for microbial strains that have health benefits. A science-based approach to the development of novel fermented foods can substantially advance their value through more secure food systems, food products that provide health-promoting microbes, and the provision of foods that improve human health.

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.

Phylogenomics reveals extensive misidentification of fungal strains from the genus Aspergillus

Modern taxonomic classification is often based on phylogenetic analyses of a few molecular markers, although single-gene studies are still common. Here, we leverage genome-scale molecular phylogenetics (phylogenomics) of species and populations to reconstruct evolutionary relationships in a dense data set of 710 fungal genomes from the biomedically and technologically important genus Aspergillus. To do so, we generated a novel set of 1,362 high-quality molecular markers specific for Aspergillus and provided profile Hidden Markov Models for each, facilitating their use by others. Examining the resulting phylogeny helped resolve ongoing taxonomic controversies, identified new ones, and revealed extensive strain misidentification (7.59% of strains were previously misidentified), underscoring the importance of population-level sampling in species classification. These findings were corroborated using the current standard, taxonomically informative loci. These findings suggest that phylogenomics of species and populations can facilitate accurate taxonomic classifications and reconstructions of the Tree of Life.IMPORTANCEIdentification of fungal species relies on the use of molecular markers. Advances in genomic technologies have made it possible to sequence the genome of any fungal strain, making it possible to use genomic data for the accurate assignment of strains to fungal species (and for the discovery of new ones). We examined the usefulness and current limitations of genomic data using a large data set of 710 publicly available genomes from multiple strains and species of the biomedically, agriculturally, and industrially important genus Aspergillus. Our evolutionary genomic analyses revealed that nearly 8% of publicly available Aspergillus genomes are misidentified. Our work highlights the usefulness of genomic data for fungal systematic biology and suggests that systematic genome sequencing of multiple strains, including reference strains (e.g., type strains), of fungal species will be required to reduce misidentification errors in public databases.

An Aspergillus flavus strain from bee bread of the Western honey bee (Apis mellifera) displays adaptations to distinctive features of the hive environment

Aspergillus fungi are ubiquitous inhabitants of colonies of the western honey bee (Apis mellifera), where they interact with bees in associations ranging from parasitism to possible mutualism. Aspergillus Flavi fungi are frequently found in bee bread (pollen processed for longterm storage) and are thought to contribute to food preparation, processing, preservation, and digestion. Conditions in the hive are challenging for fungi due, in part, to xeric and acidic properties of bee bread and the omnipresence of propolis, an antimicrobial product manufactured by bees from plant resins. We used quantitative and qualitative assays to determine whether A. flavus isolated from bee bread demonstrates tolerance for hive environmental conditions in terms of temperature, pH, osmotic pressure, and propolis exposure. Comparisons made use of three strains of A. flavus: a fungal biocontrol product not known from beehives (AF36), a strain isolated from bee bread (AFBB) in hives from central Illinois, and a pathogenic strain from a honey bee colony displaying symptoms of stonebrood (AFPA). Strain AFBB displayed higher tolerance of acidic conditions, low matric potential (simulating xeric substrate), and propolis exposure than did other strains. A genomic comparison between this new strain and the reference NRRL-3357 showed that AFBB, like AF36, might be blocked from carrying out aflatoxin biosynthesis. Sequence comparisons also revealed several missense variants in genes that encode proteins regulating osmotolerance and osmotic pressure in Aspergillus spp., including SakA, SskB, GfdA, and TcsB/Sln1. Collectively, results of our laboratory assays and genetic analyses are consistent with the suggestion that the strain isolated from bee bread is adapted to the bee bread environment and may have persisted due to a coevolutionary relationship between Aspergillus and A. mellifera. This finding bolsters recent concerns about the effects of fungicide use near bee colonies and broadens the ecological importance of highly adaptable fungal strains.

Domestication through clandestine cultivation constrained genetic diversity in magic mushrooms relative to naturalized populations

Fungi that are edible or fermentative were domesticated through selective cultivation of their desired traits. Domestication is often associated with inbreeding or selfing, which may fix traits other than those under selection, and causes an overall decrease in heterozygosity. A hallucinogenic mushroom, Psilocybe cubensis, was domesticated from its niche in livestock dung for production of psilocybin. It has caused accidental poisonings since the 1940s in Australia, which is a population hypothesized to be introduced from an unknown center of origin. We sequenced genomes of 38 isolates from Australia and compared them with 86 genomes of commercially available cultivars to determine (1) whether P. cubensis was introduced to Australia, and (2) how domestication has impacted commercial cultivars. Our analyses of genome-wide SNPs and single-copy orthologs showed that the Australian population is naturalized, having recovered its effective population size after a bottleneck when it was introduced, and it has maintained relatively high genetic diversity based on measures of nucleotide and allelic diversity. In contrast, domesticated cultivars generally have low effective population sizes and hallmarks of selfing and clonal propagation, including low genetic diversity, low heterozygosity, high linkage disequilibrium, and low allelic diversity of mating-compatibility genes. Analyses of kinship show that most cultivars are founded from related populations. Alleles in the psilocybin gene cluster are identical across most cultivars of P. cubensis with low diversity across coding sequence; however, unique allelic diversity in Australia and some cultivars may translate to differences in biosynthesis of psilocybin and its analogs.

Microbiome Assembly in Fermented Foods

For thousands of years, humans have enjoyed the novel flavors, increased shelf-life, and nutritional benefits that microbes provide in fermented foods and beverages. Recent sequencing surveys of ferments have mapped patterns of microbial diversity across space, time, and production practices. But a mechanistic understanding of how fermented food microbiomes assemble has only recently begun to emerge. Using three foods as case studies (surface-ripened cheese, sourdough starters, and fermented vegetables), we use an ecological and evolutionary framework to identify how microbial communities assemble in ferments. By combining in situ sequencing surveys with in vitro models, we are beginning to understand how dispersal, selection, diversification, and drift generate the diversity of fermented food communities. Most food producers are unaware of the ecological processes occurring in their production environments, but the theory and models of ecology and evolution can provide new approaches for managing fermented food microbiomes, from farm to ferment.

Bacterial-fungal interactions promote parallel evolution of global transcriptional regulators in a widespread Staphylococcus species

Experimental studies of microbial evolution have largely focused on monocultures of model organisms, but most microbes live in communities where interactions with other species may impact rates and modes of evolution. Using the cheese rind model microbial community, we determined how species interactions shape the evolution of the widespread food- and animal-associated bacterium Staphylococcus xylosus. We evolved S. xylosus for 450 generations alone or in co-culture with one of three microbes: the yeast Debaryomyces hansenii, the bacterium Brevibacterium aurantiacum, and the mold Penicillium solitum. We used the frequency of colony morphology mutants (pigment and colony texture phenotypes) and whole-genome sequencing of isolates to quantify phenotypic and genomic evolution. The yeast D. hansenii strongly promoted diversification of S. xylosus. By the end of the experiment, all populations co-cultured with the yeast were dominated by pigment and colony morphology mutant phenotypes. Populations of S. xylosus grown alone, with B. aurantiacum, or with P. solitum did not evolve novel phenotypic diversity. Whole-genome sequencing of individual mutant isolates across all four treatments identified numerous unique mutations in the operons for the SigB, Agr, and WalRK global regulators, but only in the D. hansenii treatment. Phenotyping and RNA-seq experiments highlighted altered pigment and biofilm production, spreading, stress tolerance, and metabolism of S. xylosus mutants. Fitness experiments revealed antagonistic pleiotropy, where beneficial mutations that evolved in the presence of the yeast had strong negative fitness effects in other biotic environments. This work demonstrates that bacterial-fungal interactions can have long-term evolutionary consequences within multispecies microbiomes by facilitating the evolution of strain diversity.

Domestication in dry-cured meat Penicillium fungi: Convergent specific phenotypes and horizontal gene transfers without strong genetic subdivision

Some fungi have been domesticated for food production, with genetic differentiation between populations from food and wild environments, and food populations often acquiring beneficial traits through horizontal gene transfers (HGTs). Studying their adaptation to human-made substrates is of fundamental and applied importance for understanding adaptation processes and for further strain improvement. We studied here the population structures and phenotypes of two distantly related Penicillium species used for dry-cured meat production, P. nalgiovense, the most common species in the dry-cured meat food industry, and P. salamii, used locally by farms. Both species displayed low genetic diversity, lacking differentiation between strains isolated from dry-cured meat and those from other environments. Nevertheless, the strains collected from dry-cured meat within each species displayed slower proteolysis and lipolysis than their wild conspecifics, and those of P. nalgiovense were whiter. Phenotypically, the non-dry-cured meat strains were more similar to their sister species than to their conspecific dry-cured meat strains, indicating an evolution of specific phenotypes in dry-cured meat strains. A comparison of available Penicillium genomes from various environments revealed HGTs, particularly between P. nalgiovense and P. salamii (representing almost 1.5 Mb of cumulative length). HGTs additionally involved P. biforme, also found in dry-cured meat products. We further detected positive selection based on amino acid changes. Our findings suggest that selection by humans has shaped the P. salamii and P. nalgiovense populations used for dry-cured meat production, which constitutes domestication. Several genetic and phenotypic changes were similar in P. salamii, P. nalgiovense and P. biforme, indicating convergent adaptation to the same human-made environment. Our findings have implications for fundamental knowledge on adaptation and for the food industry: the discovery of different phenotypes and of two mating types paves the way for strain improvement by conventional breeding, to elucidate the genomic bases of beneficial phenotypes and to generate diversity.

A new cheese population in Penicillium roqueforti and adaptation of the five populations to their ecological niche

Domestication is an excellent case study for understanding adaptation and multiple fungal lineages have been domesticated for fermenting food products. Studying domestication in fungi has thus both fundamental and applied interest. Genomic studies have revealed the existence of four populations within the blue-cheese-making fungus Penicillium roqueforti. The two cheese populations show footprints of domestication, but the adaptation of the two non-cheese populations to their ecological niches (i.e., silage/spoiled food and lumber/spoiled food) has not been investigated yet. Here, we reveal the existence of a new P. roqueforti population, specific to French Termignon cheeses, produced using small-scale traditional practices, with spontaneous blue mould colonisation. This Termignon population is genetically differentiated from the four previously identified populations, providing a novel source of genetic diversity for cheese making. The Termignon population indeed displayed substantial genetic diversity, both mating types, horizontally transferred regions previously detected in the non-Roquefort population, and intermediate phenotypes between cheese and non-cheese populations. Phenotypically, the non-Roquefort cheese population was the most differentiated, with specific traits beneficial for cheese making, in particular higher tolerance to salt, to acidic pH and to lactic acid. Our results support the view that this clonal population, used for many cheese types in multiple countries, is a domesticated lineage on which humans exerted strong selection. The lumber/spoiled food and silage/spoiled food populations were not more tolerant to crop fungicides but showed faster growth in various carbon sources (e.g., dextrose, pectin, sucrose, xylose and/or lactose), which can be beneficial in their ecological niches. Such contrasted phenotypes between P. roqueforti populations, with beneficial traits for cheese-making in the cheese populations and enhanced ability to metabolise sugars in the lumber/spoiled food population, support the inference of domestication in cheese fungi and more generally of adaptation to anthropized environments.

Composition and activity of antifungal lipopeptides produced by Bacillus spp. in daqu fermentation

Daqu is a solid-state fermentation and saccharification starter for the Chinese liquor baijou. During the daqu stage, amylolytic and proteolytic enzymes are produced by Bacillus and fungi. Bacillus spp. also produce lipopeptides with a broad spectrum of antimicrobial activities but direct evidence for their impact on community assembly in daqu is lacking. This study aimed to study the interaction between Bacillus spp. and fungi in daqu models. The antifungal activity of surfactin, fengycin, and iturin A was initially assessed in vitro. Iturin A displayed the strongest antifungal activity (MIC = 10-50 mg/L). In situ antifungal activity of B. amyloliquefaciens and B. velezensis against molds was observed in a simple daqu model inoculated with single strains of Bacillus species. Formation of lipopeptides in situ was supported by quantification of mRNA encoding for enzymes for surfactin, fengycin, and iturin A biosynthesis. In situ antifungal activity of Bacillus species was also observed in a complex daqu model that was inoculated with 8 bacterial or fungal strains plus one of the three strains of Bacillus. A relationship of lipopeptides to in situ antifungal activity was further supported by detection of the lipopeptides by liquid chromatography coupled to mass spectrometry. Both results indicated that B velezensis FUA2155 had higher antifungal activity in the daqu model, and was the only strain that produced multiple iturin A congeners in situ. Taken together, this study provides evidence that production of lipopeptides by Bacillus species in daqu may impact community assembly and hence product quality.

Trace copper-mediated asexual development via a superoxide dismutase and induction of AobrlA in Aspergillus oryzae

The filamentous fungus Aspergillus oryzae, in which sexual reproduction remains to be discovered, proliferates mainly via asexual spores (conidia). Therefore, despite its industrial importance in food fermentation and recombinant protein production, breeding beneficial strains by genetic crosses is difficult. In Aspergillus flavus, which is genetically close to A. oryzae, structures known as sclerotia are formed asexually, but they are also related to sexual development. Sclerotia are observed in some A. oryzae strains, although no sclerotia formation has been reported in most strains. A better understanding of the regulatory mechanisms underlying sclerotia formation in A. oryzae may contribute to discover its sexual development. Some factors involved in sclerotia formation have been previously identified, but their regulatory mechanisms have not been well studied in A. oryzae. In this study, we found that copper strongly inhibited sclerotia formation and induced conidiation. Deletion of AobrlA encoding a core regulator of conidiation and ecdR involved in transcriptional induction of AobrlA suppressed the copper-mediated inhibition of sclerotia formation, suggesting that AobrlA induction in response to copper leads not only to conidiation but also to inhibition of sclerotia formation. In addition, deletion of the copper-dependent superoxide dismutase (SOD) gene and its copper chaperone gene partially suppressed such copper-mediated induction of conidiation and inhibition of sclerotia formation, indicating that copper regulates asexual development via the copper-dependent SOD. Taken together, our results demonstrate that copper regulates asexual development, such as sclerotia formation and conidiation, via the copper-dependent SOD and transcriptional induction of AobrlA in A. oryzae.

Contrasting Genomic Evolution Between Domesticated and Wild Kluyveromyces lactis Yeast Populations

The process of domestication has variable consequences on genome evolution leading to different phenotypic signatures. Access to the complete genome sequences of a large number of individuals makes it possible to explore the different facets of this domestication process. Here, we sought to explore the genome evolution of Kluyveromyces lactis, a yeast species well known for its involvement in dairy processes and also present in natural environments. Using a combination of short- and long-read sequencing strategies, we investigated the genomic variability of 41 K. lactis isolates and found that the overall genetic diversity of this species is very high (θw = 3.3 × 10-2) compared with other species such as Saccharomyces cerevisiae (θw = 1.6 × 10-2). However, the domesticated dairy population shows a reduced level of diversity (θw = 1 × 10-3), probably due to a domestication bottleneck. In addition, this entire population is characterized by the introgression of the LAC4 and LAC12 genes, responsible for lactose fermentation and coming from the closely related species, Kluyveromyces marxianus, as previously described. Our results highlighted that the LAC4/LAC12 gene cluster was acquired through multiple and independent introgression events. Finally, we also identified several genes that could play a role in adaptation to dairy environments through copy number variation. These genes are involved in sugar consumption, flocculation, and drug resistance, and may play a role in dairy processes. Overall, our study illustrates contrasting genomic evolution and sheds new light on the impact of domestication processes on it.

Linolenic Acid-Derived Oxylipins Inhibit Aflatoxin Biosynthesis in Aspergillus flavus through Activation of Imizoquin Biosynthesis

Oxylipins play important signaling roles in aflatoxin (AF) biosynthesis in Aspergillus flavus. We previously showed that exogenous supply of autoxidated linolenic acid (AL) inhibited AF biosynthesis in A. flavus via oxylipins, but the molecular mechanism is still unknown. Here, we performed multiomics analyses of A. flavus grown in media with or without AL. Targeted metabolite analyses and quantitative reverse transcription (qRT)-polymerase chain reaction (PCR) showed that the imizoquin (IMQ) biosynthetic pathway was distinctly upregulated in the presence of AL. 13C-glucose labeling confirmed in parallel that the tricarboxylic acid cycle was also enhanced by AL, consistent with observed increases in mycelial growth. Moreover, we integrated thermal proteome profiling and molecular dynamics simulations to identify a potential receptor of AL; AL was found to interact with a transporter (ImqJ) located in the IMQ gene cluster, primarily through hydrophobic interactions. Further analyses of strains with an IMQ pathway transcription factor overexpressed or knocked out confirmed that this pathway was critical for AL-mediated inhibition of AF biosynthesis. Comparison of 22 assembled A. flavus and Aspergillus oryzae genomes showed that genes involved in the IMQ pathway were positively selected in A. oryzae. Taken together, the results of our study provide novel insights into oxylipin-mediated regulation of AF biosynthesis and suggest potential methods for preventing AF contamination of crops.

Convergence in domesticated fungi used for cheese and dry-cured meat maturation: beneficial traits, genomic mechanisms, and degeneration

Humans have domesticated genetically distant fungi for similar uses, the fermentation of lipid-rich and sugar-rich food to generate attractive aspects, odor and aroma, and to improve shelf life and product safety. Multiple independent domestication events also occurred within species. We review recent evidence of phenotypic convergence during the domestication of fungi for making cheese (Saccharomyces cerevisiae, Penicillium roqueforti, P. camemberti, and Geotrichum candidum) and for dry-cured meat making (P. nalgiovense and P. salamii). Convergence following adaptation to similar ecological niches involved colony aspect (fluffiness and color), lipolysis, proteolysis, volatile compound production, and competitive ability against food spoilers. We review evidence for convergence in genetic diversity loss in domesticated populations and in the degeneration of unused traits, such as toxin production and sexual reproduction. Phenotypic convergence sometimes occurred by similar mechanisms of genomic adaptation, in particular horizontal gene transfers and loss of genes.

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.

Genomic and Molecular Identification of Genes Contributing to the Caspofungin Paradoxical Effect in Aspergillus fumigatus

Aspergillus fumigatus is a deadly opportunistic fungal pathogen responsible for ~100,000 annual deaths. Azoles are the first line antifungal agent used against A. fumigatus, but azole resistance has rapidly evolved making treatment challenging. Caspofungin is an important second-line therapy against invasive pulmonary aspergillosis, a severe A. fumigatus infection. Caspofungin functions by inhibiting β-1,3-glucan synthesis, a primary and essential component of the fungal cell wall. A phenomenon termed the caspofungin paradoxical effect (CPE) has been observed in several fungal species where at higher concentrations of caspofungin, chitin replaces β-1,3-glucan, morphology returns to normal, and growth rate increases. CPE appears to occur in vivo, and it is therefore clinically important to better understand the genetic contributors to CPE. We applied genomewide association (GWA) analysis and molecular genetics to identify and validate candidate genes involved in CPE. We quantified CPE across 67 clinical isolates and conducted three independent GWA analyses to identify genetic variants associated with CPE. We identified 48 single nucleotide polymorphisms (SNPs) associated with CPE. We used a CRISPR/Cas9 approach to generate gene deletion mutants for seven genes harboring candidate SNPs. Two null mutants, ΔAfu3g13230 and ΔAfu4g07080 (dscP), resulted in reduced basal growth rate and a loss of CPE. We further characterized the dscP phosphatase-null mutant and observed a significant reduction in conidia production and extremely high sensitivity to caspofungin at both low and high concentrations. Collectively, our work reveals the contribution of Afu3g13230 and dscP in CPE and sheds new light on the complex genetic interactions governing this phenotype. IMPORTANCE This is one of the first studies to apply genomewide association (GWA) analysis to identify genes involved in an Aspergillus fumigatus phenotype. A. fumigatus is an opportunistic fungal pathogen that causes hundreds of thousands of infections and ~100,000 deaths each year, and antifungal resistance has rapidly evolved in this species. A phenomenon called the caspofungin paradoxical effect (CPE) occurs in some isolates, where high concentrations of the drug lead to increased growth rate. There is clinical relevance in understanding the genetic basis of this phenotype, since caspofungin concentrations could lead to unintended adverse clinical outcomes in certain cases. Using GWA analysis, we identified several interesting candidate polymorphisms and genes and then generated gene deletion mutants to determine whether these genes were important for CPE. Two of these mutant strains (ΔAfu3g13230 and ΔAfu4g07080/ΔdscP) displayed a loss of the CPE. This study sheds light on the genes involved in clinically important phenotype CPE.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Low-cost, specific PCR assays to identify the main aflatoxigenic species of Aspergillus section Flavi

Aflatoxins are fungal metabolites that are present as contaminants in food globally. Most aflatoxigenic species belong to Aspergillus section Flavi, and the main ones are grouped in the A. flavus clade, where many cryptic species that are difficult to discriminate are found. In this study, we investigated inter- and intraspecific diversity of the A. flavus clade to develop low-cost, species-specific PCR assays for identifying aflatoxigenic species. A total of 269 sequences of the second largest subunit of RNA polymerase II (RPB2) locus were retrieved from GenBank, and primer pairs were designed using data mining to identify A. flavus, A. parasiticus, and A. novoparasiticus. Species-specific amplicons of approximately 620, 350, and 860 bp enabled identification of target species as A. flavus, A. parasiticus, and A. novoparasiticus, respectively.

The periodic table of fermented foods: limitations and opportunities

Fermentation is one of the oldest methods of food processing and accounts for a substantial proportion of human foods, including not only staple foods such as bread, cereal porridges or fermented legumes but also fermented vegetables, meats, fish and dairy, alcoholic beverages as well as coffee, cocoa and condiments such as vinegar, soy sauce and fish sauces. Adding the regional varieties to these diverse product categories makes for an almost immeasurable diversity of fermented foods. The periodic table of fermented foods aims to map this diversity on the 118 entries of the periodic table of chemical elements. While the table fails to represent the diversity of fermented foods, it represents major fermentation substrates, product categories, fermentation processes and fermentation organisms. This communication not only addresses limitations of the graphical display on a "periodic table of fermented foods", but also identifies opportunities that relate to questions that are facilitated by this graphical presentation: on the origin and purpose of food fermentation, which fermented foods represent "indigenous" foods, differences and similarities in the assembly of microbial communities in different fermentations, differences in the global preferences for food fermentation, the link between microbial diversity, fermentation time and product properties, and opportunities of using traditional food fermentations as template for development of new products. KEY POINTS: • Fermented foods are produced in an almost immeasurable diversity. • Fermented foods were mapped on a periodic table of fermented foods. • This table facilitates identification of communalities and differences of products.

A glimpse at an early stage of microbe domestication revealed in the variable genome of Torulaspora delbrueckii, an emergent industrial yeast

Microbe domestication has a major applied relevance but is still poorly understood from an evolutionary perspective. The yeast Torulaspora delbrueckii is gaining importance for biotechnology but little is known about its population structure, variation in gene content, or possible domestication routes. Here, we show that T. delbrueckii is composed of five major clades. Among the three European clades, a lineage associated with the wild arboreal niche is sister to the two other lineages that are linked with anthropic environments, one to wine fermentations and the other to diverse sources including dairy products and bread dough (Mix- Anthropic clade). Using 64 genomes we assembled the pangenome and the variable genome of T. delbrueckii. A comparison with Saccharomyces cerevisiae indicated that the weight of the variable genome in the pangenome of T. delbrueckii is considerably smaller. An association of gene content and ecology supported the hypothesis that the Mix - Anthropic clade has the most specialized genome and indicated that some of the exclusive genes were implicated in galactose and maltose utilization. More detailed analyses traced the acquisition of a cluster of GAL genes in strains associated with dairy products and the expansion and functional diversification of MAL genes in strains isolated from bread dough. Contrary to S. cerevisiae, domestication in T. delbrueckii is not primarily driven by alcoholic fermentation but rather by adaptation to dairy and bread-production niches. This study expands our views on the processes of microbe domestication and on the trajectories leading to adaptation to anthropic niches.

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.

Regulation of gliotoxin biosynthesis and protection in Aspergillus species

Aspergillus fumigatus causes a range of human and animal diseases collectively known as aspergillosis. A. fumigatus possesses and expresses a range of genetic determinants of virulence, which facilitate colonisation and disease progression, including the secretion of mycotoxins. Gliotoxin (GT) is the best studied A. fumigatus mycotoxin with a wide range of known toxic effects that impair human immune cell function. GT is also highly toxic to A. fumigatus and this fungus has evolved self-protection mechanisms that include (i) the GT efflux pump GliA, (ii) the GT neutralising enzyme GliT, and (iii) the negative regulation of GT biosynthesis by the bis-thiomethyltransferase GtmA. The transcription factor (TF) RglT is the main regulator of GliT and this GT protection mechanism also occurs in the non-GT producing fungus A. nidulans. However, the A. nidulans genome does not encode GtmA and GliA. This work aimed at analysing the transcriptional response to exogenous GT in A. fumigatus and A. nidulans, two distantly related Aspergillus species, and to identify additional components required for GT protection. RNA-sequencing shows a highly different transcriptional response to exogenous GT with the RglT-dependent regulon also significantly differing between A. fumigatus and A. nidulans. However, we were able to observe homologs whose expression pattern was similar in both species (43 RglT-independent and 11 RglT-dependent). Based on this approach, we identified a novel RglT-dependent methyltranferase, MtrA, involved in GT protection. Taking into consideration the occurrence of RglT-independent modulated genes, we screened an A. fumigatus deletion library of 484 transcription factors (TFs) for sensitivity to GT and identified 15 TFs important for GT self-protection. Of these, the TF KojR, which is essential for kojic acid biosynthesis in Aspergillus oryzae, was also essential for virulence and GT biosynthesis in A. fumigatus, and for GT protection in A. fumigatus, A. nidulans, and A. oryzae. KojR regulates rglT, gliT, gliJ expression and sulfur metabolism in Aspergillus species. Together, this study identified conserved components required for GT protection in Aspergillus species.

A. fumigatus secretes mycotoxins that are essential for its virulence and pathogenicity. Gliotoxin (GT) is a sulfur-containing mycotoxin, which is known to impair several aspects of the human immune response. GT is also toxic to different fungal species, which have evolved several GT protection strategies. To further decipher these responses, we used transcriptional profiling aiming to compare the response to GT in the GT producer A. fumigatus and the GT non-producer A. nidulans. This analysis allowed us to identify additional genes with a potential role in GT protection. We also identified 15 transcription factors (TFs) encoded in the A. fumigatus genome that are important for conferring resistance to exogenous gliotoxin. One of these TFs, KojR, which is essential for A. oryzae kojic acid production, is also important for virulence in A. fumigatus and GT protection in A. fumigatus, A. nidulans and A. oryzae. KojR regulates the expression of genes important for gliotoxin biosynthesis and protection, and sulfur metabolism. Together, this work identified conserved components required for gliotoxin protection in Aspergillus species.

Chromolaena laevigata (Asteraceae) as a source of endophytic non-aflatoxigenic Aspergillus flavus: chemical profile in different culture conditions and biological applications

Endophytes are microorganisms that form symbiotic relationships with their host. These microorganisms can produce a variety of secondary metabolites, some of which have inhibitory effects on pests and pathogens or even act to promote plant growth. Due to these characteristics, these microorganisms are used as sources of biologically active substances for a wide range of biotechnological applications. Based on that, the aim of this study was to evaluate the production of metabolites of the endophytic Aspergillus flavus CL7 isolated from Chromolaena laevigata, in four different cultivation conditions, and to determine the antimicrobial, cytotoxic, antiviral, and antioxidant potential of these extracts. The multiphasic approach used to identify this strain was based on morphology and ITS gene sequence analysis. The chemical investigation of A. flavus using potato dextrose and minimal medium, using both stationary and agitated methods, resulted in the isolation of kojic acid, α-cyclopiazonic acid, and 20,25-dihydroxyaflavinine. Another 18 compounds in these extracts were identified by UHPLC-HRMS/MS, of which dideacetyl parasiticolide A has been described for the first time from A. flavus. Aflatoxins, important chemomarkers of A. flavus, were not detected in any of the extracts, thus indicating that the CL7 strain is non-aflatoxigenic. The biological potential of all extracts was evaluated, and the best results were observed for the extract obtained using minimal medium against Trichophyton rubrum and Mycobacterium tuberculosis.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes

•

We detected three subgroups of L. edodes with robust phenotypic differentiation.

•

The three subgroups were diverged 36,871 generations ago.

•

L. edodes cultivars in China might originate from the vicinity of Northeast China.

•

We dissected the genetic basis of adaptive evolution in L. edodes.

•

Genes related to fruiting body development are involved in adaptive evolution.

Introduction

Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi.

Objectives

This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution.

Methods

We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution.

Results

We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation.

Conclusion

Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Koji Starter and Koji World in Japan

Koji is made by culturing koji mold on grains. Koji has wide-ranging applications, for example, in alcoholic beverages and seasonings. The word 'mold' generally has a bad image, but in Japan, koji mold is valued for its usefulness, and over the years, efforts have been made to make safe, stable, and delicious food products from it. Koji mold spores, essential when making koji, are called koji starter in the industry. From the many available strains, those suitable for the production of each fermented food are chosen based on indicators such as growth rate and enzyme production capacity. In manufacturing using microorganisms, purity and yield are prioritized. However, the production of fermented foods using koji is more complex, with focus not only on the degree of decomposition of raw materials but also on factors influencing overall product design, including palatability, color, smell, and texture. Production can be facilitated by the variety of koji brought about by the diversity of koji mold combined with the solid culture method which increases the amount of enzyme production. In this report, we introduce the history of koji starter in Japan, the characteristics of koji mold in practice, and various fermented foods made from it. In addition, the factors affecting the quality of koji in solid culture are described.

Comparison of Two Aspergillus oryzae Genomes From Different Clades Reveals Independent Evolution of Alpha-Amylase Duplication, Variation in Secondary Metabolism Genes, and Differences in Primary Metabolism

Microbes (bacteria, yeasts, molds), in addition to plants and animals, were domesticated for their roles in food preservation, nutrition and flavor. Aspergillus oryzae is a domesticated filamentous fungal species traditionally used during fermentation of Asian foods and beverage, such as sake, soy sauce, and miso. To date, little is known about the extent of genome and phenotypic variation of A. oryzae isolates from different clades. Here, we used long-read Oxford Nanopore and short-read Illumina sequencing to produce a highly accurate and contiguous genome assemble of A. oryzae 14160, an industrial strain from China. To understand the relationship of this isolate, we performed phylogenetic analysis with 90 A. oryzae isolates and 1 isolate of the A. oryzae progenitor, Aspergillus flavus. This analysis showed that A. oryzae 14160 is a member of clade A, in comparison to the RIB 40 type strain, which is a member of clade F. To explore genome variation between isolates from distinct A. oryzae clades, we compared the A. oryzae 14160 genome with the complete RIB 40 genome. Our results provide evidence of independent evolution of the alpha-amylase gene duplication, which is one of the major adaptive mutations resulting from domestication. Synteny analysis revealed that both genomes have three copies of the alpha-amylase gene, but only one copy on chromosome 2 was conserved. While the RIB 40 genome had additional copies of the alpha-amylase gene on chromosomes III, and V, 14160 had a second copy on chromosome II and an third copy on chromosome VI. Additionally, we identified hundreds of lineage specific genes, and putative high impact mutations in genes involved in secondary metabolism, including several of the core biosynthetic genes. Finally, to examine the functional effects of genome variation between strains, we measured amylase activity, proteolytic activity, and growth rate on several different substrates. RIB 40 produced significantly higher levels of amylase compared to 14160 when grown on rice and starch. Accordingly, RIB 40 grew faster on rice, while 14160 grew faster on soy. Taken together, our analyses reveal substantial genome and phenotypic variation within A. oryzae.

Crucial role of the intracellular α-glucosidase MalT in the activation of the transcription factor AmyR essential for amylolytic gene expression in Aspergillus oryzae

We examined the role of the intracellular α-glucosidase gene malT, which is part of the maltose-utilizing cluster (MAL cluster) together with malR and malP, in amylolytic gene expression in Aspergillus oryzae. malT disruption severely affected fungal growth on medium containing maltose or starch. Furthermore, the transcription level of the α-amylase gene was significantly reduced by malT disruption. Given that the transcription factor AmyR responsible for amylolytic gene expression is activated by isomaltose converted from maltose incorporated into the cells, MalT may have transglycosylation activity that converts maltose to isomaltose. Indeed, transglycosylated products such as isomaltose/maltotriose and panose were generated from the substrate maltose by MalT purified from a malT-overexpressing strain. The results of this study, taken together, suggests that MalT plays a pivotal role in AmyR activation via its transglycosylation activity that converts maltose to the physiological inducer isomaltose.

Analysis of genomic characteristics and their influence on metabolism in Aspergillus luchuensis albino mutants using genome sequencing

Black Aspergillus luchuensis and its white albino mutant are essential fungi for making alcoholic beverages in Japan. A large number of industrial strains have been created using novel isolation or gene/genome mutation techniques. Such mutations influence metabolic and phenotypic characteristics in industrial strains, but few comparative studies of inter-strain mutation have been conducted. We carried out comparative genome analyses of 8 industrial strains of A. luchuensis and A. kawachii IFO 4308, the latter being the first albino strain to be isolated. Phylogenetic analysis based on 8938 concatenated genes exposed the diversity of black koji strains and uniformity among albino industrial strains, suggesting that passaged industrial albino strains have more genetic mutations compared with strain IFO 4308 and black koji strains. Comparative analysis showed that the albino strains had mutations in genes not only for conidial pigmentation but also in those that encode N-terminal acetyltransferase A and annexin XIV-like protein. The results also suggest that some mutations may have emerged through subculturing of albino strains. For example, mutations in the genes for isocitrate lyase and sugar transporters were observed only in industrial albino strains. This implies that selective pressure for increasing enzyme activity or secondary metabolites may have influenced the mutation of genes associated with environmental stress responses in A. luchuensis albino strains. Our study clarifies hitherto unknown genetic and metabolic characteristics of A. luchuensis industrial strains and provides potential applications for comparative genome analysis for breeding koji strains.

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.

Macrosynteny analysis between Lentinula edodes and Lentinula novae-zelandiae reveals signals of domestication in Lentinula edodes

The basidiomycete fungus Lentinula novae-zelandiae is endemic to New Zealand and is a sister taxon to Lentinula edodes, the second most cultivated mushroom in the world. To explore the biology of this organism, a high-quality chromosome level reference genome of L. novae-zelandiae was produced. Macrosyntenic comparisons between the genome assembly of L. novae-zelandiae, L. edodes and a set of three genome assemblies of diverse species from the Agaricomycota reveal a high degree of macrosyntenic restructuring within L. edodes consistent with signal of domestication. These results show L. edodes has undergone significant genomic change during the course of its evolutionary history, likely a result of its cultivation and domestication over the last 1000 years.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Strong effect of Penicillium roqueforti populations on volatile and metabolic compounds responsible for aromas, flavor and texture in blue cheeses

Studies of food microorganism domestication can provide important insight into adaptation mechanisms and lead to commercial applications. Penicillium roqueforti is a fungus with four genetically differentiated populations, two of which were independently domesticated for blue cheese-making, with the other two populations thriving in other environments. Most blue cheeses are made with strains from a single P. roqueforti population, whereas Roquefort cheeses are inoculated with strains from a second population. We made blue cheeses in accordance with the production specifications for Roquefort-type cheeses, inoculating each cheese with a single P. roqueforti strain, using a total of three strains from each of the four populations. We investigated differences between the cheeses made with the strains from the four P. roqueforti populations, in terms of the induced flora, the proportion of blue color, water activity and the identity and abundance of aqueous and organic metabolites as proxies for proteolysis and lipolysis as well as volatile compounds responsible for flavor and aroma. We found that the population-of-origin of the P. roqueforti strains used for inoculation had a minor impact on bacterial diversity and no effect on the abundance of the main microorganism. The cheeses produced with P. roqueforti strains from cheese populations had a higher percentage of blue area and a higher abundance of the volatile compounds typical of blue cheeses, such as methyl ketones and secondary alcohols. In particular, the Roquefort strains produced higher amounts of these aromatic compounds, partly due to more efficient proteolysis and lipolysis. The Roquefort strains also led to cheeses with a lower water availability, an important feature for preventing spoilage in blue cheeses, which is subject to controls for the sale of Roquefort cheese. The typical appearance and flavors of blue cheeses thus result from human selection on P. roqueforti, leading to the acquisition of specific features by the two cheese populations. These findings have important implications for our understanding of adaptation and domestication, and for cheese improvement.

Rapid selection response to ethanol in Saccharomyces eubayanus emulates the domestication process under brewing conditions

Although the typical genomic and phenotypic changes that characterize the evolution of organisms under the human domestication syndrome represent textbook examples of rapid evolution, the molecular processes that underpin such changes are still poorly understood. Domesticated yeasts for brewing, where short generation times and large phenotypic and genomic plasticity were attained in a few generations under selection, are prime examples. To experimentally emulate the lager yeast domestication process, we created a genetically complex (panmictic) artificial population of multiple Saccharomyces eubayanus genotypes, one of the parents of lager yeast. Then, we imposed a constant selection regime under a high ethanol concentration in 10 replicated populations during 260 generations (6 months) and compared them with propagated controls exposed solely to glucose. Propagated populations exhibited a selection differential of 60% in growth rate in ethanol, mostly explained by the proliferation of a single lineage (CL248.1) that competitively displaced all other clones. Interestingly, the outcome does not require the entire time‐course of adaptation, as four lineages monopolized the culture at generation 120. Sequencing demonstrated that de novo genetic variants were produced in all propagated lines, including SNPs, aneuploidies, INDELs and translocations. In addition, the different propagated populations showed correlated responses resembling the domestication syndrome: genomic rearrangements, faster fermentation rates, lower production of phenolic off‐flavours and lower volatile compound complexity. Expression profiling in beer wort revealed altered expression levels of genes related to methionine metabolism, flocculation, stress tolerance and diauxic shift, likely contributing to higher ethanol and fermentation stress tolerance in the evolved populations. Our study shows that experimental evolution can rebuild the brewing domestication process in ‘fast motion’ in wild yeast, and also provides a powerful tool for studying the genetics of the adaptation process in complex populations.

The International Scientific Association for Probiotics and Prebiotics (ISAPP) consensus statement on fermented foods

An expert panel was convened in September 2019 by The International Scientific Association for Probiotics and Prebiotics (ISAPP) to develop a definition for fermented foods and to describe their role in the human diet. Although these foods have been consumed for thousands of years, they are receiving increased attention among biologists, nutritionists, technologists, clinicians and consumers. Despite this interest, inconsistencies related to the use of the term 'fermented' led the panel to define fermented foods and beverages as "foods made through desired microbial growth and enzymatic conversions of food components". This definition, encompassing the many varieties of fermented foods, is intended to clarify what is (and is not) a fermented food. The distinction between fermented foods and probiotics is further clarified. The panel also addressed the current state of knowledge on the safety, risks and health benefits, including an assessment of the nutritional attributes and a mechanistic rationale for how fermented foods could improve gastrointestinal and general health. The latest advancements in our understanding of the microbial ecology and systems biology of these foods were discussed. Finally, the panel reviewed how fermented foods are regulated and discussed efforts to include them as a separate category in national dietary guidelines.

Cleaning the Cellular Factory-Deletion of McrA in Aspergillus oryzae NSAR1 and the Generation of a Novel Kojic Acid Deficient Strain for Cleaner Heterologous Production of Secondary Metabolites

The use of filamentous fungi as cellular factories, where natural product pathways can be refactored and expressed in a host strain, continues to aid the field of natural product discovery. Much work has been done to develop host strains which are genetically tractable, and for which there are multiple selectable markers and controllable expression systems. To fully exploit these strains, it is beneficial to understand their natural metabolic capabilities, as such knowledge can rule out host metabolites from analysis of transgenic lines and highlight any potential interplay between endogenous and exogenous pathways. Additionally, once identified, the deletion of secondary metabolite pathways from host strains can simplify the detection and purification of heterologous compounds. To this end, secondary metabolite production in Aspergillus oryzae strain NSAR1 has been investigated via the deletion of the newly discovered negative regulator of secondary metabolism, mcrA (multicluster regulator A). In all ascomycetes previously studied mcrA deletion led to an increase in secondary metabolite production. Surprisingly, the only detectable phenotypic change in NSAR1 was a doubling in the yields of kojic acid, with no novel secondary metabolites produced. This supports the previous claim that secondary metabolite production has been repressed in A. oryzae and demonstrates that such repression is not McrA-mediated. Strain NSAR1 was then modified by employing CRISPR-Cas9 technology to disrupt the production of kojic acid, generating the novel strain NSARΔK, which combines the various beneficial traits of NSAR1 with a uniquely clean secondary metabolite background.