Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples

A hybridization target enrichment approach for pathogen genomics

IMPORTANCE

Emerging infectious diseases require continuous pathogen monitoring. Rapid clinical diagnosis by nucleic acid amplification is limited to a small number of targets and may miss target detection due to new mutations in clinical isolates. Whole-genome sequencing (WGS) identifies genome-wide variations that may be used to determine a pathogen's drug resistance patterns and phylogenetically characterize isolates to track disease origin and transmission. WGS is typically performed using DNA isolated from cultured clinical isolates. Culturing clinical specimens increases turn-around time and may not be possible for fastidious bacteria. To overcome some of these limitations, direct sequencing of clinical specimens has been attempted using expensive capture probes to enrich the entire genomes of target pathogens. We present a method to produce a cost-effective, time-efficient, and large-scale synthesis of probes for whole-genome enrichment. We envision that our method can be used for direct clinical sequencing of a wide range of microbial pathogens for genomic epidemiology.

Revisiting Rustrela Virus: New Cases of Encephalitis and a Solution to the Capsid Enigma

Rustrela virus (RusV; species Rubivirus strelense) is a recently discovered relative of rubella virus (RuV) that has been detected in cases of encephalitis in diverse mammals. Here, we diagnosed two additional cases of fatal RusV-associated meningoencephalitis in a South American coati (Nasua nasua) and a Eurasian or European otter (Lutra lutra) that were detected in a zoological garden with history of prior RusV infections. Both animals showed abnormal movement or unusual behavior and their brains tested positive for RusV using specific reverse transcription quantitative PCR (RT-qPCR) and RNA in situ hybridization. As previous sequencing of the RusV genome proved to be very challenging, we employed a sophisticated target-specific capture enrichment with specifically designed RNA baits to generate complete RusV genome sequences from both detected encephalitic animals and apparently healthy wild yellow-necked field mice (Apodemus flavicollis). Furthermore, the technique was used to revise three previously published RusV genomes from two encephalitic animals and a wild yellow-necked field mouse. When comparing the newly generated RusV sequences to the previously published RusV genomes, we identified a previously undetected stretch of 309 nucleotides predicted to represent the intergenic region and the sequence encoding the N terminus of the capsid protein. This indicated that the original RusV sequence was likely incomplete due to misassembly of the genome at a region with an exceptionally high G+C content of >80 mol%. The new sequence data indicate that RusV has an overall genome length of 9,631 nucleotides with the longest intergenic region (290 nucleotides) and capsid protein-encoding sequence (331 codons) within the genus Rubivirus.

IMPORTANCE The detection of rustrela virus (RusV)-associated encephalitis in two carnivoran mammal species further extends the knowledge on susceptible species. Furthermore, we provide clinical and pathological data for the two new RusV cases, which were until now limited to the initial description of this fatal encephalitis. Using a sophisticated enrichment method prior to sequencing of the viral genome, we markedly improved the virus-to-background sequence ratio compared to that of standard procedures. Consequently, we were able to resolve and update the intergenic region and the coding region for the N terminus of the capsid protein of the initial RusV genome sequence. The updated putative capsid protein now resembles those of rubella and ruhugu virus in size and harbors a predicted RNA-binding domain that had not been identified in the initial RusV genome version. The newly determined complete RusV genomes strongly improve our knowledge of the genome structure of this novel rubivirus.

Decoding Transcription Regulatory Mechanisms Associated with Coccidioides immitis Phase Transition Using Total RNA

New or emerging infectious diseases are commonly caused by pathogens that cannot be readily manipulated or studied under common laboratory conditions. These limitations hinder standard experimental approaches and our abilities to define the fundamental molecular mechanisms underlying pathogenesis. The advance of capped small RNA sequencing (csRNA-seq) now enables genome-wide mapping of actively initiated transcripts from genes and other regulatory transcribed start regions (TSRs) such as enhancers at a precise moment from total RNA. As RNA is nonpathogenic and can be readily isolated from inactivated infectious samples, csRNA-seq can detect acute changes in gene regulation within or in response to a pathogen with remarkable sensitivity under common laboratory conditions. Studying valley fever (coccidioidomycosis), an emerging endemic fungal infection that increasingly impacts livestock, pet, and human health, we show how csRNA-seq can unravel transcriptional programs driving pathogenesis. Performing csRNA-seq on RNA isolated from different stages of the valley fever pathogen Coccidioides immitis revealed alternative promoter usage, connected cis-regulatory domains, and a WOPR family transcription factor, which are known regulators of virulence in other fungi, as being critical for pathogenic growth. We further demonstrate that a C. immitis WOPR homologue, CIMG_02671, activates transcription in a WOPR motif-dependent manner. Collectively, these findings provide novel insights into valley fever pathogenesis and provide a proof of principle for csRNA-seq as a powerful means to determine the genes, regulatory mechanisms, and transcription factors that control the pathogenesis of highly infectious agents. IMPORTANCE Infectious pathogens like airborne viruses or fungal spores are difficult to study; they require high-containment facilities, special equipment, and expertise. As such, establishing approaches such as genome editing or other means to identify the factors and mechanisms underlying caused diseases, and, thus, promising drug targets, is costly and time-intensive. These obstacles particularly hinder the analysis of new, emerging, or rare infectious diseases. We recently developed a method termed capped small RNA sequencing (csRNA-seq) that enables capturing acute changes in active gene expression from total RNA. Prior to csRNA-seq, such an analysis was possible only by using living cells or nuclei, in which pathogens are highly infectious. The process of RNA purification, however, inactivates pathogens and thus enables the analysis of gene expression during disease progression under standard laboratory conditions. As a proof of principle, here, we use csRNA-seq to unravel the gene regulatory programs and factors likely critical for the pathogenesis of valley fever, an emerging endemic fungal infection that increasingly impacts livestock, pet, and human health.

A novel SARS-CoV-2 related coronavirus in bats from Cambodia

Knowledge of the origin and reservoir of the coronavirus responsible for the ongoing COVID-19 pandemic is still fragmentary. To date, the closest relatives to SARS-CoV-2 have been detected in Rhinolophus bats sampled in the Yunnan province, China. Here we describe the identification of SARS-CoV-2 related coronaviruses in two Rhinolophus shameli bats sampled in Cambodia in 2010. Metagenomic sequencing identifies nearly identical viruses sharing 92.6% nucleotide identity with SARS-CoV-2. Most genomic regions are closely related to SARS-CoV-2, with the exception of a region of the spike, which is not compatible with human ACE2-mediated entry. The discovery of these viruses in a bat species not found in China indicates that SARS-CoV-2 related viruses have a much wider geographic distribution than previously reported, and suggests that Southeast Asia represents a key area to consider for future surveillance for coronaviruses.

