Virome Capture Sequencing Enables Sensitive Viral Diagnosis and Comprehensive Virome Analysis

Two Years of Viral Metagenomics in a Tertiary Diagnostics Unit: Evaluation of the First 105 Cases

Metagenomic next-generation sequencing (mNGS) can capture the full spectrum of viral pathogens in a specimen and has the potential to become an all-in-one solution for virus diagnostics. To date, clinical application is still in an early phase and limitations remain. Here, we evaluated the impact of viral mNGS for cases analyzed over two years in a tertiary diagnostics unit. High throughput mNGS was performed upon request by the treating clinician in cases where the etiology of infection remained unknown or the initial differential diagnosis was very broad. The results were compared to conventional routine testing regarding outcome and workload. In total, 163 specimens from 105 patients were sequenced. The main sample types were cerebrospinal fluid (34%), blood (33%) and throat swabs (10%). In the majority of the cases, viral encephalitis/meningitis or respiratory infection was suspected. In parallel, conventional virus diagnostic tests were performed (mean 18.5 individually probed targets/patients). mNGS detected viruses in 34 cases (32%). While often confirmatory, in multiple cases, the identified viruses were not included in the selected routine diagnostic tests. Two years of mNGS in a tertiary diagnostics unit demonstrated the advantages of a single, untargeted approach for comprehensive, rapid and efficient virus diagnostics, confirming the utility of mNGS in complementing current routine tests.

Diagnostic Testing of Neurologic Infections

Neuroinfectious diseases continue to cause morbidity and mortality worldwide, with many emerging or reemerging infections resulting in neurologic sequelae. Careful clinical evaluation coupled with appropriate laboratory investigations still forms the bedrock for making the correct etiologic diagnosis and implementing appropriate management. The treating physician needs to understand the individual test characteristics of each of the many conventional candidate-based diagnostics: culture, pathogen-specific polymerase chain reaction, antigen, antibody tests, used to diagnose the whole array of neuroinvasive infections. In addition, there is a growing need for more comprehensive, agnostic testing modalities that can identify a diversity of infections with a single assay.

Detection and characterization of a novel bat-borne coronavirus in Singapore using multiple molecular approaches

Bats are important reservoirs and vectors in the transmission of emerging infectious diseases. Many highly pathogenic viruses such as SARS-CoV and rabies-related lyssaviruses have crossed species barriers to infect humans and other animals. In this study we monitored the major roost sites of bats in Singapore, and performed surveillance for zoonotic pathogens in these bats. Screening of guano samples collected during the survey uncovered a bat coronavirus (Betacoronavirus) in Cynopterus brachyotis, commonly known as the lesser dog-faced fruit bat. Using a capture-enrichment sequencing platform, the full-length genome of the bat CoV was sequenced and found to be closely related to the bat coronavirus HKU9 species found in Leschenault’s rousette discovered in the Guangdong and Yunnan provinces.

Antibodies to Enteroviruses in Cerebrospinal Fluid of Patients with Acute Flaccid Myelitis

Acute flaccid myelitis (AFM) has caused motor paralysis in >560 children in the United States since 2014. The temporal association of enterovirus (EV) outbreaks with increases in AFM cases and reports of fever, respiratory, or gastrointestinal illness prior to AFM in >90% of cases suggest a role for infectious agents. Cerebrospinal fluid (CSF) from 14 AFM and 5 non-AFM patients with central nervous system (CNS) diseases in 2018 were investigated by viral-capture high-throughput sequencing (VirCapSeq-VERT system). These CSF and serum samples, as well as multiple controls, were tested for antibodies to human EVs using peptide microarrays. EV RNA was confirmed in CSF from only 1 adult AFM case and 1 non-AFM case. In contrast, antibodies to EV peptides were present in CSF of 11 of 14 AFM patients (79%), significantly higher than controls, including non-AFM patients (1/5 [20%]), children with Kawasaki disease (0/10), and adults with non-AFM CNS diseases (2/11 [18%]) (P = 0.023, 0.0001, and 0.0028, respectively). Six of 14 CSF samples (43%) and 8 of 11 sera (73%) from AFM patients were immunoreactive to an EV-D68-specific peptide, whereas the three control groups were not immunoreactive in either CSF (0/5, 0/10, and 0/11; P = 0.008, 0.0003, and 0.035, respectively) or sera (0/2, 0/8, and 0/5; P = 0.139, 0.002, and 0.009, respectively).IMPORTANCE The presence in cerebrospinal fluid of antibodies to EV peptides at higher levels than non-AFM controls supports the plausibility of a link between EV infection and AFM that warrants further investigation and has the potential to lead to strategies for diagnosis and prevention of disease.

Novel insights and therapeutic approaches in idiopathic multicentric Castleman disease

Castleman disease (CD) describes a heterogeneous group of hematologic disorders that share characteristic lymph node histopathology. Patients of all ages present with either a solitary enlarged lymph node (unicentric CD) or multicentric lymphadenopathy (MCD) with systemic inflammation, cytopenias, and life-threatening multiple organ dysfunction resulting from a cytokine storm often driven by interleukin 6 (IL-6). Uncontrolled human herpesvirus-8 (HHV-8) infection causes approximately 50% of MCD cases, whereas the etiology is unknown in the remaining HHV-8-negative/idiopathic MCD cases (iMCD). The limited understanding of etiology, cell types, and signaling pathways involved in iMCD has slowed development of treatments and contributed to historically poor patient outcomes. Here, recent progress for diagnosing iMCD, characterizing etio-pathogenesis, and advancing treatments are reviewed. Several clinicopathological analyses provided the evidence base for the first-ever diagnostic criteria and revealed distinct clinical subtypes: thrombocytopenia, anasarca, fever, reticulin fibrosis/renal dysfunction, organomegaly (iMCD-TAFRO) or iMCD-not otherwise specified (iMCD-NOS), which are both observed all over the world. In 2014, the anti-IL-6 therapy siltuximab became the first iMCD treatment approved by the US Food and Drug Administration, on the basis of a 34% durable response rate; consensus guidelines recommend it as front-line therapy. Recent cytokine and proteomic profiling has revealed normal IL-6 levels in many patients with iMCD and potential alternative driver cytokines. Candidate novel genomic alterations, dysregulated cell types, and signaling pathways have also been identified as candidate therapeutic targets. RNA sequencing for viral transcripts did not reveal novel viruses, HHV-8, or other viruses pathologically associated with iMCD. Despite progress, iMCD remains poorly understood. Further efforts to elucidate etiology, pathogenesis, and treatment approaches, particularly for siltuximab-refractory patients, are needed.

FDA-ARGOS is a database with public quality-controlled reference genomes for diagnostic use and regulatory science

FDA proactively invests in tools to support innovation of emerging technologies, such as infectious disease next generation sequencing (ID-NGS). Here, we introduce FDA-ARGOS quality-controlled reference genomes as a public database for diagnostic purposes and demonstrate its utility on the example of two use cases. We provide quality control metrics for the FDA-ARGOS genomic database resource and outline the need for genome quality gap filling in the public domain. In the first use case, we show more accurate microbial identification of Enterococcus avium from metagenomic samples with FDA-ARGOS reference genomes compared to non-curated GenBank genomes. In the second use case, we demonstrate the utility of FDA-ARGOS reference genomes for Ebola virus target sequence comparison as part of a composite validation strategy for ID-NGS diagnostic tests. The use of FDA-ARGOS as an in silico target sequence comparator tool combined with representative clinical testing could reduce the burden for completing ID-NGS clinical trials.

