Chaoying Cui

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.

Genetic evidence of Paleolithic colonization and Neolithic expansion of modern humans on the Tibetan Plateau

Tibetans live on the highest plateau in the world, their current population size is approximately 5 million, and most of them live at an altitude exceeding 3,500 m. Therefore, the Tibetan Plateau is a remarkable area for cultural and biological studies of human population history. However, the chronological profile of the Tibetan Plateaus colonization remains an unsolved question of human prehistory. To reconstruct the prehistoric colonization and demographic history of modern humans on the Tibetan Plateau, we systematically sampled 6,109 Tibetan individuals from 41 geographic populations across the entire region of the Tibetan Plateau and analyzed the phylogeographic patterns of both paternal (n = 2,354) and maternal (n = 6,109) lineages as well as genome-wide single nucleotide polymorphism markers (n = 50) in Tibetan populations. We found that there have been two distinct, major prehistoric migrations of modern humans into the Tibetan Plateau. The first migration was marked by ancient Tibetan genetic signatures dated to approximately 30,000 years ago, indicating that the initial peopling of the Tibetan Plateau by modern humans occurred during the Upper Paleolithic rather than Neolithic. We also found evidences for relatively young (only 7-10 thousand years old) shared Y chromosome and mitochondrial DNA haplotypes between Tibetans and Han Chinese, suggesting a second wave of migration during the early Neolithic. Collectively, the genetic data indicate that Tibetans have been adapted to a high altitude environment since initial colonization of the Tibetan Plateau in the early Upper Paleolithic, before the last glacial maximum, followed by a rapid population expansion that coincided with the establishment of farming and yak pastoralism on the Plateau in the early Neolithic.

Identification of a Tibetan-specific mutation in the hypoxic gene EGLN1 and its contribution to high-altitude adaptation

Tibetans are well adapted to high-altitude hypoxic conditions, and in recent genome-wide scans, many candidate genes have been reported involved in the physiological response to hypoxic conditions. However, the limited sequence variations analyzed in previous studies would not be sufficient to identify causal mutations. Here we conducted resequencing of the entire genomic region (59.4 kb) of the hypoxic gene EGLN1 (one of the top candidates from the genome-wide scans) in Tibetans and identified 185 sequence variations, including 13 novel variations (12 substitutions and 1 insertion or deletion). There is a nonsynonymous mutation (rs186996510, D4E) showing surprisingly deep divergence between Tibetans and lowlander populations (Fst = 0.709 between Tibetans and Han Chinese). It is highly prevalent in Tibetans (70.9% on average) but extremely rare in Han Chinese, Japanese, Europeans, and Africans (0.56-2.27%), suggesting that it might be the causal mutation of EGLN1 contributing to high-altitude hypoxic adaptation. Neutrality test confirmed the signal of Darwinian positive selection on EGLN1 in Tibetans. Haplotype network analysis revealed a Tibetan-specific haplotype, which is absent in other world populations. The estimated selective intensity (0.029 for the C allele of rs186996510) puts EGLN1 among the known genes that have undergone the strongest selection in human populations, and the onset of selection was estimated to have started at the early Neolithic (∼8,400 years ago). Finally, we detected a significant association between rs186996510 and hemoglobin levels in Tibetans, suggesting that EGLN1 contributes to the adaptively low hemoglobin level of Tibetans compared with acclimatized lowlanders at high altitude.

Genetic evidence of a recent Tibetan ancestry to Sherpas in the Himalayan region

Sherpas living around the Himalayas are renowned as high-altitude mountain climbers but when and where the Sherpa people originated from remains contentious. In this study, we collected DNA samples from 582 Sherpas living in Nepal and Tibet Autonomous Region of China to study the genetic diversity of both their maternal (mitochondrial DNA) and paternal (Y chromosome) lineages. Analysis showed that Sherpas share most of their paternal and maternal lineages with indigenous Tibetans, representing a recently derived sub-lineage. The estimated ages of two Sherpa-specific mtDNA sub-haplogroups (C4a3b1 and A15c1) indicate a shallow genetic divergence between Sherpas and Tibetans less than 1,500 years ago. These findings reject the previous theory that Sherpa and Han Chinese served as dual ancestral populations of Tibetans, and conversely suggest that Tibetans are the ancestral populations of the Sherpas, whose adaptive traits for high altitude were recently inherited from their ancestors in Tibet.

