Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

Conserved signaling modules regulate filamentous growth in fungi: a model for eukaryotic cell differentiation

Eukaryotic organisms are composed of different cell types with defined shapes and functions. Specific cell types are produced by the process of cell differentiation, which is regulated by signal transduction pathways. Signaling pathways regulate cell differentiation by sensing cues and controlling the expression of target genes whose products generate cell types with specific attributes. In studying how cells differentiate, fungi have proved valuable models because of their ease of genetic manipulation and striking cell morphologies. Many fungal species undergo filamentous growth-a specialized growth pattern where cells produce elongated tube-like projections. Filamentous growth promotes expansion into new environments, including invasion into plant and animal hosts by fungal pathogens. The same signaling pathways that regulate filamentous growth in fungi also control cell differentiation throughout eukaryotes and include highly conserved mitogen-activated protein kinase (MAPK) pathways, which is the focus of this review. In many fungal species, mucin-type sensors regulate MAPK pathways to control filamentous growth in response to diverse stimuli. Once activated, MAPK pathways reorganize cell polarity, induce changes in cell adhesion, and promote the secretion of degradative enzymes that mediate access to new environments. However, MAPK pathway regulation is complicated because related pathways can share components with each other yet induce unique responses (i.e. signal specificity). In addition, MAPK pathways function in highly integrated networks with other regulatory pathways (i.e. signal integration). Here, we discuss signal specificity and integration in several yeast models (mainly Saccharomyces cerevisiae and Candida albicans) by focusing on the filamentation MAPK pathway. Because of the strong evolutionary ties between species, a deeper understanding of the regulation of filamentous growth in established models and increasingly diverse fungal species can reveal fundamentally new mechanisms underlying eukaryotic cell differentiation.

A unique metabolic gene cluster regulates lactose and galactose metabolism in the yeast Candida intermedia

Lactose assimilation is a relatively rare trait in yeasts, and Kluyveromyces yeast species have long served as model organisms for studying lactose metabolism. Meanwhile, the metabolic strategies of most other lactose-assimilating yeasts remain unknown. In this work, we have elucidated the genetic determinants of the superior lactose-growing yeast Candida intermedia. Through genomic and transcriptomic analyses, we identified three interdependent gene clusters responsible for the metabolism of lactose and its hydrolysis product galactose: the conserved LAC cluster (LAC12, LAC4) for lactose uptake and hydrolysis, the conserved GAL cluster (GAL1, GAL7, and GAL10) for galactose catabolism through the Leloir pathway, and a "GALLAC" cluster containing the transcriptional activator gene LAC9, second copies of GAL1 and GAL10, and a XYL1 gene encoding an aldose reductase involved in carbon overflow metabolism. Bioinformatic analysis suggests that the GALLAC cluster is unique to C. intermedia and has evolved through gene duplication and divergence, and deletion mutant phenotyping proved that the cluster is indispensable for C. intermedia's growth on lactose and galactose. We also show that the regulatory network in C. intermedia, governed by Lac9 and Gal1 from the GALLAC cluster, differs significantly from the galactose and lactose regulons in Saccharomyces cerevisiae, Kluyveromyces lactis, and Candida albicans. Moreover, although lactose and galactose metabolism are closely linked in C. intermedia, our results also point to important regulatory differences.IMPORTANCEThis study paves the way to a better understanding of lactose and galactose metabolism in the non-conventional yeast C. intermedia. Notably, the unique GALLAC cluster represents a new, interesting example of metabolic network rewiring and likely helps to explain how C. intermedia has evolved into an efficient lactose-assimilating yeast. With the Leloir pathway of budding yeasts acting like a model system for understanding the function, evolution, and regulation of eukaryotic metabolism, this work provides new evolutionary insights into yeast metabolic pathways and regulatory networks. In extension, the results will facilitate future development and use of C. intermedia as a cell-factory for conversion of lactose-rich whey into value-added products.

Phylogeny and evolution of hemipteran insects based on expanded genomic and transcriptomic data

BACKGROUND

Hemiptera is the fifth species-rich order of insects and the most species-rich order of hemimetabolous insects, including numerous insect species that are of agricultural or medical significance. Despite much effort and recent advance in inferring the Hemiptera phylogeny, some high-level relationships among superfamilies remain controversial.

RESULTS

