A Quantitative Assessment of Staphylococcus aureus Community Carriage in Yuma, Arizona.

BACKGROUND: Disease control relies on pathogen identification and understanding reservoirs. Staphylococcus aureus infection prevention is based upon decades of research on colonization and infection, but diminishing returns from mitigation efforts suggest significant knowledge gaps. Existing knowledge and mitigation protocols are founded upon culture-based detection, with almost no information about pathogen quantities. METHODS: We used culture and a quantitative polymerase chain reaction assay on samples from 3 body sites to characterize colonization more comprehensively than previous studies by describing both prevalence and pathogen quantity. RESULTS: We show a much higher overall prevalence (65.9%) than previously documented, with higher quantities and prevalence associated with the nares, non-Hispanic males (86.9%), and correlating with colonization in other body sites. These results suggest that research and clinical practices likely misclassify over half of colonized persons, limiting mitigation measures and their impact. CONCLUSIONS: This work begins the process of rebuilding foundational knowledge of S aureus carriage with more accurate and wholistic approaches.

Pathogen to commensal? Longitudinal within-host population dynamics, evolution, and adaptation during a chronic >16-year Burkholderia pseudomallei infection.

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314's sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314's airways over a >16-year period, and also sampled the patient's home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host.

SaQuant: a real-time PCR assay for quantitative assessment of Staphylococcus aureus.

BACKGROUND: Molecular assays are important tools for pathogen detection but need to be periodically re-evaluated with the discovery of additional genetic diversity that may cause assays to exclude target taxa or include non-target taxa. A single well-developed assay can find broad application across research, clinical, and industrial settings. Pathogen prevalence within a population is estimated using such assays and accurate results are critical for formulating effective public health policies and guiding future research. A variety of assays for the detection of Staphylococcus aureus are currently available. The utility of commercial assays for research is limited, given proprietary signatures and lack of transparent validation. RESULTS: In silico testing of existing peer-reviewed assays show that most suffer from a lack of sensitivity and specificity. We found no assays that were specifically designed and validated for quantitative use. Here we present a qPCR assay, SaQuant, for the detection and quantification of S. aureus as might be collected on sampling swabs. Sensitivity and specificity of the assay was 95.6 and 99.9 %, respectively, with a limit of detection of between 3 and 5 genome equivalents and a limit of quantification of 8.27 genome equivalents. The presence of DNA from non-target species likely to be found in a swab sample, did not impact qualitative or quantitative abilities of the assay. CONCLUSIONS: This assay has the potential to serve as a valuable tool for the accurate detection and quantification of S. aureus collected from human body sites in order to better understand the dynamics of prevalence and transmission in community settings.

Leptospira in river and soil in a highly endemic area of Ecuador.

BACKGROUND: Leptospira are shed into the environment via urine of infected animals. Rivers are thought to be an important risk factor for transmission to humans, though much is unknown about the types of environment or characteristics that favor survival. To address this, we screened for Leptospira DNA in two rivers in rural Ecuador where Leptospirosis is endemic. RESULTS: We collected 112 longitudinal samples and recorded pH, temperature, river depth, precipitation, and dissolved oxygen. We also performed a series of three experiments designed to provide insight into Leptospira presence in the soil. In the first soil experiment, we characterized prevalence and co-occurrence of Leptospira with other bacterial taxa in the soil at dispersed sites along the rivers (n = 64). In the second soil experiment, we collected 24 river samples and 48 soil samples at three points along eight transects to compare the likelihood of finding Leptospira in the river and on the shore at different distances from the river. In a third experiment, we tested whether Leptospira presence is associated with soil moisture by collecting 25 soil samples from two different sites. In our river experiment, we found pathogenic Leptospira in only 4 (3.7%) of samples. In contrast, pathogenic Leptospira species were found in 22% of shore soil at dispersed sites, 16.7% of soil samples (compared to 4.2% of river samples) in the transects, and 40% of soil samples to test for associations with soil moisture. CONCLUSIONS: Our data are limited to two sites in a highly endemic area, but the scarcity of Leptospira DNA in the river is not consistent with the widespread contention of the importance of river water for leptospirosis transmission. While Leptospira may be shed directly into the river, onto the shores, or washed into the river from more remote sites, massive dilution and limited persistence in rivers may reduce the environmental load and therefore, the epidemiological significance of such sources. It is also possible that transmission may occur more frequently on shores where people are liable to be barefoot. Molecular studies that further explore the role of rivers and water bodies in the epidemiology of leptospirosis are needed.

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples.

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.

