Multiple recent horizontal transfers of a large genomic region in cheese making fungi

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.

Identification of quantitative trait loci (QTLs) for key cheese making phenotypes in the blue-cheese mold Penicillium roqueforti

Elucidating the genomic architecture of quantitative traits is essential for our understanding of adaptation and for breeding in domesticated organisms. Penicillium roqueforti is the mold used worldwide for the blue cheese maturation, contributing to flavors through proteolytic and lipolytic activities. The two domesticated cheese populations display very little genetic diversity, but are differentiated and carry opposite mating types. We produced haploid F1 progenies from five crosses, using parents belonging to cheese and non-cheese populations. Analyses of high-quality genome assemblies of the parental strains revealed five large translocations, two having occurred via a circular intermediate, one with footprints of Starship giant mobile elements. Offspring genotyping with genotype-by-sequencing (GBS) revealed several genomic regions with segregation distortion, possibly linked to degeneration in cheese lineages. We found transgressions for several traits relevant for cheese making, with offspring having more extreme trait values than parental strains. We identified quantitative trait loci (QTLs) for colony color, lipolysis, proteolysis, extrolite production, including mycotoxins, but not for growth rates. Some genomic regions appeared rich in QTLs for both lipid and protein metabolism, and other regions for the production of multiple extrolites, indicating that QTLs have pleiotropic effects. Some QTLs corresponded to known biosynthetic gene clusters, e.g., for the production of melanin or extrolites. F1 hybrids constitute valuable strains for cheese producers, with new traits and new allelic combinations, and allowed identifying target genomic regions for traits important in cheese making, paving the way for strain improvement. The findings further contribute to our understanding of the genetic mechanisms underlying rapid adaptation, revealing convergent adaptation targeting major gene regulators.

Comparative genomics reveals substantial divergence in metal sensitive and metal tolerant isolates of the ericoid mycorrhizal fungus Oidiodendron maius

Some heavy metal tolerant fungal isolates capable of forming ericoid mycorrhiza can also confer increased metal tolerance to the host plant. One of these fungal isolates, Oidiodendron maius Zn, has been characterized and a few molecular mechanisms underlying its metal tolerant phenotype have been identified. Here, we investigate the genomic divergences between the available genome of O. maius Zn and the genomes of metal tolerant and sensitive isolates of O. maius, with the aim of identifying genes or intergenic regions possibly involved in the display of the tolerance. The resequenced genomes of 8 tolerant and 10 sensitive isolates were mapped on the reference, O. maius Zn, yielding 357 gene models from the reference that were either missing or too polymorphic to be identified in the genomes of the sensitive isolates. These regions included genes with functions related to defense mechanisms and with unknown functions. One third of the predicted gene models turned out to be highly polymorphic, including many enriched GO terms, i.e. DNA/RNA metabolism and modification, chromosome/chromatin organization, protein biosynthesis, metabolism and function, energy consumption/transfer and mitochondrion. Overall, our findings indicate that the tolerant phenotype in O. maius likely arises from multiple genetic adaptations rather than a singular mechanism.

An Inversion Polymorphism Under Balancing Selection, Involving Giant Mobile Elements, in an Invasive Fungal Pathogen

Recombination suppression can evolve in sex or mating-type chromosomes, or in autosomal supergenes, with different haplotypes being maintained by balancing selection. In the invasive chestnut blight fungus Cryphonectria parasitica, a genomic region was suggested to lack recombination and to be partially physically linked to the mating-type (MAT) locus based on segregation analyses. Using hundreds of available C. parasitica genomes and generating new high-quality genome assemblies, we show that a ca. 1.2 Mb genomic region proximal to the mating-type locus lacks recombination, with the segregation of two highly differentiated haplotypes in balanced proportions in invasive populations. High-quality genome assemblies further revealed an inversion in one of the haplotypes in the invaded range. The two haplotypes were estimated to have diverged 1.5 million years ago, and each harboured specific genes, some of which likely belonging to Starships. These are large transposable elements, mobilized by tyrosine recombinases, able to move accessory genes, and involved in adaptation in multiple fungi. The MAT-proximal region carried genes upregulated under virus infection or vegetative incompatibility reaction. In the native range, the MAT-proximal region also appeared to have a different evolutionary history than the rest of the genome. In all continents, the MAT-Proximal region was enriched in nonsynonymous substitutions, in gene presence/absence polymorphism, in tyrosine recombinases and in transposable elements. This study thus sheds light on a case of a large nonrecombining region partially linked to a mating compatibility locus, likely maintained by balancing selection on differentiated haplotypes, possibly involved in adaptation in a devastating tree pathogen.

Mycotoxins in Cheese: Assessing Risks, Fungal Contaminants, and Control Strategies for Food Safety

According to the scientific information reviewed, cheese is highly susceptible to contamination by mycotoxin-producing fungi, primarily species from the genera Aspergillus (A. niger, A. flavus) and Penicillium (P. commune, P. solitum, P. palitans, and P. crustosum). Studies on various types of cheese made from cow's milk report an average concentration of Aflatoxin M1 (AFM1) at 13,000 ng kg-1, which is alarming since the regulatory limits for AFM1 in cheese range from 250 to 500 ng kg-1. For instance, limits set by Codex Alimentarius, the European Commission (EC), Turkey, and Iran are 250 ng kg-1. In the Netherlands, the limit is 200 ng kg-1, and in Italy, it is 450 ng kg-1. However, the concentration of mycotoxins frequently exceeds these regulatory limits, including critical mycotoxins such as ochratoxin A, citrinin, and cyclopiazonic acid, which pose significant global health concerns. Therefore, this study aims to review the mycobiota responsible for producing key mycotoxins in cheese and to assess the influence of physicochemical factors on fungal growth and mycotoxin production. By incorporating control strategies such as hygiene practices, pasteurization, and the use of preservatives, this study seeks to improve methodologies in the cheese production chain and mitigate contamination by fungi and mycotoxins.

