Step-wise evolution of azole resistance through copy number variation followed by KSR1 loss of heterozygosity in Candida albicans.

Antimicrobial drug resistance poses a global health threat, requiring a deeper understanding of the evolutionary processes that lead to its emergence in pathogens. Complex evolutionary dynamics involve multiple mutations that can result in cooperative or competitive (clonal interference) effects. Candida albicans, a major fungal pathogen, displays high rates of copy number variation (CNV) and loss of heterozygosity (LOH). CNV and LOH events involve large numbers of genes and could synergize during evolutionary adaptation. Understanding the contributions of CNV and LOH to antifungal drug adaptation is challenging, especially in the context of whole-population genome sequencing. Here, we document the sequential evolution of fluconazole tolerance and then resistance in a C. albicans isolate involving an initial CNV on chromosome 4, followed by an LOH on chromosome R that involves KSR1. Similar LOH events involving KSR1, which encodes a reductase in the sphingolipid biosynthesis pathway, were also detected in independently evolved fluconazole resistant isolates. We dissect the specific KSR1 codons that affect fluconazole resistance and tolerance. The combination of the chromosome 4 CNV and KSR1 LOH results in a >500-fold decrease in azole susceptibility relative to the progenitor, illustrating a compelling example of rapid, yet step-wise, interplay between CNV and LOH in drug resistance evolution.

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1.

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

Variations in candidalysin amino acid sequence influence toxicity and host responses.

Candida albicans causes millions of mucosal infections in humans annually. Hyphal overgrowth on mucosal surfaces is frequently associated with tissue damage caused by candidalysin, a secreted peptide toxin that destabilizes the plasma membrane of host cells thereby promoting disease and immunopathology. Candidalysin was first identified in C. albicans strain SC5314, but recent investigations have revealed candidalysin "variants" of differing amino acid sequence in isolates of C. albicans, and the related species C. dubliniensis, and C tropicalis, suggesting that sequence variation among candidalysins may be widespread in natural populations of these Candida species. Here, we analyzed ECE1 gene sequences from 182 C. albicans isolates, 10 C. dubliniensis isolates, and 78 C. tropicalis isolates and identified 10, 3, and 2 candidalysin variants in these species, respectively. Application of candidalysin variants to epithelial cells revealed differences in the ability to cause cellular damage, changes in metabolic activity, calcium influx, MAPK signalling, and cytokine secretion, while biophysical analyses indicated that variants exhibited differences in their ability to interact with and permeabilize a membrane. This study identifies candidalysin variants with differences in biological activity that are present in medically relevant Candida species. IMPORTANCE: Fungal infections are a significant burden to health. Candidalysin is a toxin produced by Candida albicans that damages host tissues, facilitating infection. Previously, we demonstrated that candidalysins exist in the related species C. dubliniensis and C. tropicalis, thereby identifying these molecules as a toxin family. Recent genomic analyses have highlighted the presence of a small number of candidalysin "variant" toxins, which have different amino acid sequences to those originally identified. Here, we screened genome sequences of isolates of C. albicans, C. dubliniensis, and C. tropicalis and identified candidalysin variants in all three species. When applied to epithelial cells, candidalysin variants differed in their ability to cause damage, activate intracellular signaling pathways, and induce innate immune responses, while biophysical analysis revealed differences in the ability of candidalysin variants to interact with lipid bilayers. These findings suggest that intraspecies variation in candidalysin amino acid sequence may influence fungal pathogenicity.

The spliceosome impacts morphogenesis in the human fungal pathogen Candida albicans.

