Assembly of a phased diploid Candida albicans genome facilitates allele-specific measurements and provides a simple model for repeat and indel structure

CRISPR-mediated Genome Editing of the Human Fungal Pathogen Candida albicans

This method describes the efficient CRISPR mediated genome editing of the diploid human fungal pathogen Candida albicans. CRISPR-mediated genome editing in C. albicans requires Cas9, guide RNA, and repair template. A plasmid expressing a yeast codon optimized Cas9 (CaCas9) has been generated. Guide sequences directly upstream from a PAM site (NGG) are cloned into the Cas9 expression vector. A repair template is then made by primer extension in vitro. Cotransformation of the repair template and vector into C. albicans leads to genome editing. Depending on the repair template used, the investigator can introduce nucleotide changes, insertions, or deletions. As C. albicans is a diploid, mutations are made in both alleles of a gene, provided that the A and B alleles do not harbor SNPs that interfere with guide targeting or repair template incorporation. Multimember gene families can be edited in parallel if suitable conserved sequences exist in all family members. The C. albicans CRISPR system described is flanked by FRT sites and encodes flippase. Upon induction of flippase, the antibiotic marker (CaCas9) and guide RNA are removed from the genome. This allows the investigator to perform subsequent edits to the genome. C. albicans CRISPR is a powerful fungal genetic engineering tool, and minor alterations to the described protocols permit the modification of other fungal species including C. glabrata, N. castellii, and S. cerevisiae.

Biofilms and vulvovaginal candidiasis

Candida species, including C. albicans, are part of the mucosal flora of most healthy women, and inhabit the gastrointestinal and genitourinary tracts. Under favourable conditions, they can colonize the vulvovaginal mucosa, giving rise to symptomatic vulvovaginal candidiasis (VVC). The mechanism by which Candida spp. produces inflammation is unknown. Both, the blastoconidia and the pseudohyphae are capable of destroying the vaginal epithelium by direct invasion. Although the symptoms are not always related to the fungal burden, in general, VVC is associated with a greater number of yeasts and pseudohyphae. Some years ago, C. albicans was the species most frequently involved in the different forms of VVC. However, infections by different species have emerged during the last two decades producing an increase in causative species of VVC such as C. glabrata, C. parapsilosis, C. krusei and C. tropicalis. Candida species are pathogenic organisms that have two forms of development: planktonic and biofilm. A biofilm is defined as a community of microorganisms attached to a surface and encompassed by an extracellular matrix. This form of presentation gives microorganisms greater resistance to antifungal agents. This review, about Candia spp. with a special emphasis on Candida albicans discusses specific areas such as biofilm structure and development, cell morphology and biofilm formation, biofilm-associated gene expression, the cell surface and adherence, the extracellular matrix, biofilm metabolism, and biofilm drug resistance in vulvovaginitis biofilms as an important virulence factor in fungi.

Gene Essentiality Analyzed by In Vivo Transposon Mutagenesis and Machine Learning in a Stable Haploid Isolate of Candida albicans

Comprehensive understanding of an organism requires that we understand the contributions of most, if not all, of its genes. Classical genetic approaches to this issue have involved systematic deletion of each gene in the genome, with comprehensive sets of mutants available only for very-well-studied model organisms. We took a different approach, harnessing the power of in vivo transposition coupled with deep sequencing to identify >500,000 different mutations, one per cell, in the prevalent human fungal pathogen Candida albicans and to map their positions across the genome. The transposition approach is efficient and less labor-intensive than classic approaches. Here, we describe the production and analysis (aided by machine learning) of a large collection of mutants and the comprehensive identification of 1,610 C. albicans genes that are essential for growth under standard laboratory conditions. Among these C. albicans essential genes, we identify those that are also essential in two distantly related model yeasts as well as those that are conserved in all four major human fungal pathogens and that are not conserved in the human genome. This list of genes with functions important for the survival of the pathogen provides a good starting point for the development of new antifungal drugs, which are greatly needed because of the emergence of fungal pathogens with elevated resistance and/or tolerance of the currently limited set of available antifungal drugs.

Knowing the full set of essential genes for a given organism provides important information about ways to promote, and to limit, its growth and survival. For many non-model organisms, the lack of a stable haploid state and low transformation efficiencies impede the use of conventional approaches to generate a genome-wide comprehensive set of mutant strains and the identification of the genes essential for growth. Here we report on the isolation and utilization of a highly stable haploid derivative of the human pathogenic fungus Candida albicans, together with a modified heterologous transposon and machine learning (ML) analysis method, to predict the degree to which all of the open reading frames are required for growth under standard laboratory conditions. We identified 1,610 C. albicans essential genes, including 1,195 with high “essentiality confidence” scores, thereby increasing the number of essential genes (currently 66 in the Candida Genome Database) by >20-fold and providing an unbiased approach to determine the degree of confidence in the determination of essentiality. Among the genes essential in C. albicans were 602 genes also essential in the model budding and fission yeasts analyzed by both deletion and transposon mutagenesis. We also identified essential genes conserved among the four major human pathogens C. albicans, Aspergillus fumigatus, Cryptococcus neoformans, and Histoplasma capsulatum and highlight those that lack homologs in humans and that thus could serve as potential targets for the design of antifungal therapies.

The Genome of the Human Pathogen Candida albicans Is Shaped by Mutation and Cryptic Sexual Recombination

