Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes.

We used whole-genome sequencing to analyze a collection of 35 fluconazole-resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two persistent clinical lineages were identified. We identified copy number variation (CNV) of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show that the annotated telomeric gene CDR1B is actually an artifactual in silico fusion of two highly similar neighboring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased the expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Genome Analysis of a Newly Discovered Yeast Species, Hanseniaspora menglaensis.

Annual surveys of Irish soil samples identified three isolates, CBS 16921 (UCD88), CBS 18246 (UCD443), and CBS 18247 (UCD483), of an apiculate yeast species within the Hanseniaspora genus. The internal transcribed spacer (ITS) and D1/D2 region of the large subunit (LSU) rRNA sequences showed that these are isolates of the recently described species Hanseniaspora menglaensis, first isolated from Southwest China. No genome sequence for H. menglaensis is currently available. The genome sequences of the three Irish isolates were determined using short-read (Illumina) sequencing, and the sequence of one isolate (CBS 16921) was assembled to chromosome level using long-read sequencing (Oxford Nanopore Technologies). Phylogenomic analysis shows that H. menglaensis belongs to the fast-evolving lineage (FEL) of Hanseniaspora. Only one MAT idiomorph (encoding MATα1) was identified in all three sequenced H. menglaensis isolates, consistent with one mating type of a heterothallic species. Genome comparisons showed that there has been a rearrangement near MATα of FEL species compared to isolates from the slowly evolving lineage (SEL).

Evaluation of outbreak persistence caused by multidrug-resistant and echinocandin-resistant Candida parapsilosis using multidimensional experimental and epidemiological approaches.

Candida parapsilosis is known to cause severe and persistent outbreaks in clinical settings. Patients infected with multidrug-resistant C. parapsilosis (MDR Cp) isolates were identified in a large Turkish hospital from 2017-2020. We subsequently identified three additional patients infected with MDR Cp isolates in 2022 from the same hospital and two echinocandin-resistant (ECR) isolates from a single patient in another hospital. The increasing number of MDR and ECR isolates contradicts the general principle that the severe fitness cost associated with these phenotypes could prevent their dominance in clinical settings. Here, we employed a multidimensional approach to systematically assess the fitness costs of MDR and ECR C. parapsilosis isolates. Whole-genome sequencing revealed a novel MDR genotype infecting two patients in 2022. Despite severe in vitro defects, the levels and tolerances of the biofilms of our ECR and MDR isolates were generally comparable to those of susceptible wild-type isolates. Surprisingly, the MDR and ECR isolates showed major alterations in their cell wall components, and some of the MDR isolates consistently displayed increased tolerance to the fungicidal activities of primary human neutrophils and were more immunoevasive during exposure to primary human macrophages. Our systemic infection mouse model showed that MDR and ECR C. parapsilosis isolates had comparable fungal burden in most organs relative to susceptible isolates. Overall, we observed a notable increase in the genotypic diversity and frequency of MDR isolates and identified MDR and ECR isolates potentially capable of causing persistent outbreaks in the future.

Draft genome sequence of the yeast Schwanniomyces capriottii strain UCD805, isolated from soil in Ireland.

Schwanniomyces capriottii is a member of the Debaryomycetaceae family in the order Saccharomycetales. Here, we present the genome sequence of S. capriottii UCD805, which was isolated from soil in Dublin, Ireland. This genome is 12.2 Mb and was assembled into 14 scaffolds plus a mitochondrial genome scaffold.

Genome sequences of two isolates of the yeast Candida zeylanoides: UCD849 from soil in Ireland, and AWD from an African wild dog.

We report genome sequences of two new isolates of the budding yeast Candida zeylanoides. Strain UCD849 from soil in Ireland was assembled into eight complete chromosomes. Strain AWD from an African Wild Dog in a US zoo was sequenced to draft level. The assemblies are 10.6 Mb and 99.57% identical.

Whole-genome sequencing confirms a persistent candidaemia clonal outbreak due to multidrug-resistant Candida parapsilosis.