At human body temperature, the fungal pathogen Candida albicans can transition from yeast to filamentous morphologies in response to host-relevant cues. Additionally, elevated temperatures encountered during febrile episodes can independently induce C. albicans filamentation. However, the underlying genetic pathways governing this developmental transition in response to elevated temperatures remain largely unexplored. Here, we conducted a functional genomic screen to unravel the genetic mechanisms orchestrating C. albicans filamentation specifically in response to elevated temperature, implicating 45% of genes associated with the spliceosome or pre-mRNA splicing in this process. Employing RNA-Seq to elucidate the relationship between mRNA splicing and filamentation, we identified greater levels of intron retention in filaments compared to yeast, which correlated with reduced expression of the affected genes. Intriguingly, homozygous deletion of a gene encoding a spliceosome component important for filamentation (PRP19) caused even greater levels of intron retention compared with wild type and displayed globally dysregulated gene expression. This suggests that intron retention is a mechanism for fine-tuning gene expression during filamentation, with perturbations of the spliceosome exacerbating this process and blocking filamentation. Overall, this study unveils a novel biological process governing C. albicans filamentation, providing new insights into the complex regulation of this key virulence trait.IMPORTANCEFungal pathogens such as Candida albicans can cause serious infections with high mortality rates in immunocompromised individuals. When C. albicans is grown at temperatures encountered during human febrile episodes, yeast cells undergo a transition to filamentous cells, and this process is key to its virulence. Here, we expanded our understanding of how C. albicans undergoes filamentation in response to elevated temperature and identified many genes involved in mRNA splicing that positively regulate filamentation. Through transcriptome analyses, we found that intron retention is a mechanism for fine-tuning gene expression in filaments, and perturbation of the spliceosome exacerbates intron retention and alters gene expression substantially, causing a block in filamentation. This work adds to the growing body of knowledge on the role of introns in fungi and provides new insights into the cellular processes that regulate a key virulence trait in C. albicans.

Metabolic reprogramming during Candida albicans planktonic-biofilm transition is modulated by the transcription factors Zcf15 and Zcf26.

Candida albicans is a commensal of the human microbiota that can form biofilms on implanted medical devices. These biofilms are tolerant to antifungals and to the host immune system. To identify novel genes modulating C. albicans biofilm formation, we performed a large-scale screen with 2,454 C. albicans doxycycline-dependent overexpression strains and identified 16 genes whose overexpression significantly hampered biofilm formation. Among those, overexpression of the ZCF15 and ZCF26 paralogs that encode transcription factors and have orthologs only in biofilm-forming species of the Candida clade, caused impaired biofilm formation both in vitro and in vivo. Interestingly, overexpression of ZCF15 impeded biofilm formation without any defect in hyphal growth. Transcript profiling, transcription factor binding, and phenotypic microarray analyses conducted upon overexpression of ZCF15 and ZCF26 demonstrated their role in reprogramming cellular metabolism by regulating central metabolism including glyoxylate and tricarboxylic acid cycle genes. Taken together, this study has identified a new set of biofilm regulators, including ZCF15 and ZCF26, that appear to control biofilm development through their specific role in metabolic remodeling.

Identification of an active RNAi pathway in Candida albicans.

RNA interference (RNAi) is a fundamental regulatory pathway with a wide range of functions, including regulation of gene expression and maintenance of genome stability. Although RNAi is widespread in the fungal kingdom, well-known species, such as the model yeast Saccharomyces cerevisiae, have lost the RNAi pathway. Until now evidence has been lacking for a fully functional RNAi pathway in Candida albicans, a human fungal pathogen considered critically important by the World Health Organization. Here, we demonstrated that the widely used C. albicans reference strain (SC5314) contains an inactivating missense mutation in the gene encoding for the central RNAi component Argonaute. In contrast, most other C. albicans isolates contain a canonical Argonaute protein predicted to be functional and RNAi-active. Indeed, using high-throughput small and long RNA sequencing combined with seamless CRISPR/Cas9-based gene editing, we demonstrate that an active C. albicans RNAi machinery represses expression of subtelomeric gene families. Thus, an intact and functional RNAi pathway exists in C. albicans, highlighting the importance of using multiple reference strains when studying this dangerous pathogen.

A gain-of-function mutation in zinc cluster transcription factor Rob1 drives Candida albicans adaptive growth in the cystic fibrosis lung environment.

