A Modified RNA-Seq Approach for Whole Genome Sequencing of RNA Viruses from Faecal and Blood Samples

To date, very large scale sequencing of many clinically important RNA viruses has been complicated by their high population molecular variation, which creates challenges for polymerase chain reaction and sequencing primer design. Many RNA viruses are also difficult or currently not possible to culture, severely limiting the amount and purity of available starting material. Here, we describe a simple, novel, high-throughput approach to Norovirus and Hepatitis C virus whole genome sequence determination based on RNA shotgun sequencing (also known as RNA-Seq). We demonstrate the effectiveness of this method by sequencing three Norovirus samples from faeces and two Hepatitis C virus samples from blood, on an Illumina MiSeq benchtop sequencer. More than 97% of reference genomes were recovered. Compared with Sanger sequencing, our method had no nucleotide differences in 14,019 nucleotides (nt) for Noroviruses (from a total of 2 Norovirus genomes obtained with Sanger sequencing), and 8 variants in 9,542 nt for Hepatitis C virus (1 variant per 1,193 nt). The three Norovirus samples had 2, 3, and 2 distinct positions called as heterozygous, while the two Hepatitis C virus samples had 117 and 131 positions called as heterozygous. To confirm that our sample and library preparation could be scaled to true high-throughput, we prepared and sequenced an additional 77 Norovirus samples in a single batch on an Illumina HiSeq 2000 sequencer, recovering >90% of the reference genome in all but one sample. No discrepancies were observed across 118,757 nt compared between Sanger and our custom RNA-Seq method in 16 samples. By generating viral genomic sequences that are not biased by primer-specific amplification or enrichment, this method offers the prospect of large-scale, affordable studies of RNA viruses which could be adapted to routine diagnostic laboratory workflows in the near future, with the potential to directly characterize within-host viral diversity.     

Dimer initiation signal of human immunodeficiency virus type 1: its role in partner selection during RNA copackaging and its effects on recombination

Frequent human immunodeficiency virus type 1 (HIV-1) recombination occurs during DNA synthesis when portions of the two copackaged RNAs are used as templates to generate a hybrid DNA copy. Therefore, the frequency of copackaging of genomic RNAs from two different viruses (heterozygous virion formation) affects the generation of genotypically different recombinants. We hypothesized that the selection of copackaged RNA partners is largely determined by Watson-Crick pairing at the dimer initiation signal (DIS), a 6-nucleotide palindromic sequence at the terminal loop of stem-loop 1 (SL1). To test our hypothesis, we examined whether heterozygous virion formation could be encouraged by manipulation of the DIS. Three pairs of viruses were generated with compensatory DIS mutations, designed so that perfect DIS base pairing could only occur between RNAs derived from different viruses, not between RNAs from the same virus. We observed that vector pairs with compensatory DIS mutations had an almost twofold increase in recombination rates compared with wild-type viruses. These data suggest that heterozygous virion formation was enhanced in viruses with compensatory DIS mutations (from 50% to more than 90% in some viral pairings). The role of the SL1 stem in heterozygous virion formation was also tested; our results indicated that the intermolecular base pairing of the stem sequences does not affect RNA partner selection. In summary, our results demonstrate that the Watson-Crick pairing of the DIS is a major determinant in the selection of the copackaged RNA partner, and altering the base pairing of the DIS can change the proportion of heterozygous viruses in a viral population. These results also strongly support the hypothesis that HIV-1 RNA dimers are formed prior to encapsidation.     

Amplification methods bias metagenomic libraries of uncultured single-stranded and double-stranded DNA viruses

Investigation of viruses in the environment often requires the amplification of viral DNA before sequencing of viral metagenomes. In this study, two of the most widely used amplification methods, the linker amplified shotgun library (LASL) and multiple displacement amplification (MDA) methods, were applied to a sample from the seawater surface. Viral DNA was extracted from viruses concentrated by tangential flow filtration and amplified by these two methods. 454 pyrosequencing was used to read the metagenomic sequences from different libraries. The resulting taxonomic classifications of the viruses, their functional assignments, and assembly patterns differed substantially depending on the amplification method. Only double-stranded DNA viruses were retrieved from the LASL, whereas most sequences in the MDA library were from single-stranded DNA viruses, and double-stranded DNA viral sequences were minorities. Thus, the two amplification methods reveal different aspects of viral diversity.     