OBJECTIVES: Although perceived as a rare clinical entity, recent studies have noted the emergence of MDR C. parapsilosis (MDR-Cp) isolates from single patients (resistant to both azole and echinocandins). We previously reported a case series of MDR-Cp isolates carrying a novel FKS1R658G mutation. Herein, we identified an echinocandin-naive patient infected with MDR-Cp a few months after the previously described isolates. WGS and CRISPR-Cas9 editing were used to explore the origin of the new MDR-Cp isolates, and to determine if the novel mutation confers echinocandin resistance. METHODS: WGS was applied to assess the clonality of these isolates and CRISPR-Cas9 editing and a Galleria mellonella model were used to examine whether FKS1R658G confers echinocandin resistance. RESULTS: Fluconazole treatment failed, and the patient was successfully treated with liposomal amphotericin B (LAMB). WGS proved that all historical and novel MDR-Cp strains were clonal and distant from the fluconazole-resistant outbreak cluster in the same hospital. CRISPR-Cas9 editing and G. mellonella virulence assays confirmed that FKS1R658G confers echinocandin resistance in vitro and in vivo. Interestingly, the FKS1R658G mutant showed a very modest fitness cost compared with the parental WT strain, consistent with the persistence of the MDR-Cp cluster in our hospital. CONCLUSIONS: Our study showcases the emergence of MDR-Cp isolates as a novel threat in clinical settings, which undermines the efficacy of the two most widely used antifungal drugs against candidiasis, leaving only LAMB as a last resort. Additionally, surveillance studies and WGS are warranted to effectively establish infection control and antifungal stewardship strategies.

Candida parapsilosis isolates carrying mutations outside FKS1 hotspot regions confer high echinocandin tolerance and facilitate the development of echinocandin resistance.

Candida parapsilosis is a significant cause of candidemia worldwide. Echinocandin-resistant (ECR) and echinocandin-tolerant (ECT) C. parapsilosis isolates have been reported in various countries but are rare. Resistance and tolerance are predominantly caused by mutations related to the hotspot (HS) regions of the FKS1 gene. A relatively high proportion of clinical C. parapsilosis isolates carrying mutations outside the HS regions has been noted in some studies, but an association with echinocandin (EC) resistance or tolerance was not explored. Herein, CRISPR-Cas9 was used and the association between amino acid substitution in FKS1 outside HS 1/2 (V595I, S745L, M1328I, F1386S, and A1422G) with EC susceptibility profile was delineated. None of the mutations conferred EC resistance, but they resulted in a significantly higher level of EC tolerance than the parental isolate, ATCC 22019. When incubated on agar plates containing ECs, specifically caspofungin and micafungin, ECR colonies were exclusively observed among ECT isolates, particularly mutants carrying V595I, S745L, and F1386S. Additionally, mutants had significantly better growth rates in yeast extract peptone dextrose (YPD) and YPD containing agents inducing membrane and oxidative stresses. The mutants had a trivial fitness cost in the Galleria mellonella model relative to ATCC 22019. Collectively, this study supports epidemiological studies to catalog mutations occurring outside the HS regions of FKS1, even if they do not confer EC resistance. These mutations are important as they potentially confer a higher level of EC tolerance and a higher propensity to develop EC resistance, therefore unveiling a novel mechanism of EC tolerance in C. parapsilosis. The identification of EC tolerance in C. parapsilosis may have direct clinical benefit in patient management.

Acquired fluconazole resistance and genetic clustering in Diutina (Candida) catenulata from clinical samples.

OBJECTIVES: Diutina (Candida) catenulata is an ascomycetous yeast isolated from environmental sources and animals, occasionally infecting humans. The aim of this study is to shed light on the in vitro antifungal susceptibility and genetic diversity of this opportunistic yeast. METHODS: Forty-five D. catenulata strains isolated from various sources (including human and environmental sources) and originating from nine countries were included. Species identification was performed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and confirmed via internal transcribed spacer ribosomal DNA barcoding. In vitro antifungal susceptibility was determined for seven systemic antifungals via the gradient strip method after 48 hours of incubation at 35°C using Etest® (Biomérieux) or Liofilchem® strips. Isolates exhibiting fluconazole minimal inhibitory concentrations (MICs) of ≥8 µg/mL were investigated for mutations in the ERG11 gene. A novel microsatellite genotyping scheme consisting of four markers was developed to assess genetic diversity. RESULTS: MIC ranges for amphotericin B, caspofungin, micafungin, isavuconazole, and posaconazole were 0.19-1 µg/mL, 0.094-0.5 µg/mL, 0.012-0.064 µg/mL, 0.003-0.047 µg/mL, and 0.006-0.032 µg/mL, respectively. By comparison, a broad range of MICs was noted for fluconazole (0.75 to >256 µg/mL) and voriconazole (0.012-0.38 mg/L), the higher values being observed among clinical strains. The Y132F amino acid substitution, associated with azole resistance in various Candida species (C. albicans, C. tropicalis, C. parapsilosis, and C. orthopsilosis), was the main substitution identified. Although microsatellite typing showed extensive genetic diversity, most strains with high fluconazole MICs clustered together, suggesting human-to-human transmission or a common source of contamination. DISCUSSION: The high rate of acquired fluconazole resistance among clinical isolates of D. catenulata is of concern. In this study, we highlight a link between the genetic diversity of D. catenulata and its antifungal resistance patterns, suggesting possible clonal transmission of resistant isolates.

Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes.

