Genetic variations in Tibetan populations and high-altitude adaptation at the Himalayas

Modern humans have occupied almost all possible environments globally since exiting Africa about 100,000 years ago. Both behavioral and biological adaptations have contributed to their success in surviving the rigors of climatic extremes, including cold, strong ultraviolet radiation, and high altitude. Among these environmental stresses, high-altitude hypoxia is the only condition in which traditional technology is incapable of mediating its effects. Inhabiting at >3,000-m high plateau, the Tibetan population provides a widely studied example of high-altitude adaptation. Yet, the genetic mechanisms underpinning long-term survival in this environmental extreme remain unknown. We performed an analysis of genome-wide sequence variations in Tibetans. In combination with the reported data, we identified strong signals of selective sweep in two hypoxia-related genes, EPAS1 and EGLN1. For these two genes, Tibetans show unusually high divergence from the non-Tibetan lowlanders (Han Chinese and Japanese) and possess high frequencies of many linked sequence variations as reflected by the Tibetan-specific haplotypes. Further analysis in seven Tibetan populations (1,334 individuals) indicates the prevalence of selective sweep across the Himalayan region. The observed indicators of natural selection on EPAS1 and EGLN1 suggest that during the long-term occupation of high-altitude areas, the functional sequence variations for acquiring biological adaptation to high-altitude hypoxia have been enriched in Tibetan populations.     

Genomic Analysis of Natural Selection and Phenotypic Variation in High-Altitude Mongolians

Deedu (DU) Mongolians, who migrated from the Mongolian steppes to the Qinghai-Tibetan Plateau approximately 500 years ago, are challenged by environmental conditions similar to native Tibetan highlanders. Identification of adaptive genetic factors in this population could provide insight into coordinated physiological responses to this environment. Here we examine genomic and phenotypic variation in this unique population and present the first complete analysis of a Mongolian whole-genome sequence. High-density SNP array data demonstrate that DU Mongolians share genetic ancestry with other Mongolian as well as Tibetan populations, specifically in genomic regions related with adaptation to high altitude. Several selection candidate genes identified in DU Mongolians are shared with other Asian groups (e.g., EDAR), neighboring Tibetan populations (including high-altitude candidates EPAS1, PKLR, and CYP2E1), as well as genes previously hypothesized to be associated with metabolic adaptation (e.g., PPARG). Hemoglobin concentration, a trait associated with high-altitude adaptation in Tibetans, is at an intermediate level in DU Mongolians compared to Tibetans and Han Chinese at comparable altitude. Whole-genome sequence from a DU Mongolian (Tianjiao1) shows that about 2% of the genomic variants, including more than 300 protein-coding changes, are specific to this individual. Our analyses of DU Mongolians and the first Mongolian genome provide valuable insight into genetic adaptation to extreme environments.     

Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.     

Genetic determinants of Tibetan high-altitude adaptation

Some highland populations have genetic adaptations that enable their successful existence in a hypoxic environment. Tibetans are protected against many of the harmful responses exhibited by non-adapted populations upon exposure to severe hypoxia, including elevated hemoglobin concentration (i.e., polycythemia). Recent studies have highlighted several genes subject to natural selection in native high-altitude Tibetans. Three of these genes, EPAS1, EGLN1 and PPARA, regulate or are regulated by hypoxia inducible factor, a principal controller of erythropoiesis and other organismal functions. Uncovering the molecular basis of hypoxic adaptation should have implications for understanding hematological and other adaptations involved in hypoxia tolerance. Because the hypoxia response involves a variety of cardiovascular, pulmonary and metabolic functions, this knowledge would improve our understanding of disease mechanisms and could ultimately be translated into targeted therapies for oxygen deprivation, cardiopulmonary and cerebral pathologies, and metabolic disorders such as diabetes and obesity.     

Mitochondrial nt3010G-nt3970C haplotype is implicated in high-altitude adaptation of Tibetans

Tibetans are well adapted to living and thriving in high-altitude environments. Mitochondria are central links to oxygen consumption, and variations in mitochondrial DNA (mtDNA) could play a role in high-altitude adaptation. Alleles at several polymorphic sites in mtDNA define common haplotypes, or haplogroups, including some that have been implicated in the risk of developing certain diseases. However, few reports have determined whether relationships exist between haplogroups and high-altitude adaptation in the Tibetan population. The D4 haplogroup is a major haplogroup of the Han Chinese. In the present study, genotypes of 12 polymorphisms were determined in members of a Tibetan population (n = 72), low altitude-Han (la-Han, n = 144), and high altitude-Han (ha-Han, n = 227) populations using polymerase chain reaction-restriction fragment length polymorphism and polymerase chain reaction-ligase detection reaction assays. The mitochondrial haplogroup D4 was negatively associated with high-altitude adaptation in Tibetans (P = 0.001 vs. la-Han, OR = 0.166, 95% CI = 0.048-0.567; P = 0.009 vs. ha-Han OR = 0.232, 95% CI = 0.069-0.778). The frequency of the nt3010G-nt3970C haplotype was significantly higher in Tibetans than in la-Han (P = 0.000) and ha-Han (P = 0.001) subjects. Findings in the present study suggest that unique mitochondrial variations determine a genetic background that is associated with high-altitude adaptation in the Tibetan population.     

