Cited by 512Gene expression and cancer classificationGenomics and Phylogenetic StudiesBioinformatics and Genomic Networks

The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance

Open access
Wang, Charles, Gong, Binsheng, Bushel, Pierre R, Thierry-Mieg, Jean, Thierry-Mieg, Danielle, Xu, Joshua, Fang, Hong, Hong, Huixiao +38 more·2014·Nature Biotechnology

Abstract

The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R(2)0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/nbt.3001

Cites

7 of 7 references verified0 citation errors · 0 not checkable

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • Among 23.2 billion reads mapped to genes other than mitochondrial or ribosomal RNAs, 85.9% mapped to RefSeq, 92.9% to GENCODE, and 97.1% to NCBI AceView.

    “Among all 23.2 billion reads that could be mapped to genes other than those encoding mitochondrial or ribosomal RNAs, 85.9% were mapped to RefSeq, 12 whereas 92.9% mapped to GENCODE, 13 and 97.1% to NCBI AceView 14 (Fig. 2a).”

    From A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control Consortium (2014) · cites this paper as [12]

    The Results (Study design; page unavailable) say the study used an AceView-based RNA-seq pipeline and refer to mapping details, but do not give the claimed read-mapping percentages. Searches for 23.2 billion reads and for the RefSeq, GENCODE, and AceView percentages found no passage in our copy reporting those figures. (LLM, conf=0.73)

    Experimental — accuracy not yet measured

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Cites

Flagged first, then by citations

Research works printed but not linked: 36

  1. [1]Hamburg MA Advancing regulatory science Science 2011 331 987 21350129
  2. [2]Chen M Zhang M Borlak J Tong W A decade of toxicogenomic research and its contribution to toxicological science Toxicological sciences: an official journal of the Society of Toxicology 2012 130 217 228 22790972
  3. [3]Shi L The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements Nature biotechnology 2006 24 1151 1161 PMC3272078 16964229
  4. [4]Guo L Rat toxicogenomic study reveals analytical consistency across microarray platforms Nature biotechnology 2006 24 1162 1169 17061323
  5. [6]Fan X Consistency of predictive signature genes and classifiers generated using different microarray platforms The pharmacogenomics journal 2010 10 247 257 20676064 PMC2920073
  6. [8]Marioni JC Mason CE Mane SM Stephens M Gilad Y RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays Genome Res 2008 18 1509 1517 18550803 PMC2527709
  7. [9]Bottomly D Evaluating gene expression in C57BL/6J and DBA/2J mouse striatum using RNA-Seq and microarrays PloS one 2011 6 e17820 21455293 PMC3063777
  8. [10]Bradford JR A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling BMC Genomics 2010 11 282 20444259 PMC2877694
  9. [12]Malone JH Oliver B Microarrays, deep sequencing and the true measure of the transcriptome BMC biology 2011 9 34 21627854 PMC3104486
  10. [13]Merrick BA RNA-Seq Profiling Reveals Novel Hepatic Gene Expression Pattern in Aflatoxin B1 Treated Rats PloS one 2013 8 e61768 23630614 PMC3632591
  11. [15]Raghavachari N A systematic comparison and evaluation of high density exon arrays and RNA-seq technology used to unravel the peripheral blood transcriptome of sickle cell disease BMC Med Genomics 2012 5 28 22747986 PMC3428653
  12. [16]Sirbu A Kerr G Crane M Ruskin HJ RNA-Seq vs dual- and single-channel microarray data: sensitivity analysis for differential expression and clustering PloS one 2012 7 e50986 23251411 PMC3518479
  13. [17]Su Z Comparing next-generation sequencing and microarray technologies in a toxicological study of the effects of aristolochic acid on rat kidneys Chem Res Toxicol 2011 24 1486 1493 21834575
  14. [18]Subramaniam S Hsiao G Gene-expression measurement: variance-modeling considerations for robust data analysis Nature immunology 2012 13 199 203 22344273 PMC4358796
  15. [19]Xiong Y RNA sequencing shows no dosage compensation of the active X-chromosome Nat Genet 2010 42 1043 1047 21102464
  16. [20]Xu W Human transcriptome array for high-throughput clinical studies Proc Natl Acad Sci U S A 2011 108 3707 3712 21317363 PMC3048146
  17. [21]Labaj PP Characterization and improvement of RNA-Seq precision in quantitative transcript expression profiling Bioinformatics 2011 27 Oxford, England 383 391 21685096 PMC3117338
  18. [22]McIntyre LM RNA-seq: technical variability and sampling BMC Genomics 2011 12 293 21645359 PMC3141664
  19. [23]Mooney M Comparative RNA-Seq and microarray analysis of gene expression changes in B-cell lymphomas of Canis familiaris PloS one 2013 8 e61088 23593398 PMC3617154
  20. [24]Sultan M A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome Science 2008 321 956 960 18599741
  21. [27]Irizarry RA Summaries of Affymetrix GeneChip probe level data Nucleic acids research 2003 31 e15 12582260 PMC150247
  22. [28]Li C Wong WH Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection Proc Natl Acad Sci U S A 2001 98 31 36 11134512 PMC14539
  23. [30]Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data Bioinformatics 2010 26 Oxford, England 139 140 19910308 PMC2796818
  24. [31]Anders S Huber W Differential expression analysis for sequence count data Genome biology 2010 11 R106 20979621 PMC3218662
  25. [32]Shi L Microarray scanner calibration curves: characteristics and implications BMC bioinformatics 2005 6 Suppl 2 S11 16026596 PMC1637029
  26. [33]Kupershmidt I Ontology-based meta-analysis of global collections of high-throughput public data PloS one 2010 5 e13066 20927376 PMC2947508
  27. [34]Lu J Bushel PR Dynamic expression of 3' UTRs revealed by Poisson hidden Markov modeling of RNA-Seq: implications in gene expression profiling Gene 2013 527 616 623 23845781 PMC3902974
  28. [35]Katz Y Wang ET Airoldi EM Burge CB Analysis and design of RNA sequencing experiments for identifying isoform regulation Nature methods 2010 7 1009 1015 21057496 PMC3037023
  29. [37]Loven J Revisiting global gene expression analysis Cell 2012 151 476 482 23101621 PMC3505597
  30. [38]Ganter B Development of a large-scale chemogenomics database to improve drug candidate selection and to understand mechanisms of chemical toxicity and action Journal of Biotechnology 2005 119 219 244 16005536
  31. [39]Liu WM Analysis of high density expression microarrays with signed-rank call algorithms Bioinformatics 2002 18 Oxford, England 1593 1599 12490443
  32. [41]Wu Z Irizarry RA Gentleman R Martinez-Murillo F Spencer F A model-based background adjustment for oligonucleotide expression arrays Journal of the American Statistical Association 2004 99 909 917
  33. [43]Wingender E The TRANSFAC system on gene expression regulation Nucleic acids research 2001 29 281 283 11125113 PMC29801
  34. [44]Wingender E TRANSFAC: an integrated system for gene expression regulation Nucleic acids research 2000 28 316 319 10592259 PMC102445
  35. [45]Breslin T Krogh M Peterson C Troein C Signal transduction pathway profiling of individual tumor samples BMC bioinformatics 2005 6 163 15987529 PMC1184060
  36. [46]Trapnell C Pachter L Salzberg SL TopHat: discovering splice junctions with RNA-Seq Bioinformatics 2009 25 Oxford, England 1105 1111 19289445 PMC2672628

