A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling.
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.
Bradford
James R
JR
Applied Computational Biology and Bioinformatics, Cancer Research UK, Paterson Institute for Cancer Research, The University of Manchester, Wilmslow Road, Manchester M20 4BX UK.
Hey
Yvonne
Y
Yates
Tim
T
Li
Yaoyong
Y
Pepper
Stuart D
SD
Miller
Crispin J
CJ
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
Nucleic Acids Res. 2008 Jan;36(Database issue):D780-6
17932061
Bioinformatics. 2005 May 1;21(9):2067-75
15657102
Methods. 2009 Jul;48(3):249-57
19336255
Nature. 2008 Jun 26;453(7199):1239-43
18488015
Bioinformatics. 2005 Sep 15;21(18):3683-5
16076888
Genome Biol. 2007;8(5):R82
17504534
Cancer Res. 2007 Jun 15;67(12):5635-42
17575129
Genome Biol. 2008;9(12):R175
19087247
Nat Methods. 2008 Jul;5(7):613-9
18516046
Nucleic Acids Res. 2009 Jan;37(Database issue):D690-7
19033362
Nat Biotechnol. 2009 May;27(5):455-7
19430453
Nat Methods. 2008 Jul;5(7):621-8
18516045
Stat Appl Genet Mol Biol. 2004;3:Article3
16646809
BMC Genomics. 2006 Dec 27;7:325
17192196
Biotechniques. 2007 Feb;42(2):181-5
17373482
Nat Methods. 2009 May;6(5):377-82
19349980
BMC Bioinformatics. 2006 Jun 02;7:276
16749918
Nat Genet. 2008 Dec;40(12):1413-5
18978789
Science. 2008 Aug 15;321(5891):956-60
18599741
Biol Direct. 2009 Apr 16;4:14
19371405
Nature. 2007 Jun 14;447(7146):799-816
17571346
Biochim Biophys Acta. 1975 Oct 20;405(2):442-51
1180967
Nucleic Acids Res. 2003 Feb 15;31(4):e15
12582260
Genome Res. 2008 Sep;18(9):1509-17
18550803
Science. 1995 Oct 20;270(5235):484-7
7570003
Science. 2008 Jun 6;320(5881):1344-9
18451266
Nat Rev Genet. 2009 Jan;10(1):57-63
19015660
Genome Biol. 2004;5(10):R80
15461798
Bioinformatics. 2010 Jan 1;26(1):139-40
19910308
Genome Res. 1997 Oct;7(10):986-95
9331369
Nucleic Acids Res. 2008 Jan;36(Database issue):D281-8
18039703