Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 10 October 2025

Sec. Antimicrobials, Resistance and Chemotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1689176

Emergence of multidrug-resistant Candidozyma auris in southern China: a multicenter study on genetic diversity and antifungal resistance

Penghao Guo&#x;Penghao Guo1Jingchun Fang&#x;Jingchun Fang2Ruizhi WangRuizhi Wang1Weihong LinWeihong Lin3Yueting JiangYueting Jiang3Jun LongJun Long4Lingjuan ChenLingjuan Chen5Xudong HuangXudong Huang6Kang Liao
Kang Liao1*Yaqin Peng
Yaqin Peng1*
  • 1Department of Clinical Laboratory, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China
  • 2Department of Clinical Laboratory, Nansha Division of The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, Guangdong, China
  • 3Department of Clinical Laboratory, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, Guangdong, China
  • 4Department of Laboratory Medicine, Microbiome Medicine Center, Zhujiang Hospital of Southern Medical University, Guangzhou, Guangdong, China
  • 5Department of Clinical Laboratory, Qingyuan People’s Hospital, The Sixth Affiliated Hospital of Guangzhou Medical University, Qingyuan, Guangdong, China
  • 6Department of Clinical Laboratory, Jieyang People’s Hospital, Jieyang, Guangdong, China

Purpose: Candidozyma auris (C. auris) is an emerging fungal pathogen that is resistant to multiple drugs and poses a serious threat to global health. This study aimed to investigate the genetic diversity and antifungal resistance profiles of C. auris in southern China.

Methods: A total of 108 clinical C. auris isolates were collected from eight hospitals in Guangdong Province between January 2023 and January 2024. All the isolates were identified via matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). Phylogenetic analysis was conducted on the basis of the sequencing results of the RPB1 and D1/D2 genes. The mechanisms of resistance to fluconazole and echinocandins were investigated through sequencing of the ERG11 and FKS1 genes. Additionally, whole-genome sequencing was performed on echinocandin-resistant and echinocandin-sensitive isolates to analyse genetic homology among the strains.

Results: All the strains were classified into two genetic clades, clade I and clade III, and all exhibited resistance to fluconazole. In the fluconazole-resistant strains, the amino acid substitutions Y132F and VF125AL were identified in the ERG11 gene. The resistance rates to caspofungin, micafungin, and anidulafungin were 7.4, 7.4, and 3.7%, respectively. Among the eight echinocandin-resistant strains, amino acid substitutions (S639Y, W691L, and S639F) were found within the HS1 hotspot region of the FKS1 gene. A phylogenetic tree was constructed on the basis of 403 SNPs and revealed two major clusters: Cluster A and Cluster B. Cluster A included the 16 isolates analyzed in this study. Cluster B consisted of 12 reference isolates retrieved from publicly available genomic databases.

Conclusion: Genetic clade I and clade III C. auris strains are prevalent in southern China and present high levels of resistance to fluconazole. Controlling the spread of C. auris in this region presents significant challenges for public health management.

1 Introduction

Candidozyma auris (formally Candida auris), an emerging multidrug-resistant fungal pathogen, was first isolated in Japan in 2009 (Satoh et al., 2009). The incidence of C. auris infections has risen dramatically around the world (Du et al., 2020; Akinbobola et al., 2023; Chakrabarti and Sood, 2021). It can cause serious infections, with an associated mortality rate ranging from 30 to 60% (Azar et al., 2017; Lockhart et al., 2017). Importantly, C. auris is resistant to most currently available antifungal medicines, and some strains are resistant to all three types of antifungal medicines (Chaabane et al., 2019). In the United States, approximately 90% of C. auris isolates are resistant to fluconazole, approximately 30% are resistant to amphotericin B, and less than 2% are resistant to echinocandins (Centers for Disease Control and Prevention, n.d.).

Additionally, C. auris can persist in hospital environments and on surfaces for extended periods, along with its high transmission ability. This poses significant challenges for clinical diagnosis, prevention, and control (Chakrabarti and Singh, 2020; Chowdhary et al., 2017). Owing to its intrinsic resistance to multiple antifungal classes and its high potential for transmission, C. auris is often referred to as a “superbug.” Consequently, the World Health Organization has included C. auris in its list of fungal priority pathogens (WHO, 2022).

In 2018, Wang et al. reported the first case of a nondrug-resistant C. auris strain at Peking University People’s Hospital in Beijing, China (Wang et al., 2018). Tian et al. identified 15 cases of C. auris colonization or infection at a general hospital in Shenyang, China (Tian et al., 2018). Recently, our team reported seven cases of C. auris bloodstream infections in two hospitals in Guangdong Province, China (Peng et al., 2024). To date, multiple healthcare facilities in southern China have reported instances of C. auris colonization and infection. To understand the epidemiological characteristics of C. auris in southern China, this study analyzed clinical isolates from eight hospitals, with a focus on geographic distribution, phylogenetic relationships, and hotspot mutations related to antifungal drug resistance. This study includes the seven previously reported bloodstream isolates.

2 Methods

2.1 Clinical data and isolates

