Copy Number Variation in Fungi and Its Implications for Wine Yeast Genetic Diversity and Adaptation

Variability and functional characterization of the Phakopsora pachyrhizi Egh16-like effectors

Effector proteins in Phakopsora pachyrhizi (Pp), the causative agent of Asian Soybean rust, are involved in the infection process. A previous study identified a rust effector Egh16-like family based expression profile during the interaction with soybean. Herein, we scrutinized available the Pp genomes to validate the predicted Egh16-like family of Pp and identify new family members. We described 22 members of the Egh16-like gene family in the Pp MT2006 genome and 18 in the UFV02 and K8108 genomes, highlighting a family expansion. Family members have a small signal peptide, conserved cysteine-rich R/Y/FxC motifs in the C-terminal region, and a virulence-related Egh16-like domain and were able to suppress PTI related responses in Benthamiana. Phylogenetic analysis placed the family members into eight clusters, with members induced during the early stages of rust infection. Members of clusters VI and VII are present in different copy numbers in Pp genomes and suppressed PAMP-related responses.

Selection for robust metabolism in domesticated yeasts is driven by adaptation to Hsp90 stress

Protein folding both promotes and constrains adaptive evolution. We uncover this surprising duality in the role of the protein-folding chaperone heat shock protein 90 (Hsp90) in maintaining the integrity of yeast metabolism amid proteotoxic stressors within industrial domestication niches. Ethanol disrupts critical Hsp90-dependent metabolic pathways and exerts strong selective pressure for redundant duplications of key genes within these pathways, yielding the classical genomic signatures of beer and bread domestication. This work demonstrates a mechanism of adaptive canalization in an ecology of major economic importance and highlights Hsp90-dependent variation as an important source of phantom heritability in complex traits.

Genetic modification of Candida maltosa, a non-pathogenic CTG species, reveals EFG1 function

Candida maltosa is closely related to important pathogenic Candida species, especially C. tropicalis and C. albicans, but it has been rarely isolated from humans. For this reason, through comparative studies, it could be a powerful model to understand the genetic underpinnings of the pathogenicity of Candida species. Here, we generated a cohesive assembly of the C. maltosa genome and developed genetic engineering tools that will facilitate studying this species at a molecular level. We used a combination of short and long-read sequencing to build a polished genomic draft composed of 14 Mbp, 45 contigs and close to 5700 genes. This assembly represents a substantial improvement from the currently available sequences that are composed of thousands of contigs. Genomic comparison with C. albicans and C. tropicalis revealed a substantial reduction in the total number of genes in C. maltosa. However, gene loss seems not to be associated to the avirulence of this species given that most genes that have been previously associated with pathogenicity were also present in C. maltosa. To be able to edit the genome of C. maltosa we generated a set of triple auxotrophic strains so that gene deletions can be performed similarly to what has been routinely done in pathogenic Candida species. As a proof of concept, we generated gene knockouts of EFG1, a gene that encodes a transcription factor that is essential for filamentation and biofilm formation in C. albicans and C. tropicalis. Characterization of these mutants showed that Efg1 also plays a role in biofilm formation and filamentous growth in C. maltosa, but it seems to be a repressor of filamentation in this species. The genome assembly and auxotrophic mutants developed here are a key step forward to start using C. maltosa for comparative and evolutionary studies at a molecular level.

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.

Genetic Diversity and Classification of Colletotrichum sublineola Pathotypes Using a Standard Set of Sorghum Differentials

Anthracnose, incited by Colletotrichum sublineola, is the most destructive foliar disease of sorghum and, under severe conditions, yield losses can exceed 80% on susceptible cultivars. The hyper-variable nature of the pathogen makes its management challenging despite the occurrence of several resistant sources. In this study, the genetic variability and pathogenicity of 140 isolates of C. sublineola, which were sequenced using restriction site-associated sequencing (RAD-Seq), resulted in 1244 quality SNPs. The genetic relationship based on the SNP data showed low to high genetic diversity based on isolates' origin. Isolates from Georgia and North Carolina were grouped into multiple clusters with some level of genetic relationships to each other. Even though some isolates from Texas formed a cluster, others clustered with isolates from Puerto Rico. The isolates from Puerto Rico showed scattered distribution, indicating the diverse nature of these isolates. A population structure and cluster analysis revealed that the genetic variation was stratified into eight populations and one admixture group. The virulence pattern of 30 sequenced isolates on 18 sorghum differential lines revealed 27 new pathotypes. SC748-5, SC112-14, and Brandes were resistant to all the tested isolates, while BTx623 was susceptible to all. Line TAM428 was susceptible to all the pathotypes, except for pathotype 26. Future use of the 18 differentials employed in this study, which contains cultivars/lines which have been used in the Americas, Asia, and Africa, could allow for better characterization of C. sublineola pathotypes at a global level, thus accelerating the development of sorghum lines with stable resistance to the anthracnose pathogen.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Ethanol Drives Evolution of Hsp90-Dependent Robustness by Redundancy in Yeast Domestication

