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 clinical and biological sources. The method is unbiased and universal; it uses random primers for cDNA synthesis and requires no prior knowledge of the viral sequence content. Before library construction, selective RNase H-based digestion is used to deplete unwanted RNA — including poly(rA) carrier and ribosomal RNA — from the viral RNA sample. Selective depletion improves both the data quality and the number of unique reads in viral RNA sequencing libraries. Moreover, a transposase-based ''tagmentation'' step is used in the protocol as it reduces overall library construction time. The protocol has enabled rapid deep sequencing of over 600 Lassa and Ebola virus samples-including collections from both blood and tissue isolates-and is broadly applicable to other microbial genomics studies.     

Individual and Bivalent Vaccines Based on Alphavirus Replicons Protect Guinea Pigs against Infection with Lassa and Ebola Viruses

Lassa and Ebola viruses cause acute, often fatal, hemorrhagic fever diseases, for which no effective vaccines are currently available. Although lethal human disease outbreaks have been confined so far to sub-Saharan Africa, they also pose significant epidemiological concern worldwide as demonstrated by several instances of accidental importation of the viruses into North America and Europe. In the present study, we developed experimental individual vaccines for Lassa virus and bivalent vaccines for Lassa and Ebola viruses that are based on an RNA replicon vector derived from an attenuated strain of Venezuelan equine encephalitis virus. The Lassa and Ebola virus genes were expressed from recombinant replicon RNAs that also encoded the replicase function and were capable of efficient intracellular self-amplification. For vaccinations, the recombinant replicons were incorporated into virus-like replicon particles. Guinea pigs vaccinated with particles expressing Lassa virus nucleoprotein or glycoprotein genes were protected from lethal challenge with Lassa virus. Vaccination with particles expressing Ebola virus glycoprotein gene also protected the animals from lethal challenge with Ebola virus. In order to evaluate a single vaccine protecting against both Lassa and Ebola viruses, we developed dual-expression particles that expressed glycoprotein genes of both Ebola and Lassa viruses. Vaccination of guinea pigs with either dual-expression particles or with a mixture of particles expressing Ebola and Lassa virus glycoprotein genes protected the animals against challenges with Ebola and Lassa viruses. The results showed that immune responses can be induced against multiple vaccine antigens coexpressed from an alphavirus replicon and suggested the possibility of engineering multivalent vaccines based upon alphavirus vectors for arenaviruses, filoviruses, and possibly other emerging pathogens.     

Lassa viral dynamics in non-human primates treated with favipiravir or ribavirin

Lassa fever is an haemorrhagic fever caused by Lassa virus (LASV). There is no vaccine approved against LASV and the only recommended antiviral treatment relies on ribavirin, despite limited evidence of efficacy. Recently, the nucleotide analogue favipiravir showed a high antiviral efficacy, with 100% survival obtained in an otherwise fully lethal non-human primate (NHP) model of Lassa fever. However the mechanism of action of the drug is not known and the absence of pharmacokinetic data limits the translation of these results to the human setting. Here we aimed to better understand the antiviral effect of favipiravir by developping the first mathematical model recapitulating Lassa viral dynamics and treatment. We analyzed the viral dynamics in 24 NHPs left untreated or treated with ribavirin or favipiravir, and we put the results in perspective with those obtained with the same drugs in the context of Ebola infection. Our model estimates favipiravir EC₅₀ in vivo to 2.89 μg.mL^-1, which is much lower than what was found against Ebola virus. The main mechanism of action of favipiravir was to decrease virus infectivity, with an efficacy of 91% at the highest dose. Based on our knowledge acquired on the drug pharmacokinetics in humans, our model predicts that favipiravir doses larger than 1200 mg twice a day should have the capability to strongly reduce the production infectious virus and provide a milestone towards a future use in humans.     

Geographic Distribution and Genetic Characterization of Lassa Virus in Sub-Saharan Mali

Background

Lassa fever is an acute viral illness characterized by multi-organ failure and hemorrhagic manifestations. Lassa fever is most frequently diagnosed in Nigeria, Sierra Leone, Liberia, and Guinea, although sporadic cases have been recorded in other West African countries, including Mali. The etiological agent of Lassa fever is Lassa virus (LASV), an Arenavirus which is maintained in nature and frequently transmitted to humans by Mastomys natalensis. The purpose of this study was to better define the geographic distribution of LASV-infected rodents in sub-Saharan Mali.

Methodologies/Principal Findings

