Expanding the Burkholderia pseudomallei Complex with the Addition of Two Novel Species: Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov.

Distinct Burkholderia strains were isolated from soil samples collected in tropical northern Australia (Northern Territory and the Torres Strait Islands, Queensland). Phylogenetic analysis of 16S rRNA and whole genome sequences revealed these strains were distinct from previously described Burkholderia species and assigned them to two novel clades within the B. pseudomallei complex (Bpc). Because average nucleotide identity and digital DNA-DNA hybridization calculations are consistent with these clades representing distinct species, we propose the names Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov. Strains assigned to B. mayonis sp. nov. include type strain BDU6T (=TSD-80; LMG 29941; ASM152374v2) and BDU8. Strains assigned to B. savannae sp. nov. include type strain MSMB266T (=TSD-82; LMG 29940; ASM152444v2), MSMB852, BDU18, and BDU19. Comparative genomics revealed unique coding regions for both putative species, including clusters of orthologous genes associated with phage. Type strains of both B. mayonis sp. nov. and B. savannae sp. nov. yielded biochemical profiles distinct from each other and from other species in the Bpc, and profiles also varied among strains within B. mayonis sp. nov. and B. savannae sp. nov. Matrix-assisted laser desorption ionization time-of-flight (MLST) analysis revealed a B. savannae sp. nov. cluster separate from other species, whereas B. mayonis sp. nov. strains did not form a distinct cluster. Neither B. mayonis sp. nov. nor B. savannae sp. nov. caused mortality in mice when delivered via the subcutaneous route. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species currently within the Bpc. IMPORTANCEBurkholderia species can be important sources of novel natural products, and new species are of interest to diverse scientific disciplines. Although many Burkholderia species are saprophytic, Burkholderia pseudomallei is the causative agent of the disease melioidosis. Understanding the genomics and virulence of the closest relatives to B. pseudomallei, i.e., the other species within the B. pseudomallei complex (Bpc), is important for identifying robust diagnostic targets specific to B. pseudomallei and for understanding the evolution of virulence in B. pseudomallei. Two proposed novel species, B. mayonis sp. nov. and B. savannae sp. nov., were isolated from soil samples collected from multiple locations in northern Australia. The two proposed species belong to the Bpc but are phylogenetically distinct from all other members of this complex. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species within this significant complex of bacteria that are available for future studies.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study.

Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patient's complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

Necrotizing cutaneous mucormycosis after a tornado in Joplin, Missouri, in 2011.

BACKGROUND: Mucormycosis is a fungal infection caused by environmentally acquired molds. We investigated a cluster of cases of cutaneous mucormycosis among persons injured during the May 22, 2011, tornado in Joplin, Missouri. METHODS: We defined a case as a soft-tissue infection in a person injured during the tornado, with evidence of a mucormycete on culture or immunohistochemical testing plus DNA sequencing. We conducted a case-control study by reviewing medical records and conducting interviews with case patients and hospitalized controls. DNA sequencing and whole-genome sequencing were performed on clinical specimens to identify species and assess strain-level differences, respectively. RESULTS: A total of 13 case patients were identified, 5 of whom (38%) died. The patients had a median of 5 wounds (range, 1 to 7); 11 patients (85%) had at least one fracture, 9 (69%) had blunt trauma, and 5 (38%) had penetrating trauma. All case patients had been located in the zone that sustained the most severe damage during the tornado. On multivariate analysis, infection was associated with penetrating trauma (adjusted odds ratio for case patients vs. controls, 8.8; 95% confidence interval [CI], 1.1 to 69.2) and an increased number of wounds (adjusted odds ratio, 2.0 for each additional wound; 95% CI, 1.2 to 3.2). Sequencing of the D1-D2 region of the 28S ribosomal DNA yielded Apophysomyces trapeziformis in all 13 case patients. Whole-genome sequencing showed that the apophysomyces isolates were four separate strains. CONCLUSIONS: We report a cluster of cases of cutaneous mucormycosis among Joplin tornado survivors that were associated with substantial morbidity and mortality. Increased awareness of fungi as a cause of necrotizing soft-tissue infections after a natural disaster is warranted.

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.

Whole-genome analysis of Exserohilum rostratum from an outbreak of fungal meningitis and other infections.