RNA viral community in human feces: prevalence of plant pathogenic viruses

The human gut is known to be a reservoir of a wide variety of microbes, including viruses. Many RNA viruses are known to be associated with gastroenteritis; however, the enteric RNA viral community present in healthy humans has not been described. Here, we present a comparative metagenomic analysis of the RNA viruses found in three fecal samples from two healthy human individuals. For this study, uncultured viruses were concentrated by tangential flow filtration, and viral RNA was extracted and cloned into shotgun viral cDNA libraries for sequencing analysis. The vast majority of the 36,769 viral sequences obtained were similar to plant pathogenic RNA viruses. The most abundant fecal virus in this study was pepper mild mottle virus (PMMV), which was found in high concentrations—up to 10⁹ virions per gram of dry weight fecal matter. PMMV was also detected in 12 (66.7%) of 18 fecal samples collected from healthy individuals on two continents, indicating that this plant virus is prevalent in the human population. A number of pepper-based foods tested positive for PMMV, suggesting dietary origins for this virus. Intriguingly, the fecal PMMV was infectious to host plants, suggesting that humans might act as a vehicle for the dissemination of certain plant viruses.     

Genome packaging within icosahedral capsids and large-scale segmentation in viral genomic sequences

The assembly and maturation of viruses with icosahedral capsids must be coordinated with icosahedral symmetry. The icosahedral symmetry imposes also the restrictions on the cooperative specific interactions between genomic RNA/DNA and coat proteins that should be reflected in quasi-regular segmentation of viral genomic sequences. Combining discrete direct and double Fourier transforms, we studied the quasi-regular large-scale segmentation in genomic sequences of different ssRNA, ssDNA, and dsDNA viruses. The particular representatives included satellite tobacco mosaic virus (STMV) and the strains of satellite tobacco necrosis virus (STNV), STNV-C, STNV-1, STNV-2, Escherichia phages MS2, ?X174, α3, and HK97, and Simian virus 40. In all their genomes, we found the significant quasi-regular segmentation of genomic sequences related to the virion assembly and the genome packaging within icosahedral capsid. We also found good correspondence between our results and available cryo-electron microscopy data on capsid structures and genome packaging in these viruses. Fourier analysis of genomic sequences provides the additional insight into mechanisms of hierarchical genome packaging and may be used for verification of the concepts of 3-fold or 5-fold intermediates in virion assembly. The results of sequence analysis should be taken into account at the choice of models and data interpretation. They also may be helpful for the development of antiviral drugs.     

The efficacy of Cot-based gene enrichment in wheat (Triticum aestivum L.).

We report the results of a study on the effectiveness of Cot filtration (CF) in the characterization of the gene space of bread wheat (Triticum aestivum L.), a large genome species (1C = 16,700 Mb) of tremendous agronomic importance. Using published Cot data as a guide, 2 genomic libraries for hexaploid wheat were constructed from the single-stranded DNA collected at Cot values > 1188 and 1639 M x s. Compared with sequences from a whole genome shotgun library from Aegilops tauschii (the D genome donor of bread wheat), the CF libraries exhibited 13.7-fold enrichment in genes, 5.8-fold enrichment in unknown low-copy sequences, and a 3-fold reduction in repetitive DNA. CF is twice as efficient as methylation filtration at enriching wheat genes. This research suggests that, with improvements, CF will be a highly useful tool in sequencing the gene space of wheat.     

Simultaneous identification of DNA and RNA viruses present in pig faeces using process-controlled deep sequencing

Background

Animal faeces comprise a community of many different microorganisms including bacteria and viruses. Only scarce information is available about the diversity of viruses present in the faeces of pigs. Here we describe a protocol, which was optimized for the purification of the total fraction of viral particles from pig faeces. The genomes of the purified DNA and RNA viruses were simultaneously amplified by PCR and subjected to deep sequencing followed by bioinformatic analyses. The efficiency of the method was monitored using a process control consisting of three bacteriophages (T4, M13 and MS2) with different morphology and genome types. Defined amounts of the bacteriophages were added to the sample and their abundance was assessed by quantitative PCR during the preparation procedure.

