Emergence and circulation of azole-resistant C. albicans, C. auris and C. parapsilosis bloodstream isolates carrying Y132F, K143R or T220L Erg11p substitutions in Colombia

Methods: Over a four-year period, 123 Candida bloodstream isolates were collected at a quaternary care hospital. The isolates were identified by MALDI-TOF MS and their fluconazole (FLC) susceptibility patterns were assessed according to CLSI guidelines. Subsequently, sequencing of ERG11, TAC1 or MRR1, and efflux pump activity were performed for resistant isolates.

Results: Out of 123 clinical strains,C. albicans accounted for 37.4%, followed by C. tropicalis 26.8%, C. parapsilosis 19.5%, C. auris 8.1%, C. glabrata 4.1%, C. krusei 2.4% and C. lusitaniae 1.6%. Resistance to FLC reached 18%; in addition, a high proportion of isolates were cross-resistant to voriconazole. Erg11 amino acid substitutions associated with FLC-resistance (Y132F, K143R, or T220L) were found in 11/19 (58%) of FLCresistant isolates. Furthermore, novel mutations were found in all genes evaluated. Regarding efflux pumps, 8/19 (42%) of FLC-resistant Candida spp strains showed significant efflux activity. Finally, 6/19 (31%) of FLC-resistant isolates neither harbored resistance-associated mutations nor showed efflux pump activity. Among FLC-resistant species, C. auris 7/10 (70%) and C. parapsilosis 6/24 (25%) displayed the highest percentages of resistance (C. albicans 6/46, 13%).

Discussion: Overall, 68% of FLC-resistant isolates exhibited a mechanism that could explain their phenotype (e.g. mutations, efflux pump activity, or both). We provide evidence that isolates from patients admitted to a Colombian hospital harbor amino acid substitutions related to resistance to one of the most commonly used molecules in the hospital setting, with Y132F being the most frequently detected.

1 Introduction

Invasive fungal infections (IFIs) due to Candida species are a frequent and life-threatening condition in hospital settings worldwide, and are often associated with high morbidity and mortality (Koehler et al., 2019). Candida albicans is the most isolated species, with close to a 30% mortality rate in candidemia (Pappas et al., 2018). However, non-albicans Candida species (NACS) such as C. glabrata and C. parapsilosis have emerged as a common cause, becoming the second or third most frequent species depending on geography, patient underlying condition, and age. For these species, the associated mortality rate is about 50% and 28%, respectively (Pemán et al., 2005; Lamoth et al., 2018). Another important NACS causing candidemia is C. auris, whose mortality rate ranges from 40% to 60% according to some studies, despite the fact its prevalence is unclear (Chen et al., 2020).

While current therapeutic options for IFIs are limited to only three classes of drugs (i.e., polyenes, azoles and echinocandins), the emergence of resistant strains to some of these molecules is even more concerning. For decades, azoles have been the most frequently used antifungal for treating Candida infections (Sheehan et al., 1999; Carradori et al., 2022). Although treatment with azoles can be effective, long-term use of fluconazole (FLC) has led to the emergence of Candida strains with decreased susceptibility. In the case of Candida spp, the molecular mechanisms behind FLC-resistance have been relatively well characterized (Morio et al., 2017; Carolus et al., 2021).

Unlike C. auris which is often resistant to FLC, C. albicans and C. parapsilosis isolates were thought to be universally susceptible to FLC, but recent studies show an increased resistance. For example, a multicenter laboratory-based survey of candidemia conducted in South Africa indicates that more than half of C. parapsilosis isolates 62% (332/531) are resistant to FLC (Govender et al., 2016). In addition, studies in Brazil, India, Kuwait, South Korea, Spain, Turkey, and the United States, as well as a recent global study, confirmed the emergence of FLC-resistance in C. parapsilosis (Berkow et al., 2015; Asadzadeh et al., 2017; Arastehfar et al., 2020; Castanheira et al., 2020; Díaz-García et al., 2022). Regarding C. albicans, this species exhibits lower levels of azole resistance. However, resistant isolates have been reported from many countries around the world, including Colombia (Costa-de-Oliveira and Rodrigues, 2020; Rojas et al., 2020; Ceballos-Garzon et al., 2020b).

The main FLC-resistance mechanisms are associated with i) up-regulation of drug transporters, ii) alteration or up-regulation of the gene encoding the enzyme being targeted, which decreases binding affinity for the drug and increases concentration of the enzyme target, iii) alterations in the ergosterol synthetic pathway and iv) activation of pathways involved in the stress response, such as the Ras/cAMP/PKA pathway, calmodulin/calcineurin pathway (CaM/CaL), and mitogen-activated protein kinase (MAPK) signaling pathways (Shapiro et al., 2011; Pristov and Ghannoum, 2019). However, the most studied of them are described below.

The ATP-binding cassette (ABC) and the major facilitator superfamily (MFS) transporters are responsible for lowering the accumulation of azoles inside the yeast cell by translocating compounds actively across the cell membrane (Costa-de-Oliveira and Rodrigues, 2020). Overexpression of genes encoding drug transporters, e.g., Cdr1/2-ABC and Mdr1/Flu1-MFS, among resistant isolates of Candida species is predominantly due to gain-of-function (GOF) mutations in genes encoding zinc cluster transcription factors, such as TAC1 (transcriptional activator of CDR genes) and MRR1 (multidrug resistance regulator). For instance, GOFs in TAC1 (T225A, R693K, A736V, H741, N972D, G980E, N997D) and MRR1 (I283R, R479K, G583R, V854A, K873N), lead to overexpression of CDR1 and MDR1 in C. albicans and C. parapsilosis, respectively (Coste et al., 2006; Liu et al., 2020).