Cited by

Sorted by citations

Authors
Wang, C., Gong, B., Bushel, P. R., Thierry-Mieg, J., Thierry-Mieg, D., Xu, J., Fang, H., Hong, H., Shen, J., Su, Z., Meehan, J., Li, X., Yang, L., Li, H., Łabaj, P. P., Kreil, D. P., Megherbi, D., Gaj, S., Caiment, F., van Delft, J., Kleinjans, J., Scherer, A., Devanarayan, V., Wang, J., Yang, Y., Qian, H.-R., Lancashire, L. J., Bessarabova, M., Nikolsky, Y., Furlanello, C., Chierici, M., Albanese, D., Jurman, G., Riccadonna, S., Filosi, M., Visintainer, R., Zhang, K. K., Li, J., Hsieh, J.-H., Svoboda, D. L., Fuscoe, J. C., Deng, Y., Shi, L., Paules, R. S., Auerbach, S. S., & Tong, W.
Affiliations
  • Loma Linda University (US)
  • National Center for Toxicological Research (US)
  • National Institutes of Health (US)
  • City Of Hope National Medical Center (US)
  • BOKU University (AT)
  • University of Warwick (GB)
  • University of Massachusetts Lowell (US)
  • Maastricht University (NL)
  • Walter and Eliza Hall Institute of Medical Research (AU)
  • AbbVie (United States) (US)
  • Eli Lilly (United States) (US)
  • Thomson Reuters (United States) (US)
  • Russian Academy of Sciences (RU)
  • Fondazione Bruno Kessler (IT)
  • Fondazione Edmund Mach (IT)
  • University of North Dakota (US)
  • National Institute of Environmental Health Sciences (US)
  • SRA International (United States) (US)
  • Rush University Medical Center (US)
Corresponding author
  • Bushel, Pierre R
  • Auerbach, Scott S
  • Tong, Weida
Venue
Nature Biotechnology
Published
2014-08-24
Volume
32
Issue
9
Pages
926-932
ISSN
  • 0733-222X (ISSN-L)
  • 1087-0156
  • 1546-1696
OpenAlex
W1989996002
Cited by
512
References
7
Citations by year
2026: 9 · 2025: 24 · 2024: 61 · 2023: 37 · 2022: 26 · 2021: 61 · 2020: 38 · 2019: 39 · 2018: 39 · 2017: 46 · 2016: 63 · 2015: 56 · 2014: 12
FWCI
24.21
Type
journal article
Language
en
Open access
Locations
  • Nature Biotechnology (publishedVersion)
  • PubMed (publishedVersion)
  • Research Publications (Maastricht University) (publishedVersion)
  • Europe PMC (PubMed Central) (acceptedVersion)
  • PubMed Central (acceptedVersion)
  • CINECA IRIS Institutional Research Information System (Fondazione Edmund Mach) (submittedVersion, other-oa)
  • Research Publications (Maastricht University) (submittedVersion)
Funders
  • Vienna Science and Technology Fund
  • National Institutes of Health
  • National Institute of Environmental Health Sciences (ZIAES102345, ES102345-04, Z01ES023026, Z01ES102345)
  • U.S. National Library of Medicine
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
Keywords
cross-platform comparison, chemical exposure, differential gene expression, mRNA levels, low-abundance transcripts, mechanism of action, hepatotoxicity, quantitative PCR validation, high-throughput compound screening
MeSH terms
Animals, RNA, Messenger, Sequence Analysis, RNA, Oligonucleotide Array Sequence Analysis, Rats

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review