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

Genetic signatures reveal high-altitude adaptation in a set of ethiopian populations

The Tibetan and Andean Plateaus and Ethiopian highlands are the largest regions to have long-term high-altitude residents. Such populations are exposed to lower barometric pressures and hence atmospheric partial pressures of oxygen. Such “hypobaric hypoxia” may limit physical functional capacity, reproductive health, and even survival. As such, selection of genetic variants advantageous to hypoxic adaptation is likely to have occurred. Identifying signatures of such selection is likely to help understanding of hypoxic adaptive processes. Here, we seek evidence of such positive selection using five Ethiopian populations, three of which are from high-altitude areas in Ethiopia. As these populations may have been recipients of Eurasian gene flow, we correct for this admixture. Using single-nucleotide polymorphism genotype data from multiple populations, we find the strongest signal of selection in BHLHE41 (also known as DEC2 or SHARP1). Remarkably, a major role of this gene is regulation of the same hypoxia response pathway on which selection has most strikingly been observed in both Tibetan and Andean populations. Because it is also an important player in the circadian rhythm pathway, BHLHE41 might also provide insights into the mechanisms underlying the recognized impacts of hypoxia on the circadian clock. These results support the view that Ethiopian, Andean, and Tibetan populations living at high altitude have adapted to hypoxia differently, with convergent evolution affecting different genes from the same pathway.

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.

EGLN1 variants influence expression and SaO2 levels to associate with high-altitude pulmonary oedema and adaptation

EGLN1 [encoding HIF (hypoxia-inducible factor)-prolyl hydroxylase 2] plays a pivotal role in the HIF pathway and has emerged as one of the most intriguing genes with respect to physiology at HA (high altitude). EGLN1, being an actual oxygen sensor, appears to have a potential role in the functional adaptation to the hypobaric hypoxic environment. In the present study, we screened 30 polymorphisms of EGLN1, evaluated its gene expression and performed association analyses. In addition, the role of allelic variants in altering TF (transcription factor)-binding sites and consequently the replacement of TFs at these loci was also investigated. The study was performed in 250 HAPE-p [HAPE (HA pulmonary oedema)-patients], 210 HAPE-f (HAPE-free controls) and 430 HLs (healthy Ladakhi highland natives). The genotypes of seven polymorphisms, rs1538664, rs479200, rs2486729, rs2790879, rs480902, rs2486736 and rs973252, differed significantly between HAPE-p and HAPE-f (P<0.008). The genotypes AA, TT, AA, GG, CC, AA and GG of rs1538664, rs479200, rs2486729, rs2790879, rs480902, rs2486736 and rs973252, prevalent in HAPE-p, were identified as risk genotypes and their counterpart homozygotes, prevalent in HLs, were identified as protective. EGLN1 expression was up-regulated 4.56-fold in HAPE-p (P=0.0084). The risk genotypes, their haplotypes and interacting genotypes were associated with up-regulated EGLN1 expression (P<0.05). Similarly, regression analysis showed that the risk alleles and susceptible haplotypes were associated with decreased SaO2 (arterial oxygen saturation) levels in the three groups. The significant inverse correlation of SaO2 levels with PASP (pulmonary artery systolic pressure) and EGLN1 expression and the association of these polymorphisms with SaO2 levels and EGLN1 expression contributed to uncovering the molecular mechanism underlying hypobaric hypoxic adaptation and maladaptation.

Regulation of erythropoiesis by hypoxia-inducible factors

A classic physiologic response to systemic hypoxia is the increase in red blood cell production. Hypoxia-inducible factors (HIFs) orchestrate this response by inducing cell-type specific gene expression changes that result in increased erythropoietin (EPO) production in kidney and liver, in enhanced iron uptake and utilization and in adjustments of the bone marrow microenvironment that facilitate erythroid progenitor maturation and proliferation. In particular HIF-2 has emerged as the transcription factor that regulates EPO synthesis in the kidney and liver and plays a critical role in the regulation of intestinal iron uptake. Its key function in the hypoxic regulation of erythropoiesis is underscored by genetic studies in human populations that live at high-altitude and by mutational analysis of patients with familial erythrocytosis. This review provides a perspective on recent insights into HIF-controlled erythropoiesis and iron metabolism, and examines cell types that have EPO-producing capability. Furthermore, the review summarizes clinical syndromes associated with mutations in the O(2)-sensing pathway and the genetic changes that occur in high altitude natives. The therapeutic potential of pharmacologic HIF activation for the treatment of anemia is discussed.