We sequenced the genomes of 64 hemipteran species from 15 superfamilies and the transcriptomes of two additional scale insect species, integrating them with existing genomic and transcriptomic data to conduct a comprehensive phylogenetic analysis of Hemiptera. Our datasets comprise an average of 1625 nuclear loci of 315 species across 27 superfamilies of Hemiptera. Our analyses supported Cicadoidea and Cercopoidea as sister groups, with Membracoidea typically positioned as the sister to Cicadoidea + Cercopoidea. In most analyses, Aleyrodoidea was recovered as the sister group of all other Sternorrhyncha. A sister-group relationship was supported between Coccoidea and Aphidoidea + Phylloxeroidea. These relationships were further supported by four-cluster likelihood mapping analyses across diverse datasets. Our ancestral state reconstruction indicates phytophagy as the primary feeding strategy for Hemiptera as a whole. However, predation likely represents an ancestral state for Heteroptera, with several phytophagous lineages having evolved from predatory ancestors. Certain lineages, like Lygaeoidea, have undergone a reversal transition from phytophagy to predation. Our divergence time estimation placed the diversification of hemipterans to be between 60 and 150 million years ago.

CONCLUSIONS

By expanding phylogenomic taxon sampling, we clarified the superfamily relationships within the infraorder Cicadomorpha. Our phylogenetic analyses supported the sister-group relationship between the superfamilies Cicadoidea and Cercopoidea, and the superfamily Membracoidea as the sister to Cicadoidea + Cercopoidea. Our divergence time estimation supported the close association of hemipteran diversification with the evolutionary success and adaptive radiation of angiosperms during the Cretaceous period.

Genetic conflicts in budding yeast: The 2μ plasmid as a model selfish element

Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2μ) plasmid as a model selfish element. The 2μ is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2μ plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

Mitochondrial ATP generation is more proteome efficient than glycolysis

Metabolic efficiency profoundly influences organismal fitness. Nonphotosynthetic organisms, from yeast to mammals, derive usable energy primarily through glycolysis and respiration. Although respiration is more energy efficient, some cells favor glycolysis even when oxygen is available (aerobic glycolysis, Warburg effect). A leading explanation is that glycolysis is more efficient in terms of ATP production per unit mass of protein (that is, faster). Through quantitative flux analysis and proteomics, we find, however, that mitochondrial respiration is actually more proteome efficient than aerobic glycolysis. This is shown across yeast strains, T cells, cancer cells, and tissues and tumors in vivo. Instead of aerobic glycolysis being valuable for fast ATP production, it correlates with high glycolytic protein expression, which promotes hypoxic growth. Aerobic glycolytic yeasts do not excel at aerobic growth but outgrow respiratory cells during oxygen limitation. We accordingly propose that aerobic glycolysis emerges from cells maintaining a proteome conducive to both aerobic and hypoxic growth.

Candida boidinii isolates from olive curation water: a promising platform for methanol-based biomanufacturing

Methanol is a promising feedstock for biomanufacturing, but the efficiency of methanol-based bioprocesses is limited by the low rate of methanol utilization pathways and methanol toxicity. Yeast diversity is an attractive biological resource to develop efficient bioprocesses since any effort with strain improvement is more deserving if applied to innate robust strains with relevant catabolic and biosynthetic potential. The present study is in line with such rational and describes the isolation and molecular identification of seven isolates of the methylotrophic species Candida boidinii from waters derived from the traditional curation of olives, in different years, and from contaminated superficial soil near fuel stations. The yeast microbiota from those habitats was also characterized. The four C. boidinii isolates obtained from the curation of olives' water exhibited significantly higher maximum specific growth rates (range 0.15-0.19 h-1), compared with the three isolates obtained from the fuel contaminated soils (range 0.05-0.06 h-1) when grown on methanol as the sole C-source (1% (v/v), in shake flasks, at 30°C). The isolates exhibit significant robustness towards methanol toxicity that increases as the cultivation temperature decreases from 30°C to 25°C. The better methanol-based growth performance exhibited by C. boidinii isolates from olives´ soaking waters could not be essentially attributed to higher methanol tolerance. These methanol-efficient catabolizing isolates are proposed as a promising platform to develop methanol-based bioprocesses.

Post-transfer adaptation of HGT-acquired genes and contribution to guanine metabolic diversification in land plants