In this study, Delaune et al., isolate and characterise a SARS-CoV-2-related coronavirus from two bats sampled in Cambodia. Their findings suggest that the geographic distribution of SARS-CoV-2-related viruses is wider than previously reported.

Genomic surveillance of enterovirus associated with aseptic meningitis cases in southern Spain, 2015-2018

New circulating Enterovirus (EV) strains often emerge through recombination. Upsurges of recombinant non-polio enteroviruses (NPEVs) associated with neurologic manifestations such as EVA71 or Echovirus 30 (E30) are a growing public health concern in Europe. Only a few complete genomes of EVs circulating in Spain are available in public databases, making it difficult to address the emergence of recombinant EVs, understand their evolutionary relatedness and the possible implication in human disease. We have used metagenomic (untargeted) NGS to generate full-length EV genomes from CSF samples of EV-positive aseptic meningitis cases in Southern Spain between 2015 and 2018. Our analyses reveal the co-circulation of multiple Enterovirus B (EV-B) types (E6, E11, E13 and E30), including a novel E13 recombinant form. We observed a genetic turnover where emergent lineages (C1 for E6 and I [tentatively proposed in this study] for E30) replaced previous lineages circulating in Spain, some concomitant with outbreaks in other parts of Europe. Metagenomic sequencing provides an effective approach for the analysis of EV genomes directly from PCR-positive CSF samples. The detection of a novel, disease-associated, recombinant form emphasizes the importance of genomic surveillance to monitor spread and evolution of EVs.

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.

Enhanced Virus Detection and Metagenomic Sequencing in Patients with Meningitis and Encephalitis

Meningitis and encephalitis are leading causes of central nervous system (CNS) disease and often result in severe neurological compromise or death. Traditional diagnostic workflows largely rely on pathogen-specific tests, sometimes over days to weeks, whereas metagenomic next-generation sequencing (mNGS) profiles all nucleic acid in a sample. In this single-center, prospective study, 68 hospitalized patients with known (n = 44) or suspected (n = 24) CNS infections underwent mNGS from RNA and DNA to identify potential pathogens and also targeted sequencing of viruses using hybrid capture. Using a computational metagenomic classification pipeline based on KrakenUniq and BLAST, we detected pathogen nucleic acid in cerebrospinal fluid (CSF) from 22 subjects, 3 of whom had no clinical diagnosis by routine workup. Among subjects diagnosed with infection by serology and/or peripheral samples, we demonstrated the utility of mNGS to detect pathogen nucleic acid in CSF, importantly for the Ixodes scapularis tick-borne pathogens Powassan virus, Borrelia burgdorferi, and Anaplasma phagocytophilum. We also evaluated two methods to enhance the detection of viral nucleic acid, hybrid capture and methylated DNA depletion. Hybrid capture nearly universally increased viral read recovery. Although results for methylated DNA depletion were mixed, it allowed the detection of varicella-zoster virus DNA in two samples that were negative by standard mNGS. Overall, mNGS is a promising approach that can test for multiple pathogens simultaneously, with efficacy similar to that of pathogen-specific tests, and can uncover geographically relevant infectious CNS disease, such as tick-borne infections in New England. With further laboratory and computational enhancements, mNGS may become a mainstay of workup for encephalitis and meningitis. IMPORTANCE Meningitis and encephalitis are leading global causes of central nervous system (CNS) disability and mortality. Current diagnostic workflows remain inefficient, requiring costly pathogen-specific assays and sometimes invasive surgical procedures. Despite intensive diagnostic efforts, 40 to 60% of people with meningitis or encephalitis have no clear cause of CNS disease identified. As diagnostic uncertainty often leads to costly inappropriate therapies, the need for novel pathogen detection methods is paramount. Metagenomic next-generation sequencing (mNGS) offers the unique opportunity to circumvent these challenges using unbiased laboratory and computational methods. Here, we performed comprehensive mNGS from 68 prospectively enrolled patients with known (n = 44) or suspected (n = 24) CNS viral infection from a single center in New England and evaluated enhanced methods to improve the detection of CNS pathogens, including those not traditionally identified in the CNS by nucleic acid detection. Overall, our work helps elucidate how mNGS can become integrated into the diagnostic toolkit for CNS infections.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

A streamlined clinical metagenomic sequencing protocol for rapid pathogen identification

Metagenomic next-generation sequencing (mNGS) holds promise as a diagnostic tool for unbiased pathogen identification and precision medicine. However, its medical utility depends largely on assay simplicity and reproducibility. In the current study, we aimed to develop a streamlined Illumina and Oxford Nanopore-based DNA/RNA library preparation protocol and rapid data analysis pipeline. The Illumina sequencing-based mNGS method was first developed and evaluated using a set of samples with known aetiology. Its sensitivity for RNA viruses (influenza A, H1N1) was < 6.4 × 10² EID50/mL, and a good correlation between viral loads and mapped reads was observed. Then, the rapid turnaround time of Nanopore sequencing was tested by sequencing influenza A virus and adenoviruses. Furthermore, 11 respiratory swabs or sputum samples pre-tested for a panel of pathogens were analysed, and the pathogens identified by Illumina sequencing showed 81.8% concordance with qPCR results. Additional sequencing of cerebrospinal fluid (CSF) samples from HIV-1-positive patients with meningitis/encephalitis detected HIV-1 RNA and Toxoplasma gondii sequences. In conclusion, we have developed a simplified protocol that realizes efficient metagenomic sequencing of a variety of clinical samples and pathogen identification in a clinically meaningful time frame.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Genomic investigation of atypical malaria cases in Kanel, northern Senegal

BACKGROUND

The diagnosis of malaria cases in regions where the malaria burden has decreased significantly and prevalence is very low is more challenging, in part because of reduced clinical presumption of malaria. The appearance of a cluster of malaria cases with atypical symptoms in Mbounguiel, a village in northern Senegal where malaria transmission is low, in September 2018 exemplifies this scenario. The collaboration between the National Malaria Control Programme (NMCP) at the Senegal Ministry of Health and the Laboratory of Parasitology and Mycology at Cheikh Anta Diop University worked together to evaluate this cluster of malaria cases using molecular and serological tools.

METHODS

Malaria cases were diagnosed primarily by rapid diagnostic test (RDT), and confirmed by photo-induced electron transfer-polymerase chain reaction (PET-PCR). 24 single nucleotide polymorphisms (SNPs) barcoding was used for Plasmodium falciparum genotyping. Unbiased metagenomic sequencing and Luminex-based multi-pathogen antibody and antigen profiling were used to assess exposure to other pathogens.

RESULTS