Small mammals were live-trapped at various locations across Mali for the purpose of identifying potential zoonotic pathogens. Serological and molecular assays were employed and determined LASV infected rodents were exclusively found in the southern Mali near the border of Côte d''Ivoire. Overall, 19.4% of Mastomys natalensis sampled in this region had evidence of LASV infection, with prevalence rates for individual villages ranging from 0 to 52%. Full-length genomic sequences were determined using high throughput sequencing methodologies for LASV isolates generated from tissue samples of rodents collected in four villages and confirmed the phylogenetic clustering of Malian LASV with strain AV.

Conclusions/Significance

The risk of human infections with LASV is greatest in villages in southern Mali. Lassa fever should be considered in the differential diagnosis for febrile individuals and appropriate diagnostic techniques need to be established to determine the incidence of infection and disease in these regions.     

Analysis of the genetic diversity of influenza A viruses using next-generation DNA sequencing

Background

Influenza viruses exist as a large group of closely related viral genomes, also called quasispecies. The composition of this influenza viral quasispecies can be determined by an accurate and sensitive sequencing technique and data analysis pipeline. We compared the suitability of two benchtop next-generation sequencers for whole genome influenza A quasispecies analysis: the Illumina MiSeq sequencing-by-synthesis and the Ion Torrent PGM semiconductor sequencing technique.

Results

We first compared the accuracy and sensitivity of both sequencers using plasmid DNA and different ratios of wild type and mutant plasmid. Illumina MiSeq sequencing reads were one and a half times more accurate than those of the Ion Torrent PGM. The majority of sequencing errors were substitutions on the Illumina MiSeq and insertions and deletions, mostly in homopolymer regions, on the Ion Torrent PGM. To evaluate the suitability of the two techniques for determining the genome diversity of influenza A virus, we generated plasmid-derived PR8 virus and grew this virus in vitro. We also optimized an RT-PCR protocol to obtain uniform coverage of all eight genomic RNA segments. The sequencing reads obtained with both sequencers could successfully be assembled de novo into the segmented influenza virus genome. After mapping of the reads to the reference genome, we found that the detection limit for reliable recognition of variants in the viral genome required a frequency of 0.5% or higher. This threshold exceeds the background error rate resulting from the RT-PCR reaction and the sequencing method. Most of the variants in the PR8 virus genome were present in hemagglutinin, and these mutations were detected by both sequencers.

Conclusions

Our approach underlines the power and limitations of two commonly used next-generation sequencers for the analysis of influenza virus gene diversity. We conclude that the Illumina MiSeq platform is better suited for detecting variant sequences whereas the Ion Torrent PGM platform has a shorter turnaround time. The data analysis pipeline that we propose here will also help to standardize variant calling in small RNA genomes based on next-generation sequencing data.     

The organisation of Ebola virus reveals a capacity for extensive, modular polyploidy

Background

Filoviruses, including Ebola virus, are unusual in being filamentous animal viruses. Structural data on the arrangement, stoichiometry and organisation of the component molecules of filoviruses has until now been lacking, partially due to the need to work under level 4 biological containment. The present study provides unique insights into the structure of this deadly pathogen.

Methodology and Principal Findings

We have investigated the structure of Ebola virus using a combination of cryo-electron microscopy, cryo-electron tomography, sub-tomogram averaging, and single particle image processing. Here we report the three-dimensional structure and architecture of Ebola virus and establish that multiple copies of the RNA genome can be packaged to produce polyploid virus particles, through an extreme degree of length polymorphism. We show that the helical Ebola virus inner nucleocapsid containing RNA and nucleoprotein is stabilized by an outer layer of VP24-VP35 bridges. Elucidation of the structure of the membrane-associated glycoprotein in its native state indicates that the putative receptor-binding site is occluded within the molecule, while a major neutralizing epitope is exposed on its surface proximal to the viral envelope. The matrix protein VP40 forms a regular lattice within the envelope, although its contacts with the nucleocapsid are irregular.

Conclusions

The results of this study demonstrate a modular organization in Ebola virus that accommodates a well-ordered, symmetrical nucleocapsid within a flexible, tubular membrane envelope.     

Optimized microRNA purification from TRIzol-treated plasma

Background

MicroRNAs (miRNAs) represent new and potentially informative diagnostic targets for disease detection and prognosis. However, little work exists documenting the effect of TRIzol, a common viral inactivation and nucleic acid extraction reagent, on miRNA purification. Here, we developed an optimized protocol for miRNA extraction from plasma samples by evaluating five different RNA extraction kits, TRIzol phase separation, purification additives, and initial plasma sample volume. This method was then used for downstream profiling of plasma miRNAs found in archived samples from one nonhuman primate (NHP) experimentally challenged with Ebola virus by the aerosol route.

Results

