<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class=""><br class="">
</div>
<div class="">
<div class="">
<div class="">Dear Ensembl dev team,</div>
<div class=""><br class="">
</div>
<div class="">I try to get rsquare between two SNPs (on the same chromosome). When I use the ldFeatureContainer function get_r_square, however, I get a warning and no return rsquare value. Here is the function signature from the documentation:</div>
<div class=""><br class="">
</div>
<div class="">
<pre class="fragment" style="font-family: monospace, fixed; font-size: 9pt; margin: 4px 8px 4px 2px; line-height: 15px; border: 1px solid rgb(196, 207, 229); background-color: rgb(251, 252, 253); padding: 4px 6px; overflow: auto; word-wrap: break-word;">$r_square = $obj->get_r_square($vf1,$vf2,$population_id);</pre>
<div class=""><br class="">
</div>
</div>
<div class="">The script I am using is the following:</div>
<div class=""><br class="">
</div>
<div class=""><font face="Courier New" class=""># START SCRIPT</font></div>
<div class=""><font face="Courier New" class="">use strict;</font></div>
<div class=""><font face="Courier New" class="">use warnings;</font></div>
<div class=""><font face="Courier New" class="">use Bio::EnsEMBL::Registry;</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">my $start_run = time();</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">my $registry = 'Bio::EnsEMBL::Registry';</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">$registry->load_registry_from_db(</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>-host => '<a href="http://ensembldb.ensembl.org" class="">ensembldb.ensembl.org</a>',</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>-user => 'anonymous'</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>);</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class=""># Connect to the databases:</font></div>
<div class=""><font face="Courier New" class="">my $variation_adaptor = $registry->get_adaptor('homo_sapiens', 'variation', 'variation');</font></div>
<div class=""><font face="Courier New" class="">my $ldfc_adaptor = $registry->get_adaptor('homo_sapiens', 'variation', 'ldfeaturecontainer');</font></div>
<div class=""><font face="Courier New" class="">my $pop_adaptor = $registry->get_adaptor('homo_sapiens', 'variation', 'population');</font></div>
<div class=""><font face="Courier New" class="">$variation_adaptor->db->use_vcf(1);</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">my $ld_population = $pop_adaptor->fetch_by_name('1000GENOMES:phase_3:CEU');</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class=""># Build tagSNP list of SNPs on same chromosome:</font></div>
<div class=""><font face="Courier New" class="">my $tagSNP1 = 'rs6678275';</font></div>
<div class=""><font face="Courier New" class="">my $tagSNP2 = 'rs6656401';</font></div>
<div class=""><font face="Courier New" class="">my @variations;</font></div>
<div class=""><font face="Courier New" class="">push (@variations, $variation_adaptor->fetch_by_name($tagSNP1));</font></div>
<div class=""><font face="Courier New" class="">push (@variations, $variation_adaptor->fetch_by_name($tagSNP2));</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">my @all_features;</font></div>
<div class=""><font face="Courier New" class="">foreach my $variation (@variations) {</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my @var_features = @{ $variation->get_all_VariationFeatures() };</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>foreach my $vf (@var_features) {</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>push (@all_features, $vf);</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>}</font></div>
<div class=""><font face="Courier New" class="">}</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">my $ldfc = $ldfc_adaptor->fetch_by_VariationFeatures(\@all_features, $ld_population);</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class=""># Print r2 between all SNPs:</font></div>
<div class=""><font face="Courier New" class="">foreach my $vf (@all_features) {</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>foreach my $vf2 (@all_features) {</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $rsid = $vf->name;</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $rsid2 = $vf2->name;</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $start = $vf->start;</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $start2 = $vf2->start;</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $region = $vf->seq_region_name;</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $region2 = $vf2->seq_region_name;</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>print "$region\t$rsid\t$start\n";</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>print "$region2\t$rsid2\t$start2\n";</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>my $r2 = $ldfc->get_r_square($vf, $vf2);</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>#print "\tr2=$r2\n";</font></div>
<div class=""><font face="Courier New" class=""><span class="Apple-tab-span" style="white-space: pre;"></span>}</font></div>
<div class=""><font face="Courier New" class="">}</font></div>
</div>
<div class=""><font face="Courier New" class=""># END SCRIPT</font></div>
<div class=""><br class="">
</div>
<div class="">I get this printout:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><font face="Courier New" class="">$ perl rsquareTagSNPs.pl out_tagsnps.txt 0.9</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6678275<span class="Apple-tab-span" style="white-space:pre"> </span>193656103</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6678275<span class="Apple-tab-span" style="white-space:pre"> </span>193656103</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">-------------------- WARNING ----------------------</font></div>
<div class=""><font face="Courier New" class="">MSG: variation features have no pairwise ld information</font></div>
<div class=""><font face="Courier New" class="">FILE: EnsEMBL/Variation/LDFeatureContainer.pm LINE: 207</font></div>
<div class=""><font face="Courier New" class="">CALLED BY: rsquareTagSNPs.pl LINE: 59</font></div>
<div class=""><font face="Courier New" class="">Date (localtime) = Mon Feb 29 09:47:44 2016</font></div>
<div class=""><font face="Courier New" class="">Ensembl API version = 83</font></div>
<div class=""><font face="Courier New" class="">---------------------------------------------------</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6678275<span class="Apple-tab-span" style="white-space:pre"> </span>193656103</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6656401<span class="Apple-tab-span" style="white-space:pre"> </span>207518704</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">-------------------- WARNING ----------------------</font></div>
<div class=""><font face="Courier New" class="">MSG: variation features have no pairwise ld information</font></div>
<div class=""><font face="Courier New" class="">FILE: EnsEMBL/Variation/LDFeatureContainer.pm LINE: 207</font></div>
<div class=""><font face="Courier New" class="">CALLED BY: rsquareTagSNPs.pl LINE: 59</font></div>
<div class=""><font face="Courier New" class="">Date (localtime) = Mon Feb 29 09:47:44 2016</font></div>
<div class=""><font face="Courier New" class="">Ensembl API version = 83</font></div>
<div class=""><font face="Courier New" class="">---------------------------------------------------</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6656401<span class="Apple-tab-span" style="white-space:pre"> </span>207518704</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6678275<span class="Apple-tab-span" style="white-space:pre"> </span>193656103</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">-------------------- WARNING ----------------------</font></div>
<div class=""><font face="Courier New" class="">MSG: variation features have no pairwise ld information</font></div>
<div class=""><font face="Courier New" class="">FILE: EnsEMBL/Variation/LDFeatureContainer.pm LINE: 207</font></div>
<div class=""><font face="Courier New" class="">CALLED BY: rsquareTagSNPs.pl LINE: 59</font></div>
<div class=""><font face="Courier New" class="">Date (localtime) = Mon Feb 29 09:47:44 2016</font></div>
<div class=""><font face="Courier New" class="">Ensembl API version = 83</font></div>
<div class=""><font face="Courier New" class="">---------------------------------------------------</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6656401<span class="Apple-tab-span" style="white-space:pre"> </span>207518704</font></div>
<div class=""><font face="Courier New" class="">1<span class="Apple-tab-span" style="white-space:pre">
</span>rs6656401<span class="Apple-tab-span" style="white-space:pre"> </span>207518704</font></div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class=""><font face="Courier New" class="">-------------------- WARNING ----------------------</font></div>
<div class=""><font face="Courier New" class="">MSG: variation features have no pairwise ld information</font></div>
<div class=""><font face="Courier New" class="">FILE: EnsEMBL/Variation/LDFeatureContainer.pm LINE: 207</font></div>
<div class=""><font face="Courier New" class="">CALLED BY: rsquareTagSNPs.pl LINE: 59</font></div>
<div class=""><font face="Courier New" class="">Date (localtime) = Mon Feb 29 09:47:44 2016</font></div>
<div class=""><font face="Courier New" class="">Ensembl API version = 83</font></div>
<div class=""><font face="Courier New" class="">---------------------------------------------------</font></div>
</div>
<div class=""><br class="">
</div>
<div class="">Am I using the function on the wrong kind of ldfc? Or have I misunderstood the documentation? The ultimate goal of my script is to find r2 between all SNPs/rsIDs I pass to the script which are on the same chromosome.</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class="">Johanne Håøy Horn</div>
<div class=""><br class="">
</div>
</div>
</body>
</html>