Nine patients, of 15 suspected cases, were evaluated, and all nine samples were found to be positive for P. falciparum only. The 24 SNPs molecular barcode showed the predominance of polygenomic infections, with identifiable strains being different from one another. All patients tested positive for the P. falciparum antigens. No other pathogenic infection was detected by either the serological panel or metagenomic sequencing.

CONCLUSIONS

This work, undertaken locally within Senegal as a collaboration between the NMCP and a research laboratory at University of Cheikh Anta Diop (UCAD) revealed that a cluster of malaria cases were caused by different strains of P. falciparum. The public health response in real time demonstrates the value of local molecular and genomics capacity in affected countries for disease control and elimination.

Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events

Analysis of 772 complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from early in the Boston-area epidemic revealed numerous introductions of the virus, a small number of which led to most cases. The data revealed two superspreading events. One, in a skilled nursing facility, led to rapid transmission and significant mortality in this vulnerable population but little broader spread, whereas other introductions into the facility had little effect. The second, at an international business conference, produced sustained community transmission and was exported, resulting in extensive regional, national, and international spread. The two events also differed substantially in the genetic variation they generated, suggesting varying transmission dynamics in superspreading events. Our results show how genomic epidemiology can help to understand the link between individual clusters and wider community spread.

Genomic surveillance of Nevada patients revealed prevalence of unique SARS-CoV-2 variants bearing mutations in the RdRp gene

Patients with signs of COVID-19 were tested through diagnostic RT-PCR for SARS-CoV-2 using RNA extracted from the nasopharyngeal/nasal swabs. To determine the variants of SARS-CoV-2 circulating in the state of Nevada, specimens from 200 COVID-19 patients were sequenced through our robust sequencing platform, which enabled sequencing of SARS-CoV-2 from specimens with even very low viral loads, without the need of culture-based amplification. High genome coverage allowed the identification of single and multi-nucleotide variants in SARS-CoV-2 in the community and their phylogenetic relationships with other variants present during the same period of the outbreak. We report the occurrence of a novel mutation at 323aa (314aa of orf1b) of nsp12 (RNA-dependent RNA polymerase) changed to phenylalanine (F) from proline (P), in the first reported isolate of SARS-CoV-2, Wuhan-Hu-1. This 323F variant was present at a very high frequency in Northern Nevada. Structural modeling determined this mutation in the interface domain, which is important for the association of accessory proteins required for the polymerase. In conclusion, we report the introduction of specific SARS-CoV-2 variants at very high frequency in distinct geographic locations, which is important for understanding the evolution and circulation of SARS-CoV-2 variants of public health importance, while it circulates in humans.

Molecular epidemiology of Plasmodium falciparum by multiplexed amplicon deep sequencing in Senegal

BACKGROUND

Molecular epidemiology can provide important information regarding the genetic diversity and transmission of Plasmodium falciparum, which can assist in designing and monitoring elimination efforts. However, malaria molecular epidemiology including understanding the genetic diversity of the parasite and performing molecular surveillance of transmission has been poorly documented in Senegal. Next Generation Sequencing (NGS) offers a practical, fast and high-throughput approach to understand malaria population genetics. This study aims to unravel the population structure of P. falciparum and to estimate the allelic diversity, multiplicity of infection (MOI), and evolutionary patterns of the malaria parasite using the NGS platform.

METHODS

Multiplex amplicon deep sequencing of merozoite surface protein 1 (PfMSP1) and merozoite surface protein 2 (PfMSP2) in fifty-three P. falciparum isolates from two epidemiologically different areas in the South and North of Senegal, was carried out.

RESULTS

A total of 76 Pfmsp1 and 116 Pfmsp2 clones were identified and 135 different alleles were found, 56 and 79 belonged to the pfmsp1 and pfmsp2 genes, respectively. K1 and IC3D7 allelic families were most predominant in both sites. The local haplotype diversity (Hd) and nucleotide diversity (π) were higher in the South than in the North for both genes. For pfmsp1, a high positive Tajima's D (TD) value was observed in the South (D = 2.0453) while negative TD value was recorded in the North (D = - 1.46045) and F-Statistic (Fst) was 0.19505. For pfmsp2, non-directional selection was found with a highly positive TD test in both areas and Fst was 0.02111. The mean MOI for both genes was 3.07 and 1.76 for the South and the North, respectively, with a statistically significant difference between areas (p = 0.001).

CONCLUSION

This study revealed a high genetic diversity of pfmsp1 and pfmsp2 genes and low genetic differentiation in P. falciparum population in Senegal. The MOI means were significantly different between the Southern and Northern areas. Findings also showed that multiplexed amplicon deep sequencing is a useful technique to investigate genetic diversity and molecular epidemiology of P. falciparum infections.

Powassan Virus Neuropathology and Genomic Diversity in Patients With Fatal Encephalitis

Background

Powassan virus (POWV) is an emerging cause of severe encephalitis; very little is known about human pathogenicity due to challenges in diagnosis and viral RNA recovery. We present 3 patients with fatal encephalitis due to POWV lineage II (deer tick virus).

Methods

We obtained 27 unique samples, including from brain biopsy and autopsy, and used metagenomic sequencing, quantitative reverse transcriptase polymerase chain reaction, and a newly developed CRISPR-based diagnostic assay to perform the first detailed characterization of POWV compartmentalization and genomics between and within human subjects.

Results

In all 3 patients, imaging and histopathology findings were notable for profound cerebellar involvement. All patients were initially diagnosed with POWV by metagenomic sequencing, and 2 of the 3 had negative clinical testing by serology. We detected POWV RNA in 13 clinical samples; levels were highest in the cerebellum, and there was very little involvement of peripheral tissue. We assembled complete POWV genomes from 8 samples, providing unique information about the strains of POWV lineage II (deer tick virus) that infect humans.

Conclusions

We demonstrate the utility of molecular assays for detecting POWV infection, including in seronegative patients, and nominate viral genomic features that may relate to human infection and neuropathogenicity. The cerebellum was identified as a key target POWV in fatal infection, by radiological and histopathological findings as well as molecular testing.

Phylogenetic analysis of SARS-CoV-2 in the Boston area highlights the role of recurrent importation and superspreading events

SARS-CoV-2 has caused a severe, ongoing outbreak of COVID-19 in Massachusetts with 111,070 confirmed cases and 8,433 deaths as of August 1, 2020. To investigate the introduction, spread, and epidemiology of COVID-19 in the Boston area, we sequenced and analyzed 772 complete SARS-CoV-2 genomes from the region, including nearly all confirmed cases within the first week of the epidemic and hundreds of cases from major outbreaks at a conference, a nursing facility, and among homeless shelter guests and staff. The data reveal over 80 introductions into the Boston area, predominantly from elsewhere in the United States and Europe. We studied two superspreading events covered by the data, events that led to very different outcomes because of the timing and populations involved. One produced rapid spread in a vulnerable population but little onward transmission, while the other was a major contributor to sustained community transmission, including outbreaks in homeless populations, and was exported to several other domestic and international sites. The same two events differed significantly in the number of new mutations seen, raising the possibility that SARS-CoV-2 superspreading might encompass disparate transmission dynamics. Our results highlight the failure of measures to prevent importation into MA early in the outbreak, underscore the role of superspreading in amplifying an outbreak in a major urban area, and lay a foundation for contact tracing informed by genetic data.