Exserohilum rostratum was the cause of most cases of fungal meningitis and other infections associated with the injection of contaminated methylprednisolone acetate produced by the New England Compounding Center (NECC). Until this outbreak, very few human cases of Exserohilum infection had been reported, and very little was known about this dematiaceous fungus, which usually infects plants. Here, we report using whole-genome sequencing (WGS) for the detection of single nucleotide polymorphisms (SNPs) and phylogenetic analysis to investigate the molecular origin of the outbreak using 22 isolates of E. rostratum retrieved from 19 case patients with meningitis or epidural/spinal abscesses, 6 isolates from contaminated NECC vials, and 7 isolates unrelated to the outbreak. Our analysis indicates that all 28 isolates associated with the outbreak had nearly identical genomes of 33.8 Mb. A total of 8 SNPs were detected among the outbreak genomes, with no more than 2 SNPs separating any 2 of the 28 genomes. The outbreak genomes were separated from the next most closely related control strain by ∼136,000 SNPs. We also observed significant genomic variability among strains unrelated to the outbreak, which may suggest the possibility of cryptic speciation in E. rostratum.

Real-time PCR assays for genotyping of Cryptococcus gattii in North America.

BACKGROUND: Cryptococcus gattii has been the cause of an ongoing outbreak starting in 1999 on Vancouver Island, British Columbia and spreading to mainland Canada and the US Pacific Northwest. In the course of the outbreak, C. gattii has been identified outside of its previously documented climate, habitat, and host disease. Genotyping of C. gattii is essential to understand the ecological and geographical expansion of this emerging pathogen. METHODS: We developed and validated a mismatch amplification mutation assay (MAMA) real-time PCR panel for genotyping C. gattii molecular types VGI-VGIV and VGII subtypes a,b,c. Subtype assays were designed based on whole-genome sequence of 20 C. gattii strains. Publically available multilocus sequence typing (MLST) data from a study of 202 strains was used for the molecular type (VGI-VGIV) assay design. All assays were validated across DNA from 112 strains of diverse international origin and sample types, including animal, environmental and human. RESULTS: Validation revealed each assay on the panel is 100% sensitive, specific and concordant with MLST. The assay panel can detect down to 0.5 picograms of template DNA. CONCLUSIONS: The (MAMA) real-time PCR panel for C. gattii accurately typed a collection of 112 diverse strains and demonstrated high sensitivity. This is a time and cost efficient method of genotyping C. gattii best suited for application in large-scale epidemiological studies.

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.].

Rapid and robust phylotyping of spa t003, a dominant MRSA clone in Luxembourg and other European countries.

BACKGROUND: spa typing is a common genotyping tool for methicillin-resistant Staphylococcus aureus (MRSA) in Europe. Given the high prevalence of dominant clones, spa-typing is proving to be limited in its ability to distinguish outbreak isolates from background isolates. New molecular tools need to be employed to improve subtyping of dominant local MRSA strains (e.g., spa type t003). METHODS: Phylogenetically critical, or canonical, SNPs (can-SNPs) were identified as subtyping targets through sequence analysis of 40 MRSA whole genomes from Luxembourg. Real-time PCR assays were designed around target SNPs and validated using a repository of 240 previously sub-typed and epidemiologically characterized Luxembourg MRSA isolates, including 153 community and hospital isolates, 69 isolates from long term care (LTC) facilities, and 21 prospectively analyzed MRSA isolates. Selected isolates were also analyzed by whole genome SNP typing (WGST) for comparison to the SNP assays and other subtyping techniques. RESULTS: Fourteen real-time PCR assays were developed and validated, including two assays to determine presence of spa t003 or t008. The other twelve assays successfully provided a high degree of resolution within the t003 subtype. WGST analysis of the LTC facility isolates provided greater resolution than other subtyping tools, identifying clusters indicative of ongoing transmission within LTC facilities. CONCLUSIONS: canSNP-based PCR assays are useful for local level MRSA phylotyping, especially in the presence of one or more dominant clones. The assays designed here can be easily adapted for investigating t003 MRSA strains in other regions in Western Europe. WGST provides substantially better resolution than other typing methods.

Genome sequence of "Candidatus Microthrix parvicella" Bio17-1, a long-chain-fatty-acid-accumulating filamentous actinobacterium from a biological wastewater treatment plant.

"Candidatus Microthrix" bacteria are deeply branching filamentous actinobacteria which occur at the water-air interface of biological wastewater treatment plants, where they are often responsible for foaming and bulking. Here, we report the first draft genome sequence of a strain from this genus: "Candidatus Microthrix parvicella" strain Bio17-1.

Molecular epidemiologic investigation of an anthrax outbreak among heroin users, Europe.

In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis-positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

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.

An attenuated strain of Bacillus anthracis (CDC 684) has a large chromosomal inversion and altered growth kinetics.