Omics-Based Comparison of Fungal Virulence Genes, Biosynthetic Gene Clusters, and Small Molecules in Penicillium expansum and Penicillium chrysogenum

Penicillium expansum is a ubiquitous pathogenic fungus that causes blue mold decay of apple fruit postharvest, and another member of the genus, Penicillium chrysogenum, is a well-studied saprophyte valued for antibiotic and small molecule production. While these two fungi have been investigated individually, a recent discovery revealed that P. chrysogenum can block P. expansum-mediated decay of apple fruit. To shed light on this observation, we conducted a comparative genomic, transcriptomic, and metabolomic study of two P. chrysogenum (404 and 413) and two P. expansum (Pe21 and R19) isolates. Global transcriptional and metabolomic outputs were disparate between the species, nearly identical for P. chrysogenum isolates, and different between P. expansum isolates. Further, the two P. chrysogenum genomes revealed secondary metabolite gene clusters that varied widely from P. expansum. This included the absence of an intact patulin gene cluster in P. chrysogenum, which corroborates the metabolomic data regarding its inability to produce patulin. Additionally, a core subset of P. expansum virulence gene homologues were identified in P. chrysogenum and were similarly transcriptionally regulated in vitro. Molecules with varying biological activities, and phytohormone-like compounds were detected for the first time in P. expansum while antibiotics like penicillin G and other biologically active molecules were discovered in P. chrysogenum culture supernatants. Our findings provide a solid omics-based foundation of small molecule production in these two fungal species with implications in postharvest context and expand the current knowledge of the Penicillium-derived chemical repertoire for broader fundamental and practical applications.

Gene acquisition by giant transposons primes eukaryotes for rapid evolution via horizontal gene transfer

Horizontal gene transfer (HGT) disseminates genetic information between species and is a powerful mechanism of adaptation. Yet, we know little about its underlying drivers in eukaryotes. Giant Starship transposons have been implicated as agents of fungal HGT, providing an unprecedented opportunity to reveal the evolutionary parameters behind this process. Here, we characterize the ssf gene cluster, which contributes to formaldehyde resistance, and use it to demonstrate how mobile element evolution shapes fungal adaptation. We found that ssf clusters have been acquired by various distantly related Starships, which each exhibit multiple instances of horizontal transfer across fungal species (at least nine events, including between different taxonomic orders). Many ssf clusters have subsequently integrated into their host's genome, illustrating how Starships shape the evolutionary trajectory of fungal hosts beyond any single transfer. Our results demonstrate the key role Starships play in mediating rapid and repeated adaptation via HGT, elevating the importance of mobile element evolution in eukaryotic biology.

Genetic diversity and population structure of Penicillium roqueforti isolates from Turkish blue cheeses

Blue cheeses, including renowned mold-ripened varieties such as Roquefort (France), Gorgonzola (Italy), Stilton (UK), Danablue (Denmark), and Cabrales (Spain), owe their distinct blue-green color and unique flavor to the fungal species Penicillium roqueforti. In Turkey, traditional cheeses similar to blue cheeses, namely mold-ripened Tulum and Civil, employ production techniques distinct from their European counterparts. Notably, mold-ripening in Turkish cheeses is spontaneous and does not involve starter cultures. Despite P. roqueforti being recognized for its distinct genetic populations sourced from various blue cheeses and non-cheese origins globally, the characteristics of the P. roqueforti population within Turkish cheeses remain unexplored. This study aimed to unravel the genetic characteristics and population structure of P. roqueforti from Turkish mold-ripened cheeses. Analysis of mold-ripened Civil (n = 22) and Tulum (n = 8) samples revealed 66 P. roqueforti isolates (76.6 % of total fungal isolates). Subsequently, these isolates (n = 66) and those from previous studies (Tulum n = 53, Golot n = 1) were used to assess genetic characteristics and mating genotypes. All 120 isolates harbored horizontal transfer regions (Wallaby and CheesyTer) and predominantly possessed the MAT1-2 mating genotype, similar to global blue cheese populations. However, most lacked the mpaC deletion associated with such populations. Analysis of the population with three polymorphic microsatellite markers revealed 36 haplotypes (HTs). Some cheeses contained isolates with different HTs or opposite mating genotypes, aligning with spontaneous fungal growth. Tulum and Civil isolates exhibited similar population diversity without forming distinct subgroups. Phylogenetic analysis of 20 selected isolates showed 75 % aligning with global blue cheese isolates, while 25 % formed unique clades. Overall, Turkish P. roqueforti isolates share genetic similarities with global populations but exhibit unique characteristics, suggesting potential new clades deserving further investigation. This research illuminates the characteristics of P. roqueforti isolates from Turkish cheeses, contributing to the knowledge of the global intraspecific diversity of the P. roqueforti species.

The Near-Gapless Penicillium fuscoglaucum Genome Enables the Discovery of Lifestyle Features as an Emerging Post-Harvest Phytopathogen

Penicillium spp. occupy many diverse biological niches that include plant pathogens, opportunistic human pathogens, saprophytes, indoor air contaminants, and those selected specifically for industrial applications to produce secondary metabolites and lifesaving antibiotics. Recent phylogenetic studies have established Penicillium fuscoglaucum as a synonym for Penicillium commune, which is an indoor air contaminant and toxin producer and can infect apple fruit during storage. During routine culturing on selective media in the lab, we obtained an isolate of P. fuscoglaucum Pf_T2 and sequenced its genome. The Pf_T2 genome is far superior to available genomic resources for the species. Our assembly exhibits a length of 35.1 Mb, a BUSCO score of 97.9% complete, and consists of five scaffolds/contigs representing the four expected chromosomes. It was determined that the Pf_T2 genome was colinear with a type specimen P. fuscoglaucum and contained a lineage-specific, intact cyclopiazonic acid (CPA) gene cluster. For comparison, a highly virulent postharvest apple pathogen, P. expansum strain TDL 12.1, was included and showed a similar growth pattern in culture to our Pf_T2 isolate but was far more aggressive in apple fruit than P. fuscoglaucum. The genome of Pf_T2 serves as a major improvement over existing resources, has superior annotation, and can inform forthcoming omics-based work and functional genetic studies to probe secondary metabolite production and disparities in aggressiveness during apple fruit decay.

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.