Candida albicans chronically colonizes the respiratory tract of patients with Cystic Fibrosis (CF). It competes with CF-associated pathogens (e.g. Pseudomonas aeruginosa) and contributes to disease severity. We hypothesize that C. albicans undergoes specific adaptation mechanisms that explain its persistence in the CF lung environment. To identify the underlying genetic and phenotypic determinants, we serially recovered 146 C. albicans clinical isolates over a period of 30 months from the sputum of 25 antifungal-naive CF patients. Multilocus sequence typing analyses revealed that most patients were individually colonized with genetically close strains, facilitating comparative analyses between serial isolates. We strikingly observed differential ability to filament and form monospecies and dual-species biofilms with P. aeruginosa among 18 serial isolates sharing the same diploid sequence type, recovered within one year from a pediatric patient. Whole genome sequencing revealed that their genomes were highly heterozygous and similar to each other, displaying a highly clonal subpopulation structure. Data mining identified 34 non-synonymous heterozygous SNPs in 19 open reading frames differentiating the hyperfilamentous and strong biofilm-former strains from the remaining isolates. Among these, we detected a glycine-to-glutamate substitution at position 299 (G299E) in the deduced amino acid sequence of the zinc cluster transcription factor ROB1 (ROB1G299E), encoding a major regulator of filamentous growth and biofilm formation. Introduction of the G299E heterozygous mutation in a co-isolated weak biofilm-former CF strain was sufficient to confer hyperfilamentous growth, increased expression of hyphal-specific genes, increased monospecies biofilm formation and increased survival in dual-species biofilms formed with P. aeruginosa, indicating that ROB1G299E is a gain-of-function mutation. Disruption of ROB1 in a hyperfilamentous isolate carrying the ROB1G299E allele abolished hyperfilamentation and biofilm formation. Our study links a single heterozygous mutation to the ability of C. albicans to better survive during the interaction with other CF-associated microbes and illuminates how adaptive traits emerge in microbial pathogens to persistently colonize and/or infect the CF-patient airways.

Secretion of the fungal toxin candidalysin is dependent on conserved precursor peptide sequences.

The opportunistic fungal pathogen Candida albicans damages host cells via its peptide toxin, candidalysin. Before secretion, candidalysin is embedded in a precursor protein, Ece1, which consists of a signal peptide, the precursor of candidalysin and seven non-candidalysin Ece1 peptides (NCEPs), and is found to be conserved in clinical isolates. Here we show that the Ece1 polyprotein does not resemble the usual precursor structure of peptide toxins. C. albicans cells are not susceptible to their own toxin, and single NCEPs adjacent to candidalysin are sufficient to prevent host cell toxicity. Using a series of Ece1 mutants, mass spectrometry and anti-candidalysin nanobodies, we show that NCEPs play a role in intracellular Ece1 folding and candidalysin secretion. Removal of single NCEPs or modifications of peptide sequences cause an unfolded protein response (UPR), which in turn inhibits hypha formation and pathogenicity in vitro. Our data indicate that the Ece1 precursor is not required to block premature pore-forming toxicity, but rather to prevent intracellular auto-aggregation of candidalysin sequences.

A CO2 sensing module modulates ß-1,3-glucan exposure in Candida albicans.

