[ensembl-dev] Wrongly documented gene_adaptor method ?
PATERSON Trevor
trevor.paterson at roslin.ed.ac.uk
Mon Jun 25 13:50:04 BST 2012
According to the Perl documentation, you should be able to get a list of genes for a given display label, using:
public Listref Bio::EnsEMBL::DBSQL::GeneAdaptor::fetch_all_by_display_label ( )
However, this method is wrongly documented in the online perl docs, and does not return an array - but a single gene object
my @genes = @{$geneAdaptor->fetch_by_display_label("TBC1D3B")};
DEBUGGER ERROR MESSAGE : Not an ARRAY reference
When you look at the real code in Release 67 API, rather than the www perldocs, the method actually returns a single gene, preferring one that is on a reference slice.
I am looking for a way of getting alternative gene models on Patch exceptions: and this method would have been promising. e.g for the gene ENSG00000224226<http://www.ensembl.org/Homo_sapiens/Gene/Summary?db=core;g=ENSG00000224226;r=HG75_PATCH:34442621-34976908> Location Chromosome 17: 34,493,061-34,503,984<http://www.ensembl.org/Homo_sapiens/Location/View?db=core;g=ENSG00000224226;r=17:34493061-34503984> how do I get its equivalent on the patch-fix ENSG00000260069<http://www.ensembl.org/Homo_sapiens/Gene/Summary?db=core;g=ENSG00000260069;r=HG75_PATCH:34442621-34976908> Location Chromosome HG75_PATCH: 34,493,132-34,503,984<http://www.ensembl.org/Homo_sapiens/Location/View?db=core;g=ENSG00000260069;r=HG75_PATCH:34493132-34503984>.
obviously I can dig into the API and call the appropriate sql code myself:
my $constraint = "x.display_label = '". $gene->display_xref()->display_id()."' AND g.is_current = 1";
my @same_labelled_genes = @{ $gene_adaptor->generic_fetch($constraint) };
but is there an alternative (simple) method to get either 'same_labelled' genes or ideally 'equivalent-exception' genes for a given gene?
the method
public Listref Bio::EnsEMBL::Gene::get_all_alt_alleles ( )
of course does something similar using the alt_alleles table , but only for alternative haplotypes not patches.
Thanks
Trevor
Trevor Paterson PhD
trevor.paterson at roslin.ed.ac.uk<mailto:trevor.paterson at roslin.ed.ac.uk>
Bioinformatics
The Roslin Institute
Royal (Dick) School of Veterinary Studies
University of Edinburgh
Easter Bush
Midlothian
EH25 9RG
Scotland UK
phone +44 (0)131 651 9157
http://bioinformatics.roslin.ed.ac.uk/
Please consider the environment before printing this e-mail
The University of Edinburgh is a charitable body, registered in Scotland with registration number SC005336
Disclaimer:This e-mail and any attachments are confidential and intended solely for the use of the recipient(s) to whom they are addressed. If you have received it in error, please destroy all copies and inform the sender.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20120625/fe287356/attachment.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20120625/fe287356/attachment.ksh>
More information about the Dev
mailing list