Search

Antimicrobial resistant lower respiratory tract infections are an increasing public health threat and an important cause of global mortality. The lung microbiome can influence susceptibility of respiratory tract infections and represents an important reservoir for exchange of antimicrobial resistance genes. Studies of the gut microbiome have found an association between age and increasing antimicrobial resistance gene burden, however, corollary studies in the lung microbiome remain absent. We performed an observational study of children and adults with acute respiratory failure admitted to the intensive care unit. From tracheal aspirate RNA sequencing data, we evaluated age-related differences in detectable antimicrobial resistance gene expression in the lung microbiome. Using a multivariable logistic regression model, we find that detection of antimicrobial resistance gene expression was significantly higher in adults compared with children after adjusting for demographic and clinical characteristics. This association remained significant after additionally adjusting for lung bacterial microbiome characteristics, and when modeling age as a continuous variable. The proportion of adults expressing beta-lactam, aminoglycoside, and tetracycline antimicrobial resistance genes was higher compared to children. Together, these findings shape our understanding of the lung resistome in critically ill patients across the lifespan, which may have implications for clinical management and global public health.

Two independent temporal-spatial clusters of hospital-onset Rhizopus infections were evaluated using whole-genome sequencing (WGS). Phylogenetic analysis confirmed that isolates within each cluster were unrelated despite epidemiological suspicion of outbreaks. The ITS1 region alone was insufficient for accurate analysis. WGS has utility for rapid rule-out of suspected nosocomial Rhizopus outbreaks.

BACKGROUNDLower respiratory tract infection (LRTI) is a leading cause of death in children worldwide. LRTI diagnosis is challenging because noninfectious respiratory illnesses appear clinically similar and because existing microbiologic tests are often falsely negative or detect incidentally carried microbes, resulting in antimicrobial overuse and adverse outcomes. Lower airway metagenomics has the potential to detect host and microbial signatures of LRTI. Whether it can be applied at scale and in a pediatric population to enable improved diagnosis and treatment remains unclear.METHODSWe used tracheal aspirate RNA-Seq to profile host gene expression and respiratory microbiota in 261 children with acute respiratory failure. We developed a gene expression classifier for LRTI by training on patients with an established diagnosis of LRTI (n = 117) or of noninfectious respiratory failure (n = 50). We then developed a classifier that integrates the host LRTI probability, abundance of respiratory viruses, and dominance in the lung microbiome of bacteria/fungi considered pathogenic by a rules-based algorithm.RESULTSThe host classifier achieved a median AUC of 0.967 by cross-validation, driven by activation markers of T cells, alveolar macrophages, and the interferon response. The integrated classifier achieved a median AUC of 0.986 and increased the confidence of patient classifications. When applied to patients with an uncertain diagnosis (n = 94), the integrated classifier indicated LRTI in 52% of cases and nominated likely causal pathogens in 98% of those.CONCLUSIONLower airway metagenomics enables accurate LRTI diagnosis and pathogen identification in a heterogeneous cohort of critically ill children through integration of host, pathogen, and microbiome features.FUNDINGSupport for this study was provided by the Eunice Kennedy Shriver National Institute of Child Health and Human Development and the National Heart, Lung, and Blood Institute (UG1HD083171, 1R01HL124103, UG1HD049983, UG01HD049934, UG1HD083170, UG1HD050096, UG1HD63108, UG1HD083116, UG1HD083166, UG1HD049981, K23HL138461, and 5R01HL155418) as well as by the Chan Zuckerberg Biohub.

The aetiology of nodding syndrome remains unclear, and comprehensive genotyping and phenotyping data from patients remain sparse. Our objectives were to characterize the phenotype of patients with nodding syndrome, investigate potential contributors to disease aetiology, and evaluate response to immunotherapy. This cohort study investigated members of a single-family unit from Lamwo District, Uganda. The participants for this study were selected by the Ugandan Ministry of Health as representative for nodding syndrome and with a conducive family structure for genomic analyses. Of the eight family members who participated in the study at the National Institutes of Health (NIH) Clinical Center, three had nodding syndrome. The three affected patients were extensively evaluated with metagenomic sequencing for infectious pathogens, exome sequencing, spinal fluid immune analyses, neurometabolic and toxicology testing, continuous electroencephalography and neuroimaging. Five unaffected family members underwent a subset of testing for comparison. A distinctive interictal pattern of sleep-activated bursts of generalized and multifocal epileptiform discharges and slowing was observed in two patients. Brain imaging showed two patients had mild generalized cerebral atrophy, and both patients and unaffected family members had excessive metal deposition in the basal ganglia. Trace metal biochemical evaluation was normal. CSF was non-inflammatory and one patient had CSF-restricted oligoclonal bands. Onchocerca volvulus-specific antibodies were present in all patients and skin snips were negative for active onchocerciasis. Metagenomic sequencing of serum and CSF revealed hepatitis B virus in the serum of one patient. Vitamin B6 metabolites were borderline low in all family members and CSF pyridoxine metabolites were normal. Mitochondrial DNA testing was normal. Exome sequencing did not identify potentially causal candidate gene variants. Nodding syndrome is characterized by a distinctive pattern of sleep-activated epileptiform activity. The associated growth stunting may be due to hypothalamic dysfunction. Extensive testing years after disease onset did not clarify a causal aetiology. A trial of immunomodulation (plasmapheresis in two patients and intravenous immunoglobulin in one patient) was given without short-term effect, but longer-term follow-up was not possible to fully assess any benefit of this intervention.

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false-negative viral PCR test results. Such tests are also susceptible to false-positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses, and nonviral conditions (n = 318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to diagnose COVID-19. We find that optimal classifiers include an interferon-stimulated gene that is strongly induced in COVID-19 compared with nonviral conditions, such as IFI6, and a second immune-response gene that is more strongly induced in other viral infections, such as GBP5. The IFI6+GBP5 classifier achieves an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n = 553). We further provide proof-of-concept demonstration that the classifier can be implemented in a clinically relevant RT-qPCR assay. Finally, we show that its performance is robust across common SARS-CoV-2 variants and is unaffected by cross-contamination, demonstrating its utility for improved accuracy of COVID-19 diagnostics. IMPORTANCE In this work, we study upper respiratory tract gene expression to develop and validate a 2-gene host-based COVID-19 diagnostic classifier and then demonstrate its implementation in a clinically practical qPCR assay. We find that the host classifier has utility for mitigating false-negative results, for example due to SARS-CoV-2 variants harboring mutations at primer target sites, and for mitigating false-positive viral PCR results due to laboratory cross-contamination. Both types of error carry serious consequences of either unrecognized viral transmission or unnecessary isolation and contact tracing. This work is directly relevant to the ongoing COVID-19 pandemic given the continued emergence of viral variants and the continued challenges of false-positive PCR assays. It also suggests the feasibility of pan-respiratory virus host-based diagnostics that would have value in congregate settings, such as hospitals and nursing homes, where unrecognized respiratory viral transmission is of particular concern.

Multidrug-resistant (MDR) Pseudomonas aeruginosa has been declared a serious threat by the United States Centers for Disease Control and Prevention. Here, we used whole genome sequencing (WGS) to investigate recurrent P. aeruginosa bloodstream infections in a severely immunocompromised patient. The infections demonstrated unusual, progressive increases in resistance to beta lactam antibiotics in the setting of active treatment with appropriate, guideline-directed agents. WGS followed by comparative genomic analysis of isolates collected over 44 days demonstrated in host evolution of a single P. aeruginosa isolate characterized by stepwise acquisition of two de-novo genetic resistance mechanisms over the course of treatment. We found a novel deletion affecting the ampC repressor ampD and neighboring gene ampE, which associated with initial cefepime treatment failure. This was followed by acquisition of a porin nonsense mutation, OprD, associated with resistance to carbapenems. This study highlights the potential for in-host evolution of P. aeruginosa during bloodstream infections in severely immunocompromised patients despite appropriate antimicrobial therapy. In addition, it demonstrates the utility of WGS for understanding unusual resistance patterns in the clinical context.