Horizontal gene transfer (HGT) is a major driving force in the evolution of prokaryotic and eukaryotic genomes. Despite recent advances in distribution and ecological importance, the extensive pattern, especially in seed plants, and post-transfer adaptation of HGT-acquired genes in land plants remain elusive. We systematically identified 1150 foreign genes in 522 land plant genomes that were likely acquired via at least 322 distinct transfers from nonplant donors and confirmed that recent HGT events were unevenly distributed between seedless and seed plants. HGT-acquired genes evolved to be more similar to native genes in terms of average intron length due to intron gains, and HGT-acquired genes containing introns exhibited higher expression levels than those lacking introns, suggesting that intron gains may be involved in the post-transfer adaptation of HGT in land plants. Functional validation of bacteria-derived gene GuaD in mosses and gymnosperms revealed that the invasion of foreign genes introduced a novel bypass of guanine degradation and resulted in the loss of native pathway genes in some gymnosperms, eventually shaping three major types of guanine metabolism in land plants. We conclude that HGT has played a critical role in land plant evolution.

Ancestral Sequence Reconstruction as a Tool to Detect and Study De Novo Gene Emergence

New protein-coding genes can evolve from previously noncoding genomic regions through a process known as de novo gene emergence. Evidence suggests that this process has likely occurred throughout evolution and across the tree of life. Yet, confidently identifying de novo emerged genes remains challenging. Ancestral sequence reconstruction is a promising approach for inferring whether a gene has emerged de novo or not, as it allows us to inspect whether a given genomic locus ancestrally harbored protein-coding capacity. However, the use of ancestral sequence reconstruction in the context of de novo emergence is still in its infancy and its capabilities, limitations, and overall potential are largely unknown. Notably, it is difficult to formally evaluate the protein-coding capacity of ancestral sequences, particularly when new gene candidates are short. How well-suited is ancestral sequence reconstruction as a tool for the detection and study of de novo genes? Here, we address this question by designing an ancestral sequence reconstruction workflow incorporating different tools and sets of parameters and by introducing a formal criterion that allows to estimate, within a desired level of confidence, when protein-coding capacity originated at a particular locus. Applying this workflow on ∼2,600 short, annotated budding yeast genes (<1,000 nucleotides), we found that ancestral sequence reconstruction robustly predicts an ancient origin for the most widely conserved genes, which constitute "easy" cases. For less robust cases, we calculated a randomization-based empirical P-value estimating whether the observed conservation between the extant and ancestral reading frame could be attributed to chance. This formal criterion allowed us to pinpoint a branch of origin for most of the less robust cases, identifying 49 genes that can unequivocally be considered de novo originated since the split of the Saccharomyces genus, including 37 Saccharomyces cerevisiae-specific genes. We find that for the remaining equivocal cases we cannot rule out different evolutionary scenarios including rapid evolution, multiple gene losses, or a recent de novo origin. Overall, our findings suggest that ancestral sequence reconstruction is a valuable tool to study de novo gene emergence but should be applied with caution and awareness of its limitations.

Reduction of Ribosomal Expansion Segments in Yeast Species of the Magnusiomyces/Saprochaete Clade

Ribosomes are ribonucleoprotein complexes highly conserved across all domains of life. The size differences of ribosomal RNAs (rRNAs) can be mainly attributed to variable regions termed expansion segments (ESs) protruding out from the ribosomal surface. The ESs were found to be involved in a range of processes including ribosome biogenesis and maturation, translation, and co-translational protein modification. Here, we analyze the rRNAs of the yeasts from the Magnusiomyces/Saprochaete clade belonging to the basal lineages of the subphylum Saccharomycotina. We find that these yeasts are missing more than 400 nt from the 25S rRNA and 150 nt from the 18S rRNAs when compared to their canonical counterparts in Saccharomyces cerevisiae. The missing regions mostly map to ESs, thus representing a shift toward a minimal rRNA structure. Despite the structural changes in rRNAs, we did not identify dramatic alterations in the ribosomal protein inventories. We also show that the size-reduced rRNAs are not limited to the species of the Magnusiomyces/Saprochaete clade, indicating that the shortening of ESs happened independently in several other lineages of the subphylum Saccharomycotina.

Phylogenomic analysis of the Candida auris-Candida haemuli clade and related taxa in the Metschnikowiaceae, and proposal of thirteen new genera, fifty-five new combinations and nine new species

