High-Altitude Andean H194R HIF2A Allele Is a Hypomorphic Allele.

For over 10,000 years, Andeans have resided at high altitude where the partial pressure of oxygen challenges human survival. Recent studies have provided evidence for positive selection acting in Andeans on the HIF2A (also known as EPAS1) locus, which encodes for a central transcription factor of the hypoxia-inducible factor pathway. However, the precise mechanism by which this allele might lead to altitude-adaptive phenotypes, if any, is unknown. By analyzing whole genome sequencing data from 46 high-coverage Peruvian Andean genomes, we confirm evidence for positive selection acting on HIF2A and a unique pattern of variation surrounding the Andean-specific single nucleotide variant (SNV), rs570553380, which encodes for an H194R amino acid substitution in HIF-2α. Genotyping the Andean-associated SNV rs570553380 in a group of 299 Peruvian Andeans from Cerro de Pasco, Peru (4,338 m), reveals a positive association with increased fraction of exhaled nitric oxide, a marker of nitric oxide biosynthesis. In vitro assays show that the H194R mutation impairs binding of HIF-2α to its heterodimeric partner, aryl hydrocarbon receptor nuclear translocator. A knockin mouse model bearing the H194R mutation in the Hif2a gene displays decreased levels of hypoxia-induced pulmonary Endothelin-1 transcripts and protection against hypoxia-induced pulmonary hypertension. We conclude the Andean H194R HIF2A allele is a hypomorphic (partial loss of function) allele.

Differential splenic responses to hyperoxic breathing at high altitude in Sherpa and lowlanders.

The human spleen contracts in response to stress-induced catecholamine secretion, resulting in a temporary rise in haemoglobin concentration ([Hb]). Recent findings highlighted enhanced splenic response to exercise at high altitude in Sherpa, possibly due to a blunted splenic response to hypoxia. To explore the potential blunted splenic contraction in Sherpas at high altitude, we examined changes in spleen volume during hyperoxic breathing, comparing acclimatized Sherpa with acclimatized individuals of lowland ancestry. Our study included 14 non-Sherpa (7 female) residing at altitude for a mean continuous duration of 3 months and 46 Sherpa (24 female) with an average of 4 years altitude exposure. Participants underwent a hyperoxic breathing test at altitude (4300 m; barrometric pressure = âˆ¼430 torr; P O 2 ${P_{{{\mathrm{O}}_{\mathrm{2}}}}}$  = âˆ¼90 torr). Throughout the test, we measured spleen volume using ultrasonography and monitored oxygen saturation ( S p O 2 ${S_{{\mathrm{p}}{{\mathrm{O}}_{\mathrm{2}}}}}$ ). During rest, Sherpa exhibited larger spleens (226 ± 70 mL) compared to non-Sherpa (165 ± 34 mL; P < 0.001; effect size (ES) = 0.95, 95% CI: 0.3-1.6). In response to hyperoxia, non-Sherpa demonstrated 22 ± 12% increase in spleen size (35 ± 17 mL, 95% CI: 20.7-48.9; P < 0.001; ES = 1.8, 95% CI: 0.93-2.66), while spleen size remained unchanged in Sherpa (-2 ± 13 mL, 95% CI: -2.4 to 7.3; P = 0.640; ES = 0.18, 95% CI: -0.10 to 0.47). Our findings suggest that Sherpa and non-Sherpas of lowland ancestry exhibit distinct variations in spleen volume during hyperoxia at high altitude, potentially indicating two distinct splenic functions. In Sherpa, this phenomenon may signify a diminished splenic response to altitude-related hypoxia at rest, potentially contributing to enhanced splenic contractions during physical stress. Conversely, non-Sherpa experienced a transient increase in spleen size during hyperoxia, indicating an active tonic contraction, which may influence early altitude acclimatization in lowlanders by raising [Hb].

Larger spleens and greater splenic contraction during exercise may be an adaptive characteristic of Nepali Sherpa at high-altitude.