Protein folding promotes and constrains adaptive evolution. We uncover this surprising duality in the role the protein-folding chaperone Hsp90 plays in mediating the interplay between proteome and the genome which acts to maintain the integrity of yeast metabolism in the face of proteotoxic stressors in anthropic niches. Of great industrial relevance, ethanol concentrations generated by fermentation in the making of beer and bread disrupt critical Hsp90-dependent nodes of metabolism and exert strong selective pressure for increased copy number of key genes encoding components of these nodes, yielding the classical genetic signatures of beer and bread domestication. This work establishes a mechanism of adaptive canalization in an ecology of major economic significance and highlights Hsp90-contingent variation as an important source of phantom heritability in complex traits.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Porcine fungal mock community analyses: Implications for mycobiome investigations

Introduction

The gut microbiome is an integral partner in host health and plays a role in immune development, altered nutrition, and pathogen prevention. The mycobiome (fungal microbiome) is considered part of the rare biosphere but is still a critical component in health. Next generation sequencing has improved our understanding of fungi in the gut, but methodological challenges remain. Biases are introduced during DNA isolation, primer design and choice, polymerase selection, sequencing platform selection, and data analyses, as fungal reference databases are often incomplete or contain erroneous sequences.

Methods

Here, we compared the accuracy of taxonomic identifications and abundances from mycobiome analyses which vary among three commonly selected target gene regions (18S, ITS1, or ITS2) and the reference database (UNITE - ITS1, ITS2 and SILVA - 18S). We analyze multiple communities including individual fungal isolates, a mixed mock community created from five common fungal isolates found in weanling piglet feces, a purchased commercial fungal mock community, and piglet fecal samples. In addition, we calculated gene copy numbers for the 18S, ITS1, and ITS2 regions of each of the five isolates from the piglet fecal mock community to determine whether copy number affects abundance estimates. Finally, we determined the abundance of taxa from several iterations of our in-house fecal community to assess the effects of community composition on taxon abundance.

Results

Overall, no marker-database combination consistently outperformed the others. Internal transcribed space markers were slightly superior to 18S in the identification of species in tested communities, but Lichtheimia corymbifera, a common member of piglet gut communities, was not amplified by ITS1 and ITS2 primers. Thus, ITS based abundance estimates of taxa in piglet mock communities were skewed while 18S marker profiles were more accurate. Kazachstania slooffiae displayed the most stable copy numbers (83-85) while L. corymbifera displayed significant variability (90-144) across gene regions.

Discussion

This study underscores the importance of preliminary studies to assess primer combinations and database choice for the mycobiome sample of interest and raises questions regarding the validity of fungal abundance estimates.

A Landscape of the Genomic Structure of Cryptococcus neoformans in Colombian Isolates

Cryptococcus neoformans species complexes are recognized as environmental fungi responsible for lethal meningoencephalitis in immunocompromised individuals. Despite the vast knowledge about the epidemiology and genetic diversity of this fungus in different regions of the world, more studies are necessary to comprehend the genomic profiles across South America, including Colombia, considered to be the second country with the highest number of Cryptococcosis. Here, we sequenced and analyzed the genomic architecture of 29 Colombian C. neoformans isolates and evaluated the phylogenetic relationship of these strains with publicly available C. neoformans genomes. The phylogenomic analysis showed that 97% of the isolates belonged to the VNI molecular type and the presence of sub-lineages and sub-clades. We evidenced a karyotype without changes, a low number of genes with copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs). Additionally, a difference in the number of SNPs between the sub-lineages/sub-clades was observed; some were involved in crucial fungi biological processes. Our study demonstrated the intraspecific divergence of C. neoformans in Colombia. These findings provide evidence that Colombian C. neoformans isolates do not probably require significant structural changes as adaptation mechanisms to the host. To the best of our knowledge, this is the first study to report the whole genome sequence of Colombian C. neoformans isolates.

