Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Applied and Environmental Science

Multiplexed Metagenomic Deep Sequencing To Analyze the Composition of High-Priority Pathogen Reagents

Michael R. Wilson, Greg Fedewa, Mark D. Stenglein, Judith Olejnik, Linda J. Rennick, Sham Nambulli, Friederike Feldmann, W. Paul Duprex, John H. Connor, Elke Mühlberger, Joseph L. DeRisi
Nicholas Chia, Editor
Michael R. Wilson
aDepartment of Neurology, University of California, San Francisco, San Francisco, California, USA
bDepartment of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Greg Fedewa
cIntegrative Program in Quantitative Biology, Bioinformatics, University of California, San Francisco, San Francisco, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Greg Fedewa
Mark D. Stenglein
dDepartment of Microbiology, Immunology, and Pathology, Colorado State University, Fort Collins, Colorado, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Judith Olejnik
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Linda J. Rennick
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sham Nambulli
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Friederike Feldmann
gRocky Mountain Veterinary Branch, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, Montana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
W. Paul Duprex
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John H. Connor
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John H. Connor
Elke Mühlberger
eDepartment of Microbiology, Boston University School of Medicine, Boston, Massachusetts, USA
fNational Emerging Infectious Diseases Laboratories, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joseph L. DeRisi
bDepartment of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, California, USA
hHoward Hughes Medical Institute, Chevy Chase, Maryland, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicholas Chia
Mayo Clinic
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00058-16
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Both the integrity and reproducibility of experiments using select agents depend in large part on unbiased validation to ensure the correct identity and purity of the species in question. Metagenomic deep sequencing (MDS) provides the required level of validation by allowing for an unbiased and comprehensive assessment of all the microbes in a laboratory stock. Laboratories studying high-priority pathogens need comprehensive methods to confirm microbial species and strains while also detecting contamination. Metagenomic deep sequencing (MDS) inventories nucleic acids present in laboratory stocks, providing an unbiased assessment of pathogen identity, the extent of genomic variation, and the presence of contaminants. Double-stranded cDNA MDS libraries were constructed from RNA extracted from in vitro-passaged stocks of six viruses (La Crosse virus, Ebola virus, canine distemper virus, measles virus, human respiratory syncytial virus, and vesicular stomatitis virus). Each library was dual indexed and pooled for sequencing. A custom bioinformatics pipeline determined the organisms present in each sample in a blinded fashion. Single nucleotide variant (SNV) analysis identified viral isolates. We confirmed that (i) each sample contained the expected microbe, (ii) dual indexing of the samples minimized false assignments of individual sequences, (iii) multiple viral and bacterial contaminants were present, and (iv) SNV analysis of the viral genomes allowed precise identification of the viral isolates. MDS can be multiplexed to allow simultaneous and unbiased interrogation of mixed microbial cultures and (i) confirm pathogen identity, (ii) characterize the extent of genomic variation, (iii) confirm the cell line used for virus propagation, and (iv) assess for contaminating microbes. These assessments ensure the true composition of these high-priority reagents and generate a comprehensive database of microbial genomes studied in each facility. MDS can serve as an integral part of a pathogen-tracking program which in turn will enhance sample security and increase experimental rigor and precision. IMPORTANCE Both the integrity and reproducibility of experiments using select agents depend in large part on unbiased validation to ensure the correct identity and purity of the species in question. Metagenomic deep sequencing (MDS) provides the required level of validation by allowing for an unbiased and comprehensive assessment of all the microbes in a laboratory stock.

INTRODUCTION

Virology laboratories must maintain constant vigilance regarding the provenance of strains in their collections. In addition to confirming that experiments are being performed on the intended virus, unbiased methods that can detect viral coinfections and/or other microbial contaminants are critical for ensuring that experiments are faithfully reporting on the biology of a particular virus and not on a polymicrobial exposure. Stock assurance is exceptionally important when working with viruses which require a high level of biocontainment, given the increased security requirements and higher costs and higher stakes of the research.

The validity and reproducibility of experiments, in addition to security and tracking concerns, demand an approach that can provide not only simple viral identification but also the entire spectrum of single nucleotide variants (SNVs) and their respective frequencies within a given stock. Molecular subtyping was critical in determining the source of the 2001 bioterrorism-associated anthrax outbreak (1). Furthermore, the recent missteps in the handling of Bacillus anthracis, influenza virus, and smallpox virus highlight the potential for inadvertent laboratory contamination and the need to improve isolate-tracking methods (2, 3).

In this study, we sequenced virus stocks of six distinct negative-sense RNA viruses, each of them with a unique passage history: (i) mycoplasma-contaminated Ebola virus (EBOV), (ii) mycoplasma-contaminated La Crosse virus (LACV) grown on cells with a known lymphocytic choriomeningitis virus (LCMV) contamination, (iii) vesicular stomatitis virus (VSV) grown on LCMV-contaminated cells, (iv) recombinant human respiratory syncytial virus (rHRSV) expressing enhanced green fluorescent protein (EGFP) from an additional transcription unit (ATU) (4), (v) recombinant canine distemper virus (rCDV) expressing Venus fluorescent protein from an ATU (5), and recombinant measles virus (rMV) expressing EGFP from an ATU (6). The investigators performing the sequence analysis were blinded to the virus isolates and the origins and known contaminations of the virus stocks. Here, we demonstrate the utility of metagenomic deep sequencing (MDS) for comprehensively assessing multiple viral isolates. We show that MDS and custom bioinformatics pipelines utilizing publicly available and free software packages can detect a wide range of bacterial and viral contaminants in high-priority virus stocks and that SNV analysis can successfully discriminate between two closely related EBOV isolates. We also demonstrate how dual-indexed barcodes dramatically decrease the false-positive assignment of sequencing reads to an input sample in multiplexed high-titer virus MDS libraries.

