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.

Genomic Epidemiology Linking Nonendemic Coccidioidomycosis to Travel.

Coccidioidomycosis is a fungal infection endemic to hot, arid regions of the western United States, northern Mexico, and parts of Central and South America. Sporadic cases outside these regions are likely travel-associated; alternatively, an infection could be acquired in as-yet unidentified newly endemic locales. A previous study of cases in nonendemic regions with patient self-reported travel history suggested that infections were acquired during travel to endemic regions. We sequenced 19 Coccidioides isolates from patients with known travel histories from that earlier investigation and performed phylogenetic analysis to identify the locations of potential source populations. Our results show that those isolates were phylogenetically linked to Coccidioides subpopulations naturally occurring in 1 of the reported travel locales, confirming that these cases were likely acquired during travel to endemic regions. Our findings demonstrate that genomic analysis is a useful tool for investigating travel-related coccidioidomycosis.

One health genomic surveillance and response to a university-based outbreak of the SARS-CoV-2 Delta AY.25 lineage, Arizona, 2021.

Genomic surveillance and wastewater tracking strategies were used to strengthen the public health response to an outbreak of the SARS-CoV-2 Delta AY.25 lineage associated with a university campus in Arizona. Epidemiologic and clinical data routinely gathered through contact tracing were matched to SARS-CoV-2 genomes belonging to an outbreak of AY.25 identified through ongoing phylogenomic analyses. Continued phylogenetic analyses were conducted to further describe the AY.25 outbreak. Wastewater collected twice weekly from sites across campus was tested for SARS-CoV-2 by RT-qPCR, and subsequently sequenced to identify variants. The AY.25 outbreak was defined by a single mutation (C18804T) and comprised 379 genomes from SARS-CoV-2 positive cases associated with the university and community. Several undergraduate student gatherings and congregate living settings on campus likely contributed to the rapid spread of COVID-19 across the university with secondary transmission into the community. The clade defining mutation was also found in wastewater samples collected from around student dormitories a week before the semester began, and 9 days before cases were identified. Genomic, epidemiologic, and wastewater surveillance provided evidence that an AY.25 clone was likely imported into the university setting just prior to the onset of the Fall 2021 semester, rapidly spread through a subset of the student population, and then subsequent spillover occurred in the surrounding community. The university and local public health department worked closely together to facilitate timely reporting of cases, identification of close contacts, and other necessary response and mitigation strategies. The emergence of new SARS-CoV-2 variants and potential threat of other infectious disease outbreaks on university campuses presents an opportunity for future comprehensive One Health genomic data driven, targeted interventions.

Phylogenomic Placement of American Southwest-Associated Clinical and Veterinary Isolates Expands Evidence for Distinct Cryptococcus gattii VGVI.

Whole-genome sequencing has advanced our understanding of the population structure of the pathogenic species complex Cryptococcus gattii, which has allowed for the phylogenomic specification of previously described major molecular type groupings and novel lineages. Recently, isolates collected in Mexico in the 1960s were determined to be genetically distant from other known molecular types and were classified as VGVI. We sequenced four clinical isolates and one veterinary isolate collected in the southwestern United States and Argentina from 2012 to 2021. Phylogenomic analysis groups these genomes with those of the Mexican VGVI isolates, expanding VGVI into a clade and establishing this molecular type as a clinically important population. These findings also potentially expand the known Cryptococcus ecological range with a previously unrecognized endemic area.

Molecular type distribution and fluconazole susceptibility of clinical Cryptococcus gattii isolates from South African laboratory-based surveillance, 2005-2013.