Tracking alternative versions of the galactose gene network in the genus Saccharomyces and their expansion after domestication

When Saccharomyces cerevisiae grows on mixtures of glucose and galactose, galactose utilization is repressed by glucose, and induction of the GAL gene network only occurs when glucose is exhausted. Contrary to reference GAL alleles, alternative alleles support faster growth on galactose, thus enabling distinct galactose utilization strategies maintained by balancing selection. Here, we report on new wild populations of Saccharomyces cerevisiae harboring alternative GAL versions and, for the first time, of Saccharomyces paradoxus alternative alleles. We also show that the non-functional GAL version found earlier in Saccharomyces kudriavzevii is phylogenetically related to the alternative versions, which constitutes a case of trans-specific maintenance of highly divergent alleles. Strains harboring the different GAL network variants show different levels of alleviation of glucose repression and growth proficiency on galactose. We propose that domestication involved specialization toward thriving in milk from a generalist ancestor partially adapted to galactose consumption in the plant niche.

Domestication signatures in the non-conventional yeast Lachancea cidri

Evaluating domestication signatures beyond model organisms is essential for a thorough understanding of the genotype-phenotype relationship in wild and human-related environments. Structural variations (SVs) can significantly impact phenotypes playing an important role in the physiological adaptation of species to different niches, including during domestication. A detailed characterization of the fitness consequences of these genomic rearrangements, however, is still limited in non-model systems, largely due to the paucity of direct comparisons between domesticated and wild isolates. Here, we used a combination of sequencing strategies to explore major genomic rearrangements in a Lachancea cidri yeast strain isolated from cider (CBS2950) and compared them to those in eight wild isolates from primary forests. Genomic analysis revealed dozens of SVs, including a large reciprocal translocation (~16 kb and 500 kb) present in the cider strain, but absent from all wild strains. Interestingly, the number of SVs was higher relative to single-nucleotide polymorphisms in the cider strain, suggesting a significant role in the strain's phenotypic variation. The set of SVs identified directly impacts dozens of genes and likely underpins the greater fermentation performance in the L. cidri CBS2950. In addition, the large reciprocal translocation affects a proline permease (PUT4) regulatory region, resulting in higher PUT4 transcript levels, which agrees with higher ethanol tolerance, improved cell growth when using proline, and higher amino acid consumption during fermentation. These results suggest that SVs are responsible for the rapid physiological adaptation of yeast to a human-related environment and demonstrate the key contribution of SVs in adaptive fermentative traits in non-model species.IMPORTANCEThe exploration of domestication signatures associated with human-related environments has predominantly focused on studies conducted on model organisms, such as Saccharomyces cerevisiae, overlooking the potential for comparisons across other non-Saccharomyces species. In our research, employing a combination of long- and short-read data, we found domestication signatures in Lachancea cidri, a non-model species recently isolated from fermentative environments in cider in France. The significance of our study lies in the identification of large array of major genomic rearrangements in a cider strain compared to wild isolates, which underly several fermentative traits. These domestication signatures result from structural variants, which are likely responsible for the phenotypic differences between strains, providing a rapid path of adaptation to human-related environments.

New colours for old in the blue-cheese fungus Penicillium roqueforti

Penicillium roqueforti is used worldwide in the production of blue-veined cheese. The blue-green colour derives from pigmented spores formed by fungal growth. Using a combination of bioinformatics, targeted gene deletions, and heterologous gene expression we discovered that pigment formation was due to a DHN-melanin biosynthesis pathway. Systematic deletion of pathway genes altered the arising spore colour, yielding white to yellow-green to red-pink-brown phenotypes, demonstrating the potential to generate new coloured strains. There was no consistent impact on mycophenolic acid production as a result of pathway interruption although levels of roquefortine C were altered in some deletants. Importantly, levels of methyl-ketones associated with blue-cheese flavour were not impacted. UV-induced colour mutants, allowed in food production, were then generated. A range of colours were obtained and certain phenotypes were successfully mapped to pathway gene mutations. Selected colour mutants were subsequently used in cheese production and generated expected new colourations with no elevated mycotoxins, offering the exciting prospect of use in future cheese manufacture.

That's no moon, it's a Starship: Giant transposons driving fungal horizontal gene transfer

To date, most reports of horizontal gene transfer (HGT) in fungi rely on genome sequence data and are therefore an indirect measure of HGT after the event has occurred. However, a novel group of class II-like transposons known as Starships may soon alter this status quo. Starships are giant transposable elements that carry dozens of genes, some of which are host-beneficial, and are linked to many recent HGT events in the fungal kingdom. These transposons remain active and mobile in many fungal genomes and their transposition has recently been shown to be driven by a conserved tyrosine-recombinase called 'Captain'. This perspective explores some of the remaining unanswered questions about how these Starship transposons move, both within a genome and between different species. We seek to outline several experimental approaches that can be used to identify the genes essential for Starship-mediated HGT and draw links to other recently discovered giant transposons outside of the fungal kingdom.

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.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

Secondary Metabolites Produced by the Blue-Cheese Ripening Mold Penicillium roqueforti; Biosynthesis and Regulation Mechanisms