BackgroundDuring the COVID-19 pandemic within the United States, much of the responsibility for diagnostic testing and epidemiologic response has relied on the action of county-level departments of public health. Here we describe the integration of genomic surveillance into epidemiologic response within Humboldt County, a rural county in northwest California.MethodsThrough a collaborative effort, 853 whole SARS-CoV-2 genomes were generated, representing ~58% of the 1,449 SARS-CoV-2-positive cases detected in Humboldt County as of March 12, 2021. Phylogenetic analysis of these data was used to develop a comprehensive understanding of SARS-CoV-2 introductions to the county and to support contact tracing and epidemiologic investigations of all large outbreaks in the county.ResultsIn the case of an outbreak on a commercial farm, viral genomic data were used to validate reported epidemiologic links and link additional cases within the community who did not report a farm exposure to the outbreak. During a separate outbreak within a skilled nursing facility, genomic surveillance data were used to rule out the putative index case, detect the emergence of an independent Spike:N501Y substitution, and verify that the outbreak had been brought under control.ConclusionsThese use cases demonstrate how developing genomic surveillance capacity within local public health departments can support timely and responsive deployment of genomic epidemiology for surveillance and outbreak response based on local needs and priorities.

We carried out integrated host and pathogen metagenomic RNA and DNA next generation sequencing (mNGS) of whole blood (n = 221) and plasma (n = 138) from critically ill patients following hospital admission. We assigned patients into sepsis groups on the basis of clinical and microbiological criteria. From whole-blood gene expression data, we distinguished patients with sepsis from patients with non-infectious systemic inflammatory conditions using a trained bagged support vector machine (bSVM) classifier (area under the receiver operating characteristic curve (AUC) = 0.81 in the training set; AUC = 0.82 in a held-out validation set). Plasma RNA also yielded a transcriptional signature of sepsis with several genes previously reported as sepsis biomarkers, and a bSVM sepsis diagnostic classifier (AUC = 0.97 training set; AUC = 0.77 validation set). Pathogen detection performance of plasma mNGS varied on the basis of pathogen and site of infection. To improve detection of virus, we developed a secondary transcriptomic classifier (AUC = 0.94 training set; AUC = 0.96 validation set). We combined host and microbial features to develop an integrated sepsis diagnostic model that identified 99% of microbiologically confirmed sepsis cases, and predicted sepsis in 74% of suspected and 89% of indeterminate sepsis cases. In summary, we suggest that integrating host transcriptional profiling and broad-range metagenomic pathogen detection from nucleic acid is a promising tool for sepsis diagnosis.

The genus Henipavirus (family Paramyxoviridae) currently comprises seven viruses, four of which have demonstrated prior evidence of zoonotic capacity. These include the biosafety level 4 agents Hendra (HeV) and Nipah (NiV) viruses, which circulate naturally in pteropodid fruit bats. Here, we describe and characterize Angavokely virus (AngV), a divergent henipavirus identified in urine samples from wild, Madagascar fruit bats. We report the nearly complete 16,740-nucleotide genome of AngV, which encodes the six major henipavirus structural proteins (nucleocapsid, phosphoprotein, matrix, fusion, glycoprotein, and L polymerase). Within the phosphoprotein (P) gene, we identify an alternative start codon encoding the AngV C protein and a putative mRNA editing site where the insertion of one or two guanine residues encodes, respectively, additional V and W proteins. In other paramyxovirus systems, C, V, and W are accessory proteins involved in antagonism of host immune responses during infection. Phylogenetic analysis suggests that AngV is ancestral to all four previously described bat henipaviruses-HeV, NiV, Cedar virus (CedV), and Ghanaian bat virus (GhV)-but evolved more recently than rodent- and shrew-derived henipaviruses, Mojiang (MojV), Gamak (GAKV), and Daeryong (DARV) viruses. Predictive structure-based alignments suggest that AngV is unlikely to bind ephrin receptors, which mediate cell entry for all other known bat henipaviruses. Identification of the AngV receptor is needed to clarify the virus's potential host range. The presence of V and W proteins in the AngV genome suggest that the virus could be pathogenic following zoonotic spillover. IMPORTANCE Henipaviruses include highly pathogenic emerging zoonotic viruses, derived from bat, rodent, and shrew reservoirs. Bat-borne Hendra (HeV) and Nipah (NiV) are the most well-known henipaviruses, for which no effective antivirals or vaccines for humans have been described. Here, we report the discovery and characterization of a novel henipavirus, Angavokely virus (AngV), isolated from wild fruit bats in Madagascar. Genomic characterization of AngV reveals all major features associated with pathogenicity in other henipaviruses, suggesting that AngV could be pathogenic following spillover to human hosts. Our work suggests that AngV is an ancestral bat henipavirus that likely uses viral entry pathways distinct from those previously described for HeV and NiV. In Madagascar, bats are consumed as a source of human food, presenting opportunities for cross-species transmission. Characterization of novel henipaviruses and documentation of their pathogenic and zoonotic potential are essential to predicting and preventing the emergence of future zoonoses that cause pandemics.

BackgroundAntimicrobial resistance (AMR) is rising at an alarming rate and complicating the management of infectious diseases including lower respiratory tract infections (LRTI). Metagenomic next-generation sequencing (mNGS) is a recently established method for culture-independent LRTI diagnosis, but its utility for predicting AMR has remained unclear. We aimed to assess the performance of mNGS for AMR prediction in bacterial LRTI and demonstrate proof of concept for epidemiological AMR surveillance and rapid AMR gene detection using Cas9 enrichment and nanopore sequencing.MethodsWe studied 88 patients with acute respiratory failure between 07/2013 and 9/2018, enrolled through a previous observational study of LRTI. Inclusion criteria were age ≥ 18, need for mechanical ventilation, and respiratory specimen collection within 72 h of intubation. Exclusion criteria were decline of study participation, unclear LRTI status, or no matched RNA and DNA mNGS data from a respiratory specimen. Patients with LRTI were identified by clinical adjudication. mNGS was performed on lower respiratory tract specimens. The primary outcome was mNGS performance for predicting phenotypic antimicrobial susceptibility and was assessed in patients with LRTI from culture-confirmed bacterial pathogens with clinical antimicrobial susceptibility testing (n = 27 patients, n = 32 pathogens). Secondary outcomes included the association between hospital exposure and AMR gene burden in the respiratory microbiome (n = 88 patients), and AMR gene detection using Cas9 targeted enrichment and nanopore sequencing (n = 10 patients).ResultsCompared to clinical antimicrobial susceptibility testing, the performance of respiratory mNGS for predicting AMR varied by pathogen, antimicrobial, and nucleic acid type sequenced. For gram-positive bacteria, a combination of RNA + DNA mNGS achieved a sensitivity of 70% (95% confidence interval (CI) 47-87%) and specificity of 95% (CI 85-99%). For gram-negative bacteria, sensitivity was 100% (CI 87-100%) and specificity 64% (CI 48-78%). Patients with hospital-onset LRTI had a greater AMR gene burden in their respiratory microbiome versus those with community-onset LRTI (p = 0.00030), or those without LRTI (p = 0.0024). We found that Cas9 targeted sequencing could enrich for low abundance AMR genes by > 2500-fold and enabled their rapid detection using a nanopore platform.ConclusionsmNGS has utility for the detection and surveillance of resistant bacterial LRTI pathogens.

SignificanceMetagenomic pathogen sequencing offers an unbiased approach to characterizing febrile illness. In resource-scarce settings with high biodiversity, it is critical to identify disease-causing pathogens in order to understand burden and to prioritize efforts for control. Here, metagenomic next-generation sequencing (mNGS) characterization of the pathogen landscape in Cambodia revealed diverse vector-borne and zoonotic pathogens irrespective of age and gender as risk factors. Identification of key pathogens led to changes in national program surveillance. This study is a "real world" example of the use of mNGS surveillance of febrile individuals, executed in-country, to identify outbreaks of vector-borne, zoonotic, and other emerging pathogens in a resource-scarce setting.