To be able to use infectious disease next generation sequencing as a diagnostic tool, appropriate reference datasets are required. Here, Sichtig et al. describe FDA-ARGOS, a reference database for high-quality microbial reference genomes, and demonstrate its utility on the example of two use cases.

Progression of Duchenne Cardiomyopathy Presenting with Chest Pain and Troponin Elevation

BACKGROUND

Improved neuromuscular and respiratory therapies have altered the natural history of Duchenne muscular dystrophy (DMD) such that the most common cause of mortality is progressive cardiomyopathy. Despite imaging evidence of progressive cardiomyopathy, troponin I (cTn) is not significantly elevated in asymptomatic DMD patients.

RESULTS

We describe eight boys with DMD evaluated for acute chest pain (ACP) and found to have acute cTn elevation with depressed left ventricular ejection fraction (LVEF). Of our eight patients, five presented with a primary complaint of ACP, while three presented with secondary myocardial injury in the context of systemic illness requiring hospitalization. Electrocardiograms showed diffuse ST changes and mean peak cTn level was 44±15.4 ng/mL (reference range <0.03 ng/mL). cTn levels normalized with only supportive care. Cardiac magnetic resonance imaging (CMR) was performed during the event on all but one patient, demonstrating increased late gadolinium enhancement (LGE) from 12.4±11.4% to 36.5±10.3% with associated deterioration of LVEF from 61±4.4% to 47.6±6.6% which remained depressed on follow-up CMR study (49.1±7.8%). All viral studies were negative. Additional investigations varied among patients, but no causative findings were demonstrated.

CONCLUSIONS

ACP with cTn elevation occurs in DMD boys and may be indicative of cardiomyopathy progression as evidenced by acute left ventricular dysfunction and development or progression of myocardial fibrosis. This clinical presentation is under recognized. These events may represent an important pathophysiological mechanism in cardiomyopathy progression.

The Intestinal Virome and Immunity

The composition of the human microbiome is considered a major source of interindividual variation in immunity and, by extension, susceptibility to diseases. Intestinal bacteria have been the major focus of research. However, diverse communities of viruses that infect microbes and the animal host cohabitate the gastrointestinal tract and collectively constitute the gut virome. Although viruses are typically investigated as pathogens, recent studies highlight a relationship between the host and animal viruses in the gut that is more akin to host-microbiome interactions and includes both beneficial and detrimental outcomes for the host. These viruses are likely sources of immune variation, both locally and extraintestinally. In this review, we describe the components of the gut virome, in particular mammalian viruses, and their ability to modulate host responses during homeostasis and disease.

Virome capture sequencing does not identify active viral infection in unicentric and idiopathic multicentric Castleman disease

Castleman disease (CD) describes a spectrum of heterogeneous disorders defined by characteristic lymph node histopathology. Enlarged lymph nodes demonstrating CD histopathology can occur in isolation (unicentric CD; UCD) sometimes accompanied by mild symptoms, or at multiple sites (multicentric CD, MCD) with systemic inflammation and cytokine-driven multi-organ dysfunction. The discovery that Kaposi sarcoma herpesvirus/human herpesvirus (HHV)-8 drives MCD in a subset of patients has led to the hypotheses that UCD and MCD patients with negative HHV-8 testing by conventional methods may represent false negatives, or that these cases are driven by another virus, known or unknown. To investigate these hypotheses, the virome capture sequencing for vertebrate viruses (VirCapSeq-VERT) platform was employed to detect RNA transcripts from known and novel viruses in fresh frozen lymph node tissue from CD patients (12 UCD, 11 HHV-8-negative MCD [idiopathic MCD; iMCD], and two HHV-8-positive MCD) and related diseases (three T cell lymphoma and three Hodgkin lymphoma). This assay detected HHV-8 in both HHV-8-positive cases; however, HHV-8 was not found in clinically HHV-8-negative iMCD or UCD cases. Additionally, no novel viruses were discovered, and no single known virus was detected with apparent association to HHV-8-negative CD cases. Herpesviridae family members, notably including Epstein-Barr virus (EBV), were detected in 7 out of 12 UCD and 5 of 11 iMCD cases with apparent correlations with markers of disease severity in iMCD. Analysis of a separate cohort of archival formalin-fixed, paraffin-embedded lymph node tissue by In situ hybridization revealed significantly fewer EBV-positive cells in UCD and iMCD compared to tissue from HHV-8-positive MCD and EBV-associated lymphoproliferative disorder. In an additional cohort, quantitative testing for EBV by PCR in peripheral blood during disease flare did not detect systemic EBV viremia, suggesting detection lymph node tissue is due to occult, local reactivation in UCD and iMCD. This study confirms that HHV-8 is not present in UCD and iMCD patients. Further, it fails to establish a clear association between any single virus, novel or known, and CD in HHV-8-negative cases. Given that distinct forms of CD exist with viral and non-viral etiological drivers, CD should be considered a group of distinct and separate diseases with heterogeneous causes worthy of further study.

Metagenomic Approaches for Public Health Surveillance of Foodborne Infections: Opportunities and Challenges

Foodborne disease surveillance in the United States is at a critical point. Clinical and diagnostic laboratories are using culture-independent diagnostic tests (CIDTs) to identify the pathogen causing foodborne illness from patient specimens. CIDTs are molecular tests that allow doctors to rapidly identify the bacteria causing illness within hours. CIDTs, unlike previous gold standard methods such as bacterial culture, do not produce an isolate that can be subtyped as part of the national molecular subtyping network for foodborne disease surveillance, PulseNet. Without subtype information, cases can no longer be linked using molecular data to identify potentially related cases that are part of an outbreak. In this review, we discuss the public health needs for a molecular subtyping approach directly from patient specimen and highlight different approaches, including amplicon and shotgun metagenomic sequencing.

Clinical metagenomics

Clinical metagenomic next-generation sequencing (mNGS), the comprehensive analysis of microbial and host genetic material (DNA and RNA) in samples from patients, is rapidly moving from research to clinical laboratories. This emerging approach is changing how physicians diagnose and treat infectious disease, with applications spanning a wide range of areas, including antimicrobial resistance, the microbiome, human host gene expression (transcriptomics) and oncology. Here, we focus on the challenges of implementing mNGS in the clinical laboratory and address potential solutions for maximizing its impact on patient care and public health.

Recent and emerging technologies for the rapid diagnosis of infection and antimicrobial resistance

The rise in antimicrobial resistance (AMR) is predicted to cause 10 million deaths per year by 2050 unless steps are taken to prevent this looming crisis. Microbiological culture is the gold standard for the diagnosis of bacterial/fungal pathogens and antimicrobial resistance and takes 48 hours or longer. Hence, antibiotic prescriptions are rarely based on a definitive diagnosis and patients often receive inappropriate treatment. Rapid diagnostic tools are urgently required to guide appropriate antimicrobial therapy, thereby improving patient outcomes and slowing AMR development. We discuss new technologies for rapid infection diagnosis including: sample-in-answer-out PCR-based tests, BioFire FilmArray and Curetis Unyvero; rapid susceptibility tests, Accelerate Pheno and microfluidic tests; and sequencing-based approaches, focusing on targeted and clinical metagenomic nanopore sequencing.

Application of a targeted-enrichment methodology for full-genome sequencing of Dengue 1-4, Chikungunya and Zika viruses directly from patient samples

