unbiased-deep-sequencing-of-rna-viruses-from-clinical-samples

Unbiased Deep Sequencing of RNA Viruses from Clinical Samples

Here we outline a next-generation RNA sequencing protocol that enables de novo assemblies and intra-host variant calls of viral genomes collected from ...

Broad Institute of MIT and Harvard

RNA interference (RNAi) is a gene-silencing phenomenon whereby double-stranded RNA (dsRNA) triggers the sequence-specific degradation of homologous mRNA. RNAi has been quickly and widely applied to discover gene functions and holds great potential to provide a new class of therapeutic agents. However, new chemistry and delivery approaches are greatly needed to silence disease-causing genes without toxic effects. We reasoned that conjugation of the cholesterol moiety to cationic lipids would enhance RNAi efficiencies and lower the toxic effects of lipid-mediated RNAi delivery. Here, we report the first design and synthesis of new cholesterol-conjugated cationic lipids for RNAi delivery using microwave-assisted quaternization (MAQ) of tertiary amines. This strategy can be employed to develop new classes of nonviral gene delivery agents under safe and fast reaction conditions.

Design and assembly of new nonviral RNAi delivery agents by microwave-assisted quaternization (MAQ) of tertiary amines.

RNA viruses are the causative agents for AIDS, influenza, SARS, and other serious health threats. Development of rapid and broadly applicable methods for complete viral genome sequencing is highly desirable to fully understand all aspects of these infectious agents as well as for surveillance of viral pandemic threats and emerging pathogens. However, traditional viral detection methods rely on prior sequence or antigen knowledge. In this study, we describe sequence-independent amplification for samples containing ultra-low amounts of viral RNA coupled with Illumina sequencing and de novo assembly optimized for viral genomes. With 5 million reads, we capture 96 to 100% of the viral protein coding region of HIV, respiratory syncytial and West Nile viral samples from as little as 100 copies of viral RNA. The methods presented here are scalable to large numbers of samples and capable of generating full or near full length viral genomes from clone and clinical samples with low amounts of viral RNA, without prior sequence information and in the presence of substantial host contamination.

Complete viral RNA genome sequencing of ultra-low copy samples by sequence-independent amplification.

Extensive genetic diversity in viral populations within infected hosts and the divergence of variants from existing reference genomes impede the analysis of deep viral sequencing data. A de novo population consensus assembly is valuable both as a single linear representation of the population and as a backbone on which intra-host variants can be accurately mapped. The availability of consensus assemblies and robustly mapped variants are crucial to the genetic study of viral disease progression, transmission dynamics, and viral evolution. Existing de novo assembly techniques fail to robustly assemble ultra-deep sequence data from genetically heterogeneous populations such as viruses into full-length genomes due to the presence of extensive genetic variability, contaminants, and variable sequence coverage.

De novo assembly of highly diverse viral populations.

The current methods for distinguishing the zygosities of twins include ultrasound scanning, which is nondefinitive, and amniocentesis, which is invasive. We explored the use of massively parallel sequencing of maternal plasma DNA for the noninvasive prenatal assessment of the zygosities of twin pregnancies.

Noninvasive prenatal determination of twin zygosity by maternal plasma DNA analysis.

This study aimed to provide an individualized assessment of fetal trisomy 21 and trisomy 18 status for twin pregnancies by maternal plasma DNA sequencing.

Noninvasive twin zygosity assessment and aneuploidy detection by maternal plasma DNA sequencing.

Human respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract disease in infants and young children and an important respiratory pathogen in the elderly and immunocompromised. While population-wide molecular epidemiology studies have shown multiple cocirculating RSV genotypes and revealed antigenic and genetic change over successive seasons, little is known about the extent of viral diversity over the course of an individual infection, the origins of novel variants, or the effect of immune pressure on viral diversity and potential immune-escape mutations. To investigate viral population diversity in the presence and absence of selective immune pressures, we studied whole-genome deep sequencing of RSV in upper airway samples from an infant with severe combined immune deficiency syndrome and persistent RSV infection. The infection continued over several months before and after bone marrow transplant (BMT) from his RSV-immune father. RSV diversity was characterized in 26 samples obtained over 78 days. Diversity increased after engraftment, as defined by T-cell presence, and populations reflected variation mostly within the G protein, the major surface antigen. Minority populations with known palivizumab resistance mutations emerged after its administration. The viral population appeared to diversify in response to selective pressures, showing a statistically significant growth in diversity in the presence of pressure from immunity. Defining escape mutations and their dynamics will be useful in the design and application of novel therapeutics and vaccines. These data can contribute to future studies of the relationship between within-host and population-wide RSV phylodynamics.