Multiple sclerosis (MS) is a heterogenous autoimmune disease in which autoreactive lymphocytes attack the myelin sheath of the central nervous system. B lymphocytes in the cerebrospinal fluid (CSF) of patients with MS contribute to inflammation and secrete oligoclonal immunoglobulins1,2. Epstein-Barr virus (EBV) infection has been epidemiologically linked to MS, but its pathological role remains unclear3. Here we demonstrate high-affinity molecular mimicry between the EBV transcription factor EBV nuclear antigen 1 (EBNA1) and the central nervous system protein glial cell adhesion molecule (GlialCAM) and provide structural and in vivo functional evidence for its relevance. A cross-reactive CSF-derived antibody was initially identified by single-cell sequencing of the paired-chain B cell repertoire of MS blood and CSF, followed by protein microarray-based testing of recombinantly expressed CSF-derived antibodies against MS-associated viruses. Sequence analysis, affinity measurements and the crystal structure of the EBNA1-peptide epitope in complex with the autoreactive Fab fragment enabled tracking of the development of the naive EBNA1-restricted antibody to a mature EBNA1-GlialCAM cross-reactive antibody. Molecular mimicry is facilitated by a post-translational modification of GlialCAM. EBNA1 immunization exacerbates disease in a mouse model of MS, and anti-EBNA1 and anti-GlialCAM antibodies are prevalent in patients with MS. Our results provide a mechanistic link for the association between MS and EBV and could guide the development of new MS therapies.

Impaired baseline lung function is associated with mortality after pediatric allogeneic hematopoietic cell transplantation (HCT), yet limited knowledge of the molecular pathways that characterize pretransplant lung function has hindered the development of lung-targeted interventions. In this study, we quantified the association between bronchoalveolar lavage (BAL) metatranscriptomes and paired pulmonary function tests performed a median of 1 to 2 weeks before allogeneic HCT in 104 children in The Netherlands. Abnormal pulmonary function was recorded in more than half the cohort, consisted most commonly of restriction and impaired diffusion, and was associated with both all-cause and lung injury-related mortality after HCT. Depletion of commensal supraglottic taxa, such as Haemophilus, and enrichment of nasal and skin taxa, such as Staphylococcus, in the BAL microbiome were associated with worse measures of lung capacity and gas diffusion. In addition, BAL gene expression signatures of alveolar epithelial activation, epithelial-mesenchymal transition, and down-regulated immunity were associated with impaired lung capacity and diffusion, suggesting a postinjury profibrotic response. Detection of microbial depletion and abnormal epithelial gene expression in BAL enhanced the prognostic utility of pre-HCT pulmonary function tests for the outcome of post-HCT mortality. These findings suggest a potentially actionable connection between microbiome depletion, alveolar injury, and pulmonary fibrosis in the pathogenesis of pre-HCT lung dysfunction.

BackgroundThe eukaryotic parasite Plasmodium falciparum causes millions of malarial infections annually while drug resistance to common anti-malarials is further confounding eradication efforts. Translation is an attractive therapeutic target that will benefit from a deeper mechanistic understanding. As the rate limiting step of translation, initiation is a primary driver of translational efficiency. It is a complex process regulated by both cis and trans acting factors, providing numerous potential targets. Relative to model organisms and humans, P. falciparum mRNAs feature unusual 5' untranslated regions suggesting cis-acting sequence complexity in this parasite may act to tune levels of protein synthesis through their effects on translational efficiency.MethodsHere, in vitro translation is deployed to compare the role of cis-acting regulatory sequences in P. falciparum and humans. Using parasite mRNAs with high or low translational efficiency, the presence, position, and termination status of upstream "AUG"s, in addition to the base composition of the 5' untranslated regions, were characterized.ResultsThe density of upstream "AUG"s differed significantly among the most and least efficiently translated genes in P. falciparum, as did the average "GC" content of the 5' untranslated regions. Using exemplars from highly translated and poorly translated mRNAs, multiple putative upstream elements were interrogated for impact on translational efficiency. Upstream "AUG"s were found to repress translation to varying degrees, depending on their position and context, while combinations of upstream "AUG"s had non-additive effects. The base composition of the 5' untranslated regions also impacted translation, but to a lesser degree. Surprisingly, the effects of cis-acting sequences were remarkably conserved between P. falciparum and humans.ConclusionsWhile translational regulation is inherently complex, this work contributes toward a more comprehensive understanding of parasite and human translational regulation by examining the impact of discrete cis-acting features, acting alone or in context.

Bats are natural reservoirs for both Alpha- and Betacoronaviruses and the hypothesized original hosts of five of seven known zoonotic coronaviruses. To date, the vast majority of bat coronavirus research has been concentrated in Asia, though coronaviruses are globally distributed; indeed, SARS-CoV and SARS-CoV-2-related Betacoronaviruses in the subgenus Sarbecovirus have been identified circulating in Rhinolophid bats in both Africa and Europe, despite the relative dearth of surveillance in these regions. As part of a long-term study examining the dynamics of potentially zoonotic viruses in three species of endemic Madagascar fruit bat (Pteropus rufus, Eidolon dupreanum, Rousettus madagascariensis), we carried out metagenomic Next Generation Sequencing (mNGS) on urine, throat, and fecal samples obtained from wild-caught individuals. We report detection of RNA derived from Betacoronavirus subgenus Nobecovirus in fecal samples from all three species and describe full genome sequences of novel Nobecoviruses in P. rufus and R. madagascariensis. Phylogenetic analysis indicates the existence of five distinct Nobecovirus clades, one of which is defined by the highly divergent ancestral sequence reported here from P. rufus bats. Madagascar Nobecoviruses derived from P. rufus and R. madagascariensis demonstrate, respectively, Asian and African phylogeographic origins, mirroring those of their fruit bat hosts. Bootscan recombination analysis indicates significant selection has taken place in the spike, nucleocapsid, and NS7 accessory protein regions of the genome for viruses derived from both bat hosts. Madagascar offers a unique phylogeographic nexus of bats and viruses with both Asian and African phylogeographic origins, providing opportunities for unprecedented mixing of viral groups and, potentially, recombination. As fruit bats are handled and consumed widely across Madagascar for subsistence, understanding the landscape of potentially zoonotic coronavirus circulation is essential for mitigation of future zoonotic threats.

BackgroundMetagenomic next-generation sequencing (mNGS) of body fluids is an emerging approach to identify occult pathogens in undiagnosed patients. We hypothesized that metagenomic testing can be simultaneously used to detect malignant neoplasms in addition to infectious pathogens.MethodsFrom two independent studies (n = 205), we used human data generated from a metagenomic sequencing pipeline to simultaneously screen for malignancies by copy number variation (CNV) detection. In the first case-control study, we analyzed body fluid samples (n = 124) from patients with a clinical diagnosis of either malignancy (positive cases, n = 65) or infection (negative controls, n = 59). In a second verification cohort, we analyzed a series of consecutive cases (n = 81) sent to cytology for malignancy workup that included malignant positives (n = 32), negatives (n = 18), or cases with an unclear gold standard (n = 31).ResultsThe overall CNV test sensitivity across all studies was 87% (55 of 63) in patients with malignancies confirmed by conventional cytology and/or flow cytometry testing and 68% (23 of 34) in patients who were ultimately diagnosed with cancer but negative by conventional testing. Specificity was 100% (95% CI 95-100%) with no false positives detected in 77 negative controls. In one example, a patient hospitalized with an unknown pulmonary illness had non-diagnostic lung biopsies, while CNVs implicating a malignancy were detectable from bronchoalveolar fluid.ConclusionsMetagenomic sequencing of body fluids can be used to identify undetected malignant neoplasms through copy number variation detection. This study illustrates the potential clinical utility of a single metagenomic test to uncover the cause of undiagnosed acute illnesses due to cancer or infection using the same specimen.

