Genotypic diversity and unrecognized antifungal resistance among populations of Candida glabrata from positive blood cultures

Microbe Profile: Candida glabrata - a master of deception

Candida glabrata is a fungal microbe associated with multiple vertebrate microbiomes and their terrestrial environments. In humans, the species has emerged as an opportunistic pathogen that now ranks as the second-leading cause of candidiasis in Europe and North America (Beardsley et al. Med Mycol 2024, 62). People at highest risk of infection include the elderly, immunocompromised individuals and/or long-term residents of hospital and assisted-living facilities. C. glabrata is intrinsically drug-resistant, metabolically versatile and able to avoid detection by the immune system. Analyses of its 12.3 Mb genome indicate a stable pangenome Marcet-Houben et al. (BMC Biol 2022, 20) and phylogenetic affinity with Saccharomyces cerevisiae. Recent phylogenetic analyses suggest reclassifying C. glabrata as Nakaseomyces glabratus Lakashima and Sugita (Med Mycol J 2022, 63: 119-132).

Nosocomial transmission of fluconazole-resistant Candida glabrata bloodstream isolates revealed by whole-genome sequencing

The clonal transmission of fluconazole-resistant Candida glabrata isolates within hospitals has seldom been analyzed by whole-genome sequencing (WGS). We performed WGS on 79 C. glabrata isolates, comprising 31 isolates from three premature infants with persistent C. glabrata bloodstream infection despite antifungal treatment in the same neonatal intensive care unit (NICU) in 2022 and 48 (27 fluconazole-resistant and 21 fluconazole-susceptible dose-dependent) bloodstream isolates from 48 patients in 15 South Korean hospitals from 2010 to 2022. Phylogenetic analysis based on WGS single-nucleotide polymorphisms (SNPs) distinguished the 79 isolates according to multilocus sequence typing (MLST) (17 sequence type [ST]3, 13 ST7, two ST22, 41 ST26, four ST55, and two ST59 isolates) and unveiled two possible clusters of nosocomial transmission among ST26 isolates. One cluster from two premature infants with overlapping NICU hospitalizations in 2022 encompassed 15 fluconazole-resistant isolates harboring pleiotropic drug-resistance transcription factor (Pdr1p) P258L (13 isolates) or N1086I (two isolates), together with 10 fluconazole-susceptible dose-dependent isolates lacking Pdr1p SNPs. The other cluster indicated unforeseen clonal transmission of fluconazole-resistant bloodstream isolates among five patients (four post-lung transplantation and one with diffuse interstitial lung disease) in the same hospital over 8 months. Among these five isolates, four obtained after exposure to azole antifungals harbored distinct Pdr1p SNPs (N1091D, E388Q, K365E, and R376Q). The findings reveal the transmission patterns of clonal bloodstream isolates of C. glabrata among patients undergoing antifungal treatment, exhibiting different levels of fluconazole susceptibility or distinct Pdr1p SNP profiles.

IMPORTANCE

The prevalence of fluconazole-resistant bloodstream infections caused by Candida glabrata is increasing globally, but the transmission of these resistant strains within hospitals has rarely been documented. Through whole-genome sequencing and epidemiological analyses, this study identified two potential clusters of C. glabrata bloodstream infections within the same hospital, revealing the transmission of clonal C. glabrata strains with different levels of fluconazole susceptibility or distinct transcription factor pleiotropic drug resistance protein 1 (Pdr1p) single-nucleotide polymorphism profiles among patients receiving antifungal therapy.

Evaluation of the Eazyplex®Candida ID LAMP Assay for the Rapid Diagnosis of Positive Blood Cultures

Rapid molecular assays can be used to identify Candida pathogens directly from positive blood cultures (BCs) in a timely manner compared to standard methods using subcultures. In this study, the eazyplex®Candida ID assay, which is based on loop-mediated amplification (LAMP) and is currently for research use only, was evaluated for the identification of the most common fungal species. A total of 190 BCs were analysed. Sensitivity and specificity were 93.88% and 99.26% for C. albicans, 89.13% and 100% for Nakaseomyces glabratus (N. glabratus), 100% and 100% for Pichia kudravzevii (P. kudriavzevii), 100% and 100% for C. tropicalis, and 100% and 99.44% for C. parapsilosis. Sample preparation took approximately 11 min and positive amplification results were obtained between 8.5 and 19 min. The eazyplex®Candida ID LAMP assay is an easy-to-use diagnostic tool that can optimise the management of patients with candidemia.