The opportunistic fungal pathogen Candida albicans lacks a conventional sexual program and is thought to evolve, at least primarily, through the clonal acquisition of genetic changes. Here, we performed an analysis of heterozygous diploid genomes from 21 clinical isolates to determine the natural evolutionary processes acting on the C. albicans genome. Mutation and recombination shaped the genomic landscape among the C. albicans isolates. Strain-specific single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) clustered across the genome. Additionally, loss-of-heterozygosity (LOH) events contributed substantially to genotypic variation, with most long-tract LOH events extending to the ends of the chromosomes suggestive of repair via break-induced replication. Consistent with a model of inheritance by descent, most polymorphisms were shared between closely related strains. However, some isolates contained highly mosaic genomes consistent with strains having experienced interclade recombination during their evolutionary history. A detailed examination of mitochondrial genomes also revealed clear examples of interclade recombination among sequenced strains. These analyses therefore establish that both (para)sexual recombination and mitotic mutational processes drive evolution of this important pathogen. To further facilitate the study of C. albicans genomes, we also introduce an online platform, SNPMap, to examine SNP patterns in sequenced isolates.IMPORTANCE Mutations introduce variation into the genome upon which selection can act. Defining the nature of these changes is critical for determining species evolution, as well as for understanding the genetic changes driving important cellular processes. The heterozygous diploid fungus Candida albicans is both a frequent commensal organism and a prevalent opportunistic pathogen. A prevailing theory is that C. albicans evolves primarily through the gradual buildup of mitotic mutations, and a pressing issue is whether sexual or parasexual processes also operate within natural populations. Here, we establish that the C. albicans genome evolves by a combination of localized mutation and both short-tract and long-tract loss-of-heterozygosity (LOH) events within the sequenced isolates. Mutations are more prevalent within noncoding and heterozygous regions and LOH increases towards chromosome ends. Furthermore, we provide evidence for genetic exchange between isolates, establishing that sexual or parasexual processes have contributed to the diversity of both nuclear and mitochondrial genomes.

Promising Antifungal Targets Against Candida albicans Based on Ion Homeostasis

In recent decades, invasive fungal infections have been increasing significantly, contributing to high incidences and mortality in immunosuppressed patients. Candida albicans (C. albicans) is the most prevalent opportunistic fungal pathogen in humans that can cause severe and often fatal bloodstream infections. Current antifungal agents have several limitations, including that only a small number of classes of antifungals are available, certain of which have severe toxicity and high cost. Moreover, the emergence of drug resistance is a new limitation to successful patient outcomes. Therefore, the development of antifungals with novel targets is an essential strategy for the efficient management of C. albicans infections. It is widely recognized that ion homeostasis is crucial for all living cells. Many studies have identified that ion-signaling and transduction networks are central to fungal survival by regulating gene expression, morphological transition, host invasion, stress response, and drug resistance. Dysregulation of ion homeostasis rapidly mediates cell death, forming the mechanistic basis of a growing number of compounds that elicit antifungal activity. Most of the potent antifungals have been widely used in the clinic, and certain of them have low toxicity, meaning that they may be expected to be used as antifungal drugs in the future. Hence, we briefly summarize the homeostasis regulation of several important ions, potential antifungal targets based on these ion-signaling networks, and antifungal compounds based on the disruption of ion homeostasis. This summary will help in designing effective drugs and identifying new targets for combating fungal diseases.

Global analysis of mutations driving microevolution of a heterozygous diploid fungal pathogen

Candida albicans is a heterozygous diploid yeast that is a commensal of the human gastrointestinal tract and a prevalent opportunistic pathogen. Here, whole-genome sequencing was performed on multiple C. albicans isolates passaged both in vitro and in vivo to characterize the complete spectrum of mutations arising in laboratory culture and in the mammalian host. We establish that, independent of culture niche, microevolution is primarily driven by de novo base substitutions and frequent short-tract loss-of-heterozygosity events. An average base-substitution rate of ∼1.2 × 10^-10 per base pair per generation was observed in vitro, with higher rates inferred during host infection. Large-scale chromosomal changes were relatively rare, although chromosome 7 trisomies frequently emerged during passaging in a gastrointestinal model and was associated with increased fitness for this niche. Multiple chromosomal features impacted mutational patterns, with mutation rates elevated in repetitive regions, subtelomeric regions, and in gene families encoding cell surface proteins involved in host adhesion. Strikingly, de novo mutation rates were more than 800-fold higher in regions immediately adjacent to emergent loss-of-heterozygosity tracts, indicative of recombination-induced mutagenesis. Furthermore, genomes showed biased patterns of mutations suggestive of extensive purifying selection during passaging. These results reveal how both cell-intrinsic and cell-extrinsic factors influence C. albicans microevolution, and provide a quantitative picture of genome dynamics in this heterozygous diploid species.

A High-Throughput Candida albicans Two-Hybrid System