Between January 2023 and January 2024, 108 C. auris isolates were collected from eight hospitals in southern China. Duplicated isolates obtained from the same patient were excluded.

All the isolates were identified as C. auris by matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) (bioMérieux, France) with the IVD knowledge base V3.2 according to the manufacturer’s instructions. Patient demographic and clinical data, such as sex, age, ward, underlying medical conditions, and antifungal treatments, were retrieved from the Hospital Information System.

2.2 Antifungal susceptibility testing

Antifungal susceptibility testing was performed using a commercial broth microdilution kit (Sensititre Yeastone™, Thermo Fisher Scientific). The antifungal agents tested included fluconazole (FLU), itraconazole (ITR), voriconazole (VOR), micafungin (MFG), caspofungin (CAS), anidulafungin (ANI), amphotericin B (AMB), 5-flucytosine (5-FC), and posaconazole (POS). Candida krusei (currently known as Pichia kudriavzevii) ATCC 6258 and Candida parapsilosis ATCC 22019 were used as quality control strains. The antifungal resistance of C. auris was interpreted according to the tentative minimum inhibitory concentration (MIC) breakpoints recommended by the US Centers for Disease Control and Prevention. Resistance was defined as an MIC ≥ 2 μg/mL for amphotericin B, ≥32 μg/mL for fluconazole, ≥2 μg/mL for caspofungin, ≥4 μg/mL for micafungin, and ≥2 μg/mL for anidulafungin.

2.3 Multilocus sequence typing

All strains were sequenced for the D1/D2 and RPB1 gene regions, and homology analysis was performed on the obtained sequences. The PCR and sequencing primers used are listed in Table 1. The procedure was carried out as follows: (1) The D1/D2 and RPB1 genes were amplified and sequenced. (2) Reference sequences of D1/D2 and RPB1 from geographically diverse strains were retrieved from the NCBI database. The gene accession numbers are provided in the supplementary material. (3) Sequences for both loci were concatenated for each strain. (4) Multiple sequence alignment and phylogenetic analysis were conducted using MEGA 7.0. A maximum likelihood phylogenetic tree was constructed with 1,000 bootstrap replicates.

Table 1
www.frontiersin.org

Table 1. List of primers used for PCR amplification and sequencing of target genes.

2.4 ERG11/FKS1 gene sequencing and analysis

The ERG11 gene was sequenced in fluconazole-resistant isolates. Additionally, the hotspot regions HS1 and HS2 of the FKS1 gene were analyzed in isolates showing reduced echinocandin susceptibility. Details of the primers used are provided in Table 1.

Sequence analysis involved assembling and aligning the obtained sequences with reference sequences from the NCBI database: FKS1 (accession numbers: OQ378936.1 and OQ378944.1) and ERG11 (accession numbers: MK294627 and MK294628). Reference sequences were derived from both fluconazole-susceptible and fluconazole-resistant strains, as well as from echinocandin-susceptible and echinocandin-resistant isolates. Potential nucleotide variations were identified using MEGA 7.0 software.

2.5 Whole-genome sequencing and phylogenetic analysis

Whole-genome sequencing was performed on 16 isolates using the TELL-seq library construction and sequencing methods. Among these strains, eight isolates exhibited echinocandin resistance, whereas the remaining eight were echinocandin susceptible. The 8 echinocandin-susceptible isolates were selected from hospitals where echinocandin-resistant strains had been isolated. During isolate selection, susceptible strains were randomly chosen according to the geographic distribution of the echinocandin-resistant strains.

TELL-seq libraries were constructed using the TELL-Seq Microbial Library Prep Kit (Universal Sequencing Technology), with 0.1 ng of genomic DNA per isolate. Library quantification was carried out using the Qubit 1X dsDNA HS Assay Kit (Thermo Fisher Scientific), followed by pooling and sequencing on an Illumina NovaSeq 6000 instrument, generating 2 × 150 paired-end reads along with 18-cycle Index 1 and 8-cycle Index 2 reads, according to the manufacturer’s protocols.

Single-nucleotide polymorphisms (SNPs) were identified using the Haplotype Caller module of the Genome Analysis Toolkit (GATK). The resulting gVCF files were merged via CombineGVCFs, and the SNPs were filtered on the basis of a variant quality score (<30) and quality by depth (<10). SNP calling was further refined using the Snippy v4.6.0 pipeline.

This study used the C. auris reference genome (GenBank accession number: GCA_014217455.1) for sequence alignment. The whole-genome sequencing data of twelve C. auris isolates (clade III) from public databases were included in the analysis. A maximum likelihood phylogenetic tree was constructed using SNP data with IQ-TREE software, with the TVMe + ASC + R3 model applied. Ultrafast bootstrap analysis was performed with 1,000 replicates.

3 Results

3.1 Clinical characteristics

The study included 108 patients with a median age of 64.3 years (range: 22–83 years). A total of 54.6% of the patients were admitted from the intensive care unit (ICU), and the male-to-female ratio was 2.2:1. The mean length of hospital stay was 53 days (range: 7–267 days). The median time from admission to the first identification of C. auris was 29.5 days (range: 4–129 days).