Candida is a polyphyletic genus of asexually reproducing yeasts in the Saccharomycotina with more than 400 species that occur in almost all families of the subclass and its name is strongly connected with the infectious disease candidiasis. During the last two decades, approximately half of the Candida species have been reassigned into more than 36 already existing genera and 14 newly proposed genera, but the polyphyletic feature of the genus largely remained. Candida auris is an important, globally emerging opportunistic pathogen that has caused life-threatening outbreaks in healthcare facilities worldwide. This species belongs to the Candida auris-Candida haemuli (CAH) clade in the Metschnikowiaceae, a clade that contains multidrug-resistant clinically relevant species, but also species isolated from natural environments. The clade is phylogenetically positioned remotely from the type species of the genus Candida that is Candida vulgaris (currently interpreted as a synonym of Candida tropicalis) and belongs to the family Debaryomycetaceae. Although previous phylogenetic and phylogenomic studies confirmed the position of C. auris in the Metschnikowiaceae, these analyses failed to resolve the position of the CAH clade within the family and its delimitation from the genera Clavispora and Metschnikowia. To resolve the position of the CAH clade, phylogenomic and comparative genomics analyses were carried out to address the phylogenetic position of C. auris and related species in the Metschnikowiaceae using several metrics, such as the average amino acid identity (AAI) values, the percentage of conserved proteins (POCP) and the presence-absence patterns of orthologs (PAPO). Based on those approaches, 13 new genera are proposed for various Candida and Hyphopichia species, including members of the CAH clade in the Metschnikowiaceae. As a result, C. auris and related species are reassigned to the genus Candidozyma. Fifty-five new combinations and nine new species are introduced and this will reduce the polyphyly of the genus Candida. Citation: Liu F, Hu Z-D, Zhao X-M, et al. 2024. Phylogenomic analysis of the Candida auris-Candida haemuli clade and related taxa in the Metschnikowiaceae, and proposal of thirteen new genera, fifty-five new combinations and nine new species. Persoonia 52: 22-43. https://doi.org/10.3767/persoonia.2024.52.02 .

Saccharomycetaceae: delineation of fungal genera based on phylogenomic analyses, genomic relatedness indices and genomics-based synapomorphies

A correct classification of fungi, including yeasts, is of prime importance to understand fungal biodiversity and to communicate about this diversity. Fungal genera are mainly defined based on phenotypic characteristics and the results of single or multigene-based phylogenetic analyses. However, because yeasts often have less phenotypic characters, their classification experienced a strong move towards DNA-based data, from short ribosomal sequences to multigene phylogenies and more recently to phylogenomics. Here, we explore the usefulness of various genomics-based parameters to circumscribe fungal genera more correctly taking the yeast domain as an example. Therefore, we compared the results of a phylogenomic analysis, average amino acid identity (AAI) values, the presence of conserved signature indels (CSIs), the percentage of conserved proteins (POCP) and the presence-absence patterns of orthologs (PAPO). These genome-based metrics were used to investigate their usefulness in demarcating 13 hitherto relatively well accepted genera in Saccharomycetaceae, namely Eremothecium, Grigorovia, Kazachstania, Kluyveromyces, Lachancea, Nakaseomyces, Naumovozyma, Saccharomyces, Tetrapisispora, Torulaspora, Vanderwaltozyma, Zygosaccharomyces and Zygotorulaspora. As a result, most of these genera are supported by the genomics-based metrics, but the genera Kazachstania, Nakaseomyces and Tetrapisispora were shown to be genetically highly diverse based on the above listed analyses. Considering the results obtained for the presently recognized genera, a range of 80-92 % POCP values and a range of 60-70 % AAI values might be valuable thresholds to discriminate genera in Saccharomycetaceae. Furthermore, the genus-specific genes identified in the PAPO analysis and the CSIs were found to be useful as synapomorphies to characterize and define genera in Saccharomycetaceae. Our results indicate that the combined monophyly-based phylogenomic analysis together with genomic relatedness indices and synapomorphies provide promising approaches to delineating yeast genera and likely those of filamentous fungi as well. The genera Kazachstania, Nakaseomyces and Tetrapisispora are revised and we propose eight new genera and 41 new combinations. Citation: Liu F, Hu Z-D, Yurkov A, et al. 2024. Saccharomycetaceae: delinaeation of fungal genera based on phylogenomic analyses, genomic relatedness indices and genomics-based synapomorphies. Persoonia 52: 1-21. https://doi.org/10.3767/persoonia.2024.52.01.

Specialization restricts the evolutionary paths available to yeast sugar transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning approximately 400 million years of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

Convergent expansions of keystone gene families drive metabolic innovation in a major eukaryotic clade

