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.

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.

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.

Delineation of the Direct Contribution of Candida auris ERG11 Mutations to Clinical Triazole Resistance

Resistance to fluconazole is one of clinical characteristics most frequently challenging the treatment of invasive Candida auris infections, and is observed among >90% of all characterized clinical isolates. In this work, the native C. auris ERG11 allele in a previously characterized fluconazole-susceptible clinical isolate was replaced with the ERG11 alleles from three highly fluconazole-resistant clinical isolates (MIC ≥256 mg/L), encoding the amino acid substitutions VF125AL, Y132F, and K143R, using Cas9-ribonucleoprotein (RNP) mediated transformation system. Reciprocally, the ERG11WT allele from the same fluconazole-susceptible clinical isolate, lacking any resistance-associated mutation, was introduced into a previously characterized fluconazole-resistant clinical isolate, replacing the native ERG11K143R allele, using the same methods. The resulting collection of strains was subjected to comprehensive triazole susceptibility testing, and the direct impact each of these clinically-derived ERG11 mutations on triazole MIC was determined. Introduction of each of the three mutant ERG11 alleles was observed to increase fluconazole and voriconazole MIC by 8- to 16-fold. The MIC for the other clinically available triazoles were not significantly impacted by any ERG11 mutation. In the fluconazole-resistant clinical isolate background, correction of the K143R encoding mutation led to a similar 16-fold decrease in fluconazole MIC, and 8-fold decrease in voriconazole MIC, while the MIC of other triazoles were minimally changed. Taken together, these findings demonstrate that mutations in C. auris ERG11 significantly contribute to fluconazole and voriconazole resistance, but alone cannot explain the substantially elevated MIC observed among clinical isolates of C. auris. IMPORTANCE Candida auris is an emerging multidrug-resistant and health care-associated pathogen of urgent clinical concern. The triazoles are the most widely prescribed antifungal agents worldwide and are commonly utilized for the treatment of invasive Candida infections. Greater than 90% of all C. auris clinical isolates are observed to be resistant to fluconazole, and nearly all fluconazole-resistant isolates of C. auris are found to have one of three mutations (encoding VF125AL, Y132F, or K143R) in the gene encoding the target of the triazoles, ERG11. However, the direct contribution of these mutations in ERG11 to fluconazole resistance and the impact these mutations may have the susceptibility of the other triazoles remains unknown. The present study seeks to address this knowledge gap and potentially inform the future application the triazole antifungals for the treatment of infections caused by C. auris.