Full Genome Virus Detection in Fecal Samples Using Sensitive Nucleic Acid Preparation,Deep Sequencing,and a Novel Iterative Sequence Classification Algorithm

We have developed a full genome virus detection process that combines sensitive nucleic acid preparation optimised for virus identification in fecal material with Illumina MiSeq sequencing and a novel post-sequencing virus identification algorithm. Enriched viral nucleic acid was converted to double-stranded DNA and subjected to Illumina MiSeq sequencing. The resulting short reads were processed with a novel iterative Python algorithm SLIM for the identification of sequences with homology to known viruses. De novo assembly was then used to generate full viral genomes. The sensitivity of this process was demonstrated with a set of fecal samples from HIV-1 infected patients. A quantitative assessment of the mammalian, plant, and bacterial virus content of this compartment was generated and the deep sequencing data were sufficient to assembly 12 complete viral genomes from 6 virus families. The method detected high levels of enteropathic viruses that are normally controlled in healthy adults, but may be involved in the pathogenesis of HIV-1 infection and will provide a powerful tool for virus detection and for analyzing changes in the fecal virome associated with HIV-1 progression and pathogenesis.     

Rapid Real-Time PCR Assay for Detection and Quantitation of Mycobacterium avium subsp. paratuberculosis DNA in Artificially Contaminated Milk

Using fluorescence resonance energy transfer technology and Lightcycler analysis, we developed a real-time PCR assay with primers and probes designed by using IS900 which allowed rapid detection of Mycobacterium avium subsp. paratuberculosis DNA in artificially contaminated milk. Initially, the PCR parameters (including primer and probe levels, assay volume, Mg²⁺ concentration, and annealing temperature) were optimized. Subsequently, the quantitative ability of the assay was tested and was found to be accurate over a broad linear range (3 × 10⁶ to 3 × 10¹ copies). The assay sensitivity when purified DNA was used was determined to be as low as five copies, with excellent reproducibility. A range of DNA isolation strategies was developed for isolating M. avium subsp. paratuberculosis DNA from spiked milk, the most effective of which involved the use of 50 mM Tris HCl, 10 mM EDTA, 2% Triton X-100, 4 M guanidinium isothiocyante, and 0.3 M sodium acetate combined with boiling, physical grinding, and nucleic acid spin columns. When this technique was used in conjunction with the real-time PCR assay, it was possible to consistently detect <100 organisms per ml of milk (equivalent to 2,000 organisms per 25 ml). Furthermore, the entire procedure (extraction and PCR) was performed in less than 3 h and was successfully adapted to quantify M. avium subsp. paratuberculosis in spiked milk from heavily and mildly contaminated samples.     

Multiple sources of contamination in samples from patients reported to have XMRV infection

Xenotropic murine leukemia virus (MLV)-related retrovirus (XMRV) was reported to be associated with prostate cancer by Urisman, et al. in 2006 and chronic fatigue syndrome (CFS) by Lombardi, et al. in 2009. To investigate this association, we independently evaluated plasma samples from 4 patients with CFS reported by Lombardi, et al. to have XMRV infection and from 5 healthy controls reported to be XMRV uninfected. We also analyzed viral sequences obtained from supernatants of cell cultures found to contain XMRV after coculture with 9 clinical samples from 8 patients. A qPCR assay capable of distinguishing XMRV from endogenous MLVs showed that the viral sequences detected in the CFS patient plasma behaved like endogenous MLVs and not XMRV. Single-genome sequences (N = 89) from CFS patient plasma were indistinguishable from endogenous MLVs found in the mouse genome that are distinct from XMRV. By contrast, XMRV sequences were detected by qPCR in 2 of the 5 plasma samples from healthy controls (sequencing of the qPCR product confirmed XMRV not MLV). Single-genome sequences (N = 234) from the 9 culture supernatants reportedly positive for XMRV were indistinguishable from XMRV sequences obtained from 22Rv1 and XMRV-contaminated 293T cell-lines. These results indicate that MLV DNA detected in the plasma samples from CFS patients evaluated in this study was from contaminating mouse genomic DNA and that XMRV detected in plasma samples from healthy controls and in cultures of patient samples was due to cross-contamination with XMRV (virus or nucleic acid).     

Polyoma virus defective DNAs. II. Physical map of a molecule with rearranged and reiterated sequences (D74).

