CoChat
Sign Up
Cited by 168Genomics and Phylogenetic StudiesEnvironmental DNA in Biodiversity StudiesMicrobial Community Ecology and Physiology

The DOI matches its Crossref record; 23 of the 24 references it cites were found in a registry or an index.

Unlocking Short Read Sequencing for Metagenomics

Open access
Rodrigue, Sébastien, Materna, Arne C., Timberlake, Sonia C., Blackburn, Matthew C., Malmstrom, Rex R., Alm, Eric J., Chisholm, Sallie W.·2010·PLoS ONE
DOI 10.1371/journal.pone.0011840PMC PMC2911387
Open this paper in CoChat — free accountOpen PDF

Abstract

BACKGROUND: Different high-throughput nucleic acid sequencing platforms are currently available but a trade-off currently exists between the cost and number of reads that can be generated versus the read length that can be achieved. METHODOLOGY/PRINCIPAL FINDINGS: We describe an experimental and computational pipeline yielding millions of reads that can exceed 200 bp with quality scores approaching that of traditional Sanger sequencing. The method combines an automatable gel-less library construction step with paired-end sequencing on a short-read instrument. With appropriately sized library inserts, mate-pair sequences can overlap, and we describe the SHERA software package that joins them to form a longer composite read. CONCLUSIONS/SIGNIFICANCE: This strategy is broadly applicable to sequencing applications that benefit from low-cost high-throughput sequencing, but require longer read lengths. We demonstrate that our approach enables metagenomic analyses using the Illumina Genome Analyzer, with low error rates, and at a fraction of the cost of pyrosequencing.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1371/journal.pone.0011840

Cites

23 of 23 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Size-dependent isolation of DNA fragments from sheared genomic DNA via dSPRI. AMPure XP SPRI beads bind DNA fragments in a size dependent manner according to the concentration of salts and polyethylene glycol (PEG) in the reaction [4]–[7], which can easily be changed by using different volume ratios of DNA to SPRI bead solutions. A two-step procedure is employed to isolate targeted DNA size fractions. Panels A to H present Bioanalyzer DNA-1000 assays showing the sheared genomic DNA used as starting material (black), the larger size DNA fragments discarded in separation 1 (red), and the size fraction purified and recovered after separation 2 (blue). Panel I is a table summarizing the conditions and results displayed in panels A to H. All Bioanalyzer DNA-1000 traces after separation 1 (panel J), and after separation 2 (panel K), are respectively displayed on a graph for the conditions presented in panels A to H. The conditions displayed in panel H were used to obtaine the Illumina composite reads discussed in the text. The wider DNA fragment size distribution from panel H allowed to better analyze the effects of shorter versus longer overlapping regions on consensus reads.
Figure 1 Size-dependent isolation of DNA fragments from sheared genomic DNA via dSPRI. AMPure XP SPRI beads bind DNA fragments in a size dependent manner according to the concentration of salts and polyethylene glycol (PEG) in the reaction [4]–[7], which can easily be changed by using different volume ratios of DNA to SPRI bead solutions. A two-step procedure is employed to isolate targeted DNA size fractions. Panels A to H present Bioanalyzer DNA-1000 assays showing the sheared genomic DNA used as starting material (black), the larger size DNA fragments discarded in separation 1 (red), and the size fraction purified and recovered after separation 2 (blue). Panel I is a table summarizing the conditions and results displayed in panels A to H. All Bioanalyzer DNA-1000 traces after separation 1 (panel J), and after separation 2 (panel K), are respectively displayed on a graph for the conditions presented in panels A to H. The conditions displayed in panel H were used to obtaine the Illumina composite reads discussed in the text. The wider DNA fragment size distribution from panel H allowed to better analyze the effects of shorter versus longer overlapping regions on consensus reads.CC BY · PubMed Central
Figure 2 Reproducibility of double-SPRI. The panel shows DNA fragment size distributions as obtained by Bioanalyzer DNA-1000 assays. The curves represent the size fractions removed during the first separation step A) or recovered after the second separation B). The two size fractions were independently reproduced in 4 or 8 separation experiments. The curves in B) represent the libraries sequenced after dSPRI based size selection, adapter ligation and PCR enrichment. While concentrations (arbitrary fluorescence units) vary between reproduced libraries the range of removed or enriched DNA fragment sizes was highly reproducible. Panel c) shows the DNA fragment size distribution recovered after the second separation when using decreasing amounts sheared genomic DNA. dSPRI allows reliable size selection in a DNA concentration independent manner.
Figure 2 Reproducibility of double-SPRI. The panel shows DNA fragment size distributions as obtained by Bioanalyzer DNA-1000 assays. The curves represent the size fractions removed during the first separation step A) or recovered after the second separation B). The two size fractions were independently reproduced in 4 or 8 separation experiments. The curves in B) represent the libraries sequenced after dSPRI based size selection, adapter ligation and PCR enrichment. While concentrations (arbitrary fluorescence units) vary between reproduced libraries the range of removed or enriched DNA fragment sizes was highly reproducible. Panel c) shows the DNA fragment size distribution recovered after the second separation when using decreasing amounts sheared genomic DNA. dSPRI allows reliable size selection in a DNA concentration independent manner.CC BY · PubMed Central
Figure 3 Quality of Illumina reads out to 143 bp. A) The mean Phred quality of single reads is shown in the solid red lines, with error bars displaying quartiles. Read quality is highly variable toward the reads' ends. Read quality as a function of base pair is worse on the second mate of the pair. B) The mean and quartiles of Phred quality by base for the average-length composite read.
Figure 3 Quality of Illumina reads out to 143 bp. A) The mean Phred quality of single reads is shown in the solid red lines, with error bars displaying quartiles. Read quality is highly variable toward the reads' ends. Read quality as a function of base pair is worse on the second mate of the pair. B) The mean and quartiles of Phred quality by base for the average-length composite read.CC BY · PubMed Central
Figure 4 High-confidence alignment yield by insert lengths. Distribution of insert lengths by aligning the original reads to the reference sequence (gray), and lengths of composite reads retained (red).
Figure 4 High-confidence alignment yield by insert lengths. Distribution of insert lengths by aligning the original reads to the reference sequence (gray), and lengths of composite reads retained (red).CC BY · PubMed Central
Figure 5 Quality of composite fragments after overlapping. The mean Phred quality (and corresponding error rate) at each position of the read. We show the original Illumina reads and the composite fragments generated by overlapping. A) 143 bp, the length of the original Illumina read; B) 180 bp, the mean length of overlapped fragment in this library and C) 250 bp, roughly the mean length generated by 454-FLX technology.
Figure 5 Quality of composite fragments after overlapping. The mean Phred quality (and corresponding error rate) at each position of the read. We show the original Illumina reads and the composite fragments generated by overlapping. A) 143 bp, the length of the original Illumina read; B) 180 bp, the mean length of overlapped fragment in this library and C) 250 bp, roughly the mean length generated by 454-FLX technology.CC BY · PubMed Central
Figure 6 Composite Illumina reads constitute a legitimate alternative to pyrosequencing for metagenomics studies. DNA from a marine metagenomics sample was sequenced with our overlapping mate-pair approach. The composite reads were directly compared to 454-FLX sequences from the exact same sample. A) Fraction of reads that could be assigned to a taxon using the composite reads (mean read length 180 bp), the entire 454-FLX dataset (mean read length 207 bp), or longer 454-FLX reads (mean read length 254 bp). B) Comparison of taxon assignments using composite reads and 454-FLX pyrosequencing reads. The top 25 represented taxa, with colored symbols next to the bracket, are listed in Table S1.
Figure 6 Composite Illumina reads constitute a legitimate alternative to pyrosequencing for metagenomics studies. DNA from a marine metagenomics sample was sequenced with our overlapping mate-pair approach. The composite reads were directly compared to 454-FLX sequences from the exact same sample. A) Fraction of reads that could be assigned to a taxon using the composite reads (mean read length 180 bp), the entire 454-FLX dataset (mean read length 207 bp), or longer 454-FLX reads (mean read length 254 bp). B) Comparison of taxon assignments using composite reads and 454-FLX pyrosequencing reads. The top 25 represented taxa, with colored symbols next to the bracket, are listed in Table S1.CC BY · PubMed Central
Figure 7 Short insertions and deletions in low-confidence composite Illumina reads. Mate-reads from the control lane (PhiX174 bacteriophage genome) were used to assess false alignments introduced by the SHERA pipeline. After constructing the composite sequences and filtering, we plotted the difference in length (if any) between a composite fragment and its insert length as predicted by MAQ by mapping the original mate-pairs to the reference.This histogram is a sum of two distributions, the overlapper software's misalignments (a broad gaussian) and a sharp peak of small (1–2 bp) indels. We used a simple and conservative linear model to remove the indels and infer the number of misalignments.
Figure 7 Short insertions and deletions in low-confidence composite Illumina reads. Mate-reads from the control lane (PhiX174 bacteriophage genome) were used to assess false alignments introduced by the SHERA pipeline. After constructing the composite sequences and filtering, we plotted the difference in length (if any) between a composite fragment and its insert length as predicted by MAQ by mapping the original mate-pairs to the reference.This histogram is a sum of two distributions, the overlapper software's misalignments (a broad gaussian) and a sharp peak of small (1–2 bp) indels. We used a simple and conservative linear model to remove the indels and infer the number of misalignments.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 24 references