ImportanceCerebrospinal fluid (CSF) cytologic testing and flow cytometry are insensitive for diagnosing neoplasms of the central nervous system (CNS). Such clinical phenotypes can mimic infectious and autoimmune causes of meningoencephalitis.ObjectiveTo ascertain whether CSF metagenomic next-generation sequencing (mNGS) can identify aneuploidy, a hallmark of malignant neoplasms, in difficult-to-diagnose cases of CNS malignant neoplasm.Design, setting, and participantsTwo case-control studies were performed at the University of California, San Francisco (UCSF). The first study used CSF specimens collected at the UCSF Clinical Laboratories between July 1, 2017, and December 31, 2019, and evaluated test performance in specimens from patients with a CNS malignant neoplasm (positive controls) or without (negative controls). The results were compared with those from CSF cytologic testing and/or flow cytometry. The second study evaluated patients who were enrolled in an ongoing prospective study between April 1, 2014, and July 31, 2019, with presentations that were suggestive of neuroinflammatory disease but who were ultimately diagnosed with a CNS malignant neoplasm. Cases of individuals whose tumors could have been detected earlier without additional invasive testing are discussed.Main outcomes and measuresThe primary outcome measures were the sensitivity and specificity of aneuploidy detection by CSF mNGS. Secondary subset analyses included a comparison of CSF and tumor tissue chromosomal abnormalities and the identification of neuroimaging characteristics that were associated with test performance.ResultsAcross both studies, 130 participants were included (median [interquartile range] age, 57.5 [43.3-68.0] years; 72 men [55.4%]). The test performance study used 125 residual laboratory CSF specimens from 47 patients with a CNS malignant neoplasm and 56 patients with other neurological diseases. The neuroinflammatory disease study enrolled 12 patients and 17 matched control participants. The sensitivity of the CSF mNGS assay was 75% (95% CI, 63%-85%), and the specificity was 100% (95% CI, 96%-100%). Aneuploidy was detected in 64% (95% CI, 41%-83%) of the patients in the test performance study with nondiagnostic cytologic testing and/or flow cytometry, and in 55% (95% CI, 23%-83%) of patients in the neuroinflammatory disease study who were ultimately diagnosed with a CNS malignant neoplasm. Of the patients in whom aneuploidy was detected, 38 (90.5%) had multiple copy number variations with tumor fractions ranging from 31% to 49%.Conclusions and relevanceThis case-control study showed that CSF mNGS, which has low specimen volume requirements, does not require the preservation of cell integrity, and was orginally developed to diagnose neurologic infections, can also detect genetic evidence of a CNS malignant neoplasm in patients in whom CSF cytologic testing and/or flow cytometry yielded negative results with a low risk of false-positive results.

Neutralizing autoantibodies against type I interferons (IFNs) have been found in some patients with critical coronavirus disease 2019 (COVID-19), the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the prevalence of these antibodies, their longitudinal dynamics across the disease severity scale, and their functional effects on circulating leukocytes remain unknown. Here, in 284 patients with COVID-19, we found type I IFN–specific autoantibodies in peripheral blood samples from 19% of patients with critical disease and 6% of patients with severe disease. We found no type I IFN autoantibodies in individuals with moderate disease. Longitudinal profiling of over 600,000 peripheral blood mononuclear cells using multiplexed single-cell epitope and transcriptome sequencing from 54 patients with COVID-19 and 26 non–COVID-19 controls revealed a lack of type I IFN–stimulated gene (ISG-I) responses in myeloid cells from patients with critical disease. This was especially evident in dendritic cell populations isolated from patients with critical disease producing type I IFN–specific autoantibodies. Moreover, we found elevated expression of the inhibitory receptor leukocyte-associated immunoglobulin-like receptor 1 (LAIR1) on the surface of monocytes isolated from patients with critical disease early in the disease course. LAIR1 expression is inversely correlated with ISG-I expression response in patients with COVID-19 but is not expressed in healthy controls. The deficient ISG-I response observed in patients with critical COVID-19 with and without type I IFN–specific autoantibodies supports a unifying model for disease pathogenesis involving ISG-I suppression through convergent mechanisms.

The immunological features that distinguish COVID-19-associated acute respiratory distress syndrome (ARDS) from other causes of ARDS are incompletely understood. Here, we report the results of comparative lower respiratory tract transcriptional profiling of tracheal aspirate from 52 critically ill patients with ARDS from COVID-19 or from other etiologies, as well as controls without ARDS. In contrast to a "cytokine storm," we observe reduced proinflammatory gene expression in COVID-19 ARDS when compared to ARDS due to other causes. COVID-19 ARDS is characterized by a dysregulated host response with increased PTEN signaling and elevated expression of genes with non-canonical roles in inflammation and immunity. In silico analysis of gene expression identifies several candidate drugs that may modulate gene expression in COVID-19 ARDS, including dexamethasone and granulocyte colony stimulating factor. Compared to ARDS due to other types of viral pneumonia, COVID-19 is characterized by impaired interferon-stimulated gene (ISG) expression. The relationship between SARS-CoV-2 viral load and expression of ISGs is decoupled in patients with COVID-19 ARDS when compared to patients with mild COVID-19. In summary, assessment of host gene expression in the lower airways of patients reveals distinct immunological features of COVID-19 ARDS.

Narnaviruses are RNA viruses detected in diverse fungi, plants, protists, arthropods, and nematodes. Though initially described as simple single-gene nonsegmented viruses encoding RNA-dependent RNA polymerase (RdRp), a subset of narnaviruses referred to as "ambigrammatic" harbor a unique genomic configuration consisting of overlapping open reading frames (ORFs) encoded on opposite strands. Phylogenetic analysis supports selection to maintain this unusual genome organization, but functional investigations are lacking. Here, we establish the mosquito-infecting Culex narnavirus 1 (CxNV1) as a model to investigate the functional role of overlapping ORFs in narnavirus replication. In CxNV1, a reverse ORF without homology to known proteins covers nearly the entire 3.2-kb segment encoding the RdRp. Additionally, two opposing and nearly completely overlapping novel ORFs are found on the second putative CxNV1 segment, the 0.8-kb "Robin" RNA. We developed a system to launch CxNV1 in a naive mosquito cell line and then showed that functional RdRp is required for persistence of both segments, and an intact reverse ORF is required on the RdRp segment for persistence. Mass spectrometry of persistently CxNV1-infected cells provided evidence for translation of this reverse ORF. Finally, ribosome profiling yielded a striking pattern of footprints for all four CxNV1 RNA strands that was distinct from actively translating ribosomes on host mRNA or coinfecting RNA viruses. Taken together, these data raise the possibility that the process of translation itself is important for persistence of ambigrammatic narnaviruses, potentially by protecting viral RNA with ribosomes, thus suggesting a heretofore undescribed viral tactic for replication and transmission. IMPORTANCE Fundamental to our understanding of RNA viruses is a description of which strand(s) of RNA are transmitted as the viral genome relative to which encode the viral proteins. Ambigrammatic narnaviruses break the mold. These viruses, found broadly in fungi, plants, and insects, have the unique feature of two overlapping genes encoded on opposite strands, comprising nearly the full length of the viral genome. Such extensive overlap is not seen in other RNA viruses and comes at the cost of reduced evolutionary flexibility in the sequence. The present study is motivated by investigating the benefits which balance that cost. We show for the first time a functional requirement for the ambigrammatic genome configuration in Culex narnavirus 1, which suggests a model for how translation of both strands might benefit this virus. Our work highlights a new blueprint for viral persistence, distinct from strategies defined by canonical definitions of the coding strand.

ObjectiveTo identify associations between microbes and host genes in cats with feline chronic gingivostomatitis (FCGS), a debilitating inflammatory oral mucosal disease with no known cause, compared with healthy cats and cats with periodontitis (control cats).Animals19 control cats and 23 cats with FCGS.ProceduresAt least 1 caudal oral mucosal swab specimen was obtained from each cat. Each specimen underwent unbiased metatranscriptomic next-generation RNA sequencing (mNGS). Filtered mNGS reads were aligned to all known genetic sequences from all organisms and to the cat transcriptome. The relative abundances of microbial and host gene read alignments were compared between FCGS-affected cats and control cats and between FCGS-affected cats that did and did not clinically respond to primary treatment. Assembled feline calicivirus (FCV) genomes were compared with reverse transcription PCR (RT-PCR) primers commonly used to identify FCV.ResultsThe only microbe strongly associated with FCGS was FCV, which was detected in 21 of 23 FCGS-affected cats but no control cats. Problematic base pair mismatches were identified between the assembled FCV genomes and RT-PCR primers. Puma feline foamy virus was detected in 9 of 13 FCGS-affected cats that were refractory to treatment and 5 healthy cats but was not detected in FCGS-affected cats that responded to tooth extractions. The most differentially expressed genes in FCGS-affected cats were those associated with antiviral activity.Conclusions and clinical relevanceResults suggested that FCGS pathogenesis has a viral component. Many FCV strains may yield false-negative results on RT-PCR-based assays. Coinfection of FCGS-affected cats with FCV and puma feline foamy virus may adversely affect response to treatment.

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.

