Adeno-associated virus type 2 in US children with acute severe hepatitis.

As of August 2022, clusters of acute severe hepatitis of unknown aetiology in children have been reported from 35 countries, including the USA1,2. Previous studies have found human adenoviruses (HAdVs) in the blood from patients in Europe and the USA3-7, although it is unclear whether this virus is causative. Here we used PCR testing, viral enrichment-based sequencing and agnostic metagenomic sequencing to analyse samples from 16 HAdV-positive cases from 1 October 2021 to 22 May 2022, in parallel with 113 controls. In blood from 14 cases, adeno-associated virus type 2 (AAV2) sequences were detected in 93% (13 of 14), compared to 4 (3.5%) of 113 controls (P < 0.001) and to 0 of 30 patients with hepatitis of defined aetiology (P < 0.001). In controls, HAdV type 41 was detected in blood from 9 (39.1%) of the 23 patients with acute gastroenteritis (without hepatitis), including 8 of 9 patients with positive stool HAdV testing, but co-infection with AAV2 was observed in only 3 (13.0%) of these 23 patients versus 93% of cases (P < 0.001). Co-infections by Epstein-Barr virus, human herpesvirus 6 and/or enterovirus A71 were also detected in 12 (85.7%) of 14 cases, with higher herpesvirus detection in cases versus controls (P < 0.001). Our findings suggest that the severity of the disease is related to co-infections involving AAV2 and one or more helper viruses.

Clinical Metagenomic Sequencing for Diagnosis of Meningitis and Encephalitis.

BACKGROUND: Metagenomic next-generation sequencing (NGS) of cerebrospinal fluid (CSF) has the potential to identify a broad range of pathogens in a single test. METHODS: In a 1-year, multicenter, prospective study, we investigated the usefulness of metagenomic NGS of CSF for the diagnosis of infectious meningitis and encephalitis in hospitalized patients. All positive tests for pathogens on metagenomic NGS were confirmed by orthogonal laboratory testing. Physician feedback was elicited by teleconferences with a clinical microbial sequencing board and by surveys. Clinical effect was evaluated by retrospective chart review. RESULTS: We enrolled 204 pediatric and adult patients at eight hospitals. Patients were severely ill: 48.5% had been admitted to the intensive care unit, and the 30-day mortality among all study patients was 11.3%. A total of 58 infections of the nervous system were diagnosed in 57 patients (27.9%). Among these 58 infections, metagenomic NGS identified 13 (22%) that were not identified by clinical testing at the source hospital. Among the remaining 45 infections (78%), metagenomic NGS made concurrent diagnoses in 19. Of the 26 infections not identified by metagenomic NGS, 11 were diagnosed by serologic testing only, 7 were diagnosed from tissue samples other than CSF, and 8 were negative on metagenomic NGS owing to low titers of pathogens in CSF. A total of 8 of 13 diagnoses made solely by metagenomic NGS had a likely clinical effect, with 7 of 13 guiding treatment. CONCLUSIONS: Routine microbiologic testing is often insufficient to detect all neuroinvasive pathogens. In this study, metagenomic NGS of CSF obtained from patients with meningitis or encephalitis improved diagnosis of neurologic infections and provided actionable information in some cases. (Funded by the National Institutes of Health and others; PDAID ClinicalTrials.gov number, NCT02910037.).

Laboratory validation of a clinical metagenomic sequencing assay for pathogen detection in cerebrospinal fluid.

