Pull out gene ontology annotation from Ensembl

[owen05nju]

Hi,

I was working on some scripts that uses Ensembl perl API to pull out newest gene build and associated annotations. However, I went into trouble when trying to retrieve gene ontology annotations associated to genes. I know you can easily have this done using BioMart in Ensembl (external references in attributes), but seems tricky to me to achieve this by perl API. Anyone got any ideas? Thanks in advance.

Bin

[flxlex]

Here is a quick copy and paste from an old script that may help or at least get you started:

#  GET THE GENE ADAPTOR FOR THIS DATABASE
my $ga = $dba->get_GeneAdaptor() ; 

#  GET THE ONTOLOGY ADAPTOR FOR THIS DATABASE
#my $go = $dba->get_OntologyAdaptor() ; 

#  GET THE TRANSLATION ADAPTOR FOR THIS DATABASE
my $tra = $dba->get_TranslationAdaptor() ; 

# Get all genes
my @genes = @{$ga->fetch_all_by_biotype('protein_coding') };

# Loop through genes
my %seen;	# prevent duplicates

foreach my $gene (@genes) {
	my $header = ">".$gene->stable_id();	# ensembl ID

	my @db_links=( $gene->get_all_DBLinks() );
    my $db_entries = shift @db_links;
    foreach my $dbe ( @{$db_entries} ) {
    	my $out="$header\t".$dbe->display_id()." ".$dbe->description."\n" if $dbe->dbname() eq "GO";
    	if ($out){
	    	$seen{$out}=1;
	        print $out;
        }
    }
	
}