The frequency of epidemics caused by Dengue viruses 1-4, Zika virus and Chikungunya viruses have been on an upward trend in recent years driven primarily by uncontrolled urbanization, mobility of human populations and geographical spread of their shared vectors, Aedes aegypti and Aedes albopictus. Infections by these viruses present with similar clinical manifestations making them challenging to diagnose; this is especially difficult in regions of the world hyperendemic for these viruses. In this study, we present a targeted-enrichment methodology to simultaneously sequence the complete viral genomes for each of these viruses directly from clinical samples. Additionally, we have also developed a customized computational tool (BaitMaker) to design these enrichment baits. This methodology is robust in its ability to capture diverse sequences and is amenable to large-scale epidemiological studies. We have applied this methodology to two large cohorts: a febrile study based in Colombo, Sri Lanka taken during the 2009-2015 dengue epidemic (n = 170) and another taken during the 2016 outbreak of Zika virus in Singapore (n = 162). Results from these studies indicate that we were able to cover an average of 97.04% ± 0.67% of the full viral genome from samples in these cohorts. We also show detection of one DENV3/ZIKV co-infected patient where we recovered full genomes for both viruses.

A Novel Orthohepadnavirus Identified in a Dead Maxwell's Duiker (Philantomba maxwellii) in Taï National Park, Côte d'Ivoire

New technologies enable viral discovery in a diversity of hosts, providing insights into viral evolution. We used one such approach, the virome capture sequencing for vertebrate viruses (VirCapSeq-VERT) platform, on 21 samples originating from six dead Maxwell’s duikers (Philantomba maxwellii) from Taï National Park, Côte d’Ivoire. We detected the presence of an orthohepadnavirus in one animal and characterized its 3128 bp genome. The highest viral copy numbers were detected in the spleen, followed by the lung, blood, and liver, with the lowest copy numbers in the kidney and heart; the virus was not detected in the jejunum. Viral copy numbers in the blood were in the range known from humans with active chronic infections leading to liver histolytic damage, suggesting this virus could be pathogenic in duikers, though many orthohepadnaviruses appear to be apathogenic in other hosts, precluding a formal test of this hypothesis. The virus was not detected in 29 other dead duiker samples from the Côte d’Ivoire and Central African Republic, suggesting either a spillover event or a low prevalence in these populations. Phylogenetic analysis placed the virus as a divergent member of the mammalian clade of orthohepadnaviruses, though its relationship to other orthohepadnaviruses remains uncertain. This represents the first orthohepadnavirus described in an artiodactyl. We have tentatively named this new member of the genus Orthohepadnavirus (family Hepadnaviridae), Taï Forest hepadnavirus. Further studies are needed to determine whether it, or some close relatives, are present in a broader range of artiodactyls, including livestock.

Emergence of the East-Central-South-African genotype of Chikungunya virus in Brazil and the city of Rio de Janeiro may have occurred years before surveillance detection

Brazil, which is hyperendemic for dengue virus (DENV), has had recent Zika (ZIKV) and (CHIKV) Chikungunya virus outbreaks. Since March 2016, CHIKV is the arbovirus infection most frequently diagnosed in Rio de Janeiro. In the analysis of 1835 syndromic patients, screened by real time RT-PCR, 56.4% of the cases were attributed to CHIKV, 29.6% to ZIKV, and 14.1% to DENV-4. Sequence analyses of CHIKV from sixteen samples revealed that the East-Central-South-African (ECSA) genotype of CHIKV has been circulating in Brazil since 2013 [95% bayesian credible interval (BCI): 03/2012-10/2013], almost a year before it was detected by arbovirus surveillance program. Brazilian cases are related to Central African Republic sequences from 1980’s. To the best of our knowledge, given the available sequence published here and elsewhere, the ECSA genotype was likely introduced to Rio de Janeiro early on 2014 (02/2014; BCI: 07/2013-08/2014) through a single event, after primary circulation in the Bahia state at the Northestern Brazil in the previous year. The observation that the ECSA genotype of CHIKV was circulating undetected underscores the need for improvements in molecular methods for viral surveillance.

The Evolution of Genotyping Strategies To Detect, Analyze, and Control Transmission of Tuberculosis

The introduction of genotypic tools to analyze Mycobacterium tuberculosis isolates has transformed our knowledge of the transmission dynamics of this pathogen. We discuss the development of the laboratory methods that have been applied in recent years to study the epidemiology of M. tuberculosis. This review integrates two approaches: on the one hand, it considers how genotyping techniques have evolved over the years; and on the other, it looks at how the way we think these techniques should be applied has changed. We begin by examining the application of fingerprinting tools to suspected outbreaks only, before moving on to universal genotyping schemes, and finally we describe the latest real-time strategies used in molecular epidemiology. We also analyze refined approaches to obtaining epidemiological data from patients and to increasing the discriminatory power of genotyping by techniques based on genomic characterization. Finally, we review the development of integrative solutions to reconcile the speed of PCR-based methods with the high discriminatory power of whole-genome sequencing in easily implemented formats adapted to low-resource settings. Our analysis of future considerations highlights the need to bring together the three key elements of high-quality surveillance of transmission in tuberculosis, namely, speed, precision, and ease of implementation.

Higher abundance of enterovirus A species in the gut of children with islet autoimmunity

Enteroviruses (EVs) are prime candidate environmental triggers of islet autoimmunity (IA), with potential as vaccine targets for type 1 diabetes prevention. However, the use of targeted virus detection methods and the selective focus on EVs by most studies increases the risk for substantial investigation bias and an overestimated association between EV and type 1 diabetes. Here we performed comprehensive virome-capture sequencing to examine all known vertebrate-infecting viruses without bias in 182 specimens (faeces and plasma) collected before or at seroconversion from 45 case children with IA and 48 matched controls. From >2.6 billion reads, 28 genera of viruses were detected and 62% of children (58/93) were positive for ≥1 vertebrate-infecting virus. We identified 129 viruses as differentially abundant between the gut of cases and controls, including 5 EV-A types significantly more abundant in the cases. Our findings further support EV’s hypothesised contribution to IA and corroborate the proposal that viral load may be an important parameter in disease pathogenesis. Furthermore, our data indicate a previously unrecognised association of IA with higher EV-A abundance in the gut of children and provide a catalog of viruses to be interrogated further to determine a causal link between virus infection and type 1 diabetes.

Albugo candida race diversity, ploidy and host-associated microbes revealed using DNA sequence capture on diseased plants in the field

Physiological races of the oomycete Albugo candida are biotrophic pathogens of diverse plant species, primarily the Brassicaceae, and cause infections that suppress host immunity to other pathogens. However, A. candida race diversity and the consequences of host immunosuppression are poorly understood in the field. We report a method that enables sequencing of DNA of plant pathogens and plant-associated microbes directly from field samples (Pathogen Enrichment Sequencing: PenSeq). We apply this method to explore race diversity in A. candida and to detect A. candida-associated microbes in the field (91 A. candida-infected plants). We show with unprecedented resolution that each host plant species supports colonization by one of 17 distinct phylogenetic lineages, each with an unique repertoire of effector candidate alleles. These data reveal the crucial role of sexual and asexual reproduction, polyploidy and host domestication in A. candida specialization on distinct plant species. Our bait design also enabled phylogenetic assignment of DNA sequences from bacteria and fungi from plants in the field. This paper shows that targeted sequencing has a great potential for the study of pathogen populations while they are colonizing their hosts. This method could be applied to other microbes, especially to those that cannot be cultured.

Distinct Gut Virome Profile of Pregnant Women With Type 1 Diabetes in the ENDIA Study

Background

The importance of gut bacteria in human physiology, immune regulation, and disease pathogenesis is well established. In contrast, the composition and dynamics of the gut virome are largely unknown; particularly lacking are studies in pregnancy. We used comprehensive virome capture sequencing to characterize the gut virome of pregnant women with and without type 1 diabetes (T1D), longitudinally followed in the Environmental Determinants of Islet Autoimmunity study.