RESULTS

Importance of dual-index barcoding.In many sequencing applications, including MDS, it is critical to be able to distinguish low-level true-positive samples from false positives that result from assignment of sequencing reads to the wrong index (“index bleed-through”). Chimeric molecules produced during amplification steps in library preparation or on the sequencer during cluster generation are particularly susceptible to such misassignment. Creating library molecules with barcode sequences on both ends of the molecule (dual indexing) has been shown to reduce the rate at which this misassignment occurs by requiring concordance between the two barcode sequences (7).

To corroborate the benefit of minimizing misassignment of sequences by using dual indexing in the context of MDS, we analyzed the rate of index bleed-through in an independently generated sequencing data set that was demultiplexed using one or both index reads. We measured the rate of bleed-through in a set of 50 samples that included two samples that were positive for a recently discovered reptarenavirus. The other 48 samples were negative by quantitative reverse transcription PCR (qRT-PCR). We aligned reads from all 50 pooled samples to a virus genome segment present in the two positive samples (8). The data sets contained a median of 2.75 × 106 read pairs, and the two positive samples contained 65,284 and 203,565 virus-mapping reads (Fig. 1). Index bleed-through was determined to have occurred when reads identified as one of the 48 virus-negative samples by barcode actually aligned to this particular reptarenavirus genome segment. The median rate of bleed-through in the data sets that were demultiplexed using a single index was 17.9 per million filtered reads (Fig. 1). The median rate in dual-index demultiplexed data sets was 0.48 per million reads. Thus, in our data set, dual indexing reduced the median rate of read misassignment by 37-fold. The absolute number of mismapping reads decreased from a median of 50 reads to one read per data set.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Dual indexing decreases the median rate of read misassignment by nearly 40-fold. Libraries from 50 samples were pooled and sequenced together. Two of the samples (indicated by arrows) were positive for snake arenavirus, and the other 48 were negative. Data sets were demultiplexed using a single index sequence (black squares) or dual index sequences (red circles), and reads from each data set were mapped with high stringency to the virus sequence. The number of virus mapping reads per million quality-filtered reads is indicated. Some dual-index-demultiplexed data sets had no misassigned reads. In these cases, red circles are not shown.

Identification of microbes.The MDS of total RNA from each sample yielded 150-nucleotide (nt) paired-end sequences (Table 1). The data were analyzed using a rapid computational pipeline developed at the University of California, San Francisco (UCSF), to classify MDS reads and identify potential pathogens by comparison to the entire NCBI nucleotide reference database as described in Materials and Methods (9).

View this table:
  • View inline
  • View popup
TABLE 1

Metagenomic deep sequencing resultsa

The results are summarized in Table 1. The EBOV and LACV samples were found to have significant contamination with Mycoplasma spp., with 3.65 and 0.02 Mycoplasma reads per viral read, respectively. The rMV sample contained more than 900 nonredundant, paired-end Epstein-Barr virus (EBV) sequences (0.013% of total reads) that were expected to come from the EBV-transformed human B-lymphoblastoid cell line (B-LCL) on which the virus was grown (10). Similarly, the rHRSV sample contained 121 paired-end sequences (0.0019% of total reads) that aligned to human papillomavirus 18 (HPV18). This was also unsurprising, since HEp-2 cells are the standard cell line used to culture HRSV and these are known to be a HeLa‑contaminated cell line; HeLa cells were recently shown to have HPV18 DNA integrated into their genome (11). The LACV sample had evidence of Syrian hamster retroviruses, and it was confirmed that the virus had been grown on baby hamster kidney (BHK) cells. The rCDV sample contained sequence reads mapping to fowlpox virus, which had been used as part of the process to generate rCDV from plasmid. Because these RNA extractions were not DNase treated, we cannot rule out that some of the EBV, HPV18, and fowlpox virus sequences were remnants of genomic DNA. Last, the LACV and VSV samples were found to be contaminated with LCMV. None of the isolates had evidence of fungal contamination.

As evidence of the minimal index bleed-through in these dual-indexed samples, only seven viral read pairs unique to LACV were found in the rCDV sample, and no other sample cross-contamination was present in any of the other samples in this set, despite being processed in parallel. Regardless, the fact that low-level index bleed-through may still occur, despite dual indexing, highlights the need for additional layers of error correction. Future iterations of this technology will likely include added features that will further reduce bleed-through of multiplexed samples, such as unique molecular identifier (UMI) barcodes, in which each cDNA molecule is uniquely indexed at the time of first-strand synthesis (12).

SNV experiment.All of the viruses contained several SNVs compared to their reference sequences that ranged in frequency, including some consensus-level SNVs, as listed in Table 2. At a frequency of greater than 0.005, rCDV had 427 SNVs (170 of those being nonsynonymous), EBOV had 143 SNVs (70 nonsynonymous), the LACV L segment had 61 SNVs (45 nonsynonymous), the LACV M segment had 31 SNVs (21 nonsynonymous), the LACV S segment had 6 SNVs (4 nonsynonymous), rMV had 122 SNVs (80 nonsynonymous), rHRSV had 115 SNVs (71 nonsynonymous), and VSV strain Indiana (VSVIN) had 109 SNVs (72 nonsynonymous). A majority of SNVs are rare variants, having a frequency of less than 0.10. All six viruses had nonsynonymous and synonymous SNVs at frequencies of greater than 0.10. Four viruses had nonsynonymous SNVs at a frequency of at least 0.90, while five viruses had synonymous SNVs at a frequency of at least 0.90. Figure 2A displays the distribution of the SNV frequency on each genomic segment.

View this table:
  • View inline
  • View popup
TABLE 2

