Background:Most seasonally circulating enteroviruses result in asymptomatic or mildly symptomatic infections. In rare cases, however, infection with some subtypes can result in paralysis or death. Of the 300 subtypes known, only poliovirus is reportable, limiting our understanding of the distribution of other enteroviruses that can cause clinical disease.Objective:The overarching objectives of this study were to: 1) describe the distribution of enteroviruses in Arizona during the late summer and fall of 2022, the time of year when they are thought to be most abundant, and 2) demonstrate the utility of viral pan-assay approaches for semi-agnostic discovery that can be followed up by more targeted assays and phylogenomics.Methods:This study utilizes pooled nasal samples collected from school-aged children and long-term care facility residents, and wastewater from multiple locations in Arizona during July-October of 2022. We used PCR to amplify and sequence a region common to all enteroviruses, followed by species-level bioinformatic characterization using the QIIME 2 platform. For Enterovirus-D68 (EV-D68), detection was carried out using RT-qPCR, followed by confirmation using near-complete whole EV-D68 genome sequencing using a newly designed tiled amplicon approach.Results:In the late summer and early fall of 2022, multiple enterovirus species were identified in Arizona wastewater, with Coxsackievirus A6, EV-D68, and Coxsackievirus A19 composing 86% of the characterized reads sequenced. While EV-D68 was not identified in pooled human nasal samples, and the only reported acute flaccid myelitis case in Arizona did not test positive for the virus, an in-depth analysis of EV-D68 in wastewater revealed that the virus was circulating from August through mid-October. A phylogenetic analysis on this relatively limited dataset revealed just a few importations into the state, with a single clade indicating local circulation.Significance:This study further supports the utility of wastewater-based epidemiology to identify potential public health threats. Our further investigations into EV-D68 shows how these data might help inform healthcare diagnoses for children presenting with concerning neurological symptoms.

Our goal for this project was to incorporate and operationalize a lean laboratory processing method that decreases costs, time, and complexity and increases throughput and accessibility for SARS-CoV-2 tiled amplicon sequencing.

Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 25, No. 1, January 2019 Increasingly, routine surveillance and monitoring of foodborne pathogens using whole-genome sequencing is creating opportunities to study foodborne illness epidemiology beyond routine outbreak investigations and case–control studies. Using a global phylogeny of Salmonella enterica serotype Typhimurium, we found that major livestock sources of the pathogen in the United States can be predicted through whole-genome sequencing data. Relatively steady rates of sequence divergence in livestock lineages enabled the inference of their recent origins. Elevated accumulation of lineage-specific pseudogenes after divergence from generalist populations and possible metabolic acclimation in a representative swine isolate indicates possible emergence of host adaptation. We developed and retrospectively applied a machine learning Random Forest classifier for genomic source prediction of Salmonella Typhimurium that correctly attributed 7 of 8 major zoonotic outbreaks in the United States during 1998–2013. We further identified 50 key genetic features that were sufficient for robust livestock source prediction.

Antimicrobial resistance (AMR) in the nosocomial pathogen,Acinetobacter baumannii, is becoming a serious public health threat. While some mechanisms of AMR have been reported, understanding novel mechanisms of resistance is critical for identifying emerging resistance. One of the first steps in identifying novel AMR mechanisms is performing genotype/phenotype association studies; however, performing these studies is complicated by the plastic nature of theA. baumanniipan-genome. In this study, we compared the antibiograms of 12 antimicrobials associated with multiple drug families for 84A. baumanniiisolates, many isolated in Arizona, USA.in silicoscreening of these genomes for known AMR mechanisms failed to identify clear correlations for most drugs. We then performed a bacterial genome wide association study (bGWAS) looking for associations between all possible 21-mers; this approach generally failed to identify mechanisms that explained the resistance phenotype. In order to decrease the genomic noise associated with population stratification, we compared four phylogenetically-related pairs of isolates with differing susceptibility profiles. RNA-Sequencing (RNA-Seq) was performed on paired isolates and differentially-expressed genes were identified. In these isolate pairs, five different potential mechanisms were identified, highlighting the difficulty of broad AMR surveillance in this species. To verify and validate differential expression, amplicon sequencing was performed. These results suggest that a diagnostic platform based on gene expression rather than genomics alone may be beneficial in certain surveillance efforts. The implementation of such advanced diagnostics coupled with increased AMR surveillance will potentially improveA. baumanniiinfection treatment and patient outcomes.

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40% among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus. To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