Many remarkable innovations have repeatedly occurred across vast evolutionary distances. When convergent traits emerge on the tree of life, they are sometimes driven by the same underlying gene families, while other times many different gene families are involved. Conversely, a gene family may be repeatedly recruited for a single trait or many different traits. To understand the general rules governing convergence at both genomic and phenotypic levels, we systematically tested associations between 56 binary metabolic traits and gene count in 14,710 gene families from 993 species of Saccharomycotina yeasts. Using a recently developed phylogenetic approach that reduces spurious correlations, we discovered that gene family expansion and contraction was significantly linked to trait gain and loss in 45/56 (80%) of traits. While 601/746 (81%) of significant gene families were associated with only one trait, we also identified several 'keystone' gene families that were significantly associated with up to 13/56 (23%) of all traits. These results indicate that metabolic innovations in yeasts are governed by a narrow set of major genetic elements and mechanisms.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Massively integrated coexpression analysis reveals transcriptional regulation, evolution and cellular implications of the yeast noncanonical translatome

BACKGROUND

Recent studies uncovered pervasive transcription and translation of thousands of noncanonical open reading frames (nORFs) outside of annotated genes. The contribution of nORFs to cellular phenotypes is difficult to infer using conventional approaches because nORFs tend to be short, of recent de novo origins, and lowly expressed. Here we develop a dedicated coexpression analysis framework that accounts for low expression to investigate the transcriptional regulation, evolution, and potential cellular roles of nORFs in Saccharomyces cerevisiae.

RESULTS

Our results reveal that nORFs tend to be preferentially coexpressed with genes involved in cellular transport or homeostasis but rarely with genes involved in RNA processing. Mechanistically, we discover that young de novo nORFs located downstream of conserved genes tend to leverage their neighbors' promoters through transcription readthrough, resulting in high coexpression and high expression levels. Transcriptional piggybacking also influences the coexpression profiles of young de novo nORFs located upstream of genes, but to a lesser extent and without detectable impact on expression levels. Transcriptional piggybacking influences, but does not determine, the transcription profiles of de novo nORFs emerging nearby genes. About 40% of nORFs are not strongly coexpressed with any gene but are transcriptionally regulated nonetheless and tend to form entirely new transcription modules. We offer a web browser interface ( https://carvunislab.csb.pitt.edu/shiny/coexpression/ ) to efficiently query, visualize, and download our coexpression inferences.

CONCLUSIONS

Our results suggest that nORF transcription is highly regulated. Our coexpression dataset serves as an unprecedented resource for unraveling how nORFs integrate into cellular networks, contribute to cellular phenotypes, and evolve.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

Candida auris from the Egyptian cobra: Role of snakes as potential reservoirs

Candida auris represents one of the most urgent threats to public health, although its ecology remains largely unknown. Because amphibians and reptiles may present favorable conditions for C. auris colonization, cloacal and blood samples (n = 68), from several snake species, were cultured and molecularly screened for C. auris using molecular amplification of glycosylphosphatidylinositol protein-encoding genes and ribosomal internal transcribed spacer sequencing. Candida auris was isolated from the cloacal swab of one Egyptian cobra (Naja haje legionis) and molecularly identified in its cloaca and blood. The isolation of C. auris from wild animals is herein reported for the first time, thus suggesting the role that these animals could play as reservoirs of this emerging pathogen. The occurrence of C. auris in blood requires further investigation, although the presence of cationic antimicrobial peptides in the plasma of reptiles could play a role in reducing the vitality of the fungus.

Discovering Fragile Clades and Causal Sequences in Phylogenomics by Evolutionary Sparse Learning

Phylogenomic analyses of long sequences, consisting of many genes and genomic segments, reconstruct organismal relationships with high statistical confidence. But, inferred relationships can be sensitive to excluding just a few sequences. Currently, there is no direct way to identify fragile relationships and the associated individual gene sequences in species. Here, we introduce novel metrics for gene-species sequence concordance and clade probability derived from evolutionary sparse learning models. We validated these metrics using fungi, plant, and animal phylogenomic datasets, highlighting the ability of the new metrics to pinpoint fragile clades and the sequences responsible. The new approach does not necessitate the investigation of alternative phylogenetic hypotheses, substitution models, or repeated data subset analyses. Our methodology offers a streamlined approach to evaluating major inferred clades and identifying sequences that may distort reconstructed phylogenies using large datasets.

Convergent reductive evolution in bee-associated lactic acid bacteria

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae (LAB) family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts towards the floral niche in LAB. In these bee-associated LAB, we observed pervasive, significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, adhE, encodes a bifunctional aldehyde-alcohol dehydrogenase associated with the evolution of fructophily, a rare phenotypic trait that was recently identified in many floral LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent loss of the same set of genes.

