Volume 21, Number 1—January 2015
Protocol for Metagenomic Virus Detection in Clinical Specimens1
Sixty percent of emerging viruses have a zoonotic origin, making transmission from animals a major threat to public health. Prompt identification and analysis of these pathogens are indispensable to taking action toward prevention and protection of the affected population. We quantifiably compared classical and modern approaches of virus purification and enrichment in theory and experiments. Eventually, we established an unbiased protocol for detection of known and novel emerging viruses from organ tissues (tissue-based universal virus detection for viral metagenomics [TUViD-VM]). The final TUViD-VM protocol was extensively validated by using real-time PCR and next-generation sequencing. We could increase the amount of detectable virus nucleic acids and improved the detection of viruses <75,000-fold compared with other tested approaches. This TUViD-VM protocol can be used in metagenomic and virome studies to increase the likelihood of detecting viruses from any biological source.
Viruses responsible for disease outbreaks in humans naturally emerge either from the human population or as zoonoses by transmission from animal hosts (1). Viruses can also emerge unnaturally, either directly (e.g., bioterrorist attacks) or accidentally (e.g., laboratory infections). Despite these possibilities of virus emergence, 60% of emerging viruses have a zoonotic origin, thus highlighting transmission from animals to humans as a major threat to public health (2). Whenever viruses emerge, prompt identification of the agent and implementation of control measures to contain the outbreak are required.
Currently, various next-generation sequencing (NGS) approaches provide solutions for detection of purified and concentrated viruses (i.e., from cell culture). However, for clinical specimens, such as blood, other fluids, or infected organ tissues, successful detection of viruses is less likely because virus-to-host genome ratios are insufficient (3–6). Use of tissues from persons with suspected infections for virus detection enables elucidation of infection directly at the site of viral replication. Although detecting viruses directly from infected organ tissue provides obvious and valuable advantages, reliable purification of viruses directly from tissues still remains a challenge.
In this study, we quantifiably and extensively compared classical and modern experimental approaches for virus purification and enrichment to finalize a protocol for unbiased detection of emerging viruses directly from organ tissues (tissue-based unbiased virus detection for viral metagenomics [TUViD-VM]) for an increased signal-to-noise ratio (ratio of virus genome to host genome) in virus detection. Use of this approach will reduce the amount of host nucleic acids required and save money and time in preparation of samples for NGS and the subsequent bioinformatic analysis.
We first describe how the protocol was developed and evaluated, We then describe the final virus purification and enrichment TUViD-VM protocol for metagenomic deep sequencing for nucleic acid from organ tissue (Figure 1).
All procedures regarding the marmoset used in this study were performed in accordance with the European Association of Zoos and Aquaria Husbandry Guidelines for Callitrichidae, 2nd ed., 2010 (http://www.marmosetcare.com/downloads/EAZA_HusbandryGuidelines.pdf), which promotes the highest possible standard for husbandry of zoo animals. The marmoset was kept in Zoo Heidelberg (Heidelberg, Germany) with other marmosets in a species-appropriate environment enriched with material for occupation and activity and adequate feeding regimens 3 times a day. The marmoset that was euthanized did not have any additional signs of illness or infection. The production of specific pathogen–free eggs (VALO BioMedia GmbH, Osterholz-Scharmbeck, Germany) was performed in accordance with guidelines of the European Pharmacopoeia (EP18.104.22.168.2) and the US Department of Agriculture Veterinary Services (Memorandum 800.65).
All procedures regarding embryonated chicken eggs were based on German Animal Protection Laws. For infection, fertilized chicken eggs at embryonation day 11 were inoculated with virus into the allantois sack or onto the chorioallantoic membrane. Development of embryos was terminated at day 17 of embryonation by cooling the eggs overnight at 4°C. No further specific approval is needed for experiments on embryonated avians before time of hatching. However, additional approval from the internal ethics advisory board of the Robert Koch Institute was obtained and is available on request.
To compare classical and modern experimental approaches of virus purification and enrichment, we designed a tissue model for internal organs of chicken, each infected with 1 of 4 viruses (poxvirus [vaccinia virus], reovirus[orthoreovirus], orthomyxovirus [influenza virus], and paramyxovirus [Sendai virus]) at low concentrations (Table 1; Table 2; Technical Appendix). Viruses were chosen on the basis of their role in emerging zoonotic diseases and their morphologic and molecular heterogeneity to obtain results for a broad range of viruses (Table 3).
Model Tissue and Protocol Development
To establish a model tissue, we inoculated specific pathogen–free embryonated chicken eggs with 1 of the prechosen viruses at different concentrations. A detailed description of egg infection and preparation of the model tissue is shown in the Technical Appendix. Reovirus (T3/Bat/Germany/342/08) (11) was chosen to represent a nonenveloped virus, orthomyxovirus (influenza A PR/8/1934) and paramyxovirus (Sendai virus) were chosen to represent enveloped viruses with an RNA genome, and poxvirus (vaccinia virus) was chosen to represent an enveloped virus with a DNA genome (Table 3). Viruses in this study were selected to optimize detection of viral zoonotic emerging diseases and possible virus bioterrorism agents.
To validate the model tissue homogeneity, we selected every ninth sample for simultaneous RNA/DNA extraction and determined copy numbers for all 4 viruses and the galTBP gene (Figure 2). Samples showed an even Gaussian distribution of virus nucleic acids per aliquot and were considered suitable for subsequent experiments.
To establish a protocol for the purification and detection of unknown viruses from animal tissue, we tested different purification techniques and their combinations, including mechanical, enzymatic, and molecular biological methods; the main aim was to eliminate as much host DNA/RNA and maintain as much virus RNA/DNA as possible to optimize random PCR amplification of unknown viruses. The novel established protocol was tested to detect any virus from lung tissue derived from a New World monkey (marmoset), which had to be euthanized because of the unknown disease-causing agent.
We compared different techniques of virus purification, enrichment, and amplification (detailed description of methods compared is shown in the Technical Appendix). In addition, complex purification techniques (digestion and ultracentrifugation) were compared by conducting experiments that had specific control factors (e.g., ultracentrifugation with different concentrations of sucrose, time and speed) (12). Organization of combinations of different control factors and their variable factors (e.g., concentration levels, duration or speed in orthogonal assays) enables conducting a minimal number of experiments. On the basis of results of all purification techniques, we developed a combined protocol to provide the maximized yield of virus RNA/DNA after purification.
Validation and Analysis of Methods Compared
All compared methods were analyzed simultaneously. Because evaluation of sample quality was ongoing, to exclude any extraction bias, an additional unprocessed control aliquot was extracted and measured with every batch. All results of 1 extraction were rigorously compared with a related control aliquot to normalize any variations caused by extraction, cDNA, and quantitative PCR (qPCR) performance.
Every result was evaluated for increasing the signal-to-noise ratio of virus to host-genome (this ratio is indicated by ). Given that ΔΔx = Δ measured – Δ control, we assume that the ratio change between virus nucleic acids and host genome is given by ΔΔCt = Δ purified – Δ unprocessed, where Ct is the cycle threshold. To visualize relative quantification (RQ), RQ (2 – ΔΔCt) was plotted against the respective methods. The RQ value indicates the x-fold change compared with that of the control aliquot (e.g., RQ value of 10 means a 10-fold higher Δbetween virus and host genomes compared with the control aliquot) (13). Per definition of the RQ method, the area of significance lays outside RQ values of 0.5 and 2 if the samples show an even Gaussian distribution. Thus, results <0.5 and >2 were considered significant.
An additional scoring system was used to evaluate different methods. For every RQ result that increased the ratio between host and virus nucleic acids, we gave 1 point (maximum +4 points if the method led to better detectability for all 4 viruses), and for every decrease, we subtracted 1 point (minimum is subsequently −4 points). Methods with the highest scores were chosen for establishment of a combined protocol that included purification of unknown viruses from any tissue source (Table 1).
Final TUViD-VM Protocol for the Enrichment and Purification of Viruses from Organ Tissue
For homogenization, a small cube of tissue (0.5–1 cm3) was placed in an autoclaved screw-cap tube (Sarstedt, Hildesheim, Germany) containing 1 mL of phosphate-buffered saline (PBS) buffer and 20–30 sterile ceramic beads. Tissue was disrupted by shaking 4 times at maximum speed at intervals of 15 s by using the FastPrep-24 Instrument (MP Biomedicals, Strasbourg, France). The duration of this procedure was ≈0.5 h.
A total of 200 μL of homogenate was placed in a 1.5-mL tube and vortexed vigorously. The homogenate was centrifuged for 5 min at 2,000 rpm in a bench top centrifuge (Eppendorf, Hamburg, Germany). The supernatant (≈170 μL) was transferred into a clean tube, and the pellet was discarded. The duration of this procedure was ≈0.25 h.
Ultracentrifugation for Virus Particle Separation
A total of 250 μL of 80% (wt/vol) sucrose solution was pipetted into a 2-3/8″ PA ultracentrifuge tube (Beckman Coulter, Krefeld, Germany) and gently overlayed with ≈3 mL of 20% (wt/vol) sucrose solution. The visibility of the phase interface between the 80% and 20% sucrose solutions was checked. The sucrose solution was gently overlayed with cleared tissue supernatant, and PBS was then added to the tubes. The tubes were centrifuged in an SW60 rotor (Beckman Coulter) at 30,000 rpm for 2 h at 4°C. The duration of this procedure was ≈2 h.
Ultracentrifugation to Pellet Virus Particles
The layer on the interface between the 20% and 80% sucrose solutions was collected and transferred into a 3-1/2″ tube (suitable for Beckmann SW32Ti rotors; Beckman Coulter). The collected layer was resuspended in ≈40 mL of PBS and mixed gently by pipetting up and down. The suspension was centrifuged for 1 h at 20,000 rpm and 4°C. The supernatant was then discarded. The duration of this procedure was ≈1 h. As an alternative method, virus particles can be precipitated overnight by using Peg-It (System Biosciences, Mountain View, CA, USA).
The pellet was resuspended in 245 μL of 1× digestion buffer (Turbo DNA Free Kit; Ambion, Darmstadt, Germany). A total of Add 5 μL of Turbo DNase (Turbo DNA Free Kit: Ambion) was added and incubated for 30 min at 37°C. The suspension was transferred to a 1.5-mL reaction tube. A total of 10 μL of stop reagent (Turbo DNA Free Kit; Ambion) wad added, incubated at room temperature for 1 min, and centrifuged at 2,000 rpm for 3 min. The supernatant was transferred to another tube, and pellet was discarded. The duration of this procedure was ≈0.75 h.
Combined TRIzol LS Extraction
A total of 750 μL of TRIzol LS (Invitrogen Life Technologies, Grand Island, NY, USA) was added to ≈250 μL of supernatant from previous procedures and homogenized by pipetting up and down 10 times. The mixture was incubated for 5 min at room temperature and centrifuged at 12,000 rpm for 10 min. The supernatant was transferred to precentrifuged phase-lock gel tube (5-Prime, Hilden, Germany). A total of 200 μL of chloroform–isoamyl alcohol was added and mixed by inverting the tube vigorously. The tube was incubated for 15 min at room temperature and centrifuged at 12,000 rpm for 15 min.
Approximately 280 μL of supernatant from the phase-lock gel tube was transferred to another tube containing 1,120 μL of AVL lysis buffer without carrier RNA (Viral RNA Mini Kit; QIAGEN, Hilden, Germany). A total of 700 μL of absolute ethanol was added and mixed by pulse vortexing. The solution was transferred in 600-μL portions to a QIAamp Mini Column, QIAGEN), centrifuged 8,000 rpm for 1 min, and the filtrate was discard. The column was placed in a new collection tube, loaded again, and centrifuged until the lysate was added to the column. A total of 500 μL of 70% (wt/vol) ethanol was added and the column was centrifuged at 8,000 rpm for 3 min.
A mixture of 10 μL of DNase and 70 µL of RDD buffer (RNase-Free DNase Set; QIAGEN) was added to the column and incubated for 15 min at room temperature, as described by the manufacturer. The column was washed with 500 μL of AW1 buffer, centrifuged at 8,000 rpm for 1 min, and the filtrate was discarded. The column was placed in a new tube, 500 μL of AW2 buffer was added, the tube was centrifuged at maximum speed for 3 min, and the filtrate was discarded. The column was then placed in a new tube, and the tube was centrifuged at maximum speed for 1 min to dry the column. A total of 30 μL of elution buffer was added to the column, incubated for 5 min at room temperature, and the column was centrifuged in a new 1.5-mL tube. A total of 30 μL of elution buffer was added to the column, incubated for 5 min at room temperature, and centrifuged in the same tube. RNA (≈60 μL) was chilled on ice. The duration of this procedure was ≈3 h.
Single-stranded cDNA was produced by using the Reverse Transcription Reagent Kit (Applied Biosystems, Foster City, CA, USA) and adapted for a 50-μL reaction containing 30 μL of RNA, 2 μL (40 μmol/L) of K8N random primer (7), 3.2 μL (25 mmol/L) of dNTPs, 4 μL 10× buffer, 9 μL (50 mmol/L) of MgCl2, 0.8 μL of RNase inhibitor, 0.6 μL of reverse transcriptase, and 0.4 μL of water). A total of 2 μL of K8N random primers and 3.2 μL of dNTPs were added to the 30 μL of RNA and heated at 95°C for 5 min before quenching on ice. The remaining contents of the mixture were heated at 42°C for 60 min before the enzyme was inactivated at 95°C for 10 min.
Double-stranded cDNA was produced by mixing 2 μL of K8N random primers, 3 μL of Klenow buffer (New England Biolabs, Ipswich, MA, USA), and 2 μL (2.5 mmol/L) of dNTPs with 19 μL of cDNA. The mixture was heated at 95°C for 2 min and cooled to 4°C. A total of 1.67 μL of Klenow fragment (New England Biolabs) was added and the mixture was at 37°C for 60 min. Double-stranded cDNA was purified by using the MSB Spin PCRapace Purification Kit (Invitek, Berlin, Germany) and an elution volume of 30 μL. Random amplification was performed by using the procedures reported by Stang and Korn (7). Successful random amplification (a 200–2,000-bp smear) was visualized by agarose gel electrophoresis of 10 μL of PCR product. The duration of this procedure was ≈4.5 h. Sequence information can be obtained by either cloning into sequencing vectors or by NGS.
RNA samples were fragmented by using the NEBNext Magnesium RNA Fragmentation Module (New England Biolabs). RNA was purified by using RNeasy MinElute (QIAGEN). For cDNA synthesis, Superscript II and Murine RNAse inhibitor (New England Biolabs) were used. Second-strand synthesis was performed by using the NEBNext mRNA Second Strand Synthesis Module (New England Biolabs) and purified by using the MinElute PCR Purification Kit (QIAGEN).
Double-stranded cDNA, DNA, and random PCR products were quantified by using the Qubit HS dsDNA Kit (Invitrogen Life Technologies). Sequencing libraries were established by Ion Xpress Plus Fragment Library Kit (without chemical fragmentation) with indices (Ion Xpress Barcode Adapters 1–16 Kit). The sequencing library was then amplified by using an emulsion-based clonal amplification PCR in the Ion OneTouch 200 Template v2 DL Kit and enriched by using an Ion OneTouch Enrichment System. Sequencing was performed on an IonTorrent PGM in the Ion PGM Sequencing 300 Kit with the Ion 318 Chip Kit (Invitrogen Life Technologies).
NGS Data Analysis
Programs used for sequence analysis were Geneious Pro R6 (Biomatters, Auckland, New Zealand) and Bowtie2align (14). The percentage of bases (Q>20) was ≈80% before length filtering (100–1,000 nt) was applied to remove shorter reads. No additional quality trimming was applied because the quality average was sufficient for our approach. Remaining reads were mapped to the whole reference genomes (or all segments of reference genome) by using Bowtie2align for paramyxovirus (Sendai virus strain Tianjin; GenBank accession no. EF679198), reovirus (T3/Bat/Germany/342/08, 10 segments; JQ412755-JQ412764), orthomyxovirus (influenza H1N1 strain A/Puerto Rico 8-SV14/1934, 8 segments; CY040170-CY040177), and poxvirus (vaccinia virus strain WR, no. AY243312). Coverage of genomes was calculated in weighted average for segmented genomes.
Development of Protocol
Every step of the TUViD protocol (homogenization of tissue, filtration, digestion, enrichment, extraction, and random amplification) was compared with alternative approaches. Results are shown in Figures 3,4,5,6,7. Each approach was tested with individual samples, which were measured by using 5 PCRs specific for viruses used and host background in 2 replicates (10 reactions/sample): Results were quantified and evaluated in qPCRs for the 4 viruses and presence of host nucleic acids (Technical Appendix; Table 4; Figures 3,4,5,6,7). A scoring system was developed to assess the optimal combination of all 4 viruses (Table 1; Figures 3,4,5,6,7). A preliminary protocol was further validated and adjusted until no host nucleic acids were detectable by qPCR. This protocol maximized the amount of amplified virus nucleic acids. Subsequently, we established an unbiased protocol for the detection of known and novel viruses in infected organ tissues (TUViD-VM).
TUViD-VM Validation by NGS
The TUViD-VM protocol was validated by NGS of 4 aliquots of the model tissue. One aliquot was prepared by using the TUViD-VM protocol developed in this study, and 3 aliquots were prepared by using other approaches commonly used for unbiased virus detection (Figure 8; Technical Appendix). We chose the Invitrogen Life Technologies platform because of its rapid run time and read length, which are crucial for diagnostic purposes. All independent runs were normalized to 1,000,000 output reads for reliable comparison (Table 5; Figure 8). NGS results confirmed the substantial increase in virus nucleic acids, as well as the decrease of host nucleic acids achieved by purification with the novel protocol. The amount of detectable virus nucleic acids was increased >1,000-fold compared with other NGS approaches (Figure 8). For example, although the best NGS approach delivered 40 reads for paramyxovirus in infected chicken tissue, the TUViD-VM protocol resulted in >60,000 reads (97.80% coverage of the complete genome) (Figures 8,9,10; Table 5).
To provide a proof of concept, we prepared lung tissue from the marmoset that was euthanized and had a natural respiratory infection with Sendai virus by using the 4 approaches and sequenced by using the Invitrogen Life Technologies protocol. Using the TUViD-VM protocol, we found that the amount of detectable virus in marmoset tissue increased 75,000-fold compared with that for other NGS approaches (>400,000 Sendai virus reads compared with 6), which represented 99.98% coverage of the Sendai virus genome and ≈50% of the total read output (Figures 8,10; Table 5).
In this study, we successfully established a purification and enrichment protocol, which shows rapid and reliable results, for detection of known and novel viruses in tissues. Likelihood of detection of RNA viruses was increased. In addition, detection of DNA incorporated in virus particles was not affected even though DNA digestion was performed. The cutoff sensitivity was 100–1,000 virus copies/mL of homogenized organ material (e.g., reovirus; Table 5). The cutoff sensitivity of compared approaches was ≥106 virus copies/mL. The TUViD-VM protocol (from solid tissue sampling to nucleic acid preparation for NGS) takes 12 h to complete. If one allows 16 h for NGS, the TUViD-VM protocol provides sequence data output within 28 h.
Current NGS techniques used for metagenomic approaches produce large amounts of sequence data, which might increase the likelihood of detection of diminutive amounts of virus in comparison with the host genome. The only limiting factor seems to be the cost required for processing 1 sample and capacities for computational analysis of results. This in silico analysis should increase the signal-to-noise ratio of relevant sequences by subtracting nonrelevant sequences, such as the host genome. However, genome sequence data for mammals are limited; only 23 sequences (0.4%) for 5,487 species (18). Just 3 genome sequences are available for bats, although they are the second most abundant mammalian species (exceeded only by rodents). There are >1,100 species of bats worldwide and they are suspected vectors of pathogenic viruses (e.g., Ebola virus, Nipah virus, Hendra virus, lyssavirus, and severe acute respiratory syndrome coronavirus). Thus, it seems inefficient to invest large amounts of time, money, and effort in obtaining large datasets, only to invest even more resources to categorize them. Furthermore, quantitative comparison of the virus-enrichment strategies described enables evaluation of multiple classical and modern approaches.
The TUViD-VM described protocol increases the signal-to-noise ratio by as much as 75,000-fold than that for compared approaches and can detect virus genomes quickly in infected tissues (Figures 9,10). Although sequencing of nucleic acid from relatively pure sources (e.g., cell culture, allantoic fluids) is well established and results in reasonable output (11,19,20), sequencing of nucleic acid clinical specimens is still challenging. Other studies reported 0.1% to <10% mammalian virus reads from clinical samples, such as tissue, guano, feces, and pharyngeal swab specimens (3,19,21–24). A method reported by Daly at al. showed promising results for detection of DNA viruses but lacked similar results for detection of RNA viruses (25). In contrast, our protocol resulted in up to 45% mammalian RNA virus reads directly from infected organ tissue (Figure 8).
After its successful and extensive validation, we highly recommend this protocol for investigation of outbreaks with unknown viral etiologic agents in humans and animals. Furthermore, this protocol can be used in metagenomic virome studies and will be beneficial whenever library construction is necessary (i.e., molecular cloning and NGS) to increase detection likelihood for viruses from any biological source. This protocol would be particularly useful for increasing the signal-to-noise ratio in virus analysis of biological samples in which levels of background nucleic acids are high, which result in difficulties in virus detection and identification. Thus, the TUViD-VM protocol described greatly increases the likelihood of detecting viruses during outbreaks of emerging infectious diseases and in metagenomic virus detection studies.
Dr Kohl is a scientist at the Centre for Biological Threats and Special Pathogens, Robert Koch Institute, Berlin Germany. Her research interests are the detection of emerging and reemerging viruses and the characterization of novel pathogens.
We thank Ute Kramer for providing technical assistance, Brunhilde Schweiger for providing influenza A virus strain PR8/38, Marc Hoferer for providing infected marmoset tissues, and Ina Smith and Ursula Erikli for copyediting this article.
This study was partially supported by a grant from the Konrad-Adenauer-Foundation and a scholarship to C.K.
- Fauci AS, Morens DM. The perpetual challenge of infectious diseases. N Engl J Med. 2012;366:454–61.
- Jones KE, Patel NGN, Levy MA, Storeygard A, Balk D, Gittleman JL, Global trends in emerging infectious diseases. Nature. 2008;451:990–3.
- He B, Li Z, Yang F, Zheng J, Feng Y, Guo H, Virome profiling of bats from Myanmar by metagenomic analysis of tissue samples reveals more novel mammalian viruses. PLoS ONE. 2013;8:e61950.
- Baker KS, Leggett RM, Bexfield NH, Alston M, Daly G, Todd S, Metagenomic study of the viruses of African straw-colored fruit bats: detection of a chiropteran poxvirus and isolation of a novel adenovirus. Virology. 2013;441:95–106.
- Nakamura S, Yang C-S, Sakon N, Ueda M, Tougan T, Yamashita A, Direct metagenomic detection of viral pathogens in nasal and fecal specimens using an unbiased high-throughput sequencing approach. PLoS ONE. 2009;4:e4219.
- Whon TW, Kim M-S, Roh SW, Shin N-R, Lee H-W, Bae J-W. Metagenomic characterization of airborne viral DNA diversity in the near-surface atmosphere. J Virol. 2012;86:8221–31.
- Stang A, Korn K. Characterization of virus isolates by particle-associated nucleic acid PCR. J Clin Microbiol. 2005;43:716–20.
- Clem AL, Sims J, Telang S, Eaton JW, Chesney J. Virus detection and identification using random multiplex (RT)-PCR with 3′-locked random primers. Virol J. 2007;4:65.
- King AM, Adams MJ, Carstens EB, Ninth L. Report of the international committee on taxonomy of viruses. San Diego: Elsevier Academic Press; 2011.
- Tidona CA, Darai G. The Springer index of viruses. Berlin: Springer; 2001.
- Kohl C, Lesnik R, Brinkmann A, Ebinger A, Radonić A, Nitsche A, Isolation and characterization of three mammalian orthoreoviruses from European bats. PLoS ONE. 2012;7:e43106.
- Taguchi G. Introduction to quality engineering : designing quality into products and processes. Tokyo: Asian Productivity Organization; 1986.
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
- Schulze M, Nitsche A, Schweiger B, Biere B. Diagnostic approach for the differentiation of the pandemic influenza A(H1N1)v virus from recent human influenza viruses by real-time PCR. PLoS ONE. 2010;5:e9966.
- Schroeder K, Nitsche A. Multicolor, multiplex real-time PCR assay for the detection of human-pathogenic poxviruses. Mol Cell Probes. 2010;24:110–3.
- Kurth A. Possible biohazard risk from infectious tissue and culture cells preserved with RNAlater. Clin Chem. 2007;53:1389–90.
- Schipper J, Chanson JS, Chiozza F, Cox NA, Hoffmann M, Katariya V, The status of the world’s land and marine mammals: diversity, threat, and knowledge. Science. 2008;322:225–30.
- Marston DA, McElhinney LM, Ellis RJ, Horton DL, Wise EL, Leech SL, Next generation sequencing of viral RNA genomes. BMC Genomics. 2013;14:444.
- Kohl C, Vidovszky MZ, Mühldorfer K, Dabrowski PW, Radonić A, Nitsche A, Genome analysis of bat adenovirus 2: indications of interspecies transmission. J Virol. 2012;86:1888–92.
- Wu Z, Ren X, Yang L, Hu Y, Yang J, He G, Virome analysis for identification of novel mammalian viruses in bat species from Chinese provinces. J Virol. 2012;86:10999–1012.
- Ge X, Li Y, Yang X, Zhang H, Zhou P, Zhang Y, Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China. J Virol. 2012;86:4620–30.
- Donaldson EF, Haskew ANA, Gates JE, Huynh J, Moore CJ, Frieman MB. Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat. J Virol. 2010;84:13004–18.
- Li L, Victoria JG, Wang C, Jones M, Fellers GM, Kunz TH, Bat guano virome: predominance of dietary viruses from insects and plants plus novel mammalian viruses. J Virol. 2010;84:6955–65.
- Daly GM, Bexfield N, Heaney J, Stubbs S, Mayer AP, Palser A, A viral discovery methodology for clinical biopsy samples utilizing massively parallel next generation sequencing. PLoS ONE. 2011;6:e28879.
TablesCite This Article
1Preliminary results of this study were presented at the Biodefense and Emerging Infectious Diseases Meeting, January 29, 2014, Washington DC, USA.