Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

MIKE: an ultrafast, assembly-, and alignment-free approach for phylogenetic tree construction

MOTIVATION

Constructing a phylogenetic tree requires calculating the evolutionary distance between samples or species via large-scale resequencing data, a process that is both time-consuming and computationally demanding. Striking the right balance between accuracy and efficiency is a significant challenge.

RESULTS

To address this, we introduce a new algorithm, MIKE (MinHash-based k-mer algorithm). This algorithm is designed for the swift calculation of the Jaccard coefficient directly from raw sequencing reads and enables the construction of phylogenetic trees based on the resultant Jaccard coefficient. Simulation results highlight the superior speed of MIKE compared to existing state-of-the-art methods. We used MIKE to reconstruct a phylogenetic tree, incorporating 238 yeast, 303 Zea, 141 Ficus, 67 Oryza, and 43 Saccharum spontaneum samples. MIKE demonstrated accurate performance across varying evolutionary scales, reproductive modes, and ploidy levels, proving itself as a powerful tool for phylogenetic tree construction.

AVAILABILITY AND IMPLEMENTATION

MIKE is publicly available on Github at https://github.com/Argonum-Clever2/mike.git.

Macroevolutionary dynamics of gene family gain and loss along multicellular eukaryotic lineages

The gain and loss of genes fluctuate over evolutionary time in major eukaryotic clades. However, the full profile of these macroevolutionary trajectories is still missing. To give a more inclusive view on the changes in genome complexity across the tree of life, here we recovered the evolutionary dynamics of gene family gain and loss ranging from the ancestor of cellular organisms to 352 eukaryotic species. We show that in all considered lineages the gene family content follows a common evolutionary pattern, where the number of gene families reaches the highest value at a major evolutionary and ecological transition, and then gradually decreases towards extant organisms. This supports theoretical predictions and suggests that the genome complexity is often decoupled from commonly perceived organismal complexity. We conclude that simplification by gene family loss is a dominant force in Phanerozoic genomes of various lineages, probably underpinned by intense ecological specializations and functional outsourcing.

The metabolic domestication syndrome of budding yeast

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Gene loss and cis-regulatory novelty shaped core histone gene evolution in the apiculate yeast Hanseniaspora uvarum

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage.

Saccharomycotina yeasts defy long-standing macroecological patterns

