Molecular phylogeny of sequenced Saccharomycetes reveals polyphyly of the alternative yeast codon usage

The yeast-human coevolution: Fungal transition from passengers, colonizers, and invaders

Fungi are the cause of more than a billion infections in humans every year, although their interactions with the host are still neglected compared to bacteria. Major systemic fungal infections are very unusual in the healthy population, due to the long history of coevolution with the human host. Humans are routinely exposed to environmental fungi and can host a commensal mycobiota, which is increasingly considered as a key player in health and disease. Here, we review the current knowledge on host-fungi coevolution and the factors that regulate their interaction. On one hand, fungi have learned to survive and inhabit the host organisms as a natural ecosystem, on the other hand, the host immune system finely tunes the response toward fungi. In turn, recognition of fungi as commensals or pathogens regulates the host immune balance in health and disease. In the human gut ecosystem, yeasts provide a fingerprint of the transient microbiota. Their status as passengers or colonizers is related to the integrity of the gut barrier and the risk of multiple disorders. Thus, the study of this less known component of the microbiota could unravel the rules of the transition from passengers to colonizers and invaders, as well as their dependence on the innate component of the host's immune response. This article is categorized under: Infectious Diseases > Environmental Factors Immune System Diseases > Environmental Factors Infectious Diseases > Molecular and Cellular Physiology.

Genome-wide profiling of piggyBac transposon insertion mutants reveals loss of the F1F0 ATPase complex causes fluconazole resistance in Candida glabrata

Invasive candidiasis caused by non-albicans species has been on the rise, with Candida glabrata emerging as the second most common etiological agent. Candida glabrata possesses an intrinsically lower susceptibility to azoles and an alarming propensity to rapidly develop high-level azole resistance during treatment. In this study, we have developed an efficient piggyBac (PB) transposon-mediated mutagenesis system in C. glabrata to conduct genome-wide genetic screens and applied it to profile genes that contribute to azole resistance. When challenged with the antifungal drug fluconazole, PB insertion into 270 genes led to significant resistance. A large subset of these genes has a role in the mitochondria, including almost all genes encoding the subunits of the F1F0 ATPase complex. We show that deleting ATP3 or ATP22 results in increased azole resistance but does not affect susceptibility to polyenes and echinocandins. The increased azole resistance is due to increased expression of PDR1 that encodes a transcription factor known to promote drug efflux pump expression. Deleting PDR1 in the atp3Δ or atp22Δ mutant resulted in hypersensitivity to fluconazole. Our results shed light on the mechanisms contributing to azole resistance in C. glabrata. This PB transposon-mediated mutagenesis system can significantly facilitate future genome-wide genetic screens.

Current advances in Candida tropicalis: Yeast overview and biotechnological applications

Candida tropicalis is a nonconventional yeast with medical and industrial significance, belonging to the CTG clade. Recent advancements in whole-genome sequencing and genetic analysis revealed its close relation to other unconventional yeasts of biotechnological importance. C. tropicalis is known for its immense potential in synthesizing various valuable biomolecules such as ethanol, xylitol, biosurfactants, lipids, enzymes, α,ω-dicarboxylic acids, single-cell proteins, and more, making it an attractive target for biotechnological applications. This review provides an update on C. tropicalis biological characteristics and its efficiency in producing a diverse range of biomolecules with industrial significance from various feedstocks. The information presented in this review contributes to a better understanding of C. tropicalis and highlights its potential for biotechnological applications and market viability.

Functional mapping of the N-terminal region of the yeast moonlighting protein Sis2/Hal3 reveals crucial residues for Ppz1 regulation