RNAseq Analysis of Rhizomania-Infected Sugar Beet Provides the First Genome Sequence of Beet Necrotic Yellow Vein Virus from the USA and Identifies a Novel Alphanecrovirus and Putative Satellite Viruses

“Rhizomania” of sugar beet is a soilborne disease complex comprised of beet necrotic yellow vein virus (BNYVV) and its plasmodiophorid vector, Polymyxa betae. Although BNYVV is considered the causal agent of rhizomania, additional viruses frequently accompany BNYVV in diseased roots. In an effort to better understand the virus cohort present in sugar beet roots exhibiting rhizomania disease symptoms, five independent RNA samples prepared from diseased beet seedlings reared in a greenhouse or from field-grown adult sugar beet plants and enriched for virus particles were subjected to RNAseq. In all but a healthy control sample, the technique was successful at identifying BNYVV and provided sequence reads of sufficient quantity and overlap to assemble > 98% of the published genome of the virus. Utilizing the derived consensus sequence of BNYVV, infectious RNA was produced from cDNA clones of RNAs 1 and 2. The approach also enabled the detection of beet soilborne mosaic virus (BSBMV), beet soilborne virus (BSBV), beet black scorch virus (BBSV), and beet virus Q (BVQ), with near-complete genome assembly afforded to BSBMV and BBSV. In one field sample, a novel virus sequence of 3682 nt was assembled with significant sequence similarity and open reading frame (ORF) organization to members within the subgenus Alphanecrovirus (genus Necrovirus; family Tombusviridae). Construction of a DNA clone based on this sequence led to the production of the novel RNA genome in vitro that was capable of inducing local lesion formation on leaves of Chenopodium quinoa. Additionally, two previously unreported satellite viruses were revealed in the study; one possessing weak similarity to satellite maize white line mosaic virus and a second possessing moderate similarity to satellite tobacco necrosis virus C. Taken together, the approach provides an efficient pipeline to characterize variation in the BNYVV genome and to document the presence of other viruses potentially associated with disease severity or the ability to overcome resistance genes used for sugar beet rhizomania disease management.

Massively multiplexed nucleic acid detection with Cas13

The great majority of globally circulating pathogens go undetected, undermining patient care and hindering outbreak preparedness and response. To enable routine surveillance and comprehensive diagnostic applications, there is a need for detection technologies that can scale to test many samples1-3 while simultaneously testing for many pathogens4-6. Here, we develop Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (CARMEN), a platform for scalable, multiplexed pathogen detection. In the CARMEN platform, nanolitre droplets containing CRISPR-based nucleic acid detection reagents7 self-organize in a microwell array8 to pair with droplets of amplified samples, testing each sample against each CRISPR RNA (crRNA) in replicate. The combination of CARMEN and Cas13 detection (CARMEN-Cas13) enables robust testing of more than 4,500 crRNA-target pairs on a single array. Using CARMEN-Cas13, we developed a multiplexed assay that simultaneously differentiates all 169 human-associated viruses with at least 10 published genome sequences and rapidly incorporated an additional crRNA to detect the causative agent of the 2020 COVID-19 pandemic. CARMEN-Cas13 further enables comprehensive subtyping of influenza A strains and multiplexed identification of dozens of HIV drug-resistance mutations. The intrinsic multiplexing and throughput capabilities of CARMEN make it practical to scale, as miniaturization decreases reagent cost per test by more than 300-fold. Scalable, highly multiplexed CRISPR-based nucleic acid detection shifts diagnostic and surveillance efforts from targeted testing of high-priority samples to comprehensive testing of large sample sets, greatly benefiting patients and public health9-11.

50 Years of Lassa Fever Research

Lassa fever was first described as a clinical entity fifty years ago. The causative agent Lassa virus was isolated from these first known cases. This chapter reviews the key publications on Lassa fever research that appeared in the scientific literature at that time and over the ensuing decades.

Novel genome sequences of cell-fusing agent virus allow comparison of virus phylogeny with the genetic structure of Aedes aegypti populations

Flaviviruses encompass not only medically relevant arthropod-borne viruses (arboviruses) but also insect-specific flaviviruses (ISFs) that are presumably maintained primarily through vertical transmission in the insect host. Interestingly, ISFs are commonly found infecting important arbovirus vectors such as the mosquito Aedes aegypti. Cell-fusing agent virus (CFAV) was the first described ISF of mosquitoes more than four decades ago. Despite evidence for widespread CFAV infections in A.aegypti populations and for CFAV potential to interfere with arbovirus transmission, little is known about CFAV evolutionary history. Here, we generated six novel CFAV genome sequences by sequencing three new virus isolates and subjecting three mosquito samples to untargeted viral metagenomics. We used these new genome sequences together with published ones to perform a global phylogenetic analysis of CFAV genetic diversity. Although there was some degree of geographical clustering among CFAV sequences, there were also notable discrepancies between geography and phylogeny. In particular, CFAV sequences from Cambodia and Thailand diverged significantly, despite confirmation that A.aegypti populations from both locations are genetically close. The apparent phylogenetic discrepancy between CFAV and its A.aegypti host in Southeast Asia indicates that other factors than host population structure shape CFAV genetic diversity.

Amplicon deep sequencing of kelch13 in Plasmodium falciparum isolates from Senegal

BACKGROUND

In 2006, the Senegalese National Malaria Control Programme recommended artemisinin-based combination therapy (ACT) with artemether-lumefantrine as the first-line treatment for uncomplicated Plasmodium falciparum malaria. To date, multiple mutations associated with artemisinin delayed parasite clearance have been described in Southeast Asia in the Pfk13 gene, such as Y493H, R539T, I543T and C580Y. Even though ACT remains clinically and parasitologically efficacious in Senegal, the spread of resistance is possible as shown by the earlier emergence of resistance to chloroquine in Southeast Asia that subsequently spread to Africa. Therefore, surveillance of artemisinin resistance in malaria endemic regions is crucial and requires the implementation of sensitive tools, such as next-generation sequencing (NGS) which can detect novel mutations at low frequency.

METHODS

Here, an amplicon sequencing approach was used to identify mutations in the Pfk13 gene in eighty-one P. falciparum isolates collected from three different regions of Senegal.