Critical Knowledge Gaps in Our Understanding of Environmental Cycling and Transmission of Leptospira spp.

Exposure to soil or water contaminated with the urine of Leptospira-infected animals is the most common way in which humans contract leptospirosis. Entire populations can be at high risk of leptospirosis while working in inundated fields, when engaging in aquatic sports, or after periods of heavy rainfall. The risk of infection after contact with these environmental sources depends on the ability of Leptospira bacteria to survive, persist, and infect new hosts. Multiple variables such as soil and water pH, temperature, and even environmental microbial communities are likely to shape the environmental conditions needed by the pathogen to persist. Here we review what is known about the environmental phase of the infectious Leptospira transmission cycle and identify knowledge gaps that will serve as a guide for future research.

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.

Whole-genome sequencing of Burkholderia pseudomallei isolates from an unusual melioidosis case identifies a polyclonal infection with the same multilocus sequence type.

Twelve Burkholderia pseudomallei isolates collected over a 32-month period from a patient with chronic melioidosis demonstrated identical multilocus sequence types (STs). However, whole-genome sequencing suggests a polyclonal infection. This study is the first to report a mixed infection with the same ST.

Estimated herd prevalence and sequence types of Coxiella burnetii in bulk tank milk samples from commercial dairies in Indiana.

BACKGROUND: Coxiella burnetii is the etiologic agent of Q fever, a zoonotic disease causing influenza-like illness, pregnancy loss, cardiovascular disease and chronic fatigue syndrome in people. C. burnetii is considered to be enzootic in ruminants, but clinical signs of infection do not always manifest. National studies have documented the presence of C. burnetii in dairy herds in Indiana. This represents an opportunity to better characterize the distribution and prevalence of C. burnetii infection at the state scale, allowing evaluation of the need for surveillance and response planning to occur at this level. A cross-sectional study was conducted to estimate the herd prevalence of C. burnetii in commercial cattle dairies in Indiana and characterize the strains of C. burnetii within these dairies. RESULTS: Bulk tank milk samples were collected between June and August of 2011 by the Indiana State Board of Animal Health (ISBOAH). A total of 316 of these samples were tested for the IS1111 transposon of C. burnetii using quantitative real time polymerase chain reaction (PCR). Single nucleotide polymorphism (SNP) genotyping was used to identify the multispacer sequence genotypes (ST) present in samples where the IS1111 transposon was identified. The geographic distribution of dairies testing positive for C. burnetii DNA and the identified STs were also evaluated. The estimated overall herd prevalence for C. burnetii DNA was 61.1 % (95 % CI 55.6-66.3 %). The highest estimated regional prevalence was 70.2 % in the Central region of Indiana. An ST was identifiable in 74 of the positive 178 samples (41.6 %) and none of the 10 negative samples tested. Of these samples, 71 (95.9 %) were identified as ST20, 2 (2.7 %) as ST8 and a combination of ST20 and ST8 was identified in a single sample. CONCLUSIONS: C. burnetii is present in dairy herds throughout Indiana. Indiana follows national trends with ST20 most commonly identified. The presence of multiple STs in a single bulk tank sample indicates that multiple strains of C. burnetii can circulate within a herd. This supports potential transmission of C. burnetii between goats and cattle, presenting the potential for a switch in the dominant genotype found in a given species.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence.

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Francisella tularensis subsp. tularensis group A.I, United States.