Cited by

Sorted by citations

See all 168

Related

Metadata

Authors
Rodrigue, S., Materna, A. C., Timberlake, S. C., Blackburn, M. C., Malmstrom, R. R., Alm, E. J., & Chisholm, S. W.
Affiliations
Massachusetts Institute of Technology (US)
Corresponding author
  • Alm, Eric J.
  • Chisholm, Sallie W.
ORCID
Venue
PLoS ONE
Publisher
Public Library of Science (PLoS)
Published
2010-07-28
Volume
5
Issue
7
Pages
e11840
ISSN
1932-6203 (ISSN-L)
DOI
10.1371/journal.pone.0011840
PMC
PMC2911387
OpenAlex
W2151671981
Cited by
168
References
24
Citations by year
2025: 5 · 2024: 2 · 2023: 6 · 2022: 2 · 2021: 1 · 2020: 6 · 2019: 3 · 2018: 10 · 2017: 7 · 2016: 13 · 2015: 21 · 2014: 32 · 2013: 29 · 2012: 16
FWCI
7.56
Type
journal article
Language
en
Open access
Licence
Locations
  • PLoS ONE (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • OSTI OAI (U.S. Department of Energy Office of Scientific and Technical Information) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • DSpace@MIT (Massachusetts Institute of Technology) (submittedVersion, cc-by-nc)
  • Figshare (submittedVersion, cc-by)
Funders
  • U.S. Department of Energy (DEFG0207ER64506)
  • Gordon and Betty Moore Foundation
  • Fonds Québécois de la Recherche sur la Nature et les Technologies
  • Natural Sciences and Engineering Research Council of Canada
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Physical Sciences › Environmental Science › Ecology › Environmental DNA in Biodiversity Studies
  • Physical Sciences › Environmental Science › Ecology › Microbial Community Ecology and Physiology
Keywords
metagenomics, short-read sequencing, long-read sequencing, paired-end sequencing, sequencing cost, read length
MeSH terms
Software, Sequence Analysis, DNA, Metagenomics

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review