Increasingly, routine surveillance and monitoring of foodborne pathogens using whole-genome sequencing is creating opportunities to study foodborne illness epidemiology beyond routine outbreak investigations and case-control studies. Using a global phylogeny of Salmonella enterica serotype Typhimurium, we found that major livestock sources of the pathogen in the United States can be predicted through whole-genome sequencing data. Relatively steady rates of sequence divergence in livestock lineages enabled the inference of their recent origins. Elevated accumulation of lineage-specific pseudogenes after divergence from generalist populations and possible metabolic acclimation in a representative swine isolate indicates possible emergence of host adaptation. We developed and retrospectively applied a machine learning Random Forest classifier for genomic source prediction of Salmonella Typhimurium that correctly attributed 7 of 8 major zoonotic outbreaks in the United States during 1998-2013. We further identified 50 key genetic features that were sufficient for robust livestock source prediction.

Strains of Staphylococcus aureus in clonal complex 8 (CC8), including USA300, USA500, and the Iberian clone, are prevalent pathogens in the United States, both inside and outside health care settings. Methods for typing CC8 strains are becoming obsolete as the strains evolve and diversify, and whole-genome sequencing has shown that some strain types fall into multiple sublineages within CC8. In this study, we attempt to clarify the strain nomenclature of CC8, classifying the major strain types based on whole-genome sequence phylogenetics using both methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) genomes. We show that isolates of the Archaic and Iberian clones from decades ago make up the most basal clade of the main CC8 lineages and that at least one successful lineage of CC8, made up mostly of MSSA, diverged before the other well-known strain types USA500 and USA300. We also show that the USA500 type includes two clades separated by the previously described "Canadian epidemic MRSA" strain CMRSA9, that one clade containing USA500 also contains the USA300 clade, and that the USA300-0114 strain type is not a monophyletic group. Additionally, we present a rapid, simple CC8 strain-typing scheme using real-time PCR assays that target single nucleotide polymorphisms (SNPs) derived from our CC8 phylogeny and show the significant benefit of using more stable genomic markers based on evolutionary lineages over traditional S. aureus typing techniques. This more accurate and accessible S. aureus typing system may improve surveillance and better inform the epidemiology of this very important pathogen. IMPORTANCE Staphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus.

Whole-genome sequencing (WGS) of bacterial isolates has become standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. NASP was developed as a reproducible method that scales well with the hundreds to thousands of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares with other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces similar, and often better, results in comparison with other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

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.

Whole genome sequencing (WGS) of bacteria is becoming standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. The Northern Arizona SNP Pipeline (NASP) was developed as a reproducible pipeline that scales well with the large amount of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares to other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces comparable, and often better, results to other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools, and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

Healthcare-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to healthcare 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 become an urgent public health crisis. In the healthcare environment, patient colonization of K. pneumoniae precedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening for K. pneumoniae colonization has cost-effective and life-saving potential. In this study, we describe the design and validation of KlebSeq, a highly informative screening tool that detects Klebsiella species and identifies clinically important strains and characteristics 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, several types of respiratory, and fecal, showing K. pneumoniae species and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool can be used to screen patients for K. pneumoniae carriage to assess risk of infection and outbreak potential, and the expansion of this tool can be used for several other HAI agents or applications.

