How Many Human Pathways Exist In Kegg Database?
6
7
Entering edit mode
13.4 years ago
Woa ★ 2.9k

How to know how many Human pathways exist in the KEGG database?

If I search THIS PAGE with 'hsa'(Homo sapiens) in the organism field and keeping the keyword field blank it returns 236. Is this the correct way to find it?

kegg pathway database • 8.3k views
ADD COMMENT
7
Entering edit mode
13.4 years ago
D-Horse ▴ 110

Hi, I think you could use the APIs in KEGG.For example to get the total numbers of genes in hsa00010.The code is written in Perl:

#!/usr/bin/env perl
use SOAP::Lite +trace => [qw(debug)];
my $serv = SOAP::Lite -> service("http://soap.genome.jp/KEGG.wsdl");
my $result = $serv->get_genes_by_pathway('path:hsa00010');
my $count = 0;
foreach (@{$result}){
    print $_,"\n";
    $count += 1;
}
print $count,"\n";

The answer is 65.

ADD COMMENT
0
Entering edit mode

Many Thanks !!!!

ADD REPLY
0
Entering edit mode

+1 for KEGG API, that's a neat resource !

ADD REPLY
6
Entering edit mode
10.6 years ago
Prakki Rama ★ 2.7k

KEGG now operates with REST instead of SOAP. Could be done in a feasible way. But I did this:

wget -O http://rest.kegg.jp/link/pathway/hsa >HSA.txt ##Saving the file
cut -d "      " -f 2 HSA.txt | cut -d ":" -f 2 | sort -u | wc -l ## counting the pathways in Human
283 (pathways in Human in KEGG at present)

~Prakki Datta

ADD COMMENT
3
Entering edit mode
13.4 years ago
Joachim ★ 2.9k

Yes, this is the correct number of pathways for human in KEGG and it is also the right way to find them.

ADD COMMENT
0
Entering edit mode

Many thanks. I've a further question: for a particular pathway like http://www.genome.jp/dbget-bin/www_bget?pathway+hsa00010, I wish to get the total no. of genes involved in that pathway. Is there an easy way to find it? Thanks in advance.

ADD REPLY
0
Entering edit mode

I guess D-Horse was faster with providing an answer.. :)

ADD REPLY
1
Entering edit mode
12.0 years ago
wchknudt ▴ 50

You are correct. The answer now is 260.

You can also use perl script to do this:

#!/usr/bin/env perl

use SOAP::Lite;
$wsdl = 'http://soap.genome.jp/KEGG.wsdl';

$results = SOAP::Lite
-> service($wsdl)
-> list_pathways("hsa");

my $file_name = "kegg_hsa_pathways.txt";

if (open (OUTFILE, '>'.$file_name))
{
    foreach $path (@{$results}) {
        print OUTFILE "$path->{entry_id}\t$path->{definition}\n";
    }
    close(OUTFILE);

    print "Content writtent to $file_name!"
}
else{
    die "Error: Error occur in  writing to file\n";
}
ADD COMMENT
0
Entering edit mode
10.6 years ago
Haluk ▴ 190

If you need their graphs, I generated them in dot format.

https://github.com/haluk/kegg

ADD COMMENT

Login before adding your answer.

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