Comparative analysis of RNA sequencing methods for degraded or low-input samples (original) (raw)

Nature Methods volume 10, pages 623–629 (2013)Cite this article

Subjects

A Corrigendum to this article was published on 30 January 2014

This article has been updated

Abstract

RNA-seq is an effective method for studying the transcriptome, but it can be difficult to apply to scarce or degraded RNA from fixed clinical samples, rare cell populations or cadavers. Recent studies have proposed several methods for RNA-seq of low-quality and/or low-quantity samples, but the relative merits of these methods have not been systematically analyzed. Here we compare five such methods using metrics relevant to transcriptome annotation, transcript discovery and gene expression. Using a single human RNA sample, we constructed and sequenced ten libraries with these methods and compared them against two control libraries. We found that the RNase H method performed best for chemically fragmented, low-quality RNA, and we confirmed this through analysis of actual degraded samples. RNase H can even effectively replace oligo(dT)-based methods for standard RNA-seq. SMART and NuGEN had distinct strengths for measuring low-quantity RNA. Our analysis allows biologists to select the most suitable methods and provides a benchmark for future method development.

This is a preview of subscription content, access via your institution

Access options

Subscribe to this journal

Receive 12 print issues and online access

$259.00 per year

only $21.58 per issue

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Additional access options:

Similar content being viewed by others

Accession codes

Primary accessions

Gene Expression Omnibus

Referenced accessions

NCBI Reference Sequence

Change history

In the version of this article initially published, in the Online Methods "RNase H libraries" section, the sentence beginning with "We added 5 μl preheated RNase H...." should have read "We added 5 μl preheated RNase H reaction mix that contains 10 U of Hybridase Thermostable RNase H (Epicentre), 0.5 μmol Tris-HCl, pH 7.5, 1 μmol NaCl and 0.2 μmol MgCl2 to the RNA and DNA oligo mix, incubated this mixture at 45 °C for 30 min and then placed it on ice." The errors have been corrected in the HTML and PDF versions of this article.

References

  1. Aviv, H. & Leder, P. Purification of biologically active globin messenger RNA by chromatography on oligothymidylic acid-cellulose. Proc. Natl. Acad. Sci. USA 69, 1408–1412 (1972).
    Article CAS PubMed PubMed Central Google Scholar
  2. Yang, L., Duff, M.O., Graveley, B.R., Carmichael, G.G. & Chen, L.L. Genomewide characterization of non-polyadenylated RNAs. Genome Biol. 12, R16 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  3. Tang, F. et al. mRNA-Seq whole-transcriptome analysis of a single cell. Nat. Methods 6, 377–382 (2009).
    Article CAS PubMed Google Scholar
  4. Ramsköld, D. et al. Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells. Nat. Biotechnol. 30, 777–782 (2012).
    Article PubMed PubMed Central Google Scholar
  5. Islam, S. et al. Characterization of the single-cell transcriptional landscape by highly multiplex RNA-seq. Genome Res. 21, 1160–1167 (2011).
    CAS PubMed PubMed Central Google Scholar
  6. Sinicropi, D. & Morlan, J. Methods for depleting RNA from nucleic acid samples. US patent application 20110111409 (2011).
  7. Morlan, J.D., Qu, K. & Sinicropi, D.V. Selective depletion of rRNA enables whole transcriptome profiling of archival fixed tissue. PLoS ONE 7, e42882 (2012).
    Article CAS PubMed PubMed Central Google Scholar
  8. Huang, R. et al. An RNA-Seq strategy to detect the complete coding and non-coding transcriptome including full-length imprinted macro ncRNAs. PLoS ONE 6, e27288 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  9. Yi, H. et al. Duplex-specific nuclease efficiently removes rRNA for prokaryotic RNA-seq. Nucleic Acids Res. 39, e140 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  10. Tariq, M.A., Kim, H.J., Jejelowo, O. & Pourmand, N. Whole-transcriptome RNAseq analysis from minute amount of total RNA. Nucleic Acids Res. 39, e120 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  11. Levin, J.Z. et al. Comprehensive comparative analysis of strand-specific RNA sequencing methods. Nat. Methods 7, 709–715 (2010).
    Article CAS PubMed PubMed Central Google Scholar
  12. DeLuca, D.S. et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 28, 1530–1532 (2012).
    Article CAS PubMed PubMed Central Google Scholar
  13. Beyer, A.L. & Osheim, Y.N. Splice site selection, rate of splicing, and alternative splicing on nascent transcripts. Genes Dev. 2, 754–765 (1988).
    Article CAS PubMed Google Scholar
  14. Yang, Y.H. et al. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 30, e15 (2002).
    Article PubMed PubMed Central Google Scholar
  15. Aird, D. et al. Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries. Genome Biol. 12, R18 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  16. Rosenkranz, R., Borodina, T., Lehrach, H. & Himmelbauer, H. Characterizing the mouse ES cell transcriptome with Illumina sequencing. Genomics 92, 187–194 (2008).
    Article CAS PubMed Google Scholar
  17. Giannoukos, G. et al. Efficient and robust RNA-seq process for cultured bacteria and complex community transcriptomes. Genome Biol. 13, R23 (2012).
    Article CAS PubMed PubMed Central Google Scholar
  18. Griffin, M., Abu-El-Haija, M., Abu-El-Haija, M., Rokhlina, T. & Uc, A. Simplified and versatile method for isolation of high-quality RNA from pancreas. Biotechniques 52, 332–334 (2012).
    Article CAS PubMed PubMed Central Google Scholar
  19. Pan, X. et al. Two methods for full-length RNA sequencing for low quantities of cells and single cells. Proc. Natl. Acad. Sci. USA 110, 594–599 (2013).
    Article CAS PubMed Google Scholar
  20. Roberts, A., Trapnell, C., Donaghey, J., Rinn, J.L. & Pachter, L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 12, R22 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  21. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).
    Article CAS PubMed PubMed Central Google Scholar
  22. Maden, B.E. et al. Clones of human ribosomal DNA containing the complete 18 S-rRNA and 28 S-rRNA genes. Characterization, a detailed map of the human ribosomal transcription unit and diversity among clones. Biochem. J. 246, 519–527 (1987).
    Article CAS PubMed PubMed Central Google Scholar
  23. Trapnell, C., Pachter, L. & Salzberg, S.L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–1111 (2009).
    Article CAS PubMed PubMed Central Google Scholar
  24. Dreszer, T.R. et al. The UCSC Genome Browser database: extensions and updates 2011. Nucleic Acids Res. 40, D918–D923 (2012).
    Article CAS PubMed Google Scholar
  25. Li, B. & Dewey, C.N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323 (2011).
    Article CAS PubMed PubMed Central Google Scholar
  26. Langmead, B., Trapnell, C., Pop, M. & Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).
    Article PubMed PubMed Central Google Scholar
  27. Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Use R!) (Springer, New York, 2009).