Virulence-Related Genes Expression in Planktonic Mixed Cultures of Candida albicans and Non-Albicans Candida Species

INTRODUCTION

Candida albicans is the most common opportunistic pathogen causing fungal infections worldwide, especially in high-risk patients. Its pathogenicity is related to virulence factors gene expression, such as hyphal growth (HWP1), cell adhesion (ALS3), and protease secretion (SAP1) during infection spreading mechanisms. In recent years, an increase in non-albicans Candida infections has been reported, which may present coinfection or competitive interactions with C. albicans, potentially aggravating the patient's condition. This study aims to evaluate the expression of genes related to virulence factors of C. albicans and non-albicans Candida during planktonic stage.

METHODS

C. albicans (ATCC MYA-3573) as well as with three clinical strains (C. albicans DCA53, C. tropicalis DCT6, and C. parapsilosis DCP1) isolated from blood samples, were grown in 24-well plates at 37°C for 20 h, either in monocultures or mixed cultures. Quantitative real-time polymerase chain reaction was used to evaluate the expression levels of the genes HWP1, ALS3, and SAP1 in cells collected during the planktonic stage. In addition, hyphal filamentation was observed using a Scanning Electron Microscope.

RESULTS

The overexpression of HWP1 and ASL3 genes in mixed growth conditions between C. albicans and non-albicans Candida species suggests a synergistic relationship as well as an increased capacity for hyphal growth and adhesion. In contrast, C. parapsilosis versus C. tropicalis interaction shows an antagonistic relationship during mixed culture, suggesting a decreased virulence profile of C. parapsilosis during initial coinfection with C. tropicalis.

CONCLUSION

The expression of HWP1, ALS3, and SAP1 genes associated with virulence factors varies under competitive conditions among species of the genus Candida during planktonic stage.

Strain heterogeneity in a non-pathogenic Aspergillus fungus highlights factors associated with virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

Implication of genotypes for prognosis of Candida glabrata bloodstream infections

BACKGROUND

Genotyping isolates of a specific pathogen may demonstrate unique patterns of antimicrobial resistance, virulence or outcomes. However, evidence for genotype-outcome association in Candida glabrata is scarce. We aimed to characterize the mycological and clinical relevance of genotypes on C. glabrata bloodstream infections (BSIs).

METHODS

Non-duplicated C. glabrata blood isolates from hospitalized adults were genotyped by MLST, and further clustered by the unweighted pair group method with arithmetic averages (UPGMA). A clonal complex (CC) was defined by UPGMA similarities of >90%. Antifungal susceptibility testing was performed by a colorimetric microdilution method and interpreted following CLSI criteria.

RESULTS

Of 48 blood isolates evaluated, 13 STs were identified. CC7 was the leading CC (n = 14; 29.2%), including 13 ST7. The overall fluconazole and echinocandin resistance rates were 6.6% and 0%, respectively. No specific resistance patterns were associated with CC7 or other CCs. Charlson comorbidity index (adjusted OR, 1.49; 95% CI, 1.05-3.11) was the only predictor for CC7. By multivariable Cox regression analyses, CC7 was independently associated with 28 day mortality [adjusted HR (aHR), 3.28; 95% CI, 1.31-8.23], even after considering potential interaction with neutropenia (aHR, 3.41; 95% CI, 1.23-9.42; P for interaction, 0.24) or limited to 34 patients with monomicrobial BSIs (aHR, 2.85; 95% CI, 1.15-7.08). Also, the Kaplan-Meier estimate showed greater mortality with CC7 (P = 0.003). Fluconazole resistance or echinocandin therapy had no significant impact on mortality.

CONCLUSIONS

Our data suggested comorbid patients were at risk of developing CC7 BSIs. Further, CC7 was independently associated with worse outcomes.

Candida albicans and Candida glabrata: global priority pathogens