Importance

Several lactic acid bacteria (LAB) species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine-learning based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference of fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

Glycosaminoglycan lyase: A new competition between bacteria and the pacific white shrimp Litopenaeus vannamei

Horizontal gene transfer (HGT) is an important evolutionary force in the formation of prokaryotic and eukaryotic genomes. In recent years, many HGT genes horizontally transferred from prokaryotes to eukaryotes have been reported, and most of them are present in arthropods. The Pacific white shrimp Litopenaeus vannamei, an important economic species of arthropod, has close relationships with bacteria, providing a platform for horizontal gene transfer (HGT). In this study, we analyzed bacteria-derived HGT based on a high-quality genome of L. vannamei via a homology search and phylogenetic analysis, and six HGT genes were identified. Among these six horizontally transferred genes, we found one gene (LOC113799989) that contains a bacterial chondroitinase AC structural domain and encodes an unknown glycosaminoglycan (GAG) lyase in L. vannamei. The real-time quantitative PCR results showed that the mRNA expression level of LOC113799989 was highest in the hepatopancreas and heart, and after stimulation by Vibrio parahaemolyticus, its mRNA expression level was rapidly up-regulated within 12 h. Furthermore, after injecting si-RNA and stimulation by V. parahaemolyticus, we found that the experimental group had a higher cumulative mortality rate in 48 h than the control group, indicating that the bacteria-derived GAG lyase can reduce the mortality of shrimp with respect to infection by V. parahaemolyticus and might be related to the resistance of shrimp to bacterial diseases. Our findings contribute to the study of the function of GAGs and provide new insights into GAG-related microbial pathogenesis and host defense mechanisms in arthropods.

Candida glabrata maintains two Hap1 homologs, Zcf27 and Zcf4, for distinct roles in ergosterol gene regulation to mediate sterol homeostasis under azole and hypoxic conditions

Candida glabrata exhibits innate resistance to azole antifungal drugs but also has the propensity to rapidly develop clinical drug resistance. Azole drugs, which target Erg11, is one of the three major classes of antifungals used to treat Candida infections. Despite their widespread use, the mechanism controlling azole-induced ERG gene expression and drug resistance in C. glabrata has primarily revolved around Upc2 and/or Pdr1. In this study, we determined the function of two zinc cluster transcription factors, Zcf27 and Zcf4, as direct but distinct regulators of ERG genes. Our phylogenetic analysis revealed C. glabrata Zcf27 and Zcf4 as the closest homologs to Saccharomyces cerevisiae Hap1. Hap1 is a known zinc cluster transcription factor in S. cerevisiae in controlling ERG gene expression under aerobic and hypoxic conditions. Interestingly, when we deleted HAP1 or ZCF27 in either S. cerevisiae or C. glabrata, respectively, both deletion strains showed altered susceptibility to azole drugs, whereas the strain deleted for ZCF4 did not exhibit azole susceptibility. We also determined that the increased azole susceptibility in a zcf27Δ strain is attributed to decreased azole-induced expression of ERG genes, resulting in decreased levels of total ergosterol. Surprisingly, Zcf4 protein expression is barely detected under aerobic conditions but is specifically induced under hypoxic conditions. However, under hypoxic conditions, Zcf4 but not Zcf27 was directly required for the repression of ERG genes. This study provides the first demonstration that Zcf27 and Zcf4 have evolved to serve distinct roles allowing C. glabrata to adapt to specific host and environmental conditions.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Saccharomycopsis fibuligera (lindner) Klocker 1907 (saccharomycetales: saccharomycopsidaceae)

Saccharomycopsis fibuligera (Lindner) Klocker 1907 is frequently employed in the fermentation of metabolites such as citric acid, ethanol, mannitol, and pyruvate. Its heat tolerance and alcohol-producing capabilities during fermentation make it a desirable option for bread and wine production. To date, the mitochondrial genome of S. fibuligera has not been sequenced. In the present study, we obtained the full mitochondrial genome of S. fibuligera, which is 57,302 bp long and has a GC content of 24.40%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 21 intronic ORFs, 25 tRNAs, and 2 rRNA genes. By utilizing the Bayesian inference phylogenetic method, we constructed phylogenetic trees for 24 Saccharomycotina fungi, which indicated that S. fibuligera is closely related to S. capsularis.

