Large Clusters of Invasive emm49 Group A Streptococcus Identified Within Arizona Health Care Facilities Through Statewide Genomic Surveillance System, 2019-2021.

A statewide genomic surveillance system for invasive Group A Streptococcus was implemented in Arizona in June 2019, resulting in 1046 isolates being submitted for genomic analysis to characterize emm types and identify transmission clusters. Eleven of the 32 identified distinct emm types comprised >80% of samples, with 29.7% of all isolates being typed as emm49 (and its genetic derivative emm151). Phylogenetic analysis initially identified an emm49 genomic cluster of 4 isolates that rapidly expanded over subsequent months (June 2019 to February 2020). Public health investigations identified epidemiologic links with 3 different long-term care facilities, resulting in specific interventions. Unbiased genomic surveillance allowed for identification and response to clusters that would have otherwise remained undetected.

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.

Suspected Locally Acquired Coccidioidomycosis in Human, Spokane, Washington, USA.

The full geographic range of coccidioidomycosis is unknown, although it is most likely expanding with environmental change. We report an apparently autochthonous coccidioidomycosis patient from Spokane, Washington, USA, a location to which Coccidioides spp. are not known to be endemic.

Coyotes as Reservoirs for Onchocerca lupi, United States, 2015-2018.

The Onchocerca lupi nematode infects dogs, cats, and humans, but whether it can be spread by coyotes has been unknown. We conducted surveillance for O. lupi nematode infection in coyotes in the southwestern United States. We identified multiple coyote populations in Arizona and New Mexico as probable reservoirs for this species.

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.

The detection of Coccidioides from ambient air in Phoenix, Arizona: Evidence of uneven distribution and seasonality.

Coccidioidomycosis is a debilitating fungal disease caused by inhalation of arthroconidia. We developed a novel approach for detection of airborne Coccidioides and used it to investigate the distribution of arthroconidia across the Phoenix, Arizona, metropolitan area. Air filters were collected daily from 21 stationary air-sampling units across the area: the first set collected before, during and after a large dust storm on August 25, 2015, and the second over the 45-day period September 25-November 8, 2016. Analysis of DNA extracted from the filters demonstrated that the day of the dust storm was not associated with increase of Coccidioides in air samples, although evidence of the low-level polymerase chain reaction (PCR) inhibition was observed in DNA extracted from samples collected on the day of the dust storm. Testing over 45 days identified uneven geographic distribution suggesting Coccidioides hot spots. In 2016, highest daily concentration of arthroconidia was observed between September 25-October 20, and only sporadic low levels were detected after that. These results provide evidence of seasonality and uneven spatial distribution of Coccidioides in the air. Our results demonstrate that routine air monitoring for arthroconidia is possible and provides an important tool for Coccidioides surveillance, which can address important questions about environmental exposure and human infection.

KlebSeq, a Diagnostic Tool for Surveillance, Detection, and Monitoring of Klebsiella pneumoniae.

Health care-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to health care costs. Multidrug-resistant and epidemic strains are a large proportion of HAI agents, and multidrug-resistant strains of Klebsiella pneumoniae, a leading HAI agent, have caused an urgent public health crisis. In the health care environment, patient colonization by K. pneumoniae precedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening for K. pneumoniae colonization has the potential to curb the number of HAIs. In this report, we describe the design and validation of KlebSeq, a highly informative screening tool that detects Klebsiella species and identifies clinically important strains and characteristics by using highly multiplexed amplicon sequencing without a live-culturing step. We demonstrate the utility of this tool on several complex specimen types, including urine, wound swabs and tissue, and several types of respiratory and fecal specimens, showing K. pneumoniae species and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool to screen patients for Klebsiella carriage could inform health care staff of the risk of infection and outbreak potential. KlebSeq also serves as a model for next-generation molecular tools for public health and health care, as expansion of this tool can be used for several other HAI agents or applications.

Understanding the exposure risk of aerosolized Coccidioides in a Valley fever endemic metropolis.

