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

Parasites have developed myriad strategies to reach and infect their specific hosts. One of the most common mechanisms for non-vector transmitted parasites to reach the internal host environment is by ingestion during feeding. In this study, we investigated the mechanisms of oyster host colonization by the alveolate Perkinsus marinus and focused on how oysters process infective waterborne P. marinus cells during feeding in order to determine the portal(s) of entry of this parasite to its host. We also compared the infectivity of freely-suspended cells of P. marinus with that of cells incorporated into marine aggregates to link changes in particle processing by the feeding organs with infection success and route. Finally, we evaluated the effect of oyster secretions (mucus) covering the feeding organs on P. marinus physiology because these host factors are involved in the processing of waterborne particles. The ensemble of results shows a unique mechanism for infection by which the parasite is mostly acquired during the feeding process, but not via ingestion. Rather, infection commonly occurs during the rejection of material as pseudofeces before reaching the mouth. The pseudofeces discharge area, a specialized area of the mantle where unwanted particles are accumulated for rejection as pseudofeces, showed significantly higher parasite loads than other host tissues including other parts of the mantle. Aggregated P. marinus cells caused significantly higher disease prevalence and infection intensities when compared to freely-suspended parasite cells. Mucus covering the mantle caused a quick and significant increase in parasite replication rates suggesting rapid impact on P. marinus physiology. A new model for P. marinus acquisition in oysters is proposed.

Early host-pathogen interactions in marine bivalves: evidence that the alveolate parasite Perkinsus marinus infects through the oyster mantle during rejection of pseudofeces.

Although several hundred regions of the human genome harbor signals of positive natural selection, few of the relevant adaptive traits and variants have been elucidated. Using full-genome sequence variation from the 1000 Genomes (1000G) Project and the composite of multiple signals (CMS) test, we investigated 412 candidate signals and leveraged functional annotation, protein structure modeling, epigenetics, and association studies to identify and extensively annotate candidate causal variants. The resulting catalog provides a tractable list for experimental follow-up; it includes 35 high-scoring nonsynonymous variants, 59 variants associated with expression levels of a nearby coding gene or lincRNA, and numerous variants associated with susceptibility to infectious disease and other phenotypes. We experimentally characterized one candidate nonsynonymous variant in Toll-like receptor 5 (TLR5) and show that it leads to altered NF-κB signaling in response to bacterial flagellin. PAPERFLICK:

Identifying recent adaptations in large-scale genomic data.

Perkinsus marinus is an important protistan parasite of the eastern oyster Crassostrea virginica. Recent findings showed that oyster pallial organs (mantle, gills) are a major portal of entry for the parasite. Therefore, mucus covering these organs represents the first host effectors encountered by P. marinus. This study consisted of several experiments designed to investigate the effect of oyster pallial mucus on the growth, protease production and infectivity of P. marinus. In each experiment, P. marinus performance in cultures supplemented with pallial mucus (mantle, gill, or both) was compared to that of parasite cells grown in unsupplemented media or in cultures supplemented with oyster plasma or digestive extracts. P. marinus grown in media supplemented with C. virginica mantle mucus showed a significantly higher growth rate than cultures enriched with the other supplemental extracts, while cultures grown in gill mucus promoted higher protease production. Conversely, P. marinus grown in cultures supplemented with pallial mucus of the non-compatible host Crassostrea gigas (Pacific oyster) were dramatically inhibited. Challenge experiments showed a significant increase in P. marinus virulence in cultures supplemented with C. virginica pallial mucus as compared to unsupplemented cultures or to those supplemented with digestive extract or plasma. These results suggest that C. virginica mucus plays a significant role in the pathogenesis of P. marinus by enhancing the proliferation and the infectivity of this devastating parasite. The contrasting results obtained with both oyster species indicate that P. marinus host specificity may begin in the mucus.

Early host-pathogen interactions in a marine bivalve: Crassostrea virginica pallial mucus modulates Perkinsus marinus growth and virulence.

