KEGG pathways associated to a gene using BiomaRt
4
5
Entering edit mode
10.4 years ago
TitoPullo ▴ 190

Hello,

is there an easy way to get the get the list of KEGG pathways associated to a specific gene using biomaRt?

I found a solution to retrieve the GO terms: How Do I Do Simple Go Term Lookup Given A Gene (Or Mrna) Identifier?.. Is it possible use something similar for KEGG pathway terms?

biomart kegg pathways R • 9.3k views
ADD COMMENT
4
Entering edit mode
10.4 years ago
Martin Morgan ★ 1.6k

The REST service is also available from the KEGGREST R / Bioconductor package.

library(KEGGREST)
keggLink("pathway", "hsa:10327")  # pathways gene 10327 is in
p = keggLink("hsa", "pathway")
map = split(names(p), unname(p))  # gene -> pathway map

Queries follow closely the API provided by KEGG.

ADD COMMENT
0
Entering edit mode

Hi,

I tried this method to link gene to pathways. But there was an error: "arguments imply differing number of rows". Tried unsplit also. Can anyone tell, how to solve this error and view the 'map' file.

ADD REPLY
1
Entering edit mode
10.4 years ago
Neilfws 49k

I don't believe that KEGG is accessible through BioMart. However, you can use the TogoWS REST service to retrieve that information. Look at the examples in the KEGG GENES section on that page.

ADD COMMENT
1
Entering edit mode
4.8 years ago
Asaf 10k

I know it's an old question but wanted to update that kegg_enzyme was added to the list of attributes in bioMart, it will return the pathway (the number portion of mapXXXXX pathway) and the enzyme EC number.

try(ensembl <-
            biomaRt::useMart("ensembl", dataset = "mmusculus_gene_ensembl"))
    try(annot <-
          biomaRt::getBM(
            c(
              "ensembl_gene_id",
              "mgi_symbol",
              "chromosome_name",
              "strand",
              "start_position",
              "end_position",
              "gene_biotype",
              "kegg_enzyme"
            ),
            mart = ensembl
          ))

Will return (for instance):

ensembl_gene_id mgi_symbol chromosome_name strand start_position end_position   gene_biotype               kegg_enzyme

1048 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            04714+2.7.11.1
1049 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            05165+2.7.11.1
1050 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            04150+2.7.11.1
1051 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            04151+2.7.11.1
1052 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            05170+2.7.11.1
1053 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            04926+2.7.11.1
1054 ENSMUSG00000030409       Dmpk               7      1       19083849     19093821 protein_coding            05163+2.7.11.1
ADD COMMENT
0
Entering edit mode

This solution is quite useless, it will give all the pathways related to the enzyme so you'll end up getting bacterial pathways related to mouse genes. Splendid

ADD REPLY
0
Entering edit mode
9.9 years ago

You can use 'graphite' R-package.

ADD COMMENT

Login before adding your answer.

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