<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Jerry,<div><br></div><div>Up to now, our main effort with RNAseq has been focussed around gene annotation. We have developed a method where we can annotate protein-coding gene models from RNAseq data. This method includes a step where we annotate introns and count the number of reads that span each intron. The intron features are stored as dna_align_features and their score column contains a count of the number of reads that support the intron. </div><div><br></div><div>The counts are also visible on our website. In the example below, we have the human MYOG gene. If you click the grey boxes for Breast intron-spanning reads and Skeletal muscle intron-spanning reads, you can see in the Read Coverage line that more reads support the skeletal muscle (51) than breast tissue (9).</div><div><br></div><div><img height="213" width="1008" apple-width="yes" apple-height="yes" id="2e542592-e38e-455b-9df7-9da20694cc4a" src="cid:C90A19FB-EC54-4C34-9283-CDEFDF3B7053@globalsuite.net"></div><div><br></div><div>In the future we would like to provide users with an expression level but we have not yet looked into this yet.</div><div><br></div><div>Cheers,</div><div>Bronwen</div><div><br></div><div><br><div><div>On 31 Aug 2012, at 02:15, Jerry Li wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Thibaut and Ensembl dev community,</div><div> My sincere gratitute for your help. I took your advice and now I'm able to map RNAseq gene ids to EntrezGene ids. </div><div> </div><div>I need to ask one more big favor from you. Has any of you tried retrieving mRNA abundance information on a per Exon or per Transcript basis from RNAseq? I have read through the Ensembl API Core documentation and tutorial but still have no clue on how to do it. </div>
<div> </div><div>Your help is greatly appreciated!!!</div><div> </div><div>Jerry</div><div><br> </div><div class="gmail_quote">On Thu, Aug 30, 2012 at 3:50 AM, th3 <span dir="ltr"><<a href="mailto:th3@sanger.ac.uk" target="_blank">th3@sanger.ac.uk</a>></span> wrote:<br>
<blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid" class="gmail_quote">Hi Jerry,<br>
<br>
We only get external reference (xref) like the EntrezGene ID for the core databases.<br>
So for every gene from the rnaseq database you need to get the gene from core that overlaps your rnaseq gene and then you can check if they have an xref of your choice.<br>
<br>
# You will need a GeneAdaptor on your core DB<br>
$coredbga = $coredb->get_GeneAdaptor;<br>
# You will need a SliceAdaptor on your core DB<br>
$coredbsa = $coredb->get_SliceAdaptor;<div class="im"><br>
foreach my $slice (@$slices) {<br>
# fetch genes from rnaseq<br>
my $genes = $slice->get_all_Genes($logic_<u></u>name);<br>
foreach my $gene (@$genes) {<br></div>
$coreslice = $coresa->fetch_by_name($gene-><u></u>feature_Slice()->name);<br>
foreach my $cgene (@{$coredbga->fetch_all_by_<u></u>Slice_and_external_dbname_<u></u>link($coreslice, undef, undef, "EntrezGene")}) {<br>
foreach my $dbes ($cgene->get_all_DBLinks("<u></u>EntrezGene")){<div class="im"><br>
foreach my $dbe (@{$dbes}){<br>
print $dbe->primary_id(), "\t";<br>
}<br>
}<br>
}<br>
}<br>
}<br>
<br></div>
Hope this helps<br>
<br>
Thibaut<div class="im"><br>
<br>
On 2012-08-30 04:24, Jerry Li wrote:<br>
</div><blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid" class="gmail_quote"><div class="im">
Hi All,<br>
Sorry to bother you all on this topic but I can really use some<br>
advice.<br>
<br>
Based on an example I found in the email archive<br>
</div><a href="http://lists.ensembl.org/pipermail/dev/2010-September/005130.html" target="_blank">http://lists.ensembl.org/<u></u>pipermail/dev/2010-September/<u></u>005130.html</a> [1]<div class="im"><br>
I put together the following code to retrieve the EntrezGene IDs based<br>
on Ensembl $gene. However, no result was returned. <br>
The database I accessed was homo_sapiens_rnaseq_68_37, and logic_name<br>
is ovary_rnaseq.<br>
Can anyone provide some pointer on how I can get this data?<br>
<br>
foreach my $slice (@$slices) {<br>
# fetch genes from rnaseq<br>
my $genes = $slice->get_all_Genes($logic_<u></u>name);<br>
foreach my $gene (@$genes) {<br>
foreach my $dbes ($gene->get_all_DBLinks("<u></u>EntrezGene")){<br>
foreach my $dbe (@{$dbes}){<br>
print $dbe->primary_id(), "t";<br>
}<br>
}<br>
}<br>
}<br>
<br>
<br>
Thanks,<br>
Jerry<br>
<br>
<br></div>
Links:<br>
------<br>
[1] <a href="http://lists.ensembl.org/pipermail/dev/2010-September/005130.html" target="_blank">http://lists.ensembl.org/<u></u>pipermail/dev/2010-September/<u></u>005130.html</a><br>
</blockquote>
<br>
<br>
______________________________<u></u>_________________<br>
Dev mailing list <a href="mailto:Dev@ensembl.org" target="_blank">Dev@ensembl.org</a><br>
List admin (including subscribe/unsubscribe): <a href="http://lists.ensembl.org/mailman/listinfo/dev" target="_blank">http://lists.ensembl.org/<u></u>mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
</blockquote></div><br>
_______________________________________________<br>Dev mailing list <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>List admin (including subscribe/unsubscribe): <a href="http://lists.ensembl.org/mailman/listinfo/dev">http://lists.ensembl.org/mailman/listinfo/dev</a><br>Ensembl Blog: <a href="http://www.ensembl.info/">http://www.ensembl.info/</a><br></blockquote></div><br></div></body></html>