The function of the Saccharomyces cerevisiae Ppz1 phosphatase is controlled by its inhibitory subunit Hal3. Hal3 is a moonlighting protein, which associates with Cab3 to form a decarboxylase involved in the CoA biosynthetic pathway. Hal3 is composed by a conserved core PD region, required for both Ppz1 regulation and CoA biosynthesis, a long N-terminal extension, and an acidic C-terminal tail. Cab3 has a similar structure, but it is not a Ppz1 inhibitor. We show here that deletion or specific mutations in a short region of the N-terminal extension of Hal3 compromise its function as a Ppz1 inhibitor in vivo and in vitro without negatively affecting its ability to interact with the phosphatase. This study defines a R-K-X₍₃₎ -VTFS- sequence whose presence explains the unexpected ability of Cab3 (but not Hal3) to regulate Ppz1 function in Candida albicans. This sequence is conserved in a subset of fungi and it could serve to estimate the relevance of Hal3 or Cab3 proteins in regulating fungal Ppz enzymes. We also show that the removal of the motif moderately affects both Ppz1 intracellular relocalization and counteraction of toxicity in cells overexpressing the phosphatase. Thus, our work contributes to our understanding of the regulation of Ppz phosphatases, which are determinants for virulence in some pathogenic fungi.

Proteogenomics analysis of CUG codon translation in the human pathogen Candida albicans

Background

Yeasts of the CTG-clade lineage, which includes the human-infecting Candida albicans, Candida parapsilosis and Candida tropicalis species, are characterized by an altered genetic code. Instead of translating CUG codons as leucine, as happens in most eukaryotes, these yeasts, whose ancestors are thought to have lost the relevant leucine-tRNA gene, translate CUG codons as serine using a serine-tRNA with a mutated anticodon, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer. Previously reported experiments have suggested that 3–5% of the CTG-clade CUG codons are mistranslated as leucine due to mischarging of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer. The mistranslation was suggested to result in variable surface proteins explaining fast host adaptation and pathogenicity.

Results

In this study, we reassess this potential mistranslation by high-resolution mass spectrometry-based proteogenomics of multiple CTG-clade yeasts, including various C. albicans strains, isolated from colonized and from infected human body sites, and C. albicans grown in yeast and hyphal forms. Our data do not support a bias towards CUG codon mistranslation as leucine. Instead, our data suggest that (i) CUG codons are mistranslated at a frequency corresponding to the normal extent of ribosomal mistranslation with no preference for specific amino acids, (ii) CUG codons are as unambiguous (or ambiguous) as the related CUU leucine and UCC serine codons, (iii) tRNA anticodon loop variation across the CTG-clade yeasts does not result in any difference of the mistranslation level, and (iv) CUG codon unambiguity is independent of C. albicans’ strain pathogenicity or growth form.

Conclusions

Our findings imply that C. albicans does not decode CUG ambiguously. This suggests that the proposed misleucylation of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathrm{tRNA}}_{\mathrm{CAG}}^{\mathrm{Ser}} $$\end{document}tRNACAGSer might be as prevalent as every other misacylation or mistranslation event and, if at all, be just one of many reasons causing phenotypic diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01197-9.

Survival Strategies of Pathogenic Candida Species in Human Blood Show Independent and Specific Adaptations

To ensure their survival, pathogens have to adapt immediately to new environments in their hosts, for example, during the transition from the gut to the bloodstream. Here, we investigated the basis of this adaptation in a group of fungal species which are among the most common causes of hospital-acquired infections, the Candida species. On the basis of a human whole-blood infection model, we studied which genes and processes are active over the course of an infection in both the host and four different Candida pathogens. Remarkably, we found that, while the human host response during the early phase of infection is predominantly uniform, the pathogens pursue largely individual strategies and each one regulates genes involved in largely disparate processes in the blood. Our results reveal that C. albicans, C. glabrata, C. parapsilosis, and C. tropicalis all have developed individual strategies for survival in the host. This indicates that their pathogenicity in humans has evolved several times independently and that genes which are central for survival in the host for one species may be irrelevant in another.

