Email updates

Keep up to date with the latest news and content from BMC Research Notes and BioMed Central.

This article is part of the series Microarray normalization and optimization.

Open Access Highly Accessed Technical Note

EMA - A R package for Easy Microarray data analysis

Nicolas Servant123*, Eleonore Gravier1234, Pierre Gestraud123, Cecile Laurent123678, Caroline Paccard123, Anne Biton1235, Isabel Brito123, Jonas Mandel123, Bernard Asselain123, Emmanuel Barillot123 and Philippe Hupé1235

Author Affiliations

1 Institut Curie, Paris F-75248, France

2 INSERM, U900, Paris F-75248, France

3 Ecole des Mines ParisTech, Fontainebleau, F-77300 France

4 Institut Curie, Departement de Transfert, Paris F-75248, France

5 CNRS, UMR144, Paris F-75248, France

6 CNRS, UMR3347, Orsay F-91405, France

7 INSERM, U1021, Orsay F-91405, France

8 Université Paris-Sud 11, Orsay F-91405, France

For all author emails, please log on.

BMC Research Notes 2010, 3:277  doi:10.1186/1756-0500-3-277

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1756-0500/3/277


Received:15 July 2010
Accepted:3 November 2010
Published:3 November 2010

© 2010 Servant et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The increasing number of methodologies and tools currently available to analyse gene expression microarray data can be confusing for non specialist users.

Findings

Based on the experience of biostatisticians of Institut Curie, we propose both a clear analysis strategy and a selection of tools to investigate microarray gene expression data. The most usual and relevant existing R functions were discussed, validated and gathered in an easy-to-use R package (EMA) devoted to gene expression microarray analysis. These functions were improved for ease of use, enhanced visualisation and better interpretation of results.

Conclusions

Strategy and tools proposed in the EMA R package could provide a useful starting point for many microarrays users. EMA is part of Comprehensive R Archive Network and is freely available at http://bioinfo.curie.fr/projects/ema/ webcite.

Findings

Numerous analysis methods and tools have been developed to study microarray, many of them being implemented as free R [1] and/or Bioconductor [2] packages. This abundance of methods makes choosing the best approach difficult for newcomers and non-specialist users.

Based on the experience of the biostatisticians of Institut Curie, we propose a clear analysis strategy combining a large variety of standard methodologies. The most usual and relevant R functions needed to perform these analyses were selected and gathered in the R package EMA (Easy Microarray data Analysis). EMA covers an entire analysis process including quality control, normalisation, exploratory analysis, unsupervised and supervised classification, functional analysis and censored data exploration. The package can be used for both one or two-colours gene expression micrarrays and for exon expression experiments.

Analysis strategy

Firstly, the quality of the data must be assessed in order to detect problematic raw probe-level data, such as spatial artifacts on the chip or poor quality hybridisation. Indeed, gene expression experiments suffer from many sources of technical and experimental variation. Removing noise and systematic biases is performed in order to both improve the biological signal and make all the arrays comparable. This is the so-called normalisation step. Secondly, we propose to discard the probesets with very low signal across the samples (i.e. genes unexpressed or below detection threshold). This filtering step leads to both a noise reduction in the data and an increase in the statistical power of the subsequent analysis.

Then, exploratory approaches are classically used to find clusters of genes (or samples) with similar profiles. Note that here, biological interpretation depends on the choice of the similarity metrics. These approaches potentially highlight outliers and/or non relevant effects (batch effect for example), which can be subsequently estimated and/or removed from the data thanks to appropriate methods.

Finally, supervised approaches aim at the identification of differentially expressed genes (DEG), or deregulated pathways by taking into account the multiple testing issues. The biological interpretation of the differential analysis results can be performed thanks to functional and gene set enrichment analyses. Sample class prediction (eg good vs poor clinical outcome) based on supervised classification methods can also be performed to highlight genes signatures.

Selected tools

For the data quality assessment, we recommend to use the arrayQualityMetrics package [3], which performs a powerful, easy-to-use and comprehensive data quality estimation as well as an automatic html report. The EMA package proposes the most famous techniques for Affymetrix GeneChip normalisation: MAS5.0 [4], RMA [5] and GCRMA [6]. We recommend to use GCRMA because it outperforms the other approaches (by ignoring the mismatch intensities and taking into account the probe sequence information) and allows an efficient filtering of irrelevant probesets thanks to its bimodal distribution of probesets expression values (Figure 1a). Other packages such as limma [7], vsn [8] or lumi [9] can be used to normalise non Affymetrix data. After this first step, the main EMA functions can be used for any type of expression data, using a simple data expression matrix as input.