The Saccharomycotina yeasts ("yeasts" hereafter) are a fungal clade of scientific, economic, and medical significance. Yeasts are highly ecologically diverse, found across a broad range of environments in every biome and continent on earth; however, little is known about what rules govern the macroecology of yeast species and their range limits in the wild. Here, we trained machine learning models on 12,816 terrestrial occurrence records and 96 environmental variables to infer global distribution maps at ~1 km2 resolution for 186 yeast species (~15% of described species from 75% of orders) and to test environmental drivers of yeast biogeography and macroecology. We found that predicted yeast diversity hotspots occur in mixed montane forests in temperate climates. Diversity in vegetation type and topography were some of the greatest predictors of yeast species richness, suggesting that microhabitats and environmental clines are key to yeast diversity. We further found that range limits in yeasts are significantly influenced by carbon niche breadth and range overlap with other yeast species, with carbon specialists and species in high-diversity environments exhibiting reduced geographic ranges. Finally, yeasts contravene many long-standing macroecological principles, including the latitudinal diversity gradient, temperature-dependent species richness, and a positive relationship between latitude and range size (Rapoport's rule). These results unveil how the environment governs the global diversity and distribution of species in the yeast subphylum. These high-resolution models of yeast species distributions will facilitate the prediction of economically relevant and emerging pathogenic species under current and future climate scenarios.

Perspectives on polarity - exploring biological asymmetry across scales

In this Perspective, Journal of Cell Science invited researchers working on cell and tissue polarity to share their thoughts on unique, emerging or open questions relating to their field. The goal of this article is to feature 'voices' from scientists around the world and at various career stages, to bring attention to innovative and thought-provoking topics of interest to the cell biology community. These voices discuss intriguing questions that consider polarity across scales, evolution, development and disease. What can yeast and protists tell us about the evolution of cell and tissue polarity in animals? How are cell fate and development influenced by emerging dynamics in cell polarity? What can we learn from atypical and extreme polarity systems? How can we arrive at a more unified biophysical understanding of polarity? Taken together, these pieces demonstrate the broad relevance of the fascinating phenomenon of cell polarization to diverse fundamental biological questions.

Evolutionary rate covariation is a reliable predictor of co-functional interactions but not necessarily physical interactions

Co-functional proteins tend to have rates of evolution that covary over time. This correlation between evolutionary rates can be measured over the branches of a phylogenetic tree through methods such as evolutionary rate covariation (ERC), and then used to construct gene networks by the identification of proteins with functional interactions. The cause of this correlation has been hypothesized to result from both compensatory coevolution at physical interfaces and nonphysical forces such as shared changes in selective pressure. This study explores whether coevolution due to compensatory mutations has a measurable effect on the ERC signal. We examined the difference in ERC signal between physically interacting protein domains within complexes compared to domains of the same proteins that do not physically interact. We found no generalizable relationship between physical interaction and high ERC, although a few complexes ranked physical interactions higher than nonphysical interactions. Therefore, we conclude that coevolution due to physical interaction is weak, but present in the signal captured by ERC, and we hypothesize that the stronger signal instead comes from selective pressures on the protein as a whole and maintenance of the general function.

Oxygenation influences xylose fermentation and gene expression in the yeast genera Spathaspora and Scheffersomyces

BACKGROUND

Cost-effective production of biofuels from lignocellulose requires the fermentation of D-xylose. Many yeast species within and closely related to the genera Spathaspora and Scheffersomyces (both of the order Serinales) natively assimilate and ferment xylose. Other species consume xylose inefficiently, leading to extracellular accumulation of xylitol. Xylitol excretion is thought to be due to the different cofactor requirements of the first two steps of xylose metabolism. Xylose reductase (XR) generally uses NADPH to reduce xylose to xylitol, while xylitol dehydrogenase (XDH) generally uses NAD+ to oxidize xylitol to xylulose, creating an imbalanced redox pathway. This imbalance is thought to be particularly consequential in hypoxic or anoxic environments.

RESULTS

We screened the growth of xylose-fermenting yeast species in high and moderate aeration and identified both ethanol producers and xylitol producers. Selected species were further characterized for their XR and XDH cofactor preferences by enzyme assays and gene expression patterns by RNA-Seq. Our data revealed that xylose metabolism is more redox balanced in some species, but it is strongly affected by oxygen levels. Under high aeration, most species switched from ethanol production to xylitol accumulation, despite the availability of ample oxygen to accept electrons from NADH. This switch was followed by decreases in enzyme activity and the expression of genes related to xylose metabolism, suggesting that bottlenecks in xylose fermentation are not always due to cofactor preferences. Finally, we expressed XYL genes from multiple Scheffersomyces species in a strain of Saccharomyces cerevisiae. Recombinant S. cerevisiae expressing XYL1 from Scheffersomyces xylosifermentans, which encodes an XR without a cofactor preference, showed improved anaerobic growth on xylose as the primary carbon source compared to S. cerevisiae strain expressing XYL genes from Scheffersomyces stipitis.

CONCLUSION

Collectively, our data do not support the hypothesis that xylitol accumulation occurs primarily due to differences in cofactor preferences between xylose reductase and xylitol dehydrogenase; instead, gene expression plays a major role in response to oxygen levels. We have also identified the yeast Sc. xylosifermentans as a potential source for genes that can be engineered into S. cerevisiae to improve xylose fermentation and biofuel production.

Evolutionary diversity of the control of the azole response by Tra1 across yeast species

Tra1 is an essential coactivator protein of the yeast SAGA and NuA4 acetyltransferase complexes that regulate gene expression through multiple mechanisms including the acetylation of histone proteins. Tra1 is a pseudokinase of the PIKK family characterized by a C-terminal PI3K domain with no known kinase activity. However, mutations of specific arginine residues to glutamine in the PI3K domains (an allele termed tra1Q3) result in reduced growth and increased sensitivity to multiple stresses. In the opportunistic fungal pathogen Candida albicans, the tra1Q3 allele reduces pathogenicity and increases sensitivity to the echinocandin antifungal drug caspofungin, which disrupts the fungal cell wall. Here, we found that compromised Tra1 function, in contrast to what is seen with caspofungin, increases tolerance to the azole class of antifungal drugs, which inhibits ergosterol synthesis. In C. albicans, tra1Q3 increases the expression of genes linked to azole resistance, such as ERG11 and CDR1. CDR1 encodes a multidrug ABC transporter associated with efflux of multiple xenobiotics, including azoles. Consequently, cells carrying tra1Q3 show reduced intracellular accumulation of fluconazole. In contrast, a tra1Q3 Saccharomyces cerevisiae strain displayed opposite phenotypes: decreased tolerance to azole, decreased expression of the efflux pump PDR5, and increased intracellular accumulation of fluconazole. Therefore, our data provide evidence that Tra1 differentially regulates the antifungal response across yeast species.

Exploring future applications of the apiculate yeast Hanseniaspora

As a metaphor, lemons get a bad rap; however the proverb 'if life gives you lemons, make lemonade' is often used in a motivational context. The same could be said of Hanseniaspora in winemaking. Despite its predominance in vineyards and grape must, this lemon-shaped yeast is underappreciated in terms of its contribution to the overall sensory profile of fine wine. Species belonging to this apiculate yeast are known for being common isolates not just on grape berries, but on many other fruits. They play a critical role in the early stages of a fermentation and can influence the quality of the final product. Their deliberate addition within mixed-culture fermentations shows promise in adding to the complexity of a wine and thus provide sensorial benefits. Hanseniaspora species are also key participants in the fermentations of a variety of other foodstuffs ranging from chocolate to apple cider. Outside of their role in fermentation, Hanseniaspora species have attractive biotechnological possibilities as revealed through studies on biocontrol potential, use as a whole-cell biocatalyst and important interactions with Drosophila flies. The growing amount of 'omics data on Hanseniaspora is revealing interesting features of the genus that sets it apart from the other Ascomycetes. This review collates the fields of research conducted on this apiculate yeast genus.

Ogataea nonmethanolica f.a, sp. nov., a novel yeast species isolated from rotting wood in Brazil and Colombia

Three yeast isolate candidates for a novel species were obtained from rotting wood samples collected in Brazil and Colombia. The Brazilian isolate differs from the Colombian isolates by one nucleotide substitution in each of the D1/D2 and small subunit (SSU) sequences. The internal transcribed spacer (ITS) and translation elongation factor 1-α gene sequences of the three isolates were identical. A phylogenetic analysis showed that this novel species belongs to the genus Ogataea. This novel species is phylogenetically related to Candida nanaspora and Candida nitratophila. The novel species differs from C. nanaspora by seven nucleotides and two indels, and by 17 nucleotides and four indels from C. nitratophila in the D1/D2 sequences. The ITS sequences of these three species differ by more than 30 nucleotides. Analyses of the sequences of the SSU and translation elongation factor 1-α gene also showed that these isolates represent a novel species of the genus Ogataea. Different from most Ogataea species, these isolates did not assimilate methanol as the sole carbon source. The name Ogataea nonmethanolica sp. nov. is proposed to accommodate these isolates. The holotype of Ogataea nonmethanolica is CBS 13485T. The MycoBank number is MB 851195.

Taxonomic revision of Geotrichum and Magnusiomyces, with the descriptions of five new Geotrichum species from China

The arthroconidial yeast-like species currently classified in the asexual genera Geotrichum and Saprochaete and the sexual genera Dipodascus, Galactomyces and Magnusiomyces are frequently associated with dairy and cosmetics production, fruit rot and human infection. However, the taxonomic system of these fungi has not been updated to accommodate the new nomenclature code adopting the "one fungus, one name" principle. Here, we performed phylogenetic analyses of these yeast-like species based on the sequences of the internal transcribed spacer (ITS) region and the D1/D2 domain of the large subunit of the rRNA gene. Two monophyletic groups were recognised from these species. One group contained Dipodascus, Galactomyces, and Geotrichum species and the other Magnusiomyces and Saprochaete species. We thus assigned the species in each group into one genus and selected the genus name Geotrichum for the first group and Magnusiomyces for the second one based on the principle of priority of publication. Five new Geotrichum species were identified from arthroconidial yeast strains recently isolated from various sources in China. The new species are described as Ge. dehoogii sp. nov., Ge. fujianense sp. nov., Ge. maricola sp. nov., Ge. smithiae sp. nov., and Ge. sinensis sp. nov.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Natural trait variation across Saccharomycotina species

Among molecular biologists, the group of fungi called Saccharomycotina is famous for its yeasts. These yeasts in turn are famous for what they have in common-genetic, biochemical, and cell-biological characteristics that serve as models for plants and animals. But behind the apparent homogeneity of Saccharomycotina species lie a wealth of differences. In this review, we discuss traits that vary across the Saccharomycotina subphylum. We describe cases of bright pigmentation; a zoo of cell shapes; metabolic specialties; and species with unique rules of gene regulation. We discuss the genetics of this diversity and why it matters, including insights into basic evolutionary principles with relevance across Eukarya.

The cell morphological diversity of Saccharomycotina yeasts

The ∼1 200 known species in subphylum Saccharomycotina are a highly diverse clade of unicellular fungi. During its lifecycle, a typical yeast exhibits multiple cell types with various morphologies; these morphologies vary across Saccharomycotina species. Here, we synthesize the evolutionary dimensions of variation in cellular morphology of yeasts across the subphylum, focusing on variation in cell shape, cell size, type of budding, and filament production. Examination of 332 representative species across the subphylum revealed that the most common budding cell shapes are ovoid, spherical, and ellipsoidal, and that their average length and width is 5.6 µm and 3.6 µm, respectively. 58.4% of yeast species examined can produce filamentous cells, and 87.3% of species reproduce asexually by multilateral budding, which does not require utilization of cell polarity for mitosis. Interestingly, ∼1.8% of species examined have not been observed to produce budding cells, but rather only produce filaments of septate hyphae and/or pseudohyphae. 76.9% of yeast species examined have sexual cycle descriptions, with most producing one to four ascospores that are most commonly hat-shaped (37.4%). Systematic description of yeast cellular morphological diversity and reconstruction of its evolution promises to enrich our understanding of the evolutionary cell biology of this major fungal lineage.

RetroGREAT signaling: The lessons we learn from yeast

The mitochondrial retrograde signaling (RTG) pathway of communication from mitochondria to the nucleus was first studied in yeast Saccharomyces cerevisiae. It rewires cellular metabolism according to the mitochondrial state by reprogramming nuclear gene expression in response to mitochondrial triggers. The main players involved in retrograde signaling are the Rtg1 and Rtg3 transcription factors, and a set of positive and negative regulators, including the Rtg2, Mks1, Lst8, and Bmh1/2 proteins. Retrograde regulation is integrated with other processes, including stress response, osmoregulation, and nutrient sensing through functional crosstalk with cellular pathways such as high osmolarity glycerol or target of rapamycin signaling. In this review, we summarize metabolic changes observed upon retrograde stimulation and analyze the progress made to uncover the mechanisms underlying the integration of regulatory circuits. Comparisons of the evolutionary adaptations of the retrograde pathway that have occurred in the different yeast groups can help to fully understand the process.

Starmerella kisarazuensis f.a., sp. nov., a novel yeast isolated from Trifolium pratense flowers

Three yeast isolates, NBRC 115909T, NBRC 115910 and NBRC 116270, were isolated from Trifolium pratense (red clover) flowers collected from Kisarazu, Chiba, Japan. Analysis of the sequences of the D1/D2 domains of the large subunit (LSU) rRNA gene and the internal transcribed spacer (ITS) regions revealed that these isolates represent a single novel species within the genus Starmerella. Also, no ascospore formation was observed. The yeast isolates were closely related to Starmerella vitae UWOPS 00-107.2T and Starmerella bombi NRRL Y-17081T. They differed from S. vitae, the most closely related species with a validly published name, by ten nucleotide substitutions with two gaps in the D1/D2 domains and 20 nucleotide substitutions in the ITS region. Moreover, the three isolates exhibited distinct phenotypic characteristics from the closely related species. Therefore, we suggest that these three isolates represent a novel species, designated as Starmerella kisarazuensis f.a., sp. nov. The holotype is NBRC 115909T (isotype: CBS 18485T).

RNA processing and degradation mechanisms shaping the mitochondrial transcriptome of budding yeasts

Yeast mitochondrial genes are expressed as polycistronic transcription units that contain RNAs from different classes and show great evolutionary variability. The promoters are simple, and transcriptional control is rudimentary. Posttranscriptional mechanisms involving RNA maturation, stability, and degradation are thus the main force shaping the transcriptome and determining the expression levels of individual genes. Primary transcripts are fragmented by tRNA excision by RNase P and tRNase Z, additional processing events occur at the dodecamer site at the 3' end of protein-coding sequences. groups I and II introns are excised in a self-splicing reaction that is supported by protein splicing factors encoded by the nuclear genes, or by the introns themselves. The 3'-to-5' exoribonucleolytic complex called mtEXO is the main RNA degradation activity involved in RNA turnover and processing, supported by an auxiliary 5'-to-3' exoribonuclease Pet127p. tRNAs and, to a lesser extent, rRNAs undergo several different base modifications. This complex gene expression system relies on the coordinated action of mitochondrial and nuclear genes and undergoes rapid evolution, contributing to speciation events. Moving beyond the classical model yeast Saccharomyces cerevisiae to other budding yeasts should provide important insights into the coevolution of both genomes that constitute the eukaryotic genetic system.

Inference of phylogenetic trees directly from raw sequencing reads using Read2Tree

Current methods for inference of phylogenetic trees require running complex pipelines at substantial computational and labor costs, with additional constraints in sequencing coverage, assembly and annotation quality, especially for large datasets. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes and bypasses traditional steps in phylogeny inference, such as genome assembly, annotation and all-versus-all sequence comparisons, while retaining accuracy. In a benchmark encompassing a broad variety of datasets, Read2Tree is 10-100 times faster than assembly-based approaches and in most cases more accurate-the exception being when sequencing coverage is high and reference species very distant. Here, to illustrate the broad applicability of the tool, we reconstruct a yeast tree of life of 435 species spanning 590 million years of evolution. We also apply Read2Tree to >10,000 Coronaviridae samples, accurately classifying highly diverse animal samples and near-identical severe acute respiratory syndrome coronavirus 2 sequences on a single tree. The speed, accuracy and versatility of Read2Tree enable comparative genomics at scale.

Dissecting an ancient stress resistance trait syndrome in the compost yeast Kluyveromyces marxianus

In the search to understand how evolution builds new traits, ancient events are often the hardest to dissect. Species-unique traits pose a particular challenge for geneticists-cases in which a character arose long ago and, in the modern day, is conserved within a species, distinguishing it from reproductively isolated relatives. In this work, we have developed the budding yeast genus Kluyveromyces as a model for mechanistic dissection of trait variation across species boundaries. Phenotypic profiling revealed robust heat and chemical-stress tolerance phenotypes that distinguished the compost yeast K. marxianus from the rest of the clade. We used culture-based, transcriptomic, and genetic approaches to characterize the metabolic requirements of the K. marxianus trait syndrome. We then generated a population-genomic resource for K. marxianus and harnessed it in molecular-evolution analyses, which found hundreds of housekeeping genes with evidence for adaptive protein variation unique to this species. Our data support a model in which, in the distant past, K. marxianus underwent a vastly complex remodeling of its proteome to achieve stress resistance. Such a polygenic architecture, involving nucleotide-level allelic variation on a massive scale, is consistent with theoretical models of the mechanisms of long-term adaptation, and suggests principles of broad relevance for interspecies trait genetics.

Changes in fungal taxonomy: mycological rationale and clinical implications

Numerous fungal species of medical importance have been recently subjected to and will likely continue to undergo nomenclatural changes as a result of the application of molecular approaches to fungal classification together with abandonment of dual nomenclature. Here, we summarize those changes affecting key groups of fungi of medical importance, explaining the mycological (taxonomic) rationale that underpinned the changes and the clinical relevance/importance (where such exists) of the key nomenclatural revisions. Potential mechanisms to mitigate unnecessary taxonomic instability are suggested, together with approaches to raise awareness of important changes to minimize potential clinical confusion.

Parallel Nonfunctionalization of CK1δ/ε Kinase Ohnologs Following a Whole-Genome Duplication Event

Whole-genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in Saccharomyces cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post-WGD. In S. cerevisiae, HRR25 encodes the casein kinase 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post-WGD and non-WGD species, and have nonconserved cellular localization, consistent with their ongoing loss of function. The analysis in Naumovozyma castellii shows that the noncomplementing ohnolog is expressed at a lower level and has become nonessential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Functional plasticity in chromosome-microtubule coupling on the evolutionary time scale

The Dam1 complex is essential for mitotic progression across evolutionarily divergent fungi. Upon analyzing amino acid (aa) sequences of Dad2, a Dam1 complex subunit, we identified a conserved 10-aa-long Dad2 signature sequence (DSS). An arginine residue (R126) in the DSS is essential for viability in Saccharomyces cerevisiae that possesses point centromeres. The corresponding arginine residues are functionally important but not essential for viability in Candida albicans and Cryptococcus neoformans; both carry several kilobases long regional centromeres. The purified recombinant Dam1 complex containing either Dad2ΔDSS or Dad2R126A failed to bind microtubules (MTs) or form any visible rings like the WT complex. Intriguingly, functional analysis revealed that the requirement of the conserved arginine residue for chromosome biorientation and mitotic progression reduced with increasing centromere length. We propose that plasticity of the invariant arginine of Dad2 in organisms with regional centromeres is achieved by conditional elevation of the kinetochore protein(s) to enable multiple kinetochore MTs to bind to each chromosome. The capacity of a chromosome to bind multiple kinetochore MTs may mask the deleterious effects of such lethal mutations.

High nitrogen concentration causes G2/M arrest in Hanseniaspora vineae

Yeasts have been widely used as a model to better understand cell cycle mechanisms and how nutritional and genetic factors can impact cell cycle progression. While nitrogen scarcity is well known to modulate cell cycle progression, the relevance of nitrogen excess for microorganisms has been overlooked. In our previous work, we observed an absence of proper entry into the quiescent state in Hanseniaspora vineae and identified a potential link between this behavior and nitrogen availability. Furthermore, the Hanseniaspora genus has gained attention due to a significant loss of genes associated with DNA repair and cell cycle. Thus, the aim of our study was to investigate the effects of varying nitrogen concentrations on H. vineae's cell cycle progression. Our findings demonstrated that nitrogen excess, regardless of the source, disrupts cell cycle progression and induces G2/M arrest in H. vineae after reaching the stationary phase. Additionally, we observed a viability decline in H. vineae cells in an ammonium-dependent manner, accompanied by increased production of reactive oxygen species, mitochondrial hyperpolarization, intracellular acidification, and DNA fragmentation. Overall, our study highlights the events of the cell cycle arrest in H. vineae induced by nitrogen excess and attempts to elucidate the possible mechanism triggering this absence of proper entry into the quiescent state.

Taxogenomic analysis of a novel yeast species isolated from soil, Pichia galeolata sp. nov

A novel budding yeast species was isolated from a soil sample collected in the United States of America. Phylogenetic analyses of multiple loci and phylogenomic analyses conclusively placed the species within the genus Pichia. Strain yHMH446 falls within a clade that includes Pichia norvegensis, Pichia pseudocactophila, Candida inconspicua, and Pichia cactophila. Whole genome sequence data were analyzed for the presence of genes known to be important for carbon and nitrogen metabolism, and the phenotypic data from the novel species were compared to all Pichia species with publicly available genomes. Across the genus, including the novel species candidate, we found that the inability to use many carbon and nitrogen sources correlated with the absence of metabolic genes. Based on these results, Pichia galeolata sp. nov. is proposed to accommodate yHMH446T (=NRRL Y-64187 = CBS 16864). This study shows how integrated taxogenomic analysis can add mechanistic insight to species descriptions.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Biological factors and statistical limitations prevent detection of most noncanonical proteins by mass spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry (MS) experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here, we leveraged recent advances in ribosome profiling and MS to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly expressed to be detected by shotgun MS at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for 4 noncanonical proteins in MS data, which were also supported by evolution and translation data. These results illustrate the power of MS to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly expressed proteins.

Taxogenomic placement of Rasporella oleae and Rasporella dianae gen. and spp. nov., two insect associated yeast species

During the course of independent studies in Europe, North America, and Africa, seven yeast strains were isolated from insect frass, decaying wood, tree flux, and olive oil sediment. Phylogenetic analysis of two barcoding DNA regions (internal transcribed spacer and the D1/D2 domain of the LSU rRNA gene) revealed that they belong to two closely related undescribed species distinct from all genera in the family Debaryomycetaceae. For reliable taxonomic placement the genomes of four strains of the two novel species and six type strains of closely related species were sequenced. Orthologous genes from 54 genomes of representatives of the Pichiomycetes and 23 outgroup taxa were concatenated to construct a fully supported phylogenetic tree. Consistent with the assumptions, we found that the two new species belong to a novel genus. In addition, the delimitation of the novel species was supported by genetic distance calculations from average nucleotide identity (ANI) and digital DNA:DNA hybridization (dDDH) values. The physiological characterization of the novel species was generally consistent with their genomic content. All strains had two alleles encoding secretory lipase in either two or three copies depending on the species. However, lipolytic activity was detected only in strains with three copies of the secretory lipase gene. Nevertheless, lipolytic activity might be related to their association with the insect gut. Based on these results, formal descriptions of the new genus Rasporella gen. nov. and of two new species Rasporella dianae sp. nov. (holotype UCDFST 68-643T , MycoBank no.: 850238) and Rasporella oleae sp. nov. (holotype ZIM 2471T , MycoBank no.: 850126) are provided.

Mitochondrial genome diversity across the subphylum Saccharomycotina

Introduction

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.

Methods

By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.

Results

Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.

Discussion

As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

Functional and evolutionary integration of a fungal gene with a bacterial operon

Siderophores are crucial for iron-scavenging in microorganisms. While many yeasts can uptake siderophores produced by other organisms, they are typically unable to synthesize siderophores themselves. In contrast, Wickerhamiella/Starmerella (W/S) clade yeasts gained the capacity to make the siderophore enterobactin following the remarkable horizontal acquisition of a bacterial operon enabling enterobactin synthesis. Yet, how these yeasts absorb the iron bound by enterobactin remains unresolved. Here, we demonstrate that Enb1 is the key enterobactin importer in the W/S-clade species Starmerella bombicola. Through phylogenomic analyses, we show that ENB1 is present in all W/S clade yeast species that retained the enterobactin biosynthetic genes. Conversely, it is absent in species that lost the ent genes, except for Starmerella stellata, making this species the only cheater in the W/S clade that can utilize enterobactin without producing it. Through phylogenetic analyses, we infer that ENB1 is a fungal gene that likely existed in the W/S clade prior to the acquisition of the ent genes and subsequently experienced multiple gene losses and duplications. Through phylogenetic topology tests, we show that ENB1 likely underwent horizontal gene transfer from an ancient W/S clade yeast to the order Saccharomycetales, which includes the model yeast Saccharomyces cerevisiae, followed by extensive secondary losses. Taken together, these results suggest that the fungal ENB1 and bacterial ent genes were cooperatively integrated into a functional unit within the W/S clade that enabled adaptation to iron-limited environments. This integrated fungal-bacterial circuit and its dynamic evolution determines the extant distribution of yeast enterobactin producers and cheaters.

Ploidy evolution in a wild yeast is linked to an interaction between cell type and metabolism

Ploidy is an evolutionarily labile trait, and its variation across the tree of life has profound impacts on evolutionary trajectories and life histories. The immediate consequences and molecular causes of ploidy variation on organismal fitness are frequently less clear, although extreme mating type skews in some fungi hint at links between cell type and adaptive traits. Here, we report an unusual recurrent ploidy reduction in replicate populations of the budding yeast Saccharomyces eubayanus experimentally evolved for improvement of a key metabolic trait, the ability to use maltose as a carbon source. We find that haploids have a substantial, but conditional, fitness advantage in the absence of other genetic variation. Using engineered genotypes that decouple the effects of ploidy and cell type, we show that increased fitness is primarily due to the distinct transcriptional program deployed by haploid-like cell types, with a significant but smaller contribution from absolute ploidy. The link between cell-type specification and the carbon metabolism adaptation can be traced to the noncanonical regulation of a maltose transporter by a haploid-specific gene. This study provides novel mechanistic insight into the molecular basis of an environment-cell type fitness interaction and illustrates how selection on traits unexpectedly linked to ploidy states or cell types can drive karyotypic evolution in fungi.

Yueomyces silvicola sp. nov., a novel ascomycetous yeast species unable to utilize ammonium, glutamate, and glutamine as sole nitrogen sources

Five yeast strains isolated from tree bark and rotten wood collected in central and southwestern China, together with four Brazilian strains (three from soil and rotting wood collected in an Amazonian rainforest biome and one from Bromeliad collected in Alagoas state) and one Costa Rican strain isolated from a flower beetle, represent a new species closely related with Yueomyces sinensis in Saccharomycetaceae, as revealed by the 26S ribosomal RNA gene D1/D2 domain and the internal transcribed spacer region sequence analysis. The name Yueomyces silvicola sp. nov. is proposed for this new species with the holotype China General Microbiological Culture Collection Center 2.6469 (= Japan Collection of Microorganisms 34885). The new species exhibits a whole-genome average nucleotide identity value of 77.8% with Y. sinensis. The two Yueomyces species shared unique physiological characteristics of being unable to utilize ammonium and the majority of the amino acids, including glutamate and glutamine, as sole nitrogen sources. Among the 20 amino acids tested, only leucine and tyrosine can be utilized by the Yueomyces species. Genome sequence comparison showed that GAT1, which encodes a GATA family protein participating in transcriptional activation of nitrogen-catabolic genes in Saccharomyces cerevisiae, is absent in the Yueomyces species. However, the failure of the Yueomyces species to utilize ammonium, glutamate, and glutamine, which are generally preferred nitrogen sources for microorganisms, implies that more complicated alterations in the central nitrogen metabolism pathway might occur in the genus Yueomyces.

Danielozyma pruni sp. nov., an asexual yeast species isolated from insect frass

Two strains (NYNU 218101 and NYNU 218104) of an asexual yeast species were isolated from insect frass collected in insect tunnels of red leaf plum trees in the Henan Province, central China. Molecular phylogenetic analysis of the D1/D2 domain of the large subunit rRNA gene and the internal transcribed spacer (ITS) region showed that these two strains belonged to the genus Danielozyma, with Danielozyma litseae as the closest known species. They differed from the type strain of D. litseae by 0.6 % substitutions (three substitutions and one gap) in the D1/D2 domain and by 5.1 % substitutions (19 substitutions and six gaps) in the ITS region, respectively. When compared with the partial ACT1, TEF1 and RPB1 gene sequences, they differed by 3 % (26 substitutions), 2.7 % (25 substitutions) and 9 %(54 substitutions) from D. litseae NRRL YB-3246T in these regions. Physiologically, they also differed from its closest known species D. litseae based on the ability to assimilate inulin and galactitol, as well as to grow in 0.1 % cycloheximide and its inability to ferment maltose and raffinose. In order to classify the two new isolates based on morphological and molecular evidence, we proposed the description of a novel species Danielozyma pruni sp. nov. with strain JCM 35735T as holotype (Mycobank MB 849101).

Biological Factors and Statistical Limitations Prevent Detection of Most Noncanonical Proteins by Mass Spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here we leveraged recent advances in ribosome profiling and mass spectrometry to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly-expressed to be detected by shotgun mass spectrometry at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for four noncanonical proteins in mass spectrometry data, which were also supported by evolution and translation data. These results illustrate the power of mass spectrometry to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly-expressed proteins.

Application of adaptive laboratory evolution for Yarrowia lipolytica: A comprehensive review

Adaptive laboratory evolution is an innovative approach utilized by researchers to enhance the characteristics of microorganisms in the field of biology. With the advancement of this technology, it is now being extended to non-model strains. Yarrowia lipolytica, an oleaginous yeast with significant industrial potential, stands out among the non-conventional fungi. However, the activity of Yarrowia lipolytica is frequently affected by specific substances and environmental factors, necessitating the development of techniques to address these challenges. This manuscript provides an overview of adaptive laboratory evolution experiments conducted on Yarrowia lipolytica, and categorizes the contents into two aspects including improving lignocellulose utilization and enhancing the production in Yarrowia lipolytica. Additionally, we selected several representative examples to illustrate how adaptive laboratory evolution can be combined with other techniques to elucidate the potential mechanisms underlying strain evolution. Lastly, we anticipate a promising future for adaptive laboratory evolution technology and Yarrowia lipolytica in tandem.

Yeast metagenomics: analytical challenges in the analysis of the eukaryotic microbiome

Even if their impact is often underestimated, yeasts and yeast-like fungi represent the most prevalent eukaryotic members of microbial communities on Earth. They play numerous roles in natural ecosystems and in association with their hosts. They are involved in the food industry and pharmaceutical production, but they can also cause diseases in other organisms, making the understanding of their biology mandatory. The ongoing loss of biodiversity due to overexploitation of environmental resources is a growing concern in many countries. Therefore, it becomes crucial to understand the ecology and evolutionary history of these organisms to systematically classify them. To achieve this, it is essential that our knowledge of the mycobiota reaches a level similar to that of the bacterial communities. To overcome the existing challenges in the study of fungal communities, the first step should be the establishment of standardized techniques for the correct identification of species, even from complex matrices, both in wet lab practices and in bioinformatic tools.

Fungi under Modified Atmosphere-The Effects of CO2 Stress on Cell Membranes and Description of New Yeast Stenotrophomyces fumitolerans gen. nov., sp. nov

High levels of carbon dioxide are known to inhibit the growth of microorganisms. A total of twenty strains of filamentous fungi and yeasts were isolated from habitats with enriched carbon dioxide concentration. Most strains were derived from modified atmosphere packed (MAP) food products or mofettes and were cultivated under an atmosphere of 20% CO2 and 80% O2. The influence of CO2 on fungal cell membrane fatty acid profiles was examined in this study. Major changes were the increase in linolenic acid (C18:3 cis 9, 12, 15) and, additionally in most strains, linoleic acid (C18:2 cis 9, 12) with a maximum of 24.8%, at the expense of oleic (C18:1 cis 9), palmitic (C16:0), palmitoleic (C16:1 cis 9) and stearic acid (C18:0). The degree of fatty acid unsaturation increased for all of the strains in the study, which consequently led to lower melting temperatures of the cell membranes after incubation with elevated levels of CO2, indicating fluidization of the membrane and a potential membrane malfunction. Growth was reduced in 18 out of 20 strains in laboratory experiments and a change in pigmentation was observed in several strains. Two of the isolated strains, strain WT5 and strain WR1, were found to represent a hitherto undescribed yeast for which the new genus and species Stenotrophomyces fumitolerans (MB# 849906) is proposed.

Genome sequencing of 15 acid-tolerant yeasts

We report draft genome sequences for 15 non-conventional Saccharomycotina yeast strains obtained from public culture repositories. Included in our collection are eight strains of Pichia with broad tolerance to dicarboxylic acids. The genome sequences of these strains will enable comparative genomics of acid-tolerant phenotypes and strain engineering of non-conventional hosts.

Isolation and sequencing of the spoilage yeast Kregervanrija delftensis

New sequencing workflows can enable the genomics of microbes isolated from craft beverages. Here, we use hybrid, short, and long-read sequencing to assemble the genome of a yeast isolated from cider, Kregervanrija delftensis NCC-J. The Kregervanrija genus has few genomes available; thus, this contributes to yeast genomics.

Extending the reach of homology by using successive computational filters to find yeast pheromone genes

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Parallel nonfunctionalization of CK1δ/ε kinase ohnologs following a whole-genome duplication event

Whole genome duplication (WGD) followed by speciation allows us to examine the parallel evolution of ohnolog pairs. In the yeast family Saccharomycetaceae, HRR25 is a rare case of repeated ohnolog maintenance. This gene has reverted to a single copy in S. cerevisiae where it is now essential, but has been maintained as pairs in at least 7 species post WGD. In S. cerevisiae, HRR25 encodes the casein kinase (CK) 1δ/ε and plays a role in a variety of functions through its kinase activity and protein-protein interactions (PPIs). We hypothesized that the maintenance of duplicated HRR25 ohnologs could be a result of repeated subfunctionalization. We tested this hypothesis through a functional complementation assay in S. cerevisiae, testing all pairwise combinations of 25 orthologs (including 7 ohnolog pairs). Contrary to our expectations, we observed no cases of pair-dependent complementation, which would have supported the subfunctionalization hypothesis. Instead, most post-WGD species have one ohnolog that failed to complement, suggesting their nonfunctionalization or neofunctionalization. The ohnologs incapable of complementation have undergone more rapid protein evolution, lost most PPIs that were observed for their functional counterparts and singletons from post and non-WGD species, and have non-conserved cellular localization, consistent with their ongoing loss of function. The analysis in N. castelli shows that the non-complementing ohnolog is expressed at a lower level and has become non-essential. Taken together, our results indicate that HRR25 orthologs are undergoing gradual nonfunctionalization.

Divergence of TORC1-mediated stress response leads to novel acquired stress resistance in a pathogenic yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

Saccharomycopsis praedatoria sp. nov., a predacious yeast isolated from soil and rotten wood in an Amazonian rainforest biome

Three yeast isolates were obtained from soil and rotting wood samples collected in an Amazonian rainforest biome in Brazil. Comparison of the intergenic spacer 5.8S region and the D1/D2 domains of the large subunit rRNA gene showed that the isolates represent a novel species of the genus Saccharomycopsis. A tree inferred from the D1/D2 sequences placed the novel species near a subclade containing Saccharomycopsis lassenensis, Saccharomycopsis fermentans, Saccharomycopsis javanensis, Saccharomycopsis babjevae, Saccharomycopsis schoenii and Saccharomycopsis oosterbeekiorum, but with low bootstrap support. In terms of sequence divergence, the novel species had the highest identity in the D1/D2 domains with Saccharomycopsis capsularis, from which it differed by 36 substitutions. In contrast, a phylogenomic analysis based on 1061 single-copy orthologs for a smaller set of Saccharomycopsis species whose whole genome sequences are available indicated that the novel species represented by strain UFMG-CM-Y6991 is phylogenetically closer to Saccharomycopsis fodiens and Saccharomycopsis sp. TF2021a (=Saccharomycopsis phalluae). The novel yeast is homothallic and produces asci with one spheroidal ascospore with an equatorial or subequatorial ledge. The name Saccharomycopsis praedatoria sp. nov. is proposed to accommodate the novel species. The holotype of Saccharomycopsis praedatoria is CBS 16589T. The MycoBank number is MB849369. S. praedatoria was able to kill cells of Saccharomyces cerevisiae by means of penetration with infection pegs, a trait common to most species of Saccharomycopsis.

Divergence of TORC1-mediated Stress Response Leads to Novel Acquired Stress Resistance in a Pathogenic Yeast

Acquired stress resistance (ASR) enables organisms to prepare for environmental changes that occur after an initial stressor. However, the genetic basis for ASR and how the underlying network evolved remain poorly understood. In this study, we discovered that a short phosphate starvation induces oxidative stress response (OSR) genes in the pathogenic yeast C. glabrata and protects it against a severe H2O2 stress; the same treatment, however, provides little benefit in the low pathogenic-potential relative, S. cerevisiae. This ASR involves the same transcription factors (TFs) as the OSR, but with different combinatorial logics. We show that Target-of-Rapamycin Complex 1 (TORC1) is differentially inhibited by phosphate starvation in the two species and contributes to the ASR via its proximal effector, Sch9. Therefore, evolution of the phosphate starvation-induced ASR involves the rewiring of TORC1's response to phosphate limitation and the repurposing of TF-target gene networks for the OSR using new regulatory logics.

Diverse signatures of convergent evolution in cacti-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently ~17 times. Using machine-learning, we further found that cactophily can be predicted with 76% accuracy from functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which is likely associated with duplication and altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved through disparate molecular mechanisms. Remarkably, multiple cactophilic lineages and their close relatives are emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-may preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

High quality de novo genome assembly of the non-conventional yeast Kazachstania bulderi describes a potential low pH production host for biorefineries

Kazachstania bulderi is a non-conventional yeast species able to grow efficiently on glucose and δ-gluconolactone at low pH. These unique traits make K. bulderi an ideal candidate for use in sustainable biotechnology processes including low pH fermentations and the production of green chemicals including organic acids. To accelerate strain development with this species, detailed information of its genetics is needed. Here, by employing long read sequencing we report a high-quality phased genome assembly for three strains of K. bulderi species, including the type strain. The sequences were assembled into 12 chromosomes with a total length of 14 Mb, and the genome was fully annotated at structural and functional levels, including allelic and structural variants, ribosomal array and mating type locus. This high-quality reference genome provides a resource to advance our fundamental knowledge of biotechnologically relevant non-conventional yeasts and to support the development of genetic tools for manipulating such strains towards their use as production hosts in biotechnological processes.