Modern day large-scale, high-density farming environments are inherently susceptible to viral outbreaks, inadvertently creating conditions that favor increased pathogen transmission and potential zoonotic spread. Metagenomic sequencing has proven to be a useful tool for characterizing the microbial burden in both people, livestock, and environmental samples. International efforts have been successful at characterizing pathogens in commercial farming environments, especially swine farms, however it is unclear whether the full extent of microbial agents have been adequately captured or is representative of farms elsewhere. To augment international efforts we performed metagenomic next-generation sequencing on nine swine slurry and three environmental samples from a United States of America (U.S.A.) farm operation, characterized the microbial composition of slurry, and identified novel viruses. We assembled a remarkable total of 1792 viral genomes, of which 554 were novel/divergent. We assembled 1637 Picobirnavirus genome segments, of which 538 are novel. In addition, we discovered 10 new viruses belonging to a novel taxon: porcine Statoviruses; which have only been previously reported in human, macaques, mouse, and cows. We assembled 3 divergent Posaviruses and 3 swine Picornaviruses. In addition to viruses described, we found other eukaryotic genera such as Entamoeba and Blastocystis, and bacterial genera such as Listeria, Treponema, Peptoclostridium and Bordetella in the slurry. Of these, two species Entamoeba histolytica and Listeria monocytogenes known to cause human disease were detected. Further, antimicrobial resistance genes such as tetracycline and MLS (macrolide, lincosamide, streptogramin) were also identified. Metagenomic surveillance in swine fecal slurry has great potential for novel and antimicrobial resistant pathogen detection.

Staphylococcus argenteus is a novel staphylococcal species associated with invasive disease. We report the first case of daptomycin/vancomycin-resistant S. argenteus, initially speciated as Staphylococcus aureus, that developed from repeated treatment with daptomycin for a complex vascular graft infection. Whole-genome sequencing of longitudinally collected isolates identified acquisition of MprF S337L, a mutation predicted to increase surface charge and repel cationic molecules.

Central nervous system B cells have several potential roles in multiple sclerosis (MS): secretors of proinflammatory cytokines and chemokines, presenters of autoantigens to T cells, producers of pathogenic antibodies, and reservoirs for viruses that trigger demyelination. To interrogate these roles, single-cell RNA sequencing (scRNA-Seq) was performed on paired cerebrospinal fluid (CSF) and blood from subjects with relapsing-remitting MS (RRMS; n = 12), other neurologic diseases (ONDs; n = 1), and healthy controls (HCs; n = 3). Single-cell immunoglobulin sequencing (scIg-Seq) was performed on a subset of these subjects and additional RRMS (n = 4), clinically isolated syndrome (n = 2), and OND (n = 2) subjects. Further, paired CSF and blood B cell subsets (RRMS; n = 7) were isolated using fluorescence activated cell sorting for bulk RNA sequencing (RNA-Seq). Independent analyses across technologies demonstrated that nuclear factor kappa B (NF-κB) and cholesterol biosynthesis pathways were activated, and specific cytokine and chemokine receptors were up-regulated in CSF memory B cells. Further, SMAD/TGF-β1 signaling was down-regulated in CSF plasmablasts/plasma cells. Clonally expanded, somatically hypermutated IgM+ and IgG1+ CSF B cells were associated with inflammation, blood-brain barrier breakdown, and intrathecal Ig synthesis. While we identified memory B cells and plasmablast/plasma cells with highly similar Ig heavy-chain sequences across MS subjects, similarities were also identified with ONDs and HCs. No viral transcripts, including from Epstein-Barr virus, were detected. Our findings support the hypothesis that in MS, CSF B cells are driven to an inflammatory and clonally expanded memory and plasmablast/plasma cell phenotype.

BackgroundTransfusion-related sepsis remains an important hospital infection control challenge. Investigation of septic transfusion events is often restricted by the limitations of bacterial culture in terms of time requirements and low yield in the setting of prior antibiotic administration.MethodsIn 3 gram-negative septic transfusion cases, we performed metagenomic next-generation sequencing (mNGS) of direct clinical blood specimens in addition to standard culture-based approaches utilized for infection control investigations. Pathogen detection leveraged IDSeq, a new open-access microbial bioinformatics portal. Phylogenetic analysis was performed to assess microbial genetic relatedness and understand transmission events.ResultsmNGS of direct clinical blood specimens afforded precision detection of pathogens responsible for each case of transfusion-related sepsis and enabled discovery of a novel Acinetobacter species in a platelet product that had become contaminated despite photochemical pathogen reduction. In each case, longitudinal assessment of pathogen burden elucidated the temporal sequence of events associated with each transfusion-transmitted infection. We found that informative data could be obtained from culture-independent mNGS of residual platelet products and leftover blood specimens that were either unsuitable or unavailable for culture or that failed to grow due to prior antibiotic administration. We additionally developed methods to enhance accuracy for detecting transfusion-associated pathogens that share taxonomic similarity to contaminants commonly found in mNGS library preparations.ConclusionsCulture-independent mNGS of blood products afforded rapid and precise assessment of pathogen identity, abundance, and genetic relatedness. Together, these challenging cases demonstrated the potential for metagenomics to advance existing methods for investigating transfusion-transmitted infections.

Paraneoplastic neurological disorders are immune-mediated diseases understood to manifest as part of a misdirected anti-tumor immune response. Paraneoplastic neurological disorder-associated autoantibodies can assist with diagnosis and enhance our understanding of tumor-associated immune processes. We designed a comprehensive library of 49-amino-acid overlapping peptides spanning the entire human proteome, including all splicing isoforms and computationally predicted coding regions. Using this library, we optimized a phage immunoprecipitation and sequencing protocol with multiple rounds of enrichment to create high-resolution epitope profiles in serum and cerebrospinal fluid (CSF) samples from patients suffering from two common paraneoplastic neurological disorders, the anti-Yo (n = 36 patients) and anti-Hu (n = 44 patients) syndromes. All (100%) anti-Yo patient samples yielded enrichment of peptides from the canonical anti-Yo (CDR2 and CDR2L) antigens, while 38% of anti-Hu patients enriched peptides deriving from the nELAVL (neuronal embryonic lethal abnormal vision like) family of proteins, the anti-Hu autoantigenic target. Among the anti-Hu patient samples that were positive for nELAVL, we noted a restricted region of immunoreactivity. To achieve single amino acid resolution, we designed a novel deep mutational scanning phage library encoding all possible single-point mutants targeting the reactive nELAVL region. This analysis revealed a distinct preference for the degenerate motif, RLDxLL, shared by ELAVL2, 3 and 4. Lastly, phage immunoprecipitation sequencing identified several known autoantigens in these same patient samples, including peptides deriving from the cancer-associated antigens ZIC and SOX families of transcription factors. Overall, this optimized phage immunoprecipitation sequencing library and protocol yielded the high-resolution epitope mapping of the autoantigens targeted in anti-Yo and anti-Hu encephalitis patients to date. The results presented here further demonstrate the utility and high-resolution capability of phage immunoprecipitation sequencing for both basic science and clinical applications and for better understanding the antigenic targets and triggers of paraneoplastic neurological disorders.