Only four species, Candida albicans, C. glabrata, C. parapsilosis, and C. tropicalis, together account for about 90% of all Candida bloodstream infections and are among the most common causes of invasive fungal infections of humans. However, virulence potential varies among these species, and the phylogenetic tree reveals that their pathogenicity may have emerged several times independently during evolution. We therefore tested these four species in a human whole-blood infection model to determine, via comprehensive dual-species RNA-sequencing analyses, which fungal infection strategies are conserved and which are recent evolutionary developments. The ex vivo infection progressed from initial immune cell interactions to nearly complete killing of all fungal cells. During the course of infection, we characterized important parameters of pathogen-host interactions, such as fungal survival, types of interacting immune cells, and cytokine release. On the transcriptional level, we obtained a predominantly uniform and species-independent human response governed by a strong upregulation of proinflammatory processes, which was downregulated at later time points after most of the fungal cells were killed. In stark contrast, we observed that the different fungal species pursued predominantly individual strategies and showed significantly different global transcriptome patterns. Among other findings, our functional analyses revealed that the fungal species relied on different metabolic pathways and virulence factors to survive the host-imposed stress. These data show that adaptation of Candida species as a response to the host is not a phylogenetic trait, but rather has likely evolved independently as a prerequisite to cause human infections.

Fungal evolution: diversity, taxonomy and phylogeny of the Fungi

The fungal kingdom comprises a hyperdiverse clade of heterotrophic eukaryotes characterized by the presence of a chitinous cell wall, the loss of phagotrophic capabilities and cell organizations that range from completely unicellular monopolar organisms to highly complex syncitial filaments that may form macroscopic structures. Fungi emerged as a 'Third Kingdom', embracing organisms that were outside the classical dichotomy of animals versus vegetals. The taxonomy of this group has a turbulent history that is only now starting to be settled with the advent of genomics and phylogenomics. We here review the current status of the phylogeny and taxonomy of fungi, providing an overview of the main defined groups. Based on current knowledge, nine phylum-level clades can be defined: Opisthosporidia, Chytridiomycota, Neocallimastigomycota, Blastocladiomycota, Zoopagomycota, Mucoromycota, Glomeromycota, Basidiomycota and Ascomycota. For each group, we discuss their main traits and their diversity, focusing on the evolutionary relationships among the main fungal clades. We also explore the diversity and phylogeny of several groups of uncertain affinities and the main phylogenetic and taxonomical controversies and hypotheses in the field.

Chromatin-Mediated Regulation of Genome Plasticity in Human Fungal Pathogens

Human fungal pathogens, such as Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans, are a public health problem, causing millions of infections and killing almost half a million people annually. The ability of these pathogens to colonise almost every organ in the human body and cause life-threating infections relies on their capacity to adapt and thrive in diverse hostile host-niche environments. Stress-induced genome instability is a key adaptive strategy used by human fungal pathogens as it increases genetic diversity, thereby allowing selection of genotype(s) better adapted to a new environment. Heterochromatin represses gene expression and deleterious recombination and could play a key role in modulating genome stability in response to environmental changes. However, very little is known about heterochromatin structure and function in human fungal pathogens. In this review, I use our knowledge of heterochromatin structure and function in fungal model systems as a road map to review the role of heterochromatin in regulating genome plasticity in the most common human fungal pathogens: Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans.

Rapid Genetic Code Evolution in Green Algal Mitochondrial Genomes

Genetic code deviations involving stop codons have been previously reported in mitochondrial genomes of several green plants (Viridiplantae), most notably chlorophyte algae (Chlorophyta). However, as changes in codon recognition from one amino acid to another are more difficult to infer, such changes might have gone unnoticed in particular lineages with high evolutionary rates that are otherwise prone to codon reassignments. To gain further insight into the evolution of the mitochondrial genetic code in green plants, we have conducted an in-depth study across mtDNAs from 51 green plants (32 chlorophytes and 19 streptophytes). Besides confirming known stop-to-sense reassignments, our study documents the first cases of sense-to-sense codon reassignments in Chlorophyta mtDNAs. In several Sphaeropleales, we report the decoding of AGG codons (normally arginine) as alanine, by tRNA(CCU) of various origins that carry the recognition signature for alanine tRNA synthetase. In Chromochloris, we identify tRNA variants decoding AGG as methionine and the synonymous codon CGG as leucine. Finally, we find strong evidence supporting the decoding of AUA codons (normally isoleucine) as methionine in Pycnococcus. Our results rely on a recently developed conceptual framework (CoreTracker) that predicts codon reassignments based on the disparity between DNA sequence (codons) and the derived protein sequence. These predictions are then validated by an evaluation of tRNA phylogeny, to identify the evolution of new tRNAs via gene duplication and loss, and structural modifications that lead to the assignment of new tRNA identities and a change in the genetic code.

