[ensembl-dev] retrieving refseq gene build has no external name?
Stef van Lieshout
mailsvl at fastmail.fm
Mon Jul 29 13:46:03 BST 2013
Nobody knows a way to do this? (retrieving refseq gene models by
gene-name)
Stef
----- Original message -----
From: Stef van Lieshout <mailsvl at fastmail.fm>
To: dev at ensembl.org
Subject: [ensembl-dev] retrieving refseq gene build has no external
name?
Date: Wed, 10 Jul 2013 13:45:01 +0200
Hi all,
It seems the refseq_human_imort gene models do not contain an external
name, is that correct?
Goal:
What I try to do (with the API) is to get all refseq transcripts (NM_*)
instead of ensembl transcripts (ENST*) for a certain gene name (eg
BRCA2).
You can do this with the "core" set by:
my $ga = $registry->get_adaptor( 'human', 'core', 'gene' );
my @genes = @{$ga->fetch_all_by_external_name( 'BRCA2' )};
And I hoped to be able to do this for refseq models with:
my $ga = $registry->get_adaptor( 'human', 'otherfeatures', 'gene' );
my @genes = @{$ga->fetch_all_by_external_name( 'BRCA2' )};
But that seems to fail because there are no external names defined (the
key is there, but set to undef) for these genes.
Is there another way to do this?
Btw, I know you can get the transcripts by using the NM_ id as input,
but this is not what I want:
my $gene = $ga->fetch_by_transcript_stable_id( 'NM_000059.3' );
/Stef
_______________________________________________
Dev mailing list Dev at ensembl.org
Posting guidelines and subscribe/unsubscribe info:
http://lists.ensembl.org/mailman/listinfo/dev
Ensembl Blog: http://www.ensembl.info/
More information about the Dev
mailing list