OBJECTIVES: The Sherpa ethnic group living at altitude in Nepal may have experienced natural selection in response to chronic hypoxia. We have previously shown that Sherpa in Kathmandu (1400 m) possess larger spleens and a greater apnea-induced splenic contraction compared to lowland Nepalis. This may be significant for exercise capacity at altitude as the human spleen responds to stress-induced catecholamine secretion by an immediate contraction, which results in transiently elevated hemoglobin concentration ([Hb]). METHODS: To investigate splenic contraction in response to exercise at high-altitude (4300 m; Pb = ~450 Torr), we recruited 63 acclimatized Sherpa (29F) and 14 Nepali non-Sherpa (7F). Spleen volume was measured before and after maximal exercise on a cycle ergometer by ultrasonography, along with [Hb] and oxygen saturation (SpO2). RESULTS: Resting spleen volume was larger in the Sherpa compared with Nepali non-Sherpa (237 ± 62 vs. 165 ± 34 mL, p < .001), as was the exercise-induced splenic contraction (Δspleen volume, 91 ± 40 vs. 38 ± 32 mL, p < .001). From rest to exercise, [Hb] increased (1.2 to 1.4 g.dl-1), SpO2 decreased (~9%) and calculated arterial oxygen content (CaO2) remained stable, but there were no significant differences between groups. In Sherpa, both resting spleen volume and the Δspleen volume were modest positive predictors of the change (Δ) in [Hb] and CaO2 with exercise (p-values from .026 to .037 and R2 values from 0.059 to 0.067 for the predictor variable). CONCLUSIONS: Larger spleens and greater splenic contraction may be an adaptive characteristic of Nepali Sherpa to increase CaO2 during exercise at altitude, but the direct link between spleen size/function and hypoxia tolerance remains unclear.

Tibetan PHD2, an allele with loss-of-function properties.

Tibetans have adapted to the chronic hypoxia of high altitude and display a distinctive suite of physiologic adaptations, including augmented hypoxic ventilatory response and resistance to pulmonary hypertension. Genome-wide studies have consistently identified compelling genetic signatures of natural selection in two genes of the Hypoxia Inducible Factor pathway, PHD2 and HIF2A The product of the former induces the degradation of the product of the latter. Key issues regarding Tibetan PHD2 are whether it is a gain-of-function or loss-of-function allele, and how it might contribute to high-altitude adaptation. Tibetan PHD2 possesses two amino acid changes, D4E and C127S. We previously showed that in vitro, Tibetan PHD2 is defective in its interaction with p23, a cochaperone of the HSP90 pathway, and we proposed that Tibetan PHD2 is a loss-of-function allele. Here, we report that additional PHD2 mutations at or near Asp-4 or Cys-127 impair interaction with p23 in vitro. We find that mice with the Tibetan Phd2 allele display augmented hypoxic ventilatory response, supporting this loss-of-function proposal. This is phenocopied by mice with a mutation in p23 that abrogates the PHD2:p23 interaction. Hif2a haploinsufficiency, but not the Tibetan Phd2 allele, ameliorates hypoxia-induced increases in right ventricular systolic pressure. The Tibetan Phd2 allele is not associated with hemoglobin levels in mice. We propose that Tibetans possess genetic alterations that both activate and inhibit selective outputs of the HIF pathway to facilitate successful adaptation to the chronic hypoxia of high altitude.

High-altitude deer mouse hypoxia-inducible factor-2α shows defective interaction with CREB-binding protein.

Numerous mammalian species have adapted to the chronic hypoxia of high altitude. Recent genomic studies have identified evidence for natural selection of genes and associated genetic changes in these species. A major gap in our knowledge is an understanding of the functional significance, if any, of these changes. Deer mice (Peromyscus maniculatus) live at both low and high altitudes in North America, providing an opportunity to identify functionally important genetic changes. High-altitude deer mice show evidence of natural selection on the Epas1 gene, which encodes for hypoxia-inducible factor-2α (Hif-2α), a central transcription factor of the hypoxia-inducible factor pathway. An SNP encoding for a T755M change in the Hif-2α protein is highly enriched in high-altitude deer mice, but its functional significance is unknown. Here, using coimmunoprecipitation and transcriptional activity assays, we show that the T755M mutation produces a defect in the interaction of Hif-2α with the transcriptional coactivator CREB-binding protein. This results in a loss of function because of decreased transcriptional activity. Intriguingly, the effect of this mutation depends on the amino acid context. Interchanges between methionine and threonine at the corresponding position in house mouse (Mus musculus) Hif-2α are without effects on CREB-binding protein binding. Furthermore, transfer of a set of deer mouse-specific Hif-2α amino acids to house mouse Hif-2α is sufficient to confer sensitivity of house mouse Hif-2α to the T755M substitution. These findings provide insight into high-altitude adaptation in deer mice and evolution at the Epas1 locus.

