Pathway Based Microarray Analysis, Utilising Enzyme Compounds and Cascade Events

Summary

Background: Pathway based microarray analysis is an effort to integrate microarray and pathway data in a holistic analytical approach, looking for coordinated changes in the expression of sets of genes forming pathways. However, it has been observed that the results produced are often cryptic, with cases of closely related genes in a pathway showing quite variable, even opposing expression.

Objectives: We propose a methodology to identify the state of activation of individual pathways, based on our hypothesis that gene members of many pathways or modules exhibit differential expression that results from their contribution to any combination of all their constituent pathways. Therefore, the observed expression of such a gene does not necessarily imply the activation state of a given pathway where its product participates, but reflects the net expression resulting from its participation in all its constituent pathways.

Methods: Firstly, in an effort to validate the hypothesis, we split the genes into two groups; single and multi-membership. We then determined and compared the proportion of differentially expressed genes in each group, for each experiment. In addition, we estimated the cumulative binomial probability of observing as many or more expressed genes in each group, in each experiment, simply by chance. Second, we propose a hill climbing methodology, aiming to maximise the agreement of gene expression per module.

Results: We detected more frequent expression of multi-membership genes and significantly lower probabilities of observing such a high proportion of differentially expressed multi-membership genes, as the one present in the dataset. The algorithm was able to correctly identify the state of activation of the KEGG glycolysis and gluconeogenesis modules, using a number of Saccharomyces cerevisiae datasets. We show that the result is equivalent to the best solution found following exhaustive search.

Conclusions: The proposed method takes into account the multi-membership nature of genes and our knowledge of the competitive nature of our exemplar modules, revealing the state of activity of a pathway.

Keywords

Microarray analysis - biochemical pathways - computational biology - computing methodologies

• References

• 1 Werner T. Bioinformatics applications for pathway analysis of microarray data. Current Opinion in Biotechnology 2008; 19 (01) 50-54.
• 2 Abelson H, Sussman GJ, Sussman J. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. PNAS 2005; 102: 15545-15550.
• 3 Keller A, Backes C, Gerasch A, Kaufman M, kohlbacher O, Meese E, Lehnof HP. A novel algorithm for detecting differentially expressed paths based on gene set enrichment analysis. Bioinformatics 2009; 25 (21) 2787-2794.
• 4 Kanehisa M, Goto S. KEGG: Kyoto encyclopaedia of genes and genomes. Nucl Acid Res 2000; 28: 27-30.
• 5 Caspi R, Altman T, Dale JM, Dreher K, Fulcher CA, Gilham F, Kaipa P, Karthikeyan AS, Kothari A, Popescu L, Pujar A, Shearer AG, Zhang P, Karp PD. The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res 2010; 38 (Suppl. 01) D473-D479.
• 6 Joshi-Tope G, Gillespie M, Vastrik I, D'Eustachio P, Schmidt E, de Bono D, Jassal B, Gopinath GR, Wu GR, Matthews L, Lewis S, Birney E, Stein L. Reactome: a knowledge base of biological pathways. Nucl Acid Res 2005; 33 Database issue D428-D432.
• 7 Cavalieri D, Castagnini C, Toti S, Maciag K, Kelder T, Gambineri L, Angioli S, Dolara P. Eu.Gene Analyzer a tool for integrating gene expression data with pathway databases. Bioinformatics 2007; 23 (19) 2631-2632.
• 8 Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR. GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways. Nat Genet 2002; 31 (01) 19-20.
• 9 Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Research 2003; 13 (11) 2498-2504.
• 10 Quadroni M, James P. Proteomics and automation. Electrophoresis 1999; 20: 664-677.
• 11 Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R. NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Research 2009; 37 (Suppl. 01) D1-D15.
• 12 Pavlidis S, Payne A, Swift S. An Improved Methodology for Pathway Based Microarray Analysis Based on Identification of Individual Pathways Responsible for Gene Regulation. IDAMAP 2008: 69-74.
• 13 Pavlidis S, Payne A, Swift S. Multi-membership Gene Regulation in Pathway Based Microarray Analysis. Algorithms for Molecular Biology 2011; 6: 22
• 14 Webb EC. Enzyme nomenclature 1992; recommendations of the Nomenclature Committee of the International Union of Biochemistry and Molecular Biology on the nomenclature and classification of enzymes. International Union of Biochemistry and Molecular Biology. Academic Press, San Diego; 1992.
• 15 Boles E, Schulte F, Miosga T, Freidel K, Schulter E, Zimmermann FK, Hollenberg CP, Heinisch JJ. Characterization of a glucose-repressed pyruvate kinase (Pyk2p) in Saccharomyces cerevisiae that is catalytically insensitive to fructose-1,6-bisphosphate. J Bacteriol 1997; 179 (09) 2987-2993.
• 16 Delgado ML, O’Connor JE, Azorin I, Renau-Piqueras J, Gil ML, Gozalbo D. The glyceraldehyde-3-phosphate dehydrogenase polypeptides encoded by the Saccharomyces cerevisiae TDH1, TDH2 and TDH3 genes are also cell wall proteins. Microbiology 2001; 147 (02) 411-417.
• 17 Champe PC, Harvey RA, Ferrier DR. Lippincotts Illustrated Reviews Biochemistry. 3rd ed.. Lippincott reverend & adventurer 2004.
• 18 Swift S, Tucker A, Vinciotti V, Martin N, Orengo C, Liu X, Kellam P. Consensus clustering and functional interpretation of gene-expression data. Genome Biology 2004; 5: R49
• 19 Michalewicz Z, Fogel DB. How to solve it: Modern heuristics. Berlin: Springer; 1998.
• 20 Herrero P, Galíndez J, Ruiz N, Martínez-Campa C, Moreno F. Transcriptional regulation of the Saccharomyces cerevisiae HXK1, HXK2 and GLK1 genes. Yeast 1995; 11 (02) 137-144.
• 21 DeRisi JL, Iyer VR, Brown PO. Exploring the metabolic and genetic control of gene expression on a genomic scale. Science 1997; 278: 680-686.
• 22 Ronen M, Botstein D. Transcriptional response of steady-state yeast cultures to transient perturbations in carbon source. Proc Natl Acad Sci U.S. A 2006; 103 (02) 389-394.
• 23 Gasch AP, Spellman PT, Kao CM, Carmel-Harel O. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell 2000; 11 (12) 4241-4257.
• 24 Bradley P, Brauer M, Rabinowitz J, Troyanskaya O. Coordinated concentration changes of transcripts and metabolites in Saccharomyces cerevisiae. PLoS Computational Biology 2006; 5 (01) e1000270