Coccidioides is the fungal causative agent of Valley fever, a primarily pulmonary disease caused by inhalation of fungal arthroconidia, or spores. Although Coccidioides has been an established pathogen for 120 years and is responsible for hundreds of thousands of infections per year, little is known about when and where infectious Coccidioides arthroconidia are present within the ambient air in endemic regions. Long-term air sampling programs provide a means to investigate these characteristics across space and time. Here we present data from > 18 months of collections from 11 air sampling sites across the Phoenix, Arizona, metropolitan area. Overall, prevalence was highly variable across space and time with no obvious spatial or temporal correlations. Several high prevalence periods were identified at select sites, with no obvious spatial or temporal associations. Comparing these data with weather and environmental factor data, wind gusts and temperature were positively associated with Coccidioides detection, while soil moisture was negatively associated with Coccidioides detection. These results provide critical insights into the frequency and distribution of airborne arthroconidia and the associated risk of inhalation and potential disease that is present across space and time in a highly endemic locale.

Comparison of carbapenem-susceptible and carbapenem-resistant Enterobacterales at nine sites in the USA, 2013-2016: a resource for antimicrobial resistance investigators.

Carbapenem-resistant Enterobacterales (CRE) are an urgent public health threat. Genomic sequencing is an important tool for investigating CRE. Through the Division of Healthcare Quality Promotion Sentinel Surveillance system, we collected CRE and carbapenem-susceptible Enterobacterales (CSE) from nine clinical laboratories in the USA from 2013 to 2016 and analysed both phenotypic and genomic sequencing data for 680 isolates. We describe the molecular epidemiology and antimicrobial susceptibility testing (AST) data of this collection of isolates. We also performed a phenotype-genotype correlation for the carbapenems and evaluated the presence of virulence genes in Klebsiella pneumoniae complex isolates. These AST and genomic sequencing data can be used to compare and contrast CRE and CSE at these sites and serve as a resource for the antimicrobial resistance research community.

Unique Genomic Epidemiology of COVID-19 in the White Mountain Apache Tribe, April to August 2020, Arizona.

The first case of coronavirus disease 2019 (COVID-19) within the White Mountain Apache Tribe (WMAT) in Arizona was diagnosed almost 1 month after community transmission was recognized in the state. Aggressive contact tracing allowed for robust genomic epidemiology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and subsequent phylogenetic analyses implicated only two virus introductions, which resulted in the spread of two unique viral lineages on the reservation. The phylogenies of these lineages reflect the nature of the introductions, the remoteness of the community, and the extraordinarily high attack rates. The timing and space-limited nature of the outbreaks validate the public health tracing efforts involved, which were illustrated by multiple short transmission chains over a period of several weeks, eventually resulting in extinction of the lineages. Comprehensive sampling and successful infection control efforts are illustrated in both the effective population size analyses and the limited mortality outcomes. The rapid spread and high attack rates of the two lineages may be due to a combination of sociological determinants of the WMAT and a seemingly enhanced transmissibility. The SARS-CoV-2 genomic epidemiology of the WMAT demonstrates a unique local history of the pandemic and highlights the extraordinary and successful efforts of their public health response. IMPORTANCE This article discusses the introduction and spread of two unique viral lineages of SARS-CoV-2 within the White Mountain Apache Tribe in Arizona. Both genomic sequencing and traditional epidemiological strategies (e.g., contract tracing) were used to understand the nature of the spread of both lineages. Beyond providing a robust genomic analysis of the epidemiology of the outbreaks, this work also highlights the successful efforts of the local public health response.

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.

Comparison of TaqMan PCR assays for detection of the melioidosis agent Burkholderia pseudomallei in clinical specimens.

Melioidosis is an emerging infectious disease caused by the soil bacterium Burkholderia pseudomallei. In diagnostic and forensic settings, molecular detection assays need not only high sensitivity with low limits of detection but also high specificity. In a direct comparison of published and newly developed TaqMan PCR assays, we found the TTS1-orf2 assay to be superior in detecting B. pseudomallei directly from clinical specimens. The YLF/BTFC multiplex assay (targeting the Yersinia-like fimbrial/Burkholderia thailandensis-like flagellum and chemotaxis region) also showed high diagnostic sensitivity and provides additional information on possible geographic origin.