Development of a Transformation Method for Metschnikowia borealis and other CUG-Serine Yeasts

Yeasts belonging to the Metschnikowia genus are particularly interesting for the unusual formation of only two needle-shaped ascospores during their mating cycle. Presently, the meiotic process that can lead to only two spores from a diploid zygote is poorly understood. The expression of fluorescent nuclear proteins should allow the meiotic process to be visualized in vivo; however, no large-spored species of Metschnikowia has ever been transformed. Accordingly, we aimed to develop a transformation method for Metschnikowia borealis, a particularly large-spored species of Metschnikowia, with the goal of enabling the genetic manipulations required to study biological processes in detail. Genetic analyses confirmed that M. borealis, and many other Metschnikowia species, are CUG-Ser yeasts. Codon-optimized selectable markers lacking CUG codons were used to successfully transform M. borealis by electroporation and lithium acetate, and transformants appeared to be the result of random integration. Mating experiments confirmed that transformed-strains were capable of generating large asci and undergoing recombination. Finally, random integration was used to transform an additional 21 yeast strains, and all attempts successfully generated transformants. The results provide a simple method to transform many yeasts from an array of different clades and can be used to study or develop many species for various applications.

CoreTracker: accurate codon reassignment prediction, applied to mitochondrial genomes

Motivation

Codon reassignments have been reported across all domains of life. With the increasing number of sequenced genomes, the development of systematic approaches for genetic code detection is essential for accurate downstream analyses. Three automated prediction tools exist so far: FACIL, GenDecoder and Bagheera; the last two respectively restricted to metazoan mitochondrial genomes and CUG reassignments in yeast nuclear genomes. These tools can only analyze a single genome at a time and are often not followed by a validation procedure, resulting in a high rate of false positives.

Results

We present CoreTracker, a new algorithm for the inference of sense-to-sense codon reassignments. CoreTracker identifies potential codon reassignments in a set of related genomes, then uses statistical evaluations and a random forest classifier to predict those that are the most likely to be correct. Predicted reassignments are then validated through a phylogeny-aware step that evaluates the impact of the new genetic code on the protein alignment. Handling simultaneously a set of genomes in a phylogenetic framework, allows tracing back the evolution of each reassignment, which provides information on its underlying mechanism. Applied to metazoan and yeast genomes, CoreTracker significantly outperforms existing methods on both precision and sensitivity.

Availability and implementation

CoreTracker is written in Python and available at https://github.com/UdeM-LBIT/CoreTracker.

Contact

mabrouk@iro.umontreal.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Endogenous Stochastic Decoding of the CUG Codon by Competing Ser- and Leu-tRNAs in Ascoidea asiatica

Although the “universal” genetic code is now known not to be universal, and stop codons can have multiple meanings, one regularity remains, namely that for a given sense codon there is a unique translation. Examining CUG usage in yeasts that have transferred CUG away from leucine, we here report the first example of dual coding: Ascoidea asiatica stochastically encodes CUG as both serine and leucine in approximately equal proportions. This is deleterious, as evidenced by CUG codons being rare, never at conserved serine or leucine residues, and predominantly in lowly expressed genes. Related yeasts solve the problem by loss of function of one of the two tRNAs. This dual coding is consistent with the tRNA-loss-driven codon reassignment hypothesis, and provides a unique example of a proteome that cannot be deterministically predicted.

Video Abstract

•

Ascoidea asiatica stochastically encodes CUG as leucine and serine

•

It is the only known example of a proteome with non-deterministic features

•