Candida albicans is a human fungal pathogen that does not follow the universal codon usage, as it translates the CUG codon into serine rather than leucine. This makes it difficult to study protein-protein interactions using the standard yeast two-hybrid (Y2H) system in the model organism Saccharomyces cerevisiae Due to the lack of adapted tools, only a small number of protein-protein interactions (PPIs) have been detected or studied using C. albicans-optimized tools despite the importance of PPIs to understand cell biology. However, with the sequencing of the whole genome of C. albicans, the availability of an ORFeome collection containing 5,099 open reading frames (ORFs) in Gateway-adapted donor vectors, and the creation of a Gateway-compatible C. albicans-specific two-hybrid (C2H) system, it became possible to study protein-protein interactions on a larger scale using C. albicans itself as the model organism. Erroneous translations are hereby eliminated compared to using the S. cerevisiae Y2H system. Here, we describe the technical adaptations and the first application of the C2H system for a high-throughput screen, thus making it possible to screen thousands of PPIs at once in C. albicans itself. This first, small-scale high-throughput screen, using Pho85 as a bait protein against 1,646 random prey proteins, yielded one interacting partner (Pcl5). The interaction found with the high-throughput setup was further confirmed with a low-throughput C2H experiment and with a coimmunoprecipitation (co-IP) experiment.IMPORTANCECandida albicans is a major fungal pathogen, and due to the rise of fungal infections and emerging resistance to the limited antifungals available, it is important to develop novel and more specific antifungals. Protein-protein interactions (PPIs) can be applied as very specific drug targets. However, because of the aberrant codon usage of C. albicans, the traditional yeast two-hybrid system in Saccharomyces cerevisiae is difficult to use, and only a limited number of PPIs have been described in C. albicans To overcome this, a C. albicans two-hybrid (C2H) system was developed in 2010. The current work describes, for the first time, the application of the C2H system in a high-throughput setup. We hereby show the usefulness of the C2H system to investigate and detect PPIs in C. albicans, making it possible to further elucidate protein networks in C. albicans, which has the potential to lead to the development of novel antifungals which specifically disrupt PPIs important for virulence.

Generating genomic platforms to study Candida albicans pathogenesis

The advent of the genomic era has made elucidating gene function on a large scale a pressing challenge. ORFeome collections, whereby almost all ORFs of a given species are cloned and can be subsequently leveraged in multiple functional genomic approaches, represent valuable resources toward this endeavor. Here we provide novel, genome-scale tools for the study of Candida albicans, a commensal yeast that is also responsible for frequent superficial and disseminated infections in humans. We have generated an ORFeome collection composed of 5099 ORFs cloned in a Gateway™ donor vector, representing 83% of the currently annotated coding sequences of C. albicans. Sequencing data of the cloned ORFs are available in the CandidaOrfDB database at http://candidaorfeome.eu. We also engineered 49 expression vectors with a choice of promoters, tags and selection markers and demonstrated their applicability to the study of target ORFs transferred from the C. albicans ORFeome. In addition, the use of the ORFeome in the detection of protein-protein interaction was demonstrated. Mating-compatible strains as well as Gateway™-compatible two-hybrid vectors were engineered, validated and used in a proof of concept experiment. These unique and valuable resources should greatly facilitate future functional studies in C. albicans and the elucidation of mechanisms that underlie its pathogenicity.

The type VI secretion system deploys antifungal effectors against microbial competitors

Interactions between bacterial and fungal cells shape many polymicrobial communities. Bacteria elaborate diverse strategies to interact and compete with other organisms, including the deployment of protein secretion systems. The type VI secretion system (T6SS) delivers toxic effector proteins into host eukaryotic cells and competitor bacterial cells, but, surprisingly, T6SS-delivered effectors targeting fungal cells have not been reported. Here we show that the 'antibacterial' T6SS of Serratia marcescens can act against fungal cells, including pathogenic Candida species, and identify the previously undescribed effector proteins responsible. These antifungal effectors, Tfe1 and Tfe2, have distinct impacts on the target cell, but both can ultimately cause fungal cell death. 'In competition' proteomics analysis revealed that T6SS-mediated delivery of Tfe2 disrupts nutrient uptake and amino acid metabolism in fungal cells, and leads to the induction of autophagy. Intoxication by Tfe1, in contrast, causes a loss of plasma membrane potential. Our findings extend the repertoire of the T6SS and suggest that antifungal T6SSs represent widespread and important determinants of the outcome of bacterial-fungal interactions.

Restriction digest screening facilitates efficient detection of site-directed mutations introduced by CRISPR in C. albicans UME6

Introduction of point mutations to a gene of interest is a powerful tool when determining protein function. CRISPR-mediated genome editing allows for more efficient transfer of a desired mutation into a wide range of model organisms. Traditionally, PCR amplification and DNA sequencing is used to determine if isolates contain the intended mutation. However, mutation efficiency is highly variable, potentially making sequencing costly and time consuming. To more efficiently screen for correct transformants, we have identified restriction enzymes sites that encode for two identical amino acids or one or two stop codons. We used CRISPR to introduce these restriction sites directly upstream of the Candida albicans UME6 Zn²⁺-binding domain, a known regulator of C. albicans filamentation. While repair templates coding for different restriction sites were not equally successful at introducing mutations, restriction digest screening enabled us to rapidly identify isolates with the intended mutation in a cost-efficient manner. In addition, mutated isolates have clear defects in filamentation and virulence compared to wild type C. albicans. Our data suggest restriction digestion screening efficiently identifies point mutations introduced by CRISPR and streamlines the process of identifying residues important for a phenotype of interest.

Gene flow contributes to diversification of the major fungal pathogen Candida albicans