Genomic characterization of polyextremotolerant black yeasts isolated from food and food production environments

Black yeasts have been isolated from acidic, low water activity, and thermally processed foods as well as from surfaces in food manufacturing plants. The genomic basis for their relative tolerance to food-relevant environmental stresses has not been well defined. In this study, we performed whole genome sequencing (WGS) on seven black yeast strains including Aureobasidium (n=5) and Exophiala (n=2) which were isolated from food or food production environments. These strains were previously characterized for their tolerance to heat, hyperosmotic pressure, high pressure processing, hypochlorite sanitizers, and ultraviolet light. Based on the WGS data, three of the strains previously identified as A. pullulans were reassigned as A. melanogenum. Both haploid and diploid A. melanogenum strains were identified in this collection. Single-locus phylogenies based on beta tubulin, RNA polymerase II, or translation elongation factor protein sequences were compared to the phylogeny produced through SNP analysis, revealing that duplication of the fungal genome in diploid strains complicates the use of single-locus phylogenetics. There was not a strong association between phylogeny and either environmental source or stress tolerance phenotype, nor were trends in the copy numbers of stress-related genes associated with extremotolerance within this collection. While there were obvious differences between the genera, the heterogenous distribution of stress tolerance phenotypes and genotypes suggests that food-relevant black yeasts may be ubiquitous rather than specialists associated with particular ecological niches. However, further evaluation of additional strains and the potential impact of gene sequence modification is necessary to confirm these findings.

Lifestyle Transitions in Fusarioid Fungi are Frequent and Lack Clear Genomic Signatures

The fungal genus Fusarium (Ascomycota) includes well-known plant pathogens that are implicated in diseases worldwide, and many of which have been genome sequenced. The genus also encompasses other diverse lifestyles, including species found ubiquitously as asymptomatic-plant inhabitants (endophytes). Here, we produced structurally annotated genome assemblies for five endophytic Fusarium strains, including the first whole-genome data for Fusarium chuoi. Phylogenomic reconstruction of Fusarium and closely related genera revealed multiple and frequent lifestyle transitions, the major exception being a monophyletic clade of mutualist insect symbionts. Differential codon usage bias and increased codon optimisation separated Fusarium sensu stricto from allied genera. We performed computational prediction of candidate secreted effector proteins (CSEPs) and carbohydrate-active enzymes (CAZymes)-both likely to be involved in the host-fungal interaction-and sought evidence that their frequencies could predict lifestyle. However, phylogenetic distance described gene variance better than lifestyle did. There was no significant difference in CSEP, CAZyme, or gene repertoires between phytopathogenic and endophytic strains, although we did find some evidence that gene copy number variation may be contributing to pathogenicity. Large numbers of accessory CSEPs (i.e., present in more than one taxon but not all) and a comparatively low number of strain-specific CSEPs suggested there is a limited specialisation among plant associated Fusarium species. We also found half of the core genes to be under positive selection and identified specific CSEPs and CAZymes predicted to be positively selected on certain lineages. Our results depict fusarioid fungi as prolific generalists and highlight the difficulty in predicting pathogenic potential in the group.

Isolation, characterization, and genome assembly of Barnettozyma botsteinii sp. nov. and novel strains of Kurtzmaniella quercitrusa isolated from the intestinal tract of the termite Macrotermes bellicosus