We used whole-genome sequence typing (WGST) to investigate an outbreak of Sarocladium kiliense bloodstream infections (BSI) associated with receipt of contaminated antinausea medication among oncology patients in Colombia and Chile during 2013-2014. Twenty-five outbreak isolates (18 from patients and 7 from medication vials) and 11 control isolates unrelated to this outbreak were subjected to WGST to elucidate a source of infection. All outbreak isolates were nearly indistinguishable (<5 single-nucleotide polymorphisms), and >21,000 single-nucleotide polymorphisms were identified from unrelated control isolates, suggesting a point source for this outbreak. S. kiliense has been previously implicated in healthcare-related infections; however, the lack of available typing methods has precluded the ability to substantiate point sources. WGST for outbreak investigation caused by eukaryotic pathogens without reference genomes or existing genotyping methods enables accurate source identification to guide implementation of appropriate control and prevention measures.

ABSTRACT Coccidioidomycosis (or valley fever) is a fungal disease with high morbidity and mortality that affects tens of thousands of people each year. This infection is caused by two sibling species, Coccidioides immitis and C. posadasii , which are endemic to specific arid locales throughout the Western Hemisphere, particularly the desert southwest of the United States. Recent epidemiological and population genetic data suggest that the geographic range of coccidioidomycosis is expanding, as new endemic clusters have been identified in the state of Washington, well outside the established endemic range. The genetic mechanisms and epidemiological consequences of this expansion are unknown and require better understanding of the population structure and evolutionary history of these pathogens. Here we performed multiple phylogenetic inference and population genomics analyses of 68 new and 18 previously published genomes. The results provide evidence of substantial population structure in C. posadasii and demonstrate the presence of distinct geographic clades in central and southern Arizona as well as dispersed populations in Texas, Mexico, South America, and Central America. Although a smaller number of C. immitis strains were included in the analyses, some evidence of phylogeographic structure was also detected in this species, which has been historically limited to California and Baja, Mexico. Bayesian analyses indicated that C. posadasii is the more ancient of the two species and that Arizona contains the most diverse subpopulations. We propose a southern Arizona-northern Mexico origin for C. posadasii and describe a pathway for dispersal and distribution out of this region. IMPORTANCE Coccidioidomycosis, or valley fever, is caused by the pathogenic fungi Coccidioides posadasii and C. immitis . The fungal species and disease are primarily found in the American desert southwest, with spotted distribution throughout the Western Hemisphere. Initial molecular studies suggested a likely anthropogenic movement of C. posadasii from North America to South America. Here we comparatively analyze eighty-six genomes of the two Coccidioides species and establish local and species-wide population structures to not only clarify the earlier dispersal hypothesis but also provide evidence of likely ancestral populations and patterns of dispersal for the known subpopulations of C. posadasii .

Background Prevention of nosocomial transmission of infections is a central responsibility in the healthcare environment, and accurate identification of transmission events presents the first challenge. Phylogenetic analysis based on whole genome sequencing provides a high-resolution approach for accurately relating isolates to one another, allowing precise identification or exclusion of transmission events and sources for nearly all cases. We sequenced 24 methicillin-resistant Staphylococcus aureus (MRSA) genomes to retrospectively investigate a suspected point source of three surgical site infections (SSIs) that occurred over a one-year period. The source of transmission was believed to be a surgical team member colonized with MRSA, involved in all surgeries preceding the SSI cases, who was subsequently decolonized. Genetic relatedness among isolates was determined using whole genome single nucleotide polymorphism (SNP) data. Results Whole genome SNP typing (WGST) revealed 283 informative SNPs between the surgical team member’s isolate and the closest SSI isolate. The second isolate was 286 and the third was thousands of SNPs different, indicating the nasal carriage strain from the surgical team member was not the source of the SSIs. Given the mutation rates estimated for S. aureus, none of the SSI isolates share a common ancestor within the past 16 years, further discounting any common point source for these infections. The decolonization procedures and resources spent on the point source infection control could have been prevented if WGST was performed at the time of the suspected transmission, instead of retrospectively. Conclusions Whole genome sequence analysis is an ideal method to exclude isolates involved in transmission events and nosocomial outbreaks, and coupling this method with epidemiological data can determine if a transmission event occurred. These methods promise to direct infection control resources more appropriately.