As is the case globally, Cryptococcus gattii is a less frequent cause of cryptococcosis than Cryptococcus neoformans in South Africa. We performed multilocus sequence typing (MLST) and fluconazole susceptibility testing of 146 isolates randomly selected from 750 South African patients with C. gattii disease identified through enhanced laboratory surveillance, 2005 to 2013. The dominant molecular type was VGIV (101/146, 70%), followed by VGI (40/146, 27%), VGII (3/146, 2%) and VGIII (2/146, 1%). Among the 146 C. gattii isolates, 99 different sequence types (STs) were identified, with ST294 (14/146, 10%) and ST155 (10/146, 7%) being most commonly observed. The fluconazole MIC50 and MIC90 values of 105 (of 146) randomly selected C. gattii isolates were 4 µg/ml and 16 µg/ml, respectively. VGIV isolates had a lower MIC50 value compared to non-VGIV isolates, but these values were within one double-dilution of each other. HIV-seropositive patients had a ten-fold increased adjusted odds of a VGIV infection compared to HIV-seronegative patients, though with small numbers (99/136; 73% vs. 2/10; 20%), the confidence interval (CI) was wide (95% CI: 1.93-55.31, p = 0.006). Whole genome phylogeny of 98 isolates of South Africa's most prevalent molecular type, VGIV, identified that this molecular type is highly diverse, with two interesting clusters of ten and six closely related isolates being identified, respectively. One of these clusters consisted only of patients from the Mpumalanga Province in South Africa, suggesting a similar environmental source. This study contributed new insights into the global population structure of this important human pathogen.

Genotype, Antifungal Susceptibility, and Virulence of Clinical South African Cryptococcus neoformans Strains from National Surveillance, 2005-2009.

In South Africa, Cryptococcus neoformans is the most common cause of adult meningitis. We performed multi locus sequence typing and fluconazole susceptibility testing of clinical C. neoformans isolates collected from 251 South African patients with cryptococcosis through national surveillance from 2005 to 2009. We examined the association between clinical characteristics of patients and genotype, and the effect of genotype on in-hospital mortality. We performed whole genome phylogenetic analysis of fifteen C. neoformans isolates with the molecular type VNB and tested their virulence in a Galleria mellonella model. Most isolates had the molecular type VNI (206/251, 82%), followed by VNII (25/251, 10%), VNB (15/251, 6%), and VNIV (5/251, 2%); 67 sequence types were identified. There were no differences in fluconazole minimum inhibitory concentration (MIC) values among molecular types and the majority of strains had low MIC values (MIC50 of 1 µg/mL and MIC90 of 4 µg/mL). Males were almost twice as likely of being infected with a non-VNI genotype (adjusted odds ratio [OR]: 1.65, 95% confidence interval [CI]: 0.25-10.99; p = 0.61). Compared to patients infected with a VNI genotype, those with a non-VNI genotype had a 50% reduced adjusted odds of dying in hospital (95% CI: 0.03-7.57; p = 0.62). However, for both these analyses, our estimates had wide confidence intervals spanning 1 with large p-values. Fifteen VNB strains were not as virulent in a G. mellonella larval model as the H99 reference strain. A majority of these VNB strains belonged to the VNBII clade and were very closely related by phylogenetic analysis.

Rhizopus microsporus Infections Associated with Surgical Procedures, Argentina, 2006-2014.

Rhizopus spp. fungi are ubiquitous in the environment and a rare but substantial cause of infection in immunosuppressed persons and surgery patients. During 2005-2017, an abnormally high number of Rhizopus infections in surgery patients, with no apparent epidemiologic links, were reported in Argentina. To determine the likelihood of a common source of the cluster, we performed whole-genome sequencing on samples collected during 2006-2014. Most isolates were separated by >60 single-nucleotide polymorphisms, and we found no evidence for recombination or nonneutral mutation accumulation; these findings do not support common source or patient-to-patient transmission. Assembled genomes of most isolates were ≈25 Mbp, and multiple isolates had substantially larger assembled genomes (43-51 Mbp), indicative of infections with strain types that underwent genome expansion. Whole-genome sequencing has become an essential tool for studying epidemiology of fungal infections. Less discriminatory techniques may miss true relationships, possibly resulting in inappropriate attribution of point source.

A unique multidrug-resistant clonal Trichophyton population distinct from Trichophyton mentagrophytes/Trichophyton interdigitale complex causing an ongoing alarming dermatophytosis outbreak in India: Genomic insights and resistance profile.

