How do I download a list of genes involved in cellular metabolism (in humans) based on a Reactome pathway?
1
0
Entering edit mode
10 months 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 • 493 views
ADD COMMENT
2
Entering edit mode
10 months ago
Pratik ▴ 800

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

Login before adding your answer.

Traffic: 2168 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