We used whole-genome sequencing to analyse a collection of 35 fluconazole resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two probable outbreak groups were identified. We identified copy number variation of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show the annotated telomeric gene CDR1B is actually an artefactual in silico fusion of two highly similar neighbouring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Identification of European isolates of the lager yeast parent Saccharomyces eubayanus.

Lager brewing first occurred in Bavaria in the 15th century, associated with restrictions of brewing to colder months. The lager yeast, Saccharomyces pastorianus, is cold tolerant. It is a hybrid between Saccharomyces cerevisiae and Saccharomyces eubayanus, and has been found only in industrial settings. Natural isolates of S. eubayanus were first discovered in Patagonia 11 years ago. They have since been isolated from China, Tibet, New Zealand, and North America, but not from Europe. Here, we describe the first European strains UCD646 and UCD650, isolated from a wooded area on a university campus in Dublin, Ireland. We generated complete chromosome level assemblies of both genomes using long- and short-read sequencing. The UCD isolates belong to the Holarctic clade. Genome analysis shows that isolates similar to the Irish strains contributed to the S. eubayanus component of S. pastorianus, but isolates from Tibet made a larger contribution.

CRISPR-Cas9 Editing Induces Loss of Heterozygosity in the Pathogenic Yeast Candida parapsilosis.

Genetic manipulation is often used to study gene function. However, unplanned genome changes (including single nucleotide polymorphisms [SNPs], aneuploidy, and loss of heterozygosity [LOH]) can affect the phenotypic traits of the engineered strains. Here, we compared the effect of classical deletion methods (replacing target alleles with selectable markers by homologous recombination) with CRISPR-Cas9 editing in the diploid human-pathogenic yeast Candida parapsilosis. We sequenced the genomes of 9 isolates that were modified using classic recombination methods and 12 that were edited using CRISPR-Cas9. As a control, the genomes of eight isolates that were transformed with a Cas9-expressing plasmid in the absence of a guide RNA were also sequenced. Following gene manipulation using classic homologous recombination, only one strain exhibited extensive LOH near the targeted gene (8.9 kb), whereas another contained multiple LOH events not associated with the intended modification. In contrast, large regions of LOH (up to >1,100 kb) were observed in most CRISPR-Cas9-edited strains. LOH most commonly occurred adjacent to the Cas9 cut site and extended to the telomere in four isolates. In two isolates, we observed LOH on chromosomes that were not targeted by CRISPR-Cas9. Among the CRISPR-edited isolates, two exhibited cysteine and methionine auxotrophy caused by LOH at a heterozygous site in MET10, approximately 11 and 157 kb downstream from the Cas9 target site, respectively. C. parapsilosis isolates have relatively low levels of heterozygosity. However, our results show that mutation complementation to confirm observed phenotypes is required when using CRISPR-Cas9. IMPORTANCE CRISPR-Cas9 has greatly streamlined gene editing and is now the gold standard and first choice for genetic engineering. However, we show that in diploid species, extra care should be taken in confirming the cause of any phenotypic changes observed. We show that the Cas9-induced double-strand break is often associated with loss of heterozygosity in the asexual diploid human fungal pathogen Candida parapsilosis. This can result in deleterious heterozygous variants (e.g., stop gain in one allele) becoming homozygous, resulting in unplanned phenotypic changes. Our results stress the importance of mutation complementation even when using CRISPR-Cas9.

