MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris

The laboratory investigation, management, and infection prevention and control of Candida auris: a narrative review to inform the 2024 national guidance update in England

The emergent fungal pathogen Candida auris is increasingly recognised as an important cause of healthcare-associated infections globally. It is highly transmissible, adaptable, and persistent, resulting in an organism with significant outbreak potential that risks devastating consequences. Progress in the ability to identify C. auris in clinical specimens is encouraging, but laboratory diagnostic capacity and surveillance systems are lacking in many countries. Intrinsic resistance to commonly used antifungals, combined with the ability to rapidly acquire resistance to therapy, substantially restricts treatment options and novel agents are desperately needed. Despite this, outbreaks can be interrupted, and mortality avoided or minimised, through the application of rigorous infection prevention and control measures with an increasing evidence base. This review provides an update on epidemiology, the impact of the COVID-19 pandemic, risk factors, identification and typing, resistance profiles, treatment, detection of colonisation, and infection prevention and control measures for C. auris. This review has informed a planned 2024 update to the United Kingdom Health Security Agency (UKHSA) guidance on the laboratory investigation, management, and infection prevention and control of Candida auris. A multidisciplinary response is needed to control C. auris transmission in a healthcare setting and should emphasise outbreak preparedness and response, rapid contact tracing and isolation or cohorting of patients and staff, strict hand hygiene and other infection prevention and control measures, dedicated or single-use equipment, appropriate disinfection, and effective communication concerning patient transfers and discharge.

Bimodal distribution of azole susceptibility in Sporothrix brasiliensis isolates in Brazil

Sporothrix brasiliensis is an emerging zoonotic fungal pathogen that can be difficult to treat. Antifungal susceptibility testing was performed on the mold phase of a convenience sample of 61 Sporothrix spp. isolates from human and cat sporotrichosis cases in Brazil using the Clinical and Laboratory Standards Institute standard M38. A bimodal distribution of azole susceptibility was observed with 50% (28/56) of S. brasiliensis isolates showing elevated itraconazole minimum inhibitory concentrations ≥16 µg/mL. Phylogenetic analysis found the in vitro resistant isolates were not clonal and were distributed across three different S. brasiliensis clades. Single nucleotide polymorphism (SNP) analysis was performed to identify potential mechanisms of in vitro resistance. Two of the 28 resistant isolates (MIC ≥16 mg/L) had a polymorphism in the cytochrome P450 gene, cyp51, corresponding to the well-known G448S substitution inducing azole resistance in Aspergillus fumigatus. SNPs corresponding to other known mechanisms of azole resistance were not identified in the remaining 26 in vitro resistant isolates.

Emergence of zoonotic sporotrichosis in Brazil: a genomic epidemiology study

BACKGROUND

Zoonotic sporotrichosis is a neglected fungal disease, whereby outbreaks are primarily driven by Sporothrix brasiliensis and linked to cat-to-human transmission. To understand the emergence and spread of sporotrichosis in Brazil, the epicentre of the current epidemic in South America, we aimed to conduct whole-genome sequencing (WGS) to describe the genomic epidemiology.

METHODS

In this genomic epidemiology study, we included Sporothrix spp isolates from sporotrichosis cases from Brazil, Colombia, and the USA. We conducted WGS using Illumina NovaSeq on isolates collected by three laboratories in Brazil from humans and cats with sporotrichosis between 2013 and 2022. All isolates that were confirmed to be Sporothrix genus by internal transcribed spacer or beta-tubulin PCR sequencing were included in this study. We downloaded eight Sporothrix genome sequences from the National Center for Biotechnology Information (six from Brazil, two from Colombia). Three Sporothrix spp genome sequences from the USA were generated by the US Centers for Disease Control and Prevention as part of this study. We did phylogenetic analyses and correlated geographical and temporal case distribution with genotypic features of Sporothrix spp isolates.

FINDINGS

72 Sporothrix spp isolates from 55 human and 17 animal sporotrichosis cases were included: 67 (93%) were from Brazil, two (3%) from Colombia, and three (4%) from the USA. Cases spanned from 1999 to 2022. Most (61 [85%]) isolates were S brasiliensis, and all were reported from Brazil. Ten (14%) were Sporothrix schenckii and were reported from Brazil, USA, and Colombia. For S schenckii isolates, two distinct clades were observed wherein isolates clustered by geography. For S brasiliensis isolates, five clades separated by more than 100 000 single-nucleotide polymorphisms were observed. Among the five S brasiliensis clades, clades A and C contained isolates from both human and cat cases, and clade A contained isolates from six different states in Brazil. Compared with S brasiliensis isolates, larger genetic diversity was observed among S schenckii isolates from animal and human cases within a clade.

