GEDI: a user-friendly toolbox for analysis of large-scale gene expression data
© Fujita et al; licensee BioMed Central Ltd. 2007
Received: 21 May 2007
Accepted: 19 November 2007
Published: 19 November 2007
Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills.
Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al.
GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.
High-throughput DNA microarray technologies yield up to tens of thousands of gene expression data, which are useful to identify differentially expressed genes, biomarkers and molecular disease profiles. In recent years, microarray platforms have become available at relatively low costs, becoming more popular among research groups which are interested in gene expression analysis. On the other hand, much effort has been spent in developing improved methods to analyze the data derived from these microarrays. These methods involve advanced mathematical and statistical models, which are quite cumbersome to biomedical researchers who attempt to implement these methods. Due to this difficulty, some of these advanced methods are often abandoned and data analysis is carried out using only the classical methods, which are implemented in popular statistical softwares. An user-friendly software could make it possible to use recently developed methods to integrate, qualify, and infer biological insights from gene expression data.
Our aim is to provide a toolbox named GEDI containing an user-friendly interface and advanced statistical methods to analyze data derived from DNA microarrays.
SVR (Support Vector Regression)  – a microarray data normalization method published in this journal, is based on a non-parametric regression, namely, Support Vector Regression, which is more robust to outliers (differentially expressed genes), therefore, this approach is superior to classical methods, such as Loess, to identify differentially expressed genes even for very lowly or very highly expressed genes, for which the expression variance is high.
DVAR (Dynamic Vector Auto Regressive Model)  – a gene expression regulatory network inference model based on time-series data, avoiding stationarity and linearity assumptions, since it is well known that different cell cycle phases involve different circuits. Hence, using DVAR, it is possible to infer different connectivities occuring during different cell cycle phases. The DVAR approach does not require model pre-specification being, therefore, unbiased. The inferred connectivities are causalities based on the Granger causality concept. This is naturally applied to networks containing cycles (feedback mechanisms).
SVAR (Sparse Vector Auto Regressive Model)  – similarly to DVAR, SVAR is also a gene expression regulatory network inference method based on time-series data. It is an extension of the VAR (Vector Autoregressive) method, consequently, it infers Granger causalities with the advantage that it is possible to infer and statistically test the connectivities under the following context: when the number of samples (microarrays) is lower than the number of parameters (genes), which is a very frequent condition nowadays.
SVR, DVAR and SVAR were available until now only as mathematical formulas, as described in the corresponding publications. Here, we introduce GEDI (Gene Expression Data Interpreter) Version 1.0 as an integrated software, providing easy access to the SVR, DVAR and SVAR algorithms, as well as to some other tools for gene expression data analysis. GEDI has an user-friendly interface and visualization capabilities to facilitate data interpretation.
The current version (GEDI 1.0) (see Additional file 1) runs on Windows and Linux operating systems and requires pre-installation of the R environment and of some R packages, which are freely available for downloading at . GEDI was entirely implemented in the R statistical language, being available upon GPL license. Source code, installation instructions, tutorials and some example input datasets are available at the GEDI's website. Since GEDI is an open source software, new tools can easily be added, allowing researchers the flexibility of implementing new functionalities, according to their own needs.
Results and discussion
Microarray data normalization: global and quantile  normalization methods are implemented in this version of GEDI. In addition, several normalization methods based on non-parametric regressions are also implemented, comprising the following methods: Loess , Splines [7, 8], Wavelets  and SVR (Support Vector Regression) . Also, for more than two microarrays, the cyclic normalization is performed as described in .
Identification of differentially expressed genes: t-test, t-test with permutation and the non-parametric Wilcoxon test with FDR (False Discovery Rate)  adjustment are available. Moreover, the recently published SAM (Significance Analysis of Microarray) is also available . Putative differentially expressed genes are listed from the lowest to the highest significant FDR-adjusted p-value.
Samples clustering and classification: often, a clinical interest requires identification of biomarkers, which may discriminate between pathological and normal samples. Therefore, GEDI has implemented the k-means clustering method, linear/quadratic Fisher discriminant analysis , hierarchical clustering  and the recently described SVM (Support Vector Machine) approach  with a cross-validation procedure.
Construction of gene regulatory networks: usually, it is of interest to identify which pathways the identified genes are related to. Unfortunately, depending on the treatment conditions, cell lines or tissues, these pathways have not yet been studied or are not yet known. GEDI offers some approaches to infer regulatory networks, based on gene expression data, with no a priori additional biological information. The methods employed are Pearson and Spearman partial correlation analysis to infer instantaneous associations and advanced methods based on Granger causality , such as VAR (Vector Autoregressive) , DVAR (Dynamic Vector Autoregressive)  and SVAR (Sparse Vector Autoregressive) . The VAR methods are of great interest because they allow infering Granger causalities from time-series gene expression data. DVAR may infer time (cell cycle)-varying connectivities, while SVAR may allow constructing large networks from only a few samples.
The input data format is very simple and independent of the microarray platform, i.e., it should consist of text files organized in a matrix, where each column is one microarray and each row is represented by one gene. To facilitate for the user, the input files have the same format for all functionalities.
The outputs are composed by graphics and numerical results. The plots may be saved as vectorial postscript files, allowing zoom without losing resolution. The numerical results may be saved in a plain tab delimited text file, which may be viewed using any text editor.
Normalized microarray data: The output is composed by graphical views of MA plots (raw data, normalized data and print-tips) and a file containing all the normalized expression data (Figure 2).
Differentially expressed genes: Given the FDR-adjusted p-value threshold, GEDI provides an ordered list from the lowest to the highest level of significance (the most differentially expressed genes) adjusted by FDR .
Samples clustering and classification: Statistics for each kind of analysis is provided, such as the number of corrected classified samples after cross-validation.
Gene expression regulatory networks: GEDI plots graphs which represent the regulatory networks (Figure 3). Each node of the graph represents the gene, and the edges represent the Granger causalities (VAR, DVAR and SVAR) and correlations (Pearson and Spearman). It also plots the time-varying connectivity graphic, time × connectivity plot, to visualize how the connectivity changes with time in the DVAR method (Figure 4).
GEDI is very user-friendly, since all that is required is to upload GEDI in the R environment, leading it to automatically start running. Moreover, one may easily add new functionalities and extend GEDI.
As perspectives, we intend to continue the development of GEDI by incorporating new functionalities as soon as new algorithms and statistical methods are developed to analyze gene expression data, allowing and facilitating the access to advanced methods by biomedical researchers.
GEDI was designed to be an integrated, user-friendly, viewer that combines the state of the art SVR, DVAR and SVAR algorithms. It facilitates the application of SVR, DVAR and SVAR, which used to be available only as cumbersome mathematical formulas, allowing to use them with visualizations for assessment of the results. Both running the statistical methods and visualizing of the results are accomplished within the graphical user interface, rendering the algorithms accessible to the large community of molecular biology researchers.
Availability and Requirements
Project name: GEDI
Project homepage: http://www.iq.usp.br/wwwdocentes/mcsoga/gedi/
Operating systems: Windows and Linux
Programming language: The software is implemented in R.
Other requirements: The R environment and some R packages (tcltk2, splines, kernlab, sna, wavethresh, MASS, pvclust, samr, affy) which are freely available at http://www.r-project.org.
License: The license is distributed under the GNU General Public License.
Any restrictions to use by non-academics: none
This research was supported by FAPESP, CAPES, CNPq, FINEP and PRP-USP.
- Fujita A, Sato JR, Rodrigues LO, Ferreira CE, Sogayar MC: Evaluating different methods of microarray data normalization. BMC Bioinformatics 2006, 7: 469. 10.1186/1471-2105-7-469PubMed CentralView ArticlePubMedGoogle Scholar
- Fujita A, Sato JR, Garay-Malpartida HM, Morettin PA, Sogayar MC, Ferreira CE: Time-varying modeling of gene expression regulatory networks using the wavelet dynamic vector autoregressive method. Bioinformatics 2007, 23: 1623–1630. 10.1093/bioinformatics/btm151View ArticlePubMedGoogle Scholar
- Fujita A, Sato JR, Garay-Malpartida HM, Yamaguchi R, Miyano S, Sogayar MC, Ferreira CE: Modeling gene expression regulatory networks with the sparse vector autoregressive model. BMC Systems Biology 2007, 1: 39. 10.1186/1752-0509-1-39PubMed CentralView ArticlePubMedGoogle Scholar
- The R project for statistical computing[http://www.r-project.org/]
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19: 185–193. 10.1093/bioinformatics/19.2.185View ArticlePubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addresing single and multiple slide systematic variation. Nucleic Acids Res 2002, 30: e15. 10.1093/nar/30.4.e15PubMed CentralView ArticlePubMedGoogle Scholar
- Baird D, Johnstone P, Wilson T: Normalization of microarray data using a spatial mixed model analysis which includes splines. Bioinformatics 2004, 17: 3196–3205. 10.1093/bioinformatics/bth384View ArticleGoogle Scholar
- Workman C, Jensen LJ, Jarmer H, Berka R, Gautier L, Nielser HB, Saxild HH, Nielsen C, Brunak S, Knudsen S: A new non-linear normalization method for reducing variability in DNA microarray experiments. Genome Biology 2002., 3:Google Scholar
- Wang J, Ma JZ, Li MD: Normalization of cDNA microarray data using wavelet regressions. Combinatorial Chemistry & High Throughput Screening 2004, 9: 783–791. 10.2174/1386207043328274View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A practical and powerful approach to multiple testing. J R Statist Soc B 1995, 57: 289–300.Google Scholar
- Tusher V, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. PNAS 2001, 98: 5116–5121. 10.1073/pnas.091062498PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher RA: The use of multiple measurements in taxonomic problems. Annals of Eugenics 1936, 7: 179–188.View ArticleGoogle Scholar
- Suzuki R, Shimodaira H: An application of multiscale bootstrap resampling to hierarchical clustering of microarray data: How accurate are these clusters? The Fifteenth International Conference on Genome Informatics 2004, 2004: 34.Google Scholar
- Yeang CH, Ramaswamy S, Tamayo P, Mukherjee S, Rifkin RM, Angelo M, Reich M, Lander E, Mesirov J, Golub T: Molecular classification of multiple tumor types. Bioinformatics 2001, 17: 316–322.View ArticleGoogle Scholar
- Granger CWJ: Investigating causal relation by econometric and cross-sectional method. Econometrica 1969, 37: 424–438. 10.2307/1912791View ArticleGoogle Scholar
- Mukhopadhyay ND, Chatterjee S: Causality and pathway search in microarray time series experiment. Bioinformatics 2007, 23: 442–449. 10.1093/bioinformatics/btl598View ArticlePubMedGoogle Scholar
This article is published under license to 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.