Draft Genome Sequence of the Yeast Blastobotrys aristata Strain UCD613, Isolated from Soil in Ireland.

Blastobotrys aristata is a member of the Trichomonascaceae family in the order Saccharomycetales. Here, we present the genome sequence of B. aristata UCD613, which was isolated from soil in Dublin, Ireland. This genome is 13.3 Mb and was assembled into 4 chromosome-size scaffolds of >2.2 Mb in size plus a mitochondrial genome scaffold.

Draft Genome Sequence of the Yeast Torulaspora quercuum Strain UCD657, Isolated from Soil in Ireland.

Torulaspora quercuum is an ascomycete yeast first isolated in 2009. Here, we present the genome sequence of T. quercuum isolate UCD657, which was isolated from soil in Ireland. This genome is 10.4 Mb and was assembled into 8 chromosome-sized scaffolds of >1 Mb in size, plus a mitochondrial genome scaffold.

Determinants of fluconazole resistance and echinocandin tolerance in C. parapsilosis isolates causing a large clonal candidemia outbreak among COVID-19 patients in a Brazilian ICU.

Patients presenting with severe COVID-19 are predisposed to acquire secondary fungal infections such as COVID-19-associated candidemia (CAC), which are associated with poor clinical outcomes despite antifungal treatment. The extreme burden imposed on clinical facilities during the COVID-19 pandemic has provided a permissive environment for the emergence of clonal outbreaks of multiple Candida species, including C. auris and C. parapsilosis. Here we report the largest clonal CAC outbreak to date caused by fluconazole resistant (FLZR) and echinocandin tolerant (ECT) C. parapsilosis. Sixty C. parapsilosis strains were obtained from 57 patients at a tertiary care hospital in Brazil, 90% of them were FLZR and ECT. Although only 35.8% of FLZR isolates contained an ERG11 mutation, all of them contained the TAC1L518F mutation and significantly overexpressed CDR1. Introduction of TAC1L518F into a susceptible background increased the MIC of fluconazole and voriconazole 8-fold and resulted in significant basal overexpression of CDR1. Additionally, FLZR isolates exclusively harboured E1939G outside of Fks1 hotspot-2, which did not confer echinocandin resistance, but significantly increased ECT. Multilocus microsatellite typing showed that 51/60 (85%) of the FLZR isolates belonged to the same cluster, while the susceptible isolates each represented a distinct lineage. Finally, biofilm production in FLZR isolates was significantly lower than in susceptible counterparts Suggesting that it may not be an outbreak determinant. In summary, we show that TAC1L518F and FKS1E1393G confer FLZR and ECT, respectively, in CAC-associated C. parapsilosis. Our study underscores the importance of antifungal stewardship and effective infection control strategies to mitigate clonal C. parapsilosis outbreaks.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance.

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Plasmid-Based CRISPR-Cas9 Editing in Multiple Candida Species.

CRISPR-Cas9 technology radically changed the approach to genetic manipulation of both medically and industrially relevant Candida species, as attested by the ever-increasing number of applications to the study of pathogenesis, drug resistance, gene expression, and host pathogen interaction and drug discovery. Here, we describe the use of plasmid-based systems for high efficiency CRISPR-Cas9 gene editing into any strain of four non-albicans Candida species, namely, Candida parapsilosis, Candida orthopsilosis, Candida metapsilosis, and Candida tropicalis. The plasmids pCP-tRNA and pCT-tRNA contain all the elements necessary for expressing the CRISPR-Cas9 machinery, and they can be used in combination with a repair template for disrupting gene function by insertion of a premature stop codon or by gene deletion. The plasmids are easily lost in the absence of selection, allowing scarless gene editing and minimizing detrimental effects of prolonged Cas9 expression.

Clinical azole cross-resistance in Candida parapsilosis is related to a novel MRR1 gain-of-function mutation.