The structure of the polyoma virus defective species D74 (74% the size of full-length polyoma virus DNA) has been determined and compared with that of polyoma virus A2 DNA. D74 appears to be composed entirely of viral DNA sequences. (No host DNA sequences have been detected.) It is made up of three DNA segments, each about 24, 24 and 27% in size. The two 24% segments appear to be identical and the 27% segment contains one copy of all the sequences found in the 24% fragments as well as a duplication of some of the sequences. When related to the physical map of A2 DNA, each segment is found to be composed of viral sequences from 1 to about 19 map units, 67 to 69 map units and 70 to 72 map units.Three features found in other polyoma virus defective species (Lund et al., 1977) are also present in D74. (1) Sequences from the region around 67 map units are linked to other (non-contiguous) viral sequences. (2) Sequences at about 72 map units are linked to sequences at 1 map unit. (3) Multiple copies of sequences from around the origin of viral DNA replication are present. From studies on other polyoma defective molecules (Griffin &; Fried, 1975; Lund et al., 1977), the origin of DNA replication for polyoma virus has been defined to lie within the sequences from 67 to 72 map units. Since D74 replicates efficiently but lacks the sequences between 69 to 70 map units, the origin of DNA replication appears to be further defined as lying within 67 and 69 map units and/or 70 to 72 map units.     

High-fidelity PCR amplification of infectious copies of the complete simian virus 40 genome from plasmids and virus-infected cell lysates

We describe here a long-polymerase chain reaction (PCR) method that can be used to amplify complete simian virus 40 (SV40) DNA with high fidelity, and we show that authentic, viable virus can be produced from molecular clones of the PCR-amplified viral DNAs. A commercial long-PCR kit that employed a combination of Taq and GB-D polymerases was used, together with a pair of overlapping primers that recognized a unique EcoRI site in the SV40 genome. Efficient amplification required linearization of the circular SV40 genomic DNAs with EcoRI. Entire SV40 genomes were successfully PCR-amplified from an SV40 plasmid and from two different SV40-infected cell lysates and were cloned into pUC-19. Three separate segments of the cloned viral genomes were DNA sequenced, and no nucleotide changes relative to the parental virus were detected, suggesting that the viral DNAs had been amplified with high fidelity. Each PCR clone was infectious, and no differences were detected in the growth characteristics of viruses derived from these clones as compared to the original viral strain. The procedure we utilized shortens and simplifies the molecular cloning of small double-stranded DNA viruses and will be useful for viral diagnostic tests and for recovery of virus from clinical samples. The results of these experiments have broad implications, as the methodology is applicable to many systems.     

Polyoma virus defective DNAs. I. Physical maps of a related set of defective molecules (D76, D91, D92).

Three related polyoma virus species, designated D92 (92% the size of full-length polyoma virus DNA), D91 (91%) and D76 (76%) have been analysed and their structures compared with that of polyoma virus A2 DNA. Three independent methods (restriction endonuclease cleavage, depurination fingerprinting and DNA-DNA hybridization) were used in the analysis.The defective DNAs appear to be: (1) entirely composed of viral sequences (no host DNA sequences were detected): (2) made up in part of long continuous sequences of DNA which appear identical to sequences of A2 DNA (D92 contains continuous sequences from 1 to 72 map units on the physical map of A2 DNA; that is, it contains the entire late region and part of the early region of the viral DNA. D91 and D76 contain those same sequences except for a 1% deletion around 18 map units): (3) made up in part of rearranged viral sequences.Several interesting features were noted about the rearranged sequences present in the defective DNAs. Sequences from the region around 67 map units were found linked to other (non-contiguous) regions of the DNA. Sequences from about 72 map units were linked to sequences from about 1 map unit. Multiple copies of sequences from 67 to 72 map units (from around the origin of DNA replication) were found (4 copies in D91 and D92, and 2 copies in D76).     

MicroRNAs and Unusual Small RNAs Discovered in Kaposi's Sarcoma-Associated Herpesvirus Virions