GCH1 plays a role in the high-altitude adaptation of Tibetans

Tibetans are well adapted to high-altitude hypoxia.Previous genome-wide scans have reported many candidate genes for this adaptation,but only a few have been studied.Here we report on a hypoxia gene (GCH1,GTP-cyclohydrolase I),involved in maintaining nitric oxide synthetase (NOS) function and normal blood pressure,that harbors many potentially adaptive variants in Tibetans.We resequenced an 80.8 kb fragment covering the entire gene region of GCH1 in 50 unrelated Tibetans.Combined with previously published data,we demonstrated many GCH1 variants showing deep divergence between highlander Tibetans and lowlander Han Chinese.Neutrality tests confirmed a signal of positive Darwinian selection on GCH1 in Tibetans.Moreover,association analysis indicated that the Tibetan version of GCH1 was significantly associated with multiple physiological traits in Tibetans,including blood nitric oxide concentration,blood oxygen saturation,and hemoglobin concentration.Taken together,we propose that GCH1 plays a role in the genetic adaptation of Tibetans to high altitude hypoxia.     

Genetic adaptation to high altitude in the Ethiopian highlands

Background

Genomic analysis of high-altitude populations residing in the Andes and Tibet has revealed several candidate loci for involvement in high-altitude adaptation, a subset of which have also been shown to be associated with hemoglobin levels, including EPAS1, EGLN1, and PPARA, which play a role in the HIF-1 pathway. Here, we have extended this work to high- and low-altitude populations living in Ethiopia, for which we have measured hemoglobin levels. We genotyped the Illumina 1M SNP array and employed several genome-wide scans for selection and targeted association with hemoglobin levels to identify genes that play a role in adaptation to high altitude.

Results

We have identified a set of candidate genes for positive selection in our high-altitude population sample, demonstrated significantly different hemoglobin levels between high- and low-altitude Ethiopians and have identified a subset of candidate genes for selection, several of which also show suggestive associations with hemoglobin levels.

Conclusions

We highlight several candidate genes for involvement in high-altitude adaptation in Ethiopia, including CBARA1, VAV3, ARNT2 and THRB. Although most of these genes have not been identified in previous studies of high-altitude Tibetan or Andean population samples, two of these genes (THRB and ARNT2) play a role in the HIF-1 pathway, a pathway implicated in previous work reported in Tibetan and Andean studies. These combined results suggest that adaptation to high altitude arose independently due to convergent evolution in high-altitude Amhara populations in Ethiopia.     

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.     

A preliminary study of copy number variation in tibetans

Genetic features of Tibetans have been broadly investigated, but the properties of copy number variation (CNV) have not been well examined. To get a preliminary view of CNV in Tibetans, we scanned 29 Tibetan genomes with the Illumina Human-1 M high-resolution genotyping microarray and identified 139 putative copy number variable regions (CNVRs), consisting of 70 deletions, 61 duplications, and 8 multi-allelic loci. Thirty-four of the 139 CNVRs showed differential allele frequencies versus other East-Asian populations, with P values <0.0001. These results indicated a distinct pattern of CNVR allele frequency distribution in Tibetans. The Tibetan CNVRs are enriched for genes in the disease class of human reproduction (such as genes from the DAZ, BPY2, CDY, and HLA-DQ and -DR gene clusters) and biological process categories of "response to DNA damage stimulus" and "DNA repair" (such as RAD51, RAD52, and MRE11A). These genes are related to the adaptive traits of high infant birth weight and darker skin tone of Tibetans, and may be attributed to recent local adaptation. Our results provide a different view of genetic diversity in Tibetans and new insights into their high-altitude adaptation.     

Population history and genomic signatures for high-altitude adaptation in Tibetan pigs

Background

The Tibetan pig is one of domestic animals indigenous to the Qinghai-Tibet Plateau. Several geographically isolated pig populations are distributed throughout the Plateau. It remained an open question if these populations have experienced different demographic histories and have evolved independent adaptive loci for the harsh environment of the Plateau. To address these questions, we herein investigated ~ 40,000 genetic variants across the pig genome in a broad panel of 678 individuals from 5 Tibetan geographic populations and 34 lowland breeds.