Overexpression of ERG11, the gene encoding lanosterol 14α-demethylase, the azole target, contributes directly to resistance as the increased abundance of the target requires higher drug doses for inhibition. Activating mutations in the gene encoding the transcription factor Upc2, which up-regulates most ergosterol biosynthesis genes, and the formation of an isochromosome with two copies of the left arm of chromosome 5 [i(5L)], or by duplication of the whole chromosome, on which ERG11 resides, are responsible for ERG11 overexpression (Cowen et al., 2014). Furthermore, point mutations in the ERG11 alter the 3D conformation of Erg11 and reduce its affinity for FLC. Some of the most frequent amino acid substitutions reported are Y132F and K143R substitutions, described in C. albicans, C. parapsilosis and C. auris (Berkow et al., 2015; Flowers et al., 2015; Chow et al., 2020).

Although Candida isolates from individual institutions may not be representative of the data of a country, such studies can provide a useful baseline snapshot of species distribution and antifungal susceptibility for candidemia in resource-limited settings (Govender et al., 2016). In Colombia, there is a lack of data about antifungal resistance and its molecular mechanisms in Candida spp. Therefore, this study aimed to investigate the prevalence of resistance and to describe the mechanisms behind FLC-resistance in a collection of bloodstream isolates.

2 Materials and methods

2.1 Ethics approval and consent to participate

The research and ethics committee of the Hospital Universitario San Ignacio (HUSI) approved this study (no. FM-CIE-8053-14). All patients are anonymized and only the code of isolates was transferred for this investigation. Therefore, no informed consent was required.

2.2 Study design

The study was a single-center retrospective analysis. One hundred twenty-three bloodstream isolates of Candida spp. obtained from 123 hospitalized patients (2016-2020) of the San Ignacio Hospital in Bogota, Colombia were included. Prior to storage at −80°C, yeast from blood cultures submitted for routine work-up to the Clinical Microbiology Laboratory were primarily identified using MicroScan (MicroScan WalkAway-96 Plus, Siemens, Deerfield, IL, USA) or VITEK 2 system (bioMérieux, Marcy-l’Etoile, France), and further characterized (this study) using the MALDI-TOF Biotyper system (Bruker Daltonik, Bremen, Germany).

2.3 MALDI-TOF MS

Isolates were streaked from a glycerol stock onto Sabouraud dextrose agar (SDA) and grown for 24–36 h at 35°C. Protein extraction was performed using formic acid/ethanol method, according to the Bruker Daltonics’ protocol. The protein mass spectra were analysed using the Flex Control software and the MALDI Biotyper version 3.1 7311 reference spectra (main spectra) (Bruker Daltonics, Bremen, Germany). MALDI-TOF MS results were obtained according to the manufacturer’s technical specifications, as follows: correct genus and species identification (≥2.0), correct genus identification (1.7–2.0), and no reliable identification (< 1.7). All clinical isolates had a score above 2.0 (Ceballos-Garzon et al., 2020a).

2.4 Antifungal susceptibility testing

Susceptibility to FLC (Sigma-Aldrich, St. Louis, MO, USA) was conducted using the Clinical and Laboratory Standards Institute broth microdilution method (CLSI-BMD), following the M27-A3 document (CLSI, 2008). Quality control was ensured by testing the CLSI-recommended strains C. parapsilosis ATCC 22019 and C. krusei ATCC 6258. For NACS isolates the CLSI breakpoints were applied (resistance to FLC was set at C. albicans, C. tropicalis, and C. parapsilosis ≥8µg/mL and C. glabrata ≥64µg/mL) (CLSI, 2018). In the case of C. auris, the FLC breakpoint recommended by the US Centers for Disease Control and Prevention (CDC) was used (≥32 µg/mL) (CDC, 2020). The MIC data obtained under routine conditions for amphotericin B (AMB), caspofungin (CAS), itraconazole (ITC), and voriconazole (VRC) by Etest^® (bioMérieux, Marcy-l’Étoile, France) and VITEK^®2 (bioMérieux) are presented in Table S1.

2.5 Sequencing analysis of Erg11, Tac1 and Mrr1-encoding genes

All isolates displaying resistance to FLC, and one susceptible isolate of each species were subjected to a single-tube PCR method to amplify and sequence the coding region of the ERG11, TAC1 or MRR1 genes (both strands) using the primers indicated in Table S2. The PCR products were purified and sequenced using a SeqStudio genetic analyzer capillary sequencer (Applied Biosystems). The sequencing results were analyzed by BLAST and compared with the published GenBank sequences: C. albicans AY856352.1 (ERG11), DQ393587 (TAC1), and C. parapsilosis GQ302972 (ERG11), HE605205 (MRR1). For C. auris, sequences download from the Candida genome database (candidagenome.org) were used, i.e., B9J08_001448 (C. auris B8441, ERG11) and B9J08_004820 (TAC1b). All sequences (FLC-susceptible and resistant clinical isolates plus reference strains) were aligned, and the dataset was used to construct a Neighbor-Joining phylogenetic tree using Maximum Composite Likelihood settings by using Molecular Evolutionary Genetics Analysis Version 11 (MEGA11) (Saitou and Nei, 1987; Tamura et al., 2021). Codon positions included were 1st + 2nd + 3rd + Noncoding. All positions containing gaps and missing data were eliminated. Evaluation of branch support was performed by Bootstrap statistical analysis with 1000 replicates (Tamura et al., 2021).

2.6 Analysis of rhodamine 6G efflux

The accumulation of R6G in growing Candida cells correlates inversely with the mRNA expression level of the ABC transporter Candida drug resistance 1 (CDR1), therefore the levels of intracellular accumulation of R6G can be used for the identification of azole-resistant strains. ABC transporter-mediated efflux was determined using rhodamine 6G (Sigma-Aldrich, USA) as previously described by Gbelska and coworkers in a step-by-step protocol (Gbelska et al., 2017). The fluorescence of the released R6G was measured at 530 nm, with an emission at 560 nm in an automated plate reader (Model 550 Microplate Reader Bio-Rad, Milan, Italy). Measurements were made before (basal) and after the addition of 20 mM glucose. Using a R6G calibration curve, the fluorescence intensity was converted into concentration.