SNVs in each virus genome segment that are at least 0.5% of populationa

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Distribution of single nucleotide variants (SNVs) for each virus and for two Ebola virus strains. For each plot, the y axis is the log10(SNV frequency) in order to display the range of low-frequency SNVs more accurately. (A) SNV analysis for each virus sample using its reference genome. Abbreviations: CDV, canine distemper virus; EBOV, Ebola virus; LACV, La Crosse virus; MV, measles virus; HRSV, human respiratory syncytial virus; VSV, vesicular stomatitis virus. (B) SNV analysis of 2,501,8050 filtered read-pairs revealed 223 SNVs in the reads that mapped to Ebola virus Mayinga with a frequency of ≥0.90, while there was only one SNV with a frequency of ≥0.90 in the reads that mapped to Ebola virus Kikwit.

The investigators were blinded to the exact isolate of EBOV since the major goal of this study was to identify a state-of-the-art approach which could precisely characterize a virus reagent with isolate specificity in an unbiased fashion. As described in Materials and Methods, SNV analysis was performed to determine the EBOV isolate. Reads were mapped to two different reference genomes: (i) Zaire ebolavirus isolate Ebola virus/Homo sapiens-tc/COD/1976/Yambuku-Mayinga, complete genome strain (GenBank accession no. NC_002549 ), and (ii) Zaire ebolavirus strain Kikwit, complete genome (GenBank accession no. JQ352763 ). Comparing rates of correctly mapped reads did not reveal isolate identity. The same 2,501,8050 filtered read-pairs mapped at virtually identical percentages, 21.03% to EBOV Mayinga and 21.05% to EBOV Kikwit. However, SNV analysis revealed 223 SNVs in the reads that mapped to EBOV Mayinga with a frequency of ≥0.90, while there was only one SNV with a frequency of ≥0.90 in the reads that mapped to EBOV Kikwit (Fig. 2B). Thus, we correctly concluded that the precise isolate sequenced in this experiment was EBOV Kikwit.

DISCUSSION

We describe the successful application of a multiplexed MDS assay and bioinformatics pipelines for comprehensively interrogating the provenance of six virus isolates in a blinded fashion. With regard to multiplexing, the frequency of reads that are assigned an incorrect barcode is a function of the abundance of those reads in the correct barcode library and the degree of multiplexing. So, for high-titer viral stock-derived data sets, where viral reads are abundant, bleed-through occurs at a higher frequency. We demonstrated that dual indexing reduces the misassignment of sequencing reads by 37-fold in a pilot experiment and saw essentially no intersample cross-contamination in our multiplexed assay despite very-high-titer virus isolates. It is important to note that most commercially available dual-indexing kits do not provide this benefit because the reduction in misassignment depends on samples being identified by unique index pairs. Most commercially available kits do not actually provide unique index pairs. Instead, they mix and match indexes to create pairs with repeated index sequences, for instance, by permuting eight i5 indexes with 12 i7 indexes to achieve 96 pairs. We also found that using dual indexing decreased the total number of reads per data set by ~6% due to the removal of unassigned read pairs. However, this is a relatively small tradeoff for the large decrease in misassigned reads.

Contaminating sequences were accurately identified, including bacteria (Mycoplasma spp.) and six viruses across the six virus stocks. Four of these viral contaminating sequences (i.e., EBV, HPV18, hamster gammaretrovirus, and LCMV [n = 2]) reflected the cells on which the viruses were cultured (10, 11, 13, 14). The fowlpox virus present in the rCDV sample was a remnant of the procedure used to generate rCDV, which uses a recombinant fowlpox virus to express T7 RNA polymerase (5). This is unsurprising given the low passage number of the rCDV.

The results from MDS of different virus stocks highlight the utility of this approach for in-depth analysis of virus sequences. There were different contaminating pathogens in each of the stocks, and these biological “fingerprints” provided precise information about their origin and passage history. Such information is invaluable for virologists embarking on time-consuming and expensive studies in biocontainment but is also useful for a broad range of microbiologists who would benefit from sample assurance. For rMV and rHRSV, there was evidence of the cell lines used to propagate the viruses, and for rCDV, there was evidence of the helper virus used for reverse genetics. For viruses with an extensive propagation history, the cell type in which the virus had been passaged could be inferred based on contaminating sequences (i.e., hamster retrovirus sequences) as well as the additional, nonmicrobial sequences that can be assembled to characterize aspects of the host cell genome. This suggests that MDS is useful for virus forensic analysis, helping to identify the manner and cell type in which the virus was cultured or recovered.

The SNV analysis identified numerous synonymous and nonsynonymous mutations present in each of the isolates. The presence and frequency of such mutations in longitudinal samples may be used to monitor inevitable virus adaptation rigorously in the laboratory setting and the quality of seed stocks (15) as well as to define particular isolates (16). Here, SNV analysis allowed the precise determination of the particular EBOV isolate (Fig. 2B). Last, detailed cataloging of SNVs for each isolate allows rapid identification and tracking of high-priority pathogens in the event of an accidental or nonaccidental release of virus into the environment and could be an integral part of any microbial source-tracking program. Viral genome sequencing including SNV analysis has become an important tool to monitor viral spread, viral evolution, and routes of transmission in an outbreak situation as exemplified by the recent EBOV disease outbreak in West Africa (17–26). If a validated reference sequence does not exist, we recommend that the first step be a de novo assembly. However, this should be followed with a remapping of all the reads back to the assembled reference and a LoFreq* or similar SNV analysis, as we have demonstrated here. The latter analysis will reveal the presence of variations or even a mixture of strains that would be obscured or lost through an assemble-and-BLAST strategy alone.

In summary, we show that a single multiplex MDS assay can comprehensively assess (i) virus and isolate identity, (ii) SNVs in virus populations, (iii) the presence of viral coinfection, and (iv) the presence of bacterial contamination and (v) can provide information about the cell line used to propagate the virus.

MATERIALS AND METHODS