Filamentous fungi are an important source of natural products. The mold Penicillium roqueforti, which is well-known for being responsible for the characteristic texture, blue-green spots, and aroma of the so-called blue-veined cheeses (French Bleu, Roquefort, Gorgonzola, Stilton, Cabrales, and Valdeón, among others), is able to synthesize different secondary metabolites, including andrastins and mycophenolic acid, as well as several mycotoxins, such as Roquefortines C and D, PR-toxin and eremofortins, Isofumigaclavines A and B, festuclavine, and Annullatins D and F. This review provides a detailed description of the biosynthetic gene clusters and pathways of the main secondary metabolites produced by P. roqueforti, as well as an overview of the regulatory mechanisms controlling secondary metabolism in this filamentous fungus.

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.

The prevalence of killer yeasts and double-stranded RNAs in the budding yeast Saccharomyces cerevisiae

Killer toxins are antifungal proteins produced by many species of "killer" yeasts, including the brewer's and baker's yeast Saccharomyces cerevisiae. Screening 1270 strains of S. cerevisiae for killer toxin production found that 50% are killer yeasts, with a higher prevalence of yeasts isolated from human clinical samples and winemaking processes. Since many killer toxins are encoded by satellite double-stranded RNAs (dsRNAs) associated with mycoviruses, S. cerevisiae strains were also assayed for the presence of dsRNAs. This screen identified that 51% of strains contained dsRNAs from the mycovirus families Totiviridae and Partitiviridae, as well as satellite dsRNAs. Killer toxin production was correlated with the presence of satellite dsRNAs but not mycoviruses. However, in most killer yeasts, whole genome analysis identified the killer toxin gene KHS1 as significantly associated with killer toxin production. Most killer yeasts had unique spectrums of antifungal activities compared to canonical killer toxins, and sequence analysis identified mutations that altered their antifungal activities. The prevalence of mycoviruses and killer toxins in S. cerevisiae is important because of their known impact on yeast fitness, with implications for academic research and industrial application of this yeast species.

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.

The Epichloë festucae Antifungal Protein Efe-AfpA Protects Creeping Bentgrass (Agrostis stolonifera) from the Plant Pathogen Clarireedia jacksonii, the Causal Agent of Dollar Spot Disease

Dollar spot disease, caused by the fungal pathogen Clarireedia jacksonii, is a major problem in many turfgrass species, particularly creeping bentgrass (Agrostis stolonifera). It is well-established that strong creeping red fescue (Festuca rubra subsp. rubra) exhibits good dollar spot resistance when infected by the fungal endophyte Epichloë festucae. This endophyte-mediated disease resistance is unique to the fine fescues and has not been observed in other grass species infected with other Epichloë spp. The mechanism underlying the unique endophyte-mediated disease resistance in strong creeping red fescue has not yet been established. We pursued the possibility that it may be due to the presence of an abundant secreted antifungal protein produced by E. festucae. Here, we compare the activity of the antifungal protein expressed in Escherichia coli, Pichia pastoris, and Penicillium chrysogenum. Active protein was recovered from all systems, with the best activity being from Pe. chrysogenum. In greenhouse assays, topical application of the purified antifungal protein to creeping bentgrass and endophyte-free strong creeping red fescue protected the plants from developing severe symptoms caused by C. jacksonii. These results support the hypothesis that Efe-AfpA is a major contributor to the dollar spot resistance observed with E. festucae-infected strong creeping red fescue in the field, and that this protein could be developed as an alternative or complement to fungicides for the management of this disease on turfgrasses.

Biosynthesis of Annullatin D in Penicillium roqueforti Implies Oxidative Lactonization between Two Hydroxyl Groups Catalyzed by a BBE-like Enzyme

Annullatins from Cordyceps annullata are alkylated aromatic polyketides including annullatin D with a fused dihydrobenzofuran lactone ring system. Here, we report the identification of a silent biosynthetic gene cluster for annullatins from Penicillium roqueforti by heterologous expression in Aspergillus nidulans, gene deletion, and feeding experiments as well as by biochemical characterization. The polyketide core structure is consecutively modified by hydroxylation and prenylation. A berberine bridge enzyme-like protein catalyzes the final step, an oxidative lactonization between two hydroxyl groups, to form (2S, 9S)-annullatin D.

High sorbic acid resistance of Penicillium roqueforti is mediated by the SORBUS gene cluster

Penicillium roqueforti is a major food-spoilage fungus known for its high resistance to the food preservative sorbic acid. Here, we demonstrate that the minimum inhibitory concentration of undissociated sorbic acid (MIC_u) ranges between 4.2 and 21.2 mM when 34 P. roqueforti strains were grown on malt extract broth. A genome-wide association study revealed that the six most resistant strains contained the 180 kbp gene cluster SORBUS, which was absent in the other 28 strains. In addition, a SNP analysis revealed five genes outside the SORBUS cluster that may be linked to sorbic acid resistance. A partial SORBUS knock-out (>100 of 180 kbp) in a resistant strain reduced sorbic acid resistance to similar levels as observed in the sensitive strains. Whole genome transcriptome analysis revealed a small set of genes present in both resistant and sensitive P. roqueforti strains that were differentially expressed in the presence of the weak acid. These genes could explain why P. roqueforti is more resistant to sorbic acid when compared to other fungi, even in the absence of the SORBUS cluster. Together, the MIC_u of 21.2 mM makes P. roqueforti among the most sorbic acid-resistant fungi, if not the most resistant fungus, which is mediated by the SORBUS gene cluster.

Chemical preservatives, such as sorbic acid, are often used in food to prevent spoilage by fungi, yet some fungi are particularly well-suited to deal with these preservatives. First, we investigated the resistance of 34 Penicillium roqueforti strains to various food preservatives. This revealed that some strains were highly resistant to sorbic acid, while others are more sensitive. Next, we used DNA sequencing to compare the genetic variation between these strains and discovered a specific genetic region (SORBUS) that is unique to the resistant strains. Through comparative analysis with other fungal species the SORBUS region was studied in more detail and with the use of genetic engineering tools we removed this unique region. Finally, the mutant lacking the SORBUS region was confirmed to have lost its sorbic acid resistance. This finding is of particular interest as it suggests that only some, not all, P. roqueforti strains are potent spoilers and that specific genetic markers could help in the identification of resistant strains.