SUMMARYA significant increase in the incidence of Candida-mediated infections has been observed in the last decade, mainly due to rising numbers of susceptible individuals. Recently, the World Health Organization published its first fungal pathogen priority list, with Candida species listed in medium, high, and critical priority categories. This review is a synthesis of information and recent advances in our understanding of two of these species-Candida albicans and Candida glabrata. Of these, C. albicans is the most common cause of candidemia around the world and is categorized as a critical priority pathogen. C. glabrata is considered a high-priority pathogen and has become an increasingly important cause of candidemia in recent years. It is now the second most common causative agent of candidemia in many geographical regions. Despite their differences and phylogenetic divergence, they are successful as pathogens and commensals of humans. Both species can cause a broad variety of infections, ranging from superficial to potentially lethal systemic infections. While they share similarities in certain infection strategies, including tissue adhesion and invasion, they differ significantly in key aspects of their biology, interaction with immune cells, host damage strategies, and metabolic adaptations. Here we provide insights on key aspects of their biology, epidemiology, commensal and pathogenic lifestyles, interactions with the immune system, and antifungal resistance.

Microevolution of Candida glabrata (Nakaseomyces glabrata) during an infection

Candida glabrata (Nakaseomyces glabrata) is an emergent and opportunistic fungal pathogen that colonizes and persists in different niches within its human host. In this work, we studied five clinical isolates from one patient (P7), that have a clonal origin, and all of which come from blood cultures except one, P7-3, obtained from a urine culture. We found phenotypic variation such as sensitivity to high temperature, oxidative stress, susceptibility to two classes of antifungal agents, and cell wall porosity. Only isolate P7-3 is highly resistant to the echinocandin caspofungin while the other four isolates from P7 are sensitive. However, this same isolate P7-3, is the only one that displays susceptibility to fluconazole (FLC), while the rest of the isolates are resistant to this antifungal. We sequenced the PDR1 gene which encodes a transcription factor required to induce the expression of several genes involved in the resistance to FLC and found that all the isolates encode for the same Pdr1 amino acid sequence except for the last isolate P7-5, which contains a single amino acid change, G1099C in the putative Pdr1 transactivation domain. Consistent with the resistance to FLC, we found that the CDR1 gene, encoding the main drug efflux pump in C. glabrata, is highly overexpressed in the FLC-resistant isolates, but not in the FLC-sensitive P7-3. In addition, the resistance to FLC observed in these isolates is dependent on the PDR1 gene. Additionally, we found that all P7 isolates have a different proportion of cell wall carbohydrates compared to our standard strains CBS138 and BG14. In P7 isolates, mannan is the most abundant cell wall component, whereas β-glucan is the most abundant component in our standard strains. Consistently, all P7 isolates have a relatively low cell wall porosity compared to our standard strains. These data show phenotypic and genotypic variability between clonal isolates from different niches within a single host, suggesting microevolution of C. glabrata during an infection.

Strain heterogeneity in a non-pathogenic fungus highlights factors contributing to virulence

Fungal pathogens exhibit extensive strain heterogeneity, including variation in virulence. Whether closely related non-pathogenic species also exhibit strain heterogeneity remains unknown. Here, we comprehensively characterized the pathogenic potentials (i.e., the ability to cause morbidity and mortality) of 16 diverse strains of Aspergillus fischeri, a non-pathogenic close relative of the major pathogen Aspergillus fumigatus. In vitro immune response assays and in vivo virulence assays using a mouse model of pulmonary aspergillosis showed that A. fischeri strains varied widely in their pathogenic potential. Furthermore, pangenome analyses suggest that A. fischeri genomic and phenotypic diversity is even greater. Genomic, transcriptomic, and metabolomic profiling identified several pathways and secondary metabolites associated with variation in virulence. Notably, strain virulence was associated with the simultaneous presence of the secondary metabolites hexadehydroastechrome and gliotoxin. We submit that examining the pathogenic potentials of non-pathogenic close relatives is key for understanding the origins of fungal pathogenicity.

Functional genetic characterization of stress tolerance and biofilm formation in Nakaseomyces (Candida) glabrata via a novel CRISPR activation system