Stochastic encoding is caused by competing tRNALeu(CAG) and tRNASer(CAG)

•

A. asiatica copes with stochastic encoding by avoiding CUG at key positions

Synthetic genome recoding: new genetic codes for new features

Full genome recoding, or rewriting codon meaning, through chemical synthesis of entire bacterial chromosomes has become feasible in the past several years. Recoding an organism can impart new properties including non-natural amino acid incorporation, virus resistance, and biocontainment. The estimated cost of construction that includes DNA synthesis, assembly by recombination, and troubleshooting, is now comparable to costs of early stage development of drugs or other high-tech products. Here, we discuss several recently published assembly methods and provide some thoughts on the future, including how synthetic efforts might benefit from the analysis of natural recoding processes and organisms that use alternative genetic codes.

[Yeasts from the CTG clade (Candida clade): Biology, impact in human health, and biotechnological applications]

Among the subdivision of Saccharomycotina (ascomycetes budding yeasts), the CTG clade (formerly the Candida clade) includes species that display a particular genetic code. In these yeasts, the CTG codon is predominantly translated as a serine instead of a leucine residue. It is now well-known that some CTG clade species have a major impact on human and its activities. Some of them are recognized as opportunistic agents of fungal infections termed candidiasis. In addition, another series of species belonging to the CTG clade draws the attention of some research groups because they exhibit a strong potential in various areas of biotechnology such as biological control, bioremediation, but also in the production of valuable biocompounds (biofuel, vitamins, sweeteners, industrial enzymes). Here we provide an overview of recent advances concerning the biology, clinical relevance, and currently tested biotechnological applications of species of the CTG clade. Future directions for scientific research on these particular yeasts are also discussed.

A CTG Clade Candida Yeast Genetically Engineered for the Genotype-Phenotype Characterization of Azole Antifungal Resistance in Human-Pathogenic Yeasts

A strain of the opportunistic pathogenic yeast Candida lusitaniae was genetically modified for use as a cellular model for assessing by allele replacement the impact of lanosterol C14α-demethylase ERG11 mutations on azole resistance. Candida lusitaniae was chosen because it is susceptible to azole antifungals, it belongs to the CTG clade of yeast, which includes most of the Candida species pathogenic for humans, and it is haploid and easily amenable to genetic transformation and molecular modeling. In this work, allelic replacement is targeted at the ERG11 locus by the reconstitution of a functional auxotrophic marker in the 3' intergenic region of ERG11 Homologous and heterologous ERG11 alleles are expressed from the resident ERG11 promoter of C. lusitaniae, allowing accurate comparison of the phenotypic change in azole susceptibility. As a proof of concept, we successfully expressed in C. lusitaniae different ERG11 alleles, either bearing or not bearing mutations retrieved from a clinical context, from two phylogenetically distant yeasts, C. albicans and Kluyveromyces marxianusCandida lusitaniae constitutes a high-fidelity expression system, giving specific Erg11p-dependent fluconazole MICs very close to those observed with the ERG11 donor strain. This work led us to characterize the phenotypic effect of two kinds of mutation: mutation conferring decreased fluconazole susceptibility in a species-specific manner and mutation conferring fluconazole resistance in several yeast species. In particular, a missense mutation affecting amino acid K143 of Erg11p in Candida species, and the equivalent position K151 in K. marxianus, plays a critical role in fluconazole resistance.

Whole Genome Sequence of the Heterozygous Clinical Isolate Candida krusei 81-B-5

Candida krusei is a diploid, heterozygous yeast that is an opportunistic fungal pathogen in immunocompromised patients. This species also is utilized for fermenting cocoa beans during chocolate production. One major concern in the clinical setting is the innate resistance of this species to the most commonly used antifungal drug fluconazole. Here, we report a high-quality genome sequence and assembly for the first clinical isolate of C. krusei, strain 81-B-5, into 11 scaffolds generated with PacBio sequencing technology. Gene annotation and comparative analysis revealed a unique profile of transporters that could play a role in drug resistance or adaptation to different environments. In addition, we show that, while 82% of the genome is highly heterozygous, a 2.0 Mb region of the largest scaffold has undergone loss of heterozygosity. This genome will serve as a reference for further genetic studies of this pathogen.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