Results

The procedure was applied to a pooled faecal sample of five pigs. From this sample, 69,613 sequence reads were generated. All of the added bacteriophages were identified by sequence analysis of the reads. In total, 7.7% of the reads showed significant sequence identities with published viral sequences. They mainly originated from bacteriophages (73.9%) and mammalian viruses (23.9%); 0.8% of the sequences showed identities to plant viruses. The most abundant detected porcine viruses were kobuvirus, rotavirus C, astrovirus, enterovirus B, sapovirus and picobirnavirus. In addition, sequences with identities to the chimpanzee stool-associated circular ssDNA virus were identified. Whole genome analysis indicates that this virus, tentatively designated as pig stool-associated circular ssDNA virus (PigSCV), represents a novel pig virus.

Conclusion

The established protocol enables the simultaneous detection of DNA and RNA viruses in pig faeces including the identification of so far unknown viruses. It may be applied in studies investigating aetiology, epidemiology and ecology of diseases. The implemented process control serves as quality control, ensures comparability of the method and may be used for further method optimization.     

Structural homology among four nodaviruses as deduced by sequencing and X-ray crystallography

The genomic RNA2s of nodaviruses encode a single gene, that of protein alpha, the precursor of virion proteins beta and gamma. We compared the sequences of the RNA2s of the nodaviruses, black beetle virus (BBV), flock house virus, boolarra virus and nodamura virus, with the objective of identifying homologies in the primary and secondary structure of these RNAs and in the structure of their encoded protein. The sequences of the four RNAs were found to be similar, so that homologous regions relating to translation and RNA replication were readily identified. However, the overall, secondary structures in solution, deduced from calculations of optimal Watson-Crick base-pairing configurations, were very different for the four RNAs. We conclude that a particular, overall, secondary structure in solution within host cells is not required for virus viability. The partially refined X-ray structure of BBV (R = 26.4% for the current model) was used as a framework for comparing the structure of the encoded proteins of the four viruses. Mapping of the four protein sequences onto the BBV capsid showed many amino acid differences on the outer surface, indicating that the exteriors of the four virions are substantially different. Mapping in the beta-barrel region showed an intermediate level of differences, indicating that some freedom in choice of amino acid residues is possible there although the basic framework of the capsids is evidently conserved. Mapping onto the interior surface of the BBV capsid showed a high degree of conservation of amino acid residues, particularly near the protein cleavage site, implying that that region is nearly identical in all four virions and has an essential role in virion maturation, and also suggests that all four capsid interior surfaces have similar surfaces exposed to the viral RNA. Apart from a small portion of the C promoter, the amino terminus of the BBV protein (residues 1 to 60) is crystallographically disordered and the amino acid residues in that region are not well conserved. The disordered portion of the BBV protein clearly projects from the capsid inner surface into the interior of the virion, the region occupied by the viral RNA. In all four viruses, residues 1 to 60 had a high proportion of basic residues, suggesting a virus-specific interaction of the amino terminus with the virion RNA.     

Comparative restriction endonuclease maps of proviral DNA of the primate type C simian sarcoma-associated virus and gibbon ape leukemia virus group.

Extrachromosomal DNA was purified from canine thymus cells acutely infected with different strains of infectious primate type C viruses of the woolly monkey (simian) sarcoma helper virus and gibbon ape leukemia virus group. All DNA preparations contained linear proviral molecules of 9.1 to 9.2 kilobases, at least some of which represent complete infectious proviral DNA. Cells infected with a replication-defective fibroblast-transforming sarcoma virus and its helper, a replication-competent nontransforming helper virus, also contained a 6.6- to 6.7-kilobase DNA. These proviral DNA molecules were digested with different restriction endonucleases, and the resultant fragments were oriented to the viral RNA by a combination of partial digestions, codigestion with more than one endonuclease, digestion of integrated proviral DNA, and hybridization with 3'- and 5'-specific viral probes. The 3'- and 5'-specific probes each hybridized to fragments from both ends of proviral DNA, indicating that, in common with those of other retroviruses, these proviruses contain a large terminal redundancy at both ends, each of which consists of sequences derived from both the 3' and 5' regions of the viral RNA. The proviral sequences are organized 3',5'-unique-3',5'. Four restriction enzymes (KpnI, SmaI, PstI, and SstI) recognized sites within the large terminal redundancies, and these sites were conserved within all the isolates tested. This suggests that both the 3' and 5' ends of the genomic RNA of these viruses are extremely closely related. In contrast, the restriction sites within the unique portion of the provirus were not strongly conserved within this group of viruses, even though they were related along most of their genomes. Whereas the 5' 60 to 70% of the RNA of these viruses was more closely related by liquid hybridization experiments than was the 3' 30 to 40%, restriction sites within this region were not preferentially conserved, suggesting that small sequence differences or point mutations or both exist throughout the entire unique portion of the genome among these viruses.     