Genetic variants in EPAS1 contribute to adaptation to high-altitude hypoxia in Sherpas

Sherpas comprise a population of Tibetan ancestry in the Himalayan region that is renowned for its mountaineering prowess. The very small amount of available genetic information for Sherpas is insufficient to explain their physiological ability to adapt to high-altitude hypoxia. Recent genetic evidence has indicated that natural selection on the endothelial PAS domain protein 1 (EPAS1) gene was occurred in the Tibetan population during their occupation in the Tibetan Plateau for millennia. Tibetan-specific variations in EPAS1 may regulate the physiological responses to high-altitude hypoxia via a hypoxia-inducible transcription factor pathway. We examined three significant tag single-nucleotide polymorphisms (SNPs, rs13419896, rs4953354, and rs4953388) in the EPAS1 gene in Sherpas, and compared these variants with Tibetan highlanders on the Tibetan Plateau as well as with non-Sherpa lowlanders. We found that Sherpas and Tibetans on the Tibetan Plateau exhibit similar patterns in three EPAS1 significant tag SNPs, but these patterns are the reverse of those in non-Sherpa lowlanders. The three SNPs were in strong linkage in Sherpas, but in weak linkage in non-Sherpas. Importantly, the haplotype structured by the Sherpa-dominant alleles was present in Sherpas but rarely present in non-Sherpas. Surprisingly, the average level of serum erythropoietin in Sherpas at 3440 m was equal to that in non-Sherpas at 1300 m, indicating a resistant response of erythropoietin to high-altitude hypoxia in Sherpas. These observations strongly suggest that EPAS1 is under selection for adaptation to the high-altitude life of Tibetan populations, including Sherpas. Understanding of the mechanism of hypoxia tolerance in Tibetans is expected to provide lights to the therapeutic solutions of some hypoxia-related human diseases, such as cardiovascular disease and cancer.

[The research progress of human genetic diversity in China]

Human genetic diversity refers to genomic variation among races, ethnic groups, isolated populations and individuals worldwide, and is one major resource and tool on discovering human evolution and migration, interaction between genetic background and environment, and factors associated with human diseases and health. China has abundant and valuable resource of human genetic diversity due to 56 ethnic groups and a large population accounting for one fifth of the total population in the world. After decades of efforts, a large number of research data on human genetic diversity have been accumulated in China, and some of outcomes reach advanced international level. This review mainly focuses on the recent progress and outcomes achieved in applying genetic markers including morphological markers, biochemical and immunological markers and DNA markers in research of genetic diversity, and the application of mitochondrial DNA, Y chromosomal DNA, HLA and others in research of the origin and relationship of Chinese ethic groups, and the origin and mi-of modern East Asian populations. This review also summarizes the advances in the research fields of preservation and utilization of Chinese genetic resource, identification of genes associated with disease selective and adaptive for natural pressure, application of whole genome association study and next generation sequencing, and Chinese human genome as well.

Toward understanding the genetic basis of adaptation to high-elevation life in poikilothermic species: a comparative transcriptomic analysis of two ranid frogs, Rana chensinensis and R. kukunoris

BACKGROUND

Understanding how organisms adapt to high-elevation environments at a genome scale provides novel insights into the process of adaptive evolution. Previous studies have mainly focused on endothermic organisms, while poikilothermic species may have evolved different mechanisms to cope with high-elevation environments. In this context, we sequenced transcriptomes of a pair of closely related anuran species, Rana chensinensis and R. kukunoris, which inhabit respective low- and high-elevation habitats. By comparing the two transcriptomes, we identified candidate genes that may be involved in high-elevation adaption in poikilothermic species.

RESULTS

Over 66 million sequence reads from each transcriptome were generated. A total of 41,858 and 39,293 transcripts for each species were obtained by de novo assembly. By comparing the orthologous transcripts, we identified 125 protein-coding genes that have likely experienced strong positive selection (Ka/Ks>1). In addition, 335 genes that may bear a signature of positive selection (1≥Ka/Ks>0.5) were also recognized. By considering their functions, fourteen candidate genes were determined to be likely involved in high-elevation adaptation, including two CYP genes, USP-1, and several others.

