To explore the feasibility of pyrosequencing for quantitative differential gene expression analysis we have performed a comparative study of the results of the sequencing experiments to those obtained by a conventional DNA microarray platform. A conclusion from our analysis is that, over a threshold of 35 normalized reads per gene, the measurements of gene expression display a good correlation with the references. The observed concordance between pyrosequencing and DNA microarray platforms beyond the threshold was of 0.8, measured as a Pearson's correlation coefficient. In differential gene expression the initial aim is the quantification the differences among transcripts when comparing experimental conditions. Thus, even in a scenario of low coverage the concordance in the measurements is quite acceptable. On the other hand, the comparatively longer read size obtained by pyrosequencing allows detecting unconventional splicing forms.
Get full access to this article
View all access options for this article.
References
1.
AlberolaT.M., Garcia-MartinezJ., AntunezO., ViladevallL., BarceloA., ArinoJ.et al.2004. A new set of DNA macrochips for the yeast Saccharomyces cerevisiae: features and uses. Int Microbiol, 7:199–206.
2.
ArinoJ.2010. Integrative responses to high pH stress in S. cerevisiae. Omics, 14:517–523.
HedbackerK., CarlsonM.2008. SNF1/AMPK pathways in yeast. Front Biosci, 13:2408–2420.
9.
HegdeP., QiR., AbernathyK., GayC., DharapS., GaspardR.et al.2000. A concise guide to cDNA microarray analysis. Biotechniques, 29:548–550552–544, 556 passim.
10.
IoannidisJ.P., AllisonD.B., BallC.A., CoulibalyI., CuiX., CulhaneA.C.et al.2009. Repeatability of published microarray gene expression analyses. Nat Genet, 41:149–155.
11.
JuneauK., PalmC., MirandaM., DavisR.W.2007. High-density yeast-tiling array reveals previously undiscovered introns and extensive regulation of meiotic splicing. Proc Natl Acad Sci USA, 104:1522–1527.
12.
LambT.M., MitchellA.P.2003. The transcription factor Rim101p governs ion tolerance and cell differentiation by direct repression of the regulatory genes NRG1 and SMP1 in Saccharomyces cerevisiae. Mol Cell Biol, 23:677–686.
13.
LambT.M., XuW., DiamondA., MitchellA.P.2001. Alkaline response genes of Saccharomyces cerevisiae and their relationship to the RIM101 pathway. J Biol Chem, 276:1850–1856.
14.
LiH., DurbinR.2010. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics, 26:589–595.
15.
MarioniJ.C., MasonC.E., ManeS.M., StephensM., GiladY.2008. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res, 18:1509–1517.
16.
MedinaI., CarbonellJ., PulidoL., MadeiraS.C., GoetzS., ConesaA.et al.2010. Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res, 38,Suppl:W210–W213.
17.
MiuraF., KawaguchiN., SeseJ., ToyodaA., HattoriM., MorishitaS.et al.2006. A large-scale full-length cDNA analysis to explore the budding yeast transcriptome. Proc Natl Acad Sci USA, 103:17846–17851.
18.
MortazaviA., WilliamsB.A., McCueK., SchaefferL., WoldB.2008. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods, 5:621–628.
19.
NagalakshmiU., WangZ., WaernK., ShouC., RahaD., GersteinM.et al.2008. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science, 320:1344–1349.
20.
OshlackA., RobinsonM.D., YoungM.D.2010. From RNA-seq reads to differential expression results. Genome Biol, 11:220.
21.
PlataraM., RuizA., SerranoR., PalominoA., MorenoF., ArinoJ.2006. The transcriptional response of the yeast Na(+)-ATPase ENA1 gene to alkaline stress involves three main signaling pathways. J Biol Chem, 281:36632–36642.
22.
RuizA., SerranoR., ArinoJ.2008. Direct regulation of genes involved in glucose utilization by the calcium/calcineurin pathway. J Biol Chem, 283:13923–13933.
23.
SerranoR., MartinH., CasamayorA., ArinoJ.2006. Signaling alkaline pH stress in the yeast Saccharomyces cerevisiae through the Wsc1 cell surface sensor and the Slt2 MAPK pathway. J Biol Chem, 281:39785–39795.
24.
SerranoR., RuizA., BernalD., ChambersJ.R., ArinoJ.2002. The transcriptional response to alkaline pH in Saccharomyces cerevisiae: evidence for calcium-mediated signalling. Mol Microbiol, 46:1319–1333.
25.
ShiL., ReidL.H., JonesW.D., ShippyR., WarringtonJ.A., BakerS.C.et al.2006. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol, 24:1151–1161.
26.
SmythG.2005. Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. GentlemanR., CareyV., DudoitS., IrizarryR., HuberW.Springer: New York, 397–420.
27.
SultanM., SchulzM.H., RichardH., MagenA., KlingenhoffA., ScherfM.et al.2008. A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science, 321:956–960.
28.
T'HoenP.A.C., AriyurekY., ThygesenH.H., VreugdenhilE., VossenR.H., de MenezesR.X.et al.2008. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res, 36:e141.
29.
TrapnellC., WilliamsB.A., PerteaG., MortazaviA., KwanG., Van BarenM.J.et al.2010. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol, 28:511–515.
30.
ViladevallL., SerranoR., RuizA., DomenechG., GiraldoJ., BarceloA.et al.2004. Characterization of the calcium-mediated response to alkaline stress in Saccharomyces cerevisiae. J Biol Chem, 279:43614–43624.
31.
WangZ., GersteinM., SnyderM.2009. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet, 10:57–63.