A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1,154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family and genome evolution are distinct from plants, animals, and filamentous ascomycetes and are characterized by small genome sizes and smaller gene numbers but larger gene family sizes. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

Phylogenetic analyses reveal insights into interdomain horizontal gene transfer of microbial lipases

Microbial lipases play a pivotal role in a wide range of biotechnological processes and in the human skin microbiome. However, their evolution remains poorly understood. Accessing the evolutionary process of lipases could contribute to future applications in health and biotechnology. We investigated genetic events associated with the evolutionary trajectory of the microbial family LIP lipases. Using phylogenetic analysis, we identified two distinct horizontal gene transfer (HGT) events from Bacteria to Fungi. Further analysis of human cutaneous mycobiome members such as the lipophilic Malassezia yeasts and CUG-Ser-1 clade (including Candida sp. and other microorganisms associated with cutaneous mycobiota) revealed recent evolutionary processes, with multiple gene duplication events. The Lid region of fungal lipases, crucial for substrate interaction, exhibits varying degrees of conservation among different groups. Our findings suggest the adaptability of the fungal LIP family in various genetic and metabolic contexts and its potential role in niche exploration.

Logistic PCA explains differences between genome-scale metabolic models in terms of metabolic pathways

Genome-scale metabolic models (GSMMs) offer a holistic view of biochemical reaction networks, enabling in-depth analyses of metabolism across species and tissues in multiple conditions. However, comparing GSMMs Against each other poses challenges as current dimensionality reduction algorithms or clustering methods lack mechanistic interpretability, and often rely on subjective assumptions. Here, we propose a new approach utilizing logisitic principal component analysis (LPCA) that efficiently clusters GSMMs while singling out mechanistic differences in terms of reactions and pathways that drive the categorization. We applied LPCA to multiple diverse datasets, including GSMMs of 222 Escherichia-strains, 343 budding yeasts (Saccharomycotina), 80 human tissues, and 2943 Firmicutes strains. Our findings demonstrate LPCA's effectiveness in preserving microbial phylogenetic relationships and discerning human tissue-specific metabolic profiles, exhibiting comparable performance to traditional methods like t-distributed stochastic neighborhood embedding (t-SNE) and Jaccard coefficients. Moreover, the subsystems and associated reactions identified by LPCA align with existing knowledge, underscoring its reliability in dissecting GSMMs and uncovering the underlying drivers of separation.

Chromosome-level genome assembly of the yeast Lodderomyces beijingensis reveals the genetic nature of metabolic adaptations and identifies subtelomeres as hotspots for amplification of mating type loci

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterized its genome sequence.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Extensive remodeling of sugar metabolism through gene loss and horizontal gene transfer in a eukaryotic lineage

BACKGROUND

In yeasts belonging to the subphylum Saccharomycotina, genes encoding components of the main metabolic pathways, like alcoholic fermentation, are usually conserved. However, in fructophilic species belonging to the floral Wickerhamiella and Starmerella genera (W/S clade), alcoholic fermentation was uniquely shaped by events of gene loss and horizontal gene transfer (HGT).

RESULTS

Because HGT and gene losses were first identified when only eight W/S-clade genomes were available, we collected publicly available genome data and sequenced the genomes of 36 additional species. A total of 63 genomes, representing most of the species described in the clade, were included in the analyses. Firstly, we inferred the phylogenomic tree of the clade and inspected the genomes for the presence of HGT-derived genes involved in fructophily and alcoholic fermentation. We predicted nine independent HGT events and several instances of secondary loss pertaining to both pathways. To investigate the possible links between gene loss and acquisition events and evolution of sugar metabolism, we conducted phenotypic characterization of 42 W/S-clade species including estimates of sugar consumption rates and fermentation byproduct formation. In some instances, the reconciliation of genotypes and phenotypes yielded unexpected results, such as the discovery of fructophily in the absence of the cornerstone gene (FFZ1) and robust alcoholic fermentation in the absence of the respective canonical pathway.

CONCLUSIONS

These observations suggest that reinstatement of alcoholic fermentation in the W/S clade triggered a surge of innovation that goes beyond the utilization of xenologous enzymes, with fructose metabolism playing a key role.

The relationship between cell density and cell count differs among Saccharomyces yeast species