Prior to positive cultures for C. auris, 73.1% of the patients had undergone mechanical ventilation, 89.7% had been administered broad-spectrum antibiotics, and 91.0% had indwelling urinary catheters. Notably, 30.8% of patients were colonized or infected at multiple anatomical sites. At discharge, 34.7% of the patients had cleared the organism. Sputum and bronchoalveolar lavage fluid (BALF) are the most common sources of clinical samples. The detailed demographic and clinical characteristics are summarized in Table 2.

Table 2
www.frontiersin.org

Table 2. Demographic and clinical characteristics of the 108 patients.

3.2 Antifungal susceptibility testing

The antifungal susceptibility profiles of all the isolates are summarized in Table 3. All the isolates were resistant to fluconazole. The resistance rates to caspofungin, micafungin, and anidulafungin were 7.4, 7.4, and 3.7%, respectively. The resistance rate to amphotericin B was 16.7%. Among the triazoles, posaconazole presented the lowest geometric mean (GM) values. Caspofungin, micafungin, and anidulafungin displayed comparable MIC50 values, GM values, and MIC distributions. Compared with the Clade III strains (10/93, 10.8% resistant), the 15 Clade I strains exhibited a higher rate of amphotericin B resistance (8/15, 53.3% resistant) (Figure 1). Notably, one strain exhibited multidrug resistance and was concurrently resistant to amphotericin B, echinocandins, and fluconazole.

Table 3
www.frontiersin.org

Table 3. In vitro antifungal susceptibility profile of all the isolates.

Figure 1
Bar graph comparing the percentage of Clade I and Clade III at different AMB MICs in micrograms per milliliter. At 0.5, Clade I is 6.7% and Clade III is 10.8%. At 1, Clade I is 40% and Clade III is 78.4%. At 2, Clade I is 53.3% and Clade III is 10.8%. Clade I is represented in white, Clade III in gray.

Figure 1. Distribution of amphotericin B MICs among clades I and III Candidozyma auris strains.

3.3 Multilocus sequence typing

On the basis of the sequencing analysis of the RPB1 and D1/D2 regions, 108 strains were grouped into two distinct genetic clades (clades I and III). Specifically, 93 strains (86.1%) belonged to clade III, whereas 15 strains (13.9%) were classified into clade I. Notably, the coexistence of clade I and clade III strains was observed in two hospitals. The phylogenetic relationships among the 108 C. auris strains were inferred using sequences from the RPB1 and D1/D2 regions and are presented in Figure 2.

Figure 2
Circular phylogenetic tree diagram showing distinct clades of C. auris. This phylogenetic tree based on the RPB1 and D1/D2 sequences. Clade labels include Clades I to V, with individual strains marked by various shapes like triangles and circles. Scale bar indicates genetic distance.

Figure 2. Phylogenetic relationships among C. auris isolates from multiple healthcare facilities. This phylogenetic tree illustrates the genetic relatedness of C. auris isolates based on the RPB1 and D1/D2 sequences. Each node represents an individual isolate. Different shapes and colors denote different healthcare facilities. The dendrogram shows that the strains are divided into two geographical clades: clade I (South Asian clade) and clade III (South African clade) and with the majority of the strains belonging to the clade III (South African clade).

3.4 ERG11 and FKS1 mutation analysis

Given that all 108 isolates exhibited fluconazole resistance, we performed ERG11 gene sequencing on all strains. All strains harbored mutations in the ERG11 gene. The VF125AL amino acid substitution was present in all the clade III strains, whereas the Y132F substitution was exclusively detected in the clade I strains. Notably, no isolate exhibited both ERG11 mutations at the same time. Additionally, eight echinocandin-resistant isolates carried mutations in the FKS1 gene. Specific amino acid substitutions (S639Y, W691L, and S639F) were detected in the HS1 hotspot region of the FKS1 gene. No mutations were detected in the HS2 hotspot region. Among these, five isolates exhibited the S639Y mutation, two carried the S639F mutation, and one contained the W691L mutation. All eight echinocandin-resistant strains belonged to clade III, with seven originating from the same hospital. Six isolates were recovered from urine cultures, one from a blood culture, and one from a catheter tip culture. The clinical characteristics and FKS1 gene mutations of the eight echinocandin-resistant strains are summarized in Figure 3.

Figure 3
Table displaying data on isolates with annotations for the FKS1 Hotspot 1 mutations at sites 639 and 691. Includes isolate identifiers, genetic clades, ward locations, sites of isolation, MIC values for ANI, CAS, and MIF, and specific FKS1mutations. Some mutants show changes such as S639F, S639Y, and W691L, with relevant wards.

Figure 3. Clinical characteristics and FKS1 mutations in eight echinocandin-resistant C. auris strains. The eight echinocandin-resistant isolates were obtained from two hospitals. All eight isolates belonged to clade III, of which six were isolated from urine specimens. The detected mutation sites, including S639Y, W691L, and S639F, were all located in the FKS1 gene.

3.5 Phylogenetic analysis of 16 echinocandin-resistant and -susceptible strains