The overexpression of genes frequently arises in Nakaseomyces (formerly Candida) glabrata via gain-of-function mutations, gene duplication, or aneuploidies, with important consequences on pathogenesis traits and antifungal drug resistance. This highlights the need to develop specific genetic tools to mimic and study genetic amplification in this important fungal pathogen. Here, we report the development, validation, and applications of the first clustered regularly interspaced short palindromic repeats (CRISPR) activation (CRISPRa) system in N. glabrata for targeted genetic overexpression. Using this system, we demonstrate the ability of CRISPRa to drive high levels of gene expression in N. glabrata, and further assess optimal guide RNA targeting for robust overexpression. We demonstrate the applications of CRISPRa to overexpress genes involved in fungal pathogenesis and drug resistance and detect corresponding phenotypic alterations in these key traits, including the characterization of novel phenotypes. Finally, we capture strain variation using our CRISPRa system in two commonly used N. glabrata genetic backgrounds. Together, this tool will expand our capacity for functional genetic overexpression in this pathogen, with numerous possibilities for future applications.IMPORTANCENakaseomyces (formerly Candida) glabrata is an important fungal pathogen that is now the second leading cause of candidiasis infections. A common strategy that this pathogen employs to resist antifungal treatment is through the upregulation of gene expression, but we have limited tools available to study this phenomenon. Here, we develop, optimize, and apply the use of CRISPRa as a means to overexpress genes in N. glabrata. We demonstrate the utility of this system to overexpress key genes involved in antifungal susceptibility, stress tolerance, and biofilm growth. This tool will be an important contribution to our ability to study the biology of this important fungal pathogen.

Measuring individual colony MICs is a more sensitive method to detect the effect of antimicrobials on antimicrobial susceptibility than the proportion of colonies resistant

The ResistAZM randomized controlled trial found that the receipt of ceftriaxone/azithromycin, compared to ceftriaxone was not associated with an increase in the proportion of oral commensal Neisseria spp. and streptococci with azithromycin resistance 14 days after treatment. We repeated the analyses by measuring the minimum inhibitory concentrations (MICs) of azithromycin and ceftriaxone for individual colonies of commensal Neisseria spp. and streptococci at day 0 and day 14 in both arms. The receipt of ceftriaxone/azithromycin but not ceftriaxone was associated with an increase in azithromycin MIC for both Neisseria spp. (P < 0.0001) and streptococci (P = 0.0076). Likewise, ceftriaxone/azithromycin but not ceftriaxone monotherapy was associated with an increase in ceftriaxone MICs in Neisseria spp. (P = 0.0035). Whereas the proportion method failed to detect an association between the receipt of azithromycin and increased macrolide resistance, the MIC distribution method detected this effect. The MIC distribution method is thus a more sensitive method to assess the effect of antimicrobials on antimicrobial susceptibility.

BACKGROUND

The ResistAZM randomized controlled trial found that the receipt of ceftriaxone/azithromycin, compared to ceftriaxone was not associated with an increase in the proportion of oral commensal Neisseria spp. and streptococci with azithromycin resistance 14 days after treatment.

METHODS

We repeated the analyses by measuring the minimum inhibitory concentrations (MICs) of azithromycin and ceftriaxone for individual colonies of commensal Neisseria spp. and streptococci at day 0 and day 14 in both arms.

RESULTS

The receipt of ceftriaxone/azithromycin but not ceftriaxone was associated with an increase in azithromycin MIC for both Neisseria spp. (P < 0.0001) and streptococci (P = 0.0076). Likewise, ceftriaxone/azithromycin but not ceftriaxone monotherapy was associated with an increase in ceftriaxone MICs in Neisseria spp. (P = 0.0035).

CONCLUSIONS

Whereas the proportion method failed to detect an association between the receipt of azithromycin and increased macrolide resistance, the MIC distribution method detected this effect. The MIC distribution method is thus a more sensitive method to assess the effect of antimicrobials on antimicrobial susceptibility.

Evaluation of Antifungal Selective Toxicity Using Candida glabrata ERG25 and Human SC4MOL Knock-In Strains

With only four classes of antifungal drugs available for the treatment of invasive systemic fungal infections, the number of resistant fungi is increasing, highlighting the urgent need for novel antifungal drugs. Ergosterol, an essential component of cell membranes, and its synthetic pathway have been targeted for antifungal drug development. Sterol-C4-methyl monooxygenase (Erg25p), which is a greater essential target than that of existing drugs, represents a promising drug target. However, the development of antifungal drugs must consider potential side effects, emphasizing the importance of evaluating their selective toxicity against fungi. In this study, we knocked in ERG25 of Candida glabrata and its human ortholog, SC4MOL, in ERG25-deleted Saccharomyces cerevisiae. Utilizing these strains, we evaluated 1181-0519, an Erg25p inhibitor, that exhibited selective toxicity against the C. glabrata ERG25 knock-in strain. Furthermore, 1181-0519 demonstrated broad-spectrum antifungal activity against pathogenic Candida species, including Candida auris. The approach of utilizing a gene that is functionally conserved between yeast and humans and subsequently screening for molecular target drugs enables the identification of selective inhibitors for both species.

Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.