There is a recent push to develop wild and non-domesticated Saccharomyces yeast strains into useful model systems for research in ecology and evolution. Yet, the variation between species and strains in important population parameters remains largely undescribed. Here, we investigated the relationship between two commonly used measures in microbiology to estimate growth rate - cell density and cell count - in 23 strains across all eight Saccharomyces species . We found that the slope of this relationship significantly differs among species and a given optical density (OD) does not translate into the same number of cells across species. We provide a cell number calculator based on our OD measurements for each strain used in this study. Surprisingly, we found a slightly positive relationship between cell size and the slope of the cell density-cell count relationship. Our results show that the strain- and species-specificity of the cell density and cell count relationship should be taken into account, for instance when running competition experiments requiring equal starting population sizes or when estimating the fitness of strains with different genetic backgrounds in experimental evolution studies.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Machine learning enables identification of an alternative yeast galactose utilization pathway

How genomic differences contribute to phenotypic differences is a major question in biology. The recently characterized genomes, isolation environments, and qualitative patterns of growth on 122 sources and conditions of 1,154 strains from 1,049 fungal species (nearly all known) in the yeast subphylum Saccharomycotina provide a powerful, yet complex, dataset for addressing this question. We used a random forest algorithm trained on these genomic, metabolic, and environmental data to predict growth on several carbon sources with high accuracy. Known structural genes involved in assimilation of these sources and presence/absence patterns of growth in other sources were important features contributing to prediction accuracy. By further examining growth on galactose, we found that it can be predicted with high accuracy from either genomic (92.2%) or growth data (82.6%) but not from isolation environment data (65.6%). Prediction accuracy was even higher (93.3%) when we combined genomic and growth data. After the GALactose utilization genes, the most important feature for predicting growth on galactose was growth on galactitol, raising the hypothesis that several species in two orders, Serinales and Pichiales (containing the emerging pathogen Candida auris and the genus Ogataea, respectively), have an alternative galactose utilization pathway because they lack the GAL genes. Growth and biochemical assays confirmed that several of these species utilize galactose through an alternative oxidoreductive D-galactose pathway, rather than the canonical GAL pathway. Machine learning approaches are powerful for investigating the evolution of the yeast genotype-phenotype map, and their application will uncover novel biology, even in well-studied traits.

Ancestral sequence reconstruction of Mic60 reveals a residue signature supporting respiration in yeast

In eukaryotes, the essential process of cellular respiration takes place in the cristae of mitochondria. The protein Mic60 is known to stabilize crista junctions; however, how the C-terminal Mitofilin domain of Mic60 mediates cristae-supported respiration remains elusive. Here, we used ancestral sequence reconstruction to generate Mitofilin ancestors up to and including the last opisthokont common ancestor (LOCA). We found that yeast-lineage derived Mitofilin ancestors as far back as the LOCA rescue respiration. By comparing Mitofilin ancestors with different respiratory phenotypes, we identify four residues that explain the difference between respiration functional yeast- and non-functional animal-derived common Mitofilin ancestors. Our results imply that Mitofilin-supported respiration in yeast stems from a conserved mechanism, and provide a foundation for investigating the divergence of candidate crista junction interactions present during the emergence of eukaryotes.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

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 determine the extant distribution of yeast enterobactin producers and cheaters.

Patterns recovered in phylogenomic analysis of Candida auris and close relatives implicate broad environmental flexibility in Candida/Clavispora clade yeasts

Fungal pathogens commonly originate from benign or non-pathogenic strains living in the natural environment. The recently emerged human pathogen, Candida auris, is one example of a fungus believed to have originated in the environment and recently transitioned into a clinical setting. To date, however, there is limited evidence about the origins of this species in the natural environment and when it began associating with humans. One approach to overcome this gap is to reconstruct phylogenetic relationships between (1) strains isolated from clinical and non-clinical environments and (2) between species known to cause disease in humans and benign environmental saprobes. C. auris belongs to the Candida/Clavispora clade, a diverse group of 45 yeast species including human pathogens and environmental saprobes. We present a phylogenomic analysis of the Candida/Clavispora clade aimed at understanding the ecological breadth and evolutionary relationships between an expanded sample of environmentally and clinically isolated yeasts. To build a robust framework for investigating these relationships, we developed a whole-genome sequence dataset of 108 isolates representing 18 species, including four newly sequenced species and 18 environmentally isolated strains. Our phylogeny, based on 619 orthologous genes, shows environmentally isolated species and strains interspersed with clinically isolated counterparts, suggesting that there have been many transitions between humans and the natural environment in this clade. Our findings highlight the breadth of environments these yeasts inhabit and imply that many clinically isolated yeasts in this clade could just as easily live outside the human body in diverse natural environments and vice versa.

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.