It is widely held that any given virus uses only one type of nucleic acid for genetic information storage. However, this consensus has been challenged slightly by several recent studies showing that many RNA species are present within a range of DNA viruses that include Kaposi''s sarcoma-associated herpesvirus (KSHV). RNAs extracted from purified DNA virus particles exhibit great diversity in terms of length, abundance, temporal expression, cellular localization, and coding capacity during viral infection. In addition to known RNA species, the current study showed that small regulatory RNAs were present in KSHV virions. A large number of viral and cellular microRNAs (miRNAs), as well as unusual small RNAs (usRNAs), were detected in KSHV virions by using deep sequencing. Both viral and host miRNAs detected in small RNAs extracted from KSHV virions were further shown to colocalize with KSHV virions directly by in situ hybridization (ISH)-electron microscopy (EM) (ISH-EM). Some of these miRNAs were differentially present in the host cells and KSHV virions, suggesting that they are not randomly present in KSHV virions. The virional miRNAs could be transported into host cells, and they are biologically functional during de novo viral infection. Our study revealed miRNAs and usRNAs as a novel group of components in KSHV virions.     

Hybridization Selection and In Vitro Translation of Autographa californica Nuclear Polyhedrosis Virus mRNA

We isolated polyadenylated RNA from the cytoplasm of cells infected with Autographa californica nuclear polyhedrosis virus late after infection (21 h postinfection). At that time intracellular protein synthesis was directed almost exclusively toward infected cell-specific proteins. The polyadenylic acid-containing RNA sequences in the cytoplasm at 21 h postinfection were radiolabeled in vitro and hybridized to A. californica nuclear polyhedrosis virus DNA restriction fragments. The polyadenylic acid-containing RNA was derived from regions representing the entire viral genome. Translation in a reticulocyte cell-free protein-synthesizing system of cytoplasmic RNA selected by hybridization to viral DNA and polyadenylic acid-containing RNA produced almost identical polypeptide patterns, suggesting that late after infection almost all of the cytoplasmic polyadenylic acid-containing RNA present in infected cells was of viral origin. Polyhedrin protein (molecular weight, 33,000) and a number of virion structural proteins were among the translation products which were identified by immunoprecipitation and by comparing molecular weights. In addition, some tentative nonstructural infected cell-specific proteins were also detected. Using the hybridization selection technique, we determined that sequences complementary to the message coding for polyhedrin were located on EcoRI fragment I of A. californica nuclear polyhedrosis virus DNA, whereas sequences coding for a putative nonstructural protein (molecular weight, 39,000) were on EcoRI fragment J.     

Characterization of the DNA and structural proteins of entomopoxviruses from Melanoplus sanguinipes,Arphia conspirsa, and Phoetaliotes nebrascensis (Orthoptera)

Entomopoxvirus (EPV) occlusion bodies were isolated from virus infected nymphs of the grasshoppers Melanoplus sanguinipes, Arphia conspirsa, and Phoetaliotes nebrascensis. Separation of the viral structural proteins by sodium dodecyl sulfate-polyacrylamide gel electrophoresis gave unique protein patterns for each of the three viruses. An occlusion body protein of approximately 100,000 MW was isolated from each virus. Cleavage of viral DNA with HinddIII and BamHI restriction endonucleases and separation of the fragments by agarose gel electrophoresis gave different DNA fragment patterns for each of the three entomopoxviruses. Molecular weight estimates of 120 × 10⁶ for M. sanguinipes EPV DNA, 129 × 10⁶ for A. conspirsa EPV DNA, and 125 × 10⁶ for P. nebrascensis EPV DNA were calculated from the sizes of the viral DNA fragments. Approximately 55% base sequence homology was detected by Southern hybridization of α-³²P-labeledM. sanguinipes EPV DNA with P. nebrascensis DNA. No base sequence homology was detected by Southern hybridization of labeled M. sanguinipes EPV DNA to Othnonius batesi EPV DNA (Coleoptera), Amsacta moorei EPV DNA (Lepidoptera), Euxoa auxiliaris EPV DNA (Lepidoptera), and vaccinia virus DNA fragments.     

Campoletis sonorensis Endoparasitic Wasps Contain Forms of C. sonorensis Virus DNA Suggestive of Integrated and Extrachromosomal Polydnavirus DNAs