CONCLUSIONS

We identified a set of candidate genes that may have promoted adaptation of R. kukunoris to its high-elevation environment. These include several genes that have previously been associated with oxygen binding, response to UV radiation, and repair of free radical injury. Detailed molecular, physiological, and phenotypic studies are required to better understand the roles of these genes in improving the performance of R. kukunoris in a high-elevation environment. We have evidence for both convergent and divergent evolution between endothermic and poikilothemic species, but additional research across a wide range of organisms will be necessary to comprehend the complexity of high-elevation adaptation.

Oxygen levels and the regulation of cell adhesion in the nervous system: a control point for morphogenesis in development, disease and evolution?

In this article, I discuss the hallmarks of hypoxia in vitro and in vivo and review work showing that many types of stem cell proliferate more robustly in lowered oxygen. I then discuss recent studies showing that alterations in the levels and the types of cell and substrate adhesion molecules are a notable response to reduced O(2) levels in both cultured primary neural stem cells and brain tissues in response to hypoxia in vivo. The ability of O(2) levels to regulate adhesion molecule expression is linked to the Wnt signaling pathway, which can control and be controlled by adhesion events. The ability of O(2) levels to influence cell adhesion also has far-reaching implications for development, ischemic trauma and neural regeneration, as well as for cancer and other diseases. Finally I discuss the possibility that the fluctuations in O(2) levels known to have occurred over evolutionary time could, by influencing adhesion systems, have contributed to early symbiotic events in unicellular organisms and to the emergence of multicellularity. It is not my intention to be exhaustive in these domains, which are far from my own field of study. Rather this article is meant to provoke and stimulate thinking about molecular evolution involving O(2) sensing and signaling during eras of geologic and atmospheric change that might inform modern studies on development and disease.

Genetic and phenotypic divergence between low- and high-altitude populations of two recently diverged cinnamon teal subspecies

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species' range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low- and high-altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P(ST) = 0.56) and exhibited significant frequency differences in a single nonsynonymous α-hemoglobin amino acid polymorphism (Asn/Ser-α9; F(ST) = 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F(ST) = 0.004-0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high-altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.

Metabolic insight into mechanisms of high-altitude adaptation in Tibetans

Recent studies have identified genes involved in high-altitude adaptation in Tibetans. Genetic variants/haplotypes within regions containing three of these genes (EPAS1, EGLN1, and PPARA) are associated with relatively decreased hemoglobin levels observed in Tibetans at high altitude, providing corroborative evidence for genetic adaptation to this extreme environment. The mechanisms that afford adaptation to high-altitude hypoxia, however, remain unclear. Considering the strong metabolic demands imposed by hypoxia, we hypothesized that a shift in fuel preference to glucose oxidation and glycolysis at the expense of fatty acid oxidation would improve adaptation to decreased oxygen availability. Correlations between serum free fatty acid and lactate concentrations in Tibetan groups living at high altitude and putatively selected haplotypes provide insight into this hypothesis. An EPAS1 haplotype that exhibits a signal of positive selection is significantly associated with increased lactate concentration, the product of anaerobic glycolysis. Furthermore, the putatively advantageous PPARA haplotype is correlated with serum free fatty acid concentrations, suggesting a possible decrease in the activity of fatty acid oxidation. Although further studies are required to assess the molecular mechanisms underlying these patterns, these associations suggest that genetic adaptation to high altitude involves alteration in energy utilization pathways.

Genetic adaptation of the hypoxia-inducible factor pathway to oxygen pressure among eurasian human populations

