Location of Repository

Error, reproducibility and sensitivity : a pipeline for data processing of Agilent oligonucleotide expression arrays

By B. M. Chain, Helen C. Bowen, John P. Hammond, Wilfried Posch, Jane Rasaiyaah, Jhen Tsang and Mahdad Noursadeghi


Background\ud Expression microarrays are increasingly used to obtain large scale transcriptomic information on a wide range of biological samples. Nevertheless, there is still much debate on the best ways to process data, to design experiments and analyse the output. Furthermore, many of the more sophisticated mathematical approaches to data analysis in the literature remain inaccessible to much of the biological research community. In this study we examine ways of extracting and analysing a large data set obtained using the Agilent long oligonucleotide transcriptomics platform, applied to a set of human macrophage and dendritic cell samples.\ud \ud Results\ud We describe and validate a series of data extraction, transformation and normalisation steps which are implemented via a new R function. Analysis of replicate normalised reference data demonstrate that intrarray variability is small (only around 2% of the mean log signal), while interarray variability from replicate array measurements has a standard deviation (SD) of around 0.5 log2 units ( 6% of mean). The common practise of working with ratios of Cy5/Cy3 signal offers little further improvement in terms of reducing error. Comparison to expression data obtained using Arabidopsis samples demonstrates that the large number of genes in each sample showing a low level of transcription reflect the real complexity of the cellular transcriptome. Multidimensional scaling is used to show that the processed data identifies an underlying structure which reflect some of the key biological variables which define the data set. This structure is robust, allowing reliable comparison of samples collected over a number of years and collected by a variety of operators.\ud \ud Conclusions\ud This study outlines a robust and easily implemented pipeline for extracting, transforming normalising and visualising transcriptomic array data from Agilent expression platform. The analysis is used to obtain quantitative estimates of the SD arising from experimental (non biological) intra- and interarray variability, and for a lower threshold for determining whether an individual gene is expressed. The study provides a reliable basis for further more extensive studies of the systems biology of eukaryotic cells.\u

Topics: QP
Publisher: BioMed Central Ltd.
Year: 2010
OAI identifier: oai:wrap.warwick.ac.uk:3339

Suggested articles



  1. (2010). Agilent Whole Human Genome Expression Arrays
  2. (2008). den Dunnen JT: Deep sequencingbased expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res doi
  3. DM: A variance-stabilizing transformation for gene-expression microarray data. Bioinformatics 2002, 18(Suppl 1):S105-S110. doi
  4. (2009). DR: Genome-wide innate immune responses in HIV-1-infected macrophages are preserved despite attenuation of the NF-kappa B activation pathway. doi
  5. (2007). Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature doi
  6. (2007). JD: A new approach to intensity-dependent normalization of two-channel microarrays. Biostatistics doi
  7. Multidimensional scaling. In Applied multivariate data analysis 2nd edition. Hodder doi
  8. (2010). RA: Frozen robust multiarray analysis (fRMA). Biostatistics doi
  9. (2007). Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics doi
  10. (2006). RD: Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project. Nat Biotechnol doi
  11. (2007). Selection and validation of normalization methods for cDNA microarrays using within-array replications. Bioinformatics doi
  12. ST, Broadley MR: A comparison of the Thlaspi caerulescens and Thlaspi arvense shoot transcriptomes. New Phytol doi
  13. (2007). Storey JD: Normalization of two-channel microarrays accounting for experimental design and intensity-dependent relationships. Genome Biol doi
  14. (2007). Systematic approaches for incorporating control spots and data quality information to improve normalization of cDNA microarray data. doi
  15. (2006). The MicroArray Quality Control (MAQC) project shows interand intraplatform reproducibility of gene expression measurements. Nat Biotechnol doi
  16. (2003). TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics doi
  17. (2009). Transcriptional and functional defects of dendritic cells derived from the MUTZ-3 leukaemia line. Immunology doi
  18. (2004). White PJ: Cesium toxicity in Arabidopsis. Plant Physiol
  19. (2002). Wolfl S: Ranking: a closer look on globalisation methods for normalisation of gene expression arrays. Nucleic Acids Res doi
  20. (2008). Yaspo ML: A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science doi

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.