We used whole-genome analysis and subsequent characterization of geographically diverse strains using new genetic signatures to identify distinct subgroups within Francisella tularensis subsp. tularensis group A.I: A.I.3, A.I.8, and A.I.12. These subgroups exhibit complex phylogeographic patterns within North America. The widest distribution was observed for A.I.12, which suggests an adaptive advantage.

High prevalence and two dominant host-specific genotypes of Coxiella burnetii in U.S. milk.

BACKGROUND: Coxiella burnetii causes Q fever in humans and Coxiellosis in animals; symptoms range from general malaise to fever, pneumonia, endocarditis and death. Livestock are a significant source of human infection as they shed C. burnetii cells in birth tissues, milk, urine and feces. Although prevalence of C. burnetii is high, few Q fever cases are reported in the U.S. and we have a limited understanding of their connectedness due to difficulties in genotyping. Here, we develop canonical SNP genotyping assays to evaluate spatial and temporal relationships among C. burnetii environmental samples and compare them across studies. Given the genotypic diversity of historical collections, we hypothesized that the current enzootic of Coxiellosis is caused by multiple circulating genotypes. We collected A) 23 milk samples from a single bovine herd, B) 134 commercial bovine and caprine milk samples from across the U.S., and C) 400 bovine and caprine samples from six milk processing plants over three years. RESULTS: We detected C. burnetii DNA in 96% of samples with no variance over time. We genotyped 88.5% of positive samples; bovine milk contained only a single genotype (ST20) and caprine milk was dominated by a second type (mostly ST8). CONCLUSIONS: The high prevalence and lack of genotypic diversity is consistent with a model of rapid spread and persistence. The segregation of genotypes between host species is indicative of species-specific adaptations or dissemination barriers and may offer insights into the relative lack of human cases and characterizing genotypes.

When outgroups fail; phylogenomics of rooting the emerging pathogen, Coxiella burnetii.

Rooting phylogenies is critical for understanding evolution, yet the importance, intricacies and difficulties of rooting are often overlooked. For rooting, polymorphic characters among the group of interest (ingroup) must be compared to those of a relative (outgroup) that diverged before the last common ancestor (LCA) of the ingroup. Problems arise if an outgroup does not exist, is unknown, or is so distant that few characters are shared, in which case duplicated genes originating before the LCA can be used as proxy outgroups to root diverse phylogenies. Here, we describe a genome-wide expansion of this technique that can be used to solve problems at the other end of the evolutionary scale: where ingroup individuals are all very closely related to each other, but the next closest relative is very distant. We used shared orthologous single nucleotide polymorphisms (SNPs) from 10 whole genome sequences of Coxiella burnetii, the causative agent of Q fever in humans, to create a robust, but unrooted phylogeny. To maximize the number of characters informative about the rooting, we searched entire genomes for polymorphic duplicated regions where orthologs of each paralog could be identified so that the paralogs could be used to root the tree. Recent radiations, such as those of emerging pathogens, often pose rooting challenges due to a lack of ingroup variation and large genomic differences with known outgroups. Using a phylogenomic approach, we created a robust, rooted phylogeny for C. burnetii. [Coxiella burnetii; paralog SNPs; pathogen evolution; phylogeny; recent radiation; root; rooting using duplicated genes.].

Genotyping of Coxiella burnetii from domestic ruminants and human in Hungary: indication of various genotypes.