Recent developments in horizontal gene transfer with the adaptive innovation of fermented foods

Horizontal gene transfer (HGT) has contributed significantly to the adaptability of bacteria, yeast and mold in fermented foods, whose evidence has been found in several fermented foods. Although not every HGT has biological significance, it plays an important role in improving the quality of fermented foods. In this review, how HGT facilitated microbial domestication and adaptive evolution in fermented foods was discussed. HGT can assist in the industrial innovation of fermented foods, and this adaptive evolution strategy can improve the quality of fermented foods. Additionally, the mechanism underlying HGT in fermented foods were analyzed. Furthermore, the critical bottlenecks involved in optimizing HGT during the production of fermented foods and strategies for optimizing HGT were proposed. Finally, the prospect of HGT for promoting the industrial innovation of fermented foods was highlighted. The comprehensive report on HGT in fermented foods provides a new trend for domesticating preferable starters for food fermentation, thus optimizing the quality and improving the industrial production of fermented foods.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

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.

A large transposable element mediates metal resistance in the fungus Paecilomyces variotii

The horizontal transfer of large gene clusters by mobile elements is a key driver of prokaryotic adaptation in response to environmental stresses. Eukaryotic microbes face similar stresses; however, a parallel role for mobile elements has not been established. A stress faced by many microorganisms is toxic metal ions in their environment. In fungi, identified mechanisms for protection against metals generally rely on genes that are dispersed within an organism's genome. Here, we discover a large (∼85 kb) region that confers tolerance to five metal/metalloid ions (arsenate, cadmium, copper, lead, and zinc) in the genomes of some, but not all, strains of a fungus, Paecilomyces variotii. We name this region HEPHAESTUS (Hφ) and present evidence that it is mobile within the P. variotii genome with features characteristic of a transposable element. HEPHAESTUS contains the greatest complement of host-beneficial genes carried by a transposable element in eukaryotes, suggesting that eukaryotic transposable elements might play a role analogous to bacteria in the horizontal transfer of large regions of host-beneficial DNA. Genes within HEPHAESTUS responsible for individual metal tolerances include those encoding a P-type ATPase transporter-PcaA-required for cadmium and lead tolerance, a transporter-ZrcA-providing tolerance to zinc, and a multicopper oxidase-McoA-conferring tolerance to copper. In addition, a subregion of Hφ confers tolerance to arsenate. The genome sequences of other fungi in the Eurotiales contain further examples of HEPHAESTUS, suggesting that it is responsible for independently assembling tolerance to a diverse array of ions, including chromium, mercury, and sodium.

Omics Approaches to Assess Flavor Development in Cheese

Cheese is characterized by a rich and complex microbiota that plays a vital role during both production and ripening, contributing significantly to the safety, quality, and sensory characteristics of the final product. In this context, it is vital to explore the microbiota composition and understand its dynamics and evolution during cheese manufacturing and ripening. Application of high-throughput DNA sequencing technologies have facilitated the more accurate identification of the cheese microbiome, detailed study of its potential functionality, and its contribution to the development of specific organoleptic properties. These technologies include amplicon sequencing, whole-metagenome shotgun sequencing, metatranscriptomics, and, most recently, metabolomics. In recent years, however, the application of multiple meta-omics approaches along with data integration analysis, which was enabled by advanced computational and bioinformatics tools, paved the way to better comprehension of the cheese ripening process, revealing significant associations between the cheese microbiota and metabolites, as well as their impact on cheese flavor and quality.

Hallstatt miners consumed blue cheese and beer during the Iron Age and retained a non-Westernized gut microbiome until the Baroque period

We subjected human paleofeces dating from the Bronze Age to the Baroque period (18th century AD) to in-depth microscopic, metagenomic, and proteomic analyses. The paleofeces were preserved in the underground salt mines of the UNESCO World Heritage site of Hallstatt in Austria. This allowed us to reconstruct the diet of the former population and gain insights into their ancient gut microbiome composition. Our dietary survey identified bran and glumes of different cereals as some of the most prevalent plant fragments. This highly fibrous, carbohydrate-rich diet was supplemented with proteins from broad beans and occasionally with fruits, nuts, or animal food products. Due to these traditional dietary habits, all ancient miners up to the Baroque period have gut microbiome structures akin to modern non-Westernized individuals whose diets are also mainly composed of unprocessed foods and fresh fruits and vegetables. This may indicate a shift in the gut community composition of modern Westernized populations due to quite recent dietary and lifestyle changes. When we extended our microbial survey to fungi present in the paleofeces, in one of the Iron Age samples, we observed a high abundance of Penicillium roqueforti and Saccharomyces cerevisiae DNA. Genome-wide analysis indicates that both fungi were involved in food fermentation and provides the first molecular evidence for blue cheese and beer consumption in Iron Age Europe.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

Genomic Analyses of Penicillium Species Have Revealed Patulin and Citrinin Gene Clusters and Novel Loci Involved in Oxylipin Production

Blue mold of apple is caused by several different Penicillium species, among which P. expansum and P. solitum are the most frequently isolated. P. expansum is the most aggressive species, and P. solitum is very weak when infecting apple fruit during storage. In this study, we report complete genomic analyses of three different Penicillium species: P. expansum R21 and P. crustosum NJ1, isolated from stored apple fruit; and P. maximae 113, isolated in 2013 from a flooded home in New Jersey, USA, in the aftermath of Hurricane Sandy. Patulin and citrinin gene cluster analyses explained the lack of patulin production in NJ1 compared to R21 and lack of citrinin production in all three strains. A Drosophila bioassay demonstrated that volatiles emitted by P. solitum SA and P. polonicum RS1 were more toxic than those from P. expansum and P. crustosum strains (R27, R11, R21, G10, and R19). The toxicity was hypothesized to be related to production of eight-carbon oxylipins. Putative lipoxygenase genes were identified in P. expansum and P. maximae strains, but not in P. crustosum. Our data will provide a better understanding of Penicillium spp. complex secondary metabolic capabilities, especially concerning the genetic bases of mycotoxins and toxic VOCs.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

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.