There has been a considerable upsurge of extensive, treatment recalcitrant, dermatophytosis presenting as tinea corporis and tinea cruris in India since the past few years. Genome analysis of Trichophytonspecies causing severe superficial dermatophytosis in North India confirmed a unique clade related to the T.mentagrophytes/interdigitale complex, seeming to belong to an early diverging clade of the complex. The Indian Trichophyton species genomes were highly related showing only up to 42 SNPs between any two isolates confirming their clonal origin. Other genetic approaches such as ITS sequencing and multigene phylogeny used in this study were contradictory or inconclusive to show the differentiation of these isolates from T. mentagrophytes/T. interdigitale. Remarkably, high rates of resistance to all three commonly used oral antifungals, i.e., 36% for terbinafine (MICs 4 to ≥32 mg/L), 39.5% for fluconazole (MIC range 32 to ≥64 mg/L) and griseofulvin (Geometric mean MIC ≥ 4 mg/L) were observed. Two important amino acid substitutions (Leu393Phe or Phe397Leu) leading to a terbinafine resistant phenotype were found in the squalene epoxidase protein of all tested terbinafine resistant isolates. All 20 examined genomes presented a high mobility group (HMG) domain transcription factor gene corresponding to mating type (+). Of these, three isolates also showed positivity for both alpha-box and HMG in the genome which might indicate hybridization or an incomplete sexual cycle. Therefore, we highlight the potential of this organism to rapidly spread alleles that might be driving antifungal resistance among its population. This new population of Trichophyton with high rates of in vitro antifungal resistance seems to be driving an ongoing outbreak of dermatophytosis in India. Our study highlights difficulties in identifying isolates from the Trichophyton mentagrophytes/interdigitale clade of the genus using currently available molecular tools. High resistance rates of terbinafine warrant further clinical studies to assess its utility in the treatment of dermatophytosis caused by this clonal strain.

Jet-Setting Koalas Spread Cryptococcus gattii VGII in Australia.

Cryptococcus gattii molecular type VGII is one of the etiologic agents of cryptococcosis, a systemic mycosis affecting a wide range of host species. Koalas (Phascolarctos cinereus) exhibit a comparatively high prevalence of cryptococcosis (clinical and subclinical) and nasal colonization, particularly in captivity. In Australia, disease associated with C. gattii VGII is typically confined to Western Australia and the Northern Territory (with sporadic cases reported in eastern Australia), occupying an enigmatic ecologic niche. A cluster of cryptococcosis in captive koalas in eastern Australia (five confirmed cases, a further two suspected), caused predominantly by C. gattii VGII, was investigated by surveying for subclinical disease, culturing koala nasal swabs and environmental samples, and genotyping cryptococcal isolates. URA5 restriction fragment length polymorphism analysis, multilocus sequence typing (MLST), and whole-genome sequencing (WGS) provided supportive evidence that the transfer of koalas from Western Australia and subsequently between several facilities in Queensland spread VGII into uncontaminated environments and environments in which C. gattii VGI was endemic. MLST identified VGII isolates as predominantly sequence type 7, while WGS further confirmed a limited genomic diversity and revealed a basal relationship with isolates from Western Australia. We hypothesize that this represents a founder effect following the introduction of a koala from Western Australia. Our findings suggest a possible competitive advantage for C. gattii VGII over VGI in the context of this captive koala environment. The ability of koalas to seed C. gattii VGII into new environments has implications for the management of captive populations and movements of koalas between zoos.IMPORTANCECryptococcus gattii molecular type VGII is one of the causes of cryptococcosis, a severe fungal disease that is acquired from the environment and affects many host species (including humans and koalas). In Australia, disease caused by C. gattii VGII is largely confined to western and central northern parts of the country, with sporadic cases reported in eastern Australia. We investigated an unusual case cluster of cryptococcosis, caused predominantly by C. gattii VGII, in a group of captive koalas in eastern Australia. This research identified that the movements of koalas between wildlife parks, including an initial transfer of a koala from Western Australia, introduced and subsequently spread C. gattii VGII in this captive environment. The spread of this pathogen by koalas could also impact other species, and these findings are significant in the implications they have for the management of koala transfers and captive environments.