Bioconversion of hemicelluloses into simpler sugars leads to the production of a significant amount of pentose sugars, such as d-xylose. However, efficient utilization of pentoses by conventional yeast production strains remains challenging. Wild yeast strains can provide new industrially relevant characteristics and efficiently utilize pentose sugars. To explore this strategy, we isolated gut-residing yeasts from the termite Macrotermes bellicosus collected in Comoé National Park, Côte d'Ivoire. The yeasts were classified through their Internal Transcribed Spacer/Large Subunit sequence, and their genomes were sequenced and annotated. We identified a novel yeast species, which we name Barnettozyma botsteinii sp. nov. 1118T (MycoBank: 833563, CBS 16679T and IBT 710) and two new strains of Kurtzmaniella quercitrusa: var. comoensis (CBS 16678, IBT 709) and var. filamentosus (CBS 16680, IBT 711). The two K. quercitrusa strains grow 15% faster on synthetic glucose medium than Saccharomyces cerevisiae CEN.PKT in acidic conditions (pH = 3.2) and both strains grow on d-xylose as the sole carbon source at a rate of 0.35 h-1. At neutral pH, the yeast form of K. quercitrusa var. filamentosus, but not var. comoensis, switched to filamentous growth in a carbon source-dependent manner. Their genomes are 11.0-13.2 Mb in size and contain between 4888 and 5475 predicted genes. Together with closely related species, we did not find any relationship between gene content and ability to grow on xylose. Besides its metabolism, K. quercitrusa var. filamentosus has a large potential as a production organism, because of its capacity to grow at low pH and to undergo a dimorphic shift.

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.

Distinct expression profiles of regulatory RNAs in the response to biocides in Staphylococcus aureus and Enterococcus faecium

The aim of the study was to characterize the antimicrobial activity of clinically-relevant biocides (chlorhexidine digluconate, benzalkonium chloride, PVP-iodine and triclosan) and to determine the sRNA expression profiles under biocide exposure in two major bacterial opportunistic pathogens, Enterococcus faecium and Staphylococcus aureus. In vitro activities were evaluated against S. aureus HG003 and E. faecium Aus0004. We determined MIC, MBC, sub-inhibitory concentrations (SIC) and growth curves under SIC conditions. sRNA expression study under SIC exposure of biocides was performed by RT-qPCR on 3 sRNAs expressed in S. aureus (RNAIII, SprD and SprX) and the first 9 sRNAs identified as expressed in E. faecium. MICs were higher against E. faecium than for S. aureus. Growth curves under increasing biocide concentrations highlighted two types of bactericidal activity: "on/off" effect for chlorhexidine, benzalkonium chloride, PVP-iodine and a "concentration-dependent" activity for triclosan. Exposure to biocide SICs led to an alteration of several sRNA expression profiles, mostly repressed. The distinct biocide activity profiles must be evaluated with other compounds and bacterial species to enrich the prediction of resistance risks associated with biocide usage. Biocide exposure induces various sRNA-mediated responses in both S. aureus and E. faecium, and further investigations are needed to decipher sRNA-driven regulatory networks.

Characterization of the Ergosterol Biosynthesis Pathway in Ceratocystidaceae

Terpenes represent the biggest group of natural compounds on earth. This large class of organic hydrocarbons is distributed among all cellular organisms, including fungi. The different classes of terpenes produced by fungi are mono, sesqui, di- and triterpenes, although triterpene ergosterol is the main sterol identified in cell membranes of these organisms. The availability of genomic data from members in the Ceratocystidaceae enabled the detection and characterization of the genes encoding the enzymes in the mevalonate and ergosterol biosynthetic pathways. Using a bioinformatics approach, fungal orthologs of sterol biosynthesis genes in nine different species of the Ceratocystidaceae were identified. Ergosterol and some of the intermediates in the pathway were also detected in seven species (Ceratocystis manginecans, C. adiposa, Huntiella moniliformis, Thielaviopsis punctulata, Bretziella fagacearum, Endoconidiophora polonica and Davidsoniella virescens), using gas chromatography-mass spectrometry analysis. The average ergosterol content differed among different genera of Ceratocystidaceae. We also identified all possible terpene related genes and possible biosynthetic clusters in the genomes used in this study. We found a highly conserved terpene biosynthesis gene cluster containing some genes encoding ergosterol biosynthesis enzymes in the analysed genomes. An additional possible terpene gene cluster was also identified in all of the Ceratocystidaceae. We also evaluated the sensitivity of the Ceratocystidaceae to a triazole fungicide that inhibits ergosterol synthesis. The results showed that different members of this family behave differently when exposed to different concentrations of triazole tebuconazole.

Large-scale sequencing and comparative analysis of oenological Saccharomyces cerevisiae strains supported by nanopore refinement of key genomes