IMA genome - F14 : Draft genome sequences of Penicillium roqueforti, Fusarium sororula, Chrysoporthe puriensis, and Chalaropsis populi

Draft genomes of Penicillium roqueforti, Fusarium sororula, Chalaropsis populi, and Chrysoporthe puriensis are presented. Penicillium roqueforti is a model fungus for genetics, physiological and metabolic studies, as well as for biotechnological applications. Fusarium sororula and Chrysoporthe puriensis are important tree pathogens, and Chalaropsis populi is a soil-borne root-pathogen. The genome sequences presented here thus contribute towards a better understanding of both the pathogenicity and biotechnological potential of these species.

Fungal natural products galaxy: Biochemistry and molecular genetics toward blockbuster drugs discovery

Secondary metabolites synthesized by fungi have become a precious source of inspiration for the design of novel drugs. Indeed, fungi are prolific producers of fascinating, diverse, structurally complex, and low-molecular-mass natural products with high therapeutic leads, such as novel antimicrobial compounds, anticancer compounds, immunosuppressive agents, among others. Given that these microorganisms possess the extraordinary capacity to secrete diverse chemical scaffolds, they have been highly exploited by the giant pharma companies to generate small molecules. This has been made possible because the isolation of metabolites from fungal natural sources is feasible and surpasses the organic synthesis of compounds, which otherwise remains a significant bottleneck in the drug discovery process. Here in this comprehensive review, we have discussed recent studies on different fungi (pathogenic, non-pathogenic, commensal, and endophytic/symbiotic) from different habitats (terrestrial and marines), the specialized metabolites they biosynthesize, and the drugs derived from these specialized metabolites. Moreover, we have unveiled the logic behind the biosynthesis of vital chemical scaffolds, such as NRPS, PKS, PKS-NRPS hybrid, RiPPS, terpenoids, indole alkaloids, and their genetic mechanisms. Besides, we have provided a glimpse of the concept behind mycotoxins, virulence factor, and host immune response based on fungal infections.

The Species-Specific Acquisition and Diversification of a K1-like Family of Killer Toxins in Budding Yeasts of the Saccharomycotina

Killer toxins are extracellular antifungal proteins that are produced by a wide variety of fungi, including Saccharomyces yeasts. Although many Saccharomyces killer toxins have been previously identified, their evolutionary origins remain uncertain given that many of these genes have been mobilized by double-stranded RNA (dsRNA) viruses. A survey of yeasts from the Saccharomyces genus has identified a novel killer toxin with a unique spectrum of activity produced by Saccharomyces paradoxus. The expression of this killer toxin is associated with the presence of a dsRNA totivirus and a satellite dsRNA. Genetic sequencing of the satellite dsRNA confirmed that it encodes a killer toxin with homology to the canonical ionophoric K1 toxin from Saccharomyces cerevisiae and has been named K1-like (K1L). Genomic homologs of K1L were identified in six non-Saccharomyces yeast species of the Saccharomycotina subphylum, predominantly in subtelomeric regions of the genome. When ectopically expressed in S. cerevisiae from cloned cDNAs, both K1L and its homologs can inhibit the growth of competing yeast species, confirming the discovery of a family of biologically active K1-like killer toxins. The sporadic distribution of these genes supports their acquisition by horizontal gene transfer followed by diversification. The phylogenetic relationship between K1L and its genomic homologs suggests a common ancestry and gene flow via dsRNAs and DNAs across taxonomic divisions. This appears to enable the acquisition of a diverse arsenal of killer toxins by different yeast species for potential use in niche competition.

[The domestication of Penicillium cheese fungi]

Domestication is the process of organism evolution under selection by humans, and as such has been a model for studying adaptation since Charles Darwin. Here we review recent studies on the genomics of adaptation and domestication syndrome in two cheese-making fungal lineages, Penicillium roqueforti used for maturing blue cheeses, and the Penicillium camemberti species complex used for making soft cheeses such as Camembert and Brie. Comparative genomics have revealed horizontal gene transfers involved in convergent adaptation to cheese. Population genomics have identified differentiated populations with contrasted traits, several populations having independently been domesticated for cheese making in both P. roqueforti and the Penicillium camemberti species complex, and having undergone bottlenecks. The different cheese populations have acquired traits beneficial for cheese making in comparison to non-cheese populations, regarding color, spore production, growth rates on cheese, salt tolerance, lipolysis, proteolysis, volatile compound or toxin production and/or competitive ability. The cheese populations also show degeneration for some unused functions such as decreased ability of sexual reproduction or of growth under harsh conditions. These recent findings have fundamental importance for our understanding of adaptation and have applied interest for strain improvement.

Domestication of the Emblematic White Cheese-Making Fungus Penicillium camemberti and Its Diversification into Two Varieties

Domestication involves recent adaptation under strong human selection and rapid diversification and therefore constitutes a good model for studies of these processes. We studied the domestication of the emblematic white mold Penicillium camemberti, used for the maturation of soft cheeses, such as Camembert and Brie, about which surprisingly little was known, despite its economic and cultural importance. Whole-genome-based analyses of genetic relationships and diversity revealed that an ancient domestication event led to the emergence of the gray-green P. biforme mold used in cheese making, by divergence from the blue-green wild P. fuscoglaucum fungus. Another much more recent domestication event led to the generation of the P. camemberti clonal lineage as a sister group to P. biforme. Penicillium biforme displayed signs of phenotypic adaptation to cheese making relative to P. fuscoglaucum, in terms of whiter color, faster growth on cheese medium under cave conditions, lower amounts of toxin production, and greater ability to prevent the growth of other fungi. The P. camemberti lineage displayed even stronger signs of domestication for all these phenotypic features. We also identified two differentiated P. camemberti varieties, apparently associated with different kinds of cheeses and with contrasted phenotypic features in terms of color, growth, toxin production, and competitive ability. We have thus identified footprints of domestication in these fungi, with genetic differentiation between cheese and wild populations, bottlenecks, and specific phenotypic traits beneficial for cheese making. This study has not only fundamental implications for our understanding of domestication but can also have important effects on cheese making.