Methods

In total, 61 pregnant women (35 with T1D and 26 without) from Australia were examined. Nucleic acid was extracted from serial fecal specimens obtained at prenatal visits, and viral genomes were sequenced by virome capture enrichment. The frequency, richness, and abundance of viruses were compared between women with and without T1D.

Results

Two viruses were more prevalent in pregnant women with T1D: picobirnaviruses (odds ratio [OR], 4.2; 95% confidence interval [CI], 1.0–17.1; P = .046) and tobamoviruses (OR, 3.2; 95% CI, 1.1–9.3; P = .037). The abundance of 77 viruses significantly differed between the 2 maternal groups (≥2-fold difference; P < .02), including 8 Enterovirus B types present at a higher abundance in women with T1D.

Conclusions

These findings provide novel insight into the composition of the gut virome during pregnancy and demonstrate a distinct profile of viruses in women with T1D.

Critical issues in application of molecular methods to environmental virology

Waterborne diseases have significant public health and socioeconomic implications worldwide. Many viral pathogens are commonly associated with water-related diseases, namely enteric viruses. Also, novel recently discovered human-associated viruses have been shown to be a causative agent of gastroenteritis or other clinical symptoms. A wide range of analytical methods is available for virus detection in environmental water samples. Viral isolation is historically carried out via propagation on permissive cell lines; however, some enteric viruses are difficult or not able to propagate on existing cell lines. Real-time polymerase chain reaction (qPCR) screening of viral nucleic acid is routinely used to investigate virus contamination in water due to the high sensitivity and specificity. Additionally, the introduction of metagenomic approaches into environmental virology has facilitated the discovery of viruses that cannot be grown in cell culture. This review (i) highlights the applications of molecular techniques in environmental virology such as PCR and its modifications to overcome the critical issues associated with the inability to discriminate between infectious viruses and nonviable viruses, (ii) outlines the strengths and weaknesses of Nucleic Acid Sequence Based Amplification (NASBA) and microarray, (iii) discusses the role of digital PCR as an emerging water quality monitoring assay and its advantages over qPCR, (iv) addresses the viral metagenomics in terms of detecting emerging viral pathogens and diversity in aquatic environment. Indeed, there are many challenges for selecting methods to detect classic and emerging viruses in environmental samples. While the existing techniques have revealed the importance and diversity of viruses in the water environment, further developments are necessary to enable more rapid and accurate methodologies for viral water quality monitoring and regulation.

Complete assembly of a dengue virus type 3 genome from a recent genotype III clade by metagenomic sequencing of serum

Background: Mosquito-borne flaviviruses, such as dengue and Japanese encephalitis virus (JEV), cause life-threatening diseases, particularly in the tropics. Methods: Here we performed unbiased metagenomic sequencing of RNA extracted from the serum of four patients and the plasma of one patient, all hospitalized at a tertiary care centre in South India with severe or prolonged febrile illness, together with the serum from one healthy control, in 2014. Results: We identified and assembled a complete dengue virus type 3 sequence from a case of severe dengue fever. We also identified a small number of JEV sequences in the serum of two adults with febrile illness, including one with severe dengue. Phylogenetic analysis revealed that the dengue sequence belonged to genotype III. It has an estimated divergence time of 13.86 years from the most highly related Indian strains. In total, 11 amino acid substitutions were predicted for this strain in the antigenic envelope protein, when compared to the parent strain used for development of the first commercial dengue vaccine.  Conclusions: We demonstrate that both genome assembly and detection of a low number of viral sequences are possible through the unbiased sequencing of clinical material. These methods may help ascertain causal agents for febrile illnesses with no known cause.

Capturing sequence diversity in metagenomes with comprehensive and scalable probe design

Metagenomic sequencing has the potential to transform microbial detection and characterization, but new tools are needed to improve its sensitivity. Here we present CATCH, a computational method to enhance nucleic acid capture for enrichment of diverse microbial taxa. CATCH designs optimal probe sets, with a specified number of oligonucleotides, that achieve full coverage of, and scale well with, known sequence diversity. We focus on applying CATCH to capture viral genomes in complex metagenomic samples. We design, synthesize, and validate multiple probe sets, including one that targets the whole genomes of the 356 viral species known to infect humans. Capture with these probe sets enriches unique viral content on average 18-fold, allowing us to assemble genomes that could not be recovered without enrichment, and accurately preserves within-sample diversity. We also use these probe sets to recover genomes from the 2018 Lassa fever outbreak in Nigeria and to improve detection of uncharacterized viral infections in human and mosquito samples. The results demonstrate that CATCH enables more sensitive and cost-effective metagenomic sequencing.

Virus genotyping by massive parallel amplicon sequencing: adenovirus and enterovirus in the Norwegian MIDIA study

OBJECTIVES

Direct genotyping of adenovirus or enterovirus from clinical material using polymerase chain reaction (PCR) followed by Sanger sequencing is often difficult due to the presence of multiple virus types in a sample, or due to varying efficacy of PCR amplifying the capsid gene on the background of foreign nucleic acids. Here we present a simple protocol for virus genotyping using massive parallel amplicon sequencing.

METHODS

The protocol utilized a set of 16 tailed degenerate primers flanking the seventh hypervariable region of the adenovirus hexon gene and 9 tailed degenerate primers targeted to the proximal portion of the enterovirus VP1 gene. Subsequent addition of dual indices enabled simultaneous sequencing of 384 different samples on an Illumina MiSeq instrument. Downstream bioinformatic analysis was based on remapping to a set of references representative of the presently known repertoire of virus types.

RESULTS

After validation with known virus types, the sequencing method was applied on 301 adenovirus-positive samples and 350 enterovirus-positive samples from a longitudinally collected series of stools from 83 children aged 3 to 36 months. We detected 7 different adenovirus types and 27 different enterovirus types. There were 37 (6.2%) samples containing more than one genotype of the same viral genus. At least one dual infection was experienced by 23 of 83 (28%) of the children observed over the 3 years' observation period.

CONCLUSIONS

Amplicon sequencing with a multiplex set of degenerate primers seems to be a rapid and reliable technical solution for genotyping of large collections of samples where simultaneous infections with multiple strains can be expected.

From Theory to Practice: Translating Whole-Genome Sequencing (WGS) into the Clinic

Hospitals worldwide are facing an increasing incidence of hard-to-treat infections. Limiting infections and providing patients with optimal drug regimens require timely strain identification as well as virulence and drug-resistance profiling. Additionally, prophylactic interventions based on the identification of environmental sources of recurrent infections (e.g., contaminated sinks) and reconstruction of transmission chains (i.e., who infected whom) could help to reduce the incidence of nosocomial infections. WGS could hold the key to solving these issues. However, uptake in the clinic has been slow. Some major scientific and logistical challenges need to be solved before WGS fulfils its potential in clinical microbial diagnostics. In this review we identify major bottlenecks that need to be resolved for WGS to routinely inform clinical intervention and discuss possible solutions.

Novel insights and therapeutic approaches in idiopathic multicentric Castleman disease