Discontinuous sequence change of human immunodeficiency virus (HIV) type 1 env sequences in plasma viral and lymphocyte-associated proviral populations in vivo: implications for models of HIV pathogenesis.

Sequence change in different hypervariable regions of the external membrane glycoprotein (gp120) of human immunodeficiency virus type 1 (HIV-1) was studied. Viral RNA associated with cell-free virus particles circulating in plasma and proviral DNA present in HIV-infected peripheral blood mononuclear cells (PBMCs) were extracted from blood samples of two currently asymptomatic hemophiliac patients over a 5-year period. HIV sequences were amplified by polymerase chain reaction to allow analysis in the V3, V4, and V5 hypervariable regions of gp120. Rapid sequence change, consisting of regular replacements by a succession of distinct viral populations, was found in both plasma virus and PBMC provirus populations. Significant differences between the frequencies of sequence variants in DNA and RNA populations within the same sample were observed, indicating that at any one time point, the predominant plasma virus variants were antigenically distinct from viruses encoded by HIV DNA sequences in PBMCs. How these findings contribute to current models of HIV pathogenesis is discussed.     

Comparative characterization of small RNAs derived from an emaravirus and a geminivirus infecting pigeonpea

High throughput sequencing technologies, supported by bioinformatics tools are employed to retrieve small RNA sequence information derived from the nucleic acids of plant infecting viruses. In addition to characterization of the small RNAs to understand the biology of the virus, the small RNA sequence can be assembled to reconstitute viral genome sequence. For the first time the semiconductor based Ion Proton sequencing technology is used to sequence the small RNAs from pigeonpea (Cajanus cajan) plants infected by two distinct viruses with RNA and DNA as their genomes. The reconstitution of the viral genome sequence revealed that the pigeonpea plant from Kalaburagi (erstwhile Gulbarga, Karnataka state) was infected by an emaravirus species Pigeonpea sterility mosaic emaravirus 1 (PPSMV-1) and another plant from New Delhi was infected by a begomovirus species Mungbean yellow mosaic India virus (MYMIV). Characterization and comparison of small RNA sequences derived from both the viruses showed vast differences in their pattern of accumulation and their size classes. In the case of PPSMV-1, the 21 nt sized siRNAs accumulated at far greater levels followed by 22 and 24 nt siRNAs. Whereas in MYMIV, the proportion of accumulation of each size class of siRNAs was similar. Further the distribution of small RNAs across the genomes of PPSMV-1 and MYMIV was mapped and the density of small RNA accumulation showed a positive correlation with the GC content of viral sequence.     

Top-heavy trophic structure within benthic viral dark matter

A better understanding of system-specific viral ecology in diverse environments is needed to predict patterns of virus–host trophic structure in the Anthropocene. This study characterised viral-host trophic structure within coral reef benthic cyanobacterial mats—a globally proliferating cause and consequence of coral reef degradation. We employed deep longitudinal multi-omic sequencing to characterise the viral assemblage (ssDNA, dsDNA, and dsRNA viruses) and profile lineage-specific host–virus interactions within benthic cyanobacterial mats sampled from Bonaire, Caribbean Netherlands. We recovered 11,012 unique viral populations spanning at least 10 viral families across the orders Caudovirales, Petitvirales, and Mindivirales. Gene-sharing network analyses provided evidence for extensive genomic novelty of mat viruses from reference and environmental viral sequences. Analysis of coverage ratios of viral sequences and computationally predicted hosts spanning 15 phyla and 21 classes revealed virus–host abundance (from DNA) and activity (from RNA) ratios consistently exceeding 1:1, suggesting a top-heavy intra-mat trophic structure with respect to virus–host interactions. Overall, our article contributes a curated database of viral sequences found in Caribbean coral reef benthic cyanobacterial mats (vMAT database) and provides multiple lines of field-based evidence demonstrating that viruses are active members of mat communities, with broader implications for mat functional ecology and demography.     

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.     

