Candida parapsilosis Mdr1B and Cdr1B Are Drivers of Mrr1-Mediated Clinical Fluconazole Resistance

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.

Adaptation of Candida albicans to specific host environments by gain-of-function mutations in transcription factors

The yeast Candida albicans is usually a harmless member of the normal microbiota in healthy persons but is also a major fungal pathogen that can colonize and infect almost every human tissue. A successful adaptation to environmental changes encountered in different host niches requires an appropriate regulation of gene expression. The zinc cluster transcription factors are the largest family of transcriptional regulators in C. albicans and are involved in the control of virtually all aspects of its biology. Under certain circumstances, mutations in these transcription factors that alter their activity and the expression of their target genes confer a selective advantage, which results in the emergence of phenotypically altered variants that are better adapted to new environmental challenges. This review describes how gain-of-function mutations in different zinc cluster transcription factors enable C. albicans to overcome antifungal therapy and to successfully establish itself in specific host niches.

High-Resolution Melting Assay to Detect the Mutations That Cause the Y132F and G458S Substitutions at the ERG11 Gene Involved in Azole Resistance in Candida parapsilosis

BACKGROUND

Candida parapsilosis is a pathogenic yeast that has reduced susceptibility to echinocandins and ranks as the second or third leading cause of candidaemia, depending on the geographical region. This yeast often causes nosocomial infections, which are frequently detected as outbreaks. In recent years, resistance to azoles in C. parapsilosis has increased globally, primarily due to the accumulation of mutations in the ERG11 gene.

OBJECTIVES

In this study, we have developed an assay based on real-time PCR and high-resolution melting (HRM) curve analysis to detect two of the most prevalent mutations at ERG11 that confer resistance to fluconazole (Y132F and G458S).

METHODS

We designed allele-specific oligonucleotides that selectively bind to either the wild type or mutated sequences and optimised the conditions to ensure amplification of the specific allele, followed by detection via high-resolution melting (HRM) analysis.

RESULTS

The designed oligonucleotides to detect the Erg11Y132F and Erg11G458S mutations produced specific amplification of either WT or mutated alleles. We conducted a duplex real-time PCR combining oligonucleotides for the wild-type sequences in one mix, and oligonucleotides for the mutated alleles in another. Following this, we performed an analysis of the HRM curve to identify the amplified allele in each case. This technique was blindly evaluated on a set of 114 C. parapsilosis isolates, all of which were unequivocally identified using our approach.

CONCLUSION

This technique offers a new method for the early detection of azole resistance mechanism in C. parapsilosis.

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.

Precise genome editing underlines the distinct contributions of mutations in ERG11, ERG3, MRR1, and TAC1 genes to antifungal resistance in Candida parapsilosis

Candida parapsilosis has recently emerged as a major threat due to the worldwide emergence of fluconazole-resistant strains causing clonal outbreaks in hospitals and poses a therapeutic challenge due to the limited antifungal armamentarium. Here, we used precise genome editing using CRISPR-Cas9 to gain further insights into the contribution of mutations in ERG11, ERG3, MRR1, and TAC1 genes and the influence of allelic dosage to antifungal resistance in C. parapsilosis. Seven of the most common amino acid substitutions previously reported in fluconazole-resistant clinical isolates (including Y132F in ERG11) were engineered in two fluconazole-susceptible C. parapsilosis lineages (ATCC 22019 and STZ5). Each mutant was then challenged in vitro against a large array of antifungals, with a focus on azoles. Any possible change in virulence was also assessed in a Galleria mellonella model. We successfully generated a total of 19 different mutants, using CRISPR-Cas9. Except for R398I (ERG11), all remaining amino acid substitutions conferred reduced susceptibility to fluconazole. However, the impact on fluconazole in vitro susceptibility varied greatly according to the engineered mutation, the stronger impact being noted for G583R acting as a gain-of-function mutation in MRR1. Cross-resistance with newer azoles, non-medical azoles, but also non-azole antifungals such as flucytosine, was occasionally noted. Posaconazole and isavuconazole remained the most active in vitro. Except for G583R, no fitness cost was associated with the acquisition of fluconazole resistance. We highlight the distinct contributions of amino acid substitutions in ERG11, ERG3, MRR1, and TAC1 genes to antifungal resistance in C. parapsilosis.

First Canadian report of transmission of fluconazole-resistant Candida parapsilosis within two hospital networks confirmed by genomic analysis