Viruses and RNA purification.All work with infectious EBOV was performed under biosafety level 4 (BSL-4) conditions at the Integrated Research Facility, Rocky Mountain Laboratories (RML), Division of Intramural Research, NIAID, NIH, Hamilton, MT. Vero E6 cells (ATCC CRL 1586) were infected with EBOV at a multiplicity of infection (MOI) of 1, and virus-containing supernatants were clarified by low-speed centrifugation at 4 days postinfection (p.i.). Total RNA from supernatant (140 µl) was purified using the Qiagen viral RNA minikit, according to the RML standard operating procedures for virus inactivation. Purified RNA was eluted in H2O. The EBOV stock used for this analysis contained a known mycoplasma contamination.

To generate rMV Khartoum, Sudan (rMVKS) stocks, B‑LCL cells were infected at an MOI of 0.01. After 4 days, when viral cytopathic effect was maximal, the stock was subjected to one freeze-thaw cycle to release the highly cell-associated measles virus into the medium. Cell debris was removed by low-speed centrifugation, and the supernatant was stored at −80°C as virus stock.

To generate rCDV Rhode Island (rCDVRI) stocks, Vero cells expressing the CDV receptor canine CD150 were infected at an MOI of 0.01. After 2 days, when viral cytopathic effect was maximal, the stock was subjected to one freeze-thaw cycle to release cell‑associated virus into the medium. Cell debris was removed by low-speed centrifugation, and the supernatant was stored at −80°C as virus stock.

To generate rHRSVB05 stocks, HEp‑2 cells were infected at an MOI of 0.01. After 3 days, when viral cytopathic effect was maximal, the medium was transferred to 50-ml tubes, and cell debris was removed by low-speed centrifugation. The supernatant was stored at −80°C as virus stock.

The LACV H78 strain was grown on BHK cells (27). BHK cells were infected with passage 3 LACV H78 at an MOI of 0.01. Virus was propagated for 48 h before medium was removed. Cell debris was removed by low-speed centrifugation, and individual aliquots of virus were frozen for future analysis. The LACV stock contained a known mycoplasma contamination.

For growth of the VSV Indiana (VSVIN) strain, BHK cells were infected with VSVIN at an MOI of 0.01. Virus was propagated for 24 h before medium was removed. Cell debris was removed by low-speed centrifugation, and individual aliquots of virus were frozen for future analysis.

For rMV, rCDV, rHRSV, LACV, and VSV, 250 µl of virus-containing supernatant was mixed with 750 µl of Trizol-LS (Ambion), and RNA was isolated according to the manufacturer’s protocol. Purified RNA was resuspended in H2O.

Comparison of single and dual indexing.A library containing a pool of 50 dual-indexed libraries was created as previously described (8). The library was sequenced on an Illumina HiSeq 2500 instrument at the UCSF Center for Advanced Technology. The 150-nt read-length, paired-end, dual-indexed data set was manually demultiplexed using either one or both index reads. To measure the rate of read misassignment, reads were first processed to remove low-quality and adapter sequences as described below. Trimmed reads were then aligned using the Bowtie2 tool to the sequence of a reptarenavirus genome segment (GenBank accession no. KP071661.1 ) actually present in two of the 50 samples (8, 28). Bowtie2 was run with parameters –local –qc-filter –score-min C,120,1. The number of aligning reads for each data set was counted.

Confirming virus identity and detecting contamination.RNA of six individual viral isolates extracted from supernatants of infected cells was provided for MDS. The investigators preparing the MDS libraries and performing the bioinformatics analysis were blinded to the presence or absence of any known microbial contaminants and also to the strain and isolate identity of each virus. The six viral isolates were LACV strain H78, EBOV, rCDVRI, rMVKS, rHRSVB05, and the VSVIN strain. Samples were processed for MDS analysis as previously described (9). Samples were randomly amplified to double-stranded cDNA using the NuGEN Ovation v.2 kit (NuGEN, San Carlos, CA), and MDS libraries were constructed using the Nextera protocol (Illumina, San Diego, CA). Each library was dual indexed in the same manner as described above. Samples were pooled into a single library before library size and concentration were determined using the Blue Pippin (Sage Science, Beverly, MA) and Kapa universal quantitative PCR (qPCR) kit (Kapa Biosystems, Woburn, MA), respectively. Samples were sequenced on an Illumina HiSeq 2500 instrument using 150-nt paired-end sequencing. The paired-end sequences were analyzed for microbes using a rapid computational pipeline for microbial detection.

Bioinformatic pipeline.Paired-end reads were quality filtered using PriceSeqFilter (version 1.2, parameters –rqf 95 0.98), a component of the paired-read iterative contig extension (PRICE) assembler (29), followed by removal of human sequences by alignment to a combined reference genome including human genome build 38 (hg38) and chimpanzee (Pan troglodytes) using the Spliced Transcripts Alignment to a Reference (STAR) aligner (30). Unaligned reads that were at least 95% identical were compressed by cd-hit-dup (v4.6.1) (31, 32). After a second Bowtie2 alignment to remove residual human sequences (using an hg38 reference database), these reads were then used as queries to search the NCBI nt database (July 2015) using gsnapl (Genentech, v2015-09-29) (33).