How tRNAs dictate nuclear codon reassignments: Only a few can capture non-cognate codons

mRNA decoding by tRNAs and tRNA charging by aminoacyl-tRNA synthetases are biochemically separated processes that nevertheless in general involve the same nucleotides. The combination of charging and decoding determines the genetic code. Codon reassignment happens when a differently charged tRNA replaces a former cognate tRNA. The recent discovery of the polyphyly of the yeast CUG sense codon reassignment challenged previous mechanistic considerations and led to the proposal of the so-called tRNA loss driven codon reassignment hypothesis. Accordingly, codon capture is caused by loss of a tRNA or by mutations in the translation termination factor, subsequent reduction of the codon frequency through reduced translation fidelity and final appearance of a new cognate tRNA. Critical for codon capture are sequence and structure of the new tRNA, which must be compatible with recognition regions of aminoacyl-tRNA synthetases. The proposed hypothesis applies to all reported nuclear and organellar codon reassignments.

Reconstructing the Backbone of the Saccharomycotina Yeast Phylogeny Using Genome-Scale Data

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.

Evolutionary genomics of yeast pathogens in the Saccharomycotina

Saccharomycotina comprises a diverse group of yeasts that includes numerous species of industrial or clinical relevance. Opportunistic pathogens within this clade are often assigned to the genus Candida but belong to phylogenetically distant lineages that also comprise non-pathogenic species. This indicates that the ability to infect humans has evolved independently several times among Saccharomycotina. Although the mechanisms of infection of the main groups of Candida pathogens are starting to be unveiled, we still lack sufficient understanding of the evolutionary paths that led to a virulent phenotype in each of the pathogenic lineages. Deciphering what genomic changes underlie the evolutionary emergence of a virulence trait will not only aid the discovery of novel virulence mechanisms but it will also provide valuable information to understand how new pathogens emerge, and what clades may pose a future danger. Here we review recent comparative genomics efforts that have revealed possible evolutionary paths to pathogenesis in different lineages, focusing on the main three agents of candidiasis worldwide: Candida albicans, C. parapsilosis and C. glabrata We will discuss what genomic traits may facilitate the emergence of virulence, and focus on two different genome evolution mechanisms able to generate drastic phenotypic changes and which have been associated to the emergence of virulence: gene family expansion and interspecies hybridization.

Comparative genomics of biotechnologically important yeasts

Ascomycete yeasts are metabolically diverse, with great potential for biotechnology. Here, we report the comparative genome analysis of 29 taxonomically and biotechnologically important yeasts, including 16 newly sequenced. We identify a genetic code change, CUG-Ala, in Pachysolen tannophilus in the clade sister to the known CUG-Ser clade. Our well-resolved yeast phylogeny shows that some traits, such as methylotrophy, are restricted to single clades, whereas others, such as l-rhamnose utilization, have patchy phylogenetic distributions. Gene clusters, with variable organization and distribution, encode many pathways of interest. Genomics can predict some biochemical traits precisely, but the genomic basis of others, such as xylose utilization, remains unresolved. Our data also provide insight into early evolution of ascomycetes. We document the loss of H3K9me2/3 heterochromatin, the origin of ascomycete mating-type switching, and panascomycete synteny at the MAT locus. These data and analyses will facilitate the engineering of efficient biosynthetic and degradative pathways and gateways for genomic manipulation.

A novel nuclear genetic code alteration in yeasts and the evolution of codon reassignment in eukaryotes