Candida parapsilosis is a common cause of non-albicans candidemia. It can be transmitted in healthcare settings resulting in serious healthcare-associated infections and can develop drug resistance to commonly used antifungal agents. Following a significant increase in the percentage of fluconazole (FLU)-nonsusceptible isolates from sterile site specimens of patients in two Ontario acute care hospital networks, we used whole genome sequence (WGS) analysis to retrospectively investigate the genetic relatedness of isolates and to assess potential in-hospital spread. Phylogenomic analysis was conducted on all 19 FLU-resistant and seven susceptible-dose dependent (SDD) isolates from the two hospital networks, as well as 13 FLU susceptible C. parapsilosis isolates from the same facilities and 20 isolates from patients not related to the investigation. Twenty-five of 26 FLU-nonsusceptible isolates (resistant or SDD) and two susceptible isolates from the two hospital networks formed a phylogenomic cluster that was highly similar genetically and distinct from other isolates. The results suggest the presence of a persistent strain of FLU-nonsusceptible C. parapsilosis causing infections over a 5.5-year period. Results from WGS were largely comparable to microsatellite typing. Twenty-seven of 28 cluster isolates had a K143R substitution in lanosterol 14-α-demethylase (ERG11) associated with azole resistance. As the first report of a healthcare-associated outbreak of FLU-nonsusceptible C. parapsilosis in Canada, this study underscores the importance of monitoring local antimicrobial resistance trends and demonstrates the value of WGS analysis to detect and characterize clusters and outbreaks. Timely access to genomic epidemiological information can inform targeted infection control measures.

Emergence of multiple fluconazole-resistant Candida parapsilosis sensu stricto clones with persistence and transmission in China

OBJECTIVES

We explored the epidemiological and molecular characteristics of Candida parapsilosis sensu stricto isolates in China, and their mechanisms of azole resistance.

METHODS

Azole susceptibilities of 2318 non-duplicate isolates were determined using CLSI broth microdilution. Isolates were genotyped by a microsatellite typing method. Molecular resistance mechanisms were also studied and functionally validated by CRISPR/Cas9-based genetic alterations.

RESULTS

Fluconazole resistance occurred in 2.4% (n = 56) of isolates, and these isolates showed a higher frequency of distribution in ICU inpatients compared with susceptible isolates (48.2%, n = 27/56 versus 27.8%, 613/2208; P = 0.019). Microsatellite-genotyping analysis yielded 29 genotypes among 56 fluconazole-resistant isolates, of which 10 genotypes, including 37 isolates, belonged to clusters, persisting and transmitting in Chinese hospitals for 1-29 months. Clusters harbouring Erg11Y132F (5/10; 50%) were predominant in China. Among these, the second most dominant cluster MT07, including seven isolates, characteristically harbouring Erg11Y132F and Mrr1Q625K, lent its carriage to being one of the strongest associations with cross-resistance and high MICs of fluconazole (>256 mg/L) and voriconazole (2-8 mg/L), causing transmission across two hospitals. Among mutations tested, Mrr1Q625K led to the highest-level increase of fluconazole MIC (32-fold), while mutations located within or near the predicted transcription factor domain of Tac1 (D440Y, T492M and L518F) conferred cross-resistance to azoles.

CONCLUSIONS

This study is the first Chinese report of persistence and transmissions of multiple fluconazole-resistant C. parapsilosis sensu stricto clones harbouring Erg11Y132F, and the first demonstration of the mutations Erg11G307A, Mrr1Q625K, Tac1L263S, Tac1D440Y and Tac1T492M as conferring resistance to azoles.

Mutations in TAC1 and ERG11 are major drivers of triazole antifungal resistance in clinical isolates of Candida parapsilosis

OBJECTIVES

The aim of this study was to determine how mutations in CpERG11 and CpTAC1 contribute to fluconazole resistance in a collection of clinical isolates.

METHODS

Sequences of CpERG11 and CpTAC1 were determined for 35 resistant Candida parapsilosis clinical isolates. A plasmid-based CRISPR-Cas9 system was used to introduce mutations leading to amino acid substitution in CpTac1 and CpErg11. Triazole susceptibility was determined by broth microdilution and E-test. Differential expression of genes mediated by CpTAC1 mutation was determined by RNA sequencing, and relative expression of individual transporter genes was assessed with RT-qPCR.

RESULTS

Six isolates carried a mutation in CpTAC1 in combination with the CpERG11 mutation, leading to the CpErg11Y132F substitution. When introduced into susceptible isolates, this CpERG11 mutation led to a 4- to 8-fold increase in fluconazole minimum inhibitory concentrations (MIC; 0.125 μg/mL vs. 0.5 μg/mL, 0.125 μg/mL vs. 0.5 μg/mL, and 0.5 μg/mL vs. 4 μg/mL). When introduced into a susceptible isolate, the CpTAC1 mutation leading to the G650E substitution resulted in an 8-fold increase in fluconazole MIC (0.25 μg/mL vs. 2 μg/mL), whereas correction of this mutation in resistant isolates led to a 16-fold reduction in MIC (32 μg/mL vs. 2 μg/mL). CpCDR1, CpCDR1B, and CpCDR1C were overexpressed in the presence CpTac1G650E. Disruption of these genes in combination resulted in a 4-fold reduction in fluconazole MIC (32 μg/mL vs. 8 μg/mL).

DISCUSSION

