Inheritance diagram for Bio::EnsEMBL::DBSQL::PredictionExonAdaptor:

Public Member Functions
protected	_tables ()

protected	_columns ()

protected	_final_clause ()

public Listref	fetch_all_by_PredictionTranscript ()

public void	store ()

public void	remove ()

public List	list_dbIDs ()

protected	_objs_from_sth ()

Public Member Functions inherited from Bio::EnsEMBL::DBSQL::BaseFeatureAdaptor
public Bio::EnsEMBL::BaseFeatureAdaptor	new ()

public Boolean	start_equals_end ()

public void	clear_cache ()

protected	_clear_slice_feature_cache ()

protected Bio::EnsEMBL::Utils::Cache	_slice_feature_cache ()

public Listref	fetch_all_by_Slice ()

public Bio::EnsEMBL::Utils::Iterator	fetch_Iterator_by_Slice_method ()

public Bio::EnsEMBL::Utils::Iterator	fetch_Iterator_by_Slice ()

public Listref	fetch_all_by_Slice_and_score ()

public Listref	fetch_all_by_Slice_constraint ()

public Listref	fetch_all_by_logic_name ()

public Listref	fetch_all_by_stable_id_list ()

protected	_create_feature ()

protected	_create_feature_fast ()

public Int	count_by_Slice_constraint ()

protected ArrayRef	_get_and_filter_Slice_projections ()

protected ArrayRef	_generate_feature_bounds ()

protected Listref	_get_by_Slice ()

protected	_slice_fetch ()

public	get_seq_region_id_external ()

public	get_seq_region_id_internal ()

protected	_pre_store ()

protected	_pre_store_userdata ()

protected	_check_start_end_strand ()

protected	_remap ()

protected	_seq_region_boundary_from_db ()

public void	store ()

public void	remove ()

public void	remove_by_Slice ()

protected	_max_feature_length ()

protected	_list_seq_region_ids ()

public	remove_by_analysis_id ()

public	remove_by_feature_id ()

public Bio::EnsEMBL::Feature	fetch_nearest_by_Feature ()

public Listref	fetch_all_by_outward_search ()

public Listref	fetch_all_nearest_by_Feature ()

public Listref	select_nearest ()

protected Int	_compute_nearest_end ()

protected Int	_compute_prime_distance ()

protected Int	_compute_midpoint ()

protected	_discard_excess_features_from_matrix ()

Public Member Functions inherited from Bio::EnsEMBL::DBSQL::BaseAdaptor
public Bio::EnsEMBL::DBSQL::BaseAdaptor	new ()

public DBI::StatementHandle	prepare ()

public Bio::EnsEMBL::DBSQL::DBAdaptor	db ()

public Bio::EnsEMBL::DBSQL::DBConnection	dbc ()

public Boolean	is_multispecies ()

public Int	species_id ()

protected	_list_dbIDs ()

protected	_straight_join ()

protected	_can_straight_join ()

public Listref	bind_param_generic_fetch ()

protected	_bind_param_generic_fetch ()

public String	generate_in_constraint ()

public Listref	generic_fetch ()

public Int	generic_count ()

protected	_generate_sql ()

public Bio::EnsEMBL::Feature	fetch_by_dbID ()

protected	_uncached_fetch_by_dbID ()

public Listref	fetch_all_by_dbID_list ()

protected	_uncached_fetch_all_by_dbID_list ()

protected ArrayRef	_uncached_fetch_all_by_id_list ()

public	fetch_all ()

public Scalar	last_insert_id ()

public	insert_ignore_clause ()

protected Bio::EnsEMBL::DBSQL::Support::BaseCache	_id_cache ()

protected Boolean	_no_id_cache ()

public Boolean	ignore_cache_override ()

public Int	schema_version ()

protected	_tables ()

protected	_columns ()

protected	_default_where_clause ()

protected	_left_join ()

protected	_final_clause ()

protected	_objs_from_sth ()

protected	_build_id_cache ()

protected	_logic_name_to_constraint ()

Detailed Description

Synopsis

$pea   = $database_adaptor->get_PredictionExonAdaptor();
$pexon = $pea->fetch_by_dbID();
 
my $slice =
  $database_adaptor->get_SliceAdaptor->fetch_by_region( 'X', 1, 1e6 );
 
my @pexons = @{ $pea->fetch_all_by_Slice($slice) };

Definition at line 19 of file PredictionExonAdaptor.pm.

Member Function Documentation

◆ _columns()

protected Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::_columns ( )

Undocumented method

Code:

click to view

◆ _final_clause()

protected Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::_final_clause ( )

Undocumented method

Code:

click to view

◆ _objs_from_sth()

protected Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::_objs_from_sth ( )

Undocumented method

Code:

click to view

