Question: Extract Pathway Topological Graphs (With Metabolites) from Reactome and WikiPathways
1
gravatar for taraeicher
4 months ago by
taraeicher30
National Institute for Advancing Translational Sciences
taraeicher30 wrote:

I am looking for a way to extract pathway-specific graphs from Reactome and WikiPathways, including both genes and metabolites in the graph. My goal is to test a novel topological pathway enrichment analysis technique on these graphs. I'm looking for something analogous to what KEGGgraph does for KEGG pathways, but for Reactome or WikiPathways. I know that the graphite package extracts graphs from various pathway databases, but it only includes genes, and the method I wish to test leverages both gene and metabolite data. Anyone know how I might go about this?

ADD COMMENTlink modified 3 months ago by dslentermsc10 • written 4 months ago by taraeicher30
1
gravatar for xanderpico
3 months ago by
xanderpico330
xanderpico330 wrote:

Reactome and WikiPathways both have apps for Cytoscape that pull their content into a graph model:

In the case of the WikiPathways app, you have the explicit option of importing a pathway as "a pathway" (meaning a complete, annotated diagram) or as "a network" (meaning a simplified graph diagram). The network version is ideal for graph analytics, merging and extension.

The app exposes a REST API via Cytoscape as well, so you can programmatically import hundreds of (or all) pathways via an R or Python script. See RCy3 and py4cytoscape packages.

By the way, all the human pathways from Reactome are also hosted by WikiPathways. So, using the WikiPathways app, you can also pull in Reactome pathways (for human) into Cytoscape "as networks" and have both in exactly the same data structure.

Caveat: Not all graph models are the same! For a topological pathway analysis, you will want to very carefully assess how pathways are modeled by these tools. For example, whether genes and metabolites interact directly or via reaction nodes, which would completely change measures of path length, diffusion, etc. In the case of WikiPathways, there is an RDF model of pathwys, for example, that is different from the Cytoscape model you get via the app. You may want to pull content in via RDF depending on your analysis plan. Or you may want to restructure the imported graphs via the app.

ADD COMMENTlink written 3 months ago by xanderpico330
1
gravatar for dslentermsc
3 months ago by
dslentermsc10
dslentermsc10 wrote:

Also check out this paper: https://dx.doi.org/10.3389%2Ffgene.2019.00059 ("Interactions with at least two annotated interaction participants (gene product, metabolite, complex) are included. Gene products have unified Ensembl (Zerbino et al., 2017) identifiers, metabolites have either Wikidata (Mietchen et al., 2015), ChEBI (Hastings et al., 2015) or HMDB identifiers (Wishart et al., 2017), and complexes have Reactome identifiers. ").

ADD COMMENTlink written 3 months ago by dslentermsc10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1293 users visited in the last hour