2.7 Statistics

Antifungal susceptibility and 6G efflux were performed in triplicate, in three independent experimental sets. The results of the efflux pump activity were analyzed statistically by the Analysis of Variance One-Way ANOVA using GraphPad Prism version 9 software. In all analyses, p values of 0.05 or less were considered statistically significant.

3 Results

3.1 Identification and Antifungal susceptibility testing of clinical isolates

The Candida species distribution from the 123 blood samples was as follows: C. albicans, 46 (37.4%); C. tropicalis, 33 (26.8%); C. parapsilosis, 24 (19.5%); C. auris, 10 (8.1%); C. glabrata, 5 (4.1%); C. krusei, 3 (2.4%); and C. lusitaniae, 2 (1.6%). Although C. albicans was the most prevalent species, accounting for 37.4%, the NACS group comprised 62.6% of the isolates identified. Concerning susceptibility, Table 1 shows the geometric MICs, the MIC ranges, the MIC50, and MIC90 distributions of FLC against 123 Candida spp strains. The MIC₉₀ values (MICs at which ≥90% of strains are inhibited) for the four most frequent species found were: C. albicans 128µg/mL, C. tropicalis 2µg/mL, C. parapsilosis 32µg/mL and C. auris >128µg/mL. As expected, C. auris presented the highest MIC₉₀ values. In addition, the range of FLC MICs was narrower for C. auris and C. glabrata than for the other species. Both, C. auris and C. glabrata had the highest geometric mean values: 39 and 42 respectively. On the other hand, C. tropicalis showed reduced susceptibility to FLC (MIC_{90 =} 2µg/mL). When the CLSI and CDC breakpoints were applied, 22 out of 123 (18%) isolates displayed in vitro resistance to FLC, among them, six of C. albicans, six of C. parapsilosis and seven of C. auris (19/22). Moreover, three isolates of C. krusei -as C. krusei is assumed to be intrinsically resistant to FLC-, were not included within the molecular study. In contrast, all isolates of C. tropicalis, C. glabrata, and C. lusitaniae were FLC-susceptible Figure 1. Regarding VRC MICs (data obtained by Etest), 10 out of 19 FLC-resistant isolates, excluding C. krusei, were VRC cross-resistant. The highest MICs were observed for C albicans (32 µg/mL) Table S3. Considering the MIC data obtained under routine conditions for AMB (mean: 0.3 µg/mL) and CAS (mean 0.06 µg/mL) all strains showed low MICs values (susceptible) Table S1.

TABLE 1

Table 1 Antifungal activity of fluconazole drug against Candida spp (n = 123) performed by CLSI.

FIGURE 1

Figure 1 Species distribution and overall fluconazole susceptibility results. (A) Distribution of the 123 species identified. (B) Percentage of fluconazole resistance. The number of isolates is indicated inside the bars. Susceptible (S) isolates are depicted in blue and resistant (R) isolates in red.

3.2 Detection of mutations in Erg11, Tac1 and Mrr1-encoding genes of 19 FLC-resistant isolates

By comparing the ERG11 coding region of C. albicans (CAAL, A-F), C. auris (CAAU, A-G) and C. parapsilosis (CAPA, A-F) FLC-resistant isolates with that of our FLC-susceptible and the published wild-type sequences, we identified 20, 16, and three mutations, respectively. As expected, some silent mutations that do not change the protein sequence were identified (data not shown). The remaining ERG11 mutations which resulted in amino acid changes are shown in, Table 2, Figure 2. Among the nonsense mutations (C. albicans 6; C. auris 5; and C. parapsilosis 2), three amino acid substitutions related with FLC-resistance (T220L, Y132F, K143R) were found from which Y132F was the most detected. Additionally, five amino acid substitutions previously described in FLC-susceptible isolates were observed (D116E, K128T, K177R, N335S, E343D). To the best of our knowledge, four amino acid substitutions (K22E, Q38T, F72V, Q77S) have not been previously reported. Overall, eight of the 19 FLC-resistant isolates did not have resistance-associated substitutions.

TABLE 2

Table 2 Amino acid substitutions found in 19 FLC-resistant isolates.

FIGURE 2

Figure 2 Phylogenetic tree and amino acid substitutions found in fluconazole-resistant isolates of (A) C albicans (CAAL), (B) C auris (CAAU) and (C) C parapsilosis (CAPA). aa., amino acids; FLC, fluconazole; MIC, minimal inhibitory concentration; CA-S., Candida FLC-susceptible. Inside the colored squares, the number of substitutions found is indicated.

On the other hand, except for isolate CAAL-A, which harbored a K128T substitution in only one ERG11 allele, other isolates were homozygous for mutations in the ERG11 allele. In the phylogenetic relationship among FLC-resistant, susceptible isolates and reference strains, a cluster of isolates carrying the substitutions Y132F in C. auris was observed, as well as in the FLC-resistant C. albicans isolates harboring K128T substitution. In addition, susceptible and resistant isolates without resistance-associated ERG11 mutations (i.e., T220L, Y132F) from these species were clustered Figures 2A, B.

Concerning TAC1 (C. albicans and C. auris) and MRR1 (C. parapsilosis) genes, no mutations previously associated with FLC-resistance were found. However, 17 Tac1 (C. albicans, n=7; C. auris, n=10) and Mrr1 (C. parapsilosis, n=1) amino acid substitutions previously described in FLC-susceptible isolates were found. Additionally, there were 37 unreported substitutions in Tac1 (C. albicans, n=7; C. auris, n=26) and Mrr1 (C. parapsilosis, n=4) Table 2.