Research into the mechanisms of human adaptation to the hypoxic environment of high altitude is of great interest to the fields of human physiology and clinical medicine. Recently, the gene EGLN1, from the hypoxia-inducible factor (HIF) pathway, was identified as being involved in the hypoxic adaptation of highland Andeans and Tibetans. Both highland Andeans and Tibetans have adapted to an extremely hypoxic habitat and less attention has been paid to populations living in normoxic conditions at sea level and mild-hypoxic environments of moderate altitude, thus, whether a common adaptive mechanism exists in response to quantitative variations of environmental oxygen pressure over a wide range of residing altitudes is unknown. Here, we first performed a genome-wide association study of 35 populations from the Human Genome Diversity-CEPH Panel who dwell at sea level to moderate altitude in Eurasia (N = 691; 0-2,500 m) to identify the genetic adaptation profile of normoxic and mild-hypoxic inhabitants. In addition, we systematically compared the results from the present study to six previously published genome-wide scans of highland Andeans and Tibetans to identify shared adaptive signals in response to quantitative variations of oxygen pressure. For normoxic and mild-hypoxic populations, the strongest adaptive signal came from the mu opioid receptor-encoding gene (OPRM1, 2.54 × 10(-9)), which has been implicated in the stimulation of respiration, while in the systematic survey the EGLN1-DISC1 locus was identified in all studies. A replication study performed with highland Tibetans (N = 733) and sea level Han Chinese (N = 748) confirmed the association between altitude and SNP allele frequencies in OPRM1 (in Tibetans only, P < 0.01) and in EGLN1-DISC1 (in Tibetans and Han Chinese, P < 0.01). Taken together, identification of the OPRM1 gene suggests that cardiopulmonary adaptation mechanisms are important and should be a focus in future studies of hypoxia adaptation. Furthermore, the identification of the EGLN1 gene from the HIF pathway suggests a common adaptive mechanism for Eurasian human populations residing at different altitudes with different oxygen pressures.

EPAS1 and EGLN1 associations with high altitude sickness in Han and Tibetan Chinese at the Qinghai-Tibetan Plateau