RESULTS

In total, 10 SNPs around the propeller domain were identified; one synonymous SNP and nine non-synonymous SNPs, and two insertions. Three of these SNPs (T478T, A578S and V637I) were located in the propeller domain. A578S, is the most frequent mutation observed in Africa, but has not previously been reported in Senegal. A previous study has suggested that A578S could disrupt the function of the Pfk13 propeller region.

CONCLUSION

As the genetic basis of possible artemisinin resistance may be distinct in Africa and Southeast Asia, further studies are necessary to assess the new SNPs reported in this study.

Spatially distinct physiology of Bacteroides fragilis within the proximal colon of gnotobiotic mice

A complex microbiota inhabits various microenvironments of the gut, with some symbiotic bacteria having evolved traits to invade the epithelial mucus layer and reside deep within intestinal tissue of animals. Whether these distinct bacterial communities across gut biogeographies exhibit divergent behaviors remains largely unknown. Global transcriptomic analysis to investigate microbial physiology in specific mucosal niches has been hampered technically by overabundance of host RNA. Herein, we employed hybrid selection RNA sequencing (hsRNA-Seq) to enable detailed spatial transcriptomic profiling of a prominent human commensal as it colonizes the colonic lumen, mucus or epithelial tissue of mice. Compared to conventional RNA-Seq, hsRNA-Seq increased reads mapping to the Bacteroides fragilis genome by 48- and 154-fold in mucus and tissue, respectively, allowing for high fidelity comparisons across biogeographic sites. Near the epithelium, B. fragilis up-regulated numerous genes involved in protein synthesis, indicating that bacteria inhabiting the mucosal niche are metabolically active. Further, a specific sulfatase (BF3086) and glycosyl hydrolase (BF3134) were highly induced in mucus and tissue compared to bacteria in the lumen. In-frame deletion of these genes impaired in vitro growth on mucus as a carbon source, as well as mucosal colonization of mice. Mutants in either B. fragilis gene displayed a fitness defect in competing for colonization against bacterial challenge, revealing the importance of site-specific gene expression for robust host-microbial symbiosis. As a versatile tool, hsRNA-Seq can be deployed to explore the in vivo spatial physiology of numerous bacterial pathogens or commensals.

Microbiome analyses of blood and tissues suggest cancer diagnostic approach

Systematic characterization of the cancer microbiome provides the opportunity to develop techniques that exploit non-human, microorganism-derived molecules in the diagnosis of a major human disease. Following recent demonstrations that some types of cancer show substantial microbial contributions1-10, we re-examined whole-genome and whole-transcriptome sequencing studies in The Cancer Genome Atlas11 (TCGA) of 33 types of cancer from treatment-naive patients (a total of 18,116 samples) for microbial reads, and found unique microbial signatures in tissue and blood within and between most major types of cancer. These TCGA blood signatures remained predictive when applied to patients with stage Ia-IIc cancer and cancers lacking any genomic alterations currently measured on two commercial-grade cell-free tumour DNA platforms, despite the use of very stringent decontamination analyses that discarded up to 92.3% of total sequence data. In addition, we could discriminate among samples from healthy, cancer-free individuals (n = 69) and those from patients with multiple types of cancer (prostate, lung, and melanoma; 100 samples in total) solely using plasma-derived, cell-free microbial nucleic acids. This potential microbiome-based oncology diagnostic tool warrants further exploration.

Current Trends in Diagnostics of Viral Infections of Unknown Etiology

Viruses are evolving at an alarming rate, spreading and inconspicuously adapting to cutting-edge therapies. Therefore, the search for rapid, informative and reliable diagnostic methods is becoming urgent as ever. Conventional clinical tests (PCR, serology, etc.) are being continually optimized, yet provide very limited data. Could high throughput sequencing (HTS) become the future gold standard in molecular diagnostics of viral infections? Compared to conventional clinical tests, HTS is universal and more precise at profiling pathogens. Nevertheless, it has not yet been widely accepted as a diagnostic tool, owing primarily to its high cost and the complexity of sample preparation and data analysis. Those obstacles must be tackled to integrate HTS into daily clinical practice. For this, three objectives are to be achieved: (1) designing and assessing universal protocols for library preparation, (2) assembling purpose-specific pipelines, and (3) building computational infrastructure to suit the needs and financial abilities of modern healthcare centers. Data harvested with HTS could not only augment diagnostics and help to choose the correct therapy, but also facilitate research in epidemiology, genetics and virology. This information, in turn, could significantly aid clinicians in battling viral infections.

Combining genomics and epidemiology to track mumps virus transmission in the United States

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks.

Programmable Inhibition and Detection of RNA Viruses Using Cas13

The CRISPR effector Cas13 could be an effective antiviral for single-stranded RNA (ssRNA) viruses because it programmably cleaves RNAs complementary to its CRISPR RNA (crRNA). Here, we computationally identify thousands of potential Cas13 crRNA target sites in hundreds of ssRNA viral species that can potentially infect humans. We experimentally demonstrate Cas13's potent activity against three distinct ssRNA viruses: lymphocytic choriomeningitis virus (LCMV); influenza A virus (IAV); and vesicular stomatitis virus (VSV). Combining this antiviral activity with Cas13-based diagnostics, we develop Cas13-assisted restriction of viral expression and readout (CARVER), an end-to-end platform that uses Cas13 to detect and destroy viral RNA. We further screen hundreds of crRNAs along the LCMV genome to evaluate how conservation and target RNA nucleotide content influence Cas13's antiviral activity. Our results demonstrate that Cas13 can be harnessed to target a wide range of ssRNA viruses and CARVER's potential broad utility for rapid diagnostic and antiviral drug development.

dsRNA-Seq: Identification of Viral Infection by Purifying and Sequencing dsRNA

RNA viruses are a major source of emerging and re-emerging infectious diseases around the world. We developed a method to identify RNA viruses that is based on the fact that RNA viruses produce double-stranded RNA (dsRNA) while replicating. Purifying and sequencing dsRNA from the total RNA isolated from infected tissue allowed us to recover dsRNA virus sequences and replicated sequences from single-stranded RNA (ssRNA) viruses. We refer to this approach as dsRNA-Seq. By assembling dsRNA sequences into contigs we identified full length or partial RNA viral genomes of varying genome types infecting mammalian culture samples, identified a known viral disease agent in laboratory infected mice, and successfully detected naturally occurring RNA viral infections in reptiles. Here, we show that dsRNA-Seq is a preferable method for identifying viruses in organisms that don’t have sequenced genomes and/or commercially available rRNA depletion reagents. In addition, a significant advantage of this method is the ability to identify replicated viral sequences of ssRNA viruses, which is useful for distinguishing infectious viral agents from potential noninfectious viral particles or contaminants.