BACKGROUND: Information about the genotypic characteristic of Coxiella burnetii from Hungary is lacking. The aim of this study is to describe the genetic diversity of C. burnetii in Hungary and compare genotypes with those found elsewhere. A total of 12 samples: (cattle, n = 6, sheep, n = 5 and human, n = 1) collected from across Hungary were studied by a 10-loci multispacer sequence typing (MST) and 6-loci multiple-locus variable-number of tandem repeat analysis (MLVA). Phylogenetic relationships among MST genotypes show how these Hungarian samples are related to others collected around the world. RESULTS: Three MST genotypes were identified: sequence type (ST) 20 has also been identified in ruminants from other European countries and the USA, ST28 was previously identified in Kazakhstan, and the proposed ST37 is novel. All MST genotypes yielded different MLVA genotypes and three different MLVA genotypes were identified within ST20 samples alone. Two novel MLVA types 0-9-5-5-6-2 (AG) and 0-8-4-5-6-2 (AF) (Ms23-Ms24-Ms27-Ms28-Ms33-Ms34) were defined in the ovine materials correlated with ST28 and ST37. Samples from different parts of the phylogenetic tree were associated with different hosts, suggesting host-specific adaptations. CONCLUSIONS: Even with the limited number of samples analysed, this study revealed high genetic diversity among C. burnetii in Hungary. Understanding the background genetic diversity will be essential in identifying and controlling outbreaks.

Staphylococcus aureus infection disparities among Hispanics and non-Hispanics in Yuma, Arizona.

Staphylococcus aureus infection patterns in Yuma, Arizona show a 2.25x higher infection rate in non-Hispanics. Males had higher infection rates in most age classes. These disparities in infection are mostly consistent with previously observed patterns in colonization, suggesting that sex and ethnicity do not differentially impact colonization and infection.

Population sequencing for diversity and transmission analyses.

Genomic diversity in a pathogen population is the foundation for evolution and adaptations in virulence, drug resistance, pathogenesis, and immune evasion. Characterizing, analyzing, and understanding population-level diversity is also essential for epidemiological and forensic tracking of sources and revealing detailed pathways of transmission and spread. For bacteria, culturing, isolating, and sequencing the large number of individual colonies required to adequately sample diversity can be prohibitively time-consuming and expensive. While sequencing directly from a mixed population will show variants among reads, they cannot be linked to reveal allele combinations associated with particular traits or phylogenetic inheritance patterns. Here, we describe the theory and method of how population sequencing directly from a mixed sample can be used in conjunction with sequencing a very small number of colonies to describe the phylogenetic diversity of a population without haplotype reconstruction. To demonstrate the utility of population sequencing in capturing phylogenetic diversity, we compared isogenic clones to population sequences of Burkholderia pseudomallei from the sputum of a single patient. We also analyzed population sequences of Staphylococcus aureus derived from different people and different body sites. Sequencing results confirm our ability to capture and characterize phylogenetic diversity in our samples. Our analyses of B. pseudomallei populations led to the surprising discovery that the pathogen population is highly structured in sputum, suggesting that for some pathogens, sputum sampling may preserve structuring in the lungs and thus present a non-invasive alternative to understanding colonization, movement, and pathogen/host interactions. Our analyses of S. aureus samples show how comparing phylogenetic diversity across populations can reveal directionality of transmission between hosts and across body sites, demonstrating the power and utility for characterizing the spread of disease and identification of reservoirs at the finest levels. We anticipate that population sequencing and analysis can be broadly applied to accelerate research in a broad range of fields reliant on a foundational understanding of population diversity.

Domestic dogs in indigenous Amazonian communities: Key players in Leptospira cycling and transmission?