Saccharomyces cerevisiae has long been part of human activities related to the production of food and wine. The industrial demand for fermented beverages with well-defined and stable characteristics boosted the isolation and selection of strains conferring a distinctive aroma profile to the final product. To uncover variants characterizing oenological strains, the sequencing of 65 new S. cerevisiae isolates, and the comparison with other 503 publicly available genomes were performed. A hybrid approach based on short Illumina and long Oxford Nanopore reads allowed the in-depth investigation of eleven genomes and the identification of putative laterally transferred regions and structural variants. A comparative analysis between clusters of strains belonging to different datasets allowed the identification of novel relevant genetic features including single nucleotide polymorphisms, insertions and structural variants. Detection of oenological single nucleotide variants shed light on the existence of different levels of modulation for the mevalonate pathway relevant for the biosynthesis of aromatic compounds.

Curing of Cocoa Beans: Fine-Scale Monitoring of the Starter Cultures Applied and Metabolomics of the Fermentation and Drying Steps

Starter culture-initiated cocoa fermentation processes can be applied to improve the quality of cured cocoa beans. However, an accurate monitoring of the microbial strains inoculated in fresh cocoa pulp-bean mass to assess their contribution to the cocoa bean curing process is still lacking. In the present study, eight different cocoa fermentation processes were carried out with Trinitario cocoa in vessels in Costa Rica to assess the contribution of two candidate yeast starter culture strains, namely Saccharomyces cerevisiae IMDO 050523 and Pichia kudriavzevii IMDO 020508, inoculated in combination with Limosilactobacillus fermentum IMDO 0611222 and Acetobacter pasteurianus IMDO 0506386. A multiphasic approach, consisting of culture-dependent selective plating and incubation, rRNA-PCR-DGGE community profiling of agar plate washes, and culture-independent high-throughput amplicon sequencing, combined with a metabolite target analysis of non-volatile and volatile organic compounds (VOCs), was performed on samples from the fermentation and/or drying steps. The different starter culture mixtures applied effectively steered the cocoa fermentation processes performed. Moreover, the use of an amplicon sequence variant (ASV) approach, aligning these ASVs to the whole-genome sequences of the inoculated strains, allowed the monitoring of these inoculated strains and their differentiation from very closely related variants naturally present in the background or spontaneous fermentation processes. Further, traits such as malolactic fermentation during the fermentation step and acetoin and tetramethylpyrazine formation during the drying step could be unraveled. Finally, the yeast strains inoculated influenced the substrate consumption and metabolite production during all starter culture-initiated fermentation processes. This had an impact on the VOC profiles of the cured cocoa beans. Whereas the P. kudriavzevii strain produced a wide range of VOCs in the cocoa pulp, the S. cerevisiae strain mostly influenced the VOC composition of the cured cocoa beans.

Genomic Evidence of an Ancient East Asian Divergence Event in Wild Saccharomyces cerevisiae

Comparative genome analyses have suggested East Asia to be the cradle of the domesticated microbe Brewer's yeast (Saccharomyces cerevisiae), used in the food and biotechnology industry worldwide. Here, we provide seven new, high-quality long-read genomes of nondomesticated yeast strains isolated from primeval forests and other natural environments in China and Taiwan. In a comprehensive analysis of our new genome assemblies, along with other long-read Saccharomycetes genomes available, we show that the newly sequenced East Asian strains are among the closest living relatives of the ancestors of the global diversity of Brewer's yeast, confirming predictions made from short-read genomic data. Three of these strains (termed the East Asian Clade IX Complex here) share a recent ancestry and evolutionary history suggesting an early divergence from other S. cerevisiae strains before the larger radiation of the species, and prior to its domestication. Our genomic analyses reveal that the wild East Asian strains contain elevated levels of structural variations. The new genomic resources provided here contribute to our understanding of the natural diversity of S. cerevisiae, expand the intraspecific genetic variation found in this heavily domesticated microbe, and provide a foundation for understanding its origin and global colonization history.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

Gene Copy Number Variation Does Not Reflect Structure or Environmental Selection in Two Recently Diverged California Populations of Suillus brevipes