We studied the host transcriptional response to SARS-CoV-2 by performing metagenomic sequencing of upper airway samples in 238 patients with COVID-19, other viral or non-viral acute respiratory illnesses (ARIs). Compared to other viral ARIs, COVID-19 was characterized by a diminished innate immune response, with reduced expression of genes involved in toll-like receptor and interleukin signaling, chemokine binding, neutrophil degranulation and interactions with lymphoid cells. Patients with COVID-19 also exhibited significantly reduced proportions of neutrophils and macrophages, and increased proportions of goblet, dendritic and B-cells, compared to other viral ARIs. Using machine learning, we built 26-, 10- and 3-gene classifiers that differentiated COVID-19 from other acute respiratory illnesses with AUCs of 0.980, 0.950 and 0.871, respectively. Classifier performance was stable at low viral loads, suggesting utility in settings where direct detection of viral nucleic acid may be unsuccessful. Taken together, our results illuminate unique aspects of the host transcriptional response to SARS-CoV-2 in comparison to other respiratory viruses and demonstrate the feasibility of COVID-19 diagnostics based on patient gene expression.

The identification of autoantigens remains a critical challenge for understanding and treating autoimmune diseases. Autoimmune polyendocrine syndrome type 1 (APS1), a rare monogenic form of autoimmunity, presents as widespread autoimmunity with T and B cell responses to multiple organs. Importantly, autoantibody discovery in APS1 can illuminate fundamental disease pathogenesis, and many of the antigens found in APS1 extend to more common autoimmune diseases. Here, we performed proteome-wide programmable phage-display (PhIP-Seq) on sera from a cohort of people with APS1 and discovered multiple common antibody targets. These novel APS1 autoantigens exhibit tissue-restricted expression, including expression in enteroendocrine cells, pineal gland, and dental enamel. Using detailed clinical phenotyping, we find novel associations between autoantibodies and organ-restricted autoimmunity, including a link between anti-KHDC3L autoantibodies and premature ovarian insufficiency, and between anti-RFX6 autoantibodies and diarrheal-type intestinal dysfunction. Our study highlights the utility of PhIP-Seq for extensively interrogating antigenic repertoires in human autoimmunity and the importance of antigen discovery for improved understanding of disease mechanisms.

New antimalarial drugs are needed. The benzoxaborole AN13762 showed excellent activity against cultured Plasmodium falciparum, against fresh Ugandan P. falciparum isolates, and in murine malaria models. To gain mechanistic insights, we selected in vitro for P. falciparum isolates resistant to AN13762. In all of 11 independent selections with 100 to 200 nM AN13762, the 50% inhibitory concentration (IC50) increased from 18-118 nM to 180-890 nM, and whole-genome sequencing of resistant parasites demonstrated mutations in prodrug activation and resistance esterase (PfPARE). The introduction of PfPARE mutations led to a similar level of resistance, and recombinant PfPARE hydrolyzed AN13762 to the benzoxaborole AN10248, which has activity similar to that of AN13762 but for which selection of resistance was not readily achieved. Parasites further selected with micromolar concentrations of AN13762 developed higher-level resistance (IC50, 1.9 to 5.0 μM), and sequencing revealed additional mutations in any of 5 genes, 4 of which were associated with ubiquitination/sumoylation enzyme cascades; the introduction of one of these mutations, in SUMO-activating enzyme subunit 2, led to a similar level of resistance. The other gene mutated in highly resistant parasites encodes the P. falciparum cleavage and specificity factor homolog PfCPSF3, previously identified as the antimalarial target of another benzoxaborole. Parasites selected for resistance to AN13762 were cross-resistant with a close analog, AN13956, but not with standard antimalarials, AN10248, or other benzoxaboroles known to have different P. falciparum targets. Thus, AN13762 appears to have a novel mechanism of antimalarial action and multiple mechanisms of resistance, including loss of function of PfPARE preventing activation to AN10248, followed by alterations in ubiquitination/sumoylation pathways or PfCPSF3.IMPORTANCE Benzoxaboroles are under study as potential new drugs to treat malaria. One benzoxaborole, AN13762, has potent activity and promising features, but its mechanisms of action and resistance are unknown. To gain insights into these mechanisms, we cultured malaria parasites with nonlethal concentrations of AN13762 and generated parasites with varied levels of resistance. Parasites with low-level resistance had mutations in PfPARE, which processes AN13762 into an active metabolite; PfPARE mutations prevented this processing. Parasites with high-level resistance had mutations in any of a number of enzymes, mostly those involved in stress responses. Parasites selected for AN13762 resistance were not resistant to other antimalarials, suggesting novel mechanisms of action and resistance for AN13762, a valuable feature for a new class of antimalarial drugs.

SARS-CoV-2 infection is characterized by peak viral load in the upper airway prior to or at the time of symptom onset, an unusual feature that has enabled widespread transmission of the virus and precipitated a global pandemic. How SARS-CoV-2 is able to achieve high titer in the absence of symptoms remains unclear. Here, we examine the upper airway host transcriptional response in patients with COVID-19 (n = 93), other viral (n = 41) or non-viral (n = 100) acute respiratory illnesses (ARIs). Compared with other viral ARIs, COVID-19 is characterized by a pronounced interferon response but attenuated activation of other innate immune pathways, including toll-like receptor, interleukin and chemokine signaling. The IL-1 and NLRP3 inflammasome pathways are markedly less responsive to SARS-CoV-2, commensurate with a signature of diminished neutrophil and macrophage recruitment. This pattern resembles previously described distinctions between symptomatic and asymptomatic viral infections and may partly explain the propensity for pre-symptomatic transmission in COVID-19. We further use machine learning to build 27-, 10- and 3-gene classifiers that differentiate COVID-19 from other ARIs with AUROCs of 0.981, 0.954 and 0.885, respectively. Classifier performance is stable across a wide range of viral load, suggesting utility in mitigating false positive or false negative results of direct SARS-CoV-2 tests.

Narnaviruses have been described as positive-sense RNA viruses with a remarkably simple genome of ~3 kb, encoding only a highly conserved RNA-dependent RNA polymerase (RdRp). Many narnaviruses, however, are 'ambigrammatic' and harbour an additional uninterrupted open reading frame (ORF) covering almost the entire length of the reverse complement strand. No function has been described for this ORF, yet the absence of stops is conserved across diverse narnaviruses, and in every case the codons in the reverse ORF and the RdRp are aligned. The >3 kb ORF overlap on opposite strands, unprecedented among RNA viruses, motivates an exploration of the constraints imposed or alleviated by the codon alignment. Here, we show that only when the codon frames are aligned can all stop codons be eliminated from the reverse strand by synonymous single-nucleotide substitutions in the RdRp gene, suggesting a mechanism for de novo gene creation within a strongly conserved amino-acid sequence. It will be fascinating to explore what implications this coding strategy has for other aspects of narnavirus biology. Beyond narnaviruses, our rapidly expanding catalogue of viral diversity may yet reveal additional examples of this broadly-extensible principle for ambigrammatic-sequence development.

The growing prevalence of deadly microbes with resistance to previously life-saving drug therapies is a dire threat to human health. Detection of low abundance pathogen sequences remains a challenge for metagenomic Next Generation Sequencing (NGS). We introduce FLASH (Finding Low Abundance Sequences by Hybridization), a next-generation CRISPR/Cas9 diagnostic method that takes advantage of the efficiency, specificity and flexibility of Cas9 to enrich for a programmed set of sequences. FLASH-NGS achieves up to 5 orders of magnitude of enrichment and sub-attomolar gene detection with minimal background. We provide an open-source software tool (FLASHit) for guide RNA design. Here we applied it to detection of antimicrobial resistance genes in respiratory fluid and dried blood spots, but FLASH-NGS is applicable to all areas that rely on multiplex PCR.

We used metagenomic next-generation sequencing to longitudinally assess the gut microbiota and antimicrobial resistomes of international travelers to clarify global exchange of resistant organisms. Travel resulted in an increase in antimicrobial resistance genes and a greater proportion of Escherichia species within gut microbial communities without impacting diversity.

BackgroundMetagenomic next-generation sequencing (NGS) of cerebrospinal fluid (CSF) has the potential to identify a broad range of pathogens in a single test.MethodsIn 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.ResultsWe 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.ConclusionsRoutine 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.).