Microbial species capable of co-existing with healthy individuals, such as the commensal fungus Candida albicans, exploit multifarious strategies to evade our immune defenses. These strategies include the masking of immunoinflammatory pathogen-associated molecular patterns (PAMPs) at their cell surface. We reported previously that C. albicans actively reduces the exposure of the proinflammatory PAMP, ß-1,3-glucan, at its cell surface in response to host-related signals such as lactate and hypoxia. Here, we show that clinical isolates of C. albicans display phenotypic variability with respect to their lactate- and hypoxia-induced ß-1,3-glucan masking. We have exploited this variability to identify responsive and non-responsive clinical isolates. We then performed RNA sequencing on these isolates to reveal genes whose expression patterns suggested potential association with lactate- or hypoxia-induced ß-1,3-glucan masking. The deletion of two such genes attenuated masking: PHO84 and NCE103. We examined NCE103-related signaling further because NCE103 has been shown previously to encode carbonic anhydrase, which promotes adenylyl cyclase-protein kinase A (PKA) signaling at low CO2 levels. We show that while CO2 does not trigger ß-1,3-glucan masking in C. albicans, the Sch9-Rca1-Nce103 signaling module strongly influences ß-1,3-glucan exposure in response to hypoxia and lactate. In addition to identifying a new regulatory module that controls PAMP exposure in C. albicans, our data imply that this module is important for PKA signaling in response to environmental inputs other than CO2.IMPORTANCEOur innate immune defenses have evolved to protect us against microbial infection in part via receptor-mediated detection of "pathogen-associated molecular patterns" (PAMPs) expressed by invading microbes, which then triggers their immune clearance. Despite this surveillance, many microbial species are able to colonize healthy, immune-competent individuals, without causing infection. To do so, these microbes must evade immunity. The commensal fungus Candida albicans exploits a variety of strategies to evade immunity, one of which involves reducing the exposure of a proinflammatory PAMP (ß-1,3-glucan) at its cell surface. Most of the ß-1,3-glucan is located in the inner layer of the C. albicans cell wall, hidden by an outer layer of mannan fibrils. Nevertheless, some ß-1,3-glucan can become exposed at the fungal cell surface. However, in response to certain specific host signals, such as lactate or hypoxia, C. albicans activates an anticipatory protective response that decreases ß-1,3-glucan exposure, thereby reducing the susceptibility of the fungus to impending innate immune attack. Here, we exploited the natural phenotypic variability of C. albicans clinical isolates to identify strains that do not display the response to ß-1,3-glucan masking signals observed for the reference isolate, SC5314. Then, using genome-wide transcriptional profiling, we compared these non-responsive isolates with responsive controls to identify genes potentially involved in ß-1,3-glucan masking. Mutational analysis of these genes revealed that a sensing module that was previously associated with CO2 sensing also modulates ß-1,3-glucan exposure in response to hypoxia and lactate in this major fungal pathogen of humans.

Loss of CorA, the primary magnesium transporter of Salmonella, is alleviated by MgtA and PhoP-dependent compensatory mechanisms.

In many Gram-negative bacteria, the stress sigma factor of RNA polymerase, σS/RpoS, remodels global gene expression to reshape the physiology of stationary phase cells and ensure their survival under non-optimal growth conditions. In the foodborne pathogen Salmonella enterica serovar Typhimurium, σS is also required for biofilm formation and virulence. We have recently shown that a ΔrpoS mutation decreases the magnesium content and expression level of the housekeeping Mg2+-transporter CorA in stationary phase Salmonella. The other two Mg2+-transporters of Salmonella are encoded by the PhoP-activated mgtA and mgtB genes and are expressed under magnesium starvation. The σS control of corA prompted us to evaluate the impact of CorA in stationary phase Salmonella cells, by using global and analytical proteomic analyses and physiological assays. The ΔcorA mutation conferred a competitive disadvantage to exit from stationary phase, and slightly impaired motility, but had no effect on total and free cellular magnesium contents. In contrast to the wild-type strain, the ΔcorA mutant produced MgtA, but not MgtB, in the presence of high extracellular magnesium concentration. Under these conditions, MgtA production in the ΔcorA mutant did not require PhoP. Consistently, a ΔmgtA, but not a ΔphoP, mutation slightly reduced the magnesium content of the ΔcorA mutant. Synthetic phenotypes were observed when the ΔphoP and ΔcorA mutations were combined, including a strong reduction in growth and motility, independently of the extracellular magnesium concentration. The abundance of several proteins involved in flagella formation, chemotaxis and secretion was lowered by the ΔcorA and ΔphoP mutations in combination, but not alone. These findings unravel the importance of PhoP-dependent functions in the absence of CorA when magnesium is sufficient. Altogether, our data pinpoint a regulatory network, where the absence of CorA is sensed by the cell and compensated by MgtA and PhoP- dependent mechanisms.