The draft genome assemblies of the 16 isolates ranged in size from 12.3 to 14.7 Mb, with an average of 23 ± 4 scaffolds. The GC content was 45.10% ± 0.03%, the scaffold N50 value was 1824.2 ± 369.1 kilobases (kb), and the coverage depth was 722 × ± 190 × on the basis of quality-trimmed reads. A phylogenetic tree was constructed on the basis of 403 SNPs and revealed two major clusters: Cluster A and Cluster B (Figure 4). Cluster A included the 16 echinocandin-resistant and -susceptible isolates analyzed in this study, all of which were collected from patients in Southern China between 2021 and 2023. Cluster B consisted of 12 reference isolates retrieved from publicly available genomic databases. These external isolates serve as a phylogenetic context, highlighting the genetic distinctness of the locally circulating strains in Cluster A.

Figure 4
Phylogenetic tree diagram showing evolutionary relationships among C. auris isolates. Strains are divided into Group A and Group B. Echinocandin-susceptible strains are highlighted in blue, and echinocandin-resistant strains are highlighted in yellow. Each branch is labeled with distance values, illustrating genetic similarities and differences among the strains.

Figure 4. Phylogenetic tree of C. auris isolates from this study and selected isolates reported in China. Cluster A comprised 16 isolates collected from patients in Southern China between 2021 and 2023. Cluster B included 12 reference isolates obtained from publicly available genomic databases.

4 Discussion

The clinical features of C. auris infections identified in this study are consistent with those reported in other studies (Du et al., 2022). First, invasive infections are relatively uncommon, with the majority of isolates obtained from nonsterile sites such as urine, sputum, or bronchoalveolar lavage fluid. Second, the majority of affected patients are elderly individuals who have been admitted to ICUs. These patients often undergo various medical interventions, including urinary catheterization, the administration of broad-spectrum antibiotics, and mechanical ventilation.

This study identified two genetic clades of C. auris (clade I and clade III) on the basis of sequencing analysis of RPB1 and D1/D2. This distribution is consistent with epidemiological reports from other regions in China (Bing et al., 2024). The first reported case of C. auris infection in Guangdong Province was associated with travel, as the patient had received medical care at a hospital in South Africa, where Clade III strains are predominant (Szekely et al., 2019). Furthermore, interhospital patient transfers have contributed to dissemination in multiple healthcare facilities. In southern China, Clade I and Clade III of C. auris are co-circulating, which is different from other regions in the country. Notably, we found both Clade I and Clade III strains present simultaneously in the same hospitals—this finding has not been reported in any published multicenter or single-center studies from mainland China. The coexistence of these two clades suggests a more complex and diverse epidemiological landscape for C. auris in southern China.

In our study, all isolates were found to be resistant to fluconazole. The resistance rates for caspofungin, micafungin, and anidulafungin were 7.4, 7.4, and 3.7%, respectively. The resistance rate to amphotericin B was 16.7%. One strain exhibited multidrug resistance, showing resistance to amphotericin B, echinocandins, and fluconazole; it was isolated from a urine sample and caused a urinary tract infection in the patient. According to MLST (multilocus sequence typing) results, this strain belongs to clade III and carries mutations in both the ERG11 and FKS1 genes. Ultimately, the patient discontinued treatment.

Consistent with the well-documented high level of fluconazole resistance (Du et al., 2022), all strains exhibited fluconazole resistance. Among the azole antifungal agents tested, posaconazole presented the lowest MIC, followed by itraconazole. In contrast, the voriconazole MICs were significantly greater than those observed for both posaconazole and itraconazole. This finding is consistent with previously published data (Gaitán et al., 2017; Larkin et al., 2017; Calvo et al., 2016; Morales-López et al., 2017; Chowdhary et al., 2018). For example, a study from Colombia reported that 23.5% of isolates were nonsusceptible to voriconazole (MIC ≥ 2 μg/mL), whereas all isolates identified in Spain and Venezuela presented voriconazole MICs of ≥2 μg/mL (Gaitán et al., 2017; Calvo et al., 2016; Morales-López et al., 2017). In our current study, 8.7% of the isolates were found to be nonsusceptible to voriconazole.

It should be noted that there are currently no C. auris-specific antifungal susceptibility breakpoints established. The breakpoints currently used are those defined by U. S. Centers for Disease Control and Prevention (CDC), based on breakpoints for closely related Candida species and expert opinion. At present, the correlation between microbiological breakpoints and clinical treatment outcomes remains unknown. Besides, for other azoles such as posaconazole, itraconazole, and voriconazole, no validated clinical breakpoints are currently available for C. auris. Consider using fluconazole susceptibility as a surrogate for second generation triazole susceptibility assessment. However, isolates that are resistant to fluconazole may respond to other triazoles occasionally. The decision to treat with another triazole will need to be made on case-by-case basis.

