How do I download a list of genes involved in cellular metabolism (in humans) based on a Reactome pathway?
1
0
Entering edit mode
2.3 years ago
emmah • 0

I am trying to get a list of genes involved in "metabolism", similar to in this paper - https://www.nature.com/articles/s42003-019-0666-1.

In this paper, they state that they "accessed information of all human metabolic pathways from the Reactome pathways database version 68. Reactome pathways are arranged into several tiers with the Reactome term “metabolism” (Reactome ID: R-HSA-1430728), encompassing 68 different metabolic pathways (see https://reactome.org/PathwayBrowser/#/R-HSA-1430728). The first-tier pathways include sixteen curated metabolic pathways which involve 2325 genes."

How do I do something similar and download a list of all genes involved in metabolism?

Also any explanation of the "tiers" thing would be great! I.e. when the authors talk about "first-tier pathways", what does this actually mean?

genes metabolism GO Reactome • 1.1k views
ADD COMMENT
3
Entering edit mode
2.3 years ago
Pratik ★ 1.0k

To answer this question: How do I do something similar and download a list of all genes involved in metabolism?

My go-to is usually the Comparative Toxicogenomics Database.

You can download the data files here in csv and tsv formats.

If for example you download the csv.gz file. You can extract it (gunzip) and then grep it for your search term using something like this:

grep R-HSA-1430728 your/file/path/here/CTD_genes_pathways.csv | cut -d ',' -f1 | sort | uniq

The command will search for R-HSA-1430728, cut the gene out (cutting the first word detecting commas as the separator), sort the genes alphabetically, and then give only the unique gene names (removes deplicates).

Note: I obtained 2171 unique genes using the command above (not 2325 genes):

grep R-HSA-1430728 your/file/path/here/CTD_genes_pathways.csv | cut -d ',' -f1 | sort | uniq | wc -l
    2171

Hope this helps!

ADD COMMENT
0
Entering edit mode

Hi Pratik, Thanks for the helpful solution above. I'm working on a similar task in Mouse species and noticed the CTD database is human-specific. I would appreciate if you have suggestions for an alternative in Mouse data.

My original question is also posted here: Mapping Mouse RNAseq Marker Genes to Reactome Pathways

ADD REPLY

Login before adding your answer.

Traffic: 2018 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6