Castleman disease (CD) describes a heterogeneous group of hematologic disorders that share characteristic lymph node histopathology. Patients of all ages present with either a solitary enlarged lymph node (unicentric CD) or multicentric lymphadenopathy (MCD) with systemic inflammation, cytopenias, and life-threatening multiple organ dysfunction resulting from a cytokine storm often driven by interleukin 6 (IL-6). Uncontrolled human herpesvirus-8 (HHV-8) infection causes approximately 50% of MCD cases, whereas the etiology is unknown in the remaining HHV-8-negative/idiopathic MCD cases (iMCD). The limited understanding of etiology, cell types, and signaling pathways involved in iMCD has slowed development of treatments and contributed to historically poor patient outcomes. Here, recent progress for diagnosing iMCD, characterizing etio-pathogenesis, and advancing treatments are reviewed. Several clinicopathological analyses provided the evidence base for the first-ever diagnostic criteria and revealed distinct clinical subtypes: thrombocytopenia, anasarca, fever, reticulin fibrosis/renal dysfunction, organomegaly (iMCD-TAFRO) or iMCD-not otherwise specified (iMCD-NOS), which are both observed all over the world. In 2014, the anti-IL-6 therapy siltuximab became the first iMCD treatment approved by the US Food and Drug Administration, on the basis of a 34% durable response rate; consensus guidelines recommend it as front-line therapy. Recent cytokine and proteomic profiling has revealed normal IL-6 levels in many patients with iMCD and potential alternative driver cytokines. Candidate novel genomic alterations, dysregulated cell types, and signaling pathways have also been identified as candidate therapeutic targets. RNA sequencing for viral transcripts did not reveal novel viruses, HHV-8, or other viruses pathologically associated with iMCD. Despite progress, iMCD remains poorly understood. Further efforts to elucidate etiology, pathogenesis, and treatment approaches, particularly for siltuximab-refractory patients, are needed.

Detection of Viruses in Clinical Samples by Use of Metagenomic Sequencing and Targeted Sequence Capture

Metagenomic shotgun sequencing (MSS) is a revolutionary approach to viral diagnostic testing that allows simultaneous detection of a broad range of viruses, detailed taxonomic assignment, and detection of mutations associated with antiviral drug resistance. To enhance sensitivity for virus detection, we previously developed ViroCap, a targeted sequence capture panel designed to enrich nucleic acid from a comprehensive set of eukaryotic viruses prior to sequencing. To demonstrate the utility of MSS with targeted sequence capture for detecting clinically important viruses and characterizing clinically important viral features, we used ViroCap to analyze clinical samples from a diagnostic virology laboratory containing a broad range of medically relevant viruses. From 26 samples, MSS with ViroCap detected all of the expected viruses and 30 additional viruses. Comparing sequencing after capture enrichment with standard MSS, we detected 13 viruses only with capture enrichment and observed a consistent increase in the number and percentage of viral sequence reads as well as the breadth and depth of coverage of the viral genomes. Compared with clinical testing, MSS enhanced taxonomic assignment for 15 viruses, and codons associated with antiviral drug resistance in influenza A virus, herpes simplex virus (HSV), human immunodeficiency virus (HIV), and hepatitis C virus (HCV) could be analyzed. Overall, in clinical samples, MSS with targeted sequence capture provides enhanced virus detection and information of clinical and epidemiologic relevance compared with clinical testing and MSS without targeted sequence capture.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.

Gut microbes as future therapeutics in treating inflammatory and infectious diseases: Lessons from recent findings

The human gut microbiota has been the interest of extensive research in recent years and our knowledge on using the potential capacity of these microbes are growing rapidly. Microorganisms colonized throughout the gastrointestinal tract of human are coevolved through symbiotic relationship and can influence physiology, metabolism, nutrition and immune functions of an individual. The gut microbes are directly involved in conferring protection against pathogen colonization by inducing direct killing, competing with nutrients and enhancing the response of the gut-associated immune repertoire. Damage in the microbiome (dysbiosis) is linked with several life-threatening outcomes viz. inflammatory bowel disease, cancer, obesity, allergy, and auto-immune disorders. Therefore, the manipulation of human gut microbiota came out as a potential choice for therapeutic intervention of the several human diseases. Herein, we review significant studies emphasizing the influence of the gut microbiota on the regulation of host responses in combating infectious and inflammatory diseases alongside describing the promises of gut microbes as future therapeutics.

Optimizing a Metatranscriptomic Next-Generation Sequencing Protocol for Bronchoalveolar Lavage Diagnostics

Compared with conventional serologic, culture-based, and molecular-based diagnostic tests, next-generation sequencing (NGS) provides sequence-evidenced detection of various microbes, without prior knowledge, and thus is becoming a novel diagnostic approach. Herein we describe an RNA-based metatranscriptomic NGS (mtNGS) protocol for culture-independent detection of potential infectious pathogens, using clinical bronchoalveolar lavage specimens as an example. We present both an optimized workflow for experimental sequence data collection and a simplified pipeline for bioinformatics sequence data processing. As shown, the whole protocol takes approximately 24 to 36 hours to detect a wide range of Gram-positive and -negative bacteria and possibly other viral and/or fungal pathogens. In particular, we introduce a spike-in RNA mix as an internal control, which plays a critical role in mitigating false-positive and false-negative results of clinical diagnostic tests. Moreover, our mtNGS method can detect antibiotic resistance genes and virulence factors; although it may not be comprehensive, such information is imperative and helpful for the clinician to make better treatment decisions. Results from our preliminary testing suggest that the mtNGS approach is a useful alterative in diagnostic detection of emerging infectious pathogens in clinical laboratories. However, further improvements are needed to achieve better sensitivity and accuracy.

Virome Analysis Reveals No Association of Head and Neck Vascular Anomalies with an Active Viral Infection

BACKGROUND/AIM

Vascular anomalies encompass different vascular malformations [arteriovenous (AVM), lymphatic (LM), venous lymphatic (VLM), venous (VM)] and vascular tumors such as hemangiomas (HA). The pathogenesis of vascular anomalies is still poorly understood. Viral infection was speculated as a possible underlying cause.

MATERIALS AND METHODS

A total of 13 human vascular anomalies and three human skin control tissues were used for viral analysis. RNA derived from AVM (n=4) and normal skin control (n=3) tissues was evaluated by RNA sequencing. The Virome Capture Sequencing Platform for Vertebrate Viruses (VirCapSeq-VERT) was deployed on 10 tissues with vascular anomalies (2×AVM, 1×HA, 1×LM, 2×VLM, 4×VM).

RESULTS

RNA sequencing did not show any correlation of AVM with viral infection. By deploying VirCapSeq-VERT, no consistent viral association was seen in the tested tissues.

CONCLUSION

The analysis does not point to the presence of an active viral infection in vascular anomalies. However, transient earlier viral infections, e.g. during pregnancy, cannot be excluded with this approach.

BacCapSeq: a Platform for Diagnosis and Characterization of Bacterial Infections

BacCapSeq is a method for differential diagnosis of bacterial infections and defining antimicrobial sensitivity profiles that has the potential to reduce morbidity and mortality, health care costs, and the inappropriate use of antibiotics that contributes to the development of antimicrobial resistance.

We report a platform that increases the sensitivity of high-throughput sequencing for detection and characterization of bacteria, virulence determinants, and antimicrobial resistance (AMR) genes. The system uses a probe set comprised of 4.2 million oligonucleotides based on the Pathosystems Resource Integration Center (PATRIC) database, the Comprehensive Antibiotic Resistance Database (CARD), and the Virulence Factor Database (VFDB), representing 307 bacterial species that include all known human-pathogenic species, known antimicrobial resistance genes, and known virulence factors, respectively. The use of bacterial capture sequencing (BacCapSeq) resulted in an up to 1,000-fold increase in bacterial reads from blood samples and lowered the limit of detection by 1 to 2 orders of magnitude compared to conventional unbiased high-throughput sequencing, down to a level comparable to that of agent-specific real-time PCR with as few as 5 million total reads generated per sample. It detected not only the presence of AMR genes but also biomarkers for AMR that included both constitutive and differentially expressed transcripts.

