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.

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.

Quantitative global studies reveal differential translational control by start codon context across the fungal kingdom.

Eukaryotic protein synthesis generally initiates at a start codon defined by an AUG and its surrounding Kozak sequence context, but the quantitative importance of this context in different species is unclear. We tested this concept in two pathogenic Cryptococcus yeast species by genome-wide mapping of translation and of mRNA 5' and 3' ends. We observed thousands of AUG-initiated upstream open reading frames (uORFs) that are a major contributor to translation repression. uORF use depends on the Kozak sequence context of its start codon, and uORFs with strong contexts promote nonsense-mediated mRNA decay. Transcript leaders in Cryptococcus and other fungi are substantially longer and more AUG-dense than in Saccharomyces. Numerous Cryptococcus mRNAs encode predicted dual-localized proteins, including many aminoacyl-tRNA synthetases, in which a leaky AUG start codon is followed by a strong Kozak context in-frame AUG, separated by mitochondrial-targeting sequence. Analysis of other fungal species shows that such dual-localization is also predicted to be common in the ascomycete mould, Neurospora crassa. Kozak-controlled regulation is correlated with insertions in translational initiation factors in fidelity-determining regions that contact the initiator tRNA. Thus, start codon context is a signal that quantitatively programs both the expression and the structures of proteins in diverse fungi.

Large-scale genome mining allows identification of neutral polymorphisms and novel resistance mutations in genes involved in Candida albicans resistance to azoles and echinocandins.

BACKGROUND: The genome of Candida albicans displays significant polymorphism. Point mutations in genes involved in resistance to antifungals may either confer phenotypic resistance or be devoid of phenotypic consequences. OBJECTIVES: To catalogue polymorphisms in azole and echinocandin resistance genes occurring in susceptible strains in order to rapidly pinpoint relevant mutations in resistant strains. METHODS: Genome sequences from 151 unrelated C. albicans strains susceptible to fluconazole and caspofungin were used to create a catalogue of non-synonymous polymorphisms in genes involved in resistance to azoles (ERG11, TAC1, MRR1 and UPC2) or echinocandins (FKS1). The potential of this catalogue to reveal putative resistance mutations was tested in 10 azole-resistant isolates, including 1 intermediate to caspofungin. Selected mutations were analysed by mutagenesis experiments or mutational prediction effect. RESULTS: In the susceptible strains, we identified 126 amino acid substitutions constituting the catalogue of phenotypically neutral polymorphisms. By excluding these neutral substitutions, we identified 22 additional substitutions in the 10 resistant strains. Among these substitutions, 10 had already been associated with resistance. The remaining 12 were in Tac1p (n = 6), Upc2p (n = 2) and Erg11p (n = 4). Four out of the six homozygous substitutions in Tac1p (H263Y, A790V, H839Y and P971S) conferred increases in azole MICs, while no effects were observed for those in Upc2p. Additionally, two homozygous substitutions (Y64H and P236S) had a predicted conformation effect on Erg11p. CONCLUSIONS: By establishing a catalogue of neutral polymorphisms occurring in genes involved in resistance to antifungal drugs, we provide a useful resource for rapid identification of mutations possibly responsible for phenotypic resistance in C. albicans.

Epidemiological investigation for grouped cases of Trichosporon asahii using whole genome and IGS1 sequencing.

BACKGROUND: Trichosporonosis is a rare invasive infection in humans mainly due to Trichosporon asahii, and especially recovered from patients having haematological malignancy. Since 2012, IGS1 region sequencing is used as a genotyping method to distinguish isolates, with high frequency of one haplotype worldwide and a geographic specificity for some haplotypes. OBJECTIVES: We compared the IGS1 genotyping method and whole genome sequencing (WGS) to study the relationship between clinical isolates involved in two grouped cases in France. METHODS: IGS1 sequencing and antifungal susceptibility testing were performed for 54 clinical isolates. Clinical data for 28 isolates included in surveillance programs were analysed. Whole genome was sequenced for 32 clinical isolates and the type strain. RESULTS: All isolates were intrinsically resistant to flucytosine, while voriconazole had the most potent in vitro activity. The majority of the isolates was recovered from patients with haematological malignancies (42.86%), with a high proportion of children (<15 yrs-old, 32.14%) and a high mortality rate at three months (46.15%). Based on the WGS analysis, isolates exhibiting IGS1 haplotype 1, 3 and 7 belonged to different clades. Five isolates recovered during the first grouped cases had the same IGS1 haplotype and shared 99% of SNPs similarity. For the second grouped cases, four isolates had 98.7% of SNPs similarity while the isolate recovered 4 years earlier was totally unlinked. CONCLUSIONS: We confirmed the usefulness of IGS1 sequencing for grouped cases infection of T. asahii. We underlined its limitation for the study of population structure and the utility of WGS analysis for the study of epidemiologically unrelated isolates.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.