Gene copy number variation across individuals has been shown to track population structure and be a source of adaptive genetic variation with significant fitness impacts. In this study, we report opposite results for both predictions based on the analysis of gene copy number variants (CNVs) of Suillus brevipes, a mycorrhizal fungus adapted to coastal and montane habitats in California. In order to assess whether gene copy number variation mirrored population structure and selection in this species, we investigated two previously studied locally adapted populations showing a highly differentiated genomic region encompassing a gene predicted to confer salt tolerance. In addition, we examined whether copy number in the genes related to salt homeostasis was differentiated between the two populations. Although we found many instances of CNV regions across the genomes of S. brevipes individuals, we also found CNVs did not recover population structure and known salt-tolerance-related genes were not under selection across the coastal population. Our results contrast with predictions of CNVs matching single-nucleotide polymorphism divergence and showed CNVs of genes for salt homeostasis are not under selection in S. brevipes.

Crossbreeding of Yeasts Domesticated for Fermentation: Infertility Challenges

Sexual reproduction is almost a universal feature of eukaryotic organisms, which allows the reproduction of new organisms by combining the genetic information from two individuals of different sexes. Based on the mechanism of sexual reproduction, crossbreeding provides an attractive opportunity to improve the traits of animals, plants, and fungi. The budding yeast Saccharomyces cerevisiae has been widely utilized in fermentative production since ancient times. Currently it is still used for many essential biotechnological processes including the production of beer, wine, and biofuels. It is surprising that many yeast strains used in the industry exhibit low rates of sporulation resulting in limited crossbreeding efficiency. Here, I provide an overview of the recent findings about infertility challenges of yeasts domesticated for fermentation along with the progress in crossbreeding technologies. The aim of this review is to create an opportunity for future crossbreeding of yeasts used for fermentation.

Integrated Analysis of Gene Expression, SNP, InDel, and CNV Identifies Candidate Avirulence Genes in Australian Isolates of the Wheat Leaf Rust Pathogen Puccinia triticina

The leaf rust pathogen, Puccinia triticina (Pt), threatens global wheat production. The deployment of leaf rust (Lr) resistance (R) genes in wheat varieties is often followed by the development of matching virulence in Pt due to presumed changes in avirulence (Avr) genes in Pt. Identifying such Avr genes is a crucial step to understand the mechanisms of wheat-rust interactions. This study is the first to develop and apply an integrated framework of gene expression, single nucleotide polymorphism (SNP), insertion/deletion (InDel), and copy number variation (CNV) analysis in a rust fungus and identify candidate avirulence genes. Using a long-read based de novo genome assembly of an isolate of Pt ('Pt104') as the reference, whole-genome resequencing data of 12 Pt pathotypes derived from three lineages Pt104, Pt53, and Pt76 were analyzed. Candidate avirulence genes were identified by correlating virulence profiles with small variants (SNP and InDel) and CNV, and RNA-seq data of an additional three Pt isolates to validate expression of genes encoding secreted proteins (SPs). Out of the annotated 29,043 genes, 2392 genes were selected as SP genes with detectable expression levels. Small variant comparisons between the isolates identified 27-40 candidates and CNV analysis identified 14-31 candidates for each Avr gene, which when combined, yielded the final 40, 64, and 69 candidates for AvrLr1, AvrLr15, and AvrLr24, respectively. Taken together, our results will facilitate future work on experimental validation and cloning of Avr genes. In addition, the integrated framework of data analysis that we have developed and reported provides a more comprehensive approach for Avr gene mining than is currently available.

Replicative aging is associated with loss of genetic heterogeneity from extrachromosomal circular DNA in Saccharomyces cerevisiae

Circular DNA can arise from all parts of eukaryotic chromosomes. In yeast, circular ribosomal DNA (rDNA) accumulates dramatically as cells age, however little is known about the accumulation of other chromosome-derived circles or the contribution of such circles to genetic variation in aged cells. We profiled circular DNA in Saccharomyces cerevisiae populations sampled when young and after extensive aging. Young cells possessed highly diverse circular DNA populations but 94% of the circular DNA were lost after ∼15 divisions, whereas rDNA circles underwent massive accumulation to >95% of circular DNA. Circles present in both young and old cells were characterized by replication origins including circles from unique regions of the genome and repetitive regions: rDNA and telomeric Y' regions. We further observed that circles can have flexible inheritance patterns: [HXT6/7circle] normally segregates to mother cells but in low glucose is present in up to 50% of cells, the majority of which must have inherited this circle from their mother. Interestingly, [HXT6/7circle] cells are eventually replaced by cells carrying stable chromosomal HXT6 HXT6/7 HXT7 amplifications, suggesting circular DNAs are intermediates in chromosomal amplifications. In conclusion, the heterogeneity of circular DNA offers flexibility in adaptation, but this heterogeneity is remarkably diminished with age.

