title: Transcriptome sequencing of the Microarray Quality Control (MAQC) RNA reference samples using next generation sequencing.
abstract: BACKGROUND
  Transcriptome sequencing using next-generation sequencing platforms will soon be competing with DNA microarray technologies for global gene expression analysis. As a preliminary evaluation of these promising technologies, we performed deep sequencing of cDNA synthesized from the Microarray Quality Control (MAQC) reference RNA samples using Roche's 454 Genome Sequencer FLX.
  
  
  RESULTS
  We generated more that 3.6 million sequence reads of average length 250 bp for the MAQC A and B samples and introduced a data analysis pipeline for translating cDNA read counts into gene expression levels. Using BLAST, 90% of the reads mapped to the human genome and 64% of the reads mapped to the RefSeq database of well annotated genes with e-values </= 10-20. We measured gene expression levels in the A and B samples by counting the numbers of reads that mapped to individual RefSeq genes in multiple sequencing runs to evaluate the MAQC quality metrics for reproducibility, sensitivity, specificity, and accuracy and compared the results with DNA microarrays and Quantitative RT-PCR (QRTPCR) from the MAQC studies. In addition, 88% of the reads were successfully aligned directly to the human genome using the AceView alignment programs with an average 90% sequence similarity to identify 137,899 unique exon junctions, including 22,193 new exon junctions not yet contained in the RefSeq database.
  
  
  CONCLUSIONS
  Using the MAQC metrics for evaluating the performance of gene expression platforms, the ExpressSeq results for gene expression levels showed excellent reproducibility, sensitivity, and specificity that improved systematically with increasing shotgun sequencing depth, and quantitative accuracy that was comparable to DNA microarrays and QRTPCR. In addition, a careful mapping of the reads to the genome using the AceView alignment programs shed new light on the complexity of the human transcriptome including the discovery of thousands of new splice variants.
authors.affiliation: Department of Biological Sciences, Virginia Tech, Blacksburg, VA 24061, USA. rvjensen@vt.edu
authors.forename: Shrinivasrao P
authors.initials: SP
authors.lastname: Mane
authors.forename: Clive
authors.initials: C
authors.lastname: Evans
authors.forename: Kristal L
authors.initials: KL
authors.lastname: Cooper
authors.forename: Oswald R
authors.initials: OR
authors.lastname: Crasta
authors.forename: Otto
authors.initials: O
authors.lastname: Folkerts
authors.forename: Stephen K
authors.initials: SK
authors.lastname: Hutchison
authors.forename: Timothy T
authors.initials: TT
authors.lastname: Harkins
authors.forename: Danielle
authors.initials: D
authors.lastname: Thierry-Mieg
authors.forename: Jean
authors.initials: J
authors.lastname: Thierry-Mieg
authors.forename: Roderick V
authors.initials: RV
authors.lastname: Jensen
mesh_terms: D018076:DNA, Complementary / Q000235:genetics; D030541:Databases, Genetic; D020869:Gene Expression Profiling / Q000379:methods*; D015723:Gene Library; D015894:Genome, Human; D006801:Humans; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D011786:Quality Control; D012015:Reference Standards; D012680:Sensitivity and Specificity; D016415:Sequence Alignment; D017423:Sequence Analysis, RNA / Q000379:methods*; D012984:Software
chemical_list: D018076:DNA, Complementary