The discovery of Bombali virus adds further support for bats as hosts of ebolaviruses

Here we describe the complete genome of a new ebolavirus, Bombali virus (BOMV) detected in free-tailed bats in Sierra Leone (little free-tailed (Chaerephon pumilus) and Angolan free-tailed (Mops condylurus)). The bats were found roosting inside houses, indicating the potential for human transmission. We show that the viral glycoprotein can mediate entry into human cells. However, further studies are required to investigate whether exposure has actually occurred or if BOMV is pathogenic in humans.

Epstein-Barr Virus-Positive Cancers Show Altered B-Cell Clonality

Around 20% of human cancers are associated with viruses. Epstein-Barr virus (EBV) contributes to gastric cancer, nasopharyngeal carcinoma, and certain lymphomas, but its role in other cancer types remains controversial. We assessed the prevalence of EBV in RNA-seq from 32 tumor types in the Cancer Genome Atlas Project (TCGA) and found EBV to be present in >5% of samples in 12 tumor types. EBV infects epithelial cells and B cells and in B cells causes proliferation. We hypothesized that the low expression of EBV in most of the tumor types was due to infiltration of B cells into the tumor. The increase in B-cell abundance and diversity in subjects where EBV was detected in the tumors strengthens this hypothesis. Overall, we found that EBV was associated with an increased and altered immune response. This result is not evidence of causality, but a potential novel biomarker for tumor immune status.

Epstein-Barr virus (EBV) is convincingly associated with gastric cancer, nasopharyngeal carcinoma, and certain lymphomas, but its role in other cancer types remains controversial. To test the hypothesis that there are additional cancer types with high prevalence of EBV, we determined EBV viral expression in all the Cancer Genome Atlas Project (TCGA) mRNA sequencing (mRNA-seq) samples (n = 10,396) from 32 different tumor types. We found that EBV was present in gastric adenocarcinoma and lymphoma, as expected, and was also present in >5% of samples in 10 additional tumor types. For most samples, EBV transcript levels were low, which suggests that EBV was likely present due to infected infiltrating B cells. In order to determine if there was a difference in the B-cell populations, we assembled B-cell receptors for each sample and found B-cell receptor abundance (P ≤ 1.4 × 10⁻²⁰) and diversity (P ≤ 8.3 × 10⁻²⁷) were significantly higher in EBV-positive samples. Moreover, diversity was independent of B-cell abundance, suggesting that the presence of EBV was associated with an increased and altered B-cell population.

IMPORTANCE Around 20% of human cancers are associated with viruses. Epstein-Barr virus (EBV) contributes to gastric cancer, nasopharyngeal carcinoma, and certain lymphomas, but its role in other cancer types remains controversial. We assessed the prevalence of EBV in RNA-seq from 32 tumor types in the Cancer Genome Atlas Project (TCGA) and found EBV to be present in >5% of samples in 12 tumor types. EBV infects epithelial cells and B cells and in B cells causes proliferation. We hypothesized that the low expression of EBV in most of the tumor types was due to infiltration of B cells into the tumor. The increase in B-cell abundance and diversity in subjects where EBV was detected in the tumors strengthens this hypothesis. Overall, we found that EBV was associated with an increased and altered immune response. This result is not evidence of causality, but a potential novel biomarker for tumor immune status.

Deep Sequencing of H7N9 Influenza A Viruses from 16 Infected Patients from 2013 to 2015 in Shanghai Reveals Genetic Diversity and Antigenic Drift

Influenza A virus (IAV) infections are a major public health concern, including annual epidemics, epizootic outbreaks, and pandemics. A significant IAV epizootic outbreak was the H7N9 avian influenza A outbreak in China, which was first detected in 2013 and which has spread over 5 waves from 2013 to 2017, causing human infections in many different Chinese provinces. Here, RNA from primary clinical throat swab samples from 20 H7N9-infected local patients with different clinical outcomes, who were admitted and treated at one hospital in Shanghai, China, from April 2013 to April 2015, was analyzed. Whole-transcriptome amplification, with positive enrichment of IAV RNA, was performed, all 20 samples were subjected to deep sequencing, and data from 16 samples were analyzed in detail. Many single-nucleotide polymorphisms, including ones not previously reported, and many nonsynonymous changes that could affect hemagglutinin head and stalk antibody binding epitopes were observed. Minor populations representing viral quasispecies, including nonsynonymous hemagglutinin changes shared by antigenically variant H7N9 clades identified in the most recent wave of H7N9 infections in 2016 to 2017, were also identified.IMPORTANCE H7N9 subtype avian influenza viruses caused infections in over 1,400 humans from 2013 to 2017 and resulted in almost 600 deaths. It is important to understand how avian influenza viruses infect and cause disease in humans and to assess their potential for efficient person-to-person transmission. In this study, we used deep sequencing of primary clinical material to assess the evolution and potential for human adaptation of H7N9 influenza viruses.

Universal Target Capture of HIV Sequences From NGS Libraries

Background: Global surveillance of viral sequence diversity is needed to keep pace with the constant evolution of HIV. Recent next generation sequencing (NGS) methods have realized the goal of sequencing circulating virus directly from patient specimens. Yet, a simple, universal approach that maximizes sensitivity and sequencing capacity remains elusive. Here we present a novel HIV enrichment strategy to yield near complete genomes from low viral load specimens. Methodology: A non-redundant biotin-labeled probe set (HIV-xGen; n = 652) was synthesized to tile all HIV-1 (groups M, N, O, and P) and HIV-2 (A and B) strains. Illumina Nextera barcoded libraries of either gene-specific or randomly primed cDNA derived from infected plasma were hybridized to probes in a single pool and unbound sequences were washed away. Captured viral cDNA was amplified by Illumina adaptor primers, sequenced on a MiSeq, and NGS reads were demultiplexed for alignment with CLC Bio software. Results: HIV-xGen probes selectively captured and amplified reads spanning the entirety of the HIV phylogenetic tree. HIV sequences clearly present in unenriched libraries of specimens but previously not observed due to high host background levels, insufficient sequencing depth or the extent of multiplexing, were now enriched by >1,000-fold. Thus, xGen selection not only substantially increased the depth of existing sequence, but also extended overall genome coverage by an average of 40%. We characterized 50 new, diverse HIV strains from clinical specimens and demonstrated a viral load cutoff of approximately log 3.5 copies/ml for full length coverage. Genome coverage was <20% for 5/10 samples with viral loads <log 3.5 copies/ml and >90% for 35/40 samples with higher viral loads. Conclusions: Characterization of >20 complete genomes at a time is now possible from a single probe hybridization and MiSeq run. With the versatility to capture all HIV strains and the sensitivity to detect low titer specimens, HIV-xGen will serve as an important tool for monitoring HIV sequence diversity.

Design and validation of a universal influenza virus enrichment probe set and its utility in deep sequence analysis of primary cloacal swab surveillance samples of wild birds

Influenza virus infections in humans and animals are major public health concerns. In the current study, a set of universal influenza enrichment probes was developed to increase the sensitivity of sequence-based virus detection and characterization for all influenza viruses. This universal influenza enrichment probe set contains 46,953 120nt RNA biotin-labeled probes designed based on all available influenza viral sequences and it can be used to enrich for influenza sequences without prior knowledge of type or subtype. Marked enrichment was demonstrated in influenza A/H1N1, influenza B, and H1-to-H16 hemagglutinin plasmids spiked into human DNA and in cultured influenza A/H2N1 virus. Furthermore, enrichment effects and mixed influenza A virus infections were revealed in wild bird cloacal swab samples. Therefore, this universal influenza virus enrichment probe system can capture and enrich influenza viral sequences selectively and effectively in different samples, especially ones with degraded RNA or containing low amount of influenza RNA.