Comparison of real-time RT-PCR results for spiked-in and endogenous miRNA sequences determined extraction efficiencies from five different RNA purification kits. These experiments showed that 50 μL plasma processed using the QIAGEN miRNeasy Mini Kit with 5 μg of glycogen as a co-precipitant yielded the highest recovery of endogenous miRNAs. Using this optimized protocol, miRNAs from archived plasma samples of one rhesus macaque challenged with aerosolized Ebola virus was profiled using a targeted real-time PCR array. A total of 519 of the 752 unique miRNAs assayed were present in the plasma samples at day 0 and day 7 (time of death) post-exposure. Statistical analyses revealed 25 sequences significantly up- or down-regulated between day 0 and day 7 post infection, validating the utility of the extraction method for plasma miRNA profiling.

Conclusions

This study contributes to the knowledgebase of circulating miRNA extraction methods and expands on the potential applications of cell-free miRNA profiling for diagnostics and pathogenesis studies. Specifically, we optimized an extraction protocol for miRNAs from TRIzol-inactivated plasma samples that can be used for highly pathogenic viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1299-5) contains supplementary material, which is available to authorized users.     

Flexible docking-based molecular dynamics simulation of natural product compounds and Ebola virus Nucleocapsid (EBOV NP): a computational approach to discover new drug for combating Ebola

Background

Ebola still remains as one of the most problematic infectious diseases in Africa with a high rate of mortality. Although this disease has been known for an almost half-century, there are no vaccines and drugs available in the market to treat Ebola. Zaire ebolavirus (EBOV), a single-stranded RNA virus which belongs to Filoviridae family and Mononegavirales order, is one of the virus causing Ebola. As one of seven proteins that EBOV encodes, Ebola virus nucleoprotein (EBOV NP) plays an imperative role in EBOV proliferation cycle. Therefore, the development of a new Ebola treatment can be targeted towards EBOV NP.

Results

In this work, we screened about 190,084 natural product compounds from ZINC15 database through in silico virtual screening and flexible docking simulation. Furthermore, the bioavailability and toxicity prediction have been conducted as well. Two best ligands according to the simulation and prediction tests were progressed into the molecular dynamics simulation.

Conclusion

In the end, we found that our proposed ligands, namely α-lipomycin (ZINC56874155) and 3-(((S)-1-amino-1,2,3,4-tetrahydroisoquinolin-5-yl)methyl)-5-((5-((5R,7S)-5,7-dihydroxy-3-oxodecyl)-2-hydroxyphenoxy) methyl)pyrrolo[3,4-b]pyrrol-5-ium (ZINC85628951), showed the promising results to be developed as a lead compounds for treating Ebola. Therefore, an experimental study is required to validate their inhibition activities against EBOV NP.

     

Ebola Virus Genome Plasticity as a Marker of Its Passaging History: A Comparison of In Vitro Passaging to Non-Human Primate Infection

To identify polymorphic sites that could be used as biomarkers of Ebola virus passage history, we repeatedly amplified Ebola virus (Kikwit variant) in vitro and in vivo and performed deep sequencing analysis of the complete genomes of the viral subpopulations. We then determined the sites undergoing selection during passage in Vero E6 cells. Four locations within the Ebola virus Kikwit genome were identified that together segregate cell culture-passaged virus and virus obtained from infected non-human primates. Three of the identified sites are located within the glycoprotein gene (GP) sequence: the poly-U (RNA editing) site at position 6925, as well as positions 6677, and 6179. One site was found in the VP24 gene at position 10833. In all cases, in vitro and in vivo, both populations (majority and minority variants) were maintained in the viral swarm, with rapid selections occurring after a few passages or infections. This analysis approach will be useful to differentiate whether filovirus stocks with unknown history have been passaged in cell culture and may support filovirus stock standardization for medical countermeasure development.     

A Bead-based Normalization for Uniform Sequencing depth (BeNUS) protocol for multi-samples sequencing exemplified by HLA-B

Background

Human leukocyte antigen (HLA) is a group of genes that are extremely polymorphic among individuals and populations and have been associated with more than 100 different diseases and adverse drug effects. HLA typing is accordingly an important tool in clinical application, medical research, and population genetics. We have previously developed a phase-defined HLA gene sequencing method using MiSeq sequencing.

Results

Here we report a simple, high-throughput, and cost-effective sequencing method that includes normalized library preparation and adjustment of DNA molar concentration. We applied long-range PCR to amplify HLA-B for 96 samples followed by transposase-based library construction and multiplex sequencing with the MiSeq sequencer. After sequencing, we observed low variation in read percentages (0.2% to 1.55%) among the 96 demultiplexed samples. On this basis, all the samples were amenable to haplotype phasing using our phase-defined sequencing method. In our study, a sequencing depth of 800x was necessary and sufficient to achieve full phasing of HLA-B alleles with reliable assignment of the allelic sequence to the 8 digit level.