Metagenomic next-generation sequencing (mNGS) for pan-pathogen detection has been successfully tested in proof-of-concept case studies in patients with acute illness of unknown etiology but to date has been largely confined to research settings. Here, we developed and validated a clinical mNGS assay for diagnosis of infectious causes of meningitis and encephalitis from cerebrospinal fluid (CSF) in a licensed microbiology laboratory. A customized bioinformatics pipeline, SURPI+, was developed to rapidly analyze mNGS data, generate an automated summary of detected pathogens, and provide a graphical user interface for evaluating and interpreting results. We established quality metrics, threshold values, and limits of detection of 0.2-313 genomic copies or colony forming units per milliliter for each representative organism type. Gross hemolysis and excess host nucleic acid reduced assay sensitivity; however, spiked phages used as internal controls were reliable indicators of sensitivity loss. Diagnostic test accuracy was evaluated by blinded mNGS testing of 95 patient samples, revealing 73% sensitivity and 99% specificity compared to original clinical test results, and 81% positive percent agreement and 99% negative percent agreement after discrepancy analysis. Subsequent mNGS challenge testing of 20 positive CSF samples prospectively collected from a cohort of pediatric patients hospitalized with meningitis, encephalitis, and/or myelitis showed 92% sensitivity and 96% specificity relative to conventional microbiological testing of CSF in identifying the causative pathogen. These results demonstrate the analytic performance of a laboratory-validated mNGS assay for pan-pathogen detection, to be used clinically for diagnosis of neurological infections from CSF.

Diagnosis of Fatal Human Case of St. Louis Encephalitis Virus Infection by Metagenomic Sequencing, California, 2016.

We used unbiased metagenomic next-generation sequencing to diagnose a fatal case of meningoencephalitis caused by St. Louis encephalitis virus in a patient from California in September 2016. This case is associated with the recent 2015-2016 reemergence of this virus in the southwestern United States.

The Structure-Function Linkage Database.