Fluconazole and other triazole antifungal agents exert their antifungal activity by competitively inhibiting sterol 14α-demethylase, a key enzyme in the ergosterol biosynthesis pathway. A common mechanism of fluconazole resistance involves mutations in the ERG11 gene, which encodes this essential enzyme. The VF125AL, Y132F, and K143R mutations have been frequently observed in fluconazole-resistant C. auris isolates (Rybak et al., 2022). Moreover, these mutations are associated with specific genetic clades (Lockhart et al., 2017). The VF125AL substitution was exclusively detected in clade III strains, whereas the Y132F substitution was exclusively detected in clade I and IV strains (Chow et al., 2020; Kwon et al., 2019). In our study, a consistent pattern was observed: all the clade III strains carried the VF125AL substitution, whereas all the clade I isolates harbored the Y132F mutation.

In contrast to the high amphotericin B resistance rate of 43% reported among C. auris strains in the United States, most isolates in China have exhibited low MICs for amphotericin B (Du et al., 2022). For example, Tian et al. reported a resistance rate of only 1.1% (1 out of 93) in Shenyang, China (Tian et al., 2021). However, the amphotericin B resistance observed in our study was greater. This discrepancy may be due to differences in the geographic clade distributions of the studied isolates. Specifically, all strains in Tian et al.’s study belonged to clade III, whereas our study included clade I and clade III strains. Emerging evidence suggests significant geographic variation in amphotericin B resistance rates among C. auris strains (Sanyaolu et al., 2022). Studies have shown that clade I has a 50% resistance rate to AMB—the highest rate in its clades. In contrast to the high amphotericin B resistance rate, clade III remains highly susceptible to AMB (da Silva et al., 2025). Previous studies suggest that putative loss-of-function mutations in ERG6, ERG3, and ERG11 can cause acquired AMB resistance in C. auris (Rybak et al., 2022; Carolus et al., 2021). However, as many AMB resistance isolates of C. auris exhibit lower levels of resistance (MIC 2–4 mg/L), and the vast majority are observed to lack mutations in ergosterol-biosynthesis genes, it appears likely that other mechanisms of AMB resistance remain to be identified (Sharma and Kadosh, 2023).

Fortunately, resistance to echinocandins remains relatively rare among C. auris isolates, suggesting that these drugs are suitable as first-line therapeutic agents. However, clinicians should be aware that echinocandin-resistant C. auris isolates have been reported. Of particular concern is the emergence of echinocandin-resistant strains in Guangdong Province within approximately 1 year of the initial detection of C. auris. The rapid emergence of echinocandin resistance in this region, occurring within approximately 1 year of the initial detection of C. auris, highlights the notable ability of this pathogen to acquire resistance. This phenomenon has been reported in studies, where strains that were initially susceptible developed echinocandin resistance (Rhodes et al., 2018; Centers for Disease Control and Prevention, 2017; Spruijtenburg et al., 2023). Bram Spruijtenburg et al. described a case in which resistance to micafungin developed after only 4 days of treatment, despite initial susceptibility (Spruijtenburg et al., 2023).

In this study, among the 8 patients with echinocandin-resistant isolates, the initial isolates from 2 patients were susceptible to echinocandins. Previous investigations into C. auris and other clinically relevant Candida species have demonstrated higher rates of echinocandin resistance among isolates recovered from urine than among those from other anatomical sources (Adams et al., 2018; Sharma et al., 2022). Consistent with these findings, the majority of the echinocandin-resistant isolates in our study were obtained from urine samples. This trend may be attributed to the low urinary concentrations of echinocandin drugs. For example, urinary excretion accounts for only 1.4% of the plasma concentration of caspofungin, whereas the rates of excretion from micafungin and anidulafungin are 0.7% and less than 0.1%, respectively (Denning, 2003). These subtherapeutic concentrations in urine are insufficient to effectively clear C. auris infection or colonization but may promote the selection and emergence of resistant strains (Spruijtenburg et al., 2023).

The primary mechanism of echinocandin resistance in Candida species involves the acquisition of mutations in the FKS genes, which encode the catalytic subunits of 1,3-β-D-glucan synthase, a key enzyme targeted by these antifungal agents. These mutations can lead to structural changes in the 1,3-β-D-glucan synthase complex, reducing its susceptibility to echinocandins. To date, three major amino acid substitutions at position S639 within the FKS1 gene have been identified in echinocandin-resistant C. auris isolates (Rybak et al., 2022; Denning, 2003; Kordalewska et al., 2018; Berkow and Lockhart, 2018; Rhodes et al., 2018). In line with these findings, our study confirms that the S639 (serine) amino acid position serves as the primary molecular determinant of echinocandin resistance in these strains.

Furthermore, whole-genome homology analysis revealed that the echinocandin- resistant strains in this study are phylogenetically more closely related to local susceptible isolates than to those from other geographic regions. The close phylogenetic clustering of echinocandin-resistant isolates within the same healthcare facility indicates a potential for nosocomial transmission or local development of resistance. However, the limited number of sequenced isolates and the absence of broader regional genomic data highlight the need for ongoing genomic surveillance.