Genome-wide identification and functional analysis of the horizontally transferred genes in Penicillium

Horizontal gene transfer (HGT) is the transmission of genetic material between different evolutionary lineages and is believed to be an important source of genomic innovation in fungi. In this study, we searched for prokaryotic-derived HGTs in 23 fully sequenced genomes using a comprehensive phylogenomic pipeline followed by manual curation. We found strong support for 60 HGT events comprising 190 genes putatively acquired from bacteria. HGT affected all Penicillium species to various degrees. Gene duplication events happened to 3 HGT genes after the transmission. Most HGT events include genes encoding a variety of enzymes, which are associated with sugar, amino acid, and lipid metabolism. Transcriptome data from 6 Penicillium species revealed that 33 of 35 HGT genes showed expression under the conditions tested and 16 genes were differentially expressed. Our results suggest an important role for inter-domain gene transfers in shaping the genome of Penicillium fungi.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Moulded by humans: The domestication of blue-veined cheese fungi

It is hard to imagine a world without food-associated microbes. The production of bread, wine, beer, salami, coffee, chocolate, cheese and many other foods and beverages all rely on specific microbes. In cheese, myriad microbial species collaborate to yield the complex organoleptic properties that are appreciated by millions of people worldwide. In the early days of cheese making, these complex communities emerged spontaneously from the natural flora associated with the raw materials, the equipment, the production environment or craftsmen involved in the production process. However, in some cases, the microbes shifted their natural habitat to the new cheese-associated environment. The most obvious cause of this is backslopping, where part of a fermented product is used to inoculate the next batch. In addition, some microbes may simply adhere to the tools used in the production process. These microbial communities gradually adapted to the novel man-made niches, a process referred to as "domestication." Domestication is associated with specific genomic and phenotypic changes and ultimately leads to lineages that are genetically and phenotypically distinct from their wild ancestors. In this issue of Molecular Ecology, Dumas et al. have investigated a prime example of cheese-associated microbes, the fungus Penicillium roqueforti. The authors identified several hallmarks of domestication in the genome and phenome of this species, allowing them to hypothesize about the origin of blue-veined cheese fungi domestication, and the specific evolutionary processes involved in adaptation to the cheese matrix.

Phenotypic and Genetic Characterization of the Cheese Ripening Yeast Geotrichum candidum

The yeast Geotrichum candidum (teleomorph Galactomyces candidus) is inoculated onto mold- and smear-ripened cheeses and plays several roles during cheese ripening. Its ability to metabolize proteins, lipids, and organic acids enables its growth on the cheese surface and promotes the development of organoleptic properties. Recent multilocus sequence typing (MLST) and phylogenetic analyses of G. candidum isolates revealed substantial genetic diversity, which may explain its strain-dependant technological capabilities. Here, we aimed to shed light on the phenotypic and genetic diversity among eight G. candidum and three Galactomyces spp. strains of environmental and dairy origin. Phenotypic tests such as carbon assimilation profiles, the ability to grow at 35°C and morphological traits on agar plates allowed us to discriminate G. candidum from Galactomyces spp. The genomes of these isolates were sequenced and assembled; whole genome comparison clustered the G. candidum strains into three subgroups and provided a reliable reference for MLST scheme optimization. Using the whole genome sequence as a reference, we optimized an MLST scheme using six loci that were proposed in two previous MLST schemes. This new MLST scheme allowed us to identify 15 sequence types (STs) out of 41 strains and revealed three major complexes named GeoA, GeoB, and GeoC. The population structure of these 41 strains was evaluated with STRUCTURE and a NeighborNet analysis of the combined six loci, which revealed recombination events between and within the complexes. These results hint that the allele variation conferring the different STs arose from recombination events. Recombination occurred for the six housekeeping genes studied, but most likely occurred throughout the genome. These recombination events may have induced an adaptive divergence between the wild strains and the cheesemaking strains, as observed for other cheese ripening fungi. Further comparative genomic studies are needed to confirm this phenomenon in G. candidum. In conclusion, the draft assembly of 11 G. candidum/Galactomyces spp. genomes allowed us to optimize a genotyping MLST scheme and, combined with the assessment of their ability to grow under different conditions, provides a reliable tool to cluster and eventually improves the selection of G. candidum strains.

A Small Horizontally Transferred Gene Cluster Contributes to the Sporulation of Alternaria alternata

Horizontal gene transfer (HGT) has been identified as an important source of genomic innovation in fungi. However, how HGT drove the evolution of Alternaria alternata, a necrotrophic fungus which can be ubiquitously isolated from soil and various plants and decaying plant materials is largely known. In this study, we identified 12 protein-encoding genes that are likely acquired from lineages outside Pezizomycotina. Phylogenetic trees and approximately unbiased comparative topology tests strongly supported the evolutionary origin of these genes. According to their predicted functions, these HGT candidates are involved in nitrogen and carbohydrate metabolism. Especially, five genes of them were likely transferred as a physically linked cluster from Tremellales (Basidiomycota). Functionally knocking out the five-gene cluster in an A. alternata isolate causing citrus brown spot resulted in an 80% decrease in asexual spore production in the deletion mutant. We further knocked out each of these five genes in this cluster and the resultant single-gene deletion mutants exhibited a various degree of reduction in spore production. Except for conidiation, functions of these genes associated with vegetative growth, stress tolerance, and virulence are very limited. Our results provide new evidence that HGT has played important roles over the course of the evolution of filamentous fungi.