Download references

Acknowledgements

We thank D. Sinicropi for alerting us to his RNase H method, Broad Genomics Platform for sequencing work, C. Nusbaum, R. Johnson, M. Yassour and V. Savova for helpful discussions, and L. Gaffney for assistance in drafting figures. Work was supported by US National Institutes of Health Pioneer Award DP1-OD003958-01, US National Human Genome Research Institute (NHGRI) 1P01HG005062-01, NHGRI Center of Excellence in Genome Science Award 1P50HG006193-01, the Howard Hughes Medical Institute, the Merkin Foundation for Stem Cell Research and the Klarman Cell Observatory at the Broad Institute (to A.R.), and by NHGRI grant HG03067. N.P. was supported by a postdoctoral research fellowship of the Fund for Scientific Research–Flanders (FWO Vlaanderen).

Author information

Author notes

  1. Xian Adiconis and Diego Borges-Rivera: These authors contributed equally to this work.

Authors and Affiliations

  1. Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA
    Xian Adiconis, Diego Borges-Rivera, Rahul Satija, David S DeLuca, Michele A Busby, Aaron M Berlin, Andrey Sivachenko, Dawn Anne Thompson, Alec Wysoker, Timothy Fennell, Andreas Gnirke, Nathalie Pochet, Aviv Regev & Joshua Z Levin
  2. Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
    Aviv Regev
  3. Howard Hughes Medical Institute, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
    Aviv Regev

Authors

  1. Xian Adiconis
    You can also search for this author inPubMed Google Scholar
  2. Diego Borges-Rivera
    You can also search for this author inPubMed Google Scholar
  3. Rahul Satija
    You can also search for this author inPubMed Google Scholar
  4. David S DeLuca
    You can also search for this author inPubMed Google Scholar
  5. Michele A Busby
    You can also search for this author inPubMed Google Scholar
  6. Aaron M Berlin
    You can also search for this author inPubMed Google Scholar
  7. Andrey Sivachenko
    You can also search for this author inPubMed Google Scholar
  8. Dawn Anne Thompson
    You can also search for this author inPubMed Google Scholar
  9. Alec Wysoker
    You can also search for this author inPubMed Google Scholar
  10. Timothy Fennell
    You can also search for this author inPubMed Google Scholar
  11. Andreas Gnirke
    You can also search for this author inPubMed Google Scholar
  12. Nathalie Pochet
    You can also search for this author inPubMed Google Scholar
  13. Aviv Regev
    You can also search for this author inPubMed Google Scholar
  14. Joshua Z Levin
    You can also search for this author inPubMed Google Scholar

Contributions

J.Z.L., X.A. and A.R. conceived the research. X.A. prepared the cDNA libraries. D.B.-R., R.S., N.P., M.A.B. and A.R. developed and performed computational analysis. D.S.D. contributed code. D.S.D., A.M.B., A.S., A.W. and T.F. helped with computational analysis. D.A.T., N.P., A.R. and J.Z.L. supervised the research. J.Z.L., X.A., D.B.-R. and A.R. wrote the paper. R.S., A.G. and D.S.D. assisted in editing the paper.

Corresponding authors

Correspondence toAviv Regev or Joshua Z Levin.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Supplementary information

Rights and permissions

About this article

Cite this article

Adiconis, X., Borges-Rivera, D., Satija, R. et al. Comparative analysis of RNA sequencing methods for degraded or low-input samples.Nat Methods 10, 623–629 (2013). https://doi.org/10.1038/nmeth.2483

Download citation

This article is cited by