The Structure-Function Linkage Database (SFLD, http://sfld.rbvi.ucsf.edu/) is a manually curated classification resource describing structure-function relationships for functionally diverse enzyme superfamilies. Members of such superfamilies are diverse in their overall reactions yet share a common ancestor and some conserved active site features associated with conserved functional attributes such as a partial reaction. Thus, despite their different functions, members of these superfamilies 'look alike', making them easy to misannotate. To address this complexity and enable rational transfer of functional features to unknowns only for those members for which we have sufficient functional information, we subdivide superfamily members into subgroups using sequence information, and lastly into families, sets of enzymes known to catalyze the same reaction using the same mechanistic strategy. Browsing and searching options in the SFLD provide access to all of these levels. The SFLD offers manually curated as well as automatically classified superfamily sets, both accompanied by search and download options for all hierarchical levels. Additional information includes multiple sequence alignments, tab-separated files of functional and other attributes, and sequence similarity networks. The latter provide a new and intuitively powerful way to visualize functional trends mapped to the context of sequence similarity.

Longitudinal comparison of the developing gut virome in infants and their mothers.

The human gut virome and its early life development are poorly understood. Prior studies have captured single-point assessments with the evolution of the infant virome remaining largely unexplored. We performed viral metagenomic sequencing on stool samples collected longitudinally from a cohort of 53 infants from age 2 weeks to 3 years (80.7 billion reads), and from their mothers (9.8 billion reads) to examine and compare viromes. The asymptomatic infant virome consisted of bacteriophages, nonhuman dietary/environmental viruses, and human-host viruses, predominantly picornaviruses. In contrast, human-host viruses were largely absent from the maternal virome. Previously undescribed, sequence-divergent vertebrate viruses were detected in the maternal but not infant virome. As infants aged, the phage component evolved to resemble the maternal virome, but by age 3, the human-host component remained dissimilar from the maternal virome. Thus, early life virome development is determined predominantly by dietary, infectious, and environmental factors rather than direct maternal acquisition.

Rapid pathogen detection by metagenomic next-generation sequencing of infected body fluids.

We developed a metagenomic next-generation sequencing (mNGS) test using cell-free DNA from body fluids to identify pathogens. The performance of mNGS testing of 182 body fluids from 160 patients with acute illness was evaluated using two sequencing platforms in comparison to microbiological testing using culture, 16S bacterial PCR and/or 28S-internal transcribed ribosomal gene spacer (28S-ITS) fungal PCR. Test sensitivity and specificity of detection were 79 and 91% for bacteria and 91 and 89% for fungi, respectively, by Illumina sequencing; and 75 and 81% for bacteria and 91 and 100% for fungi, respectively, by nanopore sequencing. In a case series of 12 patients with culture/PCR-negative body fluids but for whom an infectious diagnosis was ultimately established, seven (58%) were mNGS positive. Real-time computational analysis enabled pathogen identification by nanopore sequencing in a median 50-min sequencing and 6-h sample-to-answer time. Rapid mNGS testing is a promising tool for diagnosis of unknown infections from body fluids.

Genetic variants of human organic anion transporter 4 demonstrate altered transport of endogenous substrates.

Apical reabsorption from the urine has been shown to be important for such processes as the maintenance of critical metabolites in the blood and the excretion of nephrotoxic compounds. The solute carrier (SLC) transporter OAT4 (SLC22A11) is expressed on the apical membrane of renal proximal tubule cells and is known to mediate the transport of a variety of xenobiotic and endogenous organic anions. Functional characterization of genetic variants of apical transporters thought to mediate reabsorption, such as OAT4, may provide insight into the genetic factors influencing the complex pathways involved in drug elimination and metabolite reclamation occurring in the kidney. Naturally occurring genetic variants of OAT4 were identified in public databases and by resequencing DNA samples from 272 individuals comprising 4 distinct ethnic groups. Nine total nonsynonymous variants were identified and functionally assessed using uptake of three radiolabeled substrates. A nonsense variant, R48Stop, and three other variants (R121C, V155G, and V155M) were found at frequencies of at least 2% in an ethnic group specific fashion. The L29P, R48Stop, and H469R variants displayed a complete loss of function, and kinetic analysis identified a reduced V(max) in the common nonsynonymous variants. Plasma membrane levels of OAT4 protein were absent or reduced in the nonfunctional variants, providing a mechanistic reason for the observed loss of function. Characterization of the genetic variants of reabsorptive transporters such as OAT4 is an important step in understanding variability in tubular reabsorption with important implications in innate homeostatic processes and drug disposition.

Role of organic cation transporter 3 (SLC22A3) and its missense variants in the pharmacologic action of metformin.

OBJECTIVES: The goals of this study were to determine the role of organic cation transporter 3 (OCT3) in the pharmacological action of metformin and to identify and functionally characterize genetic variants of OCT3 with respect to the uptake of metformin and monoamines. METHODS: For pharmacological studies, we evaluated metformin-induced activation of AMP-activated protein kinase, a molecular target of metformin. We used quantitative PCR and immunostaining to localize the transporter and isotopic uptake studies in cells transfected with OCT3 and its nonsynonymous genetic variants for functional analyses. RESULTS: Quantitative PCR and immunostaining showed that OCT3 was expressed high on the plasma membrane of skeletal muscle and liver, target tissues for metformin action. Both the OCT inhibitor, cimetidine, and OCT3-specific short hairpin RNA significantly reduced the activating effect of metformin on AMP-activated protein kinase. To identify genetic variants in OCT3, we used recent data from the 1000 Genomes and the Pharmacogenomics of Membrane Transporters projects. Six novel missense variants were identified. In functional assays, using various monoamines and metformin, three variants, T44M (c.131C>T), T400I (c.1199C>T) and V423F (c.1267G>T) showed altered substrate specificity. Notably, in cells expressing T400I and V423F, the uptakes of metformin and catecholamines were significantly reduced, but the uptakes of metformin, 1-methyl-4-phenylpyridinium and histamine by T44M were significantly increased more than 50%. Structural modeling suggested that these two variants may be located in the pore lining (T400) or proximal (V423) membrane-spanning helixes. CONCLUSION: Our study suggests that OCT3 plays a role in the therapeutic action of metformin and that genetic variants of OCT3 may modulate metformin and catecholamine action.

Metagenomic sequencing of stool samples in Bangladeshi infants: virome association with poliovirus shedding after oral poliovirus vaccination.

The potential role of enteric viral infections and the developing infant virome in affecting immune responses to the oral poliovirus vaccine (OPV) is unknown. Here we performed viral metagenomic sequencing on 3 serially collected stool samples from 30 Bangladeshi infants following OPV vaccination and compared findings to stool samples from 16 age-matched infants in the United States (US). In 14 Bangladeshi infants, available post-vaccination serum samples were tested for polio-neutralizing antibodies. The abundance (p = 0.006) and richness (p = 0.013) of the eukaryotic virome increased with age and were higher than seen in age-matched US infants (p < 0.001). In contrast, phage diversity metrics remained stable and were similar to those in US infants. Non-poliovirus eukaryotic virus abundance (3.68 log10 vs. 2.25 log10, p = 0.002), particularly from potential viral pathogens (2.78log10 vs. 0.83log10, p = 0.002), and richness (p = 0.016) were inversely associated with poliovirus shedding. Following vaccination, 28.6% of 14 infants tested developed neutralizing antibodies to all three Sabin types and also exhibited higher rates of poliovirus shedding (p = 0.020). No vaccine-derived poliovirus variants were detected. These results reveal an inverse association between eukaryotic virome abundance and poliovirus shedding. Overall gut virome ecology and concurrent viral infections may impact oral vaccine responsiveness in Bangladeshi infants.

Functional genetic variation in the basal promoter of the organic cation/carnitine transporters OCTN1 (SLC22A4) and OCTN2 (SLC22A5).

The organic cation/ergothioneine transporter OCTN1 (SLC22A4) and the high-affinity carnitine transporter OCTN2 (SLC22A5), play an important role in the disposition of xenobiotics and endogenous compounds. Here, we analyzed the sequence of the proximal promoter regions of OCTN1 and OCTN2 in four ethnic groups and determined the effects of the identified genetic variants on transcriptional activities and mRNA expression. Six variants were found in the proximal promoter of OCTN1, one of which showed high allele frequency ranging from 13 to 34% in samples from individuals with ancestries in Africa, Europe, China, and Mexico. OCTN1 haplotypes had similar activities as the reference in luciferase reporter assays. For OCTN2, three of the seven variants identified in the proximal promoter showed allele frequencies greater than 29.5% in all populations, with the exception of -207C>G (rs2631367) that was monomorphic in Asian Americans. OCTN2 haplotypes containing -207G, present in all populations, were associated with a gain of function in luciferase reporter assays. Consistent with reporter assays, OCTN2 mRNA expression levels in lymphoblastoid cell lines (LCLs) from gene expression analysis were greater in samples carrying a marker for -207G. This SNP seems to contribute to racial differences in OCTN2 mRNA expression levels in LCLs. Our study with healthy subjects (n = 16) homozygous for either -207C or -207G, showed no appreciable effect of this SNP on carnitine disposition. However, there were significant effects of gender on carnitine plasma levels (p < 0.01). Further in vivo studies of OCTN2 promoter variants on carnitine disposition and variation in drug response are warranted.

Identification and characterization of proximal promoter polymorphisms in the human concentrative nucleoside transporter 2 (SLC28A2).

The human concentrative nucleoside transporter 2 (CNT2) plays an important role in the absorption, disposition, and biological effects of endogenous nucleosides and nucleoside analog drugs. We identified genetic variation in the basal promoter region of CNT2 and characterized the function of the variants. We screened DNA from an ethnically diverse population and identified five basal promoter variants in CNT2. Three major haplotypes in the CNT2 basal promoter region were identified and were found at different allele frequencies in various ethnic groups. The common promoter variants and haplotypes were constructed and characterized for their promoter activity using luciferase reporter assays. One polymorphic variant, rs2413775 (-146T>A), with an allele frequency >20% in all populations, showed a gain of function in luciferase activity. Furthermore, in vivo mouse promoter assays of these nucleotide variants using the hydrodynamic tail vein injection, leading to their expression in the liver, demonstrated similar results. Transcription factor binding site (TFBS) analysis indicated this variant alters a hepatic nuclear factor (HNF) 1 TFBS. Electrophoretic mobility shift assay demonstrated stronger binding of HNF1alpha and weaker binding of HNF1beta to the -146T and -146A regions, whereas the single nucleotide polymorphism (SNP), -146A, exhibited enhanced binding to both HNF1alpha and HNF1beta, consistent with its greater activity in reporter assays. The data collectively suggest that the common variant, -146T>A, in the proximal promoter of CNT2 may result in an enhanced transcription rate of the gene and, thus, expression levels of CNT2. This SNP may play a role in variation in the pharmacokinetics and pharmacological effects of nucleoside analogs.

Genetic variants in multidrug and toxic compound extrusion-1, hMATE1, alter transport function.

hMATE1 (human multidrug and toxin compound extrusion-1; encoded by SLC47A1) is thought to have an important function in the renal and hepatic elimination of drugs, endogenous compounds and environmental toxins. The goals of this study were to identify genetic variants of hMATE1 and to determine their effects on hMATE1 transport function. We identified four synonymous and six nonsynonymous, coding region variants in DNA samples from 272 individuals (68 Caucasians, 68 African Americans, 68 Asian Americans and 68 Mexican Americans). The overall prevalence of hMATE1 nonsynonymous variants was relatively low with three singleton variants and three variants having allele frequencies > or =2% in a specific ethnic group. The nonsynonymous hMATE1 variants were constructed and stably transfected into HEK-293 cells. Uptake studies using four known hMATE1 substrates (paraquat, metformin, tetraethylammonium and oxaliplatin) were performed in cells transfected with hMATE1 reference or variants. We found that two singleton variants, G64D and V480M, produced a complete loss of function for all four tested substrates whereas three polymorphic variants (allele frequencies > or =2%), L125F, V338I and C497S, significantly altered the transport function in a substrate-dependent manner. Confocal microscopy studies were consistent with functional studies suggesting that the altered function of the variants was due to altered localization to the plasma membrane. These data suggest that nonsynonymous variants in hMATE1 may alter drug disposition and ultimately affect clinical drug response.

The human multidrug resistance protein 4 (MRP4, ABCC4): functional analysis of a highly polymorphic gene.

ABCC4 encodes multidrug resistance protein 4 (MRP4), a member of the ATP-binding cassette family of membrane transporters involved in the efflux of endogenous and xenobiotic molecules. The aims of this study were to identify single nucleotide polymorphisms of ABCC4 and to functionally characterize selected nonsynonymous variants. Resequencing was performed in a large ethnically diverse population. Ten nonsynonymous variants were selected for analysis of transport function based on allele frequencies and evolutionary conservation. The reference and variant MRP4 cDNAs were constructed by site-directed mutagenesis and transiently transfected into human embryonic kidney cells (HEK 293T). The function of MRP4 variants was compared by measuring the intracellular accumulation of two antiviral agents, azidothymidine (AZT) and adefovir (PMEA). A total of 98 variants were identified in the coding and flanking intronic regions of ABCC4. Of these, 43 variants are in the coding region, and 22 are nonsynonymous. In a functional screen of ten variants, there was no evidence for a complete loss of function allele. However, two variants (G187W and G487E) showed a significantly reduced function compared to reference with both substrates, as evidenced by higher intracellular accumulation of AZT and PMEA compared to the reference MRP4 (43 and 69% increase in accumulation for G187W compared with the reference MRP4, with AZT and PMEA, respectively). The G187W variant also showed decreased expression following transient transfection of HEK 293T cells. Further studies are required to assess the clinical significance of this altered function and expression and to evaluate substrate specificity of this functional change.

The International Gene Trap Consortium Website: a portal to all publicly available gene trap cell lines in mouse.

Gene trapping is a method of generating murine embryonic stem (ES) cell lines containing insertional mutations in known and novel genes. A number of international groups have used this approach to create sizeable public cell line repositories available to the scientific community for the generation of mutant mouse strains. The major gene trapping groups worldwide have recently joined together to centralize access to all publicly available gene trap lines by developing a user-oriented Website for the International Gene Trap Consortium (IGTC). This collaboration provides an impressive public informatics resource comprising approximately 45 000 well-characterized ES cell lines which currently represent approximately 40% of known mouse genes, all freely available for the creation of knockout mice on a non-collaborative basis. To standardize annotation and provide high confidence data for gene trap lines, a rigorous identification and annotation pipeline has been developed combining genomic localization and transcript alignment of gene trap sequence tags to identify trapped loci. This information is stored in a new bioinformatics database accessible through the IGTC Website interface. The IGTC Website (www.genetrap.org) allows users to browse and search the database for trapped genes, BLAST sequences against gene trap sequence tags, and view trapped genes within biological pathways. In addition, IGTC data have been integrated into major genome browsers and bioinformatics sites to provide users with outside portals for viewing this data. The development of the IGTC Website marks a major advance by providing the research community with the data and tools necessary to effectively use public gene trap resources for the large-scale characterization of mammalian gene function.

Nanopore DNA Sequencing and Genome Assembly on the International Space Station.

We evaluated the performance of the MinION DNA sequencer in-flight on the International Space Station (ISS), and benchmarked its performance off-Earth against the MinION, Illumina MiSeq, and PacBio RS II sequencing platforms in terrestrial laboratories. Samples contained equimolar mixtures of genomic DNA from lambda bacteriophage, Escherichia coli (strain K12, MG1655) and Mus musculus (female BALB/c mouse). Nine sequencing runs were performed aboard the ISS over a 6-month period, yielding a total of 276,882 reads with no apparent decrease in performance over time. From sequence data collected aboard the ISS, we constructed directed assemblies of the ~4.6 Mb E. coli genome, ~48.5 kb lambda genome, and a representative M. musculus sequence (the ~16.3 kb mitochondrial genome), at 100%, 100%, and 96.7% consensus pairwise identity, respectively; de novo assembly of the E. coli genome from raw reads yielded a single contig comprising 99.9% of the genome at 98.6% consensus pairwise identity. Simulated real-time analyses of in-flight sequence data using an automated bioinformatic pipeline and laptop-based genomic assembly demonstrated the feasibility of sequencing analysis and microbial identification aboard the ISS. These findings illustrate the potential for sequencing applications including disease diagnosis, environmental monitoring, and elucidating the molecular basis for how organisms respond to spaceflight.

Comparison of methods for genomic localization of gene trap sequences.

BACKGROUND: Gene knockouts in a model organism such as mouse provide a valuable resource for the study of basic biology and human disease. Determining which gene has been inactivated by an untargeted gene trapping event poses a challenging annotation problem because gene trap sequence tags, which represent sequence near the vector insertion site of a trapped gene, are typically short and often contain unresolved residues. To understand better the localization of these sequences on the mouse genome, we compared stand-alone versions of the alignment programs BLAT, SSAHA, and MegaBLAST. A set of 3,369 sequence tags was aligned to build 34 of the mouse genome using default parameters for each algorithm. Known genome coordinates for the cognate set of full-length genes (1,659 sequences) were used to evaluate localization results. RESULTS: In general, all three programs performed well in terms of localizing sequences to a general region of the genome, with only relatively subtle errors identified for a small proportion of the sequence tags. However, large differences in performance were noted with regard to correctly identifying exon boundaries. BLAT correctly identified the vast majority of exon boundaries, while SSAHA and MegaBLAST missed the majority of exon boundaries. SSAHA consistently reported the fewest false positives and is the fastest algorithm. MegaBLAST was comparable to BLAT in speed, but was the most susceptible to localizing sequence tags incorrectly to pseudogenes. CONCLUSION: The differences in performance for sequence tags and full-length reference sequences were surprisingly small. Characteristic variations in localization results for each program were noted that affect the localization of sequence at exon boundaries, in particular.

BayGenomics: a resource of insertional mutations in mouse embryonic stem cells.

The BayGenomics gene-trap resource (http://baygenomics.ucsf.edu) provides researchers with access to thousands of mouse embryonic stem (ES) cell lines harboring characterized insertional mutations in both known and novel genes. Each cell line contains an insertional mutation in a specific gene. The identity of the gene that has been interrupted can be determined from a DNA sequence tag. Approximately 75% of our cell lines contain insertional mutations in known mouse genes or genes that share strong sequence similarities with genes that have been identified in other organisms. These cell lines readily transmit the mutation to the germline of mice and many mutant lines of mice have already been generated from this resource. BayGenomics provides facile access to our entire database, including sequence tags for each mutant ES cell line, through the World Wide Web. Investigators can browse our resource, search for specific entries, download any portion of our database and BLAST sequences of interest against our entire set of cell line sequence tags. They can then obtain the mutant ES cell line for the purpose of generating knockout mice.

Rapid metagenomic identification of viral pathogens in clinical samples by real-time nanopore sequencing analysis.

We report unbiased metagenomic detection of chikungunya virus (CHIKV), Ebola virus (EBOV), and hepatitis C virus (HCV) from four human blood samples by MinION nanopore sequencing coupled to a newly developed, web-based pipeline for real-time bioinformatics analysis on a computational server or laptop (MetaPORE). At titers ranging from 10(7)-10(8) copies per milliliter, reads to EBOV from two patients with acute hemorrhagic fever and CHIKV from an asymptomatic blood donor were detected within 4 to 10 min of data acquisition, while lower titer HCV virus (1 × 10(5) copies per milliliter) was detected within 40 min. Analysis of mapped nanopore reads alone, despite an average individual error rate of 24 % (range 8-49 %), permitted identification of the correct viral strain in all four isolates, and 90 % of the genome of CHIKV was recovered with 97-99 % accuracy. Using nanopore sequencing, metagenomic detection of viral pathogens directly from clinical samples was performed within an unprecedented <6 hr sample-to-answer turnaround time, and in a timeframe amenable to actionable clinical and public health diagnostics.

A novel outbreak enterovirus D68 strain associated with acute flaccid myelitis cases in the USA (2012-14): a retrospective cohort study.

BACKGROUND: Enterovirus D68 was implicated in a widespread outbreak of severe respiratory illness across the USA in 2014 and has also been reported sporadically in patients with acute flaccid myelitis. We aimed to investigate the association between enterovirus D68 infection and acute flaccid myelitis during the 2014 enterovirus D68 respiratory outbreak in the USA. METHODS: Patients with acute flaccid myelitis who presented to two hospitals in Colorado and California, USA, between Nov 24, 2013, and Oct 11, 2014, were included in the study. Additional cases identified from Jan 1, 2012, to Oct 4, 2014, via statewide surveillance were provided by the California Department of Public Health. We investigated the cause of these cases by metagenomic next-generation sequencing, viral genome recovery, and enterovirus D68 phylogenetic analysis. We compared patients with acute flaccid myelitis who were positive for enterovirus D68 with those with acute flaccid myelitis but negative for enterovirus D68 using the two-tailed Fisher's exact test, two-sample unpaired t test, and Mann-Whitney U test. FINDINGS: 48 patients were included: 25 with acute flaccid myelitis, two with enterovirus-associated encephalitis, five with enterovirus-D68-associated upper respiratory illness, and 16 with aseptic meningitis or encephalitis who tested positive for enterovirus. Enterovirus D68 was detected in respiratory secretions from seven (64%) of 11 patients comprising two temporally and geographically linked acute flaccid myelitis clusters at the height of the 2014 outbreak, and from 12 (48%) of 25 patients with acute flaccid myelitis overall. Phylogenetic analysis revealed that all enterovirus D68 sequences associated with acute flaccid myelitis grouped into a clade B1 strain that emerged in 2010. Of six coding polymorphisms in the clade B1 enterovirus D68 polyprotein, five were present in neuropathogenic poliovirus or enterovirus D70, or both. One child with acute flaccid myelitis and a sibling with only upper respiratory illness were both infected by identical enterovirus D68 strains. Enterovirus D68 viraemia was identified in a child experiencing acute neurological progression of his paralytic illness. Deep metagenomic sequencing of cerebrospinal fluid from 14 patients with acute flaccid myelitis did not reveal evidence of an alternative infectious cause to enterovirus D68. INTERPRETATION: These findings strengthen the putative association between enterovirus D68 and acute flaccid myelitis and the contention that acute flaccid myelitis is a rare yet severe clinical manifestation of enterovirus D68 infection in susceptible hosts. FUNDING: National Institutes of Health, University of California, Abbott Laboratories, and the Centers for Disease Control and Prevention.