3.3 Efflux pumps activity

To gain further insights into the mechanisms of azole resistance in the clinical isolates the activity of efflux pumps was evaluated using rhodamine 6G, which uses the same membrane ABC transporters (Cdr1p and Cdr2p) as FLC in Candida. Among the 19 FLC-resistant isolates, eight of them showed significant active efflux of rhodamine-6G after addition of glucose: two of them belonging to C. albicans (CAAL-C and CAAL-E); four to C. auris (CAAU-A, CAAU-B, CAAU-C and CAAU-G) and two to C. parapsilosis (CAPA-C, CAPA-F) Figure 3, Figure 1S.

FIGURE 3

Figure 3 Rhodamine 6G (R6G) efflux over time among fluconazole-resistant isolates. Colored lines indicate the concentration of R6G released after the addition of 20 mM glucose (GLU). Data are means ± SD from three experiments. *p < 0.05, **p < 0.01, ***p < 0.001. ****p≤ 0.0001. FLC, fluconazole; MIC, minimal inhibitory concentration.

Interestingly, no amino acid substitutions associated with FLC-resistance were identified in both C. albicans isolates (CAAL-C, 128µg/mL), (CAAL-D, 128µg/mL). Regarding the remaining isolates that showed efflux pump activity, C. auris (CAAU-A, CAAU-B, CAAU-C, CAAU-G) and C. parapsilosis (CAPA-C, CAPA-F) harbored the Y132F substitution.

Overall, 11/19 resistant isolates harbored an amino acid substitution associated with FLC-resistance; 8/19 displayed efflux pumps activity; 6/19 had both amino acid substitutions and efflux pump activity; and 6/19 isolates did not exhibit any of the mechanisms that could explain their resistant phenotype Figure 4.

FIGURE 4

Figure 4 Schematic representation of the results obtained in this study for the 19 FLC-resistant isolates. In purple C. albicans, light blue C. parapsilosis and yellow C. auris isolates with their respective MICs against FLC. In red font, the mutations found. The green stars indicate isolates with significant efflux pump activity. The size of the circle reflects the number of mutations, those not associated with FLC resistance (blue) and the new mutations (light green).

4 Discussion

In the region, a study that evaluated susceptibility of Candida species identified in Colombia, Ecuador, and Venezuela, found a percentage of resistance to FLC of 6.8% (De bedout et al., 2003). This shows that there is an important change in the rate of resistance, which in our study reached 18%. Nevertheless, the majority of C. albicans isolates were azole susceptible, thus the observed resistance percentage is mainly attributed to the presence of C. auris, which agrees with previous studies (Pfaller et al., 2019; Sathi et al., 2022). Furthermore, it is noteworthy that 25% of the C. parapsilosis isolates were resistant to FLC, which confirms the increase of resistance in this species (Escribano and Guinea, 2022) (see Figure 5). High azole resistance rates have been reported for this species in other single-center studies conducted in Brazil (67.9%), Italy (33%), France (9.2%), Mexico (54%), Saudi Arabia (33%), Spain (13.6%), South Africa (78%) or Turkey (26.4%) (Arastehfar et al., 2020; Magobo et al., 2020; Mesini et al., 2020; Corzo-Leon et al., 2021; Fekkar et al., 2021; Thomaz et al., 2021; Díaz-García et al., 2022).

FIGURE 5

Figure 5 Countries from which azole-resistant C. parapsilosis carrying ERG11 mutations have been reported, as of February, 2023.

Overall, the global data indicate that C. albicans remains the predominant species identified in Candida infections. In this study, C. albicans was the most predominant species, nevertheless, NACS accounted for a high proportion (62.6%). Although reported resistance rates vary from study to study, the surveillance data collected suggest that azole resistance rates for C. albicans remain low (Pfaller et al., 2019; Nishimoto et al., 2020). However, C. albicans infections in the bloodstream pose a considerable threat in immunocompromised populations and the associated high mortality remains a major problem in the clinical setting. Therefore, C. albicans should not be overlooked as a serious public health threat. Remarkably, in this study, the percentage of FLC resistance in C. albicans was 13%, which is relatively high.

The acquisition of azole resistance is a serious concern given the limited number of molecules available for the treatment of IFIs. Moreover, this is much more worrying in resource-limited regions where FLC is the only available therapy (Kneale et al., 2016). As previously described, azole resistance is mainly conferred by mutations in the ERG11 gene and by the activity of efflux pumps. The ERG11 gene is highly polymorphic and more than 140 amino acid substitutions have been reported, indicating that this protein is very permissive to conformational changes. Most substitutions occur in three amino acid hotspot regions (105-165, 266-287 and 405-488), although mutations outside these regions can also be found (Debnath and Addya, 2014; Oliveira-Carvalho and Del Negro, 2014).

In the present study, 11 of the 19 FLC-resistant isolates harbored mutations which have been previously described in resistant isolates (T220L, Y132F, K143R) (Healey et al., 2018; Arastehfar et al., 2020). Notably, the T220L substitution was observed in one C. albicans isolate. In C. auris, five isolates carried the Y132F, and in the case of C. parapsilosis four isolates harbored the Y132F and one the K143R substitution.

The substitution of lysine for threonine at position 128 (K128T) was found in CAAL-A, CAAL-B and CAAL-D isolates, which is consistent with the findings of Cernicka et al., and Peron et al., who identified this amino acid substitution in FLC-resistant strains (Cernicka and Subik, 2006; Peron et al., 2016). However, several studies refute this association because K128T substitution has been found in multiple FLC-susceptible isolates, hence, in this study, we do not consider it as a mutation associated with resistance (Morio et al., 2010; Flowers et al., 2015). Nevertheless, it might influence translation efficiency, leading to alterations in protein production as it occurs with nearby mutations such as G129R and Y132F (Kumar et al., 2020).

