{"id":32850,"date":"2011-03-20T15:42:51","date_gmt":"2011-03-20T15:42:51","guid":{"rendered":"http:\/\/euvolution.com\/futurist-transhuman-news-blog\/looking-for-relatedness-in-the-hapmap-gujaratis-gene-expression\/"},"modified":"2011-03-20T15:42:51","modified_gmt":"2011-03-20T15:42:51","slug":"looking-for-relatedness-in-the-hapmap-gujaratis-gene-expression","status":"publish","type":"post","link":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/astronomy\/looking-for-relatedness-in-the-hapmap-gujaratis-gene-expression.php","title":{"rendered":"Looking for relatedness in the HapMap Gujaratis | Gene Expression"},"content":{"rendered":"<p>Recently I was looking at a 3-D PCA animation which <a rel=\"nofollow\" target=\"_blank\" href=\"http:\/\/www.harappadna.org\/2011\/03\/ref-1-south-asians-harappa-pca\/\">Zack generated<\/a> from the <a rel=\"nofollow\" target=\"_blank\" href=\"http:\/\/www.harappadna.org\/2011\/03\/ref-1-south-asians-harappa-pca\/\">Harappa Ancestry Project<\/a> data set. Click the link and come back. Notice the outlier clusters?&nbsp;The Burusho are straightforward, they seem to have low levels of Tibetan admixture. <strong>But what about the Gujarati cluster?<\/strong> Again, we see what we&rsquo;ve seen before, the fractioning out of the Gujaratis in PCA into two groups, one a tight cluster, and the other relatively widely distributed. This prompted me to look more closely at the HapMap Gujarati sample. Today I was exploring the question with <a rel=\"nofollow\" target=\"_blank\" href=\"http:\/\/pngu.mgh.harvard.edu\/~purcell\/plink\/ibdibs.shtml\">Plink&rsquo;s identity-by-descent<\/a> feature. First I&rsquo;ll start out with a smaller data set, my family (father, mother, sibling 1, sibling 2, and myself), and an Indian (from Uttar Pradesh) and Pakistani as unrelated individuals. I merged out 23andMe derived genotypes, and with ~900,000 markers calculated pairwise IBD:<\/p><p>.\/plink --bfile IBDControl --genome<\/p><p>Here are the relevant results:<\/p><p> <strong>Individual 1<\/strong><br><strong>Individual 2<\/strong><br><strong>Z0<\/strong><br><strong>Z1<\/strong><br><strong>Z2<\/strong><br><strong>PI_HAT<\/strong><br><strong>DST<\/strong><br><strong>PPC<\/strong><br><strong>RATIO<\/strong> Indian<br>Father<br>0.768<br>0.027<br>0.205<br>0.218<br><em>0.760<\/em><br>0.160<br>1.940 Indian<br>Mother<br>0.782<br>0.010<br>0.209<br>0.214<br><em>0.759<\/em><br>0.026<br>1.886 Indian<br>Razib<br>0.767<br>0.032<br>0.202<br>0.218<br><em>0.759<\/em><br>0.500<br>2.000 Indian<br>Sibling1<br>0.769<br>0.025<br>0.206<br>0.219<br><em>0.760<\/em><br>0.198<br>1.949 Indian<br>Sibling2<br>0.766<br>0.032<br>0.203<br>0.219<br><em>0.760<\/em><br>0.685<br>2.030 Indian<br>Pakistani<br>0.781<br>0.017<br>0.203<br>0.211<br><em>0.758<\/em><br>0.533<br>2.005 Father<br>Mother<br>0.776<br>0.018<br>0.207<br>0.215<br><em>0.759<\/em><br>0.284<br>1.965 Father<br>Razib<br>0.002<br>0.777<br>0.221<br>0.610<br><em>0.851<\/em><br>1.000<br>450.800 Father<br>Sibling1<br>0.001<br>0.785<br>0.214<br>0.606<br><em>0.850<\/em><br>1.000<br>898.800 Father<br>Sibling2<br>0.002<br>0.779<br>0.220<br>0.609<br><em>0.851<\/em><br>1.000<br>643.143 Father<br>Pakistani<br>0.778<br>0.019<br>0.203<br>0.213<br><em>0.758<\/em><br>0.201<br>1.950 Mother<br>Razib<br>0.002<br>0.788<br>0.211<br>0.605<br><em>0.849<\/em><br>1.000<br>639.429 Mother<br>Sibling1<br>0.002<br>0.781<br>0.218<br>0.608<br><em>0.850<\/em><br>1.000<br>639.857 Mother<br>Sibling2<br>0.002<br>0.782<br>0.216<br>0.607<br><em>0.850<\/em><br>1.000<br>447.900 Mother<br>Pakistani<br>0.779<br>0.020<br>0.201<br>0.211<br><em>0.758<\/em><br>0.052<br>1.904 Razib<br>Sibling1<br>0.183<br>0.408<br>0.409<br>0.613<br><em>0.866<\/em><br>1.000<br>11.386 Razib<br>Sibling2<br>0.194<br>0.432<br>0.374<br>0.590<br><em>0.858<\/em><br>1.000<br>11.491 Razib<br>Pakistani<br>0.781<br>0.016<br>0.203<br>0.211<br><em>0.758<\/em><br>0.933<br>2.095 Sibling1<br>Sibling2<br>0.236<br>0.412<br>0.351<br>0.557<br><em>0.849<\/em><br>1.000<br>9.413 Sibling1<br>Pakistani<br>0.777<br>0.024<br>0.199<br>0.211<br><em>0.758<\/em><br>0.327<br>1.973 Sibling2<br>Pakistani<br>0.774<br>0.024<br>0.202<br>0.214<br><em>0.758<\/em><br>0.443<br>1.991 <\/p><p>You can infer some things without even knowing what the columns mean. Notice that there are differences between parent-child, sibling-sibling, and unrelated comparisons. The distance measure, DST, is basically exactly the same as the genome-wide comparison in 23andMe. Either the web app is running Plink, or, it&rsquo;s using the ...<img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/euvolution.com\/futurist-transhuman-news-blog\/wp-content\/plugins\/wp-o-matic\/cache\/b42dc_N82KMaJ_-KI\" height=\"1\" width=\"1\" style=\"padding-left:10px; padding-right: 10px;\"><img loading=\"lazy\" decoding=\"async\" src=\"http:\/\/euvolution.com\/futurist-transhuman-news-blog\/wp-content\/plugins\/wp-o-matic\/cache\/b42dc_a_ZlrGJdz-w\" height=\"1\" width=\"1\" style=\"padding-left:10px; padding-right: 10px;\"><\/p>","protected":false},"excerpt":{"rendered":"<p>Recently I was looking at a 3-D PCA animation which Zack generated from the Harappa Ancestry Project data set. Click the link and come back. Notice the outlier clusters?&nbsp;The Burusho are straightforward, they seem to have low levels of Tibetan &hellip; <a href=\"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/astronomy\/looking-for-relatedness-in-the-hapmap-gujaratis-gene-expression.php\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"limit_modified_date":"","last_modified_date":"","_lmt_disableupdate":"","_lmt_disable":"","footnotes":""},"categories":[21],"tags":[],"class_list":["post-32850","post","type-post","status-publish","format-standard","hentry","category-astronomy"],"modified_by":null,"_links":{"self":[{"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/posts\/32850"}],"collection":[{"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/comments?post=32850"}],"version-history":[{"count":0,"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/posts\/32850\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/media?parent=32850"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/categories?post=32850"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.euvolution.com\/futurist-transhuman-news-blog\/wp-json\/wp\/v2\/tags?post=32850"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}