Association of EGLN1 gene with high aerobic capacity of Peruvian Quechua at high altitude.

Highland native Andeans have resided at altitude for millennia. They display high aerobic capacity (VO2max) at altitude, which may be a reflection of genetic adaptation to hypoxia. Previous genomewide (GW) scans for natural selection have nominated Egl-9 homolog 1 gene (EGLN1) as a candidate gene. The encoded protein, EGLN1/PHD2, is an O2 sensor that controls levels of the Hypoxia Inducible Factor-α (HIF-α), which regulates the cellular response to hypoxia. From GW association and analysis of covariance performed on a total sample of 429 Peruvian Quechua and 94 US lowland referents, we identified 5 EGLN1 SNPs associated with higher VO2max (L⋅min-1 and mL⋅min-1⋅kg-1) in hypoxia (rs1769793, rs2064766, rs2437150, rs2491403, rs479200). For 4 of these SNPs, Quechua had the highest frequency of the advantageous (high VO2max) allele compared with 25 diverse lowland comparison populations from the 1000 Genomes Project. Genotype effects were substantial, with high versus low VO2max genotype categories differing by ∼11% (e.g., for rs1769793 SNP genotype TT = 34.2 mL⋅min-1⋅kg-1 vs. CC = 30.5 mL⋅min-1⋅kg-1). To guard against spurious association, we controlled for population stratification. Findings were replicated for EGLN1 SNP rs1769793 in an independent Andean sample collected in 2002. These findings contextualize previous reports of natural selection at EGLN1 in Andeans, and support the hypothesis that natural selection has increased the frequency of an EGLN1 causal variant that enhances O2 delivery or use during exercise at altitude in Peruvian Quechua.

Human high-altitude adaptation: forward genetics meets the HIF pathway.

Humans have adapted to the chronic hypoxia of high altitude in several locations, and recent genome-wide studies have indicated a genetic basis. In some populations, genetic signatures have been identified in the hypoxia-inducible factor (HIF) pathway, which orchestrates the transcriptional response to hypoxia. In Tibetans, they have been found in the HIF2A (EPAS1) gene, which encodes for HIF-2α, and the prolyl hydroxylase domain protein 2 (PHD2, also known as EGLN1) gene, which encodes for one of its key regulators, PHD2. High-altitude adaptation may be due to multiple genes that act in concert with one another. Unraveling their mechanism of action can offer new therapeutic approaches toward treating common human diseases characterized by chronic hypoxia.

Identifying adaptive alleles in the human genome: from selection mapping to functional validation.

The suite of phenotypic diversity across geographically distributed human populations is the outcome of genetic drift, gene flow, and natural selection throughout human evolution. Human genetic variation underlying local biological adaptations to selective pressures is incompletely characterized. With the emergence of population genetics modeling of large-scale genomic data derived from diverse populations, scientists are able to map signatures of natural selection in the genome in a process known as selection mapping. Inferred selection signals further can be used to identify candidate functional alleles that underlie putative adaptive phenotypes. Phenotypic association, fine mapping, and functional experiments facilitate the identification of candidate adaptive alleles. Functional investigation of candidate adaptive variation using novel techniques in molecular biology is slowly beginning to unravel how selection signals translate to changes in biology that underlie the phenotypic spectrum of our species. In addition to informing evolutionary hypotheses of adaptation, the discovery and functional annotation of adaptive alleles also may be of clinical significance. While selection mapping efforts in non-European populations are growing, there remains a stark under-representation of diverse human populations in current public genomic databases, of both clinical and non-clinical cohorts. This lack of inclusion limits the study of human biological variation. Identifying and functionally validating candidate adaptive alleles in more global populations is necessary for understanding basic human biology and human disease.

Complex signatures of natural selection at GYPA.