BACKGROUND: An isolate originally labeled Bacillus megaterium CDC 684 was found to contain both pXO1 and pXO2, was non-hemolytic, sensitive to gamma-phage, and produced both the protective antigen and the poly-D-glutamic acid capsule. These phenotypes prompted Ezzell et al., (J. Clin. Microbiol. 28:223) to reclassify this isolate to Bacillus anthracis in 1990. RESULTS: We demonstrate that despite these B. anthracis features, the isolate is severely attenuated in a guinea pig model. This prompted whole genome sequencing and closure. The comparative analysis of CDC 684 to other sequenced B. anthracis isolates and further analysis reveals: a) CDC 684 is a close relative of a virulent strain, Vollum A0488; b) CDC 684 defines a new B. anthracis lineage (at least 51 SNPs) that includes 15 other isolates; c) the genome of CDC 684 contains a large chromosomal inversion that spans 3.3 Mbp; d) this inversion has caused a displacement of the usual spatial orientation of the origin of replication (ori) to the termination of replication (ter) from 180° in wild-type B. anthracis to 120° in CDC 684 and e) this isolate also has altered growth kinetics in liquid media. CONCLUSIONS: We propose two alternative hypotheses explaining the attenuated phenotype of this isolate. Hypothesis 1 suggests that the skewed ori/ter relationship in CDC 684 has altered its DNA replication and/or transcriptome processes resulting in altered growth kinetics and virulence capacity. Hypothesis 2 suggests that one or more of the single nucleotide polymorphisms in CDC 684 has altered the expression of a regulatory element or other genes necessary for virulence.

Next-generation sequencing of Coccidioides immitis isolated during cluster investigation.

Next-generation sequencing enables use of whole-genome sequence typing (WGST) as a viable and discriminatory tool for genotyping and molecular epidemiologic analysis. We used WGST to confirm the linkage of a cluster of Coccidioides immitis isolates from 3 patients who received organ transplants from a single donor who later had positive test results for coccidioidomycosis. Isolates from the 3 patients were nearly genetically identical (a total of 3 single-nucleotide polymorphisms identified among them), thereby demonstrating direct descent of the 3 isolates from an original isolate. We used WGST to demonstrate the genotypic relatedness of C. immitis isolates that were also epidemiologically linked. Thus, WGST offers unique benefits to public health for investigation of clusters considered to be linked to a single source.

Cyclopropavir inhibits the normal function of the human cytomegalovirus UL97 kinase.

Cyclopropavir (CPV) is active against human cytomegalovirus (CMV), as well as both variants of human herpesvirus 6 and human herpesvirus 8. The mechanism of action of CPV against CMV is similar to that of ganciclovir (GCV) in that it is phosphorylated initially by the CMV UL97 kinase, resulting in inhibition of viral DNA synthesis. Resistance to CPV maps to the UL97 kinase but is associated primarily with H520Q mutations and thus retains good antiviral activity against most GCV-resistant isolates. An examination of CMV-infected cultures treated with CPV revealed unusual cell morphology typically associated with the absence of UL97 kinase activity. A surrogate assay for UL97 kinase activity confirmed that CPV inhibited the activity of this enzyme and that its action was similar to the inhibition seen with maribavir (MBV) in this assay. Combination studies using real-time PCR indicated that, like MBV, CPV also antagonized the efficacy of GCV and were consistent with the observed inhibition of the UL97 kinase. Deep sequencing of CPV-resistant laboratory isolates identified a frameshift mutation in UL27, presumably to compensate for a loss of UL97 enzymatic activity. We conclude that the mechanism of action of CPV against CMV is complex and involves both the inhibition of DNA synthesis and the inhibition of the normal activity of the UL97 kinase.

Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics.

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

Cost-effective interrogation of single nucleotide polymorphisms using the mismatch amplification mutation assay and capillary electrophoresis.

The ability to characterize SNPs is an important aspect of many clinical diagnostic, genetic and evolutionary studies. Here, we designed a multiplexed SNP genotyping method to survey a large number of phylogenetically informative SNPs within the genome of the bacterium Bacillus anthracis. This novel method, CE universal tail mismatch amplification mutation assay (CUMA), allows for PCR multiplexing and automatic scoring of SNP genotypes, thus providing a rapid, economical and higher throughput alternative to more expensive SNP genotyping techniques. CUMA delivered accurate B. anthracis SNP genotyping results and, when multiplexed, saved reagent costs by more than 80% compared with TaqMan real-time PCR. When real-time PCR technology and instrumentation is unavailable or the reagents are cost-prohibitive, CUMA is a powerful alternative for SNP genotyping.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

Bacillus anthracis in China and its relationship to worldwide lineages.