Down-Regulation of EPAS1 Transcription and Genetic Adaptation of Tibetans to High-Altitude Hypoxia

Abstract Tibetans are well adapted to the hypoxic environments at high altitude, yet the molecular mechanism of this adaptation remains elusive. We reported comprehensive genetic and functional analyses of EPAS1, a gene encoding hypoxia inducible factor 2α (HIF-2α) with the strongest signal of selection in previous genome-wide scans of Tibetans. We showed that the Tibetan-enriched EPAS1 variants down-regulate expression in human umbilical endothelial cells and placentas. Heterozygous EPAS1 knockout mice display blunted physiological responses to chronic hypoxia, mirroring the situation in Tibetans. Furthermore, we found that the Tibetan version of EPAS1 is not only associated with the relatively low hemoglobin level as a polycythemia protectant, but also is associated with a low pulmonary vasoconstriction response in Tibetans. We propose that the down-regulation of EPAS1 contributes to the molecular basis of Tibetans’ adaption to high-altitude hypoxia.

Physiological Responses and Evaluation of Effects of BMI, Smoking and Drinking in High Altitude Acclimatization: A Cohort Study in Chinese Han Young Males

High altitude acclimatization is a series of physiological responses taking places when subjects go to altitude. Many factors could influence these processes, such as altitude, ascending speed and individual characteristics. In this study, based on a repeated measurement design of three sequential measurements at baseline, acute phase and chronic phase, we evaluated the effect of BMI, smoking and drinking on a number of physiological responses in high altitude acclimatization by using mixed model and partial least square path model on a sample of 755 Han Chinese young males. We found that subjects with higher BMI responses were reluctant to hypoxia. The effect of smoking was not significant at acute phase. But at chronic phase, red blood cell volume increased less while respiratory function increased more for smoking subjects compared with nonsmokers. For drinking subjects, red blood cell volume increased less than nondrinkers at both acute and chronic phases, while blood pressures increased more than nondrinkers at acute phase and respiratory function, red blood cell volume and oxygen saturation increased more than nondrinkers at chronic phase. The heavy and long-term effect of smoking, drinking and other factors in high altitude acclimatization needed to be further studied.

HMOX2 Functions as a Modifier Gene for High-Altitude Adaptation in Tibetans.

Tibetans are well adapted to high‐altitude environments. Among the adaptive traits in Tibetans, the relatively low hemoglobin level is considered a blunted erythropoietic response to hypoxic challenge. Previously, EPAS1 and EGLN1, the major upstream regulators in the hypoxic pathway, were reportedly involved in the hemoglobin regulation in Tibetans. In this study, we report a downstream gene (HMOX2) involved in heme catabolism, which harbors potentially adaptive variants in Tibetans. We first resequenced the entire genomic region (45.6 kb) of HMOX2 in Tibetans, which confirmed the previously suspected signal of positive selection on HMOX2 in Tibetans. Subsequent association analyses of hemoglobin levels in two independent Tibetan populations (a total of 1,250 individuals) showed a male‐specific association between the HMOX2 variants and hemoglobin levels. Tibetan males with the derived C allele at rs4786504:T>C displayed lower hemoglobin level as compared with the T allele carriers. Furthermore, our in vitro experiments indicated that the C allele of rs4786504 could increase the expression of HMOX2, presumably leading to a more efficient breakdown of heme that may help maintain a relatively low hemoglobin level at high altitude. Collectively, we propose that HMOX2 contributes to high‐altitude adaptation in Tibetans by functioning as a modifier in the regulation of hemoglobin metabolism.

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 2 = 0.997, P = 0.002), rs7589621-G (R 2 = 0.994, P = 0.003), rs1868092-A (R 2 = 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.

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.

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.