Integrative Omics Analysis Reveals a Limited Transcriptional Shock After Yeast Interspecies Hybridization

The formation of interspecific hybrids results in the coexistence of two diverged genomes within the same nucleus. It has been hypothesized that negative epistatic interactions and regulatory interferences between the two sub-genomes may elicit a so-called genomic shock involving, among other alterations, broad transcriptional changes. To assess the magnitude of this shock in hybrid yeasts, we investigated the transcriptomic differences between a newly formed Saccharomyces cerevisiae × Saccharomyces uvarum diploid hybrid and its diploid parentals, which diverged ∼20 mya. RNA sequencing (RNA-Seq) based allele-specific expression (ASE) analysis indicated that gene expression changes in the hybrid genome are limited, with only ∼1-2% of genes significantly altering their expression with respect to a non-hybrid context. In comparison, a thermal shock altered six times more genes. Furthermore, differences in the expression between orthologous genes in the two parental species tended to be diminished for the corresponding homeologous genes in the hybrid. Finally, and consistent with the RNA-Seq results, we show a limited impact of hybridization on chromatin accessibility patterns, as assessed with assay for transposase-accessible chromatin using sequencing (ATAC-Seq). Overall, our results suggest a limited genomic shock in a newly formed yeast hybrid, which may explain the high frequency of successful hybridization in these organisms.

Into the wild: new yeast genomes from natural environments and new tools for their analysis

Genomic studies of yeasts from the wild have increased considerably in the past few years. This revolution has been fueled by advances in high-throughput sequencing technologies and a better understanding of yeast ecology and phylogeography, especially for biotechnologically important species. The present review aims to first introduce new bioinformatic tools available for the generation and analysis of yeast genomes. We also assess the accumulated genomic data of wild isolates of industrially relevant species, such as Saccharomyces spp., which provide unique opportunities to further investigate the domestication processes associated with the fermentation industry and opportunistic pathogenesis. The availability of genome sequences of other less conventional yeasts obtained from the wild has also increased substantially, including representatives of the phyla Ascomycota (e.g. Hanseniaspora) and Basidiomycota (e.g. Phaffia). Here, we review salient examples of both fundamental and applied research that demonstrate the importance of continuing to sequence and analyze genomes of wild yeasts.

Genomic features and evolution of the conditionally dispensable chromosome in the tangerine pathotype of Alternaria alternata

The tangerine pathotype of the ascomycete fungus Alternaria alternata is the causal agent of citrus brown spot, which can result in significant losses of both yield and marketability for tangerines worldwide. A conditionally dispensable chromosome (CDC), which harbours the host-selective ACT toxin gene cluster, is required for tangerine pathogenicity of A. alternata. To understand the genetic makeup and evolution of the tangerine pathotype CDC, we isolated and sequenced the CDCs of the A. alternata Z7 strain and analysed the function and evolution of their genes. The A. alternata Z7 strain has two CDCs (~1.1 and ~0.8 Mb, respectively), and the longer Z7 CDC contains all but one contig of the shorter one. Z7 CDCs contain 254 predicted protein-coding genes, which are enriched in functional categories associated with 'metabolic process' (55 genes, P = 0.037). Relatively few of the CDC genes can be classified as carbohydrate-active enzymes (CAZymes) (4) and transporters (19) and none as kinases. Evolutionary analysis of the 254 CDC proteins showed that their evolutionary conservation tends to be restricted within the genus Alternaria and that the CDC genes evolve faster than genes in the essential chromosomes, likely due to fewer selective constraints. Interestingly, phylogenetic analysis suggested that four of the 25 genes responsible for the ACT toxin production were likely transferred from Colletotrichum (Sordariomycetes). Functional experiments showed that two of them are essential for the virulence of the tangerine pathotype of A. alternata. These results provide new insights into the function and evolution of CDC genes in Alternaria.

What can be learned by scanning the genome for molecular convergence in wild populations?

