Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum

DNA damage response of major fungal pathogen Candida glabrata offers clues to explain its genetic diversity

How cells respond to DNA damage is key to maintaining genome integrity or facilitating genetic change. In fungi, DNA damage responses have been extensively characterized in the model budding yeast Saccharomyces cerevisiae, which is generally not pathogenic. However, it is not clear how closely these responses resemble those in fungal pathogens, in which genetic change plays an important role in the evolutionary arms race between pathogen and host and the evolution of antifungal drug resistance. A close relative of S. cerevisiae, Candida glabrata, is an opportunistic pathogen that displays high variability in chromosome structure among clinical isolates and rapidly evolves antifungal drug resistance. The mechanisms facilitating such genomic flexibility and evolvability in this organism are unknown. Recently we characterized the DNA damage response of C. glabrata and identified several features that distinguish it from the well characterized DNA damage response of S. cerevisiae. First, we discovered that, in contrast to the established paradigm, C. glabrata effector kinase Rad53 is not hyperphosphorylated upon DNA damage. We also uncovered evidence of an attenuated DNA damage checkpoint response, wherein in the presence of DNA damage C. glabrata cells did not accumulate in S-phase and proceeded with cell division, leading to aberrant mitoses and cell death. Finally, we identified evidence of transcriptional rewiring of the DNA damage response of C. glabrata relative to S. cerevisiae, including an upregulation of genes involved in mating and meiosis-processes that have not been reported in C. glabrata. Together, these results open new possibilities and raise tantalizing questions of how this major fungal pathogen facilitates genetic change.

Dissecting incongruence between concatenation- and quartet-based approaches in phylogenomic data

Topological conflict or incongruence is widespread in phylogenomic data. Concatenation- and coalescent-based approaches often result in incongruent topologies, but the causes of this conflict can be difficult to characterize. We examined incongruence stemming from conflict between likelihood-based signal (quantified by the difference in gene-wise log likelihood score or ΔGLS) and quartet-based topological signal (quantified by the difference in gene-wise quartet score or ΔGQS) for every gene in three phylogenomic studies in animals, fungi, and plants, which were chosen because their concatenation-based IQ-TREE (T1) and quartet-based ASTRAL (T2) phylogenies are known to produce eight conflicting internal branches (bipartitions). By comparing the types of phylogenetic signal for all genes in these three data matrices, we found that 30% - 36% of genes in each data matrix are inconsistent, that is, each of these genes has higher log likelihood score for T1 versus T2 (i.e., ΔGLS >0) whereas its T1 topology has lower quartet score than its T2 topology (i.e., ΔGQS <0) or vice versa. Comparison of inconsistent and consistent genes using a variety of metrics (e.g., evolutionary rate, gene tree topology, distribution of branch lengths, hidden paralogy, and gene tree discordance) showed that inconsistent genes are more likely to recover neither T1 nor T2 and have higher levels of gene tree discordance than consistent genes. Simulation analyses demonstrate that removal of inconsistent genes from datasets with low levels of incomplete lineage sorting (ILS) and low and medium levels of gene tree estimation error (GTEE) reduced incongruence and increased accuracy. In contrast, removal of inconsistent genes from datasets with medium and high ILS levels and high GTEE levels eliminated or extensively reduced incongruence, but the resulting congruent species phylogenies were not always topologically identical to the true species trees.

QMaker: Fast and accurate method to estimate empirical models of protein evolution