According to Chow and coworkers, mutations contributing to FLC resistance are clade-specific in C. auris, being Y132F and K143R the most predominant in clade I, F126L in clade III, and Y132F in clade IV. However, Y132F in Erg11 is the most prevalent one (Chow et al., 2020). As described for isolates of clade IV (South America), we found only one mutation (i.e., Y132F). Concerning C. parapsilosis, high prevalence of the Y132F substitution was also noted in previous studies (Grossman et al., 2015; Choi et al., 2018). In addition, the latter suggested that isolates harbouring this mutation may have a higher propensity to cause clonal transmission and to persist in nosocomial settings (Thomaz et al., 2018). Regarding our isolates, the Y132F substitution was detected in four out of six C. parapsilosis isolates. However, considering the dates of collection of the isolates, It does not seem to be a clonal spread (Table S3). Although this report is due in 2023, some resistant isolates were obtained in 2015. Therefore, prior to the reports made in some parts of the globe (Figure 5).

Although previous reports indicate that the Mrr1 substitutions I283R, R479K, G583R, A854V, K873N and L986P are associated with FLC and/or VRC resistance, none of our isolates harbored them (Branco et al., 2015). Similar to Tac1, residues located near the C terminus (760, 761, 803, 956, and 966) might contribute to azole resistance (Nishimoto et al., 2020).

This study also illustrates that acquired azole resistance commonly relies on combined molecular mechanisms in clinical isolates (Morio et al., 2013). In addition to amino acid substitutions in Erg11, six of the 19 isolates also displayed active efflux activity. Interestingly, in two isolates lacking resistance-related mutations, efflux pumps activity was observed. Considering that efflux pumps play a key role in azole resistance, this might be the explanation for the observed phenotype. However, a further study evaluating the expression of all genes involved in efflux pump activity is required.

Regarding the isolates in which no mechanism for resistance was found, namely CAAL-A,B,D, CAAU-D, CAAU-D,E and CAPA-A, six novel Tac1 or Mrr1 substitutions (G52R, V366del, F682T, F683L, T695S and P229A) were identified in three of the six isolates. Moreover, four novel Erg11 substitutions (K22E, Q38T, F72V, Q77S) in the FLC resistant C. albicans (CAAL-B, CAAL-C) and C. auris (CAAU-B) were found. The putative role of these substitutions remains to be investigated. In isolates CAAL-A-D and CAAU-E, a different mechanism to those evaluated here should confer resistance; further analysis is ongoing.

Our study has some limitations. For instance, the antifungal susceptibility data obtained for antifungals other than FLC are not complete and were not performed by broth microdilution. The molecular basis of FLC resistance was investigated only by analysis of the ERG11, TAC1, and MRR1 genes, while other mechanisms conferring resistance, such as sterol composition and gene expression, were not investigated due to lack of funds. Finally, the clinical description of the patients was missing and will be reported alongside other patients infected by Candida antifungal resistant strains elsewhere.

Recently, the WHO released the fungal priority pathogens list (WHO FPPL) which includes C. albicans, C auris (Critical group), and C. parapsilosis (high group). The WHO describes that to overcome the lack of knowledge on infections caused by these fungi, more data and evidence on fungal infections and antifungal resistance to inform and improve response to FPP is needed (World Health Organization (WHO), 2022). Although our study does not have a large number of isolates and does not include other healthcare institutions, we provide evidence of Colombian isolates harboring resistance mutations to one of the most used molecules in the hospital setting.

In conclusion, we describe for the first time the presence of Candida spp isolates harboring Erg11 FLC resistant-related substitutions (Y132F, K143R and T220L) in patients admitted to a Colombian hospital. Although FLC-resistance rates differ significantly between countries and individual health facilities, the resistance rate in our study is relatively high, emphasizing the need for active surveillance to prevent further expansion of FLC-resistant Candida spp isolates in the clinical setting.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

AP carried out the experiments. YV-C contributed to the development of efflux pumps protocol. AC-G, and AP analyzed the data. AC-G wrote the main manuscript. AP, SV-B, BA, and CP-G review and editing. and CP-G conceived the experiments and managed the resources. All authors contributed to the article and approved the submitted version.

Funding

The research office of Hospital Universitario San Ignacio and the Vice-Rectory of research at the Pontificia Universidad Javeriana in Bogotá, Colombia, supported the research (grants no.2014-52 and 20454).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1136217/full#supplementary-material

References

Arastehfar, A., Daneshnia, F., Hilmioglu-Polat, S., Fang, W., Yaşar, M., Polat, F., et al. (2020). First report of candidemia clonal outbreak caused by emerging fluconazole-resistant candida parapsilosis isolates harboring Y132F and/or Y132F+K143R in Turkey. Antimicrob. Agents Chemother. 64, e01001-20. doi: 10.1128/AAC.01001-20/SUPPL_FILE/AAC.01001-20-S0001.PDF

PubMed Abstract | CrossRef Full Text | Google Scholar

Asadzadeh, M., Ahmad, S., Al-Sweih, N., Khan, Z. (2017). Epidemiology and molecular basis of resistance to fluconazole among clinical candida parapsilosis isolates in Kuwait. Microb. Drug Resist. 23, 966–972. doi: 10.1089/MDR.2016.0336

PubMed Abstract | CrossRef Full Text | Google Scholar