thumbnailFigure 1. Graphical outputs provided by the EMA package for the class comparison study of [18]. (a) Histogram of probesets expression values across the 23 samples after GCRMA normalisation and log2 transformation. Probesets with an expression value below 3.5 (red vertical line) are discarded. (b) Individuals factor map produced by the PCA performed on the 23 filtered gene expression profiles. (c) Heatmap of the 23 gene expression profiles based on the 100 genes with the highest interquartile range (IQR) values. Sample clustering was performed using Pearson's correlation coefficient and Ward criterion. Gene clustering was performed using absolute Pearson's correlation coefficient and Ward criterion. (d) Qqplot produced by the SAM analysis on the two groups of tumours. Probesets in green are considered to be differentially expressed between the two conditions.

The EMA package provides functions to perform exploratory analyses such as Principal Component Analysis (PCA, Figure 1b), hierarchical clustering (Figure 1c) or Multiple Factor Analysis. They are based on R packages such as FactoMineR [10], cluster [11], or mostclust [12]. The use of linear model is proposed to estimate and to remove the non relevant effects potentially detected.

Various methods are proposed to perform differential analysis and their choice depends on the sample size. The multtest package provides standard approaches like Student or Mann-Whitney test associated with multiple testing correction methods. The Significance Analysis of Microarrays (SAM) approach [13] (siggenes package) is also very interesting because it both estimates the null distribution and takes into account the correlation between probesets (Figure 1d). The rank product method [14] (RankProd package) dedicated to small sample size dataset is also offered, as well as some linear model (ANOVA) functions. Alternatively, the user can apply the limma package which is a very powerful tool to assess differential expression by linear models.

The functional enrichment of the DEG list is assessed based on the GeneOntology [15], and KEGG [16] pathways annotation terms. The hyper-geometric test of the GOstats package is used to test the over-representation of the functional terms in the gene list.

For sample class prediction, we suggest to use the CMA package [17] including the most popular machine learning and gene selection algorithms. In the context of censored data, the EMA package supports Kaplan Meier and log-rank analyses using the survival package.

Example

The proposed analysis strategy was applied to the breast cancer gene expression dataset [18] comparing 12 Basal-like carcinomas (BLCs) and 11 HER2 positive carcinomas (HER2+). Some graphical outputs for data preprocessing, exploratory analysis and differential analysis steps are displayed in Figure 1. The RNA profiles were analysed using U133 plus 2.0 Affymetrix GeneChip. Three genes (P-cadherin, v-kit, FOXC1) were reported by the authors to be associated to a genes cluster over-expressed in the basal-like carcinomas and three genes (PTEN, Her2 and GRB7) to a genes cluster over-expressed in the Her2+ carcinomas. All these genes but one (v-kit) were found to be differentially expressed using the EMA package. This discrepancy is easily explained because in spite of v-kit belongs to a basal-like expression cluster, no change in v-kit expression can be observed between the two groups in this clustering analysis. This is because the hierarchical clustering was performed on genes (such as v-kit) not necessary differentially expressed between the two populations.

The R scripts used to analyse this gene expression dataset can be found in [Additional file 1]. Transcriptomic data used in this application are publicly available at Gene Expression Omnibus (Accession number: [GSE13787]) and are part of the package.

Additional file 1. R scripts applied to the breast cancer gene expression dataset [18]. R script used to analyse the breast cancer gene expression data set [18].

Format: DOC Size: 106KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Conclusions

EMA is a freely available R package which implements a complete strategy for expression microarray analysis. The package includes a vignette [Additional file 2] which describes the detailed biological/clinical analysis strategy used at Institut Curie. Most of the functions were improved for ease of use (fewer command lines, default parameters tested and chosen to be optimal). Relevant, enhanced and easy-to-interpret text and graphic outputs are offered. The package is available on The Comprehensive R Archive Network repository [19].

Additional file 2. EMA vignette. The vignette discuss the detailed biological/clinical analysis strategy used at Institut Curie and presents an application to a gene expression dataset.

Format: PDF Size: 2.8MB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Availability and requirements

• Project Name: EMA