Metagenomics for Clinical Infectious Disease Diagnostics Steps Closer to Reality

Metagenomics approaches based on shotgun next-generation sequencing hold promise for infectious disease diagnostics. Despite substantial challenges that remain, work done over the past few years justifies excitement about the potential for these approaches to transform how clinical pathogen identification and analysis are performed. In an article in this issue of the Journal of Clinical Microbiology, M. I. Ivy et al. (J Clin Microbiol 56:e00402-18, 2018, https://doi.org/10.1128/JCM.00402-18) have applied a shotgun metagenomics approach to the diagnosis of prosthetic joint infections directly from synovial fluid. The results from this work demonstrate both the potentials and challenges of this approach applied in the clinical microbiology laboratory.

Investigation of the Plasma Virome from Cases of Unexplained Febrile Illness in Tanzania from 2013 to 2014: a Comparative Analysis between Unbiased and VirCapSeq-VERT High-Throughput Sequencing Approaches

Characterization of the viruses found in the blood of febrile patients provides information pertinent to public health and diagnostic medicine. PCR and culture have historically played an important role in clinical microbiology; however, these methods require a targeted approach and may lack the capacity to identify novel or mixed viral infections. High-throughput sequencing can overcome these constraints. As the cost of running multiple samples continues to decrease, the implementation of high-throughput sequencing for diagnostic purposes is becoming more feasible. Here we present a comparative analysis of findings from an investigation of unexplained febrile illness using two strategies: unbiased high-throughput sequencing and VirCapSeq-VERT, a positive selection high-throughput sequencing system.

High-throughput sequencing can provide insights into epidemiology and medicine through comprehensive surveys of viral genetic sequences in environmental and clinical samples. Here, we characterize the plasma virome of Tanzanian patients with unexplained febrile illness by using two high-throughput sequencing methods: unbiased sequencing and VirCapSeq-VERT (a positive selection system). Sequences from dengue virus 2, West Nile virus, human immunodeficiency virus type 1, human pegivirus, and Epstein-Barr virus were identified in plasma. Both sequencing strategies recovered nearly complete genomes in samples containing multiple viruses. Whereas VirCapSeq-VERT had better sensitivity, unbiased sequencing provided better coverage of genome termini. Together, these data demonstrate the utility of high-throughput sequencing strategies in outbreak investigations.

IMPORTANCE Characterization of the viruses found in the blood of febrile patients provides information pertinent to public health and diagnostic medicine. PCR and culture have historically played an important role in clinical microbiology; however, these methods require a targeted approach and may lack the capacity to identify novel or mixed viral infections. High-throughput sequencing can overcome these constraints. As the cost of running multiple samples continues to decrease, the implementation of high-throughput sequencing for diagnostic purposes is becoming more feasible. Here we present a comparative analysis of findings from an investigation of unexplained febrile illness using two strategies: unbiased high-throughput sequencing and VirCapSeq-VERT, a positive selection high-throughput sequencing system.

The Study of Viral RNA Diversity in Bird Samples Using De Novo Designed Multiplex Genus-Specific Primer Panels

Advances in the next generation sequencing (NGS) technologies have significantly increased our ability to detect new viral pathogens and systematically determine the spectrum of viruses prevalent in various biological samples. In addition, this approach has also helped in establishing the associations of viromes with many diseases. However, unlike the metagenomic studies using 16S rRNA for the detection of bacteria, it is impossible to create universal oligonucleotides to target all known and novel viruses, owing to their genomic diversity and variability. On the other hand, sequencing the entire genome is still expensive and has relatively low sensitivity for such applications. The existing approaches for the design of oligonucleotides for targeted enrichment are usually involved in the development of primers for the PCR-based detection of particular viral species or genera, but not for families or higher taxonomic orders. In this study, we have developed a computational pipeline for designing the oligonucleotides capable of covering a significant number of known viruses within various taxonomic orders, as well as their novel variants. We have subsequently designed a genus-specific oligonucleotide panel for targeted enrichment of viral nucleic acids in biological material and demonstrated the possibility of its application for virus detection in bird samples. We have tested our panel using a number of collected samples and have observed superior efficiency in the detection and identification of viral pathogens. Since a reliable, bioinformatics-based analytical method for the rapid identification of the sequences was crucial, an NGS-based data analysis module was developed in this study, and its functionality in the detection of novel viruses and analysis of virome diversity was demonstrated.

Ultrasensitive Capture of Human Herpes Simplex Virus Genomes Directly from Clinical Samples Reveals Extraordinarily Limited Evolution in Cell Culture

Herpes simplex viruses (HSVs) are difficult to sequence due to their large DNA genome, high GC content, and the presence of repeats. To date, most HSV genomes have been recovered from culture isolates, raising concern that these genomes may not accurately represent circulating clinical strains. We report the development and validation of a DNA oligonucleotide hybridization panel to recover nearly complete HSV genomes at abundances up to 50,000-fold lower than previously reported. Using copy number information on herpesvirus and host DNA background via quantitative PCR, we developed a protocol for pooling for cost-effective recovery of more than 50 HSV-1 or HSV-2 genomes per MiSeq run. We demonstrate the ability to recover >99% of the HSV genome at >100× coverage in 72 h at viral loads that allow whole-genome recovery from latently infected ganglia. We also report a new computational pipeline for rapid HSV genome assembly and annotation. Using the above tools and a series of 17 HSV-1-positive clinical swabs sent to our laboratory for viral isolation, we show limited evolution of HSV-1 during viral isolation in human fibroblast cells compared to the original clinical samples. Our data indicate that previous studies using low-passage-number clinical isolates of herpes simplex viruses are reflective of the viral sequences present in the lesion and thus can be used in phylogenetic analyses. We also detect superinfection within a single sample with unrelated HSV-1 strains recovered from separate oral lesions in an immunosuppressed patient during a 2.5-week period, illustrating the power of direct-from-specimen sequencing of HSV.IMPORTANCE Herpes simplex viruses affect more than 4 billion people across the globe, constituting a large burden of disease. Understanding the global diversity of herpes simplex viruses is important for diagnostics and therapeutics as well as cure research and tracking transmission among humans. To date, most HSV genomics has been performed on culture isolates and DNA swabs with high quantities of virus. We describe the development of wet-lab and computational tools that enable the accurate sequencing of near-complete genomes of HSV-1 and HSV-2 directly from clinical specimens at abundances >50,000-fold lower than previously sequenced and at significantly reduced cost. We use these tools to profile circulating HSV-1 strains in the community and illustrate limited changes to the viral genome during the viral isolation process. These techniques enable cost-effective, rapid sequencing of HSV-1 and HSV-2 genomes that will help enable improved detection, surveillance, and control of this human pathogen.

Genomic Epidemiology Reconstructs the Introduction and Spread of Zika Virus in Central America and Mexico

The Zika virus (ZIKV) epidemic in the Americas established ZIKV as a major public health threat and uncovered its association with severe diseases, including microcephaly. However, genetic epidemiology in some at-risk regions, particularly Central America and Mexico, remains limited. We report 61 ZIKV genomes from this region, generated using metagenomic sequencing with ZIKV-specific enrichment, and combine phylogenetic, epidemiological, and environmental data to reconstruct ZIKV transmission. These analyses revealed multiple independent ZIKV introductions to Central America and Mexico. One introduction, likely from Brazil via Honduras, led to most infections and the undetected spread of ZIKV through the region from late 2014. Multiple lines of evidence indicate biannual peaks of ZIKV transmission in the region, likely driven by varying local environmental conditions for mosquito vectors and herd immunity. The spatial and temporal heterogeneity of ZIKV transmission in Central America and Mexico challenges arbovirus surveillance and disease control measures.

Targeted Sequencing of Respiratory Viruses in Clinical Specimens for Pathogen Identification and Genome-Wide Analysis

A large number of viruses can individually and concurrently cause various respiratory illnesses. Metagenomic sequencing using next-generation sequencing (NGS) technology is capable of identifying a variety of pathogens. Here, we describe a method using a large panel of oligo probes to enrich sequence targets of 34 respiratory DNA and RNA viruses that reduces non-viral reads in NGS data and achieves high performance of sequencing-based pathogen identification. The approach can be applied to total nucleic acids purified from respiratory swabs stored in viral transport medium. Illumina TruSeq RNA Access Library procedure is used in targeted sequencing of respiratory viruses. The samples are subjected to RNA fragmentation, random reverse transcription, random PCR amplification, and ligation with barcoded library adaptors. The libraries are pooled and subjected to two rounds of enrichments by using a large panel of oligos designed to capture whole genomes of 34 respiratory viruses. The enriched libraries are amplified and sequenced using Illumina MiSeq sequencing system and reagents. This method can achieve viral detection sensitivity comparable with molecular assay and obtain partial to complete genome sequences for each virus to allow accurate genotyping and variant analysis.

Comprehensive viral enrichment enables sensitive respiratory virus genomic identification and analysis by next generation sequencing

Next generation sequencing (NGS) technologies have revolutionized the genomics field and are becoming more commonplace for identification of human infectious diseases. However, due to the low abundance of viral nucleic acids (NAs) in relation to host, viral identification using direct NGS technologies often lacks sufficient sensitivity. Here, we describe an approach based on two complementary enrichment strategies that significantly improves the sensitivity of NGS-based virus identification. To start, we developed two sets of DNA probes to enrich virus NAs associated with respiratory diseases. The first set of probes spans the genomes, allowing for identification of known viruses and full genome sequencing, while the second set targets regions conserved among viral families or genera, providing the ability to detect both known and potentially novel members of those virus groups. Efficiency of enrichment was assessed by NGS testing reference virus and clinical samples with known infection. We show significant improvement in viral identification using enriched NGS compared to unenriched NGS. Without enrichment, we observed an average of 0.3% targeted viral reads per sample. However, after enrichment, 50%–99% of the reads per sample were the targeted viral reads for both the reference isolates and clinical specimens using both probe sets. Importantly, dramatic improvements on genome coverage were also observed following virus-specific probe enrichment. The methods described here provide improved sensitivity for virus identification by NGS, allowing for a more comprehensive analysis of disease etiology.

Complete assembly of a dengue virus type 3 genome from a recent genotype III clade by metagenomic sequencing of serum

Background: Mosquito-borne flaviviruses, such as dengue and Japanese encephalitis virus (JEV), cause life-threatening diseases, particularly in the tropics. Methods: Here we performed unbiased metagenomic sequencing of RNA extracted from the serum of four patients and the plasma of one patient, all hospitalized at a tertiary care centre in South India with severe or prolonged febrile illness, together with the serum from one healthy control, in 2014. Results: We identified and assembled a complete dengue virus type 3 sequence from a case of severe dengue fever. We also identified a small number of JEV sequences in the serum of two adults with febrile illness, including one with severe dengue. Phylogenetic analysis revealed that the dengue sequence belonged to genotype III. It has an estimated divergence time of 13.86 years from the most highly related Indian strains. In total, 11 amino acid substitutions were predicted for this strain in the antigenic envelope protein, when compared to the parent strain used for development of the first commercial dengue vaccine.  Conclusions: We demonstrate that both genome assembly and detection of a low number of viral sequences are possible through the unbiased sequencing of clinical material. These methods may help ascertain causal agents for febrile illnesses with no known cause.

Transmission, Evolution, and Endogenization: Lessons Learned from Recent Retroviral Invasions

Viruses of the subfamily Orthoretrovirinae are defined by the ability to reverse transcribe an RNA genome into DNA that integrates into the host cell genome during the intracellular virus life cycle. Exogenous retroviruses (XRVs) are horizontally transmitted between host individuals, with disease outcome depending on interactions between the retrovirus and the host organism. When retroviruses infect germ line cells of the host, they may become endogenous retroviruses (ERVs), which are permanent elements in the host germ line that are subject to vertical transmission. These ERVs sometimes remain infectious and can themselves give rise to XRVs. This review integrates recent developments in the phylogenetic classification of retroviruses and the identification of retroviral receptors to elucidate the origins and evolution of XRVs and ERVs. We consider whether ERVs may recurrently pressure XRVs to shift receptor usage to sidestep ERV interference. We discuss how related retroviruses undergo alternative fates in different host lineages after endogenization, with koala retrovirus (KoRV) receiving notable interest as a recent invader of its host germ line. KoRV is heritable but also infectious, which provides insights into the early stages of germ line invasions as well as XRV generation from ERVs. The relationship of KoRV to primate and other retroviruses is placed in the context of host biogeography and the potential role of bats and rodents as vectors for interspecies viral transmission. Combining studies of extant XRVs and "fossil" endogenous retroviruses in koalas and other Australasian species has broadened our understanding of the evolution of retroviruses and host-retrovirus interactions.

ViroFind: A novel target-enrichment deep-sequencing platform reveals a complex JC virus population in the brain of PML patients

Deep nucleotide sequencing enables the unbiased, broad-spectrum detection of viruses in clinical samples without requiring an a priori hypothesis for the source of infection. However, its use in clinical research applications is limited by low cost-effectiveness given that most of the sequencing information from clinical samples is related to the human genome, which renders the analysis of viral genomes challenging. To overcome this limitation we developed ViroFind, an in-solution target-enrichment platform for virus detection and discovery in clinical samples. ViroFind comprises 165,433 viral probes that cover the genomes of 535 selected DNA and RNA viruses that infect humans or could cause zoonosis. The ViroFind probes are used in a hybridization reaction to enrich viral sequences and therefore enhance the detection of viral genomes via deep sequencing. We used ViroFind to detect and analyze all viral populations in the brain of 5 patients with progressive multifocal leukoencephalopathy (PML) and of 18 control subjects with no known neurological disease. Compared to direct deep sequencing, by using ViroFind we enriched viral sequences present in the clinical samples up to 127-fold. We discovered highly complex polyoma virus JC populations in the PML brain samples with a remarkable degree of genetic divergence among the JC virus variants of each PML brain sample. Specifically for the viral capsid protein VP1 gene, we identified 24 single nucleotide substitutions, 12 of which were associated with amino acid changes. The most frequent (4 of 5 samples, 80%) amino acid change was D66H, which is associated with enhanced tissue tropism, and hence likely a viral fitness advantage, compared to other variants. Lastly, we also detected sparse JC virus sequences in 10 of 18 (55.5%) of control samples and sparse human herpes virus 6B (HHV6B) sequences in the brain of 11 of 18 (61.1%) control subjects. In sum, ViroFind enabled the in-depth analysis of all viral genomes in PML and control brain samples and allowed us to demonstrate a high degree of JC virus genetic divergence in vivo that has been previously underappreciated. ViroFind can be used to investigate the structure of the virome with unprecedented depth in health and disease state.