Within-host whole-genome deep sequencing and diversity analysis of human respiratory syncytial virus infection reveals dynamics of genomic diversity in the absence and presence of immune pressure.

In its largest outbreak, Ebola virus disease is spreading through Guinea, Liberia, Sierra Leone, and Nigeria. We sequenced 99 Ebola virus genomes from 78 patients in Sierra Leone to ~2000× coverage. We observed a rapid accumulation of interhost and intrahost genetic variation, allowing us to characterize patterns of viral transmission over the initial weeks of the epidemic. This West African variant likely diverged from central African lineages around 2004, crossed from Guinea to Sierra Leone in May 2014, and has exhibited sustained human-to-human transmission subsequently, with no evidence of additional zoonotic sources. Because many of the mutations alter protein sequences and other biologically meaningful targets, they should be monitored for impact on diagnostics, vaccines, and therapies critical to outbreak response.

Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak.

The 2013-2015 Ebola virus disease (EVD) epidemic is caused by the Makona variant of Ebola virus (EBOV). Early in the epidemic, genome sequencing provided insights into virus evolution and transmission and offered important information for outbreak response. Here, we analyze sequences from 232 patients sampled over 7 months in Sierra Leone, along with 86 previously released genomes from earlier in the epidemic. We confirm sustained human-to-human transmission within Sierra Leone and find no evidence for import or export of EBOV across national borders after its initial introduction. Using high-depth replicate sequencing, we observe both host-to-host transmission and recurrent emergence of intrahost genetic variants. We trace the increasing impact of purifying selection in suppressing the accumulation of nonsynonymous mutations over time. Finally, we note changes in the mucin-like domain of EBOV glycoprotein that merit further investigation. These findings clarify the movement of EBOV within the region and describe viral evolution during prolonged human-to-human transmission.

Ebola Virus Epidemiology, Transmission, and Evolution during Seven Months in Sierra Leone.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The introduction of West Nile virus (WNV) into North America in 1999 is a classic example of viral emergence in a new environment, with its subsequent dispersion across the continent having a major impact on local bird populations. Despite the importance of this epizootic, the pattern, dynamics, and determinants of WNV spread in its natural hosts remain uncertain. In particular, it is unclear whether the virus encountered major barriers to transmission, or spread in an unconstrained manner, and if specific viral lineages were favored over others indicative of intrinsic differences in fitness. To address these key questions in WNV evolution and ecology, we sequenced the complete genomes of approximately 300 avian isolates sampled across the United States between 2001 and 2012. Phylogenetic analysis revealed a relatively star-like tree structure, indicative of explosive viral spread in the United States, although with some replacement of viral genotypes through time. These data are striking in that viral sequences exhibit relatively limited clustering according to geographic region, particularly for those viruses sampled from birds, and no strong phylogenetic association with well-sampled avian species. The genome sequence data analyzed here also contain relatively little evidence for adaptive evolution, particularly of structural proteins, suggesting that most viral lineages are of similar fitness and that WNV is well adapted to the ecology of mosquito vectors and diverse avian hosts in the United States. In sum, the molecular evolution of WNV in North America depicts a largely unfettered expansion within a permissive host and geographic population with little evidence of major adaptive barriers.

Fluid Spatial Dynamics of West Nile Virus in the United States: Rapid Spread in a Permissive Host Environment.

Targeted next-generation sequencing is becoming increasingly common as a clinical diagnostic and prognostic test for patient- and tumor-specific genetic profiles as well as to optimally select targeted therapies. Here, we describe a custom-developed, next-generation sequencing test for detecting single-nucleotide variants (SNVs) and short insertions and deletions (indels) in 93 genes related to gastrointestinal cancer from routine formalin-fixed, paraffin-embedded clinical specimens. We implemented a validation strategy, based on the College of American Pathologists requirements, using reference DNA mixtures from cell lines with known genetic variants, which model a broad range of allele frequencies. Test sensitivity achieved >99% for both SNVs and indels, with allele frequencies >10%, with high specificity (97.4% for SNVs and 93.6% for indels). We further confirmed test accuracies using primary formalin-fixed, paraffin-embedded colorectal cancer specimens characterized by alternative and conventional clinical diagnostic technologies. Robust performance was observed on the formalin-fixed, paraffin-embedded specimens: sensitivity was 97.2% and specificity was 99.2%. We also observed high intrarun and inter-run reproducibility, as well as a low cross-contamination rate. Overall assessment using cell line samples and formalin-fixed, paraffin-embedded samples showed that our custom next-generation sequencing assay has consistent detection sensitivity down to 10% variant frequency.

James Qu

Broad Institute of MIT and Harvard

Unbiased Deep Sequencing of RNA Viruses from Clinical Samples