Results

Using a series of population genetic analyses, we show that Tibetan pig populations have marked genetic differentiations. Tibetan pigs appear to be 3 independent populations corresponding to the Tibetan, Gansu and Sichuan & Yunnan locations. Each population is more genetically similar to its geographic neighbors than to any of the other Tibetan populations. By applying a locus-specific branch length test, we identified both population-specific and -shared candidate genes under selection in Tibetan pigs. These genes, such as PLA2G12A, RGCC, C9ORF3, GRIN2B, GRID1 and EPAS1, are involved in high-altitude physiology including angiogenesis, pulmonary hypertension, oxygen intake, defense response and erythropoiesis. A majority of these genes have not been implicated in previous studies of highlanders and high-altitude animals.

Conclusion

Tibetan pig populations have experienced substantial genetic differentiation. Historically, Tibetan pigs likely had admixture with neighboring lowland breeds. During the long history of colonization in the Plateau, Tibetan pigs have developed a complex biological adaptation mechanism that could be different from that of Tibetans and other animals. Different Tibetan pig populations appear to have both distinct and convergent adaptive loci for the harsh environment of the Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-834) contains supplementary material, which is available to authorized users.     

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.     

Evolutionary significance of selected EDAR variants in Tibetan high-altitude adaptations

Humans have been exposed to many environmental challenges since their evolutionary origins in Africa and subsequent migrations to the rest of the world. A severe environmental challenge to human migrants was hypoxia caused by low barometric oxygen pressure at high altitudes. Several genome-wide scans have elucidated the genetic basis of human high-altitude adaptations.However, the dearth of functional variant information has led to the successful association of only a few candidate genes. In the present study, we employed a candidate gene approach and re-sequenced the EDAR locus in 45 Tibetan individuals to identify mutations involved in hypoxia adaptation. We identified 10 and five quantitative trait-associated mutations for oxygen saturation (SaO_2) and blood platelet count, respectively, at the EDAR locus. Among these, rs10865026 and rs3749110 (associated with SaO_2 and platelet count, respectively) were identified as functional candidate targets. These data demonstrate that EDAR has undergone natural selection in recent human history and indicate an important role of EDAR variants in Tibetan high-altitude adaptations.     

The EPAS1 gene influences the aerobic–anaerobic contribution in elite endurance athletes

Absract EPAS1 is a gene involved in complex oxygen sensing. It is expressed in microvascular endothelial cells, lung epithelial cells, cardiac myocytes and the brain. An association study was undertaken comparing elite endurance athletes classified into two groups according to a power–time model of performance intensity: power–time-maximum (PT-MAX; N=242, event duration 50 s to 10 min) and power–time–steady state (PT-SS; N=151, event duration ~2–10 h), with normal controls (N=444) using 12 SNPs across EPAS1. Ordinal regression analysis of allele frequencies revealed significant differences at SNPs 2 and 3 (P=0.01). Haplotype analysis revealed the presence of haplotypes involving SNPs 2–5 that significantly differentiated (P<0.05) the groups based on an ordinal ranking using the power–time classification. These same haplotypes differentiated the PT-MAX group in which a significant decrease in a haplotype (F: G-C-C-G; OR=0.57, P=0.02, 95% CI 0.36–0.92) and increase in a second haplotype (G: A-T-G-G; OR=1.75, P=0.03, 95% CI 1.05–2.91) was observed compared to controls. The PT-SS group was differentiated from the PT-MAX group by a third haplotype (H: A-T-G-A; OR=0.46, P=0.04, 95% CI 0.22–0.96). Since EPAS1 has a role as a sensor capable of integrating cardiovascular function, energetic demand, muscle activity and oxygen availability into physiological adaptation, we propose that DNA variants in EPAS1 influence the relative contribution of aerobic and anaerobic metabolism and hence the maximum sustainable metabolic power for a given event duration.     

Association between MT-CO3 haplotypes and high-altitude adaptation in Tibetan chicken