The hyper-virulent emm59 genotype of invasive group A Streptococcus was identified in northern Arizona in 2015. Eighteen isolates belonging to a genomic cluster grouped most closely with recently identified isolates in New Mexico. The continued transmission of emm59 in the southwestern United States poses a public health concern.

Multidrug-resistant Klebsiella pneumoniae producing the KPC carbapenemase have rapidly spread throughout the world, causing severe healthcare-associated infections with limited antimicrobial treatment options. Dissemination of KPC-producing K. pneumoniae is largely attributed to expansion of a single dominant strain, ST258. In this study, we explore phylogenetic relationships and evolution within ST258 and its clonal group, CG258, using whole genome sequence analysis of 167 isolates from 20 countries collected over 17 years. Our results show a common ST258 ancestor emerged from its diverse parental clonal group around 1995 and likely acquired bla KPC prior to dissemination. Over the past two decades, ST258 has remained highly clonal despite diversity in accessory elements and divergence in the capsule polysaccharide synthesis locus. Apart from the large recombination event that gave rise to ST258, few mutations set it apart from its clonal group. However, one mutation occurs in a global transcription regulator. Characterization of outer membrane protein sequences revealed a profile in ST258 that includes a truncated OmpK35 and modified OmpK37. Our work illuminates potential genomic contributors to the pathogenic success of ST258, helps us better understand the global dissemination of this strain, and identifies genetic markers unique to ST258.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

ABSTRACT Highly invasive, community-acquired Klebsiella pneumoniae infections have recently emerged, resulting in pyogenic liver abscesses. These infections are caused by hypervirulent K. pneumoniae (hvKP) isolates primarily of capsule serotype K1 or K2. Hypervirulent K1 isolates belong to clonal complex 23 (CC23), indicating that this clonal lineage has a specific genetic background conferring hypervirulence. Here, we apply whole-genome sequencing to a collection of K. pneumoniae isolates to characterize the phylogenetic background of hvKP isolates with an emphasis on CC23. Most of the hvKP isolates belonged to CC23 and grouped into a distinct monophyletic clade, revealing that CC23 is a unique clonal lineage, clearly distinct from nonhypervirulent strains. Separate phylogenetic analyses of the CC23 isolates indicated that the CC23 lineage evolved recently by clonal expansion from a single common ancestor. Limited grouping according to geographical origin was observed, suggesting that CC23 has spread globally through multiple international transmissions. Conversely, hypervirulent K2 strains clustered in genetically unrelated groups. Strikingly, homologues of a large virulence plasmid were detected in all hvKP clonal lineages, indicating a key role in K. pneumoniae hypervirulence. The plasmid encodes two siderophores, aerobactin and salmochelin, and RmpA (regulator of the mucoid phenotype); all these factors were found to be restricted to hvKP isolates. Genomic comparisons revealed additional factors specifically associated with CC23. These included a distinct variant of a genomic island encoding yersiniabactin, colibactin, and microcin E492. Furthermore, additional novel genomic regions unique to CC23 were revealed which may also be involved in the increased virulence of this important clonal lineage. IMPORTANCE During the last 3 decades, hypervirulent Klebsiella pneumoniae (hvKP) isolates have emerged, causing severe community-acquired infections primarily in the form of pyogenic liver abscesses. This syndrome has so far primarily been found in Southeast Asia, but increasing numbers of cases are being reported worldwide, indicating that the syndrome is turning into a globally emerging disease. We applied whole-genome sequencing to a collection of K. pneumoniae clinical isolates to reveal the phylogenetic background of hvKP and to identify genetic factors associated with the increased virulence. The hvKP isolates primarily belonged to clonal complex 23 (CC23), and this clonal lineage was revealed to be clearly distinct from nonhypervirulent strains. A specific virulence plasmid was found to be associated with hypervirulence, and novel genetic determinants uniquely associated with CC23 were identified. Our findings extend the understanding of the genetic background of the emergence of hvKP clones.