INTERPRETATION

Our results suggest that the ongoing epidemic driven by S brasiliensis in Brazil represents several, independent emergence events followed by animal-to-animal and animal-to human transmission within and between Brazilian states. These results describe how S brasiliensis can emerge and spread within a country.

FUNDING

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil; the São Paulo Research Foundation; Productivity in Research fellowships by the National Council for Scientific and Technological Development, and Ministry of Science and Technology of Brazil.

Candida auris detected in the oral cavity of a dog in Kansas

Candida auris is an emerging human fungal pathogen, first described in Japan in 2009, and first detected in the United States in 2016. Here, we report the first-ever description of C. auris colonizing a human pet, the first identification of C. auris in a non-human mammal in the United States and the first C. auris isolate from the state of Kansas. While analyzing the oral mycobiome of dogs from a shelter in Kansas, the oral swab from one dog was found to contain C. auris as well as three other fungal species. The presence of C. auris in a dog suggests the possibility of zoonotic transmission to humans. The isolate is a member of Clade IV, which has been found in patients in Chicago and Florida, while Clades I and III are the most prevalent in the United States. The isolate is resistant to fluconazole, terbinafine, and amphotericin B but susceptible to caspofungin, consistent with the drug-resistant characteristics of many human C. auris isolates. The source of C. auris transient colonization in this dog is unknown, and there is no evidence that it was further transmitted to humans, other dogs in the shelter, or pets in its adopted household. Isolation of C. auris from a dog in Kansas has public health implications as a potential emerging source for the zoonotic spread of this pathogenic fungus, and for the development of antifungal resistance.IMPORTANCECandida auris is an emerging fungal infection of humans and is particularly problematic because it is multi-drug resistant and difficult to treat. It is also known to be spread from person to person by contact and can remain on surfaces for long periods of time. In this report, a dog in a shelter in Kansas is found to be colonized with Candida auris. This is the first study to document the presence of Candida auris on a pet, the first to document C. auris presence on a non-human mammal in the United States, and the first to report an isolate of C. auris within the state of Kansas. The presence of C. auris in a pet dog raises the possibility of zoonotic transmission from pets to human or vice versa.

Genomic description of acquired fluconazole- and echinocandin-resistance in patients with serial Candida glabrata isolates

Candida glabrata is one of the most common causes of systemic candidiasis, often resistant to antifungal medications. To describe the genomic context of emerging resistance, we conducted a retrospective analysis of 82 serially collected isolates from 33 patients from population-based candidemia surveillance in the United States. We used whole-genome sequencing to determine the genetic relationships between isolates obtained from the same patient. Phylogenetic analysis demonstrated that isolates from 29 patients were clustered by patient. The median SNPs between isolates from the same patient was 30 (range: 7-96 SNPs), while unrelated strains infected four patients. Twenty-one isolates were resistant to echinocandins, and 24 were resistant to fluconazole. All echinocandin-resistant isolates carried a mutation either in the FKS1 or FKS2 HS1 region. Of the 24 fluconazole-resistant isolates, 17 (71%) had non-synonymous polymorphisms in the PDR1 gene, which were absent in susceptible isolates. In 11 patients, a genetically related resistant isolate was collected after recovering susceptible isolates, indicating in vivo acquisition of resistance. These findings allowed us to estimate the intra-host diversity of C. glabrata and propose an upper boundary of 96 SNPs for defining genetically related isolates, which can be used to assess donor-to-host transmission, nosocomial transmission, or acquired resistance. IMPORTANCE In our study, mutations associated to azole resistance and echinocandin resistance were detected in Candida glabrata isolates using a whole-genome sequence. C. glabrata is the second most common cause of candidemia in the United States, which rapidly acquires resistance to antifungals, in vitro and in vivo.

First Canadian report of transmission of fluconazole-resistant Candida parapsilosis within two hospital networks confirmed by genomic analysis

Candida parapsilosis is a common cause of non-albicans candidemia. It can be transmitted in healthcare settings resulting in serious healthcare-associated infections and can develop drug resistance to commonly used antifungal agents. Following a significant increase in the percentage of fluconazole (FLU)-nonsusceptible isolates from sterile site specimens of patients in two Ontario acute care hospital networks, we used whole genome sequence (WGS) analysis to retrospectively investigate the genetic relatedness of isolates and to assess potential in-hospital spread. Phylogenomic analysis was conducted on all 19 FLU-resistant and seven susceptible-dose dependent (SDD) isolates from the two hospital networks, as well as 13 FLU susceptible C. parapsilosis isolates from the same facilities and 20 isolates from patients not related to the investigation. Twenty-five of 26 FLU-nonsusceptible isolates (resistant or SDD) and two susceptible isolates from the two hospital networks formed a phylogenomic cluster that was highly similar genetically and distinct from other isolates. The results suggest the presence of a persistent strain of FLU-nonsusceptible C. parapsilosis causing infections over a 5.5-year period. Results from WGS were largely comparable to microsatellite typing. Twenty-seven of 28 cluster isolates had a K143R substitution in lanosterol 14-α-demethylase (ERG11) associated with azole resistance. As the first report of a healthcare-associated outbreak of FLU-nonsusceptible C. parapsilosis in Canada, this study underscores the importance of monitoring local antimicrobial resistance trends and demonstrates the value of WGS analysis to detect and characterize clusters and outbreaks. Timely access to genomic epidemiological information can inform targeted infection control measures.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization

Introduction

The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens.

Methods

TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens.

Results

A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances.

Discussion

TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada

A Candida auris outbreak has been ongoing in Southern Nevada since August 2021. In this manuscript we describe the sequencing of over 200 C. auris isolates from patients at several facilities. Genetically distinct subgroups of C. auris were detected from Clade I (3 distinct lineages) and III (1 lineage). Open-source bioinformatic tools were developed and implemented to aid in the epidemiological investigation. The work herein compares three methods for C. auris whole genome analysis: Nullarbor, MycoSNP and a new pipeline TheiaEuk. We also describe a novel analysis method focused on elucidating phylogenetic linkages between isolates within an ongoing outbreak. Moreover, this study places the ongoing outbreaks in a global context utilizing existing sequences provided worldwide. Lastly, we describe how the generated results were communicated to the epidemiologists and infection control to generate public health interventions.

Comparing genomic variant identification protocols for Candida auris

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Candida auris Discovery through Community Wastewater Surveillance during Healthcare Outbreak, Nevada, USA, 2022

Candida auris transmission is steadily increasing across the United States. We report culture-based detection of C. auris in wastewater and the epidemiologic link between isolated strains and southern Nevada, USA, hospitals within the sampled sewershed. Our results illustrate the potential of wastewater surveillance for containing C. auris.

Candida auris in Austria-What Is New and What Is Different

Candida auris is a novel and emerging pathogenic yeast which represents a serious global health threat. Since its first description in Japan 2009, it has been associated with large hospital outbreaks all over the world and is often resistant to more than one antifungal drug class. To date, five C. auris isolates have been detected in Austria. Morphological characterization and antifungal susceptibility profiles against echinocandins, azoles, polyenes and pyrimidines, as well as the new antifungals ibrexafungerp and manogepix, were determined. In order to assess pathogenicity of these isolates, an infection model in Galleria mellonella was performed and whole genome sequencing (WGS) analysis was conducted to determine the phylogeographic origin. We could characterize four isolates as South Asian clade I and one isolate as African clade III. All of them had elevated minimal inhibitory concentrations to at least two different antifungal classes. The new antifungal manogepix showed high in vitro efficacy against all five C. auris isolates. One isolate, belonging to the African clade III, showed an aggregating phenotype, while the other isolates belonging to South Asian clade I were non-aggregating. In the Galleria mellonella infection model, the isolate belonging to African clade III exhibited the lowest in vivo pathogenicity. As the occurrence of C. auris increases globally, it is important to raise awareness to prevent transmission and hospital outbreaks.

Genomics and metagenomics of Madurella mycetomatis, a causative agent of black grain mycetoma in Sudan

Madurella mycetomatis is one of the main causative agents of mycetoma, a debilitating neglected tropical disease. Improved understanding of the genomic diversity of the fungal and bacterial causes of mycetoma is essential to advances in diagnosis and treatment. Here, we describe a high-quality genome assembly of M. mycetomatis and results of the whole genome sequence analysis of 26 isolates from Sudan. We demonstrate evidence of at least seven genetically diverse lineages and extreme clonality among isolates within these lineages. We also performed shotgun metagenomic analysis of DNA extracted from mycetoma grains and showed that M. mycetomatis reads were detected in all sequenced samples with the average of 11,317 reads (s.d. +/- 21,269) per sample. In addition, 10 (12%) of the 81 tested grain samples contained bacterial reads including Streptococcus sp., Staphylococcus sp. and others.

Generating Complete Genome Assemblies of Candida auris

Genomic studies of Candida auris are underpinned by the generation of high-quality genome assemblies. These reference genomes have been essential for investigations of the evolution and epidemiology of this emerging fungal pathogen. In addition to genomic epidemiology studies of local outbreaks and analysis of the global emergence of this species, comparisons of genomes of isolates from the five major clades have revealed differences in gene content and genomic structure. Here, we provide a detailed protocol for generating complete genome assemblies for C. auris.