OBJECTIVES: Hereby, we describe the molecular mechanisms underlying the acquisition of azole resistance by a Candida parapsilosis isolate following fluconazole treatment due to candiduria. METHODS: A set of three consecutive C. parapsilosis isolates were recovered from the urine samples of a patient with candiduria. Whole-genome sequencing and antifungal susceptibility assays were performed. The expression of MRR1, MDR1, ERG11 and CDR1B (CPAR2_304370) was quantified by RT-qPCR. RESULTS: The initial isolate CPS-A was susceptible to all three azoles tested (fluconazole, voriconazole and posaconazole); isolate CPS-B, collected after the second cycle of treatment, exhibited a susceptible-dose-dependent phenotype to fluconazole and isolate CPS-C, recovered after the third cycle, exhibited a cross-resistance profile to fluconazole and voriconazole. Whole-genome sequencing revealed a putative resistance mechanism in isolate CPS-C, associated with a G1810A nucleotide substitution, leading to a G604R change in the Mrr1p transcription factor. Introducing this mutation into the susceptible CPS-A isolate (MRR1RI) resulted in resistance to fluconazole and voriconazole, as well as up-regulation of MRR1 and MDR1. Interestingly, the susceptible-dose-dependent phenotype exhibited by isolate CPS-B was associated with an increased copy number of the CDR1B gene. The expression of CDR1B was increased in both isolates CPS-B and CPS-C and in the MRR1RI strain, harbouring the gain-of-function mutation. CONCLUSIONS: Our results describe clinical azole cross-resistance acquisition in C. parapsilosis due to a G1810A (G604R) gain-of-function mutation, resulting in MRR1 hyperactivation and consequently, MDR1 efflux pump overexpression. We also associated amplification of the CDR1B gene with decreased fluconazole susceptibility and showed that it is a putative target of the MRR1 gain-of-function mutation.

Mating-Type Switching in Budding Yeasts, from Flip/Flop Inversion to Cassette Mechanisms.

Mating-type switching is a natural but unusual genetic control process that regulates cell identity in ascomycete yeasts. It involves physically replacing one small piece of genomic DNA by another, resulting in replacement of the master regulatory genes in the mating pathway and hence a switch of cell type and mating behavior. In this review, we concentrate on recent progress that has been made on understanding the origins and evolution of mating-type switching systems in budding yeasts (subphylum Saccharomycotina). Because of the unusual nature and the complexity of the mechanism in Saccharomyces cerevisiae, mating-type switching was assumed until recently to have originated only once or twice during yeast evolution. However, comparative genomics analysis now shows that switching mechanisms arose many times independently-at least 11 times in budding yeasts and once in fission yeasts-a dramatic example of convergent evolution. Most of these lineages switch mating types by a flip/flop mechanism that inverts a section of a chromosome and is simpler than the well-characterized 3-locus cassette mechanism (MAT/HML/HMR) used by S. cerevisiae. Mating-type switching (secondary homothallism) is one of the two possible mechanisms by which a yeast species can become self-fertile. The other mechanism (primary homothallism) has also emerged independently in multiple evolutionary lineages of budding yeasts, indicating that homothallism has been favored strongly by natural selection. Recent work shows that HO endonuclease, which makes the double-strand DNA break that initiates switching at the S. cerevisiae MAT locus, evolved from an unusual mobile genetic element that originally targeted a glycolytic gene, FBA1.

Characterization of the Candida glabrata Transcription Factor CgMar1: Role in Azole Susceptibility.

The prevalence of antifungal resistance in Candida glabrata, especially against azole drugs, results in difficult-to-treat and potentially life-threatening infections. Understanding the molecular basis of azole resistance in C. glabrata is crucial to designing more suitable therapeutic strategies. In this study, the role of the transcription factor encoded by ORF CAGL0B03421g, here denominated as CgMar1 (Multiple Azole Resistance 1), in azole susceptibility was explored. Using RNA-sequencing, CgMar1 was found to regulate 337 genes under fluconazole stress, including several related to lipid biosynthesis pathways. In this context, CgMar1 and its target CgRSB1, encoding a predicted sphingoid long-chain base efflux transporter, were found to contribute to plasma membrane sphingolipid incorporation and membrane permeability, decreasing fluconazole accumulation. CgMar1 was found to associate with the promoter of CgRSB1, which contains two instances of the CCCCTCC consensus, found to be required for CgRSB1 activation during fluconazole stress. Altogether, a regulatory pathway modulating azole susceptibility in C. glabrata is proposed, resulting from what appears to be a neofunctionalization of a Hap1-like transcription factor.