Genetic mutation in cytochrome c oxidase subunit III gene (MT-CO3) could influence the kinetics of cytochrome c oxidase (COX), which catalyzes oxygen transport capacity in oxidative phosphorylation. However, the potential relationship between MT-CO3 variants and high-altitude adaptation remains poorly understood in Tibetan chicken. Here, we sequenced MT-CO3 gene of 125 Tibetan chickens and 144 Chinese domestic chickens in areas at a low elevation (below 1000 m). Eight single nucleotide polymorphisms (SNPs) were detected; and five of them (m.10081A>G, m.10115G>A, m.10270G>A, m.10336A>G and m.10447C>T) shared by Tibetan chicken and lowland chicken with the significant difference in their respective allele frequencies. Nine haplotypes (H1–H9) were finally defined. Among them, haplotype H4 was positively associated with high-altitude adaptation whereas haplotypes H6, H7 and H8 had negative association with high-altitude adaptation. The Median-joining profile suggested that haplotype H5 had the ancestral position to the other haplotypes but had no significant relationship with high-altitude adaptation. However, there was only m.10081A>G mutation differed from haplotype H4 and H5. Results also suggested that chickens with A allele at m.10081A>G, had over 2.6 times than those with G allele in the probability of the ability to adapt hypoxia. It suggests that the synonymous mutation m.10081A>G may be a prerequisite for shaping high-altitude adaptation-specific haplotypes.     

藏族的高原适应——西藏藏族生物人类学研究回顾

藏族生活在具有世界屋脊之称的青藏高原, 特殊的生态环境和特殊的文化背景造就了藏族特殊的适应高原缺氧机制, 引起了国内外学者的广泛关注和浓厚的研究兴趣。本文根据国内外数据库的文献并结合我们的研究工作, 从高原适应的角度回顾了30多年藏族人类学研究。回顾显示, 藏族由于长期生活在高原缺氧的环境中, 不仅形态和机能发生了适应性变化, 而且体成分也表现出相应的变化, 体现了形态、机能和体成分的统一。这些变化是长期进化形成的, 与安第斯山人等有明显不同, 就是在同一高原生活的西藏、青海、四川、甘肃和云南的藏族乃至尼泊尔和印度藏族的体质也表现出地域差异, 这些差异的产生是多种因素所致, 两个关键性的基因是导致两大高原人口高原适应机制不同的最主要的原因。     

HIF2A Variants Were Associated with Different Levels of High-Altitude Hypoxia among Native Tibetans

Hypoxia inducible factors, including HIF1A and HIF2A, play central roles in response to high-altitude hypoxia and genetic variants of HIF1A or HIF2A were associated with high-altitude sickness or adaptation. However, it remains to determine whether they are associated with tolerance to different levels of high-altitude selection pressure among native Tibetans. We recruited 189 Tibetan subjects living at 2,700 meters (Low level of high altitude, LHA), 197 at 3,200 meters (Middle level of high altitude of high altitude, MHA), 249 at 3,700 meters (High level of high altitude, HHA) and 269 at 4,700 meters (Very high level of high altitude, VHA) and performed association analysis of twelve tSNPs (tagging SNPs) in HIF1A and HIF2A with high-altitude. We found (1) a increasing trend of HIF2A rs5621780-C(18.4%, 15.9%, 32.8% and 31.1%, respectively, in LHA, MHA, HHA and VHA)(P = 3.56E-9); (2) increasing trends of HIF2A rs6756667-A(68.7%, 73.4%, 79.9% and 89.6%), rs7589621- G(74.6%, 77.9%, 83.7%, and 92.1%) and rs1868092-A(64.1%, 67.3%, 75.1% and 84.4%) (P = 3.56E-9, 4.68E-16, 1.17E-13 and 7.09E-14, respectively); (3) a increasing trend of haplotype AG (68.7%, 73.1%, 79.9% and 89.6%) (P = 2.22E-7) which was constructed by rs6756667 and rs7589621; (4) a strong linear correlation between major alleles of rs6756667-A (R ² = 0.997, P = 0.002), rs7589621-G (R ² = 0.994, P = 0.003), rs1868092-A (R ² = 0.985, P = 0.008) and altitude by linear correlation test. The associations between HIF2A variants and different level of high altitude support that extremely high-altitude hypoxia challenge imposes selective effects on HIF2A variants among native Tibetans.     

EP300 contributes to high-altitude adaptation in Tibetans by regulating nitric oxide production

The genetic adaptation of Tibetans to high altitude hypoxia likely involves a group of genes in the hypoxic pathway,as suggested by earlier studies.To test the adaptive role of the previously reported candidate gene EP300 (histone acetyltransferase p300),we conducted resequencing of a 108.9 kb gene region of EP300 in 80 unrelated Tibetans.The allele-frequency and haplotype-based neutrality tests detected signals of positive Darwinian selection on EP300 in Tibetans,with a group of variants showing allelic divergence between Tibetans and lowland reference populations,including Han Chinese,Europeans,and Africans.Functional prediction suggested the involvement of multiple EP300 variants in gene expression regulation.More importantly,genetic association tests in 226 Tibetans indicated significant correlation of the adaptive EP300 variants with blood nitric oxide (NO) concentration.Collectively,we propose that EP300 harbors adaptive variants in Tibetans,which might contribute to high-altitude adaptation through regulating NO production.