SNV analysis.Illumina data from all the samples were processed in the same way for variant analysis. First reads were quality filtered using PriceSeqFilter (version 1.2) set to remove any read with less than 95% of nucleotides having a 0.98 probability of being correct (-rqf 95 0.98), any read with less than 90% called nucleotides, and any read that matched to the Illumina adapter sequences (29). Filtered high-quality reads were then aligned to reference genomes using GSNAP (version 2015-09-29) using default settings (33). Reference genome accession numbers were as follows: EBOV, JQ352763 ; LACV long (L) segment, NC_004108 ; LACV medium (M) segment, NC_004109 ; LACV short (S) segment, NC_004110 ; MV, HM439386 ; rHRSV, KF640637 ; VSVIN, J02428 . The rCDVRI sequence (unpublished) was supplied by W. Paul Duprex (Boston University). According to the recommended guidelines for the variant caller, PCR-duplicate reads were removed using Picard tools (version 2.2.4). Variants were called using LoFreq* (version 2.1.2) (34) using default settings with the exception of a conservative lower cutoff of ≥0.005 in frequency. LoFreq* is a sequencing read quality aware variant caller that models each location in the genome as a Poisson-binomial distribution of the number of nonreference bases. It also tests for common sequencing errors that may result in false positives, such as testing for the strand bias of a variant. It is also designed to take advantage of high-coverage (>500× coverage) data sets, such as those that we have (Table 2). It uses high-coverage data sets to call true-positive low-frequency variants while maintaining a low false-positive rate but still shows reasonable sensitivity at low coverage levels (50× coverage) (34). Variants were determined to be synonymous or nonsynonymous using a custom Python script.

Accession number(s).The 150-nt paired-end sequences are located in the National Center for Biotechnology Information (NCBI) Short Read Archive (SRA accession number SRP076690 ).

ACKNOWLEDGMENTS

We are grateful to H. Feldmann and members of the Laboratory of Virology for training of J.O. in the high-containment laboratory at the Rocky Mountain Laboratories, National Institute of Allergy and Infectious Diseases (NIAID), NIH, Hamilton, MT. We thank Derek Bogdanoff and Eric Chow in the UCSF Center for Advanced Technology for their expertise and assistance operating the Illumina sequencer.

FOOTNOTES

    • Received May 18, 2016.
    • Accepted June 16, 2016.
  • Copyright © 2016 Wilson et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .

REFERENCES

  1. 1.↵
    1. Rasko DA,
    2. Worsham PL,
    3. Abshire TG,
    4. Stanley ST,
    5. Bannan JD,
    6. Wilson MR,
    7. Langham RJ,
    8. Decker RS,
    9. Jiang L,
    10. Read TD,
    11. Phillippy AM,
    12. Salzberg SL,
    13. Pop M,
    14. Van Ert MN,
    15. Kenefic LJ,
    16. Keim PS,
    17. Fraser-Liggett CM,
    18. Ravel J
    . 2011. Bacillus anthracis comparative genome analysis in support of the Amerithrax investigation. Proc Natl Acad Sci U S A108:5027–5032. doi:10.1073/pnas.1016657108.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Hoffmaster AR,
    2. Fitzgerald CC,
    3. Ribot E,
    4. Mayer LW,
    5. Popovic T
    . 2002. Molecular subtyping of Bacillus anthracis and the 2001 bioterrorism-associated anthrax outbreak, United States. Emerg Infect Dis8:1111–1116. doi:10.3201/eid0810.020394.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Cohen J
    . 2014. Lab safety. Alarm over biosafety blunders. Science345:247–248. doi:10.1126/science.345.6194.247.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Lemon K,
    2. Nguyen DT,
    3. Ludlow M,
    4. Rennick LJ,
    5. Yüksel S,
    6. van Amerongen G,
    7. McQuaid S,
    8. Rima BK,
    9. de Swart RL,
    10. Duprex WP
    . 2015. Recombinant subgroup B human respiratory syncytial virus expressing enhanced green fluorescent protein efficiently replicates in primary human cells and is virulent in cotton rats. J Virol89:2849–2856. doi:10.1128/JVI.03587-14.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Ludlow M,
    2. Nguyen DT,
    3. Silin D,
    4. Lyubomska O,
    5. de Vries RD,
    6. von Messling V,
    7. McQuaid S,
    8. De Swart RL,
    9. Duprex WP
    . 2012. Recombinant canine distemper virus strain Snyder Hill expressing green or red fluorescent proteins causes meningoencephalitis in the ferret. J Virol86:7508–7519. doi:10.1128/JVI.06725-11.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Lemon K,
    2. de Vries RD,
    3. Mesman AW,
    4. McQuaid S,
    5. van Amerongen G,
    6. Yüksel S,
    7. Ludlow M,
    8. Rennick LJ,
    9. Kuiken T,
    10. Rima BK,
    11. Geijtenbeek TB,
    12. Osterhaus AD,
    13. Duprex WP,
    14. de Swart RL
    . 2011. Early target cells of measles virus after aerosol infection of non-human primates. PLoS Pathog7:e1001263. doi:10.1371/journal.ppat.1001263.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Kircher M,
    2. Sawyer S,
    3. Meyer M
    . 2012. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res40:e3. doi:10.1093/nar/gkr771.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Stenglein MD,
    2. Jacobson ER,
    3. Chang LW,
    4. Sanders C,
    5. Hawkins MG,
    6. Guzman DS,
    7. Drazenovich T,
    8. Dunker F,
    9. Kamaka EK,
    10. Fisher D,
    11. Reavill DR,
    12. Meola LF,
    13. Levens G,
    14. DeRisi JL
    . 2015. Widespread recombination, reassortment, and transmission of unbalanced compound viral genotypes in natural arenavirus infections. PLoS Pathog11:e1004900. doi:10.1371/journal.ppat.1004900.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Wilson MR,
    2. Shanbhag NM,
    3. Reid MJ,
    4. Singhal NS,
    5. Gelfand JM,
    6. Sample HA,
    7. Benkli B,
    8. O’Donovan BD,
    9. Ali IK,
    10. Keating MK,
    11. Dunnebacke TH,
    12. Wood MD,
    13. Bollen A,
    14. DeRisi JL
    . 2015. Diagnosing Balamuthia mandrillaris encephalitis with metagenomic deep sequencing. Ann Neurol78:722–730. doi:10.1002/ana.24499.
    OpenUrlCrossRef
  10. 10.↵
    1. Van Binnendijk RS,
    2. Poelen MC,
    3. de Vries P,
    4. Voorma HO,
    5. Osterhaus AD,
    6. Uytdehaag FG
    . 1989. Measles virus-specific human T cell clones. Characterization of specificity and function of CD4+ helper/cytotoxic and CD8+ cytotoxic T cell clones. J Immunol142:2847–2854.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Ogura H,
    2. Yoshinouchi M,
    3. Kudo T,
    4. Imura M,
    5. Fujiwara T,
    6. Yabe Y
    . 1993. Human papillomavirus type 18 DNA in so-called HEP-2, KB and FL cells—further evidence that these cells are HeLa cell derivatives. Cell Mol Biol (Noisy-le-grand)39:463–467.
    OpenUrlPubMedWeb of Science
  12. 12.↵
    1. Islam S,
    2. Zeisel A,
    3. Joost S,
    4. La Manno G,
    5. Zajac P,
    6. Kasper M,
    7. Lönnerberg P,
    8. Linnarsson S
    . 2014. Quantitative single-cell RNA-seq with unique molecular identifiers. Nat Methods11:163–166. doi:10.1038/nmeth.2772.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Switzer WM,
    2. Zheng H,
    3. Simmons G,
    4. Zhou Y,
    5. Tang S,
    6. Shankar A,
    7. Kapusinszky B,
    8. Delwart EL,
    9. Heneine W
    . 2011. No evidence of murine leukemia virus-related viruses in live attenuated human vaccines. PLoS One6:e29223. doi:10.1371/journal.pone.0029223.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Van der Zeijst BA,
    2. Bleumink N,
    3. Crawford LV,
    4. Swyryd EA,
    5. Stark GR
    . 1983. Viral proteins and RNAs in BHK cells persistently infected by lymphocytic choriomeningitis virus. J Virol48:262–270.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Kugelman JR,
    2. Rossi CA,
    3. Wiley MR,
    4. Ladner JT,
    5. Nagle ER,
    6. Pfeffer BP,
    7. Garcia K,
    8. Prieto K,
    9. Wada J,
    10. Kuhn JH,
    11. Palacios G
    . 2016. Informing the historical record of experimental nonhuman primate infections with Ebola virus: genomic characterization of USAMRIID Ebola virus/H. sapiens-tc/COD/1995/Kikwit-9510621 challenge stock “R4368” and its replacement “R4415.”PLoS One11:e0150919. doi:10.1371/journal.pone.0150919.
    OpenUrlCrossRef
  16. 16.↵
    1. Kuhn JH,
    2. Lofts LL,
    3. Kugelman JR,
    4. Smither SJ,
    5. Lever MS,
    6. van der Groen G,
    7. Johnson KM,
    8. Radoshitzky SR,
    9. Bavari S,
    10. Jahrling PB,
    11. Towner JS,
    12. Nichol ST,
    13. Palacios G
    . 2014. Reidentification of Ebola virus E718 and ME as Ebola virus/H.sapiens-tc/COD/1976/Yambuku-Ecran. Genome Announc2:e01178-14. doi:10.1128/genomeA.01178-14.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Ladner JT,
    2. Wiley MR,
    3. Mate S,
    4. Dudas G,
    5. Prieto K,
    6. Lovett S,
    7. Nagle ER,
    8. Beitzel B,
    9. Gilbert ML,
    10. Fakoli L,
    11. Diclaro JW II,
    12. Schoepp RJ,
    13. Fair J,
    14. Kuhn JH,
    15. Hensley LE,
    16. Park DJ,
    17. Sabeti PC,
    18. Rambaut A,
    19. Sanchez-Lockhart M,
    20. Bolay FK,
    21. Kugelman JR,
    22. Palacios G
    . 2015. Evolution and spread of Ebola virus in Liberia, 2014-2015. Cell Host Microbe18:659–669. doi:10.1016/j.chom.2015.11.008.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Hoenen T,
    2. Groseth A,
    3. Rosenke K,
    4. Fischer RJ,
    5. Hoenen A,
    6. Judson SD,
    7. Martellaro C,
    8. Falzarano D,
    9. Marzi A,
    10. Squires RB,
    11. Wollenberg KR,
    12. de Wit E,
    13. Prescott J,
    14. Safronetz D,
    15. van Doremalen N,
    16. Bushmaker T,
    17. Feldmann F,
    18. McNally K,
    19. Bolay FK,
    20. Fields B,
    21. Sealy T,
    22. Rayfield M,
    23. Nichol ST,
    24. Zoon KC,
    25. Massaquoi M,
    26. Munster VJ,
    27. Feldmann H
    . 2016. Nanopore sequencing as a rapidly deployable Ebola outbreak Tool. Emerg Infect Dis22:331–334. doi:10.3201/eid2202.151796.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Moon SJ,
    2. Kim TH,
    3. Yoon SY,
    4. Chung JH,
    5. Hwang HJ
    . 2015. Relationship between stage of chronic kidney disease and sarcopenia in Korean aged 40 years and older using the Korea National health and nutrition examination surveys (KNHANES IV-2, 3, and V-1, 2), 2008–2011. PLoS One10:e0130740. doi:10.1371/journal.pone.0130740.
    OpenUrlCrossRef
  20. 20.↵
    1. Baize S,
    2. Pannetier D,
    3. Oestereich L,
    4. Rieger T,
    5. Koivogui L,
    6. Magassouba N,
    7. Soropogui B,
    8. Sow MS,
    9. Keïta S,
    10. De Clerck H,
    11. Tiffany A,
    12. Dominguez G,
    13. Loua M,
    14. Traoré A,
    15. Kolié M,
    16. Malano ER,
    17. Heleze E,
    18. Bocquin A,
    19. Mély S,
    20. Raoul H,
    21. Caro V,
    22. Cadar D,
    23. Gabriel M,
    24. Pahlmann M,
    25. Tappe D,
    26. Schmidt-Chanasit J,
    27. Impouma B,
    28. Diallo AK,
    29. Formenty P,
    30. Van Herp M,
    31. Gunther S
    . 2014. Emergence of Zaire Ebola virus disease in Guinea. N Engl J Med371:1418–1425. doi:10.1056/NEJMoa1404505.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Gire SK,
    2. Goba A,
    3. Andersen KG,
    4. Sealfon RS,
    5. Park DJ,
    6. Kanneh L,
    7. Jalloh S,
    8. Momoh M,
    9. Fullah M,
    10. Dudas G,
    11. Wohl S,
    12. Moses LM,
    13. Yozwiak NL,
    14. Winnicki S,
    15. Matranga CB,
    16. Malboeuf CM,
    17. Qu J,
    18. Gladden AD,
    19. Schaffner SF,
    20. Yang X,
    21. Jiang PP,
    22. Nekoui M,
    23. Colubri A,
    24. Coomber MR,
    25. Fonnie M,
    26. Moigboi A,
    27. Gbakie M,
    28. Kamara FK,
    29. Tucker V,
    30. Konuwa E,
    31. Saffa S,
    32. Sellu J,
    33. Jalloh AA,
    34. Kovoma A,
    35. Koninga J,
    36. Mustapha I,
    37. Kargbo K,
    38. Foday M,
    39. Yillah M,
    40. Kanneh F,
    41. Robert W,
    42. Massally JL,
    43. Chapman SB,
    44. Bochicchio J,
    45. Murphy C,
    46. Nusbaum C,
    47. Young S,
    48. Birren BW,
    49. Grant DS,
    50. Scheiffelin JS,
    51. Lander ES,
    52. Happi C,
    53. Gevao SM,
    54. Gnirke A,
    55. Rambaut A,
    56. Garry RF,
    57. Khan SH,
    58. Sabeti PC
    . 2014. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science345:1369–1372. doi:10.1126/science.1259657.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Park DJ,
    2. Dudas G,
    3. Wohl S,
    4. Goba A,
    5. Whitmer SL,
    6. Andersen KG,
    7. Sealfon RS,
    8. Ladner JT,
    9. Kugelman JR,
    10. Matranga CB,
    11. Winnicki SM,
    12. Qu J,
    13. Gire SK,
    14. Gladden-Young A,
    15. Jalloh S,
    16. Nosamiefan D,
    17. Yozwiak NL,
    18. Moses LM,
    19. Jiang PP,
    20. Lin AE,
    21. Schaffner SF,
    22. Bird B,
    23. Towner J,
    24. Mamoh M,
    25. Gbakie M,
    26. Kanneh L,
    27. Kargbo D,
    28. Massally JL,
    29. Kamara FK,
    30. Konuwa E,
    31. Sellu J,
    32. Jalloh AA,
    33. Mustapha I,
    34. Foday M,
    35. Yillah M,
    36. Erickson BR,
    37. Sealy T,
    38. Blau D,
    39. Paddock C,
    40. Brault A,
    41. Amman B,
    42. Basile J,
    43. Bearden S,
    44. Belser J,
    45. Bergeron E,
    46. Campbell S,
    47. Chakrabarti A,
    48. Dodd K,
    49. Flint M,
    50. Gibbons A,
    51. Goodman C,
    52. Klena J,
    53. McMullan L,
    54. Morgan L,
    55. Russell B,
    56. Salzer J,
    57. Sanchez A,
    58. Wang D,
    59. Jungreis I,
    60. Tomkins-Tinch C,
    61. Kislyuk A,
    62. Lin MF,
    63. Chapman S,
    64. MacInnis B,
    65. Matthews A,
    66. Bochicchio J,
    67. Hensley LE,
    68. Kuhn JH,
    69. Nusbaum C,
    70. Schieffelin JS,
    71. Birren BW,
    72. Forget M,
    73. Nichol ST,
    74. Palacios GF,
    75. Ndiaye D,
    76. Happi C,
    77. Gevao SM,
    78. Vandi MA,
    79. Kargbo B,
    80. Holmes EC,
    81. Bedford T,
    82. Gnirke A,
    83. Ströher U,
    84. Rambaut A,
    85. Garry RF,
    86. Sabeti PC
    . 2015. Ebola virus epidemiology, transmission, and evolution during seven months in Sierra Leone. Cell161:1516–1526. doi:10.1016/j.cell.2015.06.007.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Kugelman JR,
    2. Wiley MR,
    3. Mate S,
    4. Ladner JT,
    5. Beitzel B,
    6. Fakoli L,
    7. Taweh F,
    8. Prieto K,
    9. Diclaro JW,
    10. Minogue T,
    11. Schoepp RJ,
    12. Schaecher KE,
    13. Pettitt J,
    14. Bateman S,
    15. Fair J,
    16. Kuhn JH,
    17. Hensley L,
    18. Park DJ,
    19. Sabeti PC,
    20. Sanchez-Lockhart M,
    21. Bolay FK,
    22. Palacios G,
    23. US Army Medical Research Institute of Infectious Diseases, National Institutes of Health, Integrated Research Facility—Frederick Ebola Response Team 2014–2015
    . 2015. Monitoring of Ebola virus Makona evolution through establishment of advanced genomic capability in Liberia. Emerg Infect Dis21:1135–1143. doi:10.3201/eid2107.150522.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Simon-Loriere E,
    2. Faye O,
    3. Faye O,
    4. Koivogui L,
    5. Magassouba N,
    6. Keita S,
    7. Thiberge JM,
    8. Diancourt L,
    9. Bouchier C,
    10. Vandenbogaert M,
    11. Caro V,
    12. Fall G,
    13. Buchmann JP,
    14. Matranga CB,
    15. Sabeti PC,
    16. Manuguerra JC,
    17. Holmes EC,
    18. Sall AA
    . 2015. Distinct lineages of Ebola virus in Guinea during the 2014 West African epidemic. Nature524:102–104. doi:10.1038/nature14612.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Tong YG,
    2. Shi WF,
    3. Liu D,
    4. Qian J,
    5. Liang L,
    6. Bo XC,
    7. Liu J,
    8. Ren HG,
    9. Fan H,
    10. Ni M,
    11. Sun Y,
    12. Jin Y,
    13. Teng Y,
    14. Li Z,
    15. Kargbo D,
    16. Dafae F,
    17. Kanu A,
    18. Chen CC,
    19. Lan ZH,
    20. Jiang H,
    21. Luo Y,
    22. Lu HJ,
    23. Zhang XG,
    24. Yang F,
    25. Hu Y,
    26. Cao YX,
    27. Deng YQ,
    28. Su HX,
    29. Sun Y,
    30. Liu WS,
    31. Wang Z,
    32. Wang CY,
    33. Bu ZY,
    34. Guo ZD,
    35. Zhang LB,
    36. Nie WM,
    37. Bai CQ,
    38. Sun CH,
    39. An XP,
    40. Xu PS,
    41. Zhang XL,
    42. Huang Y,
    43. Mi ZQ,
    44. Yu D,
    45. Yao HW,
    46. Feng Y,
    47. Xia ZP,
    48. Zheng XX,
    49. Yang ST,
    50. Lu B,
    51. Jiang JF,
    52. Kargbo B,
    53. He FC,
    54. Gao GF,
    55. Cao WC,
    56. China Mobile Laboratory Testing Team in Sierra Leone
    . 2015. Genetic diversity and evolutionary dynamics of Ebola virus in Sierra Leone. Nature524:93–96. doi:10.1038/nature14490.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Hoenen T,
    2. Safronetz D,
    3. Groseth A,
    4. Wollenberg KR,
    5. Koita OA,
    6. Diarra B,
    7. Fall IS,
    8. Haidara FC,
    9. Diallo F,
    10. Sanogo M,
    11. Sarro YS,
    12. Kone A,
    13. Togo AC,
    14. Traore A,
    15. Kodio M,
    16. Dosseh A,
    17. Rosenke K,
    18. de Wit E,
    19. Feldmann F,
    20. Ebihara H,
    21. Munster VJ,
    22. Zoon KC,
    23. Feldmann H,
    24. Sow S
    . 2015. Virology. Mutation rate and genotype variation of Ebola virus from Mali case sequences. Science348:117–119. doi:10.1126/science.aaa5646.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Chandler LJ,
    2. Borucki MK,
    3. Dobie DK,
    4. Wasieloski LP,
    5. Thompson WH,
    6. Gundersen CB,
    7. Case K,
    8. Beaty BJ
    . 1998. Characterization of La Crosse virus RNA in autopsied central nervous system tissues. J Clin Microbiol36:3332–3336.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with bowtie 2. Nat Methods9:357–359. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Ruby JG,
    2. Bellare P,
    3. Derisi JL
    . 2013. PRICE: software for the targeted assembly of components of (Meta) genomic sequence data. G3 (Bethesda)3:865–880. doi:10.1534/g3.113.005967.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Dobin A,
    2. Davis CA,
    3. Schlesinger F,
    4. Drenkow J,
    5. Zaleski C,
    6. Jha S,
    7. Batut P,
    8. Chaisson M,
    9. Gingeras TR
    . 2013. STAR: ultrafast universal RNA-seq aligner. Bioinformatics29:15–21. doi:10.1093/bioinformatics/bts635.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Fu L,
    2. Niu B,
    3. Zhu Z,
    4. Wu S,
    5. Li W
    . 2012. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics28:3150–3152. doi:10.1093/bioinformatics/bts565.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Li W,
    2. Godzik A
    . 2006. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics22:1658–1659. doi:10.1093/bioinformatics/btl158.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Wu TD,
    2. Nacu S
    . 2010. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics26:873–881. doi:10.1093/bioinformatics/btq057.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Wilm A,
    2. Aw PP,
    3. Bertrand D,
    4. Yeo GH,
    5. Ong SH,
    6. Wong CH,
    7. Khor CC,
    8. Petric R,
    9. Hibberd ML,
    10. Nagarajan N
    . 2012. LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res40:11189–11201. doi:10.1093/nar/gks918.
    OpenUrlCrossRefPubMedWeb of Science
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Multiplexed Metagenomic Deep Sequencing To Analyze the Composition of High-Priority Pathogen Reagents
Michael R. Wilson, Greg Fedewa, Mark D. Stenglein, Judith Olejnik, Linda J. Rennick, Sham Nambulli, Friederike Feldmann, W. Paul Duprex, John H. Connor, Elke Mühlberger, Joseph L. DeRisi
mSystems Jul 2016, 1 (4) e00058-16; DOI: 10.1128/mSystems.00058-16

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Multiplexed Metagenomic Deep Sequencing To Analyze the Composition of High-Priority Pathogen Reagents
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Multiplexed Metagenomic Deep Sequencing To Analyze the Composition of High-Priority Pathogen Reagents
Michael R. Wilson, Greg Fedewa, Mark D. Stenglein, Judith Olejnik, Linda J. Rennick, Sham Nambulli, Friederike Feldmann, W. Paul Duprex, John H. Connor, Elke Mühlberger, Joseph L. DeRisi
mSystems Jul 2016, 1 (4) e00058-16; DOI: 10.1128/mSystems.00058-16
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

NO MESH TERMS ASSIGNED AS OF 10-11-2019
metagenomics
pathogen tracking
phylogenetic analysis

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077