These results define the specific contribution made by the Y132F substitution in CpERG11 and demonstrate a role for activating mutations in CpTAC1 in triazole resistance in C. parapsilosis.

Candida parapsilosis sensu stricto Antifungal Resistance Mechanisms and Associated Epidemiology

Fungal diseases cause millions of deaths per year worldwide. Antifungal resistance has become a matter of great concern in public health. In recent years rates of non-albicans species have risen dramatically. Candida parapsilosis is now reported to be the second most frequent species causing candidemia in several countries in Europe, Latin America, South Africa and Asia. Rates of acquired azole resistance are reaching a worrisome threshold from multiple reports as in vitro susceptibility testing is now starting also to explore tolerance and heteroresistance to antifungal compounds. With this review, the authors seek to evaluate known antifungal resistance mechanisms and their worldwide distribution in Candida species infections with a specific focus on C. parapsilosis.

Candida parapsilosis Virulence and Antifungal Resistance Mechanisms: A Comprehensive Review of Key Determinants

Candida parapsilosis is the second most common Candida species isolated in Asia, Southern Europe, and Latin America and is often involved in invasive infections that seriously impact human health. This pathogen is part of the psilosis complex, which also includes Candida orthopsilosis and Candida metapsilosis. C. parapsilosis infections are particularly prevalent among neonates with low birth weights, individuals who are immunocompromised, and patients who require prolonged use of a central venous catheter or other indwelling devices, whose surfaces C. parapsilosis exhibits an enhanced capacity to adhere to and form biofilms. Despite this well-acknowledged prevalence, the biology of C. parapsilosis has not been as extensively explored as that of Candida albicans. In this paper, we describe the molecular mechanistic pathways of virulence in C. parapsilosis and show how they differ from those of C. albicans. We also describe the mode of action of antifungal drugs used for the treatment of Candida infections, namely, polyenes, echinocandins, and azoles, as well as the resistance mechanisms developed by C. parapsilosis to overcome them. Finally, we stress the importance of the ongoing search for species-specific features that may aid the development of effective control strategies and thus reduce the burden on patients and healthcare costs.

Antimicrobial Resistance and Biofilms Underlying Catheter-Related Bloodstream Coinfection by Enterobacter cloacae Complex and Candida parapsilosis

Biofilm-associated infections are a public health concern especially in the context of healthcare-associated infections such as catheter-related bloodstream infections (CRBSIs). We evaluated the biofilm formation and antimicrobials resistance (AMR) of Enterobacter cloacae complex and Candida parapsilosis co-isolated from a CRBSI patient. Antimicrobial susceptibility of central venous catheters (CVCs) and hemoculture (HC) isolates was evaluated, including whole genome sequencing (WGS) resistome analysis and evaluation of gene expression to obtain insight into their AMR determinants. Crystal violet assay was used to assess dual biofilm biomass and microscopy was used to elucidate a microorganism’s distribution within biofilms assembled on different materials. Bacteria were multidrug-resistant including resistance to colistin and beta-lactams, likely linked to the mcr-9-like phosphoethanolamine transferase and to an ACT family cephalosporin-hydrolyzing class C beta-lactamase, respectively. The R398I and Y132F mutations in the ERG11 gene and its differential expression might account for C. parapsilosis resistance to fluconazole. The phenotype of dual biofilms assembled on glass, polystyrene and polyurethane depends on the material and how biofilms were initiated by one or both pathogens. Biofilms assembled on polyurethane were denser and richer in the extracellular polymeric matrix, and microorganisms were differently distributed on the inner/outer surface of the CVC.

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

OBJECTIVES

Hereby is described the molecular mechanisms underlying the acquisition of azole resistance by a C. parapsilosis isolate following fluconazole treatment due to candiduria.

METHODS

A set of three consecutive C. parapsilosis isolates were recovered from urine samples of a patient with candiduria. Whole-genome sequencing (WGS) and antifungal susceptibility assays were performed. 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 2^nd cycle of treatment, exhibited a susceptible-dose dependent phenotype to fluconazole, while isolate CPS-C, recovered after the 3rd cycle, exhibited a cross-resistance profile to fluconazole and voriconazole. WGS 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 (MRR1_RI) resulted in resistance to fluconazole and voriconazole, as well as upregulation of MRR1 and MDR1. Interestingly, the susceptible-dose dependent phenotype exhibited by isolate CPS-B is associated with an increased copy number of the CDR1B gene. Expression of CDR1B is increased in both isolates CPS-B and CPS-C, and in the MRR1_RI strain, harboring the gain-of-function (GOF) mutation.

CONCLUSIONS

Our results describe clinical azole cross-resistance acquisition in C. parapsilosis due to a G1810A (G604R) GOF mutation resulting in MRR1 hyperactivation and consequently, MDR1 efflux pump overexpression. We also associated amplification of CDR1B gene with decreased fluconazole susceptibility and showed that it is a putative target of the MRR1 GOF mutation.