Campoletis sonorensis virus (CsV) (Polydnaviridae) previously was detected only in the calyx epithelial cells and lumen of the oviducts from female C. sonorensis (Ichneumonidae) endoparasitic wasps (Norton et al., Cell Tissue Res. 162:195-208, 1975). Using dot-blot hybridizations, we detected low amounts of CsV DNA in male and female wasp head and thorax tissues and in male abdominal tissues. Low amounts of extrachromosomal viral DNA were detected in Southern blots of undigested male wasp DNA and in male DNA purified by isopycnic centrifugation. High-molecular-weight male wasp DNA digested with any of several restriction endonucleases and hybridized with cloned viral DNAs from CsV superhelices B and Q under stringent conditions contained CsV-specific DNA fragments that differed significantly in size and number from the hybridizing fragments detected in comparably digested viral DNA. Identical offsize restriction fragments were detected in digested female head and thorax DNA. These data suggest that at least CsV DNAs B and Q are integrated in C. sonorensis cellular DNA and that the virus may be transmitted through the germline.     

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.     

Vaccinia virus Transmission through Experimentally Contaminated Milk Using a Murine Model

Bovine vaccinia (BV) is a zoonosis caused by Vaccinia virus (VACV), which affects dairy cattle and humans. Previous studies have detected the presence of viable virus particles in bovine milk samples naturally and experimentally contaminated with VACV. However, it is not known whether milk contaminated with VACV could be a route of viral transmission. However, anti-Orthopoxvirus antibodies were detected in humans from BV endemic areas, whom had no contact with affected cows, which suggest that other VACV transmission routes are possible, such as consumption of contaminated milk and dairy products. Therefore, it is important to study the possibility of VACV transmission by contaminated milk. This study aimed to examine VACV transmission, pathogenesis and shedding in mice orally inoculated with experimentally contaminated milk. Thirty mice were orally inoculated with milk containing 10⁷ PFU/ml of VACV, and ten mice were orally inoculated with uncontaminated milk. Clinical examinations were performed for 30 consecutive days, and fecal samples and oral swabs (OSs) were collected every other day. Mice were euthanized on predetermined days, and tissue and blood samples were collected. Nested-PCR, plaque reduction neutralization test (PRNT), viral isolation, histopathology, and immunohistochemistry (IHC) methods were performed on the collected samples. No clinical changes were observed in the animals. Viral DNA was detected in feces, blood, OSs and tissues, at least in one of the times tested. The lungs displayed moderate to severe interstitial lymphohistiocytic infiltrates, and only the heart, tonsils, tongue, and stomach did not show immunostaining at the IHC analysis. Neutralizing antibodies were detected at the 20^th and 30^th days post infection in 50% of infected mice. The results revealed that VACV contaminated milk could be a route of viral transmission in mice experimentally infected, showing systemic distribution and shedding through feces and oral mucosa, albeit without exhibiting any clinical signs.     

Nucleic Acid Homology of Murine Type-C Viral Genes

The nucleic acid sequence homology between various murine, endogenous type-C viruses (three host range classes of BALB/c virus, the AT-124 virus, and the CCL 52 virus) and two laboratory strains of murine leukemia virus (Rauscher and Kirsten) was determined by DNA:RNA hybridization. The viral sequences exhibit varying degrees of partial homology. DNA:DNA hybridizations were performed between [³H]DNA probes prepared from N- and X-tropic BALB/c endogenous viruses and cellular DNAs from BALB/c, NIH Swiss, and AKR inbred mouse strains as well as from California feral mice and the Asian mouse subspecies Mus musculus molossinus and M. musculus castaneus. All of these strains of mice are shown to possess multiple (six to seven per haploid genome), partially related copies of type-C virogenes in their DNAs. Thermal melting profiles of the DNA:RNA and DNA:DNA hybrids suggest that the partial homology of the viral nucleic acid sequences is the result of base alterations throughout the viral genomes, rather than the loss of discrete segments of viral sequences.     

Molecular Engineering of the Autographa californica Nuclear Polyhedrosis Virus Genome: Deletion Mutations Within the Polyhedrin Gene

We describe a method to introduce site-specific mutations into the genome of Autographa californica nuclear polyhedrosis virus. Specifically, the A. californica nuclear polyhedrosis virus gene for polyhedrin, the major protein that forms viral occlusions in infected cells, was mutagenized by introducing deletions into the cloned DNA fragment containing the gene. The mutagenized polyhedrin gene was transferred to the intact viral DNA by mixing fragment and viral DNAs, cotransfecting Spodoptera frugiperda cells, and screening for viral recombinants that had undergone allelic exchange. Recombinant viruses with mutant polyhedrin genes were obtained by selecting the progeny virus that did not produce viral occlusions in infected cells (occlusion-negative mutants). Analyses of occlusion-negative mutants demonstrated that the polyhedrin gene was not essential for the production of infectious virus and that deletion of certain sequences within the gene did not alter the control, or decrease the level of expression, of polyhedrin. An early viral protein of 25,000 molecular weight was apparently not essential for virus replication in vitro, as the synthesis of this protein was not detected in cells infected with a mutant virus.     