The human MN blood group antigens are isoforms of glycophorin A (GPA) encoded by the gene, GYPA, and are the most abundant erythrocyte sialoglycoproteins. The distribution of MN antigens has been widely studied in human populations yet the evolutionary and/or demographic factors affecting population variation remain elusive. While the primary function of GPA is yet to be discovered, it serves as the major binding site for the 175-kD erythrocyte-binding antigen (EB-175) of the malarial parasite, Plasmodium falciparum, a major selective pressure in recent human history. More specifically, exon two of GYPA encodes the receptor-binding ligand to which P. falciparum binds. Accordingly, there has been keen interest in understanding what impact, if any, natural selection has had on the distribution of variation in GYPA and exon two in particular. To this end, we resequenced GYPA in individuals sampled from both P. falciparum endemic (sub-Saharan Africa and South India) and non-endemic (Europe and East Asia) regions of the world. Observed patterns of variation suggest that GYPA has been subject to balancing selection in populations living in malaria endemic areas and in Europeans, but no such evidence was found in samples from East Asia, Oceania, and the Americas. These results are consistent with malaria acting as a selective pressure on GYPA, but also suggest that another selective force has resulted in a similar pattern of variation in Europeans. Accordingly, GYPA has perhaps a more complex evolutionary history, wherein on a global scale, spatially varying selective pressures have governed its natural history.

Estimates of continental ancestry vary widely among individuals with the same mtDNA haplogroup.

The association between a geographical region and an mtDNA haplogroup(s) has provided the basis for using mtDNA haplogroups to infer an individual's place of origin and genetic ancestry. Although it is well known that ancestry inferences using mtDNA haplogroups and those using genome-wide markers are frequently discrepant, little empirical information exists on the magnitude and scope of such discrepancies between multiple mtDNA haplogroups and worldwide populations. We compared genetic-ancestry inferences made by mtDNA-haplogroup membership to those made by autosomal SNPs in ∼940 samples of the Human Genome Diversity Panel and recently admixed populations from the 1000 Genomes Project. Continental-ancestry proportions often varied widely among individuals sharing the same mtDNA haplogroup. For only half of mtDNA haplogroups did the highest average continental-ancestry proportion match the highest continental-ancestry proportion of a majority of individuals with that haplogroup. Prediction of an individual's mtDNA haplogroup from his or her continental-ancestry proportions was often incorrect. Collectively, these results indicate that for most individuals in the worldwide populations sampled, mtDNA-haplogroup membership provides limited information about either continental ancestry or continental region of origin.

Exome sequencing as a tool for Mendelian disease gene discovery.

Exome sequencing - the targeted sequencing of the subset of the human genome that is protein coding - is a powerful and cost-effective new tool for dissecting the genetic basis of diseases and traits that have proved to be intractable to conventional gene-discovery strategies. Over the past 2 years, experimental and analytical approaches relating to exome sequencing have established a rich framework for discovering the genes underlying unsolved Mendelian disorders. Additionally, exome sequencing is being adapted to explore the extent to which rare alleles explain the heritability of complex diseases and health-related traits. These advances also set the stage for applying exome and whole-genome sequencing to facilitate clinical diagnosis and personalized disease-risk profiling.

Toll-like receptor gene variants and bacterial vaginosis among HIV-1 infected and uninfected African women.