BACKGROUND: The global pattern of distribution of 1033 B. anthracis isolates has previously been defined by a set of 12 conserved canonical single nucleotide polymorphisms (canSNP). These studies reinforced the presence of three major lineages and 12 sub-lineages and sub-groups of this anthrax-causing pathogen. Isolates that form the A lineage (unlike the B and C lineages) have become widely dispersed throughout the world and form the basis for the geographical disposition of "modern" anthrax. An archival collection of 191 different B. anthracis isolates from China provides a glimpse into the possible role of Chinese trade and commerce in the spread of certain sub-lineages of this pathogen. Canonical single nucleotide polymorphism (canSNP) and multiple locus VNTR analysis (MLVA) typing has been used to examine this archival collection of isolates. RESULTS: The canSNP study indicates that there are 5 different sub-lineages/sub-groups in China out of 12 previously described world-wide canSNP genotypes. Three of these canSNP genotypes were only found in the western-most province of China, Xinjiang. These genotypes were A.Br.008/009, a sub-group that is spread across most of Europe and Asia; A.Br.Aust 94, a sub-lineage that is present in Europe and India, and A.Br.Vollum, a lineage that is also present in Europe. The remaining two canSNP genotypes are spread across the whole of China and belong to sub-group A.Br.001/002 and the A.Br.Ames sub-lineage, two closely related genotypes. MLVA typing adds resolution to the isolates in each canSNP genotype and diversity indices for the A.Br.008/009 and A.Br.001/002 sub-groups suggest that these represent older and established clades in China. CONCLUSION: B. anthracis isolates were recovered from three canSNP sub-groups (A.Br.008/009, A.Br.Aust94, and A.Br.Vollum) in the western most portion of the large Chinese province of Xinjiang. The city of Kashi in this province appears to have served as a crossroads for not only trade but the movement of diseases such as anthrax along the ancient "silk road". Phylogenetic inference also suggests that the A.Br.Ames sub-lineage, first identified in the original Ames strain isolated from Jim Hogg County, TX, is descended from the A.Br.001/002 sub-group that has a major presence in most of China. These results suggest a genetic discontinuity between the younger Ames sub-lineage in Texas and the large Western North American sub-lineage spread across central Canada and the Dakotas.

Genomic Analyses of Acute Flaccid Myelitis Cases among a Cluster in Arizona Provide Further Evidence of Enterovirus D68 Role.

Enteroviruses are a common cause of respiratory and gastrointestinal illness, and multiple subtypes, including poliovirus, can cause neurologic disease. In recent years, enterovirus D68 (EV-D68) has been associated with serious neurologic illnesses, including acute flaccid myelitis (AFM), frequently preceded by respiratory disease. A cluster of 11 suspect cases of pediatric AFM was identified in September 2016 in Phoenix, AZ. To determine if these cases were associated with EV-D68, we performed multiple genomic analyses of nasopharyngeal (NP) swabs and cerebrospinal fluid (CSF) material from the patients, including real-time PCR and amplicon sequencing targeting the EV-D68 VP1 gene and unbiased microbiome and metagenomic sequencing. Four of the 11 patients were classified as confirmed cases of AFM, and an additional case was classified as probable AFM. Real-time PCR and amplicon sequencing detected EV-D68 virus RNA in the three AFM patients from which NP swabs were collected, as well as in a fourth patient diagnosed with acute disseminated encephalomyelitis, a disease that commonly follows bacterial or viral infections, including enterovirus. No other obvious etiological causes for AFM were identified by 16S or RNA and DNA metagenomic sequencing in these cases, strengthening the likelihood that EV-D68 is an etiological factor. Herpes simplex viral DNA was detected in the CSF of the fourth case of AFM and in one additional suspect case from the cluster. Multiple genomic techniques, such as those described here, can be used to diagnose patients with suspected EV-D68 respiratory illness, to aid in AFM diagnosis, and for future EV-D68 surveillance and epidemiology.IMPORTANCE Enteroviruses frequently result in respiratory and gastrointestinal illness; however, multiple subtypes, including poliovirus, can cause severe neurologic disease. Recent biennial increases (i.e., 2014, 2016, and 2018) in cases of non-polio acute flaccid paralysis have led to speculations that other enteroviruses, specifically enterovirus D68 (EV-D68), are emerging to fill the niche that was left from poliovirus eradication. A cluster of 11 suspect cases of pediatric acute flaccid myelitis (AFM) was identified in 2016 in Phoenix, AZ. Multiple genomic analyses identified the presence of EV-D68 in the majority of clinical AFM cases. Beyond limited detection of herpesvirus, no other likely etiologies were found in the cluster. These findings strengthen the likelihood that EV-D68 is a cause of AFM and show that the rapid molecular assays developed for this study are useful for investigations of AFM and EV-D68.