Research capacity. Enabling the genomic revolution in Africa.

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.

Limited clinical and laboratory data are available on patients with Ebola virus disease (EVD). The Kenema Government Hospital in Sierra Leone, which had an existing infrastructure for research regarding viral hemorrhagic fever, has received and cared for patients with EVD since the beginning of the outbreak in Sierra Leone in May 2014.

Clinical illness and outcomes in patients with Ebola in Sierra Leone.

We have developed a robust RNA sequencing method for generating complete de novo assemblies with intra-host variant calls of Lassa and Ebola virus genomes in clinical and biological samples. Our method uses targeted RNase H-based digestion to remove contaminating poly(rA) carrier and ribosomal RNA. This depletion step improves both the quality of data and quantity of informative reads in unbiased total RNA sequencing libraries. We have also developed a hybrid-selection protocol to further enrich the viral content of sequencing libraries. These protocols have enabled rapid deep sequencing of both Lassa and Ebola virus and are broadly applicable to other viral genomics studies.

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

Next-generation sequencing (NGS) has the potential to transform the discovery of viruses causing unexplained acute febrile illness (UAFI) because it does not depend on culturing the pathogen or a priori knowledge of the pathogen's nucleic acid sequence. More generally, it has the potential to elucidate the complete human virome, including viruses that cause no overt symptoms of disease, but may have unrecognized immunological or developmental consequences. We have used NGS to identify RNA viruses in the blood of 195 patients with UAFI and compared them with those found in 328 apparently healthy (i.e., no overt signs of illness) control individuals, all from communities in southeastern Nigeria. Among UAFI patients, we identified the presence of nucleic acids from several well-characterized pathogenic viruses, such as HIV-1, hepatitis, and Lassa virus. In our cohort of healthy individuals, however, we detected the nucleic acids of two novel rhabdoviruses. These viruses, which we call Ekpoma virus-1 (EKV-1) and Ekpoma virus-2 (EKV-2), are highly divergent, with little identity to each other or other known viruses. The most closely related rhabdoviruses are members of the genus Tibrovirus and Bas-Congo virus (BASV), which was recently identified in an individual with symptoms resembling hemorrhagic fever. Furthermore, by conducting a serosurvey of our study cohort, we find evidence for remarkably high exposure rates to the identified rhabdoviruses. The recent discoveries of novel rhabdoviruses by multiple research groups suggest that human infection with rhabdoviruses might be common. While the prevalence and clinical significance of these viruses are currently unknown, these viruses could have previously unrecognized impacts on human health; further research to understand the immunological and developmental impact of these viruses should be explored. More generally, the identification of similar novel viruses in individuals with and without overt symptoms of disease highlights the need for a broader understanding of the human virome as efforts for viral detection and discovery advance.

Discovery of novel rhabdoviruses in the blood of healthy individuals from West Africa.

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.

Although studies have identified hundreds of loci associated with human traits and diseases, pinpointing causal alleles remains difficult, particularly for non-coding variants. To address this challenge, we adapted the massively parallel reporter assay (MPRA) to identify variants that directly modulate gene expression. We applied it to 32,373 variants from 3,642 cis-expression quantitative trait loci and control regions. Detection by MPRA was strongly correlated with measures of regulatory function. We demonstrate MPRA's capabilities for pinpointing causal alleles, using it to identify 842 variants showing differential expression between alleles, including 53 well-annotated variants associated with diseases and traits. We investigated one in detail, a risk allele for ankylosing spondylitis, and provide direct evidence of a non-coding variant that alters expression of the prostaglandin EP4 receptor. These results create a resource of concrete leads and illustrate the promise of this approach for comprehensively interrogating how non-coding polymorphism shapes human biology.

Sarah Winnicki

Broad Institute of MIT and Harvard

Unbiased Deep Sequencing of RNA Viruses from Clinical Samples