Coregulation of extracellular vesicle production and fluconazole susceptibility in Cryptococcus neoformans.

Resistance to fluconazole (FLC), the most widely used antifungal drug, is typically achieved by altering the azole drug target and/or drug efflux pumps. Recent reports have suggested a link between vesicular trafficking and antifungal resistance. Here, we identified novel Cryptococcus neoformans regulators of extracellular vesicle (EV) biogenesis that impact FLC resistance. In particular, the transcription factor Hap2 does not affect the expression of the drug target or efflux pumps, yet it impacts the cellular sterol profile. Subinhibitory FLC concentrations also downregulate EV production. Moreover, in vitro spontaneous FLC-resistant colonies showed altered EV production, and the acquisition of FLC resistance was associated with decreased EV production in clinical isolates. Finally, the reversion of FLC resistance was associated with increased EV production. These data suggest a model in which fungal cells can regulate EV production in place of regulating the drug target gene expression as a first line of defense against antifungal assault in this fungal pathogen. IMPORTANCE Extracellular vesicles (EVs) are membrane-enveloped particles that are released by cells into the extracellular space. Fungal EVs can mediate community interactions and biofilm formation, but their functions remain poorly understood. Here, we report the identification of the first regulators of EV production in the major fungal pathogen Cryptococcus neoformans. Surprisingly, we uncover a novel role of EVs in modulating antifungal drug resistance. Disruption of EV production was associated with altered lipid composition and changes in fluconazole susceptibility. Spontaneous azole-resistant mutants were deficient in EV production, while loss of resistance restored initial EV production levels. These findings were recapitulated in C. neoformans clinical isolates, indicating that azole resistance and EV production are coregulated in diverse strains. Our study reveals a new mechanism of drug resistance in which cells adapt to azole stress by modulating EV production.

Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Aneuploidy and gene dosage regulate filamentation and host colonization by Candida albicans.

Aneuploidy is a frequent occurrence in fungal species where it can alter gene expression and promote adaptation to a variety of environmental cues. Multiple forms of aneuploidy have been observed in the opportunistic fungal pathogen Candida albicans, which is a common component of the human gut mycobiome but can escape this niche and cause life-threatening systemic disease. Using a barcode sequencing (Bar-seq) approach, we evaluated a set of diploid C. albicans strains and found that a strain carrying a third copy of chromosome (Chr) 7 was associated with increased fitness during both gastrointestinal (GI) colonization and systemic infection. Our analysis revealed that the presence of a Chr 7 trisomy resulted in decreased filamentation, both in vitro and during GI colonization, relative to isogenic euploid controls. A target gene approach demonstrated that NRG1, encoding a negative regulator of filamentation located on Chr 7, contributes to increased fitness of the aneuploid strain due to inhibition of filamentation in a gene dosage-dependent fashion. Together, these experiments establish how aneuploidy enables the reversible adaptation of C. albicans to its host via gene dosage-dependent regulation of morphology.

High-throughput functional profiling of the human fungal pathogen Candida albicans genome.

Candida albicans is a major fungal pathogen of humans. Although its genome has been sequenced more than two decades ago, there are still over 4300 uncharacterized C. albicans genes. We previously generated an ORFeome as well as a collection of destination vectors to facilitate overexpression of C. albicans ORFs. Here, we report the construction of ∼2500 overexpression mutants and their evaluation by in vitro spotting on rich medium and in a liquid pool experiment in rich medium, allowing the identification of genes whose overexpression has a fitness cost. The candidates were further validated at the individual strain level. This new resource allows large-scale screens in different growth conditions to be performed routinely. Altogether, based on the concept of identifying functionally related genes by cluster analysis, the availability of this overexpression mutant collection will facilitate the characterization of gene functions in C. albicans.