Berkow, E. L., Manigaba, K., Parker, J. E., Barker, K. S., Kelly, S. L., Rogers, P. D. (2015). Multidrug transporters and alterations in sterol biosynthesis contribute to azole antifungal resistance in candida parapsilosis. Antimicrob. Agents Chemother. 59, 5942–5950. doi: 10.1128/AAC.01358-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Branco, J., Silva, A. P., Silva, R. M., Silva-Dias, A., Pina-Vaz, C., Butler, G., et al. (2015). Fluconazole and voriconazole resistance in candida parapsilosis is conferred by gain-of-function mutations in MRR1 transcription factor gene. Antimicrob. Agents Chemother. 59, 6629-33. doi: 10.1128/AAC.00842-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Carolus, H., Pierson, S., Muñoz, J. F., Subotić, A., Cruz, R. B., Cuomo, C. A., et al. (2021). Genome-wide analysis of experimentally evolved candida auris reveals multiple novel mechanisms of multidrug resistance. MBio 12, e03333-20. doi: 10.1128/mBio.03333-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Carradori, S., Ammazzalorso, A., De Filippis, B., Şahin, A. F., Akdemir, A., Orekhova, A., et al. (2022). Azole-based compounds that are active against candida biofilm: In vitro, In vivo and in silico studies. Antibiot 11, 1375. doi: 10.3390/ANTIBIOTICS11101375

CrossRef Full Text | Google Scholar

Castanheira, M., Deshpande, L. M., Messer, S. A., Rhomberg, P. R., Pfaller, M. A. (2020). Analysis of global antifungal surveillance results reveals predominance of Erg11 Y132F alteration among azole-resistant candida parapsilosis and candida tropicalis and country-specific isolate dissemination. Int. J. Antimicrob. Agents 55, 105799. doi: 10.1016/J.IJANTIMICAG.2019.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

CDC (2020). Centers for disease control and prevention. antifungal susceptibility testing and interpretation to laboratorians and health professionals. Available at: https://www.cdc.gov/fungal/candida-auris/c-auris-antifungal.htmlhttps://www.cdc.gov/fungal/candida-auris/ (Accessed July 10, 2021).

Google Scholar

Ceballos-Garzon, A., Amado, D., Vélez, N., Jiménez-A, M. J., Rodríguez, C., Parra-Giraldo, C. M. (2020a). Development and validation of an in-house library of Colombian candida auris strains with MALDI-TOF MS to improve yeast identification. J. Fungi 6, 1–9. doi: 10.3390/jof6020072

CrossRef Full Text | Google Scholar

Ceballos-Garzon, A., Wintaco-Martínez, L. M., Velez, N., Hernandez-Padilla, C., de la Hoz, A., Valderrama-Beltran, S. L., et al. (2020b). Persistence of clonal azole-resistant isolates of candida albicans from a patient with chronic mucocutaneous candidiasis in Colombia. J. Glob. Infect. Dis. 12, 16–20. doi: 10.4103/jgid.jgid_74_19

PubMed Abstract | CrossRef Full Text | Google Scholar

Cernicka, J., Subik, J. (2006). Resistance mechanisms in fluconazole-resistant candida albicans isolates from vaginal candidiasis. Int. J. Antimicrob. Agents 27, 403–408. doi: 10.1016/J.IJANTIMICAG.2005.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Tian, S., Han, X., Chu, Y., Wang, Q., Zhou, B., et al. (2020). Is the superbug fungus really so scary? a systematic review and meta-analysis of global epidemiology and mortality of candida auris. BMC Infect. Dis. 20, 827. doi: 10.1186/s12879-020-05543-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, Y. J., Kim, Y. J., Yong, D., Byun, J. H., Kim, T. S., Chang, Y. S., et al. (2018). Fluconazole-resistant candida parapsilosis bloodstream isolates with Y132F mutation in ERG11 gene, south Korea. Emerg. Infect. Dis. 24, 1768. doi: 10.3201/EID2409.180625

PubMed Abstract | CrossRef Full Text | Google Scholar

Chow, N. A., Muñoz, J. F., Gade, L., Berkow, E. L., Li, X., Welsh, R. M., et al. (2020). Tracing the evolutionary history and global expansion of candida auris using population genomic analyses. MBio 11, e03364-19. doi: 10.1128/mBio.03364-19

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI (2008). M27-A3 reference method for broth dilution antifungal susceptibility testing of yeasts–third edition. Clin. Lab. Stand. Inst. 37, 1–2.

Google Scholar

CLSI (2018). M60ed1e: Performance standards for antifungal susceptibility testing of yeasts. 2nd ed (Wayne, PA, USA: Clinical and Laboratory Standards Institute).

Google Scholar

Corzo-Leon, D. E., Peacock, M., Rodriguez-Zulueta, P., Salazar-Tamayo, G. J., MacCallum, D. M. (2021). General hospital outbreak of invasive candidiasis due to azole-resistant candida parapsilosis associated with an Erg11 Y132F mutation. Med. Mycol. 59, 664–671. doi: 10.1093/MMY/MYAA098

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa-de-Oliveira, S., Rodrigues, A. G. (2020). Candida albicans antifungal resistance and tolerance in bloodstream infections: The triad yeast-Host-Antifungal. microorganisms 8, 1–19. doi: 10.3390/microorganisms8020154

CrossRef Full Text | Google Scholar

Coste, A., Turner, V., Ischer, F., Morschhäuser, J., Forche, A., Selmecki, A., et al. (2006). A mutation in Tac1p, a transcription factor regulating CDR1 and CDR2, is coupled with loss of heterozygosity at chromosome 5 to mediate antifungal resistance in candida albicans. Genetics 4, 2139–2156. doi: 10.1534/genetics.105.054767

CrossRef Full Text | Google Scholar

Cowen, L. E., Sanglard, D., Howard, S. J., Rogers, P. D., Perlin, D. S. (2014). Mechanisms of antifungal drug resistance. Cold Spring Harb. Perspect. Med. 5, 19752. doi: 10.1101/CSHPERSPECT.A019752