Conclusions

Our HLA sequencing method optimized for 96 multiplexing samples is highly time effective and cost effective and is especially suitable for automated multi-sample library preparation and sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-645) contains supplementary material, which is available to authorized users.     

A Systematic Screen of FDA-Approved Drugs for Inhibitors of Biological Threat Agents

Background

The rapid development of effective medical countermeasures against potential biological threat agents is vital. Repurposing existing drugs that may have unanticipated activities as potential countermeasures is one way to meet this important goal, since currently approved drugs already have well-established safety and pharmacokinetic profiles in patients, as well as manufacturing and distribution networks. Therefore, approved drugs could rapidly be made available for a new indication in an emergency.

Methodology/Principal Findings

A large systematic effort to determine whether existing drugs can be used against high containment bacterial and viral pathogens is described. We assembled and screened 1012 FDA-approved drugs for off-label broad-spectrum efficacy against Bacillus anthracis; Francisella tularensis; Coxiella burnetii; and Ebola, Marburg, and Lassa fever viruses using in vitro cell culture assays. We found a variety of hits against two or more of these biological threat pathogens, which were validated in secondary assays. As expected, antibiotic compounds were highly active against bacterial agents, but we did not identify any non-antibiotic compounds with broad-spectrum antibacterial activity. Lomefloxacin and erythromycin were found to be the most potent compounds in vivo protecting mice against Bacillus anthracis challenge. While multiple virus-specific inhibitors were identified, the most noteworthy antiviral compound identified was chloroquine, which disrupted entry and replication of two or more viruses in vitro and protected mice against Ebola virus challenge in vivo.

Conclusions/Significance

The feasibility of repurposing existing drugs to face novel threats is demonstrated and this represents the first effort to apply this approach to high containment bacteria and viruses.     

A new platform for ultra-high density Staphylococcus aureus transposon libraries

Background

Staphylococcus aureus readily develops resistance to antibiotics and achieving effective therapies to overcome resistance requires in-depth understanding of S. aureus biology. High throughput, parallel-sequencing methods for analyzing transposon mutant libraries have the potential to revolutionize studies of S. aureus, but the genetic tools to take advantage of the power of next generation sequencing have not been fully developed.

Results

Here we report a phage-based transposition system to make ultra-high density transposon libraries for genome-wide analysis of mutant fitness in any Φ11-transducible S. aureus strain. The high efficiency of the delivery system has made it possible to multiplex transposon cassettes containing different regulatory elements in order to make libraries in which genes are over- or under-expressed as well as deleted. By incorporating transposon-specific barcodes into the cassettes, we can evaluate how null mutations and changes in gene expression levels affect fitness in a single sequencing data set. Demonstrating the power of the system, we have prepared a library containing more than 690,000 unique insertions. Because one unique feature of the phage-based approach is that temperature-sensitive mutants are retained, we have carried out a genome-wide study of S. aureus genes involved in withstanding temperature stress. We find that many genes previously identified as essential are temperature sensitive and also identify a number of genes that, when disrupted, confer a growth advantage at elevated temperatures.

Conclusions

The platform described here reliably provides mutant collections of unparalleled genotypic diversity and will enable a wide range of functional genomic studies in S. aureus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1361-3) contains supplementary material, which is available to authorized users.     

Evaluating bias-reducing protocols for RNA sequencing library preparation

Background

Next-generation sequencing does not yield fully unbiased estimates for read abundance, which may impact on the conclusions that can be drawn from sequencing data. The ligation step in RNA sequencing library generation is a known source of bias, motivating developments in enzyme technology and library construction protocols. We present the first comparison of the standard duplex adaptor protocol supplied by Life Technologies for use on the Ion Torrent PGM with an alternate single adaptor approach involving CircLigase (CircLig protocol).A correlation between over-representation in sequenced libraries and degree of secondary structure has been reported previously, therefore we also investigated whether bias could be reduced by ligation with an enzyme that functions at a temperature not permissive for such structure.

Results

A pool of small RNA fragments of known composition was converted into a sequencing library using one of three protocols and sequenced on an Ion Torrent PGM. The CircLig protocol resulted in less over-representation of specific sequences than the standard protocol. Over-represented sequences are more likely to be predicted to have secondary structure and to co-fold with adaptor sequences. However, use of the thermostable ligase Methanobacterium thermoautotrophicum RNA ligase K97A (Mth K97A) was not sufficient to reduce bias.

Conclusions

The single adaptor CircLigase-based approach significantly reduces, but does not eliminate, bias in Ion Torrent data. Ligases that function at temperatures to remove the possible influence of secondary structure on library generation may be of value, although Mth K97A is not effective in this case.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-569) contains supplementary material, which is available to authorized users.