A phylogenetically-restricted essential cell cycle progression factor in the human pathogen Candida albicans.

Chromosomal instability caused by cell division errors is associated with antifungal drug resistance in fungal pathogens. Here, we identify potential mechanisms underlying such instability by conducting an overexpression screen monitoring chromosomal stability in the human fungal pathogen Candida albicans. Analysis of ~1000 genes uncovers six chromosomal stability (CSA) genes, five of which are related to cell division genes of other organisms. The sixth gene, CSA6, appears to be present only in species belonging to the CUG-Ser clade, which includes C. albicans and other human fungal pathogens. The protein encoded by CSA6 localizes to the spindle pole bodies, is required for exit from mitosis, and induces a checkpoint-dependent metaphase arrest upon overexpression. Thus, Csa6 is an essential cell cycle progression factor that is restricted to the CUG-Ser fungal clade, and could therefore be explored as a potential antifungal target.

Candidalysins Are a New Family of Cytolytic Fungal Peptide Toxins.

Candidalysin is the first cytolytic peptide toxin identified in any human fungal pathogen. Candidalysin is secreted by Candida albicans and is critical for driving infection and host immune responses in several model systems. However, Candida infections are also caused by non-C. albicans species. Here, we identify and characterize orthologs of C. albicans candidalysin in C. dubliniensis and C. tropicalis. The candidalysins have different amino acid sequences, are amphipathic, and adopt a predominantly α-helical secondary structure in solution. Comparative functional analysis demonstrates that each candidalysin causes epithelial damage and calcium influx and activates intracellular signaling pathways and cytokine secretion. Importantly, C. dubliniensis and C. tropicalis candidalysins have higher damaging and activation potential than C. albicans candidalysin and exhibit more rapid membrane binding and disruption, although both fungal species cause less damage to epithelial cells than C. albicans. This study identifies the first family of peptide cytolysins in human-pathogenic fungi. IMPORTANCE Pathogenic fungi kill an estimated 1.5 million people every year. Recently, we discovered that the fungal pathogen Candida albicans secretes a peptide toxin called candidalysin during mucosal infection. Candidalysin causes damage to host cells, a process that supports disease progression. However, fungal infections are also caused by Candida species other than C. albicans. In this work, we identify and characterize two additional candidalysin toxins present in the related fungal pathogens C. dubliniensis and C. tropicalis. While the three candidalysins have different amino acid sequences, all three toxins are α-helical and amphipathic. Notably, the candidalysins from C. dubliniensis and C. tropicalis are more potent at inducing cell damage, calcium influx, mitogen-activated protein kinase signaling, and cytokine responses than C. albicans candidalysin, with the C. dubliniensis candidalysin having the most rapid membrane binding kinetics. These observations identify the candidalysins as the first family of peptide toxins in human-pathogenic fungi.

Population genomic analysis of Cryptococcus Brazilian isolates reveals an African type subclade distribution.

The genomes of a large number of Cryptococcus neoformans isolates have been sequenced and analyzed in recent years. These genomes have been used to understand the global population structure of this opportunistic pathogen. However, only a small number of South American isolates have been considered in these studies, and the population structure of C. neoformans in this part of the world remains elusive. Here, we analyzed the genomic sequences of 53 Brazilian Cryptococcus isolates and deciphered the C. neoformans population structure in this country. Our data reveal an African-like structure that suggested repeated intercontinental transports from Africa to South America. We also identified a mutator phenotype in one VNBII Brazilian isolate, exemplifying how fast-evolving isolates can shape the Cryptococcus population structure. Finally, phenotypic analyses revealed wide diversity but not lineage specificity in the expression of classical virulence traits within the set of isolates.

Factors that influence bidirectional long-tract homozygosis due to double-strand break repair in Candida albicans.

