This is a test version of Biostars. For the public version, visit https://www.biostars.org.
retrieving pathways belong to my gene names

Hi there!

Sorry guys, if my question has repeated before, I don't want waste your time, but I couldn't find proper answer for my question.

I am working with numbers of genes, I want to know to which pathway of KEGG database they belong? Indeed, I have written a program by python, using bioservices module, but there isn't any way to search several pathway in a run.how can I solve it by python? if anyone could help me, I appreciate it.

python kegg bioservices
• 536 views • link
written by EFN • EFN

use kEGG rest api and post some example gene symbols

• 1 views • link
written by cpad0112 cpad0112

it is time-consuming way I think, need a faster way

can you help me?

• 1 views • link
written by EFN • EFN

post some example gene symbols

• 1 views • link
written by cpad0112 cpad0112

some of my genes are : CRK, RUNX1, PSMD11,JUP , .......

Thanks

• 1 views • link
written by EFN • EFN

It is unclear what exactly you have done so far. It sounds like you may have written a python program (that uses KEGG API?, since you used those tags) to do the search you need but you just want it done for several pathways. You may have to do some kind of a loop in that case but be aware that rapid programmatic querying may result in a IP ban from KEGG. You are supposed to buy a license (even for academic users) for extensive KEGG use.

• 1 views • link
written by GenoMax GenoMax

thanks for your response! No, I used bioservices module, I want to know if anyone can help me using API. maybe using loops is good idea.

• 1 views • link
written by EFN • EFN

Format the output as per your needs:

> from bioservices.kegg import KEGG
> k = KEGG()
> genes=["CRK","RUNX1","PSMD11","JUP"]
> for i in genes:
        j=k.get_pathway_by_gene(i,"hsa")
        for keys,values in j.items():
            print(i,keys, values, sep="\t")

output: 

CRK hsa04010    MAPK signaling pathway
CRK hsa04012    ErbB signaling pathway
CRK hsa04015    Rap1 signaling pathway
CRK hsa04062    Chemokine signaling pathway
CRK hsa04510    Focal adhesion
CRK hsa04666    Fc gamma R-mediated phagocytosis
CRK hsa04722    Neurotrophin signaling pathway
CRK hsa04810    Regulation of actin cytoskeleton
CRK hsa04910    Insulin signaling pathway
CRK hsa05100    Bacterial invasion of epithelial cells
CRK hsa05131    Shigellosis
CRK hsa05200    Pathways in cancer
CRK hsa05206    MicroRNAs in cancer
CRK hsa05211    Renal cell carcinoma
CRK hsa05220    Chronic myeloid leukemia
RUNX1   hsa04530    Tight junction
RUNX1   hsa04659    Th17 cell differentiation
RUNX1   hsa05200    Pathways in cancer
RUNX1   hsa05202    Transcriptional misregulation in cancer
RUNX1   hsa05220    Chronic myeloid leukemia
RUNX1   hsa05221    Acute myeloid leukemia
PSMD11  hsa03050    Proteasome
PSMD11  hsa05169    Epstein-Barr virus infection
JUP hsa05200    Pathways in cancer
JUP hsa05202    Transcriptional misregulation in cancer
JUP hsa05221    Acute myeloid leukemia
JUP hsa05226    Gastric cancer
JUP hsa05412    Arrhythmogenic right ventricular cardiomyopathy (ARVC)
• 1 views • link
written by cpad0112 cpad0112

thanks for your response ;)

it is really helpful

• 1 views • link
written by EFN • EFN

Hi again,

new question has arisen for me when I import a gene name that kegg cannot recognize it... the error is : if 'PATHWAY' in dic.keys(): AttributeError: 'int' object has no attribute 'keys'

I wrote an if statement, but the same error has occurred. where I'm doing wrong?

• 1 views • link
written by EFN • EFN

well, it is a generic error by bioservices. If it can't find a match, it throws a general error. This is not an error in your input/code.

• 1 views • link
written by cpad0112 cpad0112

1 answer

GeneSCF will be handy for this task (NOT python).

• 1 views • link
written by EagleEye EagleEye

Log in to answer this question.

More posts like this