title: A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling.
abstract: BACKGROUND
  RNA-Seq exploits the rapid generation of gigabases of sequence data by Massively Parallel Nucleotide Sequencing, allowing for the mapping and digital quantification of whole transcriptomes. Whilst previous comparisons between RNA-Seq and microarrays have been performed at the level of gene expression, in this study we adopt a more fine-grained approach. Using RNA samples from a normal human breast epithelial cell line (MCF-10a) and a breast cancer cell line (MCF-7), we present a comprehensive comparison between RNA-Seq data generated on the Applied Biosystems SOLiD platform and data from Affymetrix Exon 1.0ST arrays. The use of Exon arrays makes it possible to assess the performance of RNA-Seq in two key areas: detection of expression at the granularity of individual exons, and discovery of transcription outside annotated loci.
  
  
  RESULTS
  We found a high degree of correspondence between the two platforms in terms of exon-level fold changes and detection. For example, over 80% of exons detected as expressed in RNA-Seq were also detected on the Exon array, and 91% of exons flagged as changing from Absent to Present on at least one platform had fold-changes in the same direction. The greatest detection correspondence was seen when the read count threshold at which to flag exons Absent in the SOLiD data was set to t<1 suggesting that the background error rate is extremely low in RNA-Seq. We also found RNA-Seq more sensitive to detecting differentially expressed exons than the Exon array, reflecting the wider dynamic range achievable on the SOLiD platform. In addition, we find significant evidence of novel protein coding regions outside known exons, 93% of which map to Exon array probesets, and are able to infer the presence of thousands of novel transcripts through the detection of previously unreported exon-exon junctions.
  
  
  CONCLUSIONS
  By focusing on exon-level expression, we present the most fine-grained comparison between RNA-Seq and microarrays to date. Overall, our study demonstrates that data from a SOLiD RNA-Seq experiment are sufficient to generate results comparable to those produced from Affymetrix Exon arrays, even using only a single replicate from each platform, and when presented with a large genome.
authors.affiliation: Applied Computational Biology and Bioinformatics, Cancer Research UK, Paterson Institute for Cancer Research, The University of Manchester, Wilmslow Road, Manchester M20 4BX UK.
authors.forename: James R
authors.initials: JR
authors.lastname: Bradford
authors.forename: Yvonne
authors.initials: Y
authors.lastname: Hey
authors.forename: Tim
authors.initials: T
authors.lastname: Yates
authors.forename: Yaoyong
authors.initials: Y
authors.lastname: Li
authors.forename: Stuart D
authors.initials: SD
authors.lastname: Pepper
authors.forename: Crispin J
authors.initials: CJ
authors.lastname: Miller
mesh_terms: D017398:Alternative Splicing; D045744:Cell Line, Tumor; D041322:Chromosomes, Human, Y; D005091:Exons; D015870:Gene Expression; D006801:Humans; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D014158:Transcription, Genetic
references.citation: Nucleic Acids Res. 2008 Jan;36(Database issue):D780-6
references.pmid: 17932061
references.citation: Bioinformatics. 2005 May 1;21(9):2067-75
references.pmid: 15657102
references.citation: Methods. 2009 Jul;48(3):249-57
references.pmid: 19336255
references.citation: Nature. 2008 Jun 26;453(7199):1239-43
references.pmid: 18488015
references.citation: Bioinformatics. 2005 Sep 15;21(18):3683-5
references.pmid: 16076888
references.citation: Genome Biol. 2007;8(5):R82
references.pmid: 17504534
references.citation: Cancer Res. 2007 Jun 15;67(12):5635-42
references.pmid: 17575129
references.citation: Genome Biol. 2008;9(12):R175
references.pmid: 19087247
references.citation: Nat Methods. 2008 Jul;5(7):613-9
references.pmid: 18516046
references.citation: Nucleic Acids Res. 2009 Jan;37(Database issue):D690-7
references.pmid: 19033362
references.citation: Nat Biotechnol. 2009 May;27(5):455-7
references.pmid: 19430453
references.citation: Nat Methods. 2008 Jul;5(7):621-8
references.pmid: 18516045
references.citation: Stat Appl Genet Mol Biol. 2004;3:Article3
references.pmid: 16646809
references.citation: BMC Genomics. 2006 Dec 27;7:325
references.pmid: 17192196
references.citation: Biotechniques. 2007 Feb;42(2):181-5
references.pmid: 17373482
references.citation: Nat Methods. 2009 May;6(5):377-82
references.pmid: 19349980
references.citation: BMC Bioinformatics. 2006 Jun 02;7:276
references.pmid: 16749918
references.citation: Nat Genet. 2008 Dec;40(12):1413-5
references.pmid: 18978789
references.citation: Science. 2008 Aug 15;321(5891):956-60
references.pmid: 18599741
references.citation: Biol Direct. 2009 Apr 16;4:14
references.pmid: 19371405
references.citation: Nature. 2007 Jun 14;447(7146):799-816
references.pmid: 17571346
references.citation: Biochim Biophys Acta. 1975 Oct 20;405(2):442-51
references.pmid: 1180967
references.citation: Nucleic Acids Res. 2003 Feb 15;31(4):e15
references.pmid: 12582260
references.citation: Genome Res. 2008 Sep;18(9):1509-17
references.pmid: 18550803
references.citation: Science. 1995 Oct 20;270(5235):484-7
references.pmid: 7570003
references.citation: Science. 2008 Jun 6;320(5881):1344-9
references.pmid: 18451266
references.citation: Nat Rev Genet. 2009 Jan;10(1):57-63
references.pmid: 19015660
references.citation: Genome Biol. 2004;5(10):R80
references.pmid: 15461798
references.citation: Bioinformatics. 2010 Jan 1;26(1):139-40
references.pmid: 19910308
references.citation: Genome Res. 1997 Oct;7(10):986-95
references.pmid: 9331369
references.citation: Nucleic Acids Res. 2008 Jan;36(Database issue):D281-8
references.pmid: 18039703