This study represents the first comprehensive characterization of the epidemiological features of C. auris in southern China. However, the absence of whole-genome sequencing (WGS) data for all the isolates limits our ability to perform a detailed phylogenetic analysis of strain evolution. Consequently, the current findings lack sufficient molecular epidemiological evidence to delineate potential transmission routes between healthcare facilities. In addition, this study did not conduct in-depth analysis or discussion on mortality rates, clearance rates, and treatment regimens. Future studies should aim to elucidate the transmission dynamics of C. auris, investigate the mechanisms underlying resistance evolution, and assess their implications for clinical management and infection control practices.

Data availability statement

All the data generated or analyzed during this study are included in this article.

Ethics statement

This study was performed in accordance with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of the First Affiliated Hospital of Sun Yat-sen University (Ethics Number: [2022]041). Strains were collected as part of the routine clinical management of patients. Therefore, informed consent was not sought, and an informed consent waiver was approved by the institutional ethical committee of the First Affiliated Hospital of Sun Yat-sen University.

Author contributions

PG: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft. JF: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft. RW: Methodology, Writing – review & editing. WL: Resources, Writing – review & editing. YJ: Resources, Writing – review & editing. JL: Resources, Writing – review & editing. LC: Resources, Writing – review & editing. XH: Resources, Writing – review & editing. KL: Conceptualization, Methodology, Resources, Writing – review & editing. YP: Conceptualization, Data curation, Methodology, Resources, Writing – original draft.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors acknowledge Dr. Zhenxian Xiao and Dr. Yiwen Chen for their invaluable assistance and support in analyzing whole-genome data and constructing the phylogenetic tree.

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.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

References

Adams, E., Quinn, M., Tsay, S., Poirot, E., Chaturvedi, S., Southwick, K., et al. (2018). Candida auris in healthcare facilities, New York, USA, 2013-2017. Emerg. Infect. Dis. 24, 1816–1824. doi: 10.3201/eid2410.180649

PubMed Abstract | Crossref Full Text | Google Scholar

Akinbobola, A. B., Kean, R., Hanifi, S. M. A., and Quilliam, R. S. (2023). Environmental reservoirs of the drug-resistant pathogenic yeast Candida auris. PLoS Pathog. 19:e1011268. doi: 10.1371/journal.ppat.1011268

PubMed Abstract | Crossref Full Text | Google Scholar

Azar, M. M., Turbett, S. E., Fishman, J. A., and Pierce, V. M. (2017). Donor-derived transmission of Candida auris during lung transplantation. Clin. Infect. Dis. 65, 1040–1042. doi: 10.1093/cid/cix460

PubMed Abstract | Crossref Full Text | Google Scholar

Berkow, E. L., and Lockhart, S. R. (2018). Activity of CD101, a long-acting echinocandin, against clinical isolates of Candida auris. Diagn. Microbiol. Infect. Dis. 90, 196–197. doi: 10.1016/j.diagmicrobio.2017.10.021

PubMed Abstract | Crossref Full Text | Google Scholar

Bing, J., Du, H., Guo, P., Hu, T., Xiao, M., Lu, S., et al. (2024). Candida auris-associated hospitalizations and outbreaks, China, 2018-2023. Emerg. Microbes Infect. 13:2302843. doi: 10.1080/22221751.2024.2302843

PubMed Abstract | Crossref Full Text | Google Scholar

Calvo, B., Melo, A. S., Perozo-Mena, A., Hernandez, M., Francisco, E. C., Hagen, F., et al. (2016). First report of Candida auris in America: clinical and microbiological aspects of 18 episodes of candidemia. J. Infect. 73, 369–374. doi: 10.1016/j.jinf.2016.07.008

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–e03320. doi: 10.1128/mBio.03333-20

Crossref Full Text | Google Scholar

Centers for Disease Control and Prevention. Candida auris clinical update September (2017). Available online at: https://www.cdc.gov/fungal/diseases/candidiasis/c-aurisalert-09-17.html (Accessed March 26, 2019).

Google Scholar

Centers for Disease Control and Prevention. Candida auris. https://www.cdc.gov/fungal/diseases/candidiasis/candidaauris.html

Google Scholar

Chaabane, F., Graf, A., Jequier, L., and Coste, A. T. (2019). Review on antifungal resistance mechanisms in the emerging pathogen Candida auris. Front. Microbiol. 10:2788. doi: 10.3389/fmicb.2019.02788

PubMed Abstract | Crossref Full Text | Google Scholar

Chakrabarti, A., and Singh, S. (2020). Multidrug-resistant Candida auris: an epidemiological review. Expert Rev. Anti-Infect. Ther. 18, 551–562. doi: 10.1080/14787210.2020.1750368

PubMed Abstract | Crossref Full Text | Google Scholar

Chakrabarti, A., and Sood, P. (2021). On the emergence, spread and resistance of Candida auris: host, pathogen and environmental tipping points. J. Med. Microbiol. 70:001318. doi: 10.1099/jmm.0.001318

Crossref Full Text | Google Scholar

Chow, N. A., Munoz, 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

Crossref Full Text | Google Scholar

Chowdhary, A., Prakash, A., Sharma, C., Kordalewska, M., Kumar, A., Sarma, S., et al. (2018). A multicentre study of antifungal susceptibility patterns among 350 Candida auris isolates (2009-17) in India: role of the ERG11 and FKS1 genes in azole and echinocandin resistance. J. Antimicrob. Chemother. 73, 891–899. doi: 10.1093/jac/dkx480

PubMed Abstract | Crossref Full Text | Google Scholar

Chowdhary, A., Sharma, C., and Meis, J. F. (2017). Candida auris: A rapidly emerging cause of hospital-acquired multidrug-resistant fungal infections globally. PLoS Pathog. 13:e1006290. doi: 10.1371/journal.ppat.1006290

PubMed Abstract | Crossref Full Text | Google Scholar

da Silva, K. J. G., Lucini, F., Dos Santos, R. A. C., Santos, D. A., Meis, J. F., Melhem, M. S. C., et al. (2025). How does antifungal resistance vary in Candida (Candidozyma) auris and its clades? Quantitative and qualitative analyses and their clinical implications. Clin. Microbiol. Infect. 31, 1146–1156. doi: 10.1016/j.cmi.2025.04.003

PubMed Abstract | Crossref Full Text | Google Scholar

Denning, D. W. (2003). Echinocandin antifungal drugs. Lancet 362, 1142–1151. doi: 10.1016/S0140-6736(03)14472-8

PubMed Abstract | Crossref Full Text | Google Scholar

Du, H., Bing, J., Hu, T., Ennis, C. L., Nobile, C. J., and Huang, G. (2020). Candida auris: epidemiology, biology, antifungal resistance, and virulence. PLoS Pathog. 16:e1008921. doi: 10.1371/journal.ppat.1008921

PubMed Abstract | Crossref Full Text | Google Scholar

Du, H., Bing, J., Nobile, C. J., and Huang, G. (2022). Candida auris infections in China. Virulence 13, 589–591. doi: 10.1080/21505594.2022.2054120

PubMed Abstract | Crossref Full Text | Google Scholar

Gaitán, A. C. R., Moret, A., Hontangas, J. L. L., Molina, J. M., López, A. I. A., Cabezas, A. H., et al. (2017). Nosocomial fungemia by Candida auris: first four reported cases in continental Europe. Rev. Iberoam. Micol. 34, 23–27.

Google Scholar

Kordalewska, M., Lee, A., Park, S., Berrio, I., Chowdhary, A., Zhao, Y., et al. (2018). Understanding echinocandin resistance in the emerging pathogen Candida auris. Antimicrob. Agents Chemother. 62, e00238–e00218. doi: 10.1128/AAC.00238-18

Crossref Full Text | Google Scholar

Kwon, Y. J., Shin, J. H., Byun, S. A., Choi, M. J., Won, E. J., Lee, D., et al. (2019). Candida auris clinical isolates from South Korea: identification, antifungal susceptibility, and genotyping. J. Clin. Microbiol. 57, e01624–e01618. doi: 10.1128/JCM.01624-18

Crossref Full Text | Google Scholar

Larkin, E., Hager, C., Chandra, J., Mukherjee, P. K., Retuerto, M., Salem, I., et al. (2017). The emerging pathogen Candida auris: growth phenotype, virulence factors, activity of antifungals, and effect of SCY-078, a novel Glucan synthesis inhibitor, on growth morphology and biofilm formation. Antimicrob. Agents Chemother. 61, e02396–e02316. doi: 10.1128/AAC.02396-16

Crossref Full Text | Google Scholar

Lockhart, S. R., Etienne, K. A., Vallabhaneni, S., Farooqi, J., Chowdhary, A., Govender, N. P., et al. (2017). Simultaneous emergence of multidrug-resistant Candida auris on 3 continents confirmed by whole-genome sequencing and epidemiological analyses. Clin. Infect. Dis. 64, 134–140. doi: 10.1093/cid/ciw691

PubMed Abstract | Crossref Full Text | Google Scholar

Morales-López, S. E., Parra-Giraldo, C. M., Ceballos-Garzón, A., Martínez, H. P., Rodríguez, G. J., Álvarez-Moreno, C. A., et al. (2017). Invasive infections with multidrug-resistant yeast Candida auris, Colombia. Emerg. Infect. Dis. 23, 162–164. doi: 10.3201/eid2301.161497

Crossref Full Text | Google Scholar

Peng, Y., Liu, Y., and Guo, P. (2024). First report of Candida auris in Guangdong, China: clinical and microbiological characteristics of 7 episodes of candidemia. Emerg. Microbes Infect. 13:2300525. doi: 10.1080/22221751.2023.2300525

PubMed Abstract | Crossref Full Text | Google Scholar

Rhodes, J., Abdolrasouli, A., Farrer, R. A., Cuomo, C. A., Aanensen, D. M., Armstrong-James, D., et al. (2018). Genomic epidemiology of the UK outbreak of the emerging human fungal pathogen Candida auris. Emerg. Microbes Infect. 7:43. doi: 10.1038/s41426-018-0045-x

Crossref Full Text | Google Scholar

Rybak, J. M., Barker, K. S., Muñoz, J. F., Parker, J. E., Ahmad, S., Mokaddas, E., et al. (2022). In vivo emergence of high-level resistance during treatment reveals the first identified mechanism of amphotericin B resistance in Candida auris. Clin. Microbiol. Infect. 28, 838–843. doi: 10.1016/j.cmi.2021.11.024

PubMed Abstract | Crossref Full Text | Google Scholar

Rybak, J. M., Cuomo, C. A., and Rogers, P. D. (2022). The molecular and genetic basis of antifungal resistance in the emerging fungal pathogen Candida auris. Curr. Opin. Microbiol. 70:102208. doi: 10.1016/j.mib.2022.102208

PubMed Abstract | Crossref Full Text | Google Scholar

Sanyaolu, A., Okorie, C., Marinkovic, A., Abbasi, A. F., Prakash, S., Mangat, J., et al. (2022). Candida auris: an overview of the emerging drug-resistant fungal infection. Infect. Chemother. 54, 236–246. doi: 10.3947/ic.2022.0008

PubMed Abstract | Crossref Full Text | Google Scholar

Satoh, K., Makimura, K., Hasumi, Y., Nishiyama, Y., Uchida, K., and Yamaguchi, H. (2009). Candida auris sp. nov., a novel ascomycetous yeast isolated from the external ear canal of an inpatient in a Japanese hospital. Microbiol. Immunol. 53, 41–44. doi: 10.1111/j.1348-0421.2008.00083.x

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, C., and Kadosh, D. (2023). Perspective on the origin, resistance, and spread of the emerging human fungal pathogen Candida auris. PLoS Pathog. 19:e1011190. doi: 10.1371/journal.ppat.1011190

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, D., Paul, R. A., Rudramurthy, S. M., Kashyap, N., Bhattacharya, S., Soman, R., et al. (2022). Impact of FKS1 genotype on echinocandin in vitro susceptibility in Candida auris and in vivo response in a murine model of infection. Antimicrob. Agents Chemother. 66:e0165221. doi: 10.1128/AAC.01652-21

PubMed Abstract | Crossref Full Text | Google Scholar

Spruijtenburg, B., Ahmad, S., Asadzadeh, M., Alfouzan, W., Al-Obaid, I., et al. (2023). Whole genome sequencing analysis demonstrates therapy-induced echinocandin resistance in Candida auris isolates. Mycoses 66, 1079–1086. doi: 10.1111/myc.13655

PubMed Abstract | Crossref Full Text | Google Scholar

Szekely, A., Borman, A. M., and Johnson, E. M. (2019). Candida auris isolates of the southern Asian and south African lineages exhibit different phenotypic and antifungal susceptibility profiles in vitro. J. Clin. Microbiol. 57, e02055–e02018. doi: 10.1128/JCM.02055-18

Crossref Full Text | Google Scholar

Tian, S., Bing, J., Chu, Y., Chen, J., Cheng, S., Wang, Q., et al. (2021). Genomic epidemiology of Candida auris in a general hospital in Shenyang, China: a three-year surveillance study. Emerg. Microbes Infect. 10, 1088–1096. doi: 10.1080/22221751.2021.1934557

PubMed Abstract | Crossref Full Text | Google Scholar

Tian, S., Rong, C., Nian, H., Li, F., Chu, Y., Cheng, S., et al. (2018). First cases and risk factors of super yeast Candida auris infection or colonization from Shenyang, China. Emerg. Microbes Infect. 7:128. doi: 10.1038/s41426-018-0131-0

Crossref Full Text | Google Scholar

Wang, X. J., Bing, J., Zheng, Q. S., Zhang, F., Liu, J., Yue, H., et al. (2018). The first isolate of Candida auris in China: clinical and biological aspects. Emerg Microbes Infec. 7:93. doi: 10.1038/s41426-018-0095-0

Crossref Full Text | Google Scholar

WHO (2022). WHO fungal priority pathogens list to guide research, development and public health action. Available online at: https://www.who.int/publications/i/item/9789240060241

Google Scholar

Keywords: Candidozyma auris, molecular epidemiology, fungal drug resistance surveillance, antifungal resistance, drug resistance mechanisms

Citation: Guo P, Fang J, Wang R, Lin W, Jiang Y, Long J, Chen L, Huang X, Liao K and Peng Y (2025) Emergence of multidrug-resistant Candidozyma auris in southern China: a multicenter study on genetic diversity and antifungal resistance. Front. Microbiol. 16:1689176. doi: 10.3389/fmicb.2025.1689176

Received: 20 August 2025; Accepted: 29 September 2025;
Published: 10 October 2025.

Edited by:

Renátó Kovács, University of Debrecen, Hungary

Reviewed by:

Husam Eldin Salah Ismail Mohamed, Hamad Medical Corporation, Qatar
Hazrat Bilal, Jiangxi Cancer Hospital, China
Kássia Jéssica Galdino Silva, Université de Nantes, France

Copyright © 2025 Guo, Fang, Wang, Lin, Jiang, Long, Chen, Huang, Liao and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kang Liao, bGlhb2thbmcxOTcxQDE2My5jb20=; Yaqin Peng, cHlxZHJlYW1AMTYzLmNvbQ==

These authors have contributed equally to this work

Disclaimer: 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.