Genomic rearrangements have been associated with the acquisition of adaptive phenotypes, allowing organisms to efficiently generate new favorable genetic combinations. The diploid genome of Candida albicans is highly plastic, displaying numerous genomic rearrangements that are often the by-product of the repair of DNA breaks. For example, DNA double-strand breaks (DSB) repair using homologous-recombination pathways are a major source of loss-of-heterozygosity (LOH), observed ubiquitously in both clinical and laboratory strains of C. albicans. Mechanisms such as break-induced replication (BIR) or mitotic crossover (MCO) can result in long tracts of LOH, spanning hundreds of kilobases until the telomere. Analysis of I-SceI-induced BIR/MCO tracts in C. albicans revealed that the homozygosis tracts can ascend several kilobases toward the centromere, displaying homozygosis from the break site toward the centromere. We sought to investigate the molecular mechanisms that could contribute to this phenotype by characterizing a series of C. albicans DNA repair mutants, including pol32-/-, msh2-/-, mph1-/-, and mus81-/-. The impact of deleting these genes on genome stability revealed functional differences between Saccharomyces cerevisiae (a model DNA repair organism) and C. albicans. In addition, we demonstrated that ascending LOH tracts toward the centromere are associated with intrinsic features of BIR and potentially involve the mismatch repair pathway which acts upon natural heterozygous positions. Overall, this mechanistic approach to study LOH deepens our limited characterization of DNA repair pathways in C. albicans and brings forth the notion that centromere proximal alleles from DNA break sites are not guarded from undergoing LOH.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss.

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5-10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Use of CRISPR-Cas9 To Target Homologous Recombination Limits Transformation-Induced Genomic Changes in Candida albicans.

Most of our knowledge relating to molecular mechanisms of human fungal pathogenesis in Candida albicans relies on reverse genetics approaches, requiring strain engineering. DNA-mediated transformation of C. albicans has been described as highly mutagenic, potentially accentuated by the organism's genome plasticity, including the acquisition of genomic rearrangements, notably upon exposure to stress. The advent of CRISPR-Cas9 has vastly accelerated the process of genetically modifying strains, especially in diploid (such as C. albicans) and polyploid organisms. The effects of unleashing this nuclease within the genome of C. albicans are unknown, although several studies in other organisms report Cas9-associated toxicity and off-target DNA breaks. Upon the construction of a C. albicans strain collection, we took the opportunity to compare strains which were constructed using CRISPR-Cas9-free and CRISPR-Cas9-dependent transformation strategies, by quantifying and describing transformation-induced loss-of-heterozygosity and hyperploidy events. Our analysis of 57 strains highlights the mutagenic effects of transformation in C. albicans, regardless of the transformation protocol, but also underscores interesting differences in terms of genomic changes between strains obtained using different transformation protocols. Indeed, although strains constructed using the CRISPR-Cas9-free transformation method display numerous concomitant genomic changes randomly distributed throughout their genomes, the use of CRISPR-Cas9 leads to a reduced overall number of genome changes, particularly hyperploidies. Overall, in addition to facilitating strain construction by reducing the number of transformation steps, the CRISPR-Cas9-dependent transformation strategy in C. albicans appears to limit transformation-associated genome changes.IMPORTANCE Genome editing is essential to nearly all research studies aimed at gaining insight into the molecular mechanisms underlying various biological processes, including those in the opportunistic pathogen Candida albicans The adaptation of the CRISPR-Cas9 system greatly facilitates genome engineering in many organisms. However, our understanding of the effects of CRISPR-Cas9 technology on the biology of C. albicans is limited. In this study, we sought to compare the extents of transformation-induced genomic changes within strains engineered using CRISPR-Cas9-free and CRISPR-Cas9-dependent transformation methods. CRISPR-Cas9-dependent transformation allows one to simultaneously target both homologs and, importantly, appears less mutagenic in C. albicans, since strains engineered using CRISPR-Cas9 display an overall decrease in concomitant genomic changes.