Elucidating population structure and levels of genetic diversity and recombination is necessary to understand the evolution and adaptation of species. Candida albicans is the second most frequent agent of human fungal infections worldwide, causing high-mortality rates. Here we present the genomic sequences of 182 C. albicans isolates collected worldwide, including commensal isolates, as well as ones responsible for superficial and invasive infections, constituting the largest dataset to date for this major fungal pathogen. Although, C. albicans shows a predominantly clonal population structure, we find evidence of gene flow between previously known and newly identified genetic clusters, supporting the occurrence of (para)sexuality in nature. A highly clonal lineage, which experimentally shows reduced fitness, has undergone pseudogenization in genes required for virulence and morphogenesis, which may explain its niche restriction. Candida albicans thus takes advantage of both clonality and gene flow to diversify.

The fungal pathogen Candida albicans can undergo a parasexual process that may contribute to genetic diversity, but its actual relevance is unclear. Here, Ropars et al. analyse the genomic sequences of 182 C. albicans isolates collected worldwide and find evidence of gene flow and thus parasexuality in nature.

The Long History of the Diverse Roles of Short ORFs: sPEPs in Fungi

Since the completion of the genome sequence of the model eukaryote Saccharomyces cerevisiae, there have been significant advancements in the field of genome annotation, in no small part due to the availability of datasets that make large-scale comparative analyses possible. As a result, since its completion there has been a significant change in annotated ORF size distribution in this first eukaryotic genome, especially in short ORFs (sORFs) predicted to encode polypeptides less than 150 amino acids in length. Due to their small size and the difficulties associated with their study, it is only relatively recently that these genomic features and the sORF-encoded peptides (sPEPs) they encode have become a focus of many researchers. Yet while this class of peptides may seem new and exciting, the study of this part of the proteome is nothing new in S. cerevisiae, a species where the biological importance of sPEPs has been elegantly illustrated over the past 30 years. Here the authors showcase a range of different sORFs found in S. cerevisiae and the diverse biological roles of their encoded sPEPs, and provide an insight into the sORFs found in other fungal species, particularly those pathogenic to humans.

Rapid Phenotypic and Genotypic Diversification After Exposure to the Oral Host Niche in Candida albicans

In vitro studies suggest that stress may generate random standing variation and that different cellular and ploidy states may evolve more rapidly under stress. Yet this idea has not been tested with pathogenic fungi growing within their host niche in vivo Here, we analyzed the generation of both genotypic and phenotypic diversity during exposure of Candida albicans to the mouse oral cavity. Ploidy, aneuploidy, loss of heterozygosity (LOH), and recombination were determined using flow cytometry and double digest restriction site-associated DNA sequencing. Colony phenotypic changes in size and filamentous growth were evident without selection and were enriched among colonies selected for LOH of the GAL1 marker. Aneuploidy and LOH occurred on all chromosomes (Chrs), with aneuploidy more frequent for smaller Chrs and whole Chr LOH more frequent for larger Chrs. Large genome shifts in ploidy to haploidy often maintained one or more heterozygous disomic Chrs, consistent with random Chr missegregation events. Most isolates displayed several different types of genomic changes, suggesting that the oral environment rapidly generates diversity de novo In sharp contrast, following in vitro propagation, isolates were not enriched for multiple LOH events, except in those that underwent haploidization and/or had high levels of Chr loss. The frequency of events was overall 100 times higher for C. albicans populations following in vivo passage compared with in vitro These hyper-diverse in vivo isolates likely provide C. albicans with the ability to adapt rapidly to the diversity of stress environments it encounters inside the host.

New CRISPR Mutagenesis Strategies Reveal Variation in Repair Mechanisms among Fungi

We have created new vectors for clustered regularly interspaced short palindromic repeat (CRISPR) mutagenesis in Candida albicans, Saccharomyces cerevisiae, Candida glabrata, and Naumovozyma castellii These new vectors permit a comparison of the requirements for CRISPR mutagenesis in each of these species and reveal different dependencies for repair of the Cas9 double-stranded break. Both C. albicans and S. cerevisiae rely heavily on homology-directed repair, whereas C. glabrata and N. castellii use both homology-directed and nonhomologous end-joining pathways. The high efficiency of these vectors permits the creation of unmarked deletions in each of these species and the recycling of the dominant selection marker for serial mutagenesis in prototrophs. A further refinement, represented by the "Unified" Solo vectors, incorporates Cas9, guide RNA, and repair template into a single vector, thus enabling the creation of vector libraries for pooled screens. To facilitate the design of such libraries, we have identified guide sequences for each of these species with updated guide selection algorithms.IMPORTANCE CRISPR-mediated genome engineering technologies have revolutionized genetic studies in a wide range of organisms. Here we describe new vectors and guide sequences for CRISPR mutagenesis in the important human fungal pathogens C. albicans and C. glabrata, as well as in the related yeasts S. cerevisiae and N. castellii The design of these vectors enables efficient serial mutagenesis in each of these species by leaving few, if any, exogenous sequences in the genome. In addition, we describe strategies for the creation of unmarked deletions in each of these species and vector designs that permit the creation of vector libraries for pooled screens. These tools and strategies promise to advance genetic engineering of these medically and industrially important species.

Whole-Genome Sequences and Annotation of the Opportunistic Pathogen Candida albicans Strain SC5314 Grown under Two Different Environmental Conditions

