Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Geographical Genomics of Human Leukocyte Gene Expression Variation in Southern Morocco</title><author><name sortKey="Idaghdour, Youssef" sort="Idaghdour, Youssef" uniqKey="Idaghdour Y" first="Youssef" last="Idaghdour">Youssef Idaghdour</name><affiliation><nlm:aff id="A1">Department of Genetics, North Carolina State University, Raleigh NC, USA</nlm:aff></affiliation></author><author><name sortKey="Czika, Wendy" sort="Czika, Wendy" uniqKey="Czika W" first="Wendy" last="Czika">Wendy Czika</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Shianna, Kevin V" sort="Shianna, Kevin V" uniqKey="Shianna K" first="Kevin V." last="Shianna">Kevin V. Shianna</name><affiliation><nlm:aff id="A3">Institute for Genome Science and Policy, Duke University, Durham NC, USA</nlm:aff></affiliation></author><author><name sortKey="Lee, S Hong" sort="Lee, S Hong" uniqKey="Lee S" first="S. Hong" last="Lee">S. Hong Lee</name><affiliation><nlm:aff id="A4">Queensland Institute of Medical Research, Brisbane, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Visscher, Peter M" sort="Visscher, Peter M" uniqKey="Visscher P" first="Peter M." last="Visscher">Peter M. Visscher</name><affiliation><nlm:aff id="A4">Queensland Institute of Medical Research, Brisbane, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Martin, Hilary C" sort="Martin, Hilary C" uniqKey="Martin H" first="Hilary C." last="Martin">Hilary C. Martin</name><affiliation><nlm:aff id="A5">School of Biological Sciences, University of Queensland, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Miclaus, Kelci" sort="Miclaus, Kelci" uniqKey="Miclaus K" first="Kelci" last="Miclaus">Kelci Miclaus</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Jadallah, Sami J" sort="Jadallah, Sami J" uniqKey="Jadallah S" first="Sami J." last="Jadallah">Sami J. Jadallah</name><affiliation><nlm:aff id="A6">HRH Prince Sultan International Foundation for Conservation and Development of Wildlife, Agadir, Morocco</nlm:aff></affiliation></author><author><name sortKey="Goldstein, David B" sort="Goldstein, David B" uniqKey="Goldstein D" first="David B." last="Goldstein">David B. Goldstein</name><affiliation><nlm:aff id="A3">Institute for Genome Science and Policy, Duke University, Durham NC, USA</nlm:aff></affiliation></author><author><name sortKey="Wolfinger, Russell D" sort="Wolfinger, Russell D" uniqKey="Wolfinger R" first="Russell D." last="Wolfinger">Russell D. Wolfinger</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Gibson, Greg" sort="Gibson, Greg" uniqKey="Gibson G" first="Greg" last="Gibson">Greg Gibson</name><affiliation><nlm:aff id="A5">School of Biological Sciences, University of Queensland, Queensland, Australia</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">19966804</idno><idno type="pmc">2798927</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2798927</idno><idno type="RBID">PMC:2798927</idno><idno type="doi">10.1038/ng.495</idno><date when="2009">2009</date><idno type="wicri:Area/Pmc/Corpus">000225</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000225</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Geographical Genomics of Human Leukocyte Gene Expression Variation in Southern Morocco</title><author><name sortKey="Idaghdour, Youssef" sort="Idaghdour, Youssef" uniqKey="Idaghdour Y" first="Youssef" last="Idaghdour">Youssef Idaghdour</name><affiliation><nlm:aff id="A1">Department of Genetics, North Carolina State University, Raleigh NC, USA</nlm:aff></affiliation></author><author><name sortKey="Czika, Wendy" sort="Czika, Wendy" uniqKey="Czika W" first="Wendy" last="Czika">Wendy Czika</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Shianna, Kevin V" sort="Shianna, Kevin V" uniqKey="Shianna K" first="Kevin V." last="Shianna">Kevin V. Shianna</name><affiliation><nlm:aff id="A3">Institute for Genome Science and Policy, Duke University, Durham NC, USA</nlm:aff></affiliation></author><author><name sortKey="Lee, S Hong" sort="Lee, S Hong" uniqKey="Lee S" first="S. Hong" last="Lee">S. Hong Lee</name><affiliation><nlm:aff id="A4">Queensland Institute of Medical Research, Brisbane, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Visscher, Peter M" sort="Visscher, Peter M" uniqKey="Visscher P" first="Peter M." last="Visscher">Peter M. Visscher</name><affiliation><nlm:aff id="A4">Queensland Institute of Medical Research, Brisbane, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Martin, Hilary C" sort="Martin, Hilary C" uniqKey="Martin H" first="Hilary C." last="Martin">Hilary C. Martin</name><affiliation><nlm:aff id="A5">School of Biological Sciences, University of Queensland, Queensland, Australia</nlm:aff></affiliation></author><author><name sortKey="Miclaus, Kelci" sort="Miclaus, Kelci" uniqKey="Miclaus K" first="Kelci" last="Miclaus">Kelci Miclaus</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Jadallah, Sami J" sort="Jadallah, Sami J" uniqKey="Jadallah S" first="Sami J." last="Jadallah">Sami J. Jadallah</name><affiliation><nlm:aff id="A6">HRH Prince Sultan International Foundation for Conservation and Development of Wildlife, Agadir, Morocco</nlm:aff></affiliation></author><author><name sortKey="Goldstein, David B" sort="Goldstein, David B" uniqKey="Goldstein D" first="David B." last="Goldstein">David B. Goldstein</name><affiliation><nlm:aff id="A3">Institute for Genome Science and Policy, Duke University, Durham NC, USA</nlm:aff></affiliation></author><author><name sortKey="Wolfinger, Russell D" sort="Wolfinger, Russell D" uniqKey="Wolfinger R" first="Russell D." last="Wolfinger">Russell D. Wolfinger</name><affiliation><nlm:aff id="A2">SAS Institute Inc., Cary NC, USA</nlm:aff></affiliation></author><author><name sortKey="Gibson, Greg" sort="Gibson, Greg" uniqKey="Gibson G" first="Greg" last="Gibson">Greg Gibson</name><affiliation><nlm:aff id="A5">School of Biological Sciences, University of Queensland, Queensland, Australia</nlm:aff></affiliation></author></analytic><series><title level="j">Nature genetics</title><idno type="ISSN">1061-4036</idno><idno type="eISSN">1546-1718</idno><imprint><date when="2009">2009</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p id="P1">Studies of the genetics of gene expression reveal expression SNPs that explain variation in transcript abundance. Here we address the robustness of eSNP associations to environmental geography and population structure in a comparison of 194 Arab and Amazigh individuals from a city and two villages in southern Morocco. Gene expression differed between pairs of locations for up to a third of all transcripts, with notable enrichment for ribosomal biosynthesis and oxidative phosphorylation. Robust associations were observed in the leukocyte samples with <italic>cis</italic>-eSNPs (<italic>P</italic> < 10<sup>−08</sup>) for 346 genes, and <italic>trans</italic>-eSNPs (<italic>P</italic> < 10<sup>−11</sup>) with 10 genes. All of these were consistent across the three sample locations and after controlling for ethnicity and relatedness. No evidence for large-effect <italic>trans</italic>-acting mediators of the pervasive environmental influence was found and instead genetic and environmental factors acted in a largely additive manner.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Abegunde, Do" uniqKey="Abegunde D">DO Abegunde</name></author><author><name sortKey="Mathers, Cd" uniqKey="Mathers C">CD Mathers</name></author><author><name sortKey="Adam, T" uniqKey="Adam T">T Adam</name></author><author><name sortKey="Ortegon, M" uniqKey="Ortegon M">M Ortegon</name></author><author><name sortKey="Strong, K" uniqKey="Strong K">K Strong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cookson, W" uniqKey="Cookson W">W Cookson</name></author><author><name sortKey="Liang, L" uniqKey="Liang L">L Liang</name></author><author><name sortKey="Abecasis, G" uniqKey="Abecasis G">G Abecasis</name></author><author><name sortKey="Moffatt, M" uniqKey="Moffatt M">M Moffatt</name></author><author><name sortKey="Lathrop, M" uniqKey="Lathrop M">M Lathrop</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Idaghdour, Y" uniqKey="Idaghdour Y">Y Idaghdour</name></author><author><name sortKey="Storey, Jd" uniqKey="Storey J">JD Storey</name></author><author><name sortKey="Jadallah, Sj" uniqKey="Jadallah S">SJ Jadallah</name></author><author><name sortKey="Gibson, G" uniqKey="Gibson G">G Gibson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Arredi, B" uniqKey="Arredi B">B Arredi</name></author><author><name sortKey="Poloni, Es" uniqKey="Poloni E">ES Poloni</name></author><author><name sortKey="Paracchini, S" uniqKey="Paracchini S">S Paracchini</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Feezor, Rj" uniqKey="Feezor R">RJ Feezor</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Price, Al" uniqKey="Price A">AL Price</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fellay, J" uniqKey="Fellay J">J Fellay</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Biswas, S" uniqKey="Biswas S">S Biswas</name></author><author><name sortKey="Scheinfeldt, Lb" uniqKey="Scheinfeldt L">LB Scheinfeldt</name></author><author><name sortKey="Akey, Jm" uniqKey="Akey J">JM Akey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pritchard, J" uniqKey="Pritchard J">J Pritchard</name></author><author><name sortKey="Stephens, M" uniqKey="Stephens M">M Stephens</name></author><author><name sortKey="P Donnelly, P" uniqKey="P Donnelly P">P P Donnelly</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kefir, R" uniqKey="Kefir R">R Kéfir</name></author><author><name sortKey="Stevanovitch, A" uniqKey="Stevanovitch A">A Stevanovitch</name></author><author><name sortKey="Bouzaid, E" uniqKey="Bouzaid E">E Bouzaid</name></author><author><name sortKey="Beraud Colomb, E" uniqKey="Beraud Colomb E">E Béraud-Colomb</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Coudray, C" uniqKey="Coudray C">C Coudray</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ennafaa, H" uniqKey="Ennafaa H">H Ennafaa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bosch, E" uniqKey="Bosch E">E Bosch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wolfinger, Rd" uniqKey="Wolfinger R">RD Wolfinger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Royo, H" uniqKey="Royo H">H Royo</name></author><author><name sortKey="Cavaille, J" uniqKey="Cavaille J">J Cavaillé</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dixon, Al" uniqKey="Dixon A">AL Dixon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Stranger, Be" uniqKey="Stranger B">BE Stranger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cheung, Vg" uniqKey="Cheung V">VG Cheung</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Storey, Jd" uniqKey="Storey J">JD Storey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kao, Cf" uniqKey="Kao C">CF Kao</name></author><author><name sortKey="Chen, Sy" uniqKey="Chen S">SY Chen</name></author><author><name sortKey="Lee, Yh" uniqKey="Lee Y">YH Lee</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ruggero, D" uniqKey="Ruggero D">D Ruggero</name></author><author><name sortKey="Pandolfi, Pp" uniqKey="Pandolfi P">PP Pandolfi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shah, Sv" uniqKey="Shah S">SV Shah</name></author><author><name sortKey="Baliga, R" uniqKey="Baliga R">R Baliga</name></author><author><name sortKey="Rajapurkar, M" uniqKey="Rajapurkar M">M Rajapurkar</name></author><author><name sortKey="Fonseca, Va" uniqKey="Fonseca V">VA Fonseca</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goring, Hh" uniqKey="Goring H">HH Göring</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Heinzen, El" uniqKey="Heinzen E">EL Heinzen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Emilsson, V" uniqKey="Emilsson V">V Emilsson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Heap, Ga" uniqKey="Heap G">GA Heap</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hayes, Bj" uniqKey="Hayes B">BJ Hayes</name></author><author><name sortKey="Visscher, Pm" uniqKey="Visscher P">PM Visscher</name></author><author><name sortKey="Goddard, Me" uniqKey="Goddard M">ME Goddard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schadt, Ee" uniqKey="Schadt E">EE Schadt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, Ls" uniqKey="Chen L">LS Chen</name></author><author><name sortKey="Emmert Streib, F" uniqKey="Emmert Streib F">F Emmert-Streib</name></author><author><name sortKey="Storey, Jd" uniqKey="Storey J">JD Storey</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rockman, Mv" uniqKey="Rockman M">MV Rockman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Menzel, S" uniqKey="Menzel S">S Menzel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sankaran, Vg" uniqKey="Sankaran V">VG Sankaran</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sankaran, Vg" uniqKey="Sankaran V">VG Sankaran</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Melzer, D" uniqKey="Melzer D">D Melzer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pare, G" uniqKey="Pare G">G Paré</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Culverhouse, R" uniqKey="Culverhouse R">R Culverhouse</name></author><author><name sortKey="Suarez, B" uniqKey="Suarez B">B Suarez</name></author><author><name sortKey="Lin, J" uniqKey="Lin J">J Lin</name></author><author><name sortKey="Reich, T" uniqKey="Reich T">T Reich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Visscher, Pm" uniqKey="Visscher P">PM Visscher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Soranzo, N" uniqKey="Soranzo N">N Soranzo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gibson, G" uniqKey="Gibson G">G Gibson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Purcell, S" uniqKey="Purcell S">S Purcell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Visscher, Pm" uniqKey="Visscher P">PM Visscher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Idaghdour, Y" uniqKey="Idaghdour Y">Y Idaghdour</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Thomas, Pd" uniqKey="Thomas P">PD Thomas</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Okuda, S" uniqKey="Okuda S">S Okuda</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<pmc-dir>properties manuscript</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-journal-id">9216904</journal-id><journal-id journal-id-type="pubmed-jr-id">2419</journal-id><journal-id journal-id-type="nlm-ta">Nat Genet</journal-id><journal-id journal-id-type="iso-abbrev">Nat. Genet.</journal-id><journal-title-group><journal-title>Nature genetics</journal-title></journal-title-group><issn pub-type="ppub">1061-4036</issn><issn pub-type="epub">1546-1718</issn></journal-meta><article-meta><article-id pub-id-type="pmid">19966804</article-id><article-id pub-id-type="pmc">2798927</article-id><article-id pub-id-type="doi">10.1038/ng.495</article-id><article-id pub-id-type="manuscript">NIHMS157635</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Geographical Genomics of Human Leukocyte Gene Expression Variation in Southern Morocco</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Idaghdour</surname><given-names>Youssef</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Czika</surname><given-names>Wendy</given-names></name><xref ref-type="aff" rid="A2">2</xref></contrib><contrib contrib-type="author"><name><surname>Shianna</surname><given-names>Kevin V.</given-names></name><xref ref-type="aff" rid="A3">3</xref></contrib><contrib contrib-type="author"><name><surname>Lee</surname><given-names>S. Hong</given-names></name><xref ref-type="aff" rid="A4">4</xref></contrib><contrib contrib-type="author"><name><surname>Visscher</surname><given-names>Peter M.</given-names></name><xref ref-type="aff" rid="A4">4</xref></contrib><contrib contrib-type="author"><name><surname>Martin</surname><given-names>Hilary C.</given-names></name><xref ref-type="aff" rid="A5">5</xref></contrib><contrib contrib-type="author"><name><surname>Miclaus</surname><given-names>Kelci</given-names></name><xref ref-type="aff" rid="A2">2</xref></contrib><contrib contrib-type="author"><name><surname>Jadallah</surname><given-names>Sami J.</given-names></name><xref ref-type="aff" rid="A6">6</xref></contrib><contrib contrib-type="author"><name><surname>Goldstein</surname><given-names>David B.</given-names></name><xref ref-type="aff" rid="A3">3</xref></contrib><contrib contrib-type="author"><name><surname>Wolfinger</surname><given-names>Russell D.</given-names></name><xref ref-type="aff" rid="A2">2</xref></contrib><contrib contrib-type="author"><name><surname>Gibson</surname><given-names>Greg</given-names></name><xref ref-type="aff" rid="A5">5</xref><xref ref-type="corresp" rid="cor1">*</xref></contrib></contrib-group><aff id="A1"><label>1</label>Department of Genetics, North Carolina State University, Raleigh NC, USA</aff><aff id="A2"><label>2</label>SAS Institute Inc., Cary NC, USA</aff><aff id="A3"><label>3</label>Institute for Genome Science and Policy, Duke University, Durham NC, USA</aff><aff id="A4"><label>4</label>Queensland Institute of Medical Research, Brisbane, Queensland, Australia</aff><aff id="A5"><label>5</label>School of Biological Sciences, University of Queensland, Queensland, Australia</aff><aff id="A6"><label>6</label>HRH Prince Sultan International Foundation for Conservation and Development of Wildlife, Agadir, Morocco</aff><author-notes><corresp id="cor1"><label>*</label>Corresponding author: School of Biological Sciences, Goddard Building, St Lucia Campus, University of Queensland, Brisbane, QLD 4072, Australia, <email>ggibson.uq@gmail.com</email>, Ph: +61 7 3365-2194</corresp></author-notes><pub-date pub-type="nihms-submitted"><day>17</day><month>11</month><year>2009</year></pub-date><pub-date pub-type="epub"><day>06</day><month>12</month><year>2009</year></pub-date><pub-date pub-type="ppub"><month>1</month><year>2010</year></pub-date><pub-date pub-type="pmc-release"><day>01</day><month>7</month><year>2010</year></pub-date><volume>42</volume><issue>1</issue><fpage>62</fpage><lpage>67</lpage><pmc-comment>elocation-id from pubmed: 10.1038/ng.495</pmc-comment>
<permissions><license><license-p>Users may view, print, copy, download and text and data- mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use:
<ext-link ext-link-type="uri" xlink:href="http://www.nature.com/authors/editorial_policies/license.html#terms">http://www.nature.com/authors/editorial_policies/license.html#terms</ext-link></license-p></license></permissions><abstract><p id="P1">Studies of the genetics of gene expression reveal expression SNPs that explain variation in transcript abundance. Here we address the robustness of eSNP associations to environmental geography and population structure in a comparison of 194 Arab and Amazigh individuals from a city and two villages in southern Morocco. Gene expression differed between pairs of locations for up to a third of all transcripts, with notable enrichment for ribosomal biosynthesis and oxidative phosphorylation. Robust associations were observed in the leukocyte samples with <italic>cis</italic>-eSNPs (<italic>P</italic> < 10<sup>−08</sup>) for 346 genes, and <italic>trans</italic>-eSNPs (<italic>P</italic> < 10<sup>−11</sup>) with 10 genes. All of these were consistent across the three sample locations and after controlling for ethnicity and relatedness. No evidence for large-effect <italic>trans</italic>-acting mediators of the pervasive environmental influence was found and instead genetic and environmental factors acted in a largely additive manner.</p></abstract><kwd-group><kwd>Peripheral blood</kwd><kwd>eSNP</kwd><kwd>GWAS</kwd><kwd>ethnicity</kwd><kwd>relatedness</kwd><kwd>environmental geography</kwd></kwd-group></article-meta></front><body><p id="P2">The human transition from pastoral and rural to urban lifestyles has been accompanied by increased incidence of numerous chronic diseases such as asthma, diabetes and cancer<xref rid="R1" ref-type="bibr">1</xref>. Environmental contributors, likely including dietary shifts, pollution, and psychological factors, are the subject of ongoing epidemiological research. It is equally interesting to ask whether genetic influences on disease susceptibility change across environments. Since disease risk is commonly thought to often involve differential gene expression<xref rid="R2" ref-type="bibr">2</xref>, we have assessed the robustness of transcript abundance to environmental variation by performing a genome-wide association study on leukocyte gene expression profiles across two ethnicities in three locations. Our previous work had demonstrated a substantial effect of environmental geography<xref rid="R3" ref-type="bibr">3</xref> on gene expression in Moroccan Amazigh, and here we additionally add the contrast with Arabs, allowing us to test whether geography and/or ethnicity affect each of several hundred robust associations between genotypes and transcript abundance.</p><p id="P3">The Souss region in southern Morocco is home to several million people of two dominant ethnicities, living either in cities, or rural villages (<xref ref-type="fig" rid="F1">Fig. 1</xref>). The Amazigh Berbers are descendant of the first modern humans who populated North Africa 35,000 or more years ago<xref rid="R4" ref-type="bibr">4</xref>, and many still live in traditional villages in the low Atlas Mountains. The Arabs by contrast moved into southern Morocco between the 7<sup>th</sup> and 11<sup>th</sup> centuries and tend to occupy lowland villages, while both groups inhabit the cities, often retaining their linguistic and cultural identities.</p><p id="P4">We collected peripheral blood samples from 284 healthy adults in June and July of 2008 from four locations, including approximately equal numbers of men and women, and of Amazigh and Arabs. Half the sample was from two high density, low to middle-income, urban communities, Anza and Dchiera, on either side of the city of Agadir. The other half was from two rural villages near Tiznit, 120 km to the south. Boutroch is predominantly Amazigh and remains quite isolated, while Ighrem is predominantly Arab and (based on self-reported information and our observations at the collection site) many of the men, in particular, commute into the cities.</p><p id="P5">Leukocytes were isolated from serum, platelets and erythrocytes at the time of blood sampling by depletion filter technology<xref rid="R5" ref-type="bibr">5</xref> and fixed in RNALater<sup>®</sup> solution within a few minutes of blood collection. Gene expression profiles were obtained from 208 high quality RNA samples using Illumina HumanHT12 bead arrays that include 48,804 probes, of which 22,300 RefSeq probes in 16,738 genes were deemed to have signal above background. In order to minimize batch effects, all samples were processed in the same week, and the extraction, labeling and hybridization steps were all performed according to a randomized block design. Whole genome genotypes were obtained from whole blood samples using Illumina Human 610-Quad arrays. After quality control filters, 516,972 SNPs were available for 194 of the individuals who also had gene expression profiles.</p><sec sec-type="results" id="S1"><title>RESULTS</title><sec id="S2"><title>Population Structure of southern Morocco</title><p id="P6">Population structure was assessed by examining the principal components (PC) of the variance of the genotype profiles, using Eigenstrat software<xref rid="R6" ref-type="bibr">6</xref>. Initial examination revealed several clusters of siblings and other close relatives (cousins or similar) whose similarity skewed the axes; where data were available, these identities were in agreement with participant records. After removal of these relatives, analysis of 163 unrelated individuals revealed seven significant eigenvectors. None of these explain more than 5% of the variance, and PC3 through PC7 are heavily weighted by large clusters of SNPs on one or a few chromosomes. As described by others, such axes are commonly observed and do not provide reliable genome-wide estimates of population structure<xref rid="R7" ref-type="bibr">7</xref>,<xref rid="R8" ref-type="bibr">8</xref>, but it is interesting to note that PC3 distinguishes Ighrem from the other locations (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 1a online</xref>).</p><p id="P7">A plot of the first two eigenvectors (<xref ref-type="fig" rid="F2">Fig. 2a</xref>) highlights the major historical influences on population structure in southern Morocco. PC1 separates just a dozen individuals, and we inferred that this axis represents a sub-Saharan African contribution, consistent with expected levels of admixture in Morocco, by performing an analysis including 21 Yoruban individuals (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 1b online</xref>). PC2 is highly correlated with both location and self-reported ethnicity, so is inferred to capture the major component of Arab-Amazigh ancestry.</p><p id="P8">A surprising aspect of this analysis is the positioning of Ighrem Arabs between Boutroch Amazigh and half of the Agadir Arabs along PC2. This was confirmed by Structure analysis<xref rid="R9" ref-type="bibr">9</xref> of 16,000 randomly chosen autosomal SNPs assuming admixture of two ancestral populations (<xref ref-type="fig" rid="F2">Fig. 2b</xref>), which indicates that Ighrem residents tend to be a mixture, while most Amazigh are derived from one population, and only a handful of Agadir Arabs represent the other. There has thus likely been considerable admixture between these two groups over an extended period of time, possibly with movement of Arabs from other locations into Agadir recently. A slight shift of Ighrem Arabs toward the Amazigh pole of PC2, relative to Agadir Arabs, would also be consistent with some genetic exchange between the villages over 50 generations. Further sampling of villages in the region may reveal subtle population structure across southern Morocco<xref rid="R10" ref-type="bibr">10</xref>–<xref rid="R13" ref-type="bibr">13</xref>.</p></sec><sec id="S3"><title>Regional Differentiation in Gene Expression</title><p id="P9">Next we asked whether region, location and ethnicity impact gene expression profiles, and if they do so in a gender-specific manner. Since location and ethnicity are confounded in the villages, several parallel analyses were undertaken to tease apart these influences. Transcript abundance data was transformed by median centering on the log base 2 scale (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 2 online</xref>), which results in maximal overlap of profiles without altering their variance.</p><p id="P10">Gene specific analysis of variance<xref rid="R14" ref-type="bibr">14</xref> with expression as a function of region, gender, and their interaction discovered 1,521 probes significant at a 1% false discovery rate (FDR; <italic>P</italic> < 0.0007). Region, namely the rural (Boutroch plus Ighrem) versus city (Anza plus Dchiera) comparison is by far the main effect in this joint analysis. Approaching 7% of all expressed genes differentiate these individuals by this conservative criterion, whereas considerably fewer than 1% of the probes show gender differences. A full list of genes is provided in <xref ref-type="supplementary-material" rid="SD2">Supplementary Table 1 online</xref>. Among several classes of over-represented genes for this lifestyle comparison, small nucleolar RNA genes stand out: 5 of the top 8 overall and 15 of 29 members of the <italic>SNORD</italic> family are in the highly significant list, compared with just 1 of 10 <italic>SNORA</italic> genes. There is little in the literature to indicate why this is the case, or what the physiological consequences may be, but it is interesting to note that epigenetic modification has been observed for many small nucleolar RNA genes<xref rid="R15" ref-type="bibr">15</xref>.</p><p id="P11">Even more differentiation was observed when we fit analysis of variance models including location, gender, and their interaction. Since exploratory analyses indicated that the Anza and Dchiera samples are indistinguishable either for gene expression or genotypes, these were combined into a single location, Agadir, in all subsequent analyses. In the three-way comparison, 8,459 probes (38%) were significant at the 1% FDR threshold for location (<xref ref-type="supplementary-material" rid="SD3">Supplementary Table 2 online</xref>). Boutroch differs from both Ighrem and Agadir at over seven thousand probes each, with a high degree of overlap (<xref ref-type="fig" rid="F3">Fig. 3a</xref>; <xref ref-type="table" rid="T1">Table 1</xref>). Ighrem and Agadir are much more similar to one another, in part because there is considerably more diversity within the Ighrem sample that reduces the significance of the location contrast. We also noted that women are much more differentiated among locations than men (<xref ref-type="table" rid="T1">Table 1</xref>). These results confirm our previous report<xref rid="R2" ref-type="bibr">2</xref> of substantial differentiation between Bedouin nomads, urban Anza, and another remote Amazigh village, Sebt Nabor.</p><p id="P12">In order to evaluate the possible independent contribution of ethnicity more carefully, variance component analysis of the expression variation was performed. Within Agadir alone, neither ethnicity (modeled as the second eigenvector of the genotype data, gPC2) nor gender have a noteworthy impact on the principal components of the expression variation, as shown by the bar charts in <xref ref-type="fig" rid="F3">Figure 3b</xref>. However, in the total dataset there is some evidence for a contribution: <xref ref-type="fig" rid="F3">Figure 3c</xref> shows that when fit jointly with location, the ethnicity, ethnicity-and gender-by-location interaction terms make a substantial contribution to the expression profiles.</p><p id="P13">Although gender and ethnicity affect the expression of fewer genes than location, the plot of expression PC1 by PC2 for the most differentially expressed 1,500 genes in <xref ref-type="fig" rid="F4">Figure 4</xref> indicates that for many genes the interaction between these three factors is quite complex. This can also be seen in the expression profiles of characteristic individual genes (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 3 online</xref>). Boutroch and Ighrem villagers in general separate along PC1, while high values of PC2 are obtained for all Boutroch residents (cluster 1) and for Arab women in Ighrem (cluster 2). Amazigh women from Ighrem (cluster 3) and the Ighrem men (cluster 4) have lower values of PC2 similar to those observed for all Agadir residents. The simplest interpretation is that cultural or behavioral differences, likely including time spent outside the village, contribute strongly to the observed gender and ethnicity effects. Deeper sampling would be required to firmly establish whether intrinsic biological differences between the sexes and/or populations also make significant contributions to expression divergence in lymphocytes, as they appear to do for lymphoblast cell lines grown in culture<xref rid="R16" ref-type="bibr">16</xref>–<xref rid="R19" ref-type="bibr">19</xref>.</p><p id="P14">Two classes of genes stand out as significantly differentially-expressed among locations. These are ribosomal proteins of both the small and large subunits as well as the cytoplasmic and mitochondrial compartments, and proteins involved in oxidative phosphorylation, which are highly up-regulated in half of the Agadir residents (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 4a online</xref>). All of the transcripts encoding these proteins form a module of co-regulated genes, but as shown in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figure 4b online</xref>, it is noteworthy that this module is not co-expressed with the <italic>SNORD</italic>s, which tend to be relatively down-regulated in Agadir but are particularly high in the Arab women from Ighrem. Regulation of ribosomal biosynthesis may be related to response to viral infection, and it also seems to be involved in tumorigenesis in conjunction with mitochondrial activity<xref rid="R20" ref-type="bibr">20</xref>,<xref rid="R21" ref-type="bibr">21</xref>. Oxidative phosphorylation is correlated with renal health and the production or disposal of free radicals<xref rid="R22" ref-type="bibr">22</xref>, so our data suggests that deeper evaluation of health risks associated with lifestyle transitions may be revealing.</p></sec><sec id="S4"><title>Genome-wide association with gene expression variation</title><p id="P15">The genetic contribution to expression variation was evaluated by genome-wide association with expression of all 22,300 probes. Starting with a simple test of the correlation between each transcript abundance and each genotype, and filtering to retain only eSNPs with a minor allele frequency greater than 0.05, we observed 3,430 associations at <italic>P</italic> < 10<sup>−8</sup>. Further filtering of eSNPs to retain only autosomal associations with annotated genes, and imposing the additional stringency of <italic>P</italic> < 10<sup>−11</sup> for putative <italic>trans</italic> associations between an eSNP on one chromosome and a probe on another chromosome, reduced this to 1,636 associations. 1,569 (96%) of these are intra-chromosomal linkages, the vast majority within 50 kb and hence <italic>cis</italic>-acting(<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 5 online</xref>), and only 3 clearly in different chromosomal intervals. Facsimile associations were observed for 39 of the target genes represented by a second probe (37 <italic>cis</italic>, 2 <italic>trans</italic>). Reducing the dataset further to exclude linked associations within haplotype blocks leaves 346unique <italic>cis</italic> and 10 unique <italic>trans</italic> associations at the stringent genome-wide 5% significance level. These proportions are in good agreement with most other GWAS expression studies on blood or lymphocyte cell lines<xref rid="R16" ref-type="bibr">16</xref>,<xref rid="R17" ref-type="bibr">17</xref>,<xref rid="R23" ref-type="bibr">23</xref>–<xref rid="R26" ref-type="bibr">26</xref>, and a 30-fold or greater excess of <italic>cis</italic> over <italic>trans</italic> associations is also supported by 1% FDR estimates of 600 and 20 genes respectively. Complete lists of peak <italic>cis</italic> and <italic>trans</italic> associations are provided in <xref ref-type="supplementary-material" rid="SD4">Supplementary Table 3 online</xref>.</p><p id="P16">Given the high degree of population structure for gene expression, we addressed the possibility that differentiation of eSNP allele frequencies may contribute to the observed associations by calculating <italic>F<sub>ST</sub></italic> estimates for each pair-wise comparison of location for the 516,972 SNPs and 16,500 of the genes. No fixed differences were observed and plots of the <italic>F<sub>ST</sub></italic> comparisons (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 6a online</xref>) indicate only moderate overall genetic differentiation, with occasional SNPs having <italic>F<sub>ST</sub></italic> values between 0.12 and 0.3. There was no tendency for these outliers to have elevated expression differentiation and in fact almost all of the top 10% most differentially expressed genes are among the least genetically differentiated. Nor was there any correlation between <italic>F<sub>ST</sub></italic> and significance of gene expression divergence (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 6b online</xref>), confirming that the observed expression differences between locations are for the most part not attributable to gene-specific allelic frequency differences between locations.</p><p id="P17">The robustness of the 3,430 associations to environmental sources of variance and population structure was further evaluated by fitting two additional linear trend models to the data. The first included location, gender and the interaction between them. The second included two measures of ethnicity (the first three genotype eigenvectors and a four-way categorical ethnicity cluster: see methods), a matrix of relatedness based on an identity by descent measure<xref rid="R27" ref-type="bibr">27</xref>, as well as gender interactions with ethnicity cluster and genotype. <xref ref-type="fig" rid="F5">Figures 5a and 5b</xref> show the Manhattan plot of associations by chromosomal location for the second of these models, and the <italic>cis-trans</italic> plot of target against eSNP location, respectively. <xref ref-type="fig" rid="F5">Figure 5c</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 7 online</xref> show that the logarithm of the genotype significance term is highly correlated (<italic>r</italic> > 0.95) between both of these models and the original correlation test. Furthermore, <xref ref-type="fig" rid="F5">Figure 5d</xref> shows that there is no evidence for significant genotype-by-location interactions in any of the association trend tests. Neither the ethnicity nor the relatedness variance components explain an appreciable amount of the expression variation for any of the transcripts (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 8 online</xref>).</p><p id="P18">The absence of interaction effects is readily visualized by plotting expression as a function of genotype with color coding of each location, for each association. An example of a <italic>trans</italic> association in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figure 9 online</xref> shows the clear trend of increased expression of <italic>AMY1A</italic> (chromosome 1) in homozygotes for the A allele of <italic>ACTG1</italic> gamma actin (chromosome 17), consistently across the three locations despite slight overall location effects. Expression of <italic>AMY1A</italic> is highly correlated with that of <italic>AMY1B</italic> (<italic>r</italic> > 0.8) as well as of dozens of other genes in a co-expression module, but the eSNP only regulates <italic>AMY1A</italic>, because it increases expression of the gene two-fold in an additive manner. A similar plot for a representative gene that shows highly significant location and genotype effects in <italic>cis</italic>, C21ORF57, is provided in <xref ref-type="fig" rid="F6">Figure 6a</xref> and further discussed below, and further examples can be seen in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figure 9c online</xref>.</p></sec><sec id="S5"><title>Novel associations with potential disease alleles</title><p id="P19">GWAS-expression associations detected in one tissue can identify regulatory variants that may be active in other tissues that are directly engaged in the etiology of disease<xref rid="R23" ref-type="bibr">23</xref>,<xref rid="R25" ref-type="bibr">25</xref>,<xref rid="R26" ref-type="bibr">26</xref>. One example is the <italic>cis</italic>-linkages in peripheral blood with the T1D susceptibility locus at chromosome 12q13. The strongest expression association is with transcription of the <italic>RPS26</italic> ribosomal protein gene, and network analyses have been employed to argue that this is the more likely diabetes candidate gene than the initially reported <italic>ERBB3</italic><xref rid="R28" ref-type="bibr">28</xref>. However, the strongest T1D association involves a different SNP than that associated with expression and/or splicing<xref rid="R24" ref-type="bibr">24</xref> of <italic>RPS26</italic>. We further find that the same linkage group of eSNPs, centered on the rs10876864 in the <italic>SUOX</italic> gene 35kb from <italic>RPS26</italic>, is also associated in <italic>trans</italic> with half a dozen other RP26 paralogs (probably due to cross-hybridization), and with <italic>CCDC4</italic> on chromosome 4, albeit at the suggestive significance level of<italic>P</italic> = 3.5×10<sup>−10</sup>. Intriguingly, expression of <italic>RPS26</italic> is only weakly correlated with that of the module of ribosomal proteins that differentiate locations (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 4b online</xref>), so this association does not contribute to the environmental effect on transcription of ribosomal protein genes.</p><p id="P20">Another <italic>trans</italic>-association of interest involves rs11987927 in <italic>MYOM2</italic> at 8p23 with <italic>ZNF71</italic> at 19q13, but also with its own <italic>MYOM2</italic> transcript. Logic suggests that the <italic>cis</italic>-association likely affects the abundance of the MYOM2 myomesin protein, which in turn regulates <italic>ZNF71</italic>, but the <italic>trans</italic> association is actually significantly stronger and conditional dependence analysis<xref rid="R29" ref-type="bibr">29</xref>,<xref rid="R30" ref-type="bibr">30</xref> points in the opposite direction, namely that the <italic>MYOM2</italic> regulatory site influences <italic>ZNF71</italic>, which then feeds back on the <italic>MYOM2</italic> transcript (<xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 10 online</xref>). This example may be a cautionary tale concerning the interpretation of conditional dependence results. It is worth mentioning that four of the seven strongest <italic>trans</italic> associations involve regulation by loci that include genes that encode structural proteins, the others being the <italic>LAMA5</italic> laminin (20q13) with <italic>OSBPL2</italic>, and the <italic>PLEKHM1</italic> plekstrin homology domain protein (17q21) with <italic>MAPK8IP1</italic>.</p><p id="P21">One further <italic>trans</italic>-association is of particular interest. Prolongation of fetal gamma hemoglobin expression in adults is often observed in thalassemia patients. We found association of two probes that detect both <italic>HBG1</italic> and <italic>HBG2</italic> transcripts from 11p15 with rs766432 in the second intron of the <italic>BCL11A</italic> zinc finger proto-oncogene at 2p16. This same SNP has previously been associated with the fraction of erythrocytes that contain measurable fetal hemoglobin<xref rid="R31" ref-type="bibr">31</xref>, and alteration of BCL11A activity was recently shown drive differences in globin switching between mice and humans<xref rid="R32" ref-type="bibr">32</xref>. Another SNP in <italic>BCL11A</italic>, rs4671393 has been associated with abundance of two <italic>BCL11A</italic> transcript isoforms in CEU and YRI HapMap lymphoblast cell lines<xref rid="R33" ref-type="bibr">33</xref>, but is not associated with <italic>BCL11A</italic> transcript abundance in our leukocyte data, suggesting that regulation of BCL11A translation or protein activity is more likely to be affecting <italic>HBG</italic> expression in our sample.</p><p id="P22">Numerous <italic>cis</italic>-associations are also likely to be of interest. We scanned the GWAS association database for overlap between our study and established disease associations at <italic>p</italic><10<sup>−5</sup>. Of 1,628 entries, 10 involve <italic>cis</italic> associations observed in our dataset that explain between 15 and 55% of the transcript variance (<xref ref-type="supplementary-material" rid="SD5">Supplementary Table 4 online</xref>). Five of the associations are with disease conditions (rheumatoid arthritis, celiac disease, T1D, ulcerative colitis, and SLE) and five are with endophenotypes (PAFAH1B2 and ICAM-1 protein levels, triglycerides, LDL cholesterol, and hip bone mineral density). The two serum protein associations<xref rid="R34" ref-type="bibr">34</xref>,<xref rid="R35" ref-type="bibr">35</xref> are with the same SNPs as we detect and hence suggest that protein abundance is largely regulated at the transcriptional level.</p></sec></sec><sec sec-type="discussion" id="S6"><title>DISCUSSION</title><sec id="S7"><title>The genetic and environmental contributions to expression variation</title><p id="P23">Our geographical genomic survey of gene expression variation in southern Morocco has highlighted two parallel and for the most part non-overlapping insights. On the one hand, it is evident that as much as half of the transcriptome is influenced by the environment in a highly coordinated manner such that where a person lives explains up to a quarter of the variation for a substantial fraction of the transcripts. The environmental influences are likely a combination of biotic and abiotic factors, as well as cultural and behavioral ones, while genetic differences between the two North African ethnicities are relatively minor. On the other hand, the genome is littered with strong genetic associations, mainly in <italic>cis</italic>, that explain between 15 and 60 percent of the variance of 5% of the transcripts. Impressive as these associations are, particularly since they are discovered in a sample of just under 200 individuals, they have essentially no bearing on the vast majority of the transcriptional variation, and are not informative of the genetic basis of the environmental response.</p><p id="P24">The robustness of the observed associations to the environmental effect raises the issue of whether genotype-by-environment interactions influence the peripheral blood transcriptome at all. Genome-wide significant interaction effects are generally unlikely to occur in the absence of significant main genotype effects<xref rid="R36" ref-type="bibr">36</xref>. The only circumstances in which they will are if the genotype effect is in the opposite direction in two locations, and if the genetic effect in these locations is at least the same magnitude as the main effects detected in this GWAS, namely explaining over 30% of the variance of a particular transcript. While a few such interactions may exist, it would take a study comparing several thousand individuals from each location to reveal weaker genotype-by-environment interactions. If the genetic architecture of transcription is generally similar to that of visible phenotypes like height and body mass<xref rid="R37" ref-type="bibr">37</xref>,<xref rid="R38" ref-type="bibr">38</xref>, even such a study will be underpowered to explain the vast majority of transcriptional variance.</p><p id="P25">A related question is whether or not genotype-by-environment interactions at the level of transcription are necessary to explain genotype-by-environment interactions for disease. It is possible the small interactions beneath the level of detection of GWAS are prevalent, or alternatively that disease arises primarily as a result of rare alleles of major effect, whose penetrance may be modulated in an environment-specific manner. However, transcriptional interactions are not required to explain the increased incidence of chronic disease. It is not difficult to imagine that individuals that fall into the major categories of transcriptome profiles (such as those implicated in <xref ref-type="fig" rid="F4">Fig. 4</xref> and <xref ref-type="supplementary-material" rid="SD1">Supplementary Fig. 4 online</xref>) have different distributions of disease susceptibility that alter the genotype-disease association matrix genome-wide, thereby inducing environment-by-genotype interactions for disease. Transcription of some genes that contribute to this expression component may also correlate with disease directly, effectively uncovering cryptic variation and resulting in environment-specific eSNP disease associations without any interaction effect at the level of transcription (<xref ref-type="fig" rid="F6">Fig. 6</xref>)<xref rid="R39" ref-type="bibr">39</xref>. A corollary of this is that gene expression profiling might be used to stratify individuals at elevated risk for disease, thereby increasing the resolution of genome-wide association studies by focusing attention on the subset of individuals where genetic effects on disease are most pronounced.</p></sec></sec><sec sec-type="supplementary-material" id="SM"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="SD1"><label>1</label><media xlink:href="NIHMS157635-supplement-1.pdf" mimetype="application" mime-subtype="pdf" xlink:type="simple" id="d37e772" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD2"><label>2</label><media xlink:href="NIHMS157635-supplement-2.xls" mimetype="application" mime-subtype="vnd.ms-excel" xlink:type="simple" id="d37e776" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD3"><label>3</label><media xlink:href="NIHMS157635-supplement-3.xls" mimetype="application" mime-subtype="vnd.ms-excel" xlink:type="simple" id="d37e780" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD4"><label>4</label><media xlink:href="NIHMS157635-supplement-4.xls" mimetype="application" mime-subtype="vnd.ms-excel" xlink:type="simple" id="d37e784" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD5"><label>5</label><media xlink:href="NIHMS157635-supplement-5.xls" mimetype="application" mime-subtype="vnd.ms-excel" xlink:type="simple" id="d37e788" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD6"><label>6</label><media xlink:href="NIHMS157635-supplement-6.xls" mimetype="application" mime-subtype="vnd.ms-excel" xlink:type="simple" id="d37e792" position="anchor"></media></supplementary-material><supplementary-material content-type="local-data" id="SD7"><label>01</label><media xlink:href="NIHMS157635-supplement-01.pdf" mimetype="application" mime-subtype="pdf" xlink:type="simple" id="d37e796" position="anchor"></media></supplementary-material></sec></body><back><fn-group><fn id="FN1"><p id="P26"><bold>ACCESSION NUMBERS</bold></p><p id="P27">Accession codes. NCBI GEO: Gene expression data from this study have been deposited under series GSE17065.</p></fn><fn id="FN2" fn-type="con"><p id="P28"><bold>AUTHOR CONTRIBUTIONS</bold></p><p id="P29">YI collected the samples with the assistance of SJ and processed them with KS and DG; KM, HL, DG, PV and RW provided statistical and conceptual support for analysis of the data by YI, WC, HM and GG; and YI and GG conceived the study and wrote the paper.</p></fn><fn id="FN3"><p id="P30"><bold>COMPETING INTEREST:</bold> Co-authors Wolfinger, Czika and Miclaus are all employees of SAS, Inc, which are the producers of commercial JMP Genomics software used in the analysis of the data.</p></fn></fn-group><ack id="S8"><title>ACKNOWLEDGMENTS</title><p id="P31">We would like to thank all of the study participants in Agadir, Ighrem and Boutroch, as well as numerous individuals who facilitated sample collection, in particular Malika, Aicha and Ahmed Idaghdour. Dongliang Ge and Alison Motsinger-Reif provided timely computational support, and we also thank Shameek Biswas and Josh Akey for providing HapMap genotypes. Funding for the study was provided by the University of Queensland, with YI supported by a Fulbright Fellowship and GG by an ARC Australian Professorial Fellowship.</p></ack><ref-list><title>References</title><ref id="R1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abegunde</surname><given-names>DO</given-names></name><name><surname>Mathers</surname><given-names>CD</given-names></name><name><surname>Adam</surname><given-names>T</given-names></name><name><surname>Ortegon</surname><given-names>M</given-names></name><name><surname>Strong</surname><given-names>K</given-names></name></person-group><article-title>The burden and costs of chronic diseases in low-incomeand middle-income countries</article-title><source>Lancet</source><year>2007</year><volume>370</volume><fpage>1929</fpage><lpage>1938</lpage><pub-id pub-id-type="pmid">18063029</pub-id></element-citation></ref><ref id="R2"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cookson</surname><given-names>W</given-names></name><name><surname>Liang</surname><given-names>L</given-names></name><name><surname>Abecasis</surname><given-names>G</given-names></name><name><surname>Moffatt</surname><given-names>M</given-names></name><name><surname>Lathrop</surname><given-names>M</given-names></name></person-group><article-title>Mapping complex disease traits with global gene expression</article-title><source>Nat Rev Genet</source><year>2009</year><volume>10</volume><fpage>184</fpage><lpage>194</lpage><pub-id pub-id-type="pmid">19223927</pub-id></element-citation></ref><ref id="R3"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Idaghdour</surname><given-names>Y</given-names></name><name><surname>Storey</surname><given-names>JD</given-names></name><name><surname>Jadallah</surname><given-names>SJ</given-names></name><name><surname>Gibson</surname><given-names>G</given-names></name></person-group><article-title>A genome-wide gene expression signature of environmental geography in leukocytes of Moroccan Amazighs</article-title><source>PLoS Genet</source><year>2008</year><volume>4</volume><fpage>e52</fpage></element-citation></ref><ref id="R4"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arredi</surname><given-names>B</given-names></name><name><surname>Poloni</surname><given-names>ES</given-names></name><name><surname>Paracchini</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>A predominantly Neolithic origin for Y-chromosomal DNA variation in North Africa</article-title><source>Am. J. Hum. Genet</source><year>2004</year><volume>75</volume><fpage>338</fpage><lpage>345</lpage><pub-id pub-id-type="pmid">15202071</pub-id></element-citation></ref><ref id="R5"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Feezor</surname><given-names>RJ</given-names></name><etal></etal></person-group><article-title>Whole blood and leukocyte RNA isolation for gene expression analyses</article-title><source>Physiol Genom</source><year>2004</year><volume>19</volume><fpage>247</fpage><lpage>254</lpage></element-citation></ref><ref id="R6"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Price</surname><given-names>AL</given-names></name><etal></etal></person-group><article-title>Principal components analysis corrects for stratification in genome-wide association studies</article-title><source>Nat Genet</source><year>2006</year><volume>38</volume><fpage>904</fpage><lpage>909</lpage><pub-id pub-id-type="pmid">16862161</pub-id></element-citation></ref><ref id="R7"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fellay</surname><given-names>J</given-names></name><etal></etal></person-group><article-title>A whole-genome association study of major determinants for host control of HIV-1</article-title><source>Science</source><year>2007</year><volume>317</volume><fpage>944</fpage><lpage>947</lpage><pub-id pub-id-type="pmid">17641165</pub-id></element-citation></ref><ref id="R8"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Biswas</surname><given-names>S</given-names></name><name><surname>Scheinfeldt</surname><given-names>LB</given-names></name><name><surname>Akey</surname><given-names>JM</given-names></name></person-group><article-title>Genome-wide insights into the patterns and determinants of fine-scale population structure in humans</article-title><source>Am. J. Hum. Genet</source><year>2009</year><volume>84</volume><fpage>641</fpage><lpage>650</lpage><pub-id pub-id-type="pmid">19442770</pub-id></element-citation></ref><ref id="R9"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pritchard</surname><given-names>J</given-names></name><name><surname>Stephens</surname><given-names>M</given-names></name><name><surname>P Donnelly</surname><given-names>P</given-names></name></person-group><article-title>Inference of population structure using multilocus genotype data</article-title><source>Genetics</source><year>2000</year><volume>155</volume><fpage>945</fpage><lpage>959</lpage><pub-id pub-id-type="pmid">10835412</pub-id></element-citation></ref><ref id="R10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kéfir</surname><given-names>R</given-names></name><name><surname>Stevanovitch</surname><given-names>A</given-names></name><name><surname>Bouzaid</surname><given-names>E</given-names></name><name><surname>Béraud-Colomb</surname><given-names>E</given-names></name></person-group><article-title>Diversité mitochondriale de la population de Taforalt (12.000 ans bp - Maroc): Une approche génétique à l'étude du peuplement de l'Afrique du nord</article-title><source>Anthropologie</source><year>2005</year><volume>43</volume><fpage>1</fpage><lpage>11</lpage></element-citation></ref><ref id="R11"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coudray</surname><given-names>C</given-names></name><etal></etal></person-group><article-title>Population genetic data of 15 tetrameric short tandem repeats (STRs) in Berbers from Morocco</article-title><source>Forensic Sci. Int</source><year>2007</year><volume>167</volume><fpage>81</fpage><lpage>86</lpage><pub-id pub-id-type="pmid">16464552</pub-id></element-citation></ref><ref id="R12"><label>12</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ennafaa</surname><given-names>H</given-names></name><etal></etal></person-group><article-title>Mitochondrial DNA haplogroup H structure in North Africa</article-title><source>BMC Genet</source><year>2009</year><volume>10</volume><fpage>8</fpage><pub-id pub-id-type="pmid">19243582</pub-id></element-citation></ref><ref id="R13"><label>13</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bosch</surname><given-names>E</given-names></name><etal></etal></person-group><article-title>Population history of North Africa: evidence from classical genetic markers</article-title><source>Hum. Biol</source><year>1997</year><volume>69</volume><fpage>295</fpage><lpage>311</lpage><pub-id pub-id-type="pmid">9164042</pub-id></element-citation></ref><ref id="R14"><label>14</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wolfinger</surname><given-names>RD</given-names></name><etal></etal></person-group><article-title>Assessing gene significance from cDNA microarray expression data via mixed models</article-title><source>J. Comput. Biol</source><year>2001</year><volume>8</volume><fpage>625</fpage><lpage>637</lpage><pub-id pub-id-type="pmid">11747616</pub-id></element-citation></ref><ref id="R15"><label>15</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Royo</surname><given-names>H</given-names></name><name><surname>Cavaillé</surname><given-names>J</given-names></name></person-group><article-title>Non-coding RNAs in imprinted gene clusters</article-title><source>Biol. Cell</source><year>2008</year><volume>100</volume><fpage>149</fpage><lpage>166</lpage><pub-id pub-id-type="pmid">18271756</pub-id></element-citation></ref><ref id="R16"><label>16</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dixon</surname><given-names>AL</given-names></name><etal></etal></person-group><article-title>A genome-wide association study of global gene expression</article-title><source>Nat. Genet</source><year>2007</year><volume>39</volume><fpage>1202</fpage><lpage>1207</lpage><pub-id pub-id-type="pmid">17873877</pub-id></element-citation></ref><ref id="R17"><label>17</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stranger</surname><given-names>BE</given-names></name><etal></etal></person-group><article-title>Population genomics of human gene expression</article-title><source>Nat. Genet</source><year>2007</year><volume>39</volume><fpage>1217</fpage><lpage>1224</lpage><pub-id pub-id-type="pmid">17873874</pub-id></element-citation></ref><ref id="R18"><label>18</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cheung</surname><given-names>VG</given-names></name><etal></etal></person-group><article-title>Mapping determinants of human gene expression by regional and genome-wide association</article-title><source>Nature</source><year>2005</year><volume>437</volume><fpage>1365</fpage><lpage>1369</lpage><pub-id pub-id-type="pmid">16251966</pub-id></element-citation></ref><ref id="R19"><label>19</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Storey</surname><given-names>JD</given-names></name><etal></etal></person-group><article-title>Gene-expression variation within and among human populations</article-title><source>Am. J. Hum. Genet</source><year>2007</year><volume>80</volume><fpage>502</fpage><lpage>509</lpage><pub-id pub-id-type="pmid">17273971</pub-id></element-citation></ref><ref id="R20"><label>20</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kao</surname><given-names>CF</given-names></name><name><surname>Chen</surname><given-names>SY</given-names></name><name><surname>Lee</surname><given-names>YH</given-names></name></person-group><article-title>Activation of RNA polymerase I transcription by hepatitis C virus core protein</article-title><source>J. Biomed. Sci</source><year>2004</year><volume>11</volume><fpage>72</fpage><lpage>94</lpage><pub-id pub-id-type="pmid">14730212</pub-id></element-citation></ref><ref id="R21"><label>21</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ruggero</surname><given-names>D</given-names></name><name><surname>Pandolfi</surname><given-names>PP</given-names></name></person-group><article-title>Does the ribosome translate cancer?</article-title><source>Nat. Rev. Cancer</source><year>2003</year><volume>3</volume><fpage>179</fpage><lpage>192</lpage><pub-id pub-id-type="pmid">12612653</pub-id></element-citation></ref><ref id="R22"><label>22</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shah</surname><given-names>SV</given-names></name><name><surname>Baliga</surname><given-names>R</given-names></name><name><surname>Rajapurkar</surname><given-names>M</given-names></name><name><surname>Fonseca</surname><given-names>VA</given-names></name></person-group><article-title>Oxidants in chronic kidney disease</article-title><source>J. Am. Soc. Nephrol</source><year>2007</year><volume>18</volume><fpage>16</fpage><lpage>28</lpage><pub-id pub-id-type="pmid">17167116</pub-id></element-citation></ref><ref id="R23"><label>23</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Göring</surname><given-names>HH</given-names></name><etal></etal></person-group><article-title>Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytes</article-title><source>Nat. Genet</source><year>2007</year><volume>39</volume><fpage>1208</fpage><lpage>1216</lpage><pub-id pub-id-type="pmid">17873875</pub-id></element-citation></ref><ref id="R24"><label>24</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Heinzen</surname><given-names>EL</given-names></name><etal></etal></person-group><article-title>Tissue-specific genetic control of splicing: implications for the study of complex traits</article-title><source>PLoS Biol</source><year>2008</year><volume>6</volume><fpage>e1000001</fpage></element-citation></ref><ref id="R25"><label>25</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Emilsson</surname><given-names>V</given-names></name><etal></etal></person-group><article-title>Genetics of gene expression and its effect on disease</article-title><source>Nature</source><year>2008</year><volume>452</volume><fpage>423</fpage><lpage>428</lpage><pub-id pub-id-type="pmid">18344981</pub-id></element-citation></ref><ref id="R26"><label>26</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Heap</surname><given-names>GA</given-names></name><etal></etal></person-group><article-title>Complex nature of SNP genotype effects on gene expression in primary human leucocytes</article-title><source>BMC Med Genomics</source><year>2009</year><volume>2</volume><fpage>1</fpage><pub-id pub-id-type="pmid">19128478</pub-id></element-citation></ref><ref id="R27"><label>27</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hayes</surname><given-names>BJ</given-names></name><name><surname>Visscher</surname><given-names>PM</given-names></name><name><surname>Goddard</surname><given-names>ME</given-names></name></person-group><article-title>Increased accuracy of artificial selection by using the realized relationship matrix</article-title><source>Genet. Res</source><year>2009</year><volume>91</volume><fpage>47</fpage><lpage>60</lpage></element-citation></ref><ref id="R28"><label>28</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schadt</surname><given-names>EE</given-names></name><etal></etal></person-group><article-title>Mapping the genetic architecture of gene expression in human liver</article-title><source>PLoS Biol</source><year>2009</year><volume>6</volume><fpage>e107</fpage><pub-id pub-id-type="pmid">18462017</pub-id></element-citation></ref><ref id="R29"><label>29</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>LS</given-names></name><name><surname>Emmert-Streib</surname><given-names>F</given-names></name><name><surname>Storey</surname><given-names>JD</given-names></name></person-group><article-title>Harnessing naturally randomized transcription to infer regulatory relationships among genes</article-title><source>Genome Biol</source><year>2007</year><volume>8</volume><fpage>R219</fpage><pub-id pub-id-type="pmid">17931418</pub-id></element-citation></ref><ref id="R30"><label>30</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rockman</surname><given-names>MV</given-names></name></person-group><article-title>Reverse engineering the genotype-phenotype map with natural genetic variation</article-title><source>Nature</source><year>2008</year><volume>456</volume><fpage>738</fpage><lpage>744</lpage><pub-id pub-id-type="pmid">19079051</pub-id></element-citation></ref><ref id="R31"><label>31</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Menzel</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>A QTL influencing F cell production maps to a gene encoding a zinc-finger protein on chromosome 2p15</article-title><source>Nat. Genet</source><year>2007</year><volume>39</volume><fpage>1197</fpage><lpage>1199</lpage><pub-id pub-id-type="pmid">17767159</pub-id></element-citation></ref><ref id="R32"><label>32</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sankaran</surname><given-names>VG</given-names></name><etal></etal></person-group><article-title>Developmental and species-divergent globin switching are driven by BCL11A</article-title><source>Nature</source><year>2009</year><volume>460</volume><fpage>1093</fpage><lpage>1097</lpage><pub-id pub-id-type="pmid">19657335</pub-id></element-citation></ref><ref id="R33"><label>33</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sankaran</surname><given-names>VG</given-names></name><etal></etal></person-group><article-title>Human fetal hemoglobin expression is regulated by the developmental stage-specific repressor BCL11A</article-title><source>Science</source><year>2008</year><volume>322</volume><fpage>1839</fpage><lpage>1842</lpage><pub-id pub-id-type="pmid">19056937</pub-id></element-citation></ref><ref id="R34"><label>34</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Melzer</surname><given-names>D</given-names></name><etal></etal></person-group><article-title>A Genome-Wide Association Study Identifies Protein Quantitative Trait Loci (pQTLs)</article-title><source>PLoS Genet</source><year>2008</year><volume>4</volume><fpage>1000072</fpage></element-citation></ref><ref id="R35"><label>35</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Paré</surname><given-names>G</given-names></name><etal></etal></person-group><article-title>Novel association of ABO histo-blood group antigen with soluble ICAM-1: results of a genome-wide association study of 6,578 women</article-title><source>PLoS Genet</source><year>2008</year><volume>4</volume><fpage>e1000118</fpage><pub-id pub-id-type="pmid">18604267</pub-id></element-citation></ref><ref id="R36"><label>36</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Culverhouse</surname><given-names>R</given-names></name><name><surname>Suarez</surname><given-names>B</given-names></name><name><surname>Lin</surname><given-names>J</given-names></name><name><surname>Reich</surname><given-names>T</given-names></name></person-group><article-title>A perspective on epistasis: limits of models displaying no main effect</article-title><source>Am. J. Hum. Genet</source><year>2002</year><volume>70</volume><fpage>461</fpage><lpage>471</lpage><pub-id pub-id-type="pmid">11791213</pub-id></element-citation></ref><ref id="R37"><label>37</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Visscher</surname><given-names>PM</given-names></name></person-group><article-title>Sizing up human height variation</article-title><source>Nat Genet</source><year>2008</year><volume>40</volume><fpage>489</fpage><lpage>490</lpage><pub-id pub-id-type="pmid">18443579</pub-id></element-citation></ref><ref id="R38"><label>38</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Soranzo</surname><given-names>N</given-names></name><etal></etal></person-group><article-title>Meta-analysis of genome-wide scans for human adult stature identifies novel Loci and associations with measures of skeletal frame size</article-title><source>PLoS Genet</source><year>2009</year><volume>5</volume><fpage>e1000445</fpage><pub-id pub-id-type="pmid">19343178</pub-id></element-citation></ref><ref id="R39"><label>39</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gibson</surname><given-names>G</given-names></name></person-group><article-title>Decanalization and the origin of complex disease</article-title><source>Nat. Rev. Genet</source><year>2009</year><volume>10</volume><fpage>134</fpage><lpage>140</lpage><pub-id pub-id-type="pmid">19119265</pub-id></element-citation></ref><ref id="R40"><label>40</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Purcell</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>PLINK: a tool set for whole-genome association and population-based linkage analyses</article-title><source>Am. J. Hum. Genet</source><year>2007</year><volume>81</volume><fpage>559</fpage><lpage>575</lpage><pub-id pub-id-type="pmid">17701901</pub-id></element-citation></ref><ref id="R41"><label>41</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Visscher</surname><given-names>PM</given-names></name><etal></etal></person-group><article-title>Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings</article-title><source>PLoS Genet</source><year>2006</year><volume>2</volume><fpage>e41</fpage><pub-id pub-id-type="pmid">16565746</pub-id></element-citation></ref><ref id="R42"><label>42</label><element-citation publication-type="thesis"><person-group person-group-type="author"><name><surname>Idaghdour</surname><given-names>Y</given-names></name></person-group><source>Genetic and Environmental Components of Human Leukocyte Gene Expression Variation in Morocco. PhD thesis</source><comment>available from the NCSU electronic library at <ext-link ext-link-type="uri" xlink:href="http://www.lib.ncsu.edu/ETD-db/ETD-browse/browse?first_letter=I">http://www.lib.ncsu.edu/ETD-db/ETD-browse/browse?first_letter=I</ext-link></comment></element-citation></ref><ref id="R43"><label>43</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thomas</surname><given-names>PD</given-names></name><etal></etal></person-group><article-title>PANTHER: a library of protein families and subfamilies indexed by function</article-title><source>Genome Res</source><year>2003</year><volume>13</volume><fpage>2129</fpage><lpage>2141</lpage><pub-id pub-id-type="pmid">12952881</pub-id></element-citation></ref><ref id="R44"><label>44</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Okuda</surname><given-names>S</given-names></name><etal></etal></person-group><article-title>KEGG Atlas mapping for global analysis of metabolic pathways</article-title><source>Nucleic Acids Res</source><year>2008</year><volume>36</volume><issue>Web Server issue</issue><fpage>W423</fpage><lpage>W426</lpage><pub-id pub-id-type="pmid">18477636</pub-id></element-citation></ref></ref-list></back><floats-group><fig id="F1" position="float"><label>Figure 1</label><caption><title>Map of the Souss region of southern Morocco</title><p id="P32">showing the location of the two rural villages, Boutroch and Ighrem, near the town of Tiznit, relative to the urban locations Anza and Dchiera north and south of the city of Agadir, respectively.</p></caption><graphic xlink:href="nihms157635f1"></graphic></fig><fig id="F2" position="float"><label>Figure 2</label><caption><title>Population structure in southern Morocco</title><p id="P33">(<bold>a</bold>)Eigenstrat principal component analysis of 579,144 SNPs reveals 7 significant eigenvectors, the first two of which, explaining just 1.3 and 0.8 % of the genotypic variance respectively, are plotted here. By self-report, Boutroch Amazigh are blue squares, Agadir Amazigh green triangles, Agadir Arabs green plus symbols, Ighrem Arabsred circles, and Ighrem Amazigh red triangles. 3 individuals with uncertain ethnicity possibly including sub-Saharan heritage, are indicated as gray spots, and have high values of PC1, which is characteristic of Yoruban ancestry as shown in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figure 1b online</xref>. (<bold>b</bold>) Structure analysis of 16,000 autosomal SNPs, with k=3 and employing the admixture model with correlated allele frequencies, highlights the same individuals with large PC1 values (brown bars) and shows that Boutroch Amazigh are predominantly derived from one population group (pale blue) while all other samples are a mixture of the two populations represented by pale red and blue bars.</p></caption><graphic xlink:href="nihms157635f2"></graphic></fig><fig id="F3" position="float"><label>Figure 3</label><caption><title>Location impacts gene expression transcriptome-wide</title><p id="P34">(<bold>a</bold>) Venn diagram of the number of genes significant at 1% FDR for ANOVA of the three pair-wise comparisons indicated. Variance components of expression variation (<bold>b</bold>) just in the 118 residents of Agadir (excluding 9 individuals with strongly positive gPC1 scores, and including reassignment of ethnicity according to gPC2 for just 11 individuals relative to self-report, <xref ref-type="supplementary-material" rid="SD6">Supplementary Table 5</xref>),where Ethnicity is modeled as the PC2 of the genotype variation as shown in <xref ref-type="fig" rid="F1">Figure 1a</xref>, or (<bold>c</bold>) for all 22,300 probes in the full sample of 208 individuals.</p></caption><graphic xlink:href="nihms157635f3"></graphic></fig><fig id="F4" position="float"><label>Figure 4</label><caption><title>Principal component plot for the most differentially expressed genes</title><p id="P35">The two major principal components of the expression of the 1,500 most significant genes shows significant separation of individuals by location (PC1 and PC2) and gender (PC2) (all <italic>P</italic> < 0.0001) as described in the text. Individuals from Boutroch are blue, Ighrem red, and Agadir green. Arabs are indicated with solid spots, Amazigh open circles, and males are lighter symbols for each color. Boutroch and Arab women from Ighrem (clusters 1 and 2) separate from Amazigh women and Arab men from Ighrem (clusters 3 and 4) who are closer to Agadir residents. If Boutroch residents and Ighrem Arab women are grouped and contrasted with Agadir residents, Ighrem Amazigh women, and Ighrem men, 8,239 genes are significantly differentially expressed at the 1% FDR rate, more than any pair-wise comparison of locations. A similar plot for all genes is shown in <xref ref-type="supplementary-material" rid="SD1">Supplementary Figure 11</xref>.</p></caption><graphic xlink:href="nihms157635f4"></graphic></fig><fig id="F5" position="float"><label>Figure 5</label><caption><title>Genome-wide association with transcript abundance</title><p id="P36">(<bold>a</bold>)Manhattan plot of all 1,636 genome-wide associations at <italic>P</italic> < 10<sup>−8</sup>(NLP > 8) for model 3, which includes control for genotype-determined ethnicity, location, relatedness, and gender. Each chromosome is indicated by a different color. The horizontal red line indicates the genome-wide significance threshold (NLP > 11.4) for <italic>trans</italic> associations. Note the excess of peaks at the MHC complex on chromosome 6 due to multiple cis-eSNPs. (<bold>b</bold>) Cis-<italic>Trans</italic> plot showing target transcript location against eSNP location indicating that most eSNPs are in <italic>cis</italic> to the regulated transcript, while just 13 <italic>trans</italic> associations at NLP > 11.4are visible. (<bold>c</bold>) High correlation of significance measures for all eSNPs detected by simple correlation of genotype with expression (model 1) or robust control for ethnicity, gender and location (model 3). (<bold>d</bold>) Absence of genome-wide significance for the Genotype-by-Location interaction effect, which is not correlated with the Genotype effect.</p></caption><graphic xlink:href="nihms157635f5"></graphic></fig><fig id="F6" position="float"><label>Figure 6</label><caption><title>The relationship between genotype, expression, and phenotype</title><p id="P37">(<bold>a</bold>) A typical example of a transcript (encoding C21ORF57, a putative metallo-proteinase) that shows both a significant difference between locations (<italic>P</italic> < 10<sup>−5</sup>) and a <italic>cis</italic>-eSNP association, with rs1556337 (<italic>P</italic> < 10<sup>−13</sup>) but no interaction effect in an additive model on the log scale. Expression is lower in Boutroch (blue points and line), while genotype has a consistent effect across all three locations (Ighrem, red; Agadir, green). (<bold>b</bold>) The Actual vs Predicted plot separates the genotypes by location for clarity. Suppose that a disease or phenotype is only seen in individuals with transcript abundance less than 1.0 (on a relative log2 scale), indicated by the gray area. Then in Agadir and Ighrem (green and red respectively) almost all affected are AA homozygotes, whereas in Boutroch (blue) heterozygotes and some GG homozygotes are also affected. There is thus a G×E interaction for the phenotype in the absence of a G×E interaction for transcription, because the environment shifts more individuals into the susceptible zone. Similar arguments would apply for phenotypes with high expression values, and for graded rather than threshold-dependent traits.</p></caption><graphic xlink:href="nihms157635f6"></graphic></fig><table-wrap id="T1" position="float" orientation="landscape"><label>Table 1</label><caption><p id="P38">Number of transcripts significant at 1% FDR</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left" rowspan="1" colspan="1">Location</th><th align="right" rowspan="1" colspan="1"></th><th align="right" rowspan="1" colspan="1"></th><th align="left" rowspan="1" colspan="1">Gender</th><th align="right" rowspan="1" colspan="1"></th><th align="right" rowspan="1" colspan="1"></th><th align="left" rowspan="1" colspan="1">Interaction</th><th align="right" rowspan="1" colspan="1"></th><th align="right" rowspan="1" colspan="1"></th></tr><tr><th align="left" rowspan="1" colspan="1">ANOVA</th><th align="right" colspan="2" rowspan="1">ANCOVA</th><th align="left" rowspan="1" colspan="1"></th><th align="right" rowspan="1" colspan="1">ANOVA</th><th align="right" rowspan="1" colspan="1">ANCOVA</th><th align="left" rowspan="1" colspan="1"></th><th align="right" rowspan="1" colspan="1">ANOVA</th><th align="right" rowspan="1" colspan="1">ANCOVA</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">3-way</td><td align="right" rowspan="1" colspan="1">8459</td><td align="right" rowspan="1" colspan="1">7057</td><td align="left" rowspan="1" colspan="1">Male : Female</td><td align="right" rowspan="1" colspan="1">151</td><td align="right" rowspan="1" colspan="1">233</td><td align="left" rowspan="1" colspan="1">Location*Gender</td><td align="right" rowspan="1" colspan="1">133</td><td align="right" rowspan="1" colspan="1">203</td></tr><tr><td align="left" rowspan="1" colspan="1">Aga : Bou</td><td align="right" rowspan="1" colspan="1">6744</td><td align="right" rowspan="1" colspan="1">4974</td><td align="left" rowspan="1" colspan="1">In Agadir</td><td align="right" rowspan="1" colspan="1">24</td><td align="right" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">Fem (Aga : Bou)</td><td align="right" rowspan="1" colspan="1">4830</td><td align="right" rowspan="1" colspan="1">3791</td></tr><tr><td align="left" rowspan="1" colspan="1">Aga : Igh</td><td align="right" rowspan="1" colspan="1">635</td><td align="right" rowspan="1" colspan="1">651</td><td align="left" rowspan="1" colspan="1">In Boutroch</td><td align="right" rowspan="1" colspan="1">13</td><td align="right" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">Fem (Aga : Igh)</td><td align="right" rowspan="1" colspan="1">1451</td><td align="right" rowspan="1" colspan="1">1467</td></tr><tr><td align="left" rowspan="1" colspan="1">Bou : Igh</td><td align="right" rowspan="1" colspan="1">7339</td><td align="right" rowspan="1" colspan="1">6286</td><td align="left" rowspan="1" colspan="1">In Ighrem</td><td align="right" rowspan="1" colspan="1">589</td><td align="right" rowspan="1" colspan="1">890</td><td align="left" rowspan="1" colspan="1">Mal (Aga : Bou)<xref ref-type="table-fn" rid="TFN2">†</xref></td><td align="right" rowspan="1" colspan="1">407</td><td align="right" rowspan="1" colspan="1">806</td></tr><tr><td align="left" rowspan="1" colspan="1">Aga : Rural</td><td align="right" rowspan="1" colspan="1">1521</td><td align="right" rowspan="1" colspan="1">607</td><td align="left" rowspan="1" colspan="1"></td><td align="right" rowspan="1" colspan="1"></td><td align="right" rowspan="1" colspan="1"></td><td align="left" rowspan="1" colspan="1">Mal (Aga : Igh)</td><td align="right" rowspan="1" colspan="1">8</td><td align="right" rowspan="1" colspan="1">8</td></tr></tbody></table><table-wrap-foot><fn id="TFN1"><p id="P39">ANOVA includes terms for Location, Gender, and Location*Gender interaction. The False Discovery Rate was evaluated using the conservative Benjamin and Hochberg method. The left hand columns show the number of genes significant at the 1% FDR threshold for Location effects (either in the 3-way comparison of Agadir (Aga), Boutroch (Bou) and Ighrem (Igh); between pairs of locations, or between Agadir (Aga) and the two rural sites combined). The central columns contrast gender (male versus female) effects, either in the total sample or each location individually. The right hand columns show interaction effects, either in the total sample or showing the indicated contrast between Agadir and either village, for females or males separately. ANCOVA is the same model with an additional continuous covariate for ethnicity, genotypic PC2.</p></fn><fn id="TFN2"><label>†</label><p id="P40">Significance of this contrast was reduced by the small sample of Boutroch males (12, cf 26 females).</p></fn></table-wrap-foot></table-wrap></floats-group></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Linguistique/explor/TamazightV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000225 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000225 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Linguistique
|area= TamazightV2
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |