Molecular epidemiology and antimicrobial resistance of vaginal Candida glabrata isolates in Namibia

Abstract Candida glabrata is the most common non-albicans Candida species that causes vulvovaginal candidiasis (VVC). Given the intrinsically low susceptibility of C. glabrata to azole drugs, investigations into C. glabrata prevalence, fungal susceptibility profile, and molecular epidemiology are necessary to optimise the treatment of VVC. This molecular epidemiological study was conducted to determine antifungal drug profile, single nucleotide polymorphisms (SNPs) associated with phenotypic antifungal resistance and epidemic diversity of C. glabrata isolates from women with VVC in Namibia. Candida glabrata isolates were identified using phenotypic and molecular methods. Antifungal susceptibility of strains was determined for fluconazole, itraconazole, amphotericin B, and anidulafungin. Whole genome sequencing was used to determine SNPs in antifungal resistance genes and sequence type (ST) allocation. Among C. glabrata isolates, all (20/20; 100%) exhibited phenotypic resistance to the azole class antifungal drug, (fluconazole), and phenotypic susceptibility to the polyene class (amphotericin B), and the echinocandins (anidulafungin). Non-synonymous SNPs were identified in antifungal resistance genes of all fluconazole-resistant C. glabrata isolates including ERG6 (15%), ERG7 (15%), CgCDR1 (25%), CgPDR1 (60%), SNQ2 (10%), FKS1 (5.0%), FKS2 (5.0%), CgFPS1 (5.0%), and MSH2 (15%). ST15 (n = 8/20, 40%) was predominant. This study provides important insight into phenotypic and genotypic antifungal resistance across C. glabrata isolates from women with VVC in Namibia. In this study, azole resistance is determined by an extensive range of SNPs, while the observed polyene and echinocandin resistance-associated SNPs despite phenotypic susceptibility require further investigation.