Amino acid substitution models play a crucial role in phylogenetic analyses. Maximum likelihood (ML) methods have been proposed to estimate amino acid substitution models; however, they are typically complicated and slow. In this article, we propose QMaker, a new ML method to estimate a general time-reversible \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Q$\end{document} matrix from a large protein data set consisting of multiple sequence alignments. QMaker combines an efficient ML tree search algorithm, a model selection for handling the model heterogeneity among alignments, and the consideration of rate mixture models among sites. We provide QMaker as a user-friendly function in the IQ-TREE software package (http://www.iqtree.org) supporting the use of multiple CPU cores so that biologists can easily estimate amino acid substitution models from their own protein alignments. We used QMaker to estimate new empirical general amino acid substitution models from the current Pfam database as well as five clade-specific models for mammals, birds, insects, yeasts, and plants. Our results show that the new models considerably improve the fit between model and data and in some cases influence the inference of phylogenetic tree topologies.[Amino acid replacement matrices; amino acid substitution models; maximum likelihood estimation; phylogenetic inferences.]

Engineering heterologous molybdenum-cofactor-biosynthesis and nitrate-assimilation pathways enables nitrate utilization by Saccharomyces cerevisiae

Metabolic capabilities of cells are not only defined by their repertoire of enzymes and metabolites, but also by availability of enzyme cofactors. The molybdenum cofactor (Moco) is widespread among eukaryotes but absent from the industrial yeast Saccharomyces cerevisiae. No less than 50 Moco-dependent enzymes covering over 30 catalytic activities have been described to date, introduction of a functional Moco synthesis pathway offers interesting options to further broaden the biocatalytic repertoire of S. cerevisiae. In this study, we identified seven Moco biosynthesis genes in the non-conventional yeast Ogataea parapolymorpha by SpyCas9-mediated mutational analysis and expressed them in S. cerevisiae. Functionality of the heterologously expressed Moco biosynthesis pathway in S. cerevisiae was assessed by co-expressing O. parapolymorpha nitrate-assimilation enzymes, including the Moco-dependent nitrate reductase. Following two-weeks of incubation, growth of the engineered S. cerevisiae strain was observed on nitrate as sole nitrogen source. Relative to the rationally engineered strain, the evolved derivatives showed increased copy numbers of the heterologous genes, increased levels of the encoded proteins and a 5-fold higher nitrate-reductase activity in cell extracts. Growth at nM molybdate concentrations was enabled by co-expression of a Chlamydomonas reinhardtii high-affinity molybdate transporter. In serial batch cultures on nitrate-containing medium, a non-engineered S. cerevisiae strain was rapidly outcompeted by the spoilage yeast Brettanomyces bruxellensis. In contrast, an engineered and evolved nitrate-assimilating S. cerevisiae strain persisted during 35 generations of co-cultivation. This result indicates that the ability of engineered strains to use nitrate may be applicable to improve competitiveness of baker's yeast in industrial processes upon contamination with spoilage yeasts.

Tracking gene expression, metabolic profiles, and biochemical analysis in the halotolerant basidiomycetous yeast Rhodotorula mucilaginosa EXF-1630 during benzo[a]pyrene and phenanthrene biodegradation under hypersaline conditions

Polyaromatic phenanthrene (Phe) and benzo[a]pyrene (BaP) are highly toxic, mutagenic, and carcinogenic contaminants widely dispersed in nature, including saline environments. Polyextremotolerant Rhodotorula mucilaginosa EXF-1630, isolated from Arctic sea ice, was grown on a huge concentration range -10 to 500 ppm- of Phe and BaP as sole carbon sources at hypersaline conditions (1 M NaCl). Selected polycyclic aromatic hydrocarbons (PAHs) supported growth as well as glucose, even at high PAH concentrations. Initially, up to 40% of Phe and BaP were adsorbed, followed by biodegradation, resulting in 80% removal in 10 days. While extracellular laccase, peroxidase, and un-specific peroxygenase activities were not detected, NADPH-cytochrome c reductase activity peaked at 4 days. The successful removal of PAHs and the absence of toxic metabolites were confirmed by toxicological tests on moss Physcomitrium patens, bacterium Aliivibrio fischeri, human erythrocytes, and pulmonary epithelial cells (A549). Metabolic profiles were determined at the midpoint of the biodegradation exponential phase, with added Phe and BaP (100 ppm) and 1 M NaCl. Different hydroxylated products were found in the culture medium, while the conjugative metabolite 1-phenanthryl-β-D-glucopyranose was detected in the medium and in the cells. Transcriptome analysis resulted in 870 upregulated and 2,288 downregulated transcripts on PAHs, in comparison to glucose. Genomic mining of 61 available yeast genomes showed a widespread distribution of 31 xenobiotic degradation pathways in different yeast lineages. Two distributions with similar metabolic capacities included black yeasts and mainly members of the Sporidiobolaceae family (including EXF-1630), respectively. This is the first work describing a metabolic profile and transcriptomic analysis of PAH degradation by yeast.

PhyKIT: a broadly applicable UNIX shell toolkit for processing and analyzing phylogenomic data

MOTIVATION

Diverse disciplines in biology process and analyze multiple sequence alignments (MSAs) and phylogenetic trees to evaluate their information content, infer evolutionary events and processes, and predict gene function. However, automated processing of MSAs and trees remains a challenge due to the lack of a unified toolkit. To fill this gap, we introduce PhyKIT, a toolkit for the UNIX shell environment with 30 functions that process MSAs and trees, including but not limited to estimation of mutation rate, evaluation of sequence composition biases, calculation of the degree of violation of a molecular clock, and collapsing bipartitions (internal branches) with low support.

RESULTS

To demonstrate the utility of PhyKIT, we detail three use cases: (1) summarizing information content in MSAs and phylogenetic trees for diagnosing potential biases in sequence or tree data; (2) evaluating gene-gene covariation of evolutionary rates to identify functional relationships, including novel ones, among genes; and (3) identify lack of resolution events or polytomies in phylogenetic trees, which are suggestive of rapid radiation events or lack of data. We anticipate PhyKIT will be useful for processing, examining, and deriving biological meaning from increasingly large phylogenomic datasets.

AVAILABILITY

PhyKIT is freely available on GitHub (https://github.com/JLSteenwyk/PhyKIT), PyPi (https://pypi.org/project/phykit/), and the Anaconda Cloud (https://anaconda.org/JLSteenwyk/phykit) under the MIT license with extensive documentation and user tutorials (https://jlsteenwyk.com/PhyKIT).

SUPPLEMENTARY INFORMATION

Supplementary data are available on figshare (doi: 10.6084/m9.figshare.13118600) and are available at Bioinformatics online.

Systematic Comparison of Cell Wall-Related Proteins of Different Yeasts

Yeast cell walls have two major roles, to preserve physical integrity of the cell, and to ensure communication with surrounding molecules and cells. While the first function requires evolutionary conserved polysaccharide network synthesis, the second needs to be flexible and provide adaptability to different habitats and lifestyles. In this study, the comparative in silico analysis of proteins required for cell wall biosynthesis and functions containing 187 proteins of 92 different yeasts was performed in order to assess which proteins were broadly conserved among yeasts and which were more species specific. Proteins were divided into several groups according to their role and localization. As expected, many Saccharomyces cerevisiae proteins involved in protein glycosylation, glycosylphosphatidylinositol (GPI) synthesis and the synthesis of wall polysaccharides had orthologues in most other yeasts. Similarly, a group of GPI anchored proteins involved in cell wall biosynthesis (Gas proteins and Dfg5p/Dcw1p) and other non-GPI anchored cell wall proteins involved in the wall synthesis and remodeling were highly conserved. However, GPI anchored proteins involved in flocculation, aggregation, cell separation, and those of still unknown functions were not highly conserved. The proteins localized in the cell walls of various yeast species were also analyzed by protein biotinylation and blotting. Pronounced differences were found both in the patterns, as well as in the overall amounts of different groups of proteins. The amount of GPI-anchored proteins correlated with the mannan to glucan ratio of the wall. Changes of the wall proteome upon temperature shift to 42 °C were detected.

Novakomyces olei sp. nov., the First Member of a Novel Taphrinomycotina Lineage

Taphrinomycotina is the smallest subphylum of the phylum Ascomycota. It is an assemblage of distantly related early diverging lineages of the phylum, comprising organisms with divergent morphology and ecology; however, phylogenomic analyses support its monophyly. In this study, we report the isolation of a yeast strain, which could not be assigned to any of the currently recognised five classes of Taphrinomycotina. The strain of the novel budding species was recovered from extra virgin olive oil and characterised phenotypically by standard methods. The ultrastructure of the cell wall was investigated by transmission electron microscopy. Comparisons of barcoding DNA sequences indicated that the investigated strain is not closely related to any known organism. Tentative phylogenetic placement was achieved by maximum-likelihood analysis of the D1/D2 domain of the nuclear LSU rRNA gene. The genome of the investigated strain was sequenced, assembled, and annotated. Phylogenomic analyses placed it next to the fission Schizosaccharomyces species. To accommodate the novel species, Novakomyces olei, a novel genus Novakomyces, a novel family Novakomycetaceae, a novel order Novakomycetales, and a novel class Novakomycetes is proposed as well. Functional analysis of genes missing in N. olei in comparison to Schizosaccharomyces pombe revealed that they are biased towards biosynthesis of complex organic molecules, regulation of mRNA, and the electron transport chain. Correlating the genome content and physiology among species of Taphrinomycotina revealed some discordance between pheno- and genotype. N. olei produced ascospores in axenic culture preceded by conjugation between two cells. We confirmed that N. olei is a primary homothallic species lacking genes for different mating types.

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

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

Ancient balancing selection maintains incompatible versions of the galactose pathway in yeast

Metabolic pathways differ across species but are expected to be similar within a species. We discovered two functional, incompatible versions of the galactose pathway in Saccharomyces cerevisiae We identified a three-locus genetic interaction for growth in galactose, and used precisely engineered alleles to show that it arises from variation in the galactose utilization genes GAL2, GAL1/10/7, and phosphoglucomutase (PGM1), and that the reference allele of PGM1 is incompatible with the alternative alleles of the other genes. Multiloci balancing selection has maintained the two incompatible versions of the pathway for millions of years. Strains with alternative alleles are found primarily in galactose-rich dairy environments, and they grow faster in galactose but slower in glucose, revealing a trade-off on which balancing selection may have acted.

Translational landscape and protein biogenesis demands of the early secretory pathway in Komagataella phaffii

Background

Eukaryotes use distinct networks of biogenesis factors to synthesize, fold, monitor, traffic, and secrete proteins. During heterologous expression, saturation of any of these networks may bottleneck titer and yield. To understand the flux through various routes into the early secretory pathway, we quantified the global and membrane-associated translatomes of Komagataella phaffii.

Results

By coupling Ribo-seq with long-read mRNA sequencing, we generated a new annotation of protein-encoding genes. By using Ribo-seq with subcellular fractionation, we quantified demands on co- and posttranslational translocation pathways. During exponential growth in rich media, protein components of the cell-wall represent the greatest number of nascent chains entering the ER. Transcripts encoding the transmembrane protein PMA1 sequester more ribosomes at the ER membrane than any others. Comparison to Saccharomyces cerevisiae reveals conservation in the resources allocated by gene ontology, but variation in the diversity of gene products entering the secretory pathway.

Conclusion

A subset of host proteins, particularly cell-wall components, impose the greatest biosynthetic demands in the early secretory pathway. These proteins are potential targets in strain engineering aimed at alleviating bottlenecks during heterologous protein production.

Komagataella phaffii as Emerging Model Organism in Fundamental Research

Komagataella phaffii (Pichia pastoris) is one of the most extensively applied yeast species in pharmaceutical and biotechnological industries, and, therefore, also called the biotech yeast. However, thanks to more advanced strain engineering techniques, it recently started to gain attention as model organism in fundamental research. So far, the most studied model yeast is its distant cousin, Saccharomyces cerevisiae. While these data are of great importance, they limit our knowledge to one organism only. Since the divergence of the two species 250 million years ago, K. phaffii appears to have evolved less rapidly than S. cerevisiae, which is why it remains more characteristic of the common ancient yeast ancestors and shares more features with metazoan cells. This makes K. phaffii a valuable model organism for research on eukaryotic molecular cell biology, a potential we are only beginning to fully exploit. As methylotrophic yeast, K. phaffii has the intriguing property of being able to efficiently assimilate methanol as a sole source of carbon and energy. Therefore, major efforts have been made using K. phaffii as model organism to study methanol assimilation, peroxisome biogenesis and pexophagy. Other research topics covered in this review range from yeast genetics including mating and sporulation behavior to other cellular processes such as protein secretion, lipid biosynthesis and cell wall biogenesis. In this review article, we compare data obtained from K. phaffii with S. cerevisiae and other yeasts whenever relevant, elucidate major differences, and, most importantly, highlight the big potential of using K. phaffii in fundamental research.

Expanding the Knowledge on the Skillful Yeast Cyberlindnera jadinii

Cyberlindnera jadinii is widely used as a source of single-cell protein and is known for its ability to synthesize a great variety of valuable compounds for the food and pharmaceutical industries. Its capacity to produce compounds such as food additives, supplements, and organic acids, among other fine chemicals, has turned it into an attractive microorganism in the biotechnology field. In this review, we performed a robust phylogenetic analysis using the core proteome of C. jadinii and other fungal species, from Asco- to Basidiomycota, to elucidate the evolutionary roots of this species. In addition, we report the evolution of this species nomenclature over-time and the existence of a teleomorph (C. jadinii) and anamorph state (Candida utilis) and summarize the current nomenclature of most common strains. Finally, we highlight relevant traits of its physiology, the solute membrane transporters so far characterized, as well as the molecular tools currently available for its genomic manipulation. The emerging applications of this yeast reinforce its potential in the white biotechnology sector. Nonetheless, it is necessary to expand the knowledge on its metabolism, regulatory networks, and transport mechanisms, as well as to develop more robust genetic manipulation systems and synthetic biology tools to promote the full exploitation of C. jadinii.

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides. High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

Giant GAL gene clusters for the melibiose-galactose pathway in Torulaspora

In many yeast species, the three genes at the centre of the galactose catabolism pathway, GAL1, GAL10 and GAL7, are neighbours in the genome and form a metabolic gene cluster. We report here that some yeast strains in the genus Torulaspora have much larger GAL clusters that include genes for melibiase (MEL1), galactose permease (GAL2), glucose transporter (HGT1), phosphoglucomutase (PGM1) and the transcription factor GAL4, in addition to GAL1, GAL10, and GAL7. Together, these eight genes encode almost all the steps in the pathway for catabolism of extracellular melibiose (a disaccharide of galactose and glucose). We show that a progenitor 5-gene cluster containing GAL 7-1-10-4-2 was likely present in the common ancestor of Torulaspora and Zygotorulaspora. It added PGM1 and MEL1 in the ancestor of most Torulaspora species. It underwent further expansion in the T. pretoriensis clade, involving the fusion of three progenitor clusters in tandem and the gain of HGT1. These giant GAL clusters are highly polymorphic in structure, and subject to horizontal transfers, pseudogenization and gene losses. We identify recent horizontal transfers of complete GAL clusters from T. franciscae into one strain of T. delbrueckii, and from a relative of T. maleeae into one strain of T. globosa. The variability and dynamic evolution of GAL clusters in Torulaspora indicates that there is strong natural selection on the GAL pathway in this genus.

Sequencing platform shifts provide opportunities but pose challenges for combining genomic data sets

Technological advances in DNA sequencing over the last decade now permit the production and curation of large genomic data sets in an increasing number of nonmodel species. Additionally, these new data provide the opportunity for combining data sets, resulting in larger studies with a broader taxonomic range. Whilst the development of new sequencing platforms has been beneficial, resulting in a higher throughput of data at a lower per-base cost, shifts in sequencing technology can also pose challenges for those wishing to combine new sequencing data with data sequenced on older platforms. Here, we outline the types of studies where the use of curated data might be beneficial, and highlight potential biases that might be introduced by combining data from different sequencing platforms. As an example of the challenges associated with combining data across sequencing platforms, we focus on the impact of the shift in Illumina's base calling technology from a four-channel system to a two-channel system. We caution that when data are combined from these two systems, erroneous guanine base calls that result from the two-channel chemistry can make their way through a bioinformatic pipeline, eventually leading to inaccurate and potentially misleading conclusions. We also suggest solutions for dealing with such potential artefacts, which make samples sequenced on different sequencing platforms appear more differentiated from one another than they really are. Finally, we stress the importance of archiving tissue samples and the associated sequences for the continued reproducibility and reusability of sequencing data in the face of ever-changing sequencing platform technology.

Step-by-Step Evolution of Telomeres: Lessons from Yeasts

In virtually every eukaryotic species, the ends of nuclear chromosomes are protected by telomeres, nucleoprotein structures counteracting the end-replication problem and suppressing recombination and undue DNA repair. Although in most cases, the primary structure of telomeric DNA is conserved, there are several exceptions to this rule. One is represented by the telomeric repeats of ascomycetous yeasts, which encompass a great variety of sequences, whose evolutionary origin has been puzzling for several decades. At present, the key questions concerning the driving force behind their rapid evolution and the means of co-evolution of telomeric repeats and telomere-binding proteins remain largely unanswered. Previously published studies addressed mostly the general concepts of the evolutionary origin of telomeres, key properties of telomeric proteins as well as the molecular mechanisms of telomere maintenance; however, the evolutionary process itself has not been analyzed thoroughly. Here, we aimed to inspect the evolution of telomeres in ascomycetous yeasts from the subphyla Saccharomycotina and Taphrinomycotina, with special focus on the evolutionary origin of species-specific telomeric repeats. We analyzed the sequences of telomeric repeats from 204 yeast species classified into 20 families and as a result, we propose a step-by-step model, which integrates the diversity of telomeric repeats, telomerase RNAs, telomere-binding protein complexes and explains a propensity of certain species to generate the repeat heterogeneity within a single telomeric array.

A Noncanonical DNA Damage Checkpoint Response in a Major Fungal Pathogen

In order to preserve genome integrity, all cells must mount appropriate responses to DNA damage, including slowing down or arresting the cell cycle to give the cells time to repair the damage and changing gene expression, for example to induce genes involved in DNA repair. The Rad53 protein kinase is a conserved central mediator of these responses in eukaryotic cells, and its extensive phosphorylation upon DNA damage is necessary for its activation and subsequent activity.

DNA damage checkpoints are key guardians of genome integrity. Eukaryotic cells respond to DNA damage by triggering extensive phosphorylation of Rad53/CHK2 effector kinase, whereupon activated Rad53/CHK2 mediates further aspects of checkpoint activation, including cell cycle arrest and transcriptional changes. Budding yeast Candida glabrata, closely related to model eukaryote Saccharomyces cerevisiae, is an opportunistic pathogen characterized by high genetic diversity and rapid emergence of drug-resistant mutants. However, the mechanisms underlying this genetic variability are unclear. We used Western blotting and mass spectrometry to show that, unlike S. cerevisiae, C. glabrata cells exposed to DNA damage did not induce C. glabrata Rad53 (CgRad53) phosphorylation. Furthermore, flow cytometry analysis showed that, unlike S. cerevisiae, C. glabrata cells did not accumulate in S phase upon DNA damage. Consistent with these observations, time-lapse microscopy showed C. glabrata cells continuing to divide in the presence of DNA damage, resulting in mitotic errors and cell death. Finally, transcriptome sequencing (RNAseq) analysis revealed transcriptional rewiring of the DNA damage response in C. glabrata and identified several key protectors of genome stability upregulated by DNA damage in S. cerevisiae but downregulated in C. glabrata, including proliferating cell nuclear antigen (PCNA). Together, our results reveal a noncanonical fungal DNA damage response in C. glabrata, which may contribute to rapidly generating genetic change and drug resistance.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

An investigation of irreproducibility in maximum likelihood phylogenetic inference

Phylogenetic trees are essential for studying biology, but their reproducibility under identical parameter settings remains unexplored. Here, we find that 3515 (18.11%) IQ-TREE-inferred and 1813 (9.34%) RAxML-NG-inferred maximum likelihood (ML) gene trees are topologically irreproducible when executing two replicates (Run1 and Run2) for each of 19,414 gene alignments in 15 animal, plant, and fungal phylogenomic datasets. Notably, coalescent-based ASTRAL species phylogenies inferred from Run1 and Run2 sets of individual gene trees are topologically irreproducible for 9/15 phylogenomic datasets, whereas concatenation-based phylogenies inferred twice from the same supermatrix are reproducible. Our simulations further show that irreproducible phylogenies are more likely to be incorrect than reproducible phylogenies. These results suggest that a considerable fraction of single-gene ML trees may be irreproducible. Increasing reproducibility in ML inference will benefit from providing analyses’ log files, which contain typically reported parameters (e.g., program, substitution model, number of tree searches) but also typically unreported ones (e.g., random starting seed number, number of threads, processor type).

Replicate runs of maximum likelihood phylogenetic analyses can generate different tree topologies due to differences in parameters, such as random seeds. Here, Shen et al. demonstrate that replicate runs can generate substantially different tree topologies even with identical data and parameters.

Evolutionary Overview of Molecular Interactions and Enzymatic Activities in the Yeast Cell Walls

Fungal cell walls are composed of a polysaccharide network that serves as a scaffold in which different glycoproteins are embedded. Investigation of fungal cell walls, besides simple identification and characterization of the main cell wall building blocks, covers the pathways and regulations of synthesis of each individual component of the wall and biochemical reactions by which they are cross-linked and remodeled in response to different growth phase and environmental signals. In this review, a survey of composition and organization of so far identified and characterized cell wall components of different yeast genera including Saccharomyces, Candida, Kluyveromyces, Yarrowia, and Schizosaccharomyces are presented with the focus on their cell wall proteomes.

A yeast living ancestor reveals the origin of genomic introgressions

Genome introgressions drive evolution across the animal¹, plant² and fungal³ kingdoms. Introgressions initiate from archaic admixtures followed by repeated backcrossing to one parental species. However, how introgressions arise in reproductively isolated species, such as yeast⁴, has remained unclear. Here we identify a clonal descendant of the ancestral yeast hybrid that founded the extant Saccharomyces cerevisiae Alpechin lineage⁵, which carries abundant Saccharomyces paradoxus introgressions. We show that this clonal descendant, hereafter defined as a 'living ancestor', retained the ancestral genome structure of the first-generation hybrid with contiguous S. cerevisiae and S. paradoxus subgenomes. The ancestral first-generation hybrid underwent catastrophic genomic instability through more than a hundred mitotic recombination events, mainly manifesting as homozygous genome blocks generated by loss of heterozygosity. These homozygous sequence blocks rescue hybrid fertility by restoring meiotic recombination and are the direct origins of the introgressions present in the Alpechin lineage. We suggest a plausible route for introgression evolution through the reconstruction of extinct stages and propose that genome instability allows hybrids to overcome reproductive isolation and enables introgressions to emerge.

Genomic and phenotypic analysis of COVID-19-associated pulmonary aspergillosis isolates of Aspergillus fumigatus

The ongoing global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) first described from Wuhan, China. A subset of COVID-19 patients has been reported to have acquired secondary infections by microbial pathogens, such as fungal opportunistic pathogens from the genus Aspergillus . To gain insight into COVID-19 associated pulmonary aspergillosis (CAPA), we analyzed the genomes and characterized the phenotypic profiles of four CAPA isolates of Aspergillus fumigatus obtained from patients treated in the area of North Rhine-Westphalia, Germany. By examining the mutational spectrum of single nucleotide polymorphisms, insertion-deletion polymorphisms, and copy number variants among 206 genes known to modulate A. fumigatus virulence, we found that CAPA isolate genomes do not exhibit major differences from the genome of the Af293 reference strain. By examining virulence in an invertebrate moth model, growth in the presence of osmotic, cell wall, and oxidative stressors, and the minimum inhibitory concentration of antifungal drugs, we found that CAPA isolates were generally, but not always, similar to A. fumigatus reference strains Af293 and CEA17. Notably, CAPA isolate D had more putative loss of function mutations in genes known to increase virulence when deleted (e.g., in the FLEA gene, which encodes a lectin recognized by macrophages). Moreover, CAPA isolate D was significantly more virulent than the other three CAPA isolates and the A. fumigatus reference strains tested. These findings expand our understanding of the genomic and phenotypic characteristics of isolates that cause CAPA.

Genome-scale phylogeny and contrasting modes of genome evolution in the fungal phylum Ascomycota

Ascomycota, the largest and most well-studied phylum of fungi, contains three subphyla: Saccharomycotina (budding yeasts), Pezizomycotina (filamentous fungi), and Taphrinomycotina (fission yeasts). Despite its importance, we lack a comprehensive genome-scale phylogeny or understanding of the similarities and differences in the mode of genome evolution within this phylum. By examining 1107 genomes from Saccharomycotina (332), Pezizomycotina (761), and Taphrinomycotina (14) species, we inferred a robust genome-wide phylogeny that resolves several contentious relationships and estimated that the Ascomycota last common ancestor likely originated in the Ediacaran period. Comparisons of genomic properties revealed that Saccharomycotina and Pezizomycotina differ greatly in their genome properties and enabled inference of the direction of evolutionary change. The Saccharomycotina typically have smaller genomes, lower guanine-cytosine contents, lower numbers of genes, and higher rates of molecular sequence evolution compared with Pezizomycotina. These results provide a robust evolutionary framework for understanding the diversity and ecological lifestyles of the largest fungal phylum.

Candida metrosideri pro tempore sp. nov. and Candida ohialehuae pro tempore sp. nov., two antifungal-resistant yeasts associated with Metrosideros polymorpha flowers in Hawaii

Flowers produce an array of nutrient-rich exudates in which microbes can thrive, making them hotspots for microbial abundance and diversity. During a diversity study of yeasts inhabiting the flowers of Metrosideros polymorpha (Myrtaceae) in the Hawai’i Volcanoes National Park (HI, USA), five isolates were found to represent two novel species. Morphological and physiological characterization, and sequence analysis of the small subunit ribosomal RNA (rRNA) genes, the D1/D2 domains of the large subunit rRNA genes, the internal transcribed spacer (ITS) regions, and the genes encoding the largest and second largest subunits of the RNA polymerase II (RPB1 and RPB2, respectively), classified both species in the family Metschnikowiaceae, and we propose the names Candida metrosideri pro tempore sp. nov. (JK22^T = CBS 16091 = MUCL 57821) and Candida ohialehuae pro tempore sp. nov. (JK58.2^T = CBS 16092 = MUCL 57822) for such new taxa. Both novel Candida species form a well-supported subclade in the Metschnikowiaceae containing species associated with insects, flowers, and a few species of clinical importance. The ascosporic state of the novel species was not observed. The two novel yeast species showed elevated minimum inhibitory concentrations to the antifungal drug amphotericin B (>4 μg/mL). The ecology and phylogenetic relationships of C. metrosideri and C. ohialehuae are also discussed.

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

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

Rediscovering a Forgotten System of Symbiosis: Historical Perspective and Future Potential

While the majority of symbiosis research is focused on bacteria, microbial eukaryotes play important roles in the microbiota and as pathogens, especially the incredibly diverse Fungi kingdom. The recent emergence of widespread pathogens in wildlife (bats, amphibians, snakes) and multidrug-resistant opportunists in human populations (Candida auris) has highlighted the importance of better understanding animal-fungus interactions. Regardless of their prominence there are few animal-fungus symbiosis models, but modern technological advances are allowing researchers to utilize novel organisms and systems. Here, I review a forgotten system of animal-fungus interactions: the beetle-fungus symbioses of Drugstore and Cigarette beetles with their symbiont Symbiotaphrina. As pioneering systems for the study of mutualistic symbioses, they were heavily researched between 1920 and 1970, but have received only sporadic attention in the past 40 years. Several features make them unique research organisms, including (1) the symbiont is both extracellular and intracellular during the life cycle of the host, and (2) both beetle and fungus can be cultured in isolation. Specifically, fungal symbionts intracellularly infect cells in the larval and adult beetle gut, while accessory glands in adult females harbor extracellular fungi. In this way, research on the microbiota, pathogenesis/infection, and mutualism can be performed. Furthermore, these beetles are economically important stored-product pests found worldwide. In addition to providing a historical perspective of the research undertaken and an overview of beetle biology and their symbiosis with Symbiotaphrina, I performed two analyses on publicly available genomic data. First, in a preliminary comparative genomic analysis of the fungal symbionts, I found striking differences in the pathways for the biosynthesis of two B vitamins important for the host beetle, thiamine and biotin. Second, I estimated the most recent common ancestor for Drugstore and Cigarette beetles at 8.8-13.5 Mya using sequence divergence (CO1 gene). Together, these analyses demonstrate that modern methods and data (genomics, transcriptomes, etc.) have great potential to transform these beetle-fungus systems into model systems again.

Unlocking a signal of introgression from codons in Lachancea kluyveri using a mutation-selection model

Background

For decades, codon usage has been used as a measure of adaptation for translational efficiency and translation accuracy of a gene’s coding sequence. These patterns of codon usage reflect both the selective and mutational environment in which the coding sequences evolved. Over this same period, gene transfer between lineages has become widely recognized as an important biological phenomenon. Nevertheless, most studies of codon usage implicitly assume that all genes within a genome evolved under the same selective and mutational environment, an assumption violated when introgression occurs. In order to better understand the effects of introgression on codon usage patterns and vice versa, we examine the patterns of codon usage in Lachancea kluyveri, a yeast which has experienced a large introgression. We quantify the effects of mutation bias and selection for translation efficiency on the codon usage pattern of the endogenous and introgressed exogenous genes using a Bayesian mixture model, ROC SEMPPR, which is built on mechanistic assumptions about protein synthesis and grounded in population genetics.

Results

We find substantial differences in codon usage between the endogenous and exogenous genes, and show that these differences can be largely attributed to differences in mutation bias favoring A/T ending codons in the endogenous genes while favoring C/G ending codons in the exogenous genes. Recognizing the two different signatures of mutation bias and selection improves our ability to predict protein synthesis rate by 42% and allowed us to accurately assess the decaying signal of endogenous codon mutation and preferences. In addition, using our estimates of mutation bias and selection, we identify Eremothecium gossypii as the closest relative to the exogenous genes, providing an alternative hypothesis about the origin of the exogenous genes, estimate that the introgression occurred ∼6×10⁸ generation ago, and estimate its historic and current selection against mismatched codon usage.

Conclusions

Our work illustrates how mechanistic, population genetic models like ROC SEMPPR can separate the effects of mutation and selection on codon usage and provide quantitative estimates from sequence data.

In Kluyveromyces lactis a Pair of Paralogous Isozymes Catalyze the First Committed Step of Leucine Biosynthesis in Either the Mitochondria or the Cytosol

Divergence of paralogous pairs, resulting from gene duplication, plays an important role in the evolution of specialized or novel gene functions. Analysis of selected duplicated pairs has elucidated some of the mechanisms underlying the functional diversification of Saccharomyces cerevisiae (S. cerevisiae) paralogous genes. Similar studies of the orthologous pairs extant in pre-whole genome duplication yeast species, such as Kluyveromyces lactis (K. lactis) remain to be addressed. The genome of K. lactis, an aerobic yeast, includes gene pairs generated by sporadic duplications. The genome of this organism comprises the KlLEU4 and KlLEU4BIS paralogous pair, annotated as putative α-isopropylmalate synthases (α-IPMSs), considered to be the orthologs of the S. cerevisiae ScLEU4/ScLEU9 paralogous genes. The enzymes encoded by the latter two genes are mitochondrially located, differing in their sensitivity to leucine allosteric inhibition resulting in ScLeu4-ScLeu4 and ScLeu4-ScLeu9 sensitive dimers and ScLeu9-ScLeu9 relatively resistant homodimers. Previous work has shown that, in a Scleu4Δ mutant, ScLEU9 expression is increased and assembly of ScLeu9-ScLeu9 leucine resistant homodimers results in loss of feedback regulation of leucine biosynthesis, leading to leucine accumulation and decreased growth rate. Here we report that: (i) K. lactis harbors a sporadic gene duplication, comprising the KlLEU4, syntenic with S. cerevisiae ScLEU4 and ScLEU9, and the non-syntenic KlLEU4BIS, arising from a pre-WGD event. (ii) That both, KlLEU4 and KlLEU4BIS encode leucine sensitive α-IPMSs isozymes, located in the mitochondria (KlLeu4) and the cytosol (KlLeu4BIS), respectively. (iii) That both, KlLEU4 or KlLEU4BIS complement the Scleu4Δ Scleu9Δ leucine auxotrophic phenotype and revert the enhanced ScLEU9 transcription observed in a Scleu4Δ ScLEU9 mutant. The Scleu4Δ ScLEU9 growth mutant phenotype is only fully complemented when transformed with the syntenic KlLEU4 mitochondrial isoform. KlLEU4 and KlLEU4BIS underwent a different diversification pathways than that leading to ScLEU4/ScLEU9. KlLEU4 could be considered as the functional ortholog of ScLEU4, since its encoded isozyme can complement both the Scleu4Δ Scleu9Δ leucine auxotrophy and the Scleu4Δ ScLEU9 complex phenotype.

Suhomyces drosophilae sp. nov., isolated from Drosophila flies feeding on a stinkhorn mushroom

The majority of Suhomyces species have been isolated from fungus-feeding insects and particularly from the gut of beetles. In the present study, seven yeast strains were isolated from the gut of Drosophila species feeding on gleba, the spore-bearing inner mass, of a stinkhorn mushroom belonging to the family Phallaceae. Based on phenotypic, biochemical characterization and sequence analysis of the D1/D2 region of the large subunit rRNA gene and the internal transcribed spacer (ITS) region, two of these yeast strains, DGY3 and DGY4, represented a novel species of the genus Suhomyces. The novel species is closely related to an undescribed species of Candida ST-370 (DQ404513) and with Suhomyces canberraensis, wherein, the novel species differs from S. canberraensis by 40 nucleotide substitutions and three gaps (7.7 % sequence variation) in the D1/D2 region and 50 nucleotide substitutions and seven gaps (13.7 % sequence variation) in the ITS region. Several morphological and physiological differences were also observed between S. canberraensis and the strains obtained during this study. These data support the proposal of Suhomyces drosophilae as a novel species, with DGY3^T as the holotype and CBS 16329^T and MCC 1871^T as ex-type strains.

Towards yeast taxogenomics: lessons from novel species descriptions based on complete genome sequences

In recent years, ‘multi-omic’ sciences have affected all aspects of fundamental and applied biological research. Yeast taxonomists, though somewhat timidly, have begun to incorporate complete genomic sequences into the description of novel taxa, taking advantage of these powerful data to calculate more reliable genetic distances, construct more robust phylogenies, correlate genotype with phenotype and even reveal cryptic sexual behaviors. However, the use of genomic data in formal yeast species descriptions is far from widespread. The present review examines published examples of genome-based species descriptions of yeasts, highlights relevant bioinformatic approaches, provides recommendations for new users and discusses some of the challenges facing the genome-based systematics of yeasts.

Toward a Fully Resolved Fungal Tree of Life

In this review, we discuss the current status and future challenges for fully elucidating the fungal tree of life. In the last 15 years, advances in genomic technologies have revolutionized fungal systematics, ushering the field into the phylogenomic era. This has made the unthinkable possible, namely access to the entire genetic record of all known extant taxa. We first review the current status of the fungal tree and highlight areas where additional effort will be required. We then review the analytical challenges imposed by the volume of data and discuss methods to recover the most accurate species tree given the sea of gene trees. Highly resolved and deeply sampled trees are being leveraged in novel ways to study fungal radiations, species delimitation, and metabolic evolution. Finally, we discuss the critical issue of incorporating the unnamed and uncultured dark matter taxa that represent the vast majority of fungal diversity.

Hybridization and the origin of new yeast lineages

Hybrids originate from the mating of two diverged organisms, resulting in novel lineages that have chimeric genomes. Hybrids may exhibit unique phenotypic traits that are not necessarily intermediate between those present in the progenitors. These unique traits may enable them to thrive in new environments. Many hybrid lineages have been discovered among yeasts in the Saccharomycotina, of which many have industrial or clinical relevance, but this might reflect a bias toward investigating species with relevance to humans. Hybridization has also been proposed to be at the root of the whole-genome duplication in the lineage leading to Saccharomyces cerevisiae. Thus, hybridization seems to have played a prominent role in the evolution of Saccharomycotina yeasts, although it is still unclear how common this evolutionary process has been during the evolution of this and other fungal clades. Similarly, the evolutionary aftermath of hybridization, including implications at the genomic, transcriptional, physiological or ecological levels, remains poorly understood. In this review, I survey recent findings from genomic analysis of yeast hybrids of industrial or clinical relevance, and discuss the evolutionary implications of genomic hybridization for the origin of new lineages, including when such hybridization results in a whole-genome duplication.

Kluyveromyces osmophilus is not a synonym of Zygosaccharomyces mellis; reinstatement as Zygosaccharomyces osmophilus comb. nov

Kluyveromyces osmophilus, a single-strain species isolated from Mozambique sugar, has been treated a synonym of Zygosaccharomyces mellis. Analyses of D1/D2 LSU rRNA gene sequences confirmed that the species belongs to the genus Zygosaccharomyces but showed it to be distinct from strains of Z. mellis. During studies of yeasts associated with stingless bees in Brazil, nine additional isolates of the species were obtained from unripe and ripe honey and pollen of Scaptotrigona cfr. bipunctata, as well as ripe honey of Tetragonisca angustula. The D1/D2 sequences of the Brazilian isolates were identical to those of the type strain of K. osmophilus CBS 5499 (=ATCC 22027), indicating that they represent the same species. Phylogenomic analyses using 4038 orthologous genes support the reinstatement of K. osmophilus as a member of the genus Zygosaccharomyces. We, therefore, propose the name Zygosaccharomyces osmophilus comb. nov. (lectotype ATCC 22027; MycoBank no. MB 833739).

Starmera foglemanii sp. nov. and Starmera ilhagrandensis sp. nov., two novel yeast species isolated from ephemeral plant substrates

Four isolates of two novel ascosporogenous species belonging to the clade Starmera were obtained from cactus tissues and rotting wood in Brazil. Results of analyses of the sequences of the ITS and D1/D2 domains of the large subunit rRNA gene indicated that the two isolates of the cactophilic species are related to Starmera caribaea and Starmera pilosocereana, yeasts that are associated with cacti and require an organic source of sulfur for growth. We propose the novel species Starmera foglemanii sp. nov. (CBS 16113T; MycoBank number: MB 834400) to accommodate these isolates. The other two isolates are phylogenetically related to Candida dendrica, Candida laemsonensis and Candida berthetii, also in the Starmera clade. The novel species name Starmera ilhagrandensis sp. nov. (CBS 16316T; MycoBank number: MB 834402) is proposed for this species.

The Methylotroph Gene Order Browser (MGOB) reveals conserved synteny and ancestral centromere locations in the yeast family Pichiaceae

The yeast family Pichiaceae, also known as the 'methylotrophs clade', is a relatively little studied group of yeasts despite its economic and clinical relevance. To explore the genome evolution and synteny relationships within this family, we developed the Methylotroph Gene Order Browser (MGOB, http://mgob.ucd.ie) similar to our previous gene order browsers for other yeast families. The dataset contains genome sequences from nine Pichiaceae species, including our recent reference sequence of Pichia kudriavzevii. As an example, we demonstrate the conservation of synteny around the MOX1 locus among species both containing and lacking the MOX1 gene for methanol assimilation. We found ancient clusters of genes that are conserved as adjacent between Pichiaceae and Saccharomycetaceae. Surprisingly, we found evidence that the locations of some centromeres have been conserved among Pichiaceae species, and between Pichiaceae and Saccharomycetaceae, even though the centromeres fall into different structural categories-point centromeres, inverted repeats and retrotransposon cluster centromeres.

Takashi Nakase's last tweet: what is the current direction of microbial taxonomy research?

During the last few decades, type strains of most yeast species have been barcoded using the D1/D2 domain of their LSU rRNA gene and internal transcribed spacer (ITS) region. Species identification using DNA sequences regarding conspecificity in yeasts has also been studied. Most yeast species can be identified according to the sequence divergence of their ITS region or a combination of the D1/D2 and ITS regions. Studies that have examined intraspecific diversity have used multilocus sequence analyses, whereas the marker regions used in this analysis vary depending upon taxa. D1/D2 domain and ITS region sequences have been used as barcodes to develop primers suitable for the detection of the biological diversity of environmental DNA and the microbiome. Using these barcode sequences, it is possible to identify relative lineages and infer their gene products and function, and how they adapt to their environment. If barcode sequence was not variable enough to identify a described species, one could investigate the other biological traits of these yeasts, considering geological distance, environmental circumstances and isolation of reproduction. This article is dedicated to late Dr Takashi Nakase (1939-2018).

Osh6 requires Ist2 for localization to ER-PM contacts and efficient phosphatidylserine transport in budding yeast

Osh6 and Osh7 are lipid transfer proteins (LTPs) that move phosphatidylserine (PS) from the endoplasmic reticulum (ER) to the plasma membrane (PM). High PS levels at the PM are key for many cellular functions. Intriguingly, Osh6 and Osh7 localize to ER-PM contact sites, although they lack membrane-targeting motifs, in contrast to multidomain LTPs that both bridge membranes and convey lipids. We show that Osh6 localization to contact sites depends on its interaction with the cytosolic tail of the ER-PM tether Ist2, a homolog of TMEM16 proteins. We identify a motif in the Ist2 tail, conserved in yeasts, as the Osh6-binding region, and we map an Ist2-binding surface on Osh6. Mutations in the Ist2 tail phenocopy osh6Δ osh7Δ deletion: they decrease cellular PS levels and block PS transport to the PM. Our study unveils an unexpected partnership between a TMEM16-like protein and a soluble LTP, which together mediate lipid transport at contact sites.This article has an associated First Person interview with the first author of the paper.

Performance, genomic rearrangements, and signatures of adaptive evolution: Lessons from fermentative yeasts

The capacity of some yeasts to extract energy from single sugars, generating CO2 and ethanol (=fermentation), even in the presence of oxygen, is known as the Crabtree effect. This phenomenon represents an important adaptation as it allowed the utilization of the ecological niche given by modern fruits, an abundant source of food that emerged in the terrestrial environment in the Cretaceous. However, identifying the evolutionary events that triggered fermentative capacity in Crabtree-positive species is challenging, as microorganisms do not leave fossil evidence. Thus, key innovations should be inferred based only on traits measured under culture conditions. Here, we reanalyzed data from a common garden experiment where several proxies of fermentative capacity were recorded in Crabtree-positive and Crabtree-negative species, representing yeast phylogenetic diversity. In particular, we applied the "lasso-OU" algorithm which detects points of adaptive shifts, using traits that are proxies of fermentative performance. We tested whether multiple events or a single event explains the actual fermentative capacity of yeasts. According to the lasso-OU procedure, evolutionary changes in the three proxies of fermentative capacity that we considered (i.e., glycerol production, ethanol yield, and respiratory quotient) are consistent with a single evolutionary episode (a whole-genomic duplication, WGD), instead of a series of small genomic rearrangements. Thus, the WGD appears as the key event behind the diversification of fermentative yeasts, which by increasing gene dosage, and maximized their capacity of energy extraction for exploiting the new ecological niche provided by single sugars.

Divergence of Peroxisome Membrane Gene Sequence and Expression Between Yeast Species

Large population-genomic sequencing studies can enable highly-powered analyses of sequence signatures of natural selection. Genome repositories now available for Saccharomyces yeast make it a premier model for studies of the molecular mechanisms of adaptation. We mined the genomes of hundreds of isolates of the sister species S. cerevisiae and S. paradoxus to identify sequence hallmarks of adaptive divergence between the two. From the top hits we focused on a set of genes encoding membrane proteins of the peroxisome, an organelle devoted to lipid breakdown and other specialized metabolic pathways. In-depth population- and comparative-genomic sequence analyses of these genes revealed striking divergence between S. cerevisiae and S. paradoxus. And from transcriptional profiles we detected non-neutral, directional cis-regulatory variation at the peroxisome membrane genes, with overall high expression in S. cerevisiae relative to S. paradoxus. Taken together, these data support a model in which yeast species have differentially tuned the expression of peroxisome components to boost their fitness in distinct niches.

Spatial inter-centromeric interactions facilitated the emergence of evolutionary new centromeres

Centromeres of Candida albicans form on unique and different DNA sequences but a closely related species, Candida tropicalis, possesses homogenized inverted repeat (HIR)-associated centromeres. To investigate the mechanism of centromere type transition, we improved the fragmented genome assembly and constructed a chromosome-level genome assembly of C. tropicalis by employing PacBio sequencing, chromosome conformation capture sequencing (3C-seq), chromoblot, and genetic analysis of engineered aneuploid strains. Further, we analyzed the 3D genome organization using 3C-seq data, which revealed spatial proximity among the centromeres as well as telomeres of seven chromosomes in C. tropicalis. Intriguingly, we observed evidence of inter-centromeric translocations in the common ancestor of C. albicans and C. tropicalis. Identification of putative centromeres in closely related Candida sojae, Candida viswanathii and Candida parapsilosis indicates loss of ancestral HIR-associated centromeres and establishment of evolutionary new centromeres (ENCs) in C. albicans. We propose that spatial proximity of the homologous centromere DNA sequences facilitated karyotype rearrangements and centromere type transitions in human pathogenic yeasts of the CUG-Ser1 clade.

Draft Genome Sequence of the Yeast Kazachstania telluris CBS 16338 Isolated from Forest Soil in Ireland

Draft genomes of several Kazachstania species, a large group of ascomycetous budding yeasts, have been recently published. However, there is none yet available for pathogenic species from the K. telluris clade, including K. telluris sensu stricto, an opportunistic yeast that has been isolated from various niches and human clinical samples. Here we provide the first draft genome sequence of the strain K. telluris CBS 16338, that we isolated from forest soil in Ireland.

Stress-induced expression is enriched for evolutionarily young genes in diverse budding yeasts

The Saccharomycotina subphylum (budding yeasts) spans 400 million years of evolution and includes species that thrive in diverse environments. To study niche-adaptation, we identify changes in gene expression in three divergent yeasts grown in the presence of various stressors. Duplicated and non-conserved genes are significantly more likely to respond to stress than genes that are conserved as single-copy orthologs. Next, we develop a sorting method that considers evolutionary origin and duplication timing to assign an evolutionary age to each gene. Subsequent analysis reveals that genes that emerged in recent evolutionary time are enriched amongst stress-responsive genes for each species. This gene expression pattern suggests that budding yeasts share a stress adaptation mechanism, whereby selective pressure leads to functionalization of young genes to improve growth in adverse conditions. Further characterization of young genes from species that thrive in harsh environments can inform the design of more robust strains for biotechnology.

Fermentation parameters of industrial processes are often not the ideal growth conditions for industrial microbes. Here, the authors reveal that young genes are more responsive to environmental stress than ancient genes using a new gene age assignment method and provide targeted genes for metabolic engineering.

The yeast mating-type switching endonuclease HO is a domesticated member of an unorthodox homing genetic element family

The mating-type switching endonuclease HO plays a central role in the natural life cycle of Saccharomyces cerevisiae, but its evolutionary origin is unknown. HO is a recent addition to yeast genomes, present in only a few genera close to Saccharomyces. Here we show that HO is structurally and phylogenetically related to a family of unorthodox homing genetic elements found in Torulaspora and Lachancea yeasts. These WHO elements home into the aldolase gene FBA1, replacing its 3' end each time they integrate. They resemble inteins but they operate by a different mechanism that does not require protein splicing. We show that a WHO protein cleaves Torulaspora delbrueckii FBA1 efficiently and in an allele-specific manner, leading to DNA repair by gene conversion or NHEJ. The DNA rearrangement steps during WHO element homing are very similar to those during mating-type switching, and indicate that HO is a domesticated WHO-like element.

In the same way as a sperm from a male and an egg from a female join together to form an embryo in most animals, yeast cells have two sexes that coordinate how they reproduce. These are called “mating types” and, rather than male or female, an individual yeast cell can either be mating type “a” or “alpha”. Every yeast cell contains the genes for both mating types, and each cell’s mating type is determined by which of those genes it has active.

Only one mating type gene can be ‘on’ at a time, but some yeast species can swap mating type on demand by switching the corresponding genes ‘on’ or ‘off’. This switch is unusual. Rather than simply activate one of the genes it already has, the yeast cell keeps an inactive version of each mating type gene tucked away, makes a copy of the gene it wants to be active and pastes that copy into a different location in its genome. To do all of this yeast need another gene called HO. This gene codes for an enzyme that cuts the DNA at the location of the active mating type gene. This makes an opening that allows the cell to replace the ‘a’ gene with the ‘alpha’ gene, or vice versa. This system allows yeast cells to continue mating even if all the cells in a colony start off as the same mating type. But, cutting into the DNA is risky, and can damage the health of the cell. So, why did yeast cells evolve a system that could cause them harm?

To find out where the HO gene came from, Coughlan et al. searched through all the available genomes from yeast species for other genes with similar sequences and identified a cluster which they nicknamed “weird HO” genes, or WHO genes for short. Testing these genes revealed that they also code for enzymes that make cuts in the yeast genome, but the way the cell repairs the cuts is different. The WHO genes are jumping genes. When the enzyme encoded by a WHO gene makes a cut in the genome, the yeast cell copies the gene into the gap, allowing the gene to ‘jump’ from one part of the genome to another. It is possible that this was the starting point for the evolution of the HO gene. Changes to a WHO gene could have allowed it to cut into the mating type region of the yeast genome, giving the yeast an opportunity to ‘domesticate’ it. Over time, the yeast cell stopped the WHO gene from jumping into the gap and started using the cut to change its mating type.

Understanding how cells adapt genes for different purposes is a key question in evolutionary biology. There are many other examples of domesticated jumping genes in other organisms, including in the human immune system. Understanding the evolution of HO not only sheds light on how yeast mating type switching evolved, but on how other species might harness and adapt their genes.