Pathway analysis explained
Pathway is the term from molecular biology for a curated schematic representation of a well characterized segment of the molecular physiological machinery, such as a metabolic pathway describing an enzymatic process within a cell or tissue or a signaling pathway model representing a regulatory process that might, in its turn, enable a metabolic or another regulatory process downstream. A typical pathway model starts with an extracellular signaling molecule that activates a specific receptor, thus triggering a chain of molecular interactions.[1] A pathway is most often represented as a relatively small graph with gene, protein, and/or small molecule nodes connected by edges of known functional relations. While a simpler pathway might appear as a chain,[2] complex pathway topologies with loops and alternative routes are much more common. Computational analyses employ special formats of pathway representation.[3] [4] In the simplest form, however, a pathway might be represented as a list of member molecules with order and relations unspecified. Such a representation, generally called Functional Gene Set (FGS), can also refer to other functionally characterised groups such as protein families, Gene Ontology (GO) and Disease Ontology (DO) terms etc. In bioinformatics, methods of pathway analysis might be used to identify key genes/proteins within a previously known pathway in relation to a particular experiment / pathological condition or building a pathway de novo from proteins that have been identified as key affected elements. By examining changes in e.g. gene expression in a pathway, its biological activity can be explored. However most frequently, pathway analysis refers to a method of initial characterization and interpretation of an experimental (or pathological) condition that was studied with omics tools or genome-wide association study.[5] Such studies might identify long lists of altered genes. A visual inspection is then challenging and the information is hard to summarize, since the altered genes map to a broad range of pathways, processes, and molecular functions (with a large gene fraction lacking any annotation). In such situations, the most productive way of exploring the list is to identify enrichment of specific s in it. The general approach of enrichment analyses is to identify FGSs, members of which were most frequently or most strongly altered in the given condition, in comparison to a gene set sampled by chance. In other words, enrichment can map canonical prior knowledge structured in the form of FGSs to the condition represented by altered genes.
Use
The data for pathway analysis come from high throughput biology. This includes high throughput sequencing data and microarray data. Before pathway analysis can be done, each gene's alteration should be evaluated using the omics dataset in either quantitative (differential expression analysis) or qualitative (detection of somatic point mutations or mapping neighbor genes to a disease-associated SNP). It is also possible to combine datasets from different research groups or multiple omics platform with a meta-analysis and cross-platform regularization.[6] [7] Further, a list where gene identifiers are accompanied by the alteration attributes is subjected to a pathway analysis. By using pathway analysis software, researchers can determine which s are enriched with the altered experimental genes[8] [9] For example, pathway analysis of several independent microarray experiments (meta-analysis) helped to discover potential biomarkers in a single pathway important for fast-to-slow switch fiber type transition in Duchenne muscular dystrophy.[10] In another study meta-analysis identified two biomarkers in blood of patients with Parkinson's disease, which can be useful for monitoring the disease.[11] Candidate gene alleles causative of Alzheimer's disease and elderly dementia where first discovered via genome-wide association study and further validated with network enrichment analysis against consisting of known Alzheimer's genes.[12] [13]
Databases
Pathway collections and interaction networks constitute the knowledge base required for a pathway analysis. Pathway content, structure, format, and functionality vary between different database resources such as KEGG,[14] WikiPathways, or Reactome.[15] Also exist proprietary pathways collections used by e.g. Pathway Studio[16] and Ingenuity Pathway Analysis[17] tools. Public online tools can provide pre-compiled and ready-to-go menus of pathways and networks from different open sources (e.g. EviNet).
Methods and software
Pathway analysis software can be found in the form of desktop programs, web-based applications, or packages coded in such languages as R and Python and shared openly through the BioConductor[18] and GitHub[19] projects. The methodology of pathway analysis evolves fast and the classification is still discussable,[20] with the following main categories of pathway enrichment analysis applicable to high-throughput data:
Over-representation analysis (ORA)
This method measures the overlap between, on the one hand, a set of genes (or proteins) in an and, on the other hand, a list of most altered genes generally called Altered Gene Sets (AGS). A typical AGS example is a list of top N differentially expressed genes from an RNA-Seq assay. The basic assumption behind ORA is that a biologically relevant pathway can be identified by excess of genes in it compared to the number expected by chance. The aim of ORA is to identify such enriched pathways, judging by statistical significance of the overlap between FGS and AGS as determined either by an appropriate statistic, such as Jaccard index or by a statistical test producing p-values (Fisher's exact test or the test using hypergeometric distribution).
Functional class scoring (FCS)
This method identifies by considering their relative positions in the full list of genes studied in the experiment. This full list should be therefore ranked in advance by a statistic (such as mRNA expression fold-change, Student's t-test etc.) or a p-value - while watching the direction of fold change, since p-values are non-directional. Thus FCS takes into account every FGS gene regardless of its statistical significance and does not require pre-compiled . One of the first and most popular methods deploying the FCS approach was the Gene Set Enrichment Analysis (GSEA).
Pathway topology analysis (PTA)
Similarly to, PTA accounts for high-throughput data for every gene.[21] In addition, specific topological information is used about role, position, and interaction directions of the pathway genes. This requires additional input data from a pathway database in a pre-specified format, such as KEGG Markup Language (KGML). Using this information, PTA estimates a pathway significance by considering how much each individual gene alteration might have affected the whole pathway. Multiple alteration types can be used in parallel (somatic copy-number variations, point mutations etc.) when available.[22] The set of PTA methods includes the Impact Analysis,[23] [24] EnrichNet,[25] GGEA,[26] and TopoGSA.[27]
Network enrichment analysis (NEA)
Network enrichment analysis (NEA) has been an extension of gene-set enrichment analysis to the domain of global gene networks[28] [29] [30] [31] The major principle of NEA can be understood in comparison with, where enrichment of in genes of the is determined by how many genes are directly shared by AGS and FGS. In NEA, on the contrary, the global network is searched for network edges that connect any genes of AGS with any genes of FGS. Since enrichment significance is influenced by the highly variable node degrees of individual AGS and FGS genes, it should be determined by a dedicated statistical test, which compares the observed number of network edges to the number expected by chance in the same network context. Some valuable properties of NEA are that:
- it is more robust to biological and technical variability between sample replicates;[7] [32]
- genes may not necessarily be annotated as pathway members;[33]
- members do not have to be altered themselves, but still are accounted for due to possessing network links to AGS genes.[34]
Commercial solutions
Beyond open-source tools, such as STRING or Cytoscape, a number of companies sell licensed software products to analyse gene sets. While most of the publicly available solutions use online and public pathway collections, the commercial products mostly promote own, proprietary pathways and networks. The choice of such products might be driven by customers' skills, financial and time resources, and needs. Ingenuity, for example, maintains a knowledge base for comparative analysis of gene expression data.[35] Pathways Studio[36] is commercial software which allows searching for biologically relevant facts, analyze experiments, and create pathways. Pathways Studio Viewer[37] is a free resource from the same company for presenting the Pathway Studio interactive pathway collection and database. Two commercial solutions offer : iPathwayGuide from Advaita Corporation and MetaCore from Thomson Reuters.[38] Advaita uses the peer reviewed Impact Analysis method while the MetaCore method is unpublished. Correlation Engine uses the Running Fisher algorithm for gene set enrichment within its Pathway Enrichment application.[39]
Limitations
Lack of annotations
Application of pathway analysis methods depends on annotations found in existing databases, such as gene set membership in pathways, pathway topology, presence of genes in the global network etc. These annotations, however, are far from being complete and have highly variable degrees of confidence. In addition, such information is usually general, i.e. deprived of e.g. cell type, compartment, or developmental context. Therefore, interpretation of pathway analysis results for omics datasets should be done with caution Partially, the problem can be addressed by analysing larger gene sets in a more, such as big pathway collections or global interaction networks.[40]
See also
Notes and References
- Book: Berg JM, Tymoczko JL, Stryer L . Biochemistry . 2002 . W.H. Freeman . New York . 978-0-7167-3051-4 . 5th.
- Ohlrogge J, Browse J . Lipid biosynthesis . The Plant Cell . 7 . 7 . 957–70 . July 1995 . 7640528 . 10.1105/tpc.7.7.957 . 160893 . 219201001 . free .
- Web site: Main Page - SBML.caltech.edu . sbml.org . en.
- Web site: KGML (KEGG Markup Language) . www.genome.jp.
- García-Campos MA, Espinal-Enríquez J, Hernández-Lemus E . Pathway Analysis: State of the Art . Frontiers in Physiology . 6 . 383 . 2015 . 26733877 . 4681784 . 10.3389/fphys.2015.00383 . free .
- Walsh CJ, Hu P, Batt J, Santos CC . Microarray Meta-Analysis and Cross-Platform Normalization: Integrative Genomics for Robust Biomarker Discovery . Microarrays . 4 . 3 . 389–406 . August 2015 . 27600230 . 4996376 . 10.3390/microarrays4030389 . free .
- Suo C, Hrydziuszko O, Lee D, Pramana S, Saputra D, Joshi H, Calza S, Pawitan Y . Integration of somatic mutation, expression and functional data reveals potential driver genes predictive of breast cancer survival . Bioinformatics . 31 . 16 . 2607–13 . August 2015 . 25810432 . 10.1093/bioinformatics/btv164 . free .
- Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM . Systematic determination of genetic network architecture . Nature Genetics . 22 . 3 . 281–5 . July 1999 . 10391217 . 10.1038/10343 . 14688842 .
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP . Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles . Proceedings of the National Academy of Sciences of the United States of America . 102 . 43 . 15545–50 . October 2005 . 16199517 . 1239896 . 10.1073/pnas.0506580102 . 2005PNAS..10215545S . free .
- Kotelnikova E, Shkrob MA, Pyatnitskiy MA, Ferlini A, Daraselia N . Novel approach to meta-analysis of microarray datasets reveals muscle remodeling-related drug targets and biomarkers in Duchenne muscular dystrophy . PLOS Computational Biology . 8 . 2 . e1002365 . February 2012 . 22319435 . 3271016 . 10.1371/journal.pcbi.1002365 . 2012PLSCB...8E2365K . free .
- Santiago JA, Potashkin JA . Judith Potashkin. Network-based metaanalysis identifies HNF4A and PTBP1 as longitudinally dynamic biomarkers for Parkinson's disease . Proceedings of the National Academy of Sciences of the United States of America . 112 . 7 . 2257–62 . February 2015 . 25646437 . 4343174 . 10.1073/pnas.1423573112 . 2015PNAS..112.2257S . free .
- Reynolds CA, Hong MG, Eriksson UK, Blennow K, Wiklund F, Johansson B, Malmberg B, Berg S, Alexeyenko A, Grönberg H, Gatz M, Pedersen NL, Prince JA . Analysis of lipid pathway genes indicates association of sequence variation near SREBF1/TOM1L2/ATPAF2 with dementia risk . Human Molecular Genetics . 19 . 10 . 2068–78 . May 2010 . 20167577 . 2860895 . 10.1093/hmg/ddq079 .
- Bennet AM, Reynolds CA, Eriksson UK, Hong MG, Blennow K, Gatz M, Alexeyenko A, Pedersen NL, Prince JA . Genetic association of sequence variants near AGER/NOTCH4 and dementia . Journal of Alzheimer's Disease . 24 . 3 . 475–84 . 1 January 2011 . 21297263 . 3477600 . 10.3233/jad-2011-101848 .
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M . KEGG: Kyoto Encyclopedia of Genes and Genomes . Nucleic Acids Research . 27 . 1 . 29–34 . January 1999 . 9847135 . 148090 . 10.1093/nar/27.1.29 .
- Vastrik I, D'Eustachio P, Schmidt E, Joshi-Tope G, Gopinath G, Croft D, de Bono B, Gillespie M, Jassal B, Lewis S, Matthews L, Wu G, Birney E, Stein L . Reactome: a knowledge base of biologic pathways and processes . Genome Biology . 8 . 3 . R39 . 2007 . 17367534 . 1868929 . 10.1186/gb-2007-8-3-r39 . free .
- https://mammalcedfx.pathwaystudio.com/app/search Pathway Studio Pathways
- https://www.qiagen.com/us/shop/genes-and-pathways/pathway-central/ Pathway Central
- 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 Biology . 5 . 10 . R80 . 2004 . 15461798 . 545600 . 10.1186/gb-2004-5-10-r80 . Sandrine Dudoit . free .
- Book: Dabbish L, Stuart C, Tsay J, Herbsleb J . Social coding in GitHub: transparency and collaboration in an open software repository. . Proceedings of the ACM 2012 Conference on Computer Supported Cooperative Work . February 2012 . 1277–1286 . https://www.cs.cmu.edu/afs/cs/Web/People/xia/resources/Documents/cscw2012_Github-paper-FinalVersion-1.pdf . Association for Computing Machinery . New York . 10.1145/2145204.21453 . 2024-07-21 .
- Henderson-Maclennan NK, Papp JC, Talbot CC, McCabe ER, Presson AP . Pathway analysis software: annotation errors and solutions . Molecular Genetics and Metabolism . 101 . 2–3 . 134–40 . 2010 . 20663702 . 2950253 . 10.1016/j.ymgme.2010.06.005 .
- Emmert-Streib F, Dehmer M . Networks for systems biology: conceptual connection of data and function . IET Systems Biology . 5 . 3 . 185–207 . May 2011 . 21639592 . 10.1049/iet-syb.2010.0025 .
- Khatri P, Sirota M, Butte AJ . Ten years of pathway analysis: current approaches and outstanding challenges . PLOS Computational Biology . 8 . 2 . e1002375 . 23 February 2012 . 22383865 . 3285573 . 10.1371/journal.pcbi.1002375 . 2012PLSCB...8E2375K . free .
- Draghici S, Khatri P, Tarca AL, Amin K, Done A, Voichita C, Georgescu C, Romero R . A systems biology approach for pathway level analysis . Genome Research . 17 . 10 . 1537–45 . October 2007 . 17785539 . 1987343 . 10.1101/gr.6202607 .
- Tarca AL, Draghici S, Khatri P, Hassan SS, Mittal P, Kim JS, Kim CJ, Kusanovic JP, Romero R . A novel signaling pathway impact analysis . Bioinformatics . 25 . 1 . 75–82 . January 2009 . 18990722 . 2732297 . 10.1093/bioinformatics/btn577 .
- Glaab E, Baudot A, Krasnogor N, Schneider R, Valencia A . EnrichNet: network-based gene set enrichment analysis . Bioinformatics . 28 . 18 . i451–i457 . September 2012 . 22962466 . 3436816 . 10.1093/bioinformatics/bts389 .
- Geistlinger L, Csaba G, Küffner R, Mulder N, Zimmer R . From sets to graphs: towards a realistic enrichment analysis of transcriptomic systems . Bioinformatics . 27 . 13 . i366-73 . July 2011 . 21685094 . 3117393 . 10.1093/bioinformatics/btr228 .
- Glaab E, Baudot A, Krasnogor N, Valencia A . TopoGSA: network topological gene set analysis . Bioinformatics . 26 . 9 . 1271–2 . May 2010 . 20335277 . 2859135 . 10.1093/bioinformatics/btq131 .
- Shojaie A, Michailidis G . Network enrichment analysis in complex experiments . Statistical Applications in Genetics and Molecular Biology . 9 . 1 . Article22 . 22 May 2010 . 20597848 . 2898649 . 10.2202/1544-6115.1483 .
- Huttenhower C, Haley EM, Hibbs MA, Dumeaux V, Barrett DR, Coller HA, Troyanskaya OG . Exploring the human genome with functional maps . Genome Research . 19 . 6 . 1093–106 . June 2009 . 19246570 . 10.1101/gr.082214.108 . 2694471 . free .
- Alexeyenko A, Lee W, Pernemalm M, Guegan J, Dessen P, Lazar V, Lehtiö J, Pawitan Y . Network enrichment analysis: extension of gene-set enrichment analysis to gene networks . BMC Bioinformatics . 13 . 226 . September 2012 . 22966941 . 3505158 . 10.1186/1471-2105-13-226 . free .
- Signorelli M, Vinciotti V, Wit EC . NEAT: an efficient network enrichment analysis test . BMC Bioinformatics . 17 . 1 . 352 . September 2016 . 27597310 . 10.1186/s12859-016-1203-6 . 5011912 . 1604.01210 . 2274758 . free .
- Jeggari A, Alexeyenko A . NEArender: an R package for functional interpretation of 'omics' data via network enrichment analysis . BMC Bioinformatics . 18 . Suppl 5 . 118 . March 2017 . 28361684 . 5374688 . 10.1186/s12859-017-1534-y . free .
- Hong MG, Alexeyenko A, Lambert JC, Amouyel P, Prince JA . Genome-wide pathway analysis implicates intracellular transmembrane protein transport in Alzheimer disease . Journal of Human Genetics . 55 . 10 . 707–9 . October 2010 . 20668461 . 10.1038/jhg.2010.92 . 27020289 . free .
- Jeggari A, Alekseenko Z, Petrov I, Dias JM, Ericson J, Alexeyenko A . EviNet: a web platform for network enrichment analysis with flexible definition of gene sets . Nucleic Acids Research . 46 . W1 . W163–W170 . July 2018 . 29893885 . 6030852 . 10.1093/nar/gky485 .
- Web site: Ingenuity IPA - Integrate and Understand Complex 'omics Data. . Ingenuity . 8 April 2015 .
- http://www.pathwaystudio.com/ Pathway Studio
- https://mammalcedfx.pathwaystudio.com/app/search Pathway Studio Viewer
- Mitrea C, Taghavi Z, Bokanizad B, Hanoudi S, Tagett R, Donato M, Voichiţa C, Drăghici S . Methods and approaches in the topology-based analysis of biological pathways . Frontiers in Physiology . 4 . 278 . October 2013 . 24133454 . 3794382 . 10.3389/fphys.2013.00278 . free .
- Kupershmidt I, Su QJ, Grewal A, Sundaresh S, Halperin I, Flynn J, Shekar M, Wang H, Park J, Cui W, Wall GD, Wisotzkey R, Alag S, Akhtari S, Ronaghi M . Ontology-based meta-analysis of global collections of high-throughput public data . PLOS ONE . 5 . 9 . e13066 . September 2010 . 20927376 . 2947508 . 10.1371/journal.pone.0013066 . free . 2010PLoSO...513066K . Aziz RK .
- Franco M, Jeggari A, Peuget S, Böttger F, Selivanova G, Alexeyenko A . Prediction of response to anti-cancer drugs becomes robust via network integration of molecular data . Scientific Reports . 9 . 1 . 2379 . February 2019 . 30787419 . 6382934 . 10.1038/s41598-019-39019-2 . 2019NatSR...9.2379F .