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

Blood lead levels in Peruvian adults are associated with proximity to mining and DNA methylation

Background: Inorganic lead (Pb) is common in the environment, and is toxic to neurological, renal, and cardiovascular systems. Pb exposure influences the epigenome with documented effects on DNA methylation (DNAm). We assessed the impact of low levels of Pb exposure on DNAm among non-miner individuals from two locations in Peru: Lima, the capital, and Cerro de Pasco, a highland mining town, to study the effects of Pb exposure on physiological outcomes and DNAm. Methods: Pb levels were measured in whole blood (n = 305). Blood leukocyte DNAm was determined for 90 DNA samples using the Illumina MethylationEPIC chip. An epigenome-wide association study was performed to assess the relationship between Pb and DNAm. Results: Individuals from Cerro de Pasco had higher Pb than individuals from Lima (p-value = 2.00E-16). Males had higher Pb than females (p-value = 2.36E-04). Pb was positively associated with hemoglobin (p-value = 8.60E-04). In Cerro de Pasco, blood Pb decreased with the distance from the mine (p-value = 0.04), and association with soil Pb was approaching significance (p-value = 0.08). We identified differentially methylated positions (DMPs) associated with genes SOX18, ZMIZ1, and KDM1A linked to neurological function. We also found 45 differentially methylated regions (DMRs), seven of which were associated with genes involved in metal ion binding and nine to neurological function and development. Conclusions: Our results demonstrate that even low levels of Pb can have a significant impact on the body including changes to DNAm. We report associations between Pb and hemoglobin, Pb and distance from mining, and between blood and soil Pb. We also report associations between loci- and region-specific DNAm and Pb

LINE-1 and EPAS1 DNA methylation associations with high-altitude exposure

Recent discoveries indicate a genetic basis for high-altitude adaptation among human groups who have resided at high altitude for millennia, including Andeans, Tibetans, and Ethiopians. Yet, genetics alone does not explain the extent of variation in altitude-adaptive phenotypes. Current and past environments may also play a role, and one way to determine the effect of the environment is through the epigenome. To characterize if Andean adaptive responses to high altitude have an epigenetic component, we analyzed DNA methylation of the promoter region of EPAS1 and LINE-1 repetitive element among 572 Quechua individuals from high- (4,388 m) and low-altitude (0 m) in Peru. Participants recruited at high altitude had lower EPAS1 DNA methylation and higher LINE-1 methylation. Altitude of birth was associated with higher LINE-1 methylation, not with EPAS1 methylation. The number of years lived at high altitude was negatively associated with EPAS1 methylation and positively associated with LINE-1 methylation. We found four one-carbon metabolism SNPs (MTHFD1 rs2236225, TYMS rs502396, FOLH1 rs202676, GLDC rs10975681) that cumulatively explained 11.29% of the variation in average LINE-1 methylation. And identified an association between LINE-1 methylation and genome-wide SNP principal component 1 that distinguishes European from Indigenous American ancestry suggesting that European admixture decreases LINE-1 methylation. Our results indicate that both current and lifetime exposure to high-altitude hypoxia have an effect on EPAS1 and LINE-1 methylation among Andean Quechua, suggesting that epigenetic modifications may play a role in high-altitude adaptation

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

Angiotensin-Converting Enzyme Genotype and Arterial Oxygen Saturation at High Altitude in Peruvian Quechua

Bigham, Abigail W., Melisa Kiyamu, Fabiola León-Verlarde, Esteban J. Parra, Maria Rivera-Ch, Mark D. Shriver, and Tom D. Brutsaert. Angiotensin-converting enzyme genotype and arterial oxygen saturation at high altitude in Peruvian Quechua. High Alt. Med. Biol. 9:167–178, 2008.—The I-allele of the angiotensin-converting enzyme (ACE) gene insertion/deletion (I/D) polymorphism has been associated with performance benefits at high altitude (HA). In n = 142 young males and females of largely Quechua origins in Peru, we evaluated 3 specific hypotheses with regard to the HA benefits of the I-allele: (1) the I-allele is associated with higher arterial oxygen saturation (\documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document}) at HA, (2) the I-allele effect depends on the acclimatization state of the subjects, and (3) the putative I-allele effect on \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} is mediated by the isocapnic hypoxic ventilatory response (HVR, \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$$1 / {\rm min}^{- 1} / \%{\rm Sa}_{\rm O_2}{- 1}$$\end{document}). The subject participants comprised two different study groups including BLA subjects (born at low altitude) who were lifelong sea-level residents transiently exposed to hypobaric hypoxia (<24 h) and BHA subjects (born at HA) who were lifelong residents of HA. To control for the possibility of population stratification, Native American ancestry proportion (NAAP) was estimated as a covariate for each individual using a panel of 70 ancestry-informative molecular markers (AIMS). At HA, resting and exercise \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} was strongly associated with the ACE genotype, p = 0.008 with ∼4% of the total variance in \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} attributed to ACE genotype. Moreover, I/I individuals maintained ∼2.3 percentage point higher \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} compared to I/D and D/D. This I-allele effect was evident in both BLA and BHA groups, suggesting that acclimatization state has little influence on the phenotypic expression of the ACE gene. Finally, ACE genotype was not associated with the isocapnic HVR, although HVR had a strong independent effect on \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} (p = 0.001). This suggests that the I-allele effect on \documentclass{aastex}\usepackage{amsbsy}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{bm}\usepackage{mathrsfs}\usepackage{pifont}\usepackage{stmaryrd}\usepackage{textcomp}\usepackage{portland,xspace}\usepackage{amsmath,amsxtra}\pagestyle{empty}\DeclareMathSizes{10}{9}{7}{6}\begin{document}$${\rm Sa}_{\rm O_2}$$\end{document} is not mediated by the peripheral control of breathing, but rather by some other central cardiopulmonary effect of the ACE gene on the renin–angiotensin–aldosterone system (RAAS)

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