sub _objs_from_sth {
  my ($self, $sth, $mapper, $dest_slice) = @_;
 
  #
  # This code is ugly because an attempt has been made to remove as many
  # function calls as possible for speed purposes.  Thus many caches and
  # a fair bit of gymnastics is used.
  #
  my $sa = $self->db()->get_SliceAdaptor();
 
  my @exons;
  my %slice_hash;
  my %sr_name_hash;
  my %sr_cs_hash;
 
  my(
    $prediction_exon_id, $seq_region_id,
    $seq_region_start,   $seq_region_end, $seq_region_strand,
    $start_phase,        $score,          $p_value);
 
  $sth->bind_columns(\(
                      $prediction_exon_id, $seq_region_id,
                      $seq_region_start,   $seq_region_end, $seq_region_strand,
                      $start_phase,        $score,          $p_value));
 
  my $dest_slice_start;
  my $dest_slice_end;
  my $dest_slice_strand;
  my $dest_slice_length;
  my $dest_slice_cs;
  my $dest_slice_sr_name;
  my $dest_slice_sr_id;
  my $asma;
 
  if ($dest_slice) {
    $dest_slice_start   = $dest_slice->start();
    $dest_slice_end     = $dest_slice->end();
    $dest_slice_strand  = $dest_slice->strand();
    $dest_slice_length  = $dest_slice->length();
    $dest_slice_cs      = $dest_slice->coord_system();
    $dest_slice_sr_name = $dest_slice->seq_region_name();
    $dest_slice_sr_id   = $dest_slice->get_seq_region_id();
    $asma               = $self->db->get_AssemblyMapperAdaptor();
  }
 
  FEATURE: while($sth->fetch()) {
 
    #need to get the internal_seq_region, if present
    $seq_region_id = $self->get_seq_region_id_internal($seq_region_id);
    my $slice = $slice_hash{"ID:".$seq_region_id};
 
    if (!$slice) {
      $slice                            = $sa->fetch_by_seq_region_id($seq_region_id);
      $slice_hash{"ID:".$seq_region_id} = $slice;
      $sr_name_hash{$seq_region_id}     = $slice->seq_region_name();
      $sr_cs_hash{$seq_region_id}       = $slice->coord_system();
    }
 
    #obtain a mapper if none was defined, but a dest_seq_region was
    if(!$mapper && $dest_slice && !$dest_slice_cs->equals($slice->coord_system)) {
      $mapper = $asma->fetch_by_CoordSystems($dest_slice_cs, $slice->coord_system);
    }
 
    my $sr_name = $sr_name_hash{$seq_region_id};
    my $sr_cs   = $sr_cs_hash{$seq_region_id};
 
    #
    # remap the feature coordinates to another coord system
    # if a mapper was provided
    #
 
    if ($mapper) {
 
      if (defined $dest_slice && $mapper->isa('Bio::EnsEMBL::ChainedAssemblyMapper') ) {
        ($seq_region_id, $seq_region_start, $seq_region_end, $seq_region_strand) =
         $mapper->map($sr_name, $seq_region_start, $seq_region_end, $seq_region_strand, $sr_cs, 1, $dest_slice);
 
      } else {
        ($seq_region_id, $seq_region_start, $seq_region_end, $seq_region_strand) =
         $mapper->fastmap($sr_name, $seq_region_start, $seq_region_end, $seq_region_strand, $sr_cs);
      }
 
      #skip features that map to gaps or coord system boundaries
      next FEATURE if (!defined($seq_region_id));
 
      #get a slice in the coord system we just mapped to
      $slice = $slice_hash{"ID:".$seq_region_id} ||= $sa->fetch_by_seq_region_id($seq_region_id);
    }
 
    #
    # If a destination slice was provided convert the coords.
    #
    if (defined($dest_slice)) {
      my $seq_region_len = $dest_slice->seq_region_length();
 
      if ( $dest_slice_strand == 1 ) {
        $seq_region_start = $seq_region_start - $dest_slice_start + 1;
        $seq_region_end   = $seq_region_end - $dest_slice_start + 1;
 
        if ( $dest_slice->is_circular ) {
        # Handle circular chromosomes.
 
          if ( $seq_region_start > $seq_region_end ) {
            # Looking at a feature overlapping the chromosome origin.
 
            if ( $seq_region_end > $dest_slice_start ) {
              # Looking at the region in the beginning of the chromosome
              $seq_region_start -= $seq_region_len;
            }
            if ( $seq_region_end < 0 ) {
              $seq_region_end += $seq_region_len;
            }
          } else {
            if ($dest_slice_start > $dest_slice_end && $seq_region_end < 0) {
              # Looking at the region overlapping the chromosome
              # origin and a feature which is at the beginning of the
              # chromosome.
              $seq_region_start += $seq_region_len;
              $seq_region_end   += $seq_region_len;
            }
          }
        }
      } else {
 
        my $start = $dest_slice_end - $seq_region_end + 1;
        my $end = $dest_slice_end - $seq_region_start + 1;
 
        if ($dest_slice->is_circular()) {
 
          if ($dest_slice_start > $dest_slice_end) {
            # slice spans origin or replication
 
            if ($seq_region_start >= $dest_slice_start) {
              $end += $seq_region_len;
              $start += $seq_region_len if $seq_region_end > $dest_slice_start;
 
            } elsif ($seq_region_start <= $dest_slice_end) {
              # do nothing
            } elsif ($seq_region_end >= $dest_slice_start) {
              $start += $seq_region_len;
              $end += $seq_region_len;
 
            } elsif ($seq_region_end <= $dest_slice_end) {
              $end += $seq_region_len if $end < 0;
 
            } elsif ($seq_region_start > $seq_region_end) {
              $end += $seq_region_len;
            }
 
          } else {
 
            if ($seq_region_start <= $dest_slice_end and $seq_region_end >= $dest_slice_start) {
              # do nothing
            } elsif ($seq_region_start > $seq_region_end) {
              if ($seq_region_start <= $dest_slice_end) {
                $start -= $seq_region_len;
              } elsif ($seq_region_end >= $dest_slice_start) {
                $end += $seq_region_len;
              }
            }
          }
        }
 
        $seq_region_start = $start;
        $seq_region_end = $end;
        $seq_region_strand *= -1;
 
      } ## end else [ if ( $dest_slice_strand...)]
 
      # Throw away features off the end of the requested slice or on
      # different seq_region.
      if ($seq_region_end < 1
          || $seq_region_start > $dest_slice_length
          || ($dest_slice_sr_id != $seq_region_id)) {
        next FEATURE;
      }
      $slice = $dest_slice;
    }
 
    # Finally, create the new PredictionExon.
    push( @exons,
          $self->_create_feature( 'Bio::EnsEMBL::PredictionExon', {
                                    '-start'   => $seq_region_start,
                                    '-end'     => $seq_region_end,
                                    '-strand'  => $seq_region_strand,
                                    '-adaptor' => $self,
                                    '-slice'   => $slice,
                                    '-dbID'    => $prediction_exon_id,
                                    '-phase'   => $start_phase,
                                    '-score'   => $score,
                                    '-p_value' => $p_value
                                  } ) );
 
  }
 
  return \@exons;
}