Autologous Antibody Capture to Enrich Immunogenic Viruses for Viral Discovery

Discovery of new viruses has been boosted by novel deep sequencing technologies. Currently, many viruses can be identified by sequencing without knowledge of the pathogenicity of the virus. However, attributing the presence of a virus in patient material to a disease in the patient can be a challenge. One approach to meet this challenge is identification of viral sequences based on enrichment by autologous patient antibody capture. This method facilitates identification of viruses that have provoked an immune response within the patient and may increase the sensitivity of the current virus discovery techniques. To demonstrate the utility of this method, virus discovery deep sequencing (VIDISCA-454) was performed on clinical samples from 19 patients: 13 with a known respiratory viral infection and 6 with a known gastrointestinal viral infection. Patient sera was collected from one to several months after the acute infection phase. Input and antibody capture material was sequenced and enrichment was assessed. In 18 of the 19 patients, viral reads from immunogenic viruses were enriched by antibody capture (ranging between 1.5x to 343x in respiratory material, and 1.4x to 53x in stool). Enriched reads were also determined in an identity independent manner by using a novel algorithm Xcompare. In 16 of the 19 patients, 21% to 100% of the enriched reads were derived from infecting viruses. In conclusion, the technique provides a novel approach to specifically identify immunogenic viral sequences among the bulk of sequences which are usually encountered during virus discovery metagenomics.     

Rat sequences of the Kirsten and Harvey murine sarcoma virus genomes: nature, origin, and expression in rat tumor RNA.

Two murine sarcoma viruses, the Kirsten and the Harvey, were isolated by passage of mouse type C leukemia viruses through rats. These sarcoma viruses have genomes containing portions of their parental type C mouse leukemia virus genomes, in stable association with specific rat cellular sequences that we find to be quite likely not those of a rat type C leukemia virus. To determine if these murine sarcoma viruses provide a model relevant to the events occurring in spontaneous tumors, we have hybridized DNA and RNA prepared from rat tumors and normal rat tissues to [3H]DNA prepared from the Kirsten murine sarcoma virus. We have also hybridized these rat tissue nucleic acids to [3H]DNA prepared from a respresentative endogenous rat type C leukemia virus, the WFU (Wistar-Furth). Sarcoma-viral rat cellular sequences and endogenous rat leukemia viral sequences were detected in the DNA of both tumor and normal tissues, with no evidence of either gene amplification or additional sequences being present in tumor DNA. Sarcoma-viral rat cellular sequences and endogenous rat leukemia viral sequences were detected at elevated concentrations in the RNA of many rat tumors and in specific groups of normal tissues.     

New DNA viruses identified in patients with acute viral infection syndrome

A sequence-independent PCR amplification method was used to identify viral nucleic acids in the plasma samples of 25 individuals presenting with symptoms of acute viral infection following high-risk behavior for human immunodeficiency virus type 1 transmission. GB virus C/hepatitis G virus was identified in three individuals and hepatitis B virus in one individual. Three previously undescribed DNA viruses were also detected, a parvovirus and two viruses related to TT virus (TTV). Nucleic acids in human plasma that were distantly related to bacterial sequences or with no detectable similarities to known sequences were also found. Nearly complete viral genome sequencing and phylogenetic analysis confirmed the presence of a new parvovirus distinct from known human and animal parvoviruses and of two related TTV-like viruses highly divergent from both the TTV and TTV-like minivirus groups. The detection of two previously undescribed viral species in a small group of individuals presenting acute viral syndrome with unknown etiology indicates that a rich yield of new human viruses may be readily identifiable using simple methods of sequence-independent nucleic acid amplification and limited sequencing.