The genetic variability of the opportunistic pathogen Candida albicans is an important adaptive mechanism. Here, we present the whole-genome sequences of the C. albicans SC5314 strain under two different growth conditions, providing useful information for comparative genomic studies and further intraspecific analysis.

Accessibility and contribution to glucan masking of natural and genetically tagged versions of yeast wall protein 1 of Candida albicans

Yeast wall protein 1 (Ywp1) is an abundant glycoprotein of the cell wall of the yeast form of Candida albicans, the most prevalent fungal pathogen of humans. Antibodies that bind to the polypeptide backbone of isolated Ywp1 show little binding to intact yeast cells, presumably because the Ywp1 epitopes are masked by the polysaccharides of the mannoproteins that form the outer layer of the cell wall. Rare cells do exhibit much greater anti-Ywp1 binding, however, and one of these was isolated and characterized. No differences were seen in its Ywp1, but it exhibited greater adhesiveness, sensitivity to wall perturbing agents, and exposure of its underlying β-1,3-glucan layer to external antibodies. The molecular basis for this greater epitope accessibility has not been determined, but has facilitated exploration of how these properties change as a function of cell growth and morphology. In addition, previously engineered strains with reduced quantities of Ywp1 in their cell walls were also found to have greater β-1,3-glucan exposure, indicating that Ywp1 itself contributes to the masking of wall epitopes, which may be important for understanding the anti-adhesive effect of Ywp1. Ectopic production of Ywp1 by hyphae, which reduces the adhesivity of these filamentous forms of C. albicans, was similarly found to reduce exposure of the β-1,3-glucan in their walls. To monitor Ywp1 in the cell wall irrespective of its accessibility, green fluorescent protein (Gfp) was genetically inserted into wall-anchored Ywp1 using a bifunctional cassette that also allowed production from a single transfection of a soluble, anchor-free version. The wall-anchored Ywp1-Gfp-Ywp1 accumulated in the wall of the yeast forms but not hyphae, and appeared to have properties similar to native Ywp1, including its adhesion-inhibiting effect. Some pseudohyphal walls also detectably accumulated this probe. Strains of C. albicans with tandem hemagglutinin (HA) epitopes inserted into wall-anchored Ywp1 were previously created by others, and were further explored here. As above, rare cells with much greater accessibility of the HA epitopes were isolated, and also found to exhibit greater exposure of Ywp1 and β-1,3-glucan. The placement of the HA cassette inhibited the normal N-glycosylation and propeptide cleavage of Ywp1, but the wall-anchored Ywp1-HA-Ywp1 still accumulated in the cell wall of yeast forms. Bifunctional transformation cassettes were used to additionally tag these molecules with Gfp, generating soluble Ywp1-HA-Gfp and wall-anchored Ywp1-HA-Gfp-Ywp1 molecules. The former revealed unexpected electrophoretic properties caused by the HA insertion, while the latter further highlighted differences between the presence of a tagged Ywp1 molecule (as revealed by Gfp fluorescence) and its accessibility in the cell wall to externally applied antibodies specific for HA, Gfp and Ywp1, with accessibility being greatest in the rapidly expanding walls of budding daughter cells. These strains and results increase our understanding of cell wall properties and how C. albicans masks itself from recognition by the human immune system.

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.

Harnessing Whole Genome Sequencing in Medical Mycology

Purpose of Review

Comparative genome sequencing studies of human fungal pathogens enable identification of genes and variants associated with virulence and drug resistance. This review describes current approaches, resources, and advances in applying whole genome sequencing to study clinically important fungal pathogens.

Recent Findings

Genomes for some important fungal pathogens were only recently assembled, revealing gene family expansions in many species and extreme gene loss in one obligate species. The scale and scope of species sequenced is rapidly expanding, leveraging technological advances to assemble and annotate genomes with higher precision. By using iteratively improved reference assemblies or those generated de novo for new species, recent studies have compared the sequence of isolates representing populations or clinical cohorts. Whole genome approaches provide the resolution necessary for comparison of closely related isolates, for example, in the analysis of outbreaks or sampled across time within a single host.

Summary

Genomic analysis of fungal pathogens has enabled both basic research and diagnostic studies. The increased scale of sequencing can be applied across populations, and new metagenomic methods allow direct analysis of complex samples.

Haplotyping a Non-meiotic Diploid Fungal Pathogen Using Induced Aneuploidies and SNP/CGH Microarray Analysis

The generation of haplotype information has recently become very attractive due to its utility for identifying mutations associated with human disease and for the development of personalized medicine. Haplotype information also is crucial for studying recombination mechanisms and genetic diversity, and for analyzing allele-specific gene expression. Classic haplotyping methods require the analysis of hundreds of meiotic progeny. To facilitate haplotyping in the non-meiotic human fungal pathogen Candida albicans, we exploited trisomic heterozygous chromosomes generated via the UAU1 selection strategy. Using this system, we obtained phasing information from allelic biases, detected by SNP/CGH microarray analysis. This strategy has the potential to be applicable to other diploid, asexual Candida species that are important causes of human disease.

The Candida Genome Database (CGD): incorporation of Assembly 22, systematic identifiers and visualization of high throughput sequencing data