◆ _tables()

protected Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::_tables ( )

Undocumented method

Code:

click to view

◆ fetch_all_by_PredictionTranscript()

public Listref Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::fetch_all_by_PredictionTranscript ( )

  Arg [1]    : Bio::EnsEMBL::PredcitionTranscript $transcript
  Example    :

none

  Description: Retrieves all Exons for the Transcript in 5-3 order
  Returntype : listref Bio::EnsEMBL::Exon on Transcript slice 
  Exceptions : throws if transcript does not have a slice
  Caller     : Transcript->get_all_Exons()
  Status     : Stable

Code:

click to view

◆ list_dbIDs()

public List Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::list_dbIDs ( )

  Arg [1]    : none
  Example    :

@exon_ids = @{$exon_adaptor->list_dbIDs()};

  Description: Gets an array of internal ids for all exons in the current db
  Arg[1]     : <optional> int. not 0 for the ids to be sorted by the seq_region.
  Returntype : list of ints
  Exceptions : none
  Caller     : ?
  Status     : Stable

Code:

click to view

◆ remove()

public void Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::remove ( )

  Arg [1]    : Bio::EnsEMBL::PredictionExon $exon
               the exon to remove from the database 
  Example    :

$exon_adaptor->remove($exon);

  Description: Removes an exon from the database
  Returntype : none
  Exceptions : none
  Caller     : general
  Status     : Stable

Code:

click to view

◆ store()

public void Bio::EnsEMBL::DBSQL::PredictionExonAdaptor::store ( )

  Arg [1]    : Bio::EnsEMBL::PredictionExon $exon
               The exon to store in this database
  Arg [2]    : int $prediction_transcript_id
               The internal identifier of the prediction exon that that this
               exon is associated with.
  Arg [3]    : int $rank
               The rank of the exon in the transcript (starting at 1)
  Example    :

$pexon_adaptor->store($pexon, 1211, 2);

  Description: Stores a PredictionExon in the database
  Returntype : none
  Exceptions : thrown if exon does not have a slice attached
               or if $exon->start, $exon->end, $exon->strand, or $exon->phase 
               are not defined or if $exon is not a Bio::EnsEMBL::PredictionExon 
  Caller     : general
  Status     : Stable

Code:

click to view

The documentation for this class was generated from the following file:

ensembl/modules/Bio/EnsEMBL/DBSQL/PredictionExonAdaptor.pm

Public Member Functions

Detailed Description

Synopsis

Member Function Documentation

◆ _columns()

◆ _final_clause()

◆ _objs_from_sth()

◆ _tables()

◆ fetch_all_by_PredictionTranscript()

◆ list_dbIDs()

◆ remove()

◆ store()