The utility and perspectives of NGS-based methods in BSL-3 and BSL-4 laboratory - sequencing and analysis strategies

Modern diagnostics is in general based on molecular biology methods. Nowadays sequencing-based methods, especially whole genome sequencing, are becoming increasingly important. Implementation of such methods into routine diagnostic of highly dangerous pathogens, like Bacillus anthracis, Francisella tularensis, Yersinia pestis, Ebola virus, MERS, Lassa virus etc. would be very helpful. The best diagnostic strategy would be the metagenomic sequencing directly from the clinical sample. Implementation of majority of currently available WGS platforms inside the BSL-3 or 4 laboratory is impractical because of the size of the equipment and time consuming wet lab part (e.g. library preparation). Nowadays there is a possibility to implement pocket size MinION - real time whole genome sequencer into BSL-3 and 4 laboratory for rapid and precise diagnostic purposes.

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.

Introduction of Ebola virus into a remote border district of Sierra Leone, 2014: use of field epidemiology and RNA sequencing to describe chains of transmission

In early October 2014, 7 months after the 2014-2015 Ebola epidemic in West Africa began, a cluster of reported deaths in Koinadugu, a remote district of Sierra Leone, was the first evidence of Ebola virus disease (Ebola) in the district. Prior to this event, geographic isolation was thought to have prevented the introduction of Ebola to this area. We describe our initial investigation of this cluster of deaths and subsequent public health actions after Ebola was confirmed, and present challenges to our investigation and methods of overcoming them. We present a transmission tree and results of whole genome sequencing of selected isolates to identify the source of infection in Koinadugu and demonstrate transmission between its villages. Koinadugu's experience highlights the danger of assuming that remote location and geographic isolation can prevent the spread of Ebola, but also demonstrates how deployment of rapid field response teams can help limit spread once Ebola is detected.

Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea)

Bats are unique mammals, exhibit distinctive life history traits and have unique immunological approaches to suppression of viral diseases upon infection. High-throughput next-generation sequencing has been used in characterizing the virome of different bat species. The cave nectar bat, Eonycteris spelaea, has a broad geographical range across Southeast Asia, India and southern China, however, little is known about their involvement in virus transmission. Here we investigate the diversity and abundance of viral communities from a colony of Eonycteris spelaea residing in Singapore. Our results detected 47 and 22 different virus families from bat fecal and urine samples, respectively. Among these, we identify a large number of virus families including Adenoviridae, Flaviviridae, Reoviridae, Papillomaviridae, Paramyxoviridae, Parvoviridae, Picornaviridae, and Polyomaviridae. In most cases, viral sequences from Eonycteris spelaea are genetically related to a group of bat viruses from other bat genera (e.g., Eidolon, Miniopterus, Rhinolophus and Rousettus). The results of this study improve our knowledge of the host range, spread and evolution of several important viral pathogens. More significantly, our findings provide a baseline to study the temporal patterns of virus shedding and how they correlate with bat phenological trends.

Identification and genetic characterization of a novel Orthobunyavirus species by a straightforward high-throughput sequencing-based approach

Identification and characterization of novel unknown viruses is of great importance. The introduction of high-throughput sequencing (HTS)-based methods has paved the way for genomics-based detection of pathogens without any prior assumptions about the characteristics of the organisms. However, the use of HTS for the characterization of viral pathogens from clinical samples remains limited. Here, we report the identification of a novel Orthobunyavirus species isolated from horse plasma. The identification was based on a straightforward HTS approach. Following enrichment in cell culture, RNA was extracted from the growth medium and rapid library preparation, HTS and primary bioinformatic analyses were performed in less than 12 hours. Taxonomical profiling of the sequencing reads did not reveal sequence similarities to any known virus. Subsequent application of de novo assembly tools to the sequencing reads produced contigs, of which three showed some similarity to the L, M, and S segments of viruses belonging to the Orthobunyavirus genus. Further refinement of these contigs resulted in high-quality, full-length genomic sequences of the three genomic segments (L, M and S) of a novel Orthobunyavirus. Characterization of the genomic sequence, including the prediction of open reading frames and the inspection of consensus genomic termini and phylogenetic analysis, further confirmed that the novel virus is indeed a new species, which we named Ness Ziona virus.

Marburg virus disease outbreak in Kween District Uganda, 2017: Epidemiological and laboratory findings

INTRODUCTION

In October 2017, a blood sample from a resident of Kween District, Eastern Uganda, tested positive for Marburg virus. Within 24 hour of confirmation, a rapid outbreak response was initiated. Here, we present results of epidemiological and laboratory investigations.

METHODS

A district task force was activated consisting of specialised teams to conduct case finding, case management and isolation, contact listing and follow up, sample collection and testing, and community engagement. An ecological investigation was also carried out to identify the potential source of infection. Virus isolation and Next Generation sequencing were performed to identify the strain of Marburg virus.

RESULTS

Seventy individuals (34 MVD suspected cases and 36 close contacts of confirmed cases) were epidemiologically investigated, with blood samples tested for MVD. Only four cases met the MVD case definition; one was categorized as a probable case while the other three were confirmed cases. A total of 299 contacts were identified; during follow- up, two were confirmed as MVD. Of the four confirmed and probable MVD cases, three died, yielding a case fatality rate of 75%. All four cases belonged to a single family and 50% (2/4) of the MVD cases were female. All confirmed cases had clinical symptoms of fever, vomiting, abdominal pain and bleeding from body orifices. Viral sequences indicated that the Marburg virus strain responsible for this outbreak was closely related to virus strains previously shown to be circulating in Uganda.

CONCLUSION

This outbreak of MVD occurred as a family cluster with no additional transmission outside of the four related cases. Rapid case detection, prompt laboratory testing at the Uganda National VHF Reference Laboratory and presence of pre-trained, well-prepared national and district rapid response teams facilitated the containment and control of this outbreak within one month, preventing nationwide and global transmission of the disease.

Mutations present in a low-passage Zika virus isolate result in attenuated pathogenesis in mice

Zika virus (ZIKV) infection can result in neurological disorders including Congenital Zika Syndrome in infants exposed to the virus in utero. Pregnant women can be infected by mosquito bite as well as by sexual transmission from infected men. Herein, the variants of ZIKV within the male reproductive tract and ejaculates were assessed in inoculated mice. We identified two non-synonymous variants at positions E-V330L and NS1-W98G. These variants were also present in the passage three PRVABC59 isolate and infectious clone relative to the patient serum PRVABC59 sequence. In subsequent studies, ZIKV E-330L was less pathogenic in mice than ZIKV E-330V as evident by increased average survival times. In Vero cells, ZIKV E-330L/NS1-98G outcompeted ZIKV E-330V/NS1-98W within 3 passages. These results suggest that the E-330L/NS1-98G variants are attenuating in mice and were enriched during cell culture passaging. Cell culture propagation of ZIKV could significantly affect animal model development and vaccine efficacy studies.

