Repositories of genome-wide expression studies such as ArrayExpress  have been growing rapidly over the last few years and continue to do so. The more experimental data are deposited into these repositories, the more likely it becomes that some of them can provide a meaningful biological context to aid in the planning and analysis of new studies. Retrieval of experiments based on their textual description and experimental design has several shortcomings. First of all, textual description of an experiment or its results is not as information-rich as the actual data itself. Secondly, information about the experimental design alone is only of limited use in retrieving biologically relevant data because it does not reflect the results, which contain the bulk of the information and may reveal unexpected relationships. We introduce novel retrieval methods that incorporate the actual gene expression measurements into the search process, along with visualization tools for interpreting and exploring the results .
We developed a two-stage procedure, first identifying differentially active gene sets in each experiment using a recent nonparametric statistical method , and then combining gene set activation patterns into higher-level structures, so-called biological topics, using a state-of-the-art probabilistic model . The probabilistic formulation enables the use of a natural and rigorous metric for assessing the similarity between two experiments. For interpreting and exploring retrieval results, we have developed visualization methods that also provide insight into the model used to perform the retrieval.
We show that gene sets corresponding to each biological topic form highly coherent and holistic components. Several case studies performed on a subset of ArrayExpress show that our method can retrieve experiments relevant to a biological question, as long as sufficient amounts of data are available, and highlight relations between experiments, either because the same biological questions were targeted, or because of unexpected relationships that were confirmed in the literature. The visualization methods allow us to both efficiently interpret the model and put retrieval results in the context of the whole set of experiments (see Figure 1 for an example).
Figure 1. 2D NeRV projection of retrieval results when the model is queried with a malignant melanoma experiment. Each experiment is represented as a striped glyph. Colors indicate biological topics. Stripe widths indicate the predominance of each biological topic in each experiment. Glyph size indicates relevance to the malignant melanoma query experiment .
Using a combination of existing and novel methods for modeling and visualizing a heterogeneous collection of gene expression experiments, we were able to decompose and relate experiments via biologically meaningful components. Our approach allows search within a gene expression database to be driven by actual measurement data.
This work was supported by TEKES (grant no. 40101/07). JC, AF and SK are additionally partially supported by PASCAL 2 Network of Excellence, ICT 216886. JC is additionally supported by a doctoral grant from the Portuguese Foundation for Science and Technology (FCT). NG is supported by a PhD fellowship of the European Molecular Biology Laboratory (EMBL).
Parkinson H, Kapushesky M, Kolesnikov N, Rustici G, Shojatalab M, Abeygunawardena N, Brube H, Dylab M, Emam I, Farne A: ArrayExpress update – from an archive of functional genomics experiments to the atlas of gene expression.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gilette MA, Paulovich A, Pomeroy SL, Golub TR, et al.: Gene set enrichment analysis – a knowledge-based approach for interpreting genome-wide expression profiles.
J Mach Learn Res 2003, 3:993-1022. Publisher Full Text