[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