Convergent evolution, where independent lineages evolve similar phenotypes in response to similar challenges, can provide valuable insight into how selection operates and the limitations it encounters. However, it has only recently become possible to explore how convergent evolution is reflected at the genomic level. The overlapping outlier approach (OOA), where genome scans of multiple independent lineages are used to find outliers that overlap and therefore identify convergently evolving loci, is becoming popular. Here, we present a quantitative analysis of 34 studies that used this approach across many sampling designs, taxa, and sampling intensities. We found that OOA studies with increased biological sampling power within replicates have increased likelihood of finding overlapping, "convergent" signals of adaptation between them. When identifying convergent loci as overlapping outliers, it is tempting to assume that any false-positive outliers derived from individual scans will fail to overlap across replicates, but this cannot be guaranteed. We highlight how population demographics and genomic context can contribute toward both true convergence and false positives in OOA studies. We finish with an exploration of emerging methods that couple genome scans with phenotype and environmental measures, leveraging added information from genome data to more directly test hypotheses of the likelihood of convergent evolution.

Polar bear evolution is marked by rapid changes in gene copy number in response to dietary shift

Polar bear (Ursus maritimus) and brown bear (Ursus arctos) are recently diverged species that inhabit vastly differing habitats. Thus, analysis of the polar bear and brown bear genomes represents a unique opportunity to investigate the evolutionary mechanisms and genetic underpinnings of rapid ecological adaptation in mammals. Copy number (CN) differences in genomic regions between closely related species can underlie adaptive phenotypes and this form of genetic variation has not been explored in the context of polar bear evolution. Here, we analyzed the CN profiles of 17 polar bears, 9 brown bears, and 2 black bears (Ursus americanus). We identified an average of 318 genes per individual that showed evidence of CN variation (CNV). Nearly 200 genes displayed species-specific CN differences between polar bear and brown bear species. Principal component analysis of gene CN provides strong evidence that CNV evolved rapidly in the polar bear lineage and mainly resulted in CN loss. Olfactory receptors composed 47% of CN differentiated genes, with the majority of these genes being at lower CN in the polar bear. Additionally, we found significantly fewer copies of several genes involved in fatty acid metabolism as well as AMY1B, the salivary amylase-encoding gene in the polar bear. These results suggest that natural selection shaped patterns of CNV in response to the transition from an omnivorous to primarily carnivorous diet during polar bear evolution. Our analyses of CNV shed light on the genomic underpinnings of ecological adaptation during polar bear evolution.

Extensive loss of cell-cycle and DNA repair genes in an ancient lineage of bipolar budding yeasts

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

A population genomic characterization of copy number variation in the opportunistic fungal pathogen Aspergillus fumigatus

Aspergillus fumigatus is a potentially deadly opportunistic fungal pathogen. Molecular studies have shaped our understanding of the genes, proteins, and molecules that contribute to A. fumigatus pathogenicity, but few studies have characterized genome-wide patterns of genetic variation at the population level. Of A. fumigatus genomic studies to-date, most focus mainly on single nucleotide polymorphisms and large structural variants, while overlooking the contribution of copy number variation (CNV). CNV is a class of small structural variation defined as loci that vary in their number of copies between individuals due to duplication, gain, or deletion. CNV can influence phenotype, including fungal virulence. In the present study, we characterized the population genomic patterns of CNV in a diverse collection of 71 A. fumigatus isolates using publicly available sequencing data. We used genome-wide single nucleotide polymorphisms to infer the population structure of these isolates and identified three populations consisting of at least 8 isolates. We then computationally predicted genome-wide CNV profiles for each isolate and conducted analyses at the species-, population-, and individual levels. Our results suggest that CNV contributes to genetic variation in A. fumigatus, with ~10% of the genome being CN variable. Our analysis indicates that CNV is non-randomly distributed across the A. fumigatus genome, and is overrepresented in subtelomeric regions. Analysis of gene ontology categories in genes that overlapped CN variants revealed an enrichment of genes related to transposable element and secondary metabolism functions. We further identified 72 loci containing 33 genes that showed divergent copy number profiles between the three A. fumigatus populations. Many of these genes encode proteins that interact with the cell surface or are involved in pathogenicity. Our results suggest that CNV is an important source of genetic variation that could account for some of the phenotypic differences between A. fumigatus populations and isolates.