The Candida Genome Database (CGD, http://www.candidagenome.org/) is a freely available online resource that provides gene, protein and sequence information for multiple Candida species, along with web-based tools for accessing, analyzing and exploring these data. The mission of CGD is to facilitate and accelerate research into Candida pathogenesis and biology, by curating the scientific literature in real time, and connecting literature-derived annotations to the latest version of the genomic sequence and its annotations. Here, we report the incorporation into CGD of Assembly 22, the first chromosome-level, phased diploid assembly of the C. albicans genome, coupled with improvements that we have made to the assembly using additional available sequence data. We also report the creation of systematic identifiers for C. albicans genes and sequence features using a system similar to that adopted by the yeast community over two decades ago. Finally, we describe the incorporation of JBrowse into CGD, which allows online browsing of mapped high throughput sequencing data, and its implementation for several RNA-Seq data sets, as well as the whole genome sequencing data that was used in the construction of Assembly 22.

Analysis of Repair Mechanisms following an Induced Double-Strand Break Uncovers Recessive Deleterious Alleles in the Candida albicans Diploid Genome

The diploid genome of the yeast Candida albicans is highly plastic, exhibiting frequent loss-of-heterozygosity (LOH) events. To provide a deeper understanding of the mechanisms leading to LOH, we investigated the repair of a unique DNA double-strand break (DSB) in the laboratory C. albicans SC5314 strain using the I-SceI meganuclease. Upon I-SceI induction, we detected a strong increase in the frequency of LOH events at an I-SceI target locus positioned on chromosome 4 (Chr4), including events spreading from this locus to the proximal telomere. Characterization of the repair events by single nucleotide polymorphism (SNP) typing and whole-genome sequencing revealed a predominance of gene conversions, but we also observed mitotic crossover or break-induced replication events, as well as combinations of independent events. Importantly, progeny that had undergone homozygosis of part or all of Chr4 haplotype B (Chr4B) were inviable. Mining of genome sequencing data for 155 C. albicans isolates allowed the identification of a recessive lethal allele in the GPI16 gene on Chr4B unique to C. albicans strain SC5314 which is responsible for this inviability. Additional recessive lethal or deleterious alleles were identified in the genomes of strain SC5314 and two clinical isolates. Our results demonstrate that recessive lethal alleles in the genomes of C. albicans isolates prevent the occurrence of specific extended LOH events. While these and other recessive lethal and deleterious alleles are likely to accumulate in C. albicans due to clonal reproduction, their occurrence may in turn promote the maintenance of corresponding nondeleterious alleles and, consequently, heterozygosity in the C. albicans species.

IMPORTANCE

Recessive lethal alleles impose significant constraints on the biology of diploid organisms. Using a combination of an I-SceI meganuclease-mediated DNA DSB, a fluorescence-activated cell sorter (FACS)-optimized reporter of LOH, and a compendium of 155 genome sequences, we were able to unmask and identify recessive lethal and deleterious alleles in isolates of Candida albicans, a diploid yeast and the major fungal pathogen of humans. Accumulation of recessive deleterious mutations upon clonal reproduction of C. albicans could contribute to the maintenance of heterozygosity despite the high frequency of LOH events in this species.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum

The malaria parasite Plasmodium falciparum has a great capacity for evolutionary adaptation to evade host immunity and develop drug resistance. Current understanding of parasite evolution is impeded by the fact that a large fraction of the genome is either highly repetitive or highly variable and thus difficult to analyze using short-read sequencing technologies. Here, we describe a resource of deep sequencing data on parents and progeny from genetic crosses, which has enabled us to perform the first genome-wide, integrated analysis of SNP, indel and complex polymorphisms, using Mendelian error rates as an indicator of genotypic accuracy. These data reveal that indels are exceptionally abundant, being more common than SNPs and thus the dominant mode of polymorphism within the core genome. We use the high density of SNP and indel markers to analyze patterns of meiotic recombination, confirming a high rate of crossover events and providing the first estimates for the rate of non-crossover events and the length of conversion tracts. We observe several instances of meiotic recombination within copy number variants associated with drug resistance, demonstrating a mechanism whereby fitness costs associated with resistance mutations could be compensated and greater phenotypic plasticity could be acquired.

Phenotypic Consequences of a Spontaneous Loss of Heterozygosity in a Common Laboratory Strain of Candida albicans

By testing the susceptibility to DNA damaging agents of several Candida albicans mutant strains derived from the commonly used laboratory strain, CAI4, we uncovered sensitivity to methyl methanesulfonate (MMS) in CAI4 and its derivatives, but not in CAF2-1. This sensitivity is not a result of URA3 disruption because the phenotype was not restored after URA3 reintroduction. Rather, we found that homozygosis of a short region of chromosome 3R (Chr3R), which is naturally heterozygous in the MMS-resistant-related strains CAF4-2 and CAF2-1, confers MMS sensitivity and modulates growth polarization in response to MMS. Furthermore, induction of homozygosity in this region in CAF2-1 or CAF4-2 resulted in MMS sensitivity. We identified 11 genes by SNP/comparative genomic hybridization containing only the a alleles in all the MMS-sensitive strains. Four candidate genes, SNF5, POL1, orf19.5854.1, and MBP1, were analyzed by generating hemizygous configurations in CAF2-1 and CAF4-2 for each allele of all four genes. Only hemizygous MBP1a/mbp1b::SAT1-FLIP strains became MMS sensitive, indicating that MBP1a in the homo- or hemizygosis state was sufficient to account for the MMS-sensitive phenotype. In yeast, Mbp1 regulates G1/S genes involved in DNA repair. A second region of homozygosis on Chr2L increased MMS sensitivity in CAI4 (Chr3R homozygous) but not CAF4-2 (Chr3R heterozygous). This is the first example of sign epistasis in C. albicans.

How to Use the Candida Genome Database

Studying Candida biology requires access to genomic sequence data in conjunction with experimental information that provides functional context to genes and proteins. The Candida Genome Database (CGD) integrates functional information about Candida genes and their products with a set of analysis tools that facilitate searching for sets of genes and exploring their biological roles. This chapter describes how the various types of information available at CGD can be searched, retrieved, and analyzed. Starting with the guided tour of the CGD Home page and Locus Summary page, this unit shows how to navigate the various assemblies of the C. albicans genome, how to use Gene Ontology tools to make sense of large-scale data, and how to access the microarray data archived at CGD.

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.

A comprehensive Candida albicans PeptideAtlas build enables deep proteome coverage

To provide new and expanded proteome documentation of the opportunistically pathogen Candida albicans, we have developed new protein extraction and analysis routines to provide a new, extended and enhanced version of the C. albicans PeptideAtlas. Two new datasets, resulting from experiments consisting of exhaustive subcellular fractionations and different growing conditions, plus two additional datasets from previous experiments on the surface and the secreted proteomes, have been incorporated to increase the coverage of the proteome. High resolution precursor mass spectrometry (MS) and ion trap tandem MS spectra were analyzed with three different search engines using a database containing allele-specific sequences. This approach, novel for a large-scale C. albicans proteomics project, was combined with the post-processing and filtering implemented in the Trans Proteomic Pipeline consistently used in the PeptideAtlas project and resulted in 49,372 additional peptides (3-fold increase) and 1630 more proteins (1.6-fold increase) identified in the new C. albicans PeptideAtlas with respect to the previous build. A total of 71,310 peptides and 4174 canonical (minimal non-redundant set) proteins (4115 if one protein per pair of alleles is considered) were identified representing 66% of the 6218 proteins in the predicted proteome. This makes the new PeptideAtlas build the most comprehensive C. albicans proteomics resource available and the only large-scale one with detections of individual alleles.

Budding off: bringing functional genomics to Candida albicans

Candida species are the most prevalent human fungal pathogens, with Candida albicans being the most clinically relevant species. Candida albicans resides as a commensal of the human gastrointestinal tract but is a frequent cause of opportunistic mucosal and systemic infections. Investigation of C. albicans virulence has traditionally relied on candidate gene approaches, but recent advances in functional genomics have now facilitated global, unbiased studies of gene function. Such studies include comparative genomics (both between and within Candida species), analysis of total RNA expression, and regulation and delineation of protein-DNA interactions. Additionally, large collections of mutant strains have begun to aid systematic screening of clinically relevant phenotypes. Here, we will highlight the development of functional genomics in C. albicans and discuss the use of these approaches to addressing both commensalism and pathogenesis in this species.

The mitochondrial protein Mcu1 plays important roles in carbon source utilization, filamentation, and virulence in Candida albicans

The fungus Candida albicans is both a pathogen and a commensal in humans. The ability to utilize different carbon sources available in diverse host niches is vital for both commensalism and pathogenicity. N-acetylglucosamine (GlcNAc) is an important signaling molecule as well as a carbon source in C. albicans. Here, we report the discovery of a novel gene MCU1 essential for GlcNAc utilization. Mcu1 is located in mitochondria and associated with multiple energy- and metabolism-related proteins including Por1, Atp1, Pet9, and Mdh1. Consistently, inactivating Por1 impaired GlcNAc utilization as well. Deletion of MCU1 also caused defects in utilizing non-fermentable carbon sources and amino acids. Furthermore, MCU1 is required for filamentation in several inducing conditions and virulence in a mouse systemic infection model. We also deleted TGL99 and GUP1, two genes adjacent to MCU1, and found that the gup1/gup1 mutant exhibited mild defects in the utilization of several carbon sources including GlcNAc, maltose, galactose, amino acids, and ethanol. Our results indicate that MCU1 exists in a cluster of genes involved in the metabolism of carbon sources. Given its importance in metabolism and lack of a homolog in humans, Mcu1 could be a potential target for developing antifungal agents.

A FACS-optimized screen identifies regulators of genome stability in Candida albicans

Loss of heterozygosity (LOH) plays important roles in genome dynamics, notably, during tumorigenesis. In the fungal pathogen Candida albicans, LOH contributes to the acquisition of antifungal resistance. In order to investigate the mechanisms that regulate LOH in C. albicans, we have established a novel method combining an artificial heterozygous locus harboring the blue fluorescent protein and green fluorescent protein markers and flow cytometry to detect LOH events at the single-cell level. Using this fluorescence-based method, we have confirmed that elevated temperature, treatment with methyl methanesulfonate, and inactivation of the Mec1 DNA damage checkpoint kinase triggered an increase in the frequency of LOH. Taking advantage of this system, we have searched for C. albicans genes whose overexpression triggered an increase in LOH and identified four candidates, some of which are known regulators of genome dynamics with human homologues contributing to cancer progression. Hence, the approach presented here will allow the implementation of new screens to identify genes that are important for genome stability in C. albicans and more generally in eukaryotic cells.

Genetic and phenotypic intra-species variation in Candida albicans

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

YMAP: a pipeline for visualization of copy number variation and loss of heterozygosity in eukaryotic pathogens

The design of effective antimicrobial therapies for serious eukaryotic pathogens requires a clear understanding of their highly variable genomes. To facilitate analysis of copy number variations, single nucleotide polymorphisms and loss of heterozygosity events in these pathogens, we developed a pipeline for analyzing diverse genome-scale datasets from microarray, deep sequencing, and restriction site associated DNA sequence experiments for clinical and laboratory strains of Candida albicans, the most prevalent human fungal pathogen. The Y_MAP pipeline (http://lovelace.cs.umn.edu/Ymap/) automatically illustrates genome-wide information in a single intuitive figure and is readily modified for the analysis of other pathogens with small genomes.

Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.

The Candida pathogenic species complex

Candida species are the most common causes of fungal infection. Approximately 90% of infections are caused by five species: Candida albicans, Candida glabrata, Candida tropicalis, Candida parapsilosis, and Candida krusei. Three (C. albicans, C. tropicalis, and C. parapsilosis) belong to the CTG clade, in which the CTG codon is translated as serine and not leucine. C. albicans remains the most commonly isolated but is decreasing relative to the other species. The increasing incidence of C. glabrata is related to its reduced susceptibility to azole drugs. Genome analysis suggests that virulence in the CTG clade is associated with expansion of gene families, particularly of cell wall genes. Similar independent processes took place in the C. glabrata species group. Gene loss and expansion in an ancestor of C. glabrata may have resulted in preadaptations that enabled pathogenicity.

Rapid mechanisms for generating genome diversity: whole ploidy shifts, aneuploidy, and loss of heterozygosity

Human fungal pathogens can exist in a variety of ploidy states, including euploid and aneuploid forms. Ploidy change has a major impact on phenotypic properties, including the regulation of interactions with the human host. In addition, the rapid emergence of drug-resistant isolates is often associated with the formation of specific supernumerary chromosomes. Pathogens such as Candida albicans and Cryptococcus neoformans appear particularly well adapted for propagation in multiple ploidy states with novel pathways driving ploidy variation. In both species, heterozygous cells also readily undergo loss of heterozygosity (LOH), leading to additional phenotypic changes such as altered drug resistance. Here, we examine the sexual and parasexual cycles that drive ploidy variation in human fungal pathogens and discuss ploidy and LOH events with respect to their far-reaching roles in fungal adaptation and pathogenesis.

Extensive and coordinated control of allele-specific expression by both transcription and translation in Candida albicans

Though sequence differences between alleles are often limited to a few polymorphisms, these differences can cause large and widespread allelic variation at the expression level. Such allele-specific expression (ASE) has been extensively explored at the level of transcription but not translation. Here we measured ASE in the diploid yeast Candida albicans at both the transcriptional and translational levels using RNA-seq and ribosome profiling, respectively. Since C. albicans is an obligate diploid, our analysis isolates ASE arising from cis elements in a natural, nonhybrid organism, where allelic effects reflect evolutionary forces. Importantly, we find that ASE arising from translation is of a similar magnitude as transcriptional ASE, both in terms of the number of genes affected and the magnitude of the bias. We further observe coordination between ASE at the levels of transcription and translation for single genes. Specifically, reinforcing relationships—where transcription and translation favor the same allele—are more frequent than expected by chance, consistent with selective pressure tuning ASE at multiple regulatory steps. Finally, we parameterize alleles based on a range of properties and find that SNP location and predicted mRNA-structure stability are associated with translational ASE in cis. Since this analysis probes more than 4000 allelic pairs spanning a broad range of variations, our data provide a genome-wide view into the relative impact of cis elements that regulate translation.

Molecular epidemiology, phylogeny and evolution of Candida albicans

A small number of Candida species form part of the normal microbial flora of mucosal surfaces in humans and may give rise to opportunistic infections when host defences are impaired. Candida albicans is by far the most prevalent commensal and pathogenic Candida species. Several different molecular typing approaches including multilocus sequence typing, multilocus microsatellite typing and DNA fingerprinting using C. albicans-specific repetitive sequence-containing DNA probes have yielded a wealth of information regarding the epidemiology and population structure of this species. Such studies revealed that the C. albicans population structure consists of multiple major and minor clades, some of which exhibit geographical or phenotypic enrichment and that C. albicans reproduction is predominantly clonal. Despite this, losses of heterozygosity by recombination, the existence of a parasexual cycle, toleration of a wide range of aneuploidies and the recent description of viable haploid strains have all demonstrated the extensive plasticity of the C. albicans genome. Recombination and gross chromosomal rearrangements are more common under stressful environmental conditions, and have played a significant role in the evolution of this opportunistic pathogen. Surprisingly, Candida dubliniensis, the closest relative of C. albicans exhibits more karyotype variability than C. albicans, but is significantly less adaptable to unfavourable environments. This disparity most likely reflects the evolutionary processes that occurred during or soon after the divergence of both species from their common ancestor. Whilst C. dubliniensis underwent significant gene loss and pseudogenisation, C. albicans expanded gene families considered to be important in virulence. It is likely that technological developments in whole genome sequencing and data analysis in coming years will facilitate its routine use for population structure, epidemiological investigations, and phylogenetic analyses of Candida species. These are likely to reveal more minor C. albicans clades and to enhance our understanding of the population biology of this versatile organism.