Precision epidemiology for infectious disease control

Advances in genomics and computing are transforming the capacity for the characterization of biological systems, and researchers are now poised for a precision-focused transformation in the way they prepare for, and respond to, infectious diseases. This includes the use of genome-based approaches to inform molecular diagnosis and individual-level treatment regimens. In addition, advances in the speed and granularity of pathogen genome generation have improved the capability to track and understand pathogen transmission, leading to potential improvements in the design and implementation of population-level public health interventions. In this Perspective, we outline several trends that are driving the development of precision epidemiology of infectious disease and their implications for scientists' ability to respond to outbreaks.

Diagnostics for Nipah virus: a zoonotic pathogen endemic to Southeast Asia

Nipah virus (NiV) is an emerging pathogen that, unlike other priority pathogens identified by WHO, is endemic to Southeast Asia. It is most commonly transmitted through exposure to saliva or excrement from the Pteropus fruit bat, or direct contact with intermediate animal hosts, such as pigs. NiV infection causes severe febrile encephalitic disease and/or respiratory disease; treatment options are limited to supportive care. A number of in-house diagnostic assays for NiV using serological and nucleic acid amplification techniques have been developed for NiV and are used in laboratory settings, including some early multiplex panels for differentiation of NiV infection from other febrile diseases. However, given the often rural and remote nature of NiV outbreak settings, there remains a need for rapid diagnostic tests that can be implemented at the point of care. Additionally, more reliable assays for surveillance of communities and livestock will be vital to achieving a better understanding of the ecology of the fruit bat host and transmission risk to other intermediate hosts, enabling implementation of a ‘One Health’ approach to outbreak prevention and the management of this zoonotic disease. An improved understanding of NiV viral diversity and infection kinetics or dynamics will be central to the development of new diagnostics, and access to clinical specimens must be improved to enable effective validation and external quality assessments. Target product profiles for NiV diagnostics should be refined to take into account these outstanding needs.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

Metagenomic Sequencing of HIV-1 in the Blood and Female Genital Tract Reveals Little Quasispecies Diversity during Acute Infection

Due to error-prone replication, HIV-1 generates a diverse population of viruses within a chronically infected individual. When HIV-1 is transmitted to a new individual, one or a few viruses establish the new infection, leading to a genetic bottleneck in the virus population. Understanding the timing and nature of this bottleneck may provide insight into HIV-1 vaccine design and other preventative strategies. We examined the HIV-1 population in three women enrolled in a unique prospective cohort in South Africa who were followed closely during the earliest stages of HIV-1 infection. We found very little HIV-1 diversity in the blood and female genital tract during the first 2 weeks after virus was detected in the bloodstream. These results are compatible with a very early HIV-1 population bottleneck, suggesting the need to study the HIV-1 population in the female genital tract before virus is detectable in the bloodstream.

Heterosexual transmission of human immunodeficiency virus type 1 (HIV-1) is associated with a significant bottleneck in the viral quasispecies population, yet the timing of that bottleneck is poorly understood. We characterized HIV-1 diversity in the blood and female genital tract (FGT) within 2 weeks after detection of infection in three women enrolled in a unique prospective cohort in South Africa. We assembled full-length HIV-1 genomes from matched cervicovaginal lavage (CVL) samples and plasma. Deep sequencing allowed us to identify intrahost single-nucleotide variants (iSNVs) and to characterize within-sample HIV-1 diversity. Our results demonstrated very little HIV-1 diversity in the FGT and plasma by the time viremia was detectable. Within each subject, the consensus HIV-1 sequences were identical in plasma and CVL fluid. No iSNV was present at >6% frequency. One subject had 77 low-frequency iSNVs across both CVL fluid and plasma, another subject had 14 iSNVs in only CVL fluid from the earliest time point, and the third subject had no iSNVs in CVL fluid or plasma. Overall, the small amount of diversity that we detected was greater in the FGT than in plasma and declined over the first 2 weeks after viremia was detectable, compatible with a very early HIV-1 transmission bottleneck. To our knowledge, our study represents the earliest genomic analysis of HIV-1 in the FGT after transmission. Further, the use of metagenomic sequencing allowed us to characterize other organisms in the FGT, including commensal bacteria and sexually transmitted infections, highlighting the utility of the method to sequence both HIV-1 and its metagenomic environment.

IMPORTANCE Due to error-prone replication, HIV-1 generates a diverse population of viruses within a chronically infected individual. When HIV-1 is transmitted to a new individual, one or a few viruses establish the new infection, leading to a genetic bottleneck in the virus population. Understanding the timing and nature of this bottleneck may provide insight into HIV-1 vaccine design and other preventative strategies. We examined the HIV-1 population in three women enrolled in a unique prospective cohort in South Africa who were followed closely during the earliest stages of HIV-1 infection. We found very little HIV-1 diversity in the blood and female genital tract during the first 2 weeks after virus was detected in the bloodstream. These results are compatible with a very early HIV-1 population bottleneck, suggesting the need to study the HIV-1 population in the female genital tract before virus is detectable in the bloodstream.

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.

Tracking virus outbreaks in the twenty-first century

Emerging viruses have the potential to impose substantial mortality, morbidity and economic burdens on human populations. Tracking the spread of infectious diseases to assist in their control has traditionally relied on the analysis of case data gathered as the outbreak proceeds. Here, we describe how many of the key questions in infectious disease epidemiology, from the initial detection and characterization of outbreak viruses, to transmission chain tracking and outbreak mapping, can now be much more accurately addressed using recent advances in virus sequencing and phylogenetics. We highlight the utility of this approach with the hypothetical outbreak of an unknown pathogen, ‘Disease X’, suggested by the World Health Organization to be a potential cause of a future major epidemic. We also outline the requirements and challenges, including the need for flexible platforms that generate sequence data in real-time, and for these data to be shared as widely and openly as possible.

This Review Article describes how recent advances in viral genome sequencing and phylogenetics have enabled key issues associated with outbreak epidemiology to be more accurately addressed, and highlights the requirements and challenges for generating, sharing and using such data when tackling a viral outbreak.

miRTrace reveals the organismal origins of microRNA sequencing data

