Question: How to download an unspliced transcript for a given Ensembl ID using the perl API
0
gravatar for Joseph Hughes
17 months ago by
Joseph Hughes2.7k
Scotland, UK
Joseph Hughes2.7k wrote:

This is the code snippet I have so far:

my $gene_adaptor = $registry->get_adaptor( 'Human', 'Core', 'Gene' );

# Get the 'ensemblID' gene from human
my $gene = $gene_adaptor->fetch_by_gene_stable_id( 'ENSG00000099889');
print "GENE ", $gene->stable_id(), "\n";
print_DBEntries( $gene->get_all_DBEntries() );

foreach my $transcript ( @{ $gene->get_all_Transcripts() } ) {
  print "cDNA: ", $transcript->spliced_seq(), "\n";}

But how do you get the unspliced transcript (i.e. the Exons+Introns without 5' and 3' UTR)?

transcript unspliced ensembl • 463 views
ADD COMMENTlink modified 17 months ago by Emily_Ensembl18k • written 17 months ago by Joseph Hughes2.7k
4
gravatar for Emily_Ensembl
17 months ago by
Emily_Ensembl18k
EMBL-EBI
Emily_Ensembl18k wrote:

There's no direct method for this. I would get the transcript slice, then the sequence of the slice.

print "unspliced: ", $transcript->slice->seq, "/n";
ADD COMMENTlink written 17 months ago by Emily_Ensembl18k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 986 users visited in the last hour