Bacterial vaginosis (BV) is a common vaginal syndrome associated with altered microflora that increases the risk of preterm delivery and acquisition of sexually transmitted diseases. The cause of BV is unknown although toll-like receptors (TLRs), that are central to innate immune responses, may be important. We evaluated associations between TLR SNPs and BV among HIV-1 infected and uninfected African women. Logistic regression was used to assess associations between SNPs (N=99) in TLRs 2-4, 7-9 and BV (as classified by Nugent's criteria). Among HIV-1 uninfected women, TLR7 rs5743737 and TLR7 rs1634323 were associated with a decreased risk of BV, whereas TLR7 rs179012 was associated with an increased risk. TLR2 SNP rs3804099 was associated with a decreased risk of BV among HIV-1 infected women. Our findings indicate that there may be differences in TLR association with BV among HIV-1 infected and HIV-1 uninfected women.

Defective Tibetan PHD2 binding to p23 links high altitude adaption to altered oxygen sensing.

The Tibetan population has adapted to the chronic hypoxia of high altitude. Tibetans bear a genetic signature in the prolyl hydroxylase domain protein 2 (PHD2/EGLN1) gene, which encodes for the central oxygen sensor of the hypoxia-inducible factor (HIF) pathway. Recent studies have focused attention on two nonsynonymous coding region substitutions, D4E and C127S, both of which are markedly enriched in the Tibetan population. These amino acids reside in a region of PHD2 that harbors a zinc finger, which we have previously discovered binds to a Pro-Xaa-Leu-Glu (PXLE) motif in the HSP90 cochaperone p23, thereby recruiting PHD2 to the HSP90 pathway to facilitate HIF-α hydroxylation. We herein report that the Tibetan PHD2 haplotype (D4E/C127S) strikingly diminishes the interaction of PHD2 with p23, resulting in impaired PHD2 down-regulation of the HIF pathway. The defective binding to p23 depends on both the D4E and C127S substitutions. We also identify a PXLE motif in HSP90 itself that can mediate binding to PHD2 but find that this interaction is maintained with the D4E/C127S PHD2 haplotype. We propose that the Tibetan PHD2 variant is a loss of function (hypomorphic) allele, leading to augmented HIF activation to facilitate adaptation to high altitude.

Targeted capture and massively parallel sequencing of 12 human exomes.

Genome-wide association studies suggest that common genetic variants explain only a modest fraction of heritable risk for common diseases, raising the question of whether rare variants account for a significant fraction of unexplained heritability. Although DNA sequencing costs have fallen markedly, they remain far from what is necessary for rare and novel variants to be routinely identified at a genome-wide scale in large cohorts. We have therefore sought to develop second-generation methods for targeted sequencing of all protein-coding regions ('exomes'), to reduce costs while enriching for discovery of highly penetrant variants. Here we report on the targeted capture and massively parallel sequencing of the exomes of 12 humans. These include eight HapMap individuals representing three populations, and four unrelated individuals with a rare dominantly inherited disorder, Freeman-Sheldon syndrome (FSS). We demonstrate the sensitive and specific identification of rare and common variants in over 300 megabases of coding sequence. Using FSS as a proof-of-concept, we show that candidate genes for Mendelian disorders can be identified by exome sequencing of a small number of unrelated, affected individuals. This strategy may be extendable to diseases with more complex genetics through larger sample sizes and appropriate weighting of non-synonymous variants by predicted functional impact.

Toll-like receptor polymorphism associations with HIV-1 outcomes among sub-Saharan Africans.

OBJECTIVE: We evaluated Toll-like receptors (TLRs) single nucleotide polymorphisms (SNPs) for associations with HIV-1 acquisition, set-point and disease progression in African couples. METHODS: Seven candidate and 116 haplotype-tagging SNPs (tagSNPs) were genotyped in 504 HIV-1 infected cases, and 343 seronegative controls. RESULTS: TLR9 1635A/G was associated with reduced HIV-1 acquisition among HIV-seronegative controls with high but not low HIV-1 exposure (odds ratio [OR] = 0.7; P = .03 and OR = 0.9, P = .5, respectively). TLR7 rs179012 and TLR2 597C/T reduced set-point; the latter modified by time since HIV-1 acquisition. TLR8 1A/G reduced disease progression. CONCLUSIONS: TLR SNPs impact HIV-1 outcomes with epidemiologic factors modifying these relationships.

Maternal PRKAA1 and EDNRA genotypes are associated with birth weight, and PRKAA1 with uterine artery diameter and metabolic homeostasis at high altitude.

Low birth weight and intrauterine growth restriction (IUGR) increase the risk of mortality and morbidity during the perinatal period as well as in adulthood. Environmental and genetic factors contribute to IUGR, but the influence of maternal genetic variation on birth weight is largely unknown. We implemented a gene-by-environment study wherein we utilized the growth restrictive effects of high altitude. Multigenerational high-altitude residents (Andeans) are protected from altitude-associated IUGR compared with recent migrants (Europeans). Using a combined cohort of low- and high-altitude European and Andean women, we tested 63 single nucleotide polymorphisms (SNPs) from 16 natural selection-nominated candidate gene regions for associations with infant birth weight. We identified significant SNP associations with birth weight near coding regions for two genes involved in oxygen sensing and vascular control, PRKAA1 and EDNRA, respectively. Next, we identified a significant association for the PRKAA1 SNP with an intermediate phenotype, uterine artery diameter, which has been shown to be related to Andean protection from altitude-associated reductions in fetal growth. To explore potential functional relationships for the effect of maternal SNP genotype on birth weight, we evaluated the relationship between maternal PRKAA1 SNP genotype and gene expression patterns in general and, in particular, of key pathways involved in metabolic homeostasis that have been proposed to play a role in the pathophysiology of IUGR. Our observations suggest that maternal genetic variation within genes that regulate oxygen sensing, metabolic homeostasis, and vascular control influence fetal growth and birth weight outcomes and hence Andean adaptation to high altitude.

Andean and Tibetan patterns of adaptation to high altitude.

OBJECTIVES: High-altitude hypoxia, or decreased oxygen levels caused by low barometric pressure, challenges the ability of humans to live and reproduce. Despite these challenges, human populations have lived on the Andean Altiplano and the Tibetan Plateau for millennia and exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. We and others have identified natural selection candidate genes and gene regions for these adaptations using dense genome scan data. One gene previously known to be important in cellular oxygen sensing, egl nine homolog 1 (EGLN1), shows evidence of positive selection in both Tibetans and Andeans. Interestingly, the pattern of variation for this gene differs between the two populations. Continued research among Tibetan populations has identified statistical associations between hemoglobin concentration and single nucleotide polymorphism (SNP) genotype at EGLN1 and a second gene, endothelial PAS domain protein 1 (EPAS1). METHODS: To measure for the effects of EGLN1 and EPAS1 altitude genotypes on hemoglobin concentration among Andean highlanders, we performed a multiple linear regression analysis of 10 candidate SNPs in or near these two genes. RESULTS: Our analysis did not identify significant associations between EPAS1 or EGLN1 SNP genotypes and hemoglobin concentration in Andeans. CONCLUSIONS: These results contribute to our understanding of the unique set of adaptations developed in different highland groups to the hypoxia of high altitude. Overall, the results provide key insights into the patterns of genetic adaptation to high altitude in Andean and Tibetan populations.

The Legacy of Infectious Disease Exposure on the Genomic Diversity of Indigenous Southern Mexicans.

To characterize host risk factors for infectious disease in Mesoamerican populations, we interrogated 857,481 SNPs assayed using the Affymetrix 6.0 genotyping array for signatures of natural selection in immune response genes. We applied three statistical tests to identify signatures of natural selection: locus-specific branch length (LSBL), the cross-population extended haplotype homozygosity (XP-EHH), and the integrated haplotype score (iHS). Each of the haplotype tests (XP-EHH and iHS) were paired with LSBL and significance was determined at the 1% level. For the paired analyses, we identified 95 statistically significant windows for XP-EHH/LSBL and 63 statistically significant windows for iHS/LSBL. Among our top immune response loci, we found evidence of recent directional selection associated with the major histocompatibility complex (MHC) and the peroxisome proliferator-activated receptor gamma (PPAR-γ) signaling pathway. These findings illustrate that Mesoamerican populations' immunity has been shaped by exposure to infectious disease. As targets of selection, these variants are likely to encode phenotypes that manifest themselves physiologically and therefore may contribute to population-level variation in immune response. Our results shed light on past selective events influencing the host response to modern diseases, both pathogenic infection as well as autoimmune disorders.

OPRM1 and EGFR contribute to skin pigmentation differences between Indigenous Americans and Europeans.

Contemporary variation in skin pigmentation is the result of hundreds of thousands years of human evolution in new and changing environments. Previous studies have identified several genes involved in skin pigmentation differences among African, Asian, and European populations. However, none have examined skin pigmentation variation among Indigenous American populations, creating a critical gap in our understanding of skin pigmentation variation. This study investigates signatures of selection at 76 pigmentation candidate genes that may contribute to skin pigmentation differences between Indigenous Americans and Europeans. Analysis was performed on two samples of Indigenous Americans genotyped on genome-wide SNP arrays. Using four tests for natural selection--locus-specific branch length (LSBL), ratio of heterozygosities (lnRH), Tajima's D difference, and extended haplotype homozygosity (EHH)--we identified 14 selection-nominated candidate genes (SNCGs). SNPs in each of the SNCGs were tested for association with skin pigmentation in 515 admixed Indigenous American and European individuals from regions of the Americas with high ground-level ultraviolet radiation. In addition to SLC24A5 and SLC45A2, genes previously associated with European/non-European differences in skin pigmentation, OPRM1 and EGFR were associated with variation in skin pigmentation in New World populations for the first time.