Loss of Polyoma Virus Infectivity as a Result of a Single Amino Acid Change in a Region of Polyoma Virus Large T-Antigen Which Has Extensive Amino Acid Homology with Simian Virus 40 Large T-Antigen

The polyoma virus (Py) transformed cell line 7axB, selected by in vivo passage of an in vitro transformed cell, contains an integrated tandem array of 2.4 genomes and produces the large, middle, and small Py T-antigen species, with molecular weights of 100,000, 55,000, and 22,000, respectively (Hayday et al., J. Virol. 44:67-77, 1982; Lania et al., Cold Spring Harbor Symp. Quant. Biol. 44:597-603, 1980). The integrated viral and adjacent host DNA sequences have been molecularly cloned as three EcoRI fragments (Hayday et al.). One of these fragments (7B-M), derived from within the tandem viral sequences, is equivalent to an EcoRI viral linear molecule. Fragment 7B-M has been found to be transformation competent but incapable of producing infectious virus after DNA transfection (Hayday et al.). By constructing chimerae between 7B-M and Py DNA and by direct DNA sequencing, the mutation responsible for the loss of infectivity has been located to a single base change (adenine to guanine) at nucleotide 2503. This results in a conversion of an aspartic acid to a glycine in the C-terminal region of the Py large T-antigen but does not appear to affect the binding of the Py large T-antigen to Py DNA at the putative DNA replication and autoregulation binding sites. The mutation is located within a 21-amino acid homology region shared by the simian virus 40 large T-antigen (Friedmann et al., Cell 17:715-724, 1979). These results suggest that the mutation in the 7axB large T-antigen may be involved in the active site of the protein for DNA replication.     

Partial characterization of DNA from five entomopoxviruses

Extensive genomic heterogeneity was detected in the restriction endonuclease cleavage patterns of DNA from five entomopoxvirus isolates and vaccinia virus, strain WR. An 8.2 kilobase pair extra-chromosomal element was detected in Amsacta moorei entomopoxvirus and a 22 kilobase pair extra-chromosomal DNA element was isolated from Choristoneura biennis EPV. The extent of DNA base sequence homology was determined by Southern hybridization of HindIII and BamHI DNA restriction fragments of C. biennis EPV DNA and A. moorei EPV DNA with (α³²P)-labeledA. moorei EPV DNA. Methylation of 5′-C^mCGG-3′ sequences was not detected in the DNA of A. moorei, C. biennis, E. auxiliaris, M. sanguinipes, and A. conspersa entomopoxviruses after cleavage of the viral DNAs with MspI and HpaII restriction endonucleases. Based upon the DNA base sequence homology data presented here, the five entomopoxviruses used in this study appear to be unrelated.     

Construction of an Expression Vector for Production and Purification of Human Somatostatin in Escherichia coli

Cassava mosaic disease, caused by cassava mosaic geminiviruses are transmitted by Bemisia tabaci. The B. tabaci adults from colonies reared on virus free cassava plant produced from apical meristem culture was studied to determine their ability to transmit Indian cassava mosaic virus (ICMV) and Sri Lankan cassava mosaic virus (SLCMV) from cassava to cassava. Virus free plants were confirmed by polymerase chain reaction (PCR) using geminivirus degenerate primers. The virus acquisition access period (AAP) of 48 h on virus infected cassava leaves and 48 h virus inoculation access periods on virus free healthy leaves were investigated. Both ICMV and SLCMV were absolutely transmitted by whiteflies reared on cassava. Virus specific primers were designed in the replicase region and used to detect virus in B. tabaci after different AAP. The PCR amplified replicase genes from virus transmitted cassava leaves were cloned the plasmid DNA was isolated from a recombinant colony of E. coli DH5α after their confirmation by colony PCR and sequenced them. The nucleotide sequences obtained from automated DNA sequencing were confirmed as ICMV and SLCMV replicase gene after homology searching by BLAST and found to be a new isolates. The nucleotide sequences of new isolates were submitted in GenBank (accession number JN652126 and JN595785).