Comparative genomics applied to Mucor species with different lifestyles

BACKGROUND

Despite a growing number of investigations on early diverging fungi, the corresponding lineages have not been as extensively characterized as Ascomycota or Basidiomycota ones. The Mucor genus, pertaining to one of these lineages is not an exception. To this date, a restricted number of Mucor annotated genomes is publicly available and mainly correspond to the reference species, Mucor circinelloides, and to medically relevant species. However, the Mucor genus is composed of a large number of ubiquitous species as well as few species that have been reported to specifically occur in certain habitats. The present study aimed to expand the range of Mucor genomes available and identify potential genomic imprints of adaptation to different environments and lifestyles in the Mucor genus.

RESULTS

In this study, we report four newly sequenced genomes of Mucor isolates collected from non-clinical environments pertaining to species with contrasted lifestyles, namely Mucor fuscus and Mucor lanceolatus, two species used in cheese production (during ripening), Mucor racemosus, a recurrent cheese spoiler sometimes described as an opportunistic animal and human pathogen, and Mucor endophyticus, a plant endophyte. Comparison of these new genomes with those previously available for six Mucor and two Rhizopus (formerly identified as M. racemosus) isolates allowed global structural and functional description such as their TE content, core and species-specific genes and specialized genes. We proposed gene candidates involved in iron metabolism; some of these genes being known to be involved in pathogenicity; and described patterns such as a reduced number of CAZymes in the species used for cheese ripening as well as in the endophytic isolate that might be related to adaptation to different environments and lifestyles within the Mucor genus.

CONCLUSIONS

This study extended the descriptive data set for Mucor genomes, pointed out the complexity of obtaining a robust phylogeny even with multiple genes families and allowed identifying contrasting potentially lifestyle-associated gene repertoires. The obtained data will allow investigating further the link between genetic and its biological data, especially in terms of adaptation to a given habitat.

Independent domestication events in the blue-cheese fungus Penicillium roqueforti

Domestication provides an excellent framework for studying adaptive divergence. Using population genomics and phenotypic assays, we reconstructed the domestication history of the blue cheese mould Penicillium roqueforti. We showed that this fungus was domesticated twice independently. The population used in Roquefort originated from an old domestication event associated with weak bottlenecks and exhibited traits beneficial for pre‐industrial cheese production (slower growth in cheese and greater spore production on bread, the traditional multiplication medium). The other cheese population originated more recently from the selection of a single clonal lineage, was associated with all types of blue cheese worldwide except Roquefort, and displayed phenotypes more suited for industrial cheese production (high lipolytic activity, efficient cheese cavity colonization ability and salt tolerance). We detected genomic regions affected by recent positive selection and putative horizontal gene transfers. This study sheds light on the processes of rapid adaptation and raises questions about genetic resource conservation.

see also the Perspective by Brigida Gallone, Jan Steensels and Kevin J. Verstrepen.

Penicillium camemberti galacturonate reductase: C-1 oxidation/reduction of uronic acids and substrate inhibition mitigation by aldonic acids

The enzyme galacturonate oxidoreductase PcGOR from Penicillium camemberti reduces the C-1 carbon of D-glucuronate and C-4 epimer D-galacturonate to their corresponding aldonic acids, important reactions in both pectin catabolism and ascorbate biosynthesis. PcGOR was active on both glucuronic acid and galacturonic acid, with similar substrate specificities (k_cat/K_m) using the preferred co-substrate NADPH. Substrate acceptance extended to lactone congeners, and D-glucurono-3,6-lactone was converted to L-gulono-1,4-lactone, an immediate precursor of ascorbate. Reaction with glucuronate showed only minor substrate inhibition, and the product L-gulonate and L-gulono-1,4-lactone were both found to be competitive inhibitors with K_i in the low mM range. In contrast, reaction with C-4 epimer galacturonate displayed marked substrate inhibition. Moreover, the product L-galactonate and L-galactono-1,4-lactone were observed to mitigate substrate inhibition by galacturonate, with the lactone having a greater effect than the acid.

Rapid Phenotypic and Metabolomic Domestication of Wild Penicillium Molds on Cheese

Fermented foods provide novel ecological opportunities for natural populations of microbes to evolve through successive recolonization of resource-rich substrates. Comparative genomic data have reconstructed the evolutionary histories of microbes adapted to food environments, but experimental studies directly demonstrating the process of domestication are lacking for most fermented food microbes. Here, we show that during adaptation to cheese, phenotypic and metabolomic traits of wild Penicillium molds rapidly change to produce domesticated phenotypes with properties similar to those of the industrial cultures used to make Camembert and other bloomy rind cheeses. Over a period of just a few weeks, populations of wild Penicillium strains serially passaged on cheese had reduced pigment, spore, and mycotoxin production. Domesticated strains also had a striking change in volatile metabolite production, shifting from production of earthy or musty volatile compounds (e.g., geosmin) to fatty and cheesy volatiles (e.g., 2-nonanone, 2-undecanone). RNA sequencing demonstrated a significant decrease in expression of 356 genes in domesticated strains, with an enrichment of many secondary metabolite production pathways in these downregulated genes. By manipulating the presence of neighboring microbial species and overall resource availability, we demonstrate that the limited competition and high nutrient availability of the cheese environment promote rapid trait evolution of Penicillium molds.IMPORTANCE Industrial cultures of filamentous fungi are used to add unique aesthetics and flavors to cheeses and other microbial foods. How these microbes adapted to live in food environments is generally unknown as most microbial domestication is unintentional. Our work demonstrates that wild molds closely related to the starter culture Penicillium camemberti can readily lose traits and quickly shift toward producing desirable aroma compounds. In addition to experimentally demonstrating a putative domestication pathway for P. camemberti, our work suggests that wild Penicillium isolates could be rapidly domesticated to produce new flavors and aesthetics in fermented foods.