CrossRef Full Text | Google Scholar

De bedout, C., Ayabaca, J., Vega, R., Mendez, M., Santiago, A., Pabon, M., et al. (2003). Evaluation of candida species’ susceptibility to fluconazole with the disk diffusion method. Biomedica 1, 31–37.

Google Scholar

Debnath, S., Addya, S. (2014). Structural basis for heterogeneous phenotype of ERG11 dependent azole resistance in c.albicans clinical isolates. Springerplus 3, 660. doi: 10.1186/2193-1801-3-660

PubMed Abstract | CrossRef Full Text | Google Scholar

Díaz-García, J., Gómez, A., Alcalá, L., Reigadas, E., Sánchez-Carrillo, C., Pérez-Ayala, A., et al. (2022). Evidence of fluconazole-resistant candida parapsilosis genotypes spreading across hospitals located in Madrid, Spain and harboring the Y132F ERG11p substitution. Antimicrob. Agents Chemother. 66, e0071022. doi: 10.1128/AAC.00710-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Escribano, P., Guinea, J. (2022). Fluconazole-resistant candida parapsilosis: A new emerging threat in the fungi arena. Front. Fungal Biol. 0. doi: 10.3389/FFUNB.2022.1010782

CrossRef Full Text | Google Scholar

Fekkar, A., Blaize, M., Bouglé, A., Normand, A. C., Raoelina, A., Kornblum, D., et al. (2021). Hospital outbreak of fluconazole-resistant candida parapsilosis: Arguments for clonal transmission and long-term persistence. Antimicrob. Agents Chemother. 65, e02036-20. doi: 10.1128/AAC.02036-20/SUPPL_FILE/AAC.02036-20-S000S1.PDF

CrossRef Full Text | Google Scholar

Flowers, S. A., Colón, B., Whaley, S. G., Schuler, M. A., David Rogers, P. (2015). Contribution of clinically derived mutations in ERG11 to azole resistance in candida albicans. Antimicrob. Agents Chemother. 59, 450–460. doi: 10.1128/AAC.03470-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Gbelska, Y., Toth Hervay, N., Dzugasova, V., Konecna, A. (2017). Measurement of energy-dependent rhodamine 6G efflux in yeast species. BIO-PROTOCOL 7, 1–5. doi: 10.21769/BIOPROTOC.2428

CrossRef Full Text | Google Scholar

Govender, N. P., Patel, J., Magobo, R. E., Naicker, S., Wadula, J., Whitelaw, A., et al. (2016). Emergence of azole-resistant candida parapsilosis causing bloodstream infection: results from laboratory-based sentinel surveillance in south Africa. J. Antimicrob. Chemother. 71, 1994–2004. doi: 10.1093/JAC/DKW091

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossman, N. T., Pham, C. D., Cleveland, A. A., Lockhart, S. R. (2015). Molecular mechanisms of fluconazole resistance in candida parapsilosis isolates from a U.S. surveillance system. Antimicrob. Agents Chemother. 59, 1030–1037. doi: 10.1128/AAC.04613-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Healey, K. R., Kordalewska, M., Ortigosa, C. J., Singh, A., Berrío, I., Chowdhary, A., et al. (2018). Limited ERG11 mutations identified in isolates of candida auris directly contribute to reduced azole susceptibility. Antimicrob. Agents Chemother. 62, e01427-18. doi: 10.1128/AAC.01427-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Kneale, M., Bartholomew, J. S., Davies, E., Denning, D. W. (2016). Global access to antifungal therapy and its variable cost. J. Antimicrob. Chemother. 71, 3599–3606. doi: 10.1093/JAC/DKW325

PubMed Abstract | CrossRef Full Text | Google Scholar

Koehler, P., Stecher, M., Cornely, O. A., Koehler, D., Vehreschild, M. J. G. T., Bohlius, J., et al. (2019). Morbidity and mortality of candidaemia in Europe: an epidemiologic meta-analysis. Clin. Microbiol. Infect. 25, 1200–1212. doi: 10.1016/J.CMI.2019.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, A., Nair, R., Kumar, M., Banerjee, A., Chakrabarti, A., Rudramurthy, S. M., et al. (2020). Assessment of antifungal resistance and associated molecular mechanism in candida albicans isolates from different cohorts of patients in north Indian state of haryana. Folia Microbiol. 65 (4), 747–754. doi: 10.1007/S12223-020-00785-6

CrossRef Full Text | Google Scholar

Lamoth, F., Lockhart, S. R., Berkow, E. L., Calandra, T. (2018). Changes in the epidemiological landscape of invasive candidiasis. J. Antimicrob. Chemother. 73, i4–i13. doi: 10.1093/jac/dkx444

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J.-Y., Wei, B., Wang, Y., Shi, C., Li, W.-J., Zhao, Y., et al. (2020). The H741D mutation in Tac1p contributes to the upregulation of CDR1 and CDR2 expression in candida albicans. Braz. J. Microbiol. 51, 1553–1561. doi: 10.1007/s42770-020-00336-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Magobo, R. E., Lockhart, S. R., Govender, N. P. (2020). Fluconazole-resistant candida parapsilosis strains with a Y132F substitution in the ERG11 gene causing invasive infections in a neonatal unit, south Africa. Mycoses 63, 471–477. doi: 10.1111/MYC.13070

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesini, A., Mikulska, M., Giacobbe, D. R., Del Puente, F., Gandolfo, N., Codda, G., et al. (2020). Changing epidemiology of candidaemia: Increase in fluconazole-resistant candida parapsilosis. Mycoses 63, 361–368. doi: 10.1111/MYC.13050

PubMed Abstract | CrossRef Full Text | Google Scholar