BACKGROUND: Leptospirosis is the world's most common zoonotic disease. Mitigation and control rely on pathogen identification and understanding the roles of potential reservoirs in cycling and transmission. Underreporting and misdiagnosis obscure the magnitude of the problem and confound efforts to understand key epidemiological components. Difficulties in culturing hamper the use of serological diagnostics and delay the development of DNA detection methods. As a result, especially in complex ecosystems, we know very little about the importance of different mammalian host species in cycling and transmission to humans. METHODOLOGY/PRINCIPAL FINDINGS: We sampled dogs from five indigenous Kichwa communities living in the Yasuní National Park in the Ecuadorian Amazon basin. Blood and urine samples from domestic dogs were collected to assess the exposure of these animals to Leptospira and to identify the circulating species. Microscopic Agglutination Tests with a panel of 22 different serovars showed anti-leptospira antibodies in 36 sampled dogs (75%), and 7 serogroups were detected. Two DNA-based detection assays revealed pathogenic Leptospira DNA in 18 of 19 dog urine samples (94.7%). Amplicon sequencing and phylogenetic analysis of 16S rRNA and SecY genes from 15 urine samples revealed genetic diversity within two of three different Leptospira species: noguchii (n = 7), santarosai (n = 7), and interrogans (n = 1). CONCLUSIONS/SIGNIFICANCE: The high prevalence of antibodies and Leptospira DNA provides strong evidence for high rates of past and current infections. Such high prevalence has not been previously reported for dogs. These dogs live in the peridomestic environment in close contact with humans, yet they are free-ranging animals that interact with wildlife. This complex web of interactions may explain the diverse types of pathogenic Leptospira observed in this study. Our results suggest that domestic dogs are likely to play an important role in the cycling and transmission of Leptospira. Future studies in areas with complex ecoepidemiology will enable better parsing of the significance of genotypic, environmental, and host characteristics.

Natural reversion promotes LPS elongation in an attenuated Coxiella burnetii strain.

Lipopolysaccharide (LPS) phase variation is a critical aspect of virulence in many Gram-negative bacteria. It is of particular importance to Coxiella burnetii, the biothreat pathogen that causes Q fever, as in vitro propagation of this organism leads to LPS truncation, which is associated with an attenuated and exempted from select agent status (Nine Mile II, NMII). Here, we demonstrate that NMII was recovered from the spleens of infected guinea pigs. Moreover, these strains exhibit a previously unrecognized form of elongated LPS and display increased virulence in comparison with the initial NMII strain. The reversion of a 3-bp mutation in the gene cbu0533 directly leads to LPS elongation. To address potential safety concerns, we introduce a modified NMII strain unable to produce elongated LPS.

Host population structure and rare dispersal events drive leptospirosis transmission patterns among Rattus norvegicus in Boston, Massachusetts, US.

Leptospirosis (caused by pathogenic bacteria in the genus Leptospira ) is prevalent worldwide but more common in tropical and subtropical regions. Transmission can occur following direct exposure to infected urine from reservoir hosts, such as rats, or a urine-contaminated environment, which then can serve as an infection source for additional rats and other mammals, including humans. The brown rat, Rattus norvegicus , is an important reservoir of leptospirosis in urban settings. We investigated leptospirosis among brown rats in Boston, Massachusetts and hypothesized that rat dispersal in this urban setting influences the movement, persistence, and diversity of Leptospira . We analyzed DNA from 328 rat kidney samples collected from 17 sites in Boston over a seven-year period (2016-2022); 59 rats representing 12 of 17 sites were positive for Leptospira . We used 21 neutral microsatellite loci to genotype 311 rats and utilized the resulting data to investigate genetic connectivity among sampling sites. We generated whole genome sequences for 28 Leptospira isolates obtained from frozen and fresh tissue from some of the 59 Leptospira -positive rat kidneys. When isolates were not obtained, we attempted Leptospira genomic DNA capture and enrichment, which yielded 14 additional Leptospira genomes from rats. We also generated an enriched Leptospira genome from a 2018 human case in Boston. We found evidence of high genetic structure and limited dispersal among rat populations that is likely influenced by major roads and/or other unknown dispersal barriers, resulting in distinct rat population groups within the city; at certain sites these groups persisted for multiple years. We identified multiple distinct phylogenetic clades of L. interrogans among rats, with specific clades tightly linked to distinct rat populations. This pattern suggests L. interrogans persists in local rat populations and movement of leptospirosis in this urban rat community is driven by rat dispersal. Finally, our genomic analyses of the 2018 human leptospirosis case in Boston suggests a link to rats as the source. These findings will be useful for guiding rat control and human leptospirosis mitigation efforts in this and other urban settings.

Pan-Enterovirus Characterization Reveals Cryptic Circulation of Clinically Relevant Subtypes in Arizona Wastewater.

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.