BackgroundDespite improved diagnostics, pulmonary pathogens in immunocompromised children frequently evade detection, leading to significant mortality. Therefore, we aimed to develop a highly sensitive metagenomic next-generation sequencing (mNGS) assay capable of evaluating the pulmonary microbiome and identifying diverse pathogens in the lungs of immunocompromised children.MethodsWe collected 41 lower respiratory specimens from 34 immunocompromised children undergoing evaluation for pulmonary disease at 3 children's hospitals from 2014-2016. Samples underwent mechanical homogenization, parallel RNA/DNA extraction, and metagenomic sequencing. Sequencing reads were aligned to the National Center for Biotechnology Information nucleotide reference database to determine taxonomic identities. Statistical outliers were determined based on abundance within each sample and relative to other samples in the cohort.ResultsWe identified a rich cross-domain pulmonary microbiome that contained bacteria, fungi, RNA viruses, and DNA viruses in each patient. Potentially pathogenic bacteria were ubiquitous among samples but could be distinguished as possible causes of disease by parsing for outlier organisms. Samples with bacterial outliers had significantly depressed alpha-diversity (median, 0.61; interquartile range [IQR], 0.33-0.72 vs median, 0.96; IQR, 0.94-0.96; P < .001). Potential pathogens were detected in half of samples previously negative by clinical diagnostics, demonstrating increased sensitivity for missed pulmonary pathogens (P < .001).ConclusionsAn optimized mNGS assay for pulmonary microbes demonstrates significant inoculation of the lower airways of immunocompromised children with diverse bacteria, fungi, and viruses. Potential pathogens can be identified based on absolute and relative abundance. Ongoing investigation is needed to determine the pathogenic significance of outlier microbes in the lungs of immunocompromised children with pulmonary disease.

Here, we report the coding-complete genome sequence of an avian metapneumovirus from a monk parakeet (Myiopsitta monachus), identified by metagenomic next-generation sequencing during an investigation into a disease outbreak in a captive parrot breeding facility. Based on divergence from known strains, this sequence represents a new subgroup of avian metapneumovirus.

Lower respiratory tract infections (LRTIs) lead to more deaths each year than any other infectious disease category. Despite this, etiologic LRTI pathogens are infrequently identified due to limitations of existing microbiologic tests. In critically ill patients, noninfectious inflammatory syndromes resembling LRTIs further complicate diagnosis. To address the need for improved LRTI diagnostics, we performed metagenomic next-generation sequencing (mNGS) on tracheal aspirates from 92 adults with acute respiratory failure and simultaneously assessed pathogens, the airway microbiome, and the host transcriptome. To differentiate pathogens from respiratory commensals, we developed a rules-based model (RBM) and logistic regression model (LRM) in a derivation cohort of 20 patients with LRTIs or noninfectious acute respiratory illnesses. When tested in an independent validation cohort of 24 patients, both models achieved accuracies of 95.5%. We next developed pathogen, microbiome diversity, and host gene expression metrics to identify LRTI-positive patients and differentiate them from critically ill controls with noninfectious acute respiratory illnesses. When tested in the validation cohort, the pathogen metric performed with an area under the receiver-operating curve (AUC) of 0.96 (95% CI, 0.86-1.00), the diversity metric with an AUC of 0.80 (95% CI, 0.63-0.98), and the host transcriptional classifier with an AUC of 0.88 (95% CI, 0.75-1.00). Combining these achieved a negative predictive value of 100%. This study suggests that a single streamlined protocol offering an integrated genomic portrait of pathogen, microbiome, and host transcriptome may hold promise as a tool for LRTI diagnosis.

Here we present a compendium of single-cell transcriptomic data from the model organism Mus musculus that comprises more than 100,000 cells from 20 organs and tissues. These data represent a new resource for cell biology, reveal gene expression in poorly characterized cell populations and enable the direct and controlled comparison of gene expression in cell types that are shared between tissues, such as T lymphocytes and endothelial cells from different anatomical locations. Two distinct technical approaches were used for most organs: one approach, microfluidic droplet-based 3'-end counting, enabled the survey of thousands of cells at relatively low coverage, whereas the other, full-length transcript analysis based on fluorescence-activated cell sorting, enabled the characterization of cell types with high sensitivity and coverage. The cumulative data provide the foundation for an atlas of transcriptomic cell biology.

