Hi Graham,<br><br>Thanks for your fast answer. This approach however doesn't seem to work for me.. I tried both with your given example of a SNP and some SNPs that were annotated as "REGULATORY_REGION" in previous versions of Ensembl e.g. rs13076120. <br>
<br>When I run the code I don't get any error messages, but I don't get any
annotation to regulatory regions either. I tried with all the 3
methods, but in all cases the result was the same. Do you know what
could be a possible reason for that?<br><br>I have tried it both with a local installation of a database, and connecting to your server. I also can't see a way to extract this information using Biomart.<br><br><div class="gmail_quote">
<div>Best regards,<br><br>Agata<br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><font color="#888888"></font><div><div class="h5">
<div class="gmail_quote">On Thu, May 12, 2011 at 5:07 PM, Graham Ritchie <span dir="ltr"><<a href="mailto:grsr@ebi.ac.uk" target="_blank">grsr@ebi.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

Hi Agata,<br>
<br>
We changed the consequence type calculation code for release 62 and the way in which we integrate with the Regulation databases has changed significantly. We used to only include consequences on regulatory regions that were associated with transcripts somehow, but now you can find out if your variants lie in any of the various types of feature stored in the regulation database. There are 3 new methods in the VariationFeature class that let you find out if your variant overlaps a feature from the regulation database: get_all_RegulatoryFeatureVariations, get_all_ExternalFeatureVariations and get_all_MotifFeatureVariations, the objects returned from these methods have a similar interface as for TranscriptVariations, so you can print out the consequence term with, for example, the display_consequence method. Please refer to the API documentation for more details, but here's a little bit of example code:<br>


<br>
my $vdba = $reg->get_DBAdaptor('human', 'variation');<br>
<br>
my $vfa = $vdba->get_VariationFeatureAdaptor;<br>
my $va = $vdba->get_VariationAdaptor;<br>
<br>
my $v = $va->fetch_by_name('rs10054597');<br>
<br>
my $vf = $vfa->fetch_all_by_Variation($v)->[0];<br>
<br>
for my $rfv (@{ $vf->get_all_RegulatoryFeatureVariations }) {<br>
    print $rfv->display_consequence, "\n";<br>
}<br>
<br>
Regarding your second question, the regulation team store miRNA target sites from miRanda in the regulation database as ExternalFeatures, so you can check if your variant lies in one of these by calling get_all_ExternalFeatureVariations and if any are returned checking if the associated ExternalFeatureVariation->feature->feature_set->name is 'miRanda miRNA targets'.<br>


<br>
Cheers,<br>
<font color="#888888"><br>
Graham<br>
</font><div><div></div><div><br>
<br>
On 12 May 2011, at 15:02, Agata Wesolowska wrote:<br>
<br>
> Hi,<br>
><br>
> I was extracting known SNPs together with their consequences on transcripts from a list of genes using Ensembl API. I noticed that the consequence "REGULATORY_REGION" is not present anymore among the possible transcript consequences. Is there any reason for that? And is there any way to obtain this information?<br>


><br>
> Also another question - is it possible to find SNPs in microRNA target sites using API?<br>
><br>
> Best regards,<br>
><br>
> Agata Wesolowska<br>
><br>
><br>
><br>
><br>
> _______________________________________________<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/mailman/listinfo/dev</a><br>
> Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
<br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>