annotate variant_effect_predictor/Bio/EnsEMBL/Variation/ExternalFeatureVariation.pm @ 3:d30fa12e4cc5 default tip

Merge heads 2:a5976b2dce6f and 1:09613ce8151e which were created as a result of a recently fixed bug.
author devteam <devteam@galaxyproject.org>
date Mon, 13 Jan 2014 10:38:30 -0500
parents 1f6dce3d34e0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
1 =head1 LICENSE
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
2
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
3 Copyright (c) 1999-2012 The European Bioinformatics Institute and
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
4 Genome Research Limited. All rights reserved.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
5
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
6 This software is distributed under a modified Apache license.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
7 For license details, please see
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
8
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
9 http://www.ensembl.org/info/about/code_licence.html
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
10
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
11 =head1 CONTACT
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
12
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
13 Please email comments or questions to the public Ensembl
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
14 developers list at <dev@ensembl.org>.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
15
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
16 Questions may also be sent to the Ensembl help desk at
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
17 <helpdesk@ensembl.org>.
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
18
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
19 =cut
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
20
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
21 package Bio::EnsEMBL::Variation::ExternalFeatureVariation;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
22
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
23 use strict;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
24 use warnings;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
25
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
26 use Bio::EnsEMBL::Variation::ExternalFeatureVariationAllele;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
27 use Bio::EnsEMBL::Variation::Utils::VariationEffect qw(overlap);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
28
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
29 use base qw(Bio::EnsEMBL::Variation::RegulationVariation);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
30
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
31 sub new {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
32 my $class = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
33
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
34 my %args = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
35
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
36 # swap a '-external_feature' argument for a '-feature' one for the superclass
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
37
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
38 for my $arg (keys %args) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
39 if (lc($arg) eq '-external_feature') {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
40 $args{'-feature'} = delete $args{$arg};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
41 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
42 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
43
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
44 # call the superclass constructor
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
45 my $self = $class->SUPER::new(%args) || return undef;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
46
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
47 # rebless the alleles from vfoas to efvas
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
48 map { bless $_, 'Bio::EnsEMBL::Variation::ExternalFeatureVariationAllele' }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
49 @{ $self->get_all_ExternalFeatureVariationAlleles };
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
50
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
51 return $self;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
52 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
53
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
54 sub external_feature_stable_id {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
55 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
56 return $self->SUPER::feature_stable_id(@_);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
57 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
58
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
59 sub external_feature {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
60 my ($self, $ef) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
61 return $self->SUPER::feature($ef, 'ExternalFeature');
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
62 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
63
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
64 sub add_ExternalFeatureVariationAllele {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
65 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
66 return $self->SUPER::add_VariationFeatureOverlapAllele(@_);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
67 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
68
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
69 sub get_reference_ExternalFeatureVariationAllele {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
70 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
71 return $self->SUPER::get_reference_VariationFeatureOverlapAllele(@_);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
72 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
73
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
74 sub get_all_alternate_ExternalFeatureVariationAlleles {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
75 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
76 return $self->SUPER::get_all_alternate_VariationFeatureOverlapAlleles(@_);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
77 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
78
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
79 sub get_all_ExternalFeatureVariationAlleles {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
80 my $self = shift;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
81 return $self->SUPER::get_all_VariationFeatureOverlapAlleles(@_);
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
82 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
83
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
84 sub target_feature_stable_id {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
85 my ($self, $target_feature_stable_id) = @_;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
86
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
87 $self->{target_feature_stable_id} = $target_feature_stable_id if $target_feature_stable_id;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
88
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
89 unless ($self->{target_feature_stable_id}) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
90
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
91 # try to fetch it from the funcgen database
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
92
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
93 if (my $species = $self->{adaptor}->db->species) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
94 for my $entry (
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
95 @{ $self->external_feature->get_all_DBEntries($species.'_core_Transcript') },
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
96 @{ $self->external_feature->get_all_DBEntries($species.'_core_Gene') } ) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
97 if (my $id = $entry->primary_id) {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
98
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
99 $self->{target_feature_stable_id} = $id;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
100
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
101 # there should never be more than one, so we last out of the loop
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
102
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
103 last;
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
104 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
105 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
106 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
107 else {
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
108 warn "Failed to get species from adaptor";
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
109 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
110 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
111
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
112 return $self->{target_feature_stable_id};
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
113 }
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
114
1f6dce3d34e0 Uploaded
mahtabm
parents:
diff changeset
115 1;