RationaleQuantification of type 2 inflammation provided a molecular basis for heterogeneity in asthma. Non-type 2 pathways that contribute to asthma pathogenesis are not well understood.ObjectivesTo identify dysregulated pathways beyond type 2 inflammation.MethodsWe applied RNA sequencing to airway epithelial brushings obtained from subjects with stable mild asthma not on corticosteroids (n = 19) and healthy control subjects (n = 16). Sequencing reads were mapped to human and viral genomes. In the same cohort, and in a separate group with severe asthma (n = 301), we profiled blood gene expression with microarrays.Measurements and main resultsIn airway brushings from mild asthma on inhaled corticosteroids, RNA sequencing yielded 1,379 differentially expressed genes (false discovery rate

ObjectivesThe Stopping Cavities Trial investigated effectiveness and safety of 38% silver diamine fluoride in arresting caries lesions.Materials and methodsThe study was a double-blind randomized placebo-controlled superiority trial with 2 parallel groups. The sites were Oregon preschools. Sixty-six preschool children with ≥1 lesion were enrolled. Silver diamine fluoride (38%) or placebo (blue-tinted water), applied topically to the lesion. The primary endpoint was caries arrest (lesion inactivity, Nyvad criteria) 14-21days post intervention. Dental plaque was collected from all children, and microbial composition was assessed by RNA sequencing from 2 lesions and 1 unaffected surface before treatment and at follow-up for 3 children from each group.Results and conclusionAverage proportion of arrested caries lesions in the silver diamine fluoride group was higher (0.72; 95% CI; 0.55, 0.84) than in the placebo group (0.05; 95% CI; 0.00, 0.16). Confirmatory analysis using generalized estimating equation log-linear regression, based on the number of arrested lesions and accounting for the number of treated surfaces and length of follow-up, indicates the risk of arrested caries was significantly higher in the treatment group (relative risk, 17.3; 95% CI: 4.3 to 69.4). No harms were observed. RNA sequencing analysis identified no consistent changes in relative abundance of caries-associated microbes, nor emergence of antibiotic or metal resistance gene expression. Topical 38% silver diamine fluoride is effective and safe in arresting cavities in preschool children.Clinical significanceThe treatment is applicable to primary care practice and may reduce the burden of untreated tooth decay in the population.

IntroductionCurrently, the detection of pathogens or mutations associated with intraocular lymphomas heavily relies on prespecified, directed PCRs. With metagenomic deep sequencing (MDS), an unbiased high-throughput sequencing approach, all pathogens as well as all mutations present in the host's genome can be detected in the same small amount of ocular fluid.MethodsIn this cross-sectional case series, aqueous fluid samples from two patients were submitted to MDS to identify pathogens as well as common and rare cancer mutations.ResultsMDS of aqueous fluid from the first patient with vitreal lymphoma revealed the presence of both Epstein-Barr virus (HHV-4/EBV) and human herpes virus 8 (HHV-8) RNA. Aqueous fluid from the second patient with intraocular B-cell lymphoma demonstrated a less common mutation in the MYD88 gene associated with B-cell lymphoma.ConclusionMDS detects pathogens that, in some instances, may drive the development of intraocular lymphomas. Moreover, MDS is able to identify both common and rare mutations associated with lymphomas.

BackgroundThe widespread availability of next generation genome sequencing technologies has enabled a wide range of variant detection applications, especially in cancer and inborn genetic disorders. For model systems and microorganisms, the same technology may be used to discover the causative mutations for any phenotype, including those generated in response to chemical perturbation. In the case of pathogenic organisms, these approaches have allowed the determination of drug targets by means of resistance selection followed by genome sequencing.ResultsMinorityReport is open source software written in python that facilitates the comparison of any two sets of genome alignments for the purpose of rapidly identifying the spectrum of nonsynonymous changes, insertions or deletions, and copy number variations in a presumed mutant relative to its parent. Specifically, MinorityReport relates mapped sequence reads in SAM format output from any alignment tool for both the mutant and parent genome, relative to a reference genome, and produces the set of variants that distinguishes the mutant from the parent, all presented in an intuitive, straightforward report format. MinorityReport features tunable parameters for evaluating evidence and a scoring system that prioritizes reported variants based on relative proportions of read counts supporting the variant in the mutant versus parent data sets. The utility of MinorityReport is demonstrated using previously published publicly available data sets to find the determinants of resistance for novel anti-malarial drugs.ConclusionsMinorityReport is readily available (github: JeremyHorst/MinorityReport) to identify the genetic mechanisms of drug resistance in Plasmodium, genotype-phenotype relationships in human diads, or genomic variations between any two related organisms.

Here, we report the genome sequence of a divergent human rhinovirus C isolate identified from an infant with a severe community-acquired respiratory infection. RNA sequencing performed on an Illumina platform identified reads aligning to human rhinovirus species, which were de novo assembled to produce a coding-complete genome sequence.

Botanical, mycological, zoological, and prokaryotic species names follow the Linnaean format, consisting of an italicized Latinized binomen with a capitalized genus name and a lower case species epithet (e.g., Homo sapiens). Virus species names, however, do not follow a uniform format, and, even when binomial, are not Linnaean in style. In this thought exercise, we attempted to convert all currently official names of species included in the virus family Arenaviridae and the virus order Mononegavirales to Linnaean binomials, and to identify and address associated challenges and concerns. Surprisingly, this endeavor was not as complicated or time-consuming as even the authors of this article expected when conceiving the experiment. [Arenaviridae; binomials; ICTV; International Committee on Taxonomy of Viruses; Mononegavirales; virus nomenclature; virus taxonomy.].

Genetic manipulation of the deadly malaria parasite Plasmodium falciparum remains challenging, but the rise of CRISPR/Cas9-based genome editing tools is increasing the feasibility of altering this parasite's genome in order to study its biology. Of particular interest is the investigation of drug targets and drug resistance mechanisms, which have major implications for fighting malaria. We present a new method for introducing drug resistance mutations in P. falciparum without the use of plasmids or the need for cloning homologous recombination templates. We demonstrate this method by introducing edits into the sodium efflux channel PfATP4 by transfection of a purified CRISPR/Cas9-guide RNA ribonucleoprotein complex and a 200-nucleotide single-stranded oligodeoxynucleotide (ssODN) repair template. Analysis of whole genome sequencing data with the variant-finding program MinorityReport confirmed that only the intended edits were made, and growth inhibition assays confirmed that these mutations confer resistance to the antimalarial SJ733. The method described here is ideally suited for the introduction of mutations that confer a fitness advantage under selection conditions, and the novel finding that an ssODN can function as a repair template in P. falciparum could greatly simplify future editing attempts regardless of the nuclease used or the delivery method.

Objectives: The Stopping Cavities Trial investigated effectiveness and safety of 38% silver diamine fluoride in arresting caries lesions. Materials and Methods: Double-blind randomized placebo-controlled superiority trial with 2 parallel groups. Oregon preschools. 66 preschool children with ≥1 lesion. 38% silver diamine fluoride or placebo (blue-tinted water), applied topically to the lesion. The primary endpoint was caries arrest (lesion inactivity, Nyvad criteria) 14- 21 days post intervention. Dental plaque was collected from all children, and microbial composition was assessed by RNA sequencing from 2 lesions and 1 unaffected surface before treatment and at follow-up for 3 children from each group. Results and Conclusion: Mean fraction of arrested caries lesions in the silver diamine fluoride group was higher (0.72; 95% CI; 0.55, 0.84) than in the placebo group (0.05; 95% CI; 0.00, 0.16). Confirmatory analysis using generalized estimating equation log-linear regression, accounting for the number of treated surfaces and length of follow-up, indicated the fraction of arrested caries was significantly higher in the treatment group (relative risk, 17.3; 95% CI: 4.3 to 69.4). No harms were observed. RNA sequencing analysis identified no consistent changes in relative abundance of caries-associated microbes, nor emergence of antibiotic or metal resistance gene expression. Topical 38% silver diamine fluoride was effective and safe in arresting cavities in preschool children. The treatment is applicable to primary care practice and may reduce the burden of untreated tooth decay in the population. Trial Registration: ClinicalTrials.gov NCT02536040 . Clinical Significance In this clinical trial, 72% of caries lesions were arrested by silver diamine fluoride, with no harms. Contrary to the presumed antibacterial mechanism, lesion bacterial composition changed negligibly. This simple topical treatment is applicable to primary care practice and may reduce the burden of untreated tooth decay in the population.

It has come to our attention that there is an error in Fig. 3a for this article [1]. The correct version of Fig. 3a can be found below. The red markers now reflect the sequence differences. The text is correct. There was also a row omitted in Additional file 1: Table S1. The revised version can be found below.

BackgroundOcular infections remain a major cause of blindness and morbidity worldwide. While prognosis is dependent on the timing and accuracy of diagnosis, the etiology remains elusive in ~50 % of presumed infectious uveitis cases. The objective of this study is to determine if unbiased metagenomic deep sequencing (MDS) can accurately detect pathogens in intraocular fluid samples of patients with uveitis.MethodsThis is a proof-of-concept study, in which intraocular fluid samples were obtained from five subjects with known diagnoses, and one subject with bilateral chronic uveitis without a known etiology. Samples were subjected to MDS, and results were compared with those from conventional diagnostic tests. Pathogens were identified using a rapid computational pipeline to analyze the non-host sequences obtained from MDS.ResultsUnbiased MDS of intraocular fluid produced results concordant with known diagnoses in subjects with (n = 4) and without (n = 1) uveitis. Samples positive for Cryptococcus neoformans, Toxoplasma gondii, and herpes simplex virus 1 as tested by a Clinical Laboratory Improvement Amendments-certified laboratory were correctly identified with MDS. Rubella virus was identified in one case of chronic bilateral idiopathic uveitis. The subject's strain was most closely related to a German rubella virus strain isolated in 1992, one year before he developed a fever and rash while living in Germany. The pattern and the number of viral identified mutations present in the patient's strain were consistent with long-term viral replication in the eye.ConclusionsMDS can identify fungi, parasites, and DNA and RNA viruses in minute volumes of intraocular fluid samples. The identification of chronic intraocular rubella virus infection highlights the eye's role as a long-term pathogen reservoir, which has implications for virus eradication and emerging global epidemics.

Laboratories studying high-priority pathogens need comprehensive methods to confirm microbial species and strains while also detecting contamination. Metagenomic deep sequencing (MDS) inventories nucleic acids present in laboratory stocks, providing an unbiased assessment of pathogen identity, the extent of genomic variation, and the presence of contaminants. Double-stranded cDNA MDS libraries were constructed from RNA extracted from in vitro-passaged stocks of six viruses (La Crosse virus, Ebola virus, canine distemper virus, measles virus, human respiratory syncytial virus, and vesicular stomatitis virus). Each library was dual indexed and pooled for sequencing. A custom bioinformatics pipeline determined the organisms present in each sample in a blinded fashion. Single nucleotide variant (SNV) analysis identified viral isolates. We confirmed that (i) each sample contained the expected microbe, (ii) dual indexing of the samples minimized false assignments of individual sequences, (iii) multiple viral and bacterial contaminants were present, and (iv) SNV analysis of the viral genomes allowed precise identification of the viral isolates. MDS can be multiplexed to allow simultaneous and unbiased interrogation of mixed microbial cultures and (i) confirm pathogen identity, (ii) characterize the extent of genomic variation, (iii) confirm the cell line used for virus propagation, and (iv) assess for contaminating microbes. These assessments ensure the true composition of these high-priority reagents and generate a comprehensive database of microbial genomes studied in each facility. MDS can serve as an integral part of a pathogen-tracking program which in turn will enhance sample security and increase experimental rigor and precision. IMPORTANCE Both the integrity and reproducibility of experiments using select agents depend in large part on unbiased validation to ensure the correct identity and purity of the species in question. Metagenomic deep sequencing (MDS) provides the required level of validation by allowing for an unbiased and comprehensive assessment of all the microbes in a laboratory stock.