Morio, F., Jensen, R. H., Le Pape, P., Arendrup, M. C. (2017). Molecular basis of antifungal drug resistance in yeasts. Int. J. Antimicrob. Agents 50, 599–606. doi: 10.1016/j.ijantimicag.2017.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Morio, F., Loge, C., Besse, B., Hennequin, C., Le Pape, P. (2010). Screening for amino acid substitutions in the candida albicans Erg11 protein of azole-susceptible and azole-resistant clinical isolates: new substitutions and a review of the literature. Diagn. Microbiol. Infect. Dis. 66, 373–384. doi: 10.1016/j.diagmicrobio.2009.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Morio, F., Pagniez, F., Besse, M., Gay-Andrieu, F., Miegeville, M., Le Pape, P. (2013). Deciphering azole resistance mechanisms with a focus on transcription factor-encoding genes TAC1, MRR1 and UPC2 in a set of fluconazole-resistant clinical isolates of candida albicans. Int. J. Antimicrob. Agents 42, 410-5. doi: 10.1016/j.ijantimicag.2013.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishimoto, A. T., Sharma, C., Rogers, P. D. (2020). Molecular and genetic basis of azole antifungal resistance in the opportunistic pathogenic fungus candida albicans. J. Antimicrob. Chemother. 75, 257–270. doi: 10.1093/JAC/DKZ400

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveira-Carvalho, V., Del Negro, G. M. B. (2014). Is the S405F mutation in candida albicans ERG11 gene sufficient to confer resistance to fluconazole? J. Mycol. Med. 24, 241–242. doi: 10.1016/J.MYCMED.2014.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Pappas, P. G., Lionakis, M. S., Arendrup, M. C., Ostrosky-Zeichner, L., Kullberg, B. J. (2018). Invasive candidiasis. Nat. Rev. Dis. Prim. 4, 1–20. doi: 10.1038/nrdp.2018.26

CrossRef Full Text | Google Scholar

Pemán, J., Cantón, E., Gobernado, M. (2005). Epidemiology and antifungal susceptibility of candida species isolated from blood: results of a 2-year multicentre study in Spain. Eur. J. Clin. Microbiol. Infect. Dis. 24, 23–30. doi: 10.1007/S10096-004-1267-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Peron, I. H., Reichert-Lima, F., Busso-Lopes, A. F., Nagasako, C. K., Lyra, L., Moretti, M. L., et al. (2016). Resistance surveillance in candida albicans: A five-year antifungal susceptibility evaluation in a Brazilian university hospital. PloS One 11, e0158126. doi: 10.1371/journal.pone.0158126

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaller, M. A., Diekema, D. J., Turnidge, J. D., Castanheira, M., Jones, R. N. (2019). Twenty years of the SENTRY antifungal surveillance program: Results for candida species from 1997-2016. Open Forum Infect. Dis. 6, S79–S94. doi: 10.1093/ofid/ofy358

PubMed Abstract | CrossRef Full Text | Google Scholar

Pristov, K. E., Ghannoum, M. A. (2019). Resistance of candida to azoles and echinocandins worldwide. Clin. Microbiol. Infect. 25, 792–798. doi: 10.1016/J.CMI.2019.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojas, A. E., Pérez, J. E., Hernández, J. S., Zapata, Y., Rojas, A. E., Pérez, J. E., et al. (2020). Análisis cuantitativo de la expresión de genes de resistencia a fluconazol en cepas de candida albicans aisladas al ingreso de adultos mayores a una unidad de cuidados intensivos de manizales, Colombia. Biomédica 40, 153–165. doi: 10.7705/BIOMEDICA.4723

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425. doi: 10.1093/OXFORDJOURNALS.MOLBEV.A040454

PubMed Abstract | CrossRef Full Text | Google Scholar

Sathi, F. A., Paul, S. K., Ahmed, S., Alam, M. M., Nasreen, S. A., Haque, N., et al. (2022). Prevalence and antifungal susceptibility of clinically relevant candida species, identification of candida auris and kodamaea ohmeri in Bangladesh. Trop. Med. Infect. Dis. 7, 211. doi: 10.3390/TROPICALMED7090211/S1

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapiro, R. S., Robbins, N., Cowen, L. E. (2011). Regulatory circuitry governing fungal development, drug resistance, and disease. Microbiol. Mol. Biol. Rev. 75, 213–267. doi: 10.1128/MMBR.00045-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheehan, D. J., Hitchcock, C. A., Sibley, C. M. (1999). Current and emerging azole antifungal agents. Clin. Microbiol. Rev. 12, 40–79. doi: 10.1128/CMR.12.1.40

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Stecher, G., Kumar, S. (2021). MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38, 3022–3027. doi: 10.1093/MOLBEV/MSAB120

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomaz, D. Y., De Almeida, J. N., Lima, G. M. E., De Oliveira Nunes, M., Camargo, C. H., De Carvalho Grenfell, R., et al. (2018). An azole-resistant candida parapsilosis outbreak: Clonal persistence in the intensive care unit of a Brazilian teaching hospital. Front. Microbiol. 92997. doi: 10.3389/FMICB.2018.02997/BIBTEX

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomaz, D. Y., de Almeida, J. N., Sejas, O. N. E., Del Negro, G. M. B., Carvalho, G. O. M. H., Gimenes, V. M. F., et al. (2021). Environmental clonal spread of azole-resistant candida parapsilosis with Erg11-Y132F mutation causing a Large candidemia outbreak in a Brazilian cancer referral center. J. Fungi 2021 Vol. 7 Page 259 7, 259. doi: 10.3390/JOF7040259

CrossRef Full Text | Google Scholar

World Health Organization (WHO) (2022). Fungal priority pathogens list to guide research, development and public health action. Available at: https://www.who.int/publications/i/item/9789240060241 (Accessed November 3, 2022).

Google Scholar