We present here miRTrace, the first algorithm to trace microRNA sequencing data back to their taxonomic origins. This is a challenge with profound implications for forensics, parasitology, food control, and research settings where cross-contamination can compromise results. miRTrace accurately (> 99%) assigns real and simulated data to 14 important animal and plant groups, sensitively detects parasitic infection in mammals, and discovers the primate origin of single cells. Applying our algorithm to over 700 public datasets, we find evidence that over 7% are cross-contaminated and present a novel solution to clean these computationally, even after sequencing has occurred. miRTrace is freely available at https://github.com/friedlanderlab/mirtrace .

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.

Using noninvasive metagenomics to characterize viral communities from wildlife

Microbial communities play an important role in organismal and ecosystem health. While high-throughput metabarcoding has revolutionized the study of bacterial communities, generating comparable viral communities has proven elusive, particularly in wildlife samples where the diversity of viruses and limited quantities of viral nucleic acid present distinctive challenges. Metagenomic sequencing is a promising solution for studying viral communities, but the lack of standardized methods currently precludes comparisons across host taxa or localities. Here, we developed an untargeted shotgun metagenomic sequencing protocol to generate comparable viral communities from noninvasively collected faecal and oropharyngeal swabs. Using samples from common vampire bats (Desmodus rotundus), a key species for virus transmission to humans and domestic animals, we tested how different storage media, nucleic acid extraction procedures and enrichment steps affect viral community detection. Based on finding viral contamination in foetal bovine serum, we recommend storing swabs in RNAlater or another nonbiological medium. We recommend extracting nucleic acid directly from swabs rather than from supernatant or pelleted material, which had undetectable levels of viral RNA. Results from a low-input RNA library preparation protocol suggest that ribosomal RNA depletion and light DNase treatment reduce host and bacterial nucleic acid, and improve virus detection. Finally, applying our approach to twelve pooled samples from seven localities in Peru, we showed that detected viral communities saturated at the attained sequencing depth, allowing unbiased comparisons of viral community composition. Future studies using the methods outlined here will elucidate the determinants of viral communities across host species, environments and time.

Characterization of Unknown Orthobunya-Like Viruses from India

Next-generation sequencing (NGS) of agents causing idiopathic human diseases has been crucial in the identification of novel viruses. This study describes the isolation and characterization of two novel orthobunyaviruses obtained from a jungle myna and a paddy bird from Karnataka State, India. Using an NGS approach, these isolates were classified as Cat Que and Balagodu viruses belonging to the Manzanilla clade of the Simbu serogroup. Closely related viruses in the Manzanilla clade have been isolated from mosquitos, humans, birds, and pigs across a wide geographic region. Since Orthobunyaviruses exhibit high reassortment frequency and can cause acute, self-limiting febrile illness, these data suggest that human and livestock infections of the Oya/Cat Que/Manzanilla virus may be more widespread and/or under-reported than anticipated. It therefore becomes imperative to identify novel and unknown viruses in order to understand their role in human and animal pathogenesis. The current study is a step forward in this regard and would act as a prototype method for isolation, identification and detection of several other emerging viruses.

Phylogeny of Yellow Fever Virus, Uganda, 2016

In April 2016, a yellow fever outbreak was detected in Uganda. Removal of contaminating ribosomal RNA in a clinical sample improved the sensitivity of next-generation sequencing. Molecular analyses determined the Uganda yellow fever outbreak was distinct from the concurrent yellow fever outbreak in Angola, improving our understanding of yellow fever epidemiology.

Field-deployable viral diagnostics using CRISPR-Cas13

Mitigating global infectious disease requires diagnostic tools that are sensitive, specific, and rapidly field deployable. In this study, we demonstrate that the Cas13-based SHERLOCK (specific high-sensitivity enzymatic reporter unlocking) platform can detect Zika virus (ZIKV) and dengue virus (DENV) in patient samples at concentrations as low as 1 copy per microliter. We developed HUDSON (heating unextracted diagnostic samples to obliterate nucleases), a protocol that pairs with SHERLOCK for viral detection directly from bodily fluids, enabling instrument-free DENV detection directly from patient samples in <2 hours. We further demonstrate that SHERLOCK can distinguish the four DENV serotypes, as well as region-specific strains of ZIKV from the 2015-2016 pandemic. Finally, we report the rapid (<1 week) design and testing of instrument-free assays to detect clinically relevant viral single-nucleotide polymorphisms.

Negative Selection by Spiral Inertial Microfluidics Improves Viral Recovery and Sequencing from Blood

In blood samples from patients with viral infection, it is often important to separate viral particles from human cells, for example, to minimize background in performing viral whole genome sequencing. Here, we present a microfluidic device that uses spiral inertial microfluidics with continuous circulation to separate host cells from viral particles and free nucleic acid. We demonstrate that this device effectively reduces white blood cells, red blood cells, and platelets from both whole blood and plasma samples with excellent recovery of viral nucleic acid. Furthermore, microfluidic separation leads to greater viral genome coverage and depth, highlighting an important application of this device in processing clinical samples for viral genome sequencing.

Rapid Detection of Powassan Virus in a Patient With Encephalitis by Metagenomic Sequencing

We describe a patient with severe and progressive encephalitis of unknown etiology. We performed rapid metagenomic sequencing from cerebrospinal fluid and identified Powassan virus, an emerging tick-borne flavivirus that has been increasingly detected in the United States.

Accurate RNA Sequencing From Formalin-Fixed Cancer Tissue To Represent High-Quality Transcriptome From Frozen Tissue

PURPOSE

Accurate transcriptional sequencing (RNA-seq) from formalin-fixation and paraffin-embedding (FFPE) tumor samples presents an important challenge for translational research and diagnostic development. In addition, there are now several different protocols to prepare a sequencing library from total RNA. We evaluated the accuracy of RNA-seq data generated from FFPE samples in terms of expression profiling.

METHODS

We designed a biospecimen study to directly compare gene expression results from different protocols to prepare libraries for RNA-seq from human breast cancer tissues, with randomization to fresh-frozen (FF) or FFPE conditions. The protocols were compared using multiple computational methods to assess alignment of reads to reference genome, and the uniformity and continuity of coverage; as well as the variance and correlation, of overall gene expression and patterns of measuring coding sequence, phenotypic patterns of gene expression, and measurements from representative multigene signatures.

RESULTS

The principal determinant of variance in gene expression was use of exon capture probes, followed by the conditions of preservation (FF versus FFPE), and phenotypic differences between breast cancers. One protocol, with RNase H-based rRNA depletion, exhibited least variability of gene expression measurements, strongest correlation between FF and FFPE samples, and was generally representative of the transcriptome from standard FF RNA-seq protocols.

CONCLUSION

Method of RNA-seq library preparation from FFPE samples had marked effect on the accuracy of gene expression measurement compared to matched FF samples. Nevertheless, some protocols produced highly concordant expression data from FFPE RNA-seq data, compared to RNA-seq results from matched frozen samples.