The genetic code is the cellular translation table for the conversion of nucleotide sequences into amino acid sequences. Changes to the meaning of sense codons would introduce errors into almost every translated message and are expected to be highly detrimental. However, reassignment of single or multiple codons in mitochondria and nuclear genomes, although extremely rare, demonstrates that the code can evolve. Several models for the mechanism of alteration of nuclear genetic codes have been proposed (including “codon capture,” “genome streamlining,” and “ambiguous intermediate” theories), but with little resolution. Here, we report a novel sense codon reassignment in Pachysolen tannophilus, a yeast related to the Pichiaceae. By generating proteomics data and using tRNA sequence comparisons, we show that Pachysolen translates CUG codons as alanine and not as the more usual leucine. The Pachysolen tRNA_CAG is an anticodon-mutated tRNA^Ala containing all major alanine tRNA recognition sites. The polyphyly of the CUG-decoding tRNAs in yeasts is best explained by a tRNA loss driven codon reassignment mechanism. Loss of the CUG-tRNA in the ancient yeast is followed by gradual decrease of respective codons and subsequent codon capture by tRNAs whose anticodon is not part of the aminoacyl-tRNA synthetase recognition region. Our hypothesis applies to all nuclear genetic code alterations and provides several testable predictions. We anticipate more codon reassignments to be uncovered in existing and upcoming genome projects.

The birth of a deadly yeast: tracing the evolutionary emergence of virulence traits in Candida glabrata

The yeast Candida glabrata is an opportunistic human fungal pathogen whose incidence has increased in the last two decades. Despite its name, this yeast is only distantly related to the model fungal pathogen C. albicans, and more closely related to Saccharomyces cerevisiae and other yeasts that underwent an ancient whole-genome duplication. Understanding what specific traits make C. glabrata a successful opportunistic pathogen within a clade of mostly innocuous yeasts, and how these compare to virulence traits in distant pathogens such as C. albicans is a focus of intense research. From an evolutionary perspective, uncovering how the ability to infect humans has emerged multiple, independent times in different lineages may reveal new disease mechanisms and provide us with the capacity to predict which genomic features in a clade may confer a higher potential to develop virulence against humans.

Candida glabrata is an opportunistic human pathogen; genomics analyses have revealed its evolutionary path to virulence.

The transcriptome of Candida albicans mitochondria and the evolution of organellar transcription units in yeasts

Background

Yeasts show remarkable variation in the organization of their mitochondrial genomes, yet there is little experimental data on organellar gene expression outside few model species. Candida albicans is interesting as a human pathogen, and as a representative of a clade that is distant from the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. Unlike them, it encodes seven Complex I subunits in its mtDNA. No experimental data regarding organellar expression were available prior to this study.

Methods

We used high-throughput RNA sequencing and traditional RNA biology techniques to study the mitochondrial transcriptome of C. albicans strains BWP17 and SN148.

Results

The 14 protein-coding genes, two ribosomal RNA genes, and 24 tRNA genes are expressed as eight primary polycistronic transcription units. We also found transcriptional activity in the noncoding regions, and antisense transcripts that could be a part of a regulatory mechanism. The promoter sequence is a variant of the nonanucleotide identified in other yeast mtDNAs, but some of the active promoters show significant departures from the consensus. The primary transcripts are processed by a tRNA punctuation mechanism into the monocistronic and bicistronic mature RNAs. The steady state levels of various mature transcripts exhibit large differences that are a result of posttranscriptional regulation. Transcriptome analysis allowed to precisely annotate the positions of introns in the RNL (2), COB (2) and COX1 (4) genes, as well as to refine the annotation of tRNAs and rRNAs. Comparative study of the mitochondrial genome organization in various Candida species indicates that they undergo shuffling in blocks usually containing 2–3 genes, and that their arrangement in primary transcripts is not conserved. tRNA genes with their associated promoters, as well as GC-rich sequence elements play an important role in these evolutionary events.

Conclusions

The main evolutionary force shaping the mitochondrial genomes of yeasts is the frequent recombination, constantly breaking apart and joining genes into novel primary transcription units. The mitochondrial transcription units are constantly rearranged in evolution shaping the features of gene expression, such as the presence of secondary promoter sites that are inactive, or act as “booster” promoters, simplified transcriptional regulation and reliance on posttranscriptional mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2078-z) contains supplementary material, which is available to authorized users.