• Project home page:

http://bioinfo.curie.fr/projects/ema/ webcite

http://cran.r-project.org/ webcite

• Operating systems: Linux, Windows

• Programming language: R

• Other requirements: R version ≥ 2.10. R packages: cluster, Hmisc, heatmap.plus, FactoMineR, GOstats, survival, multtest, affy, gcrma, rgl, GSA, RankProd, siggenes, MASS, hgu133plus2.db, xtable, biomaRt.

• License: GNU GPL

• Any restrictions to use by non-academics: none

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NS and EG discussed the choice of the strategy and tools, participated to the development of the EMA package and wrote the paper. PG, CL, CP, AB, IB, JM discussed the choice of the strategy and tools and participated to the development of the EMA package. BA, EB and PH discussed the choice of the strategy and tools and supervised the work group. All authors read and approved the final manuscript.

Acknowledgements

We thank all the persons from the ABCIS methodology group for their collaboration and fruitful discussions. This work was supported by a grant from the Institut Curie "Programmes incitatifs et coopératifs Bioinformatique et Biostatistique". EB and PH are members of the team Systems Biology of Cancer "équipes labellisées par la Ligue Nationale Contre le Cancer".

References

  1. R Development Core Team: R: A Language and Environment for Statistical Computing. [http://www.R-project.org] webcite

    R Foundation for Statistical Computing, Vienna, Austria; 2009.

    [ISBN 3-900051-07-0]

  2. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics.

    Genome Biol 2004, 5:R80. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  3. Kauffmann A, Gentleman R, Huber W: arrayQualityMetrics-a bioconductor package for quality assessment of microarray data.

    Bioinformatics 2009, 25:415-6. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Statistical Algorithms Description Document [http://www.affymetrix.com/support/technical/whitepapers.affx] webcite

  5. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data.

    Biostatistics 2003, 4:249-64. PubMed Abstract | Publisher Full Text OpenURL

  6. Wu Z, Irizarry R, Gentleman R, Martinez-Murillo F, Spencer F: A Model-Based Background Adjustment for Oligonucleotide Expression Arrays.

    Journal of the American Statistical Association 2004, 99(468):909-917. Publisher Full Text OpenURL

  7. K SG, P ST: Normalization of cDNA microarray data.

    Methods 2003, 31:265-273. PubMed Abstract | Publisher Full Text OpenURL

  8. A HWVHASHP, M V: Variance stabilization applied to microarray data calibration and to the quantification of differential expression.

    Bioinformatics 2002, 18:S96-S104. PubMed Abstract | Publisher Full Text OpenURL

  9. Du P LS Kibbe WA: Lumi:a pipeline for processing Illumina microarray.

    Bioinformatics 2008, 24:1547-1548. PubMed Abstract | Publisher Full Text OpenURL

  10. Lê S, Josse J, Husson F: FactoMineR: an R package for multivariate analysis.

    Journal of statistical software 2008, 25:1-18. OpenURL

  11. Kaufman L, Rousseeuw P: Finding groups in data - An introduction to cluster analysis. New York: John Wiley & Sons; 1990. OpenURL

  12. Bertoni A, Valentini G: Model order selection for bio-molecular data clustering.

    BMC Bioinformatics 2007, 8(Suppl 2):S7. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response.

    Proc Natl Acad Sci USA 2001, 98:5116-21. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis.

    Bioinformatics 2006, 22:2825-7. PubMed Abstract | Publisher Full Text OpenURL

  15. Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la Cruz N, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource.

    Nucleic Acids Res 2004, 32:D258-61. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes.

    Nucleic Acids Res 1999, 27:29-34. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Slawski M, Daumer M, Boulesteix AL: CMA: a comprehensive Bioconductor package for supervised classification with high dimensional data.

    BMC Bioinformatics 2008, 9:439. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Marty B, Maire V, Gravier E, Rigaill G, Vincent-Salomon A, Kappler M, Lebigot I, Djelti F, Tourdes A, Gestraud P, Hupe P, Barillot E, Cruzalegui F, Tucker GC, Stern MH, Thiery JP, Hickman JA, Dubois T: Frequent PTEN genomic alterations and activated phosphatidylinositol 3-kinase pathway in basal-like breast cancer cells.

    Breast Cancer Res 2008, 10:R101. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. The Comprehensive R Archive Network [http://cran.r-project.org/] webcite