Dominance of multidrug resistant CC271 clones in macrolide-resistant streptococcus pneumoniae in Arizona.

BACKGROUND: Rates of resistance to macrolide antibiotics in Streptococcus pneumoniae are rising around the world due to the spread of mobile genetic elements harboring mef(E) and erm(B) genes and post-vaccine clonal expansion of strains that carry them. RESULTS: Characterization of 592 clinical isolates collected in Arizona over a 10 year period shows 23.6% are macrolide resistant. The largest portion of the macrolide-resistant population, 52%, is dual mef(E)/erm(B)-positive. All dual-positive isolates are multidrug-resistant clonal lineages of Taiwan19F-14, mostly multilocus sequence type 320, carrying the recently described transposon Tn2010. The remainder of the macrolide resistant S. pneumoniae collection includes 31% mef(E)-positive, and 9% erm(B)-positive strains. CONCLUSIONS: The dual-positive, multidrug-resistant S. pneumoniae clones have likely expanded by switching to non-vaccine serotypes after the heptavalent pneumococcal conjugate vaccine release, and their success limits therapy options. This upsurge could have a considerable clinical impact in Arizona.

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.

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.

Normalization of SARS-CoV-2 viral load via RT-qPCR provides higher-resolution data for comparison across time and between patients.

The 2020 pandemic has transformed the world and elicited thousands of studies to better understand the SARS-CoV-2 virus. Viral load has been a common measure to monitor treatment therapies and associate viral dynamics with patient outcomes; however, methods associated with viral load have varied across studies. These variations have the potential to sacrifice the accuracy of findings as they often do not account for inter-assay variation or variation across samples. In a retrospective study of nasopharyngeal samples, we found a significant amount of variation within the DNA and RNA targets; for example, across time within a single patient, there was an average of a 32-fold change. Further, we explore the impacts of host normalization on 94 clinical samples using the TGen Quantitative SARS-CoV-2 assay, finding that without host normalization samples with the same viral concentration can have up to 100-fold variation in the viral load.

Applying Genomic Epidemiology to Characterize a COVID-19 Outbreak in a Developmentally Disabled Adult Group Home Setting, Arizona.

Individuals living in congregate settings, including those in group homes, have been disproportionately impacted by COVID-19 and may be at increased risk of exposure or infection due to underlying illness. In mid-May 2020, local public health officials responded to an outbreak of COVID-19 among staff and residents associated with a multi-residential group home that provides care for adults with intellectual and developmental disabilities. Samples were collected at 16 of the homes. In four of the homes all the residents tested positive, and in the remaining 12 houses where samples were collected, all residents tested negative. Of the 152 individuals tested, 15/58 (25.9%) residents and 27/94 (28.7%) staff were positive for SARS-CoV-2, including eight hospitalizations and four deaths. Phylogenetic analysis of genomes from this outbreak in the context of genomes from Northern Arizona shows that very few mutations separate the samples from this outbreak. A potential transmission network was developed to illustrate person-place epidemiologic linkages and further demonstrates the dynamic connections between staff and residents with respect to each group home location. Epidemiologic and genomic evidence correlate, and suggest that asymptomatic infected staff likely introduced and spread COVID-19 in this setting. Implementation of public health prevention measures alongside rapid genomic analysis can help guide policy development and guide management efforts to prevent and mitigate future outbreaks.

Methods for sequencing the pandemic: benefits of rapid or high-throughput processing.

Genomic epidemiology has proven successful for real-time and retrospective monitoring of small and large-scale outbreaks. Here, we report two genomic sequencing and analysis strategies for rapid-turnaround or high-throughput processing of metagenomic samples. The rapid-turnaround method was designed to provide a quick phylogenetic snapshot of samples at the heart of active outbreaks, and has a total turnaround time of <48 hours from raw sample to analyzed data. The high-throughput method was designed for semi-retrospective data analysis, and is both cost effective and highly scalable. Though these methods were developed and utilized for the SARS-CoV-2 pandemic response in Arizona, U.S, and we envision their use for infectious disease epidemiology in the 21 st Century.