High altitude sickness (HAS) occurs among humans visiting or inhabiting high altitude environments. Genetic differences in the EPAS1 and EGLN1 genes have been found between lowland (Han) and highland (Tibetan) Chinese. Three SNPs within EPAS1 and EGLN1 were evaluated in Han and Tibetan patients with acute mountain sickness (AMS) and chronic mountain sickness (CMS). We compared 85 patients with AMS to 79 Han unaffected with mountain sickness (MS) as well as 45 CMS patients to 34 unaffected Tibetan subjects. The three SNPs studied were EPAS1 [ch2: 46441523 (hg18], EGLN1 (rs480902) and (rs516651). Direct sequencing was used to identify individual genotypes for the three SNPs. Age was found to be significantly associated with the EPAS1 SNP in the CMS patients while heart rate (HR) and oxygen saturation level of hemoglobin (SaO(2)) were found to be significantly associated with the EGLN1 (rs480902) SNP in the Han patients with AMS. The individuals with CMS were found to diverge significantly for the EPAS1 SNP compared to their Tibetan control group as measured by genetic distance (0.123) indicating positive selection of the EPAS-G allele with age and illness. The EGLN1 (rs480902) SNP had a significant correlation with hematocrit (HCT), HR and SaO(2) in AMS patients. AMS and CMS were found to be significantly associated with the EPAS1 and EGLN1 SNPs compared to their Han and Tibetan control groups, respectively, indicating these nucleotide alterations have a physiological effect for the development of high altitude sickness.

Transcriptome and network changes in climbers at extreme altitudes

Extreme altitude can induce a range of cellular and systemic responses. Although it is known that hypoxia underlies the major changes and that the physiological responses include hemodynamic changes and erythropoiesis, the molecular mechanisms and signaling pathways mediating such changes are largely unknown. To obtain a more complete picture of the transcriptional regulatory landscape and networks involved in extreme altitude response, we followed four climbers on an expedition up Mount Xixiabangma (8,012 m), and collected blood samples at four stages during the climb for mRNA and miRNA expression assays. By analyzing dynamic changes of gene networks in response to extreme altitudes, we uncovered a highly modular network with 7 modules of various functions that changed in response to extreme altitudes. The erythrocyte differentiation module is the most prominently up-regulated, reflecting increased erythrocyte differentiation from hematopoietic stem cells, probably at the expense of differentiation into other cell lineages. These changes are accompanied by coordinated down-regulation of general translation. Network topology and flow analyses also uncovered regulators known to modulate hypoxia responses and erythrocyte development, as well as unknown regulators, such as the OCT4 gene, an important regulator in stem cells and assumed to only function in stem cells. We predicted computationally and validated experimentally that increased OCT4 expression at extreme altitude can directly elevate the expression of hemoglobin genes. Our approach established a new framework for analyzing the transcriptional regulatory network from a very limited number of samples.

Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude

We performed a mitochondrial whole-genome comparison study in 40 Tibetan and 50 Han Chinese. All subjects could be classified into 13 haplogroups pertained to the Macrohaplogroup M and N that pitched different quadrants by principal component analysis. We observed a difference in the M9 haplogroup and identified 18 significant variants by comparing whole sequences between Tibetan and Han populations. Variants in ND2, COX2, tRNA alanine and 12S rRNA were predicted to confer increased protein stability in Tibetans. We compared the base substitutions of nonsynonymous (NS) versus synonymous (S) of 13 protein-encoding genes and found the NS/S values of the ATP6, ATP8, and Cyt b genes were larger (>1) in Tibetans than that in Han population. Our findings provide clues for the existence of adaptive selection for the ATP6, ATP8, Cyt b, ND2, COX2, tRNA alanine and 12S rRNA genes in Tibetans which likely contributed to adaptation to their specific geographic environment, such as high 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.

High altitude adaptation in Daghestani populations from the Caucasus

We have surveyed 15 high-altitude adaptation candidate genes for signals of positive selection in North Caucasian highlanders using targeted re-sequencing. A total of 49 unrelated Daghestani from three ethnic groups (Avars, Kubachians, and Laks) living in ancient villages located at around 2,000 m above sea level were chosen as the study population. Caucasian (Adygei living at sea level, N = 20) and CEU (CEPH Utah residents with ancestry from northern and western Europe; N = 20) were used as controls. Candidate genes were compared with 20 putatively neutral control regions resequenced in the same individuals. The regions of interest were amplified by long-PCR, pooled according to individual, indexed by adding an eight-nucleotide tag, and sequenced using the Illumina GAII platform. 1,066 SNPs were called using false discovery and false negative thresholds of ~6%. The neutral regions provided an empirical null distribution to compare with the candidate genes for signals of selection. Two genes stood out. In Laks, a non-synonymous variant within HIF1A already known to be associated with improvement in oxygen metabolism was rediscovered, and in Kubachians a cluster of 13 SNPs located in a conserved intronic region within EGLN1 showing high population differentiation was found. These variants illustrate both the common pathways of adaptation to high altitude in different populations and features specific to the Daghestani populations, showing how even a mildly hypoxic environment can lead to genetic adaptation.

Matrix metalloproteinase (MMP)-9 induced by Wnt signaling increases the proliferation and migration of embryonic neural stem cells at low O2 levels

Recent studies have shown that various neural and embryonic stem cells cultured in 1-8% oxygen (O(2)), levels lower than those typically used in cell culture (20.9%), displayed increased rates of proliferation; however, the molecular mechanisms underlying these changes are largely undefined. In this study, using rigorously controlled O(2) levels, we found that neural stem cells (NSCs) from embryonic day 15 rat cortex increased their rate of proliferation and migration in 1% O(2) relative to 20% O(2) without changes in viability. We sought to identify molecular changes in NSCs grown in 1% O(2) that might account for these increases. In 1% O(2), levels of the hypoxia-inducible transcription factor HIF-1α were transiently increased. Reduced adherence of NSCs in 1% O(2) to basement membrane-coated plates was observed, and quantitative RT-PCR analysis confirmed that the levels of mRNA for an assortment of cell adhesion and extracellular matrix molecules were altered. Most notable was a 5-fold increase in matrix metalloproteinase (MMP)-9 mRNA. Specific inhibition of MMP-9 activity, verified using a fluorescent substrate assay, prevented the increase in proliferation and migration in 1% O(2). The canonical Wnt pathway was recently shown to be activated in stem cells in low O(2) via HIF-1α. Inhibition of Wnt signaling by DKK-1 also prevented the increase in proliferation, migration, and MMP-9 expression. Thus, MMP-9 is a key molecular effector, downstream of HIF-1α and Wnt activation, responsible for increased rates of NSC proliferation and migration in 1% O(2).