[ensembl-dev] retrieving refseq gene build has no external name?

Stef van Lieshout mailsvl at fastmail.fm
Wed Jul 10 12:45:01 BST 2013


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




More information about the Dev mailing list