Yong-Bo Guo, Yao-Xi He, Chao-Ying Cui, Ouzhuluobu, Baimakangzhuo, Duojizhuoma, Dejiquzong, Bianba, Yi Peng, Cai-juan Bai, Gonggalanzi, Yong-Yue Pan, Qula, Kangmin, Cirenyangji, Baimayangji, Wei Guo, Yangla, Hui Zhang, Xiao-Ming Zhang, Wang-Shan Zheng,, Shu-Hua Xu, Hua Chen, Sheng-Guo Zhao, Yuan Cai, Shi-Ming Liu, Tian-Yi Wu, Xue-Bin Qi,*, Bing Su,*

1College of Animal Science and Technology,Gansu Agricultural University,Lanzhou Gansu730070,China

2State Key Laboratory of Genetic Resources and Evolution,Kunming Institute of Zoology,Chinese Academy of Sciences,Kunming Yunnan650223,China

3High Altitude Medical Research Center,School of Medicine,Tibetan University,Lhasa Tibet850000,China

4Kunming College of Life Science,University of Chinese Academy of Sciences,Kunming Yunnan 650204,China

5Chinese Academy of Sciences Key Laboratory of Computational Biology,Max Planck Independent Research Group on Population Genomics, CAS-MPG Partner Institute for Computational Biology (PICB),Shanghai Institutes for Biological Sciences,Chinese Academy of Sciences, Shanghai200031,China

6Center for Computational Genomics,Beijing Institute of Genomics,Chinese Academy of Sciences,Beijing100101,China

7National Key Laboratory of High Altitude Medicine,High Altitude Medical Research Institute,XiningQinghai810012,China

8School of Life Science and Technology,Shanghai Tech University,Shanghai200031,China

9Collaborative Innovation Center of Genetics and Development,Shanghai200438,China

ZOOLOGICAL RESEARCH

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

Yong-Bo Guo1,2,#, Yao-Xi He2,4,#, Chao-Ying Cui3,#, Ouzhuluobu3, Baimakangzhuo3, Duojizhuoma3, Dejiquzong3, Bianba3, Yi Peng2, Cai-juan Bai3, Gonggalanzi3, Yong-Yue Pan3, Qula3, Kangmin3, Cirenyangji3, Baimayangji3, Wei Guo3, Yangla3, Hui Zhang2, Xiao-Ming Zhang2, Wang-Shan Zheng1,2, Shu-Hua Xu5,8,9, Hua Chen6, Sheng-Guo Zhao1, Yuan Cai1, Shi-Ming Liu7, Tian-Yi Wu7, Xue-Bin Qi2,*, Bing Su2,*

1College of Animal Science and Technology,Gansu Agricultural University,Lanzhou Gansu730070,China

2State Key Laboratory of Genetic Resources and Evolution,Kunming Institute of Zoology,Chinese Academy of Sciences,Kunming Yunnan650223,China

3High Altitude Medical Research Center,School of Medicine,Tibetan University,Lhasa Tibet850000,China

4Kunming College of Life Science,University of Chinese Academy of Sciences,Kunming Yunnan 650204,China

5Chinese Academy of Sciences Key Laboratory of Computational Biology,Max Planck Independent Research Group on Population Genomics, CAS-MPG Partner Institute for Computational Biology (PICB),Shanghai Institutes for Biological Sciences,Chinese Academy of Sciences, Shanghai200031,China

6Center for Computational Genomics,Beijing Institute of Genomics,Chinese Academy of Sciences,Beijing100101,China

7National Key Laboratory of High Altitude Medicine,High Altitude Medical Research Institute,XiningQinghai810012,China

8School of Life Science and Technology,Shanghai Tech University,Shanghai200031,China

9Collaborative Innovation Center of Genetics and Development,Shanghai200438,China

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 ofGCH1in 50 unrelated Tibetans. Combined with previously published data, we demonstrated manyGCH1variants showing deep divergence between highlander Tibetans and lowlander Han Chinese. Neutrality tests confirmed a signal of positive Darwinian selection onGCH1in Tibetans. Moreover, association analysis indicated that the Tibetan version ofGCH1was significantly associated with multiple physiological traits in Tibetans, including blood nitric oxide concentration, blood oxygen saturation, and hemoglobin concentration. Taken together, we propose thatGCH1plays a role in the genetic adaptation of Tibetans to high altitude hypoxia.

GCH1; Positive selection; Tibetan; Hypoxia adaptation; Nitric oxide; Hemoglobin; Oxygen saturation

INTRODUCTION

Tibetans are a well-known example of successful adaptation to extreme environments at high-altitude. Compared to lowlanders moving to high altitude, Tibetans show greater lung capacity, function, diffusion, and ventilation, as well as lower hemoglobin (Hb) levels, better blood oxygen saturation, low hypoxic pulmonary vasoconstriction, high nitric oxide (NO) concentrations, and lower incidence of reduced birth weight (Beall, 2006; Beall et al., 1997; Erzurum et al., 2007; Wu & Kayser, 2006). These traits have been acquired during a long period of naturalselection at high altitude after the ancestors of modern Tibetans permanently settled on the Qinghai-Tibetan Plateau during the early Upper Paleolithic period (Qi et al., 2013).

Various research groups have compared the genetic differences between Tibetan and Han Chinese using genome-wide scans, with two key genes identified (hypoxia-inducible factor 2α,HIF2α, also calledEPAS1, andEGLN1) in the hypoxic pathway showing deep between-population divergence (Beall et al., 2010; Bigham et al., 2010; Peng et al., 2011; Simonson et al., 2010; Xu et al., 2011; Yi et al., 2010). For these two genes, Tibetan-specific haplotypes have been found, which are highly enriched in Tibetans (～80%), but rare or absent in other world populations (Lorenzo et al., 2014; Peng et al., 2011, 2017; Xiang et al., 2013). The function of selection on these two genes have been shown to cause blunted physiological responses under high altitude hypoxic conditions (Lorenzo et al., 2014; Peng et al., 2017; Xiang et al., 2013). For example, the Tibetan versions ofEPAS1andEGLN1protect Tibetans from the over production of red blood cells (polycythemia), a common deleterious physiological response when lowlanders migrate to high altitude areas (Beall et al., 2010; Erzurum et al., 2007; Lorenzo et al., 2014; Peng et al., 2011, 2017; Petousi et al., 2014; Xiang et al., 2013; Xu et al., 2011; Zhang et al., 2010).

Many genes are likely involved in complex traits like high altitude adaptation. In addition toEPAS1andEGLN1, previous genome-wide studies have identified other candidate genes that might also contribute to genetic adaptation in Tibetans (Beall et al., 2010; Bigham et al., 2010; Peng et al., 2011; Simonson et al., 2010; Xu et al., 2011; Yi et al., 2010). The GTP-cyclohydrolase I (GCH1) gene is a reported candidate that shows relatively deep allelic divergence between Tibetans and Han Chinese, a sign of positive Darwinian selection (Peng et al., 2011).GCH1is located on human chromosome 14 (14q22.1-q22.2), spanning 60.8 kb with six exons. Furthermore,GCH1is a rate-limiting enzyme in thede novosynthesis of tetrahydrobiopterin (BH4). It has been reported that under hypoxia,GCH1can promote cancer growth, and its expression and that of endothelial nitric oxide synthetase (eNOS) is upregulated (Pickert et al., 2013).GCH1is considered as a major factor in maintaining nitric oxide synthetase (NOS) function and normal blood pressure, and its inhibition can increase blood pressure due to NOS uncoupling, which is found in many cardiovascular diseases such as hypertension and atherosclerosis (Antoniades et al., 2006). Hence, the known functional role ofGCH1also makes it a candidate for high-altitude adaptation in Tibetans.

We resequenced an 80.8 kb fragment covering the entire gene region ofGCH1in 50 unrelated Tibetans. Combined with published data, we found signals of positive selection onGCH1in Tibetans, with multiple sequence variants showing deep genetic divergence between highlander Tibetans and lowlander Han Chinese. Genetic association analysis detected significant correlation of theGCH1variants with multiple physiological traits of Tibetans, including blood nitric oxide (NO) concentration, blood oxygen saturation level, and hemoglobin concentration. Hence,GCH1might play an essential role in high-altitude adaptation in Tibetans.

MATERIALS AND METHODS

DNA sample collection and resequencing of GCH1 gene fragment

The 50 unrelated Tibetan samples were obtained from previous research (Peng et al., 2011). We resequenced an 80.8 kb fragment covering the gene region ofGCH1and its flanking sequences (10 kb up- and down-stream ofGCH1). For association analysis, we collected blood samples and extracted DNA from 226 unrelated adult Tibetans, whose physiological data were also collected with written informed consent. The protocol of this study was evaluated and approved by the Internal Review Board of Kunming Institute of China, Chinese Academy of Sciences.

Detection of selection on GCH1 in Tibetans

The initially identified sequence variants were filtered by removing single nucleotide polymorphisms (SNPs) showing a significant deviation from the Hardy-Weinberg Equilibrium (HWE<0.000 1), and those with an excessive missing genotype rate (MGR>0.05). Four methods were used for selection testing, including two allele-frequency-based and two haplotype-based tests. Locus specificFSTwas calculated between 83 Tibetans and three reference populations (103 Han Chinese, CHB; 99 Europeans, CEU; and 108 Africans, YRI) following Weir & Cockerham (1984). The Tajima’sD-test was also performed following the published procedure (Tajima, 1989).

The two haplotype-based tests included the iHS and XP-EHH tests (Sabeti et al., 2002). The iHS score was calculated for each site using selscan (Szpiech & Hernandez, 2014) based on the phased haplotypes, and only those allelic loci whose ancestral alleles were known with certainty were included in the analysis (Voight et al., 2006). The XP-EHH (Sabeti et al., 2007) analysis was used to detect the extended haplotypes due to positive selection. Han Chinese were used as the reference population in the XP-EHH test. We computed XP-EHH scores using selscan (Szpiech & Hernandez, 2014) based on phased haplotypes of Tibetans and Han Chinese. The XP-EHH score of each SNP was standardized by the mean XP-EHH and the standard deviation (SD) over the genome.

Functional prediction of GCH1 candidate SNPs

Functional enrichment analyses of the candidate variants were performed using the Combined Annotation Dependent Depletion (CADD) database (http://krishna.gs.washington.edu/ download/CADD/v1.3/1000G_phase3_inclAnno.tsv.gz), which incorporates data from ENCODE (http://genome.ucsc.edu/ ENCODE/) and NIH Roadmap Epigenomics using ChromHMM (https://sites.google.com/site/anshulkundaje/projects/epigenom e roadmap#TOC-Core-Integrative-chromatin-state-maps-127-Epigenomes-) (Ernst & Kellis, 2012).

Evolutionary constraint is an indication of functional importance. We used Genome Evolutionary Rate Profiling (GERP) to evaluate how conserved a test SNP was compared with SNP-containing sequences from different species (http:// mendel.stanford.edu/SidowLab/downloads/gerp/). The GERP++ method was used to calculate site-specific RS scores (Davydovet al., 2010). A positive GERP++ score indicates evolutionary constraint, and the greater the score, the greater the level of evolutionary constraint inferred to be acting on that site.

The H3K4Me1 value indicates the maximum ENCODE H3K4 methylation level (maximum value observed across 16 ENCODE cell lines at a given position), suggestive of an enhancer or other regulatory activities. The H3K4Me3 value indicates the maximum ENCODE H3K4 trimethylation level, and is an indication of a promoter. The DNasePvalue indicates evidence for open chromatin. The transcription factor binding site (TFBS) is indicated by the number of different overlapping ChIP transcription factor binding sites. In addition, the splice site results indicate if the tested variant is within an ACCEPTOR or a DONOR (Supplementary Table S1).

Measurements of physiological traits

Physiological data and blood samples were collected from 226 unrelated Tibetans permanently residing in Bange county (n=135, 37.41±3.8 years old) at an elevation of 4 700 m and Lhasa city (n=91, 35.33±6.8 years old) at an elevation of 3 600 m. Written informed consent from all participants was obtained. For physiological parameters, we collected hemoglobin (Hb) concentration, arterial oxygen saturation (SaO2) level, and blood nitric oxide concentration, representing key adaptive physiological traits in Tibetans (Wu & Kayser, 2006).

The Hb level in blood was measured using a HemoCue Hb 201+ analyzer (Angelholm, Sweden). The SaO2level was measured using fingertip blood with a hand-held pulse oximeter (Nellcor NPB-40, CA, USA). Blood NO was measured using a nitric oxide analyzer (Sievers Model-280, GE Analytical Instruments; Boulder, CO, USA).

SNP genotyping and association analysis

We selected nine tag SNPs, covering the entire gene region ofGCH1(Table 1). Genotyping was conducted by SNaPshot on an ABI 3130 sequencer (Applied Bio-Systems, Forster City, CA, USA). Genetic association analysis was executed using PLINK 1.07 (Purcell et al., 2007). We used an additive model in the association analysis as all candidate SNPs were non-coding and likely influenced the level of gene expression. Permutations (100 000 times for each test) were performed for statistical assessment and correction for multiple tests.

Table 1 Association of nine GCH1 variants with three physiological traits in Tibetans

RESULTS

Resequencing of GCH1 in Tibetans and tests of selectionPrevious DNA array-based genome-wide studies have only covered a limited number ofGCH1sequence variants (Peng et al., 2011). To obtain complete sequence data ofGCH1, we first resequenced an 80.8 kb fragment covering the entire gene region ofGCH1(60.8 kb) as well as the flanking sequences (10 kb from upstream and 10 kb from downstream regions). In total, we sequenced 50 unrelated Tibetan individuals, as reported previously (Peng et al., 2011). In addition, we obtained theGCH1sequences of 33 Tibetans from recently published whole genome sequencing data (Lu et al., 2016), for a final sample size of 83.

We identified a total of 384GCH1sequence variants (SNPs) in the 83 unrelated Tibetans, among which 245 were shared between Tibetans and the three lowland reference populations from the 1000 Genomes Project (http://www.1000genomes.org) (103 Han Chinese, CHB; 99 Europeans, CEU and 108 Africans, YRI). The remaining 139 SNPs were rare in Tibetans (<1.0%), and therefore not included in our analysis. The 245 SNPs were all located in the non-coding regions ofGCH1.

To detect whetherGCH1was under selection in Tibetans, we conducted four different tests of selection, including two allelefrequency-based tests (FSTand Tajima’sD) and two haplotypebased tests (iHS and XP-EHH). We identified 49GCH1SNPs with large between-population (Tibetan vs. Han) divergence (FST>0.2), much larger than the genome average (FST=0.03). These high-FSTvariants also showed high iHS and XP-EHH values (iHS>0.2 and/or XP-EHH>0.2) (Figure 1). They were aggregated in a relatively short region (7.5 kb) covering intron-1 and intron-2 ofGCH1. These results suggest a clear signal of positive Darwinian selection onGCH1in Tibetans (Supplementary Table S1)

Functional prediction and genetic association analysis

To determine the potential function of the SNPs with selection signals, we performed functional predictions using evolutionary constraint (GERP), transcription factor binding sites (TFBS), splicing motif, H3K4Me1/H3K4Me3 sites, and DNase-I hypersensitive sites. Results showed that there might be multiple functional sites, reflected by the consistent signals of different functional predictions (Supplementary Table S1).

To test whether these candidateGCH1variants contribute to the adaptive physiological traits of Tibetans, we measured blood NO concentrations, SaO2levels, and Hb concentrations in 226 unrelated adult Tibetans (refer to methods for details). For genotyping, we selected nine tag SNPs to represent the entire gene region ofGCH1. Notably, four SNPs (rs7148266, rs17128004, rs4411417, and rs10136972) showed significant association with blood NO levels, with the adaptive alleles exhibiting decreased NO. Each SNP accounted for 1.88%–3.76% of NO variance (Figure 2). Similar results were observed when males and females were analyzed separately (Table 1). Additionally, two SNPs (rs146540091 and rs137863726) showed association with hemoglobin concentration, with the adaptive alleles having lower Hb levels, consistent with previous results forEPAS1andEGLN1(Peng et al., 2011, 2017; Xiang et al., 2013). Another SNP (rs137863726) was associated with blood oxygen saturation, with the adaptive allele exhibiting higher SaO2levels.

For the four SNPs (rs7148266, rs17128004, rs4411417, and rs10136972) showing significant association with blood NO levels, functional prediction analysis indicated that they were located in theGCH1intron regions with peak signals for H3K4Me1, H3K4Me3, DNase-I, and TFBS. For example, the H3K4Me1 peak values for rs10136972 and rs4411417 were 8.2 and 5.9, respectively, indicating their potential role in gene expression regulation ofGCH1.

DISCUSSION

Hypoxia serves as a key stress in high altitude environments. For lowlanders, prolonged exposure to high altitude hypoxia can cause chronic mountain sickness, reflected by the over production of red blood cells (polycythemia) as well as other deleterious physiological changes (Hackett & Roach, 2001; Macinnis et al., 2010). Tibetans are genetically adapted to high altitude hypoxia, and exhibit blunted physiological responses, e.g., relatively low hemoglobin levels (Lorenzo et al., 2014; Peng et al., 2011, 2017; Xiang et al., 2013). Key hypoxic pathway genesEPAS1andEGLN1were reported to be responsible for these blunted physiological responses (Lorenzo et al., 2014; Peng et al., 2017; Xiang et al., 2013). However, while the Tibetan version of these two genes provide protection against polycythemia, they do not explain all physiological changes in Tibetans, suggesting there might be other genes involved given the complexity of high altitude adaptation.

Based on resequencing and population analysis, we confirmed a signal of selection onGCH1in Tibetans. We identified more than 40 variants showing deep allelic divergence between highlander Tibetans and lowlander Han Chinese, with some having potential functional effects based on prediction using ENCODE.GCH1is a rate-limiting enzyme, acting as a crucial factor for maintaining normal NO synthetase function and blood pressure. Inhibition ofGCH1activity is related to several cardiovascular diseases, withGCH1found to prevent hypoxia-induced pulmonary hypertension (Khoo et al., 2005). Hence, the function of selection onGCH1in Tibetans is expected to help maintain proper cardiovascular function at high altitude.

Figure 2 Genetic association of seven GCH1 variants with three physiological traits (blood NO concentration, blood oxygen saturation level, and hemoglobin concentration)

Hypoxic pulmonary vasoconstriction and pulmonary vascular structural remodeling are dominant pathophysiological characteristics of hypoxic pulmonary hypertension (Galiè et al., 2016; McLaughlin & McGoon, 2006; Schermuly et al., 2011). When lowlanders move to high altitudes, pulmonary hypertension usually occurs within a few weeks (Wilkins et al., 2015; Wu & Kayser, 2006); however, Tibetans rarely develop this condition. We showed thatGCH1SNPs were associated with NO levels in the blood.GCH1is involved in the synthesis of tetrahydrobiopterin (BH4), a vital regulator of eNOS, the endothelial-form enzyme producing NO, an important molecule for vasodilation, which is considered the main reason for the superior blood flow and pulmonary pressure in Tibetans (Pickertet al., 2013). In hph-1 mice, deficiency of BH4 causes hypoxiainduced pulmonary hypertension even under normoxic conditions (Khoo et al., 2005). The overexpression ofGCH1in mice could prevent hypoxia-induced pulmonary hypertension due to the augmentation of BH4 (Khoo et al., 2005). Hence, it is possible thatGCH1regulates pulmonary vasoconstriction responses in Tibetans by influencing NO production in the blood. We observed fourGCH1variants showing significant association with blood NO levels. As these variants are located in theGCH1intron regions with peak enhancer and/or promoter activity signals, they are likely involved in the regulation ofGCH1expression and eventually affect blood NO production, which needs further investigation. We also observed associations ofGCH1SNPs with oxygen saturation and hemoglobin; however, the underlying molecular mechanisms are yet to be studied.

In summary, we demonstrated thatGCH1has been under positive selection in Tibetans. We identified many variants with deep allelic divergence between Tibetans and lowlanders. The association and known function results suggest the potential involvement ofGCH1in the regulation of multiple physiological traits in Tibetans.

ACKNOWLEDGEMENTS

We are grateful to all the volunteers participated in this study.

REFERENCES

Antoniades C, Shirodaria C, Warrick N, Cai SJ, de Bono J, Lee J, Leeson P, Neubauer S, Ratnatunga C, Pillai R, Refsum H, Channon KM. 2006. 5-methyltetrahydrofolate rapidly improves endothelial function and decreases superoxide production in human vessels: effects on vascular tetrahydrobiopterin availability and endothelial nitric oxide synthase coupling.Circulation,114(11): 1193-1201.

Beall CM, Strohl KP, Blangero J, Williams-Blangero S, Almasy LA, Decker MJ, Worthman CM, Goldstein MC, Vargas E, Villena M, Soria R, Alarcon AM, Gonzales C. 1997. Ventilation and hypoxic ventilatory response of Tibetan and Aymara high altitude natives.American Journal of Physical Anthropology,104(4): 427-447.

Beall CM. 2006. Andean, Tibetan, and Ethiopian patterns of adaptation to high-altitude hypoxia.Integrative and Comparative Biology,46(1): 18-24.

Beall CM, Cavalleri GL, Deng LB, Elston RC, Gao Y, Knight J, Li CH, Li JC, Liang Y, McCormack M, Montgomery HE, Pan H, Robbins PA, Shianna KV, Tam SC, Tsering N, Veeramah KR, Wang W, Wangdui P, Weale ME, Xu YM, Xu Z, Yang L, Zaman MJ, Zeng CQ, Zhang L, Zhang XL, Zhaxi P, Zheng YT. 2010. Natural selection onEPAS1(HIF2α) associated with low hemoglobin concentration in Tibetan highlanders.Proceedings of the National Academy of Sciences of the United States of America,107(25): 11459-11464.

Bigham A, Bauchet M, Pinto D, Mao XY, Akey JM, Mei R, Scherer SW, Julian CG, Wilson MJ, Herráez DL, Brutsaert T, Parra EJ, Moore LG, Shriver MD. 2010. Identifying signatures of natural selection in Tibetan and Andean populations using dense genome scan data.PLoS Genetics,6(9): e1001116.

Davydov EV, Goode DL, Sirota M, Cooper GM, Sidow A, Batzoglou S. 2010. Identifying a high fraction of the human genome to be under selective constraint using GERP++.PLoS Computational Biology,6(12): e1001025.

Ernst J, Kellis M. 2012. ChromHMM: automating chromatin-state discovery and characterization.Nature Methods,9(3): 215-216.

Erzurum SC, Ghosh S, Janocha AJ, Xu W, Bauer S, Bryan NS, Tejero J, Hemann C, Hille R, Stuehr DJ, Feelisch M, Beall CM. 2007. Higher blood flow and circulating NO products offset high-altitude hypoxia among Tibetans.Proceedings of the National Academy of Sciences of the United States of America,104(45): 17593-17598.

Galiè N, Humbert M, Vachiery JL, Gibbs S, Lang I, Torbicki A, Simonneau G, Peacock A, Noordegraaf AV, Beghetti M, Ghofrani A, Sanchez MAG, Hansmann G, Klepetko W, Lancellotti P, Matucci M, McDonagh T, Pierard LA, Trindade PT, Zompatori M, Hoeper M, Aboyans V, Vaz Carneiro A, Achenbach S, Agewall S, Allanore Y, Asteggiano R, Badano LP, Barberà JA, Bouvaist H, Bueno H, Byrne RA, Carerj S, Castro G, Erol C, Falk V, Funck-Brentano C, Gorenflo M, Granton J, Iung B, Kiely DG, Kirchhof P, Kjellstrom B, Landmesser U, Lekakis J, Lionis C, Lip GY, Orfanos SE, Park MH, Piepoli MF, Ponikowski P, Revel MP, Rigau D, Rosenkranz S, V?ller H, Zamorano JL. 2016. 2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: the joint task force for the diagnosis and treatment of pulmonary hypertension of the european society of cardiology (ESC) and the european respiratory society (ERS): endorsed by: association for european paediatric and congenital cardiology (AEPC), international society for heart and lung transplantation (ISHLT).European Heart Journal,37(1): 67-119.

Hackett PH, Roach RC. 2001. High-altitude illness.The New England Journal of Medicine,345(2): 107-114.

Khoo JP, Zhao L, Alp NJ, Bendall JK, Nicoli T, Rockett K, Wilkins MR, Channon KM. 2005. Pivotal role for endothelial tetrahydrobiopterin in pulmonary hypertension.Circulation,111(16): 2126-2133.

Lorenzo FR, Huff C, Myllym?ki M, Olenchock B, Swierczek S, Tashi T, Gordeuk V, Wuren T, Ri-Li G, McClain DA, Khan TM, Koul PA, Guchhait P, Salama ME, Xing JC, Semenza GL, Liberzon E, Wilson A, Simonson TS, Jorde LB, Kaelin WG, Jr., Koivunen P, Prchal JT. 2014. A genetic mechanism for Tibetan high-altitude adaptation.Nature Genetics,46(9): 951-956.

Lu DS, Lou HY, Yuan K, Wang XJ, Wang YC, Zhang C, Lu Y, Yang X, Deng L, Zhou Y, Feng QD, Hu Y, Ding QL, Yang YJ, Li SL, Jin L, Guan YQ, Su B, Kang LL, Xu SH. 2016. Ancestral origins and genetic history of tibetan highlanders.The American Journal of Human Genetics,99(3): 580-594.

Macinnis MJ, Koehle MS, Rupert JL. 2010. Evidence for a genetic basis for altitude illness: 2010 update.High Altitude Medicine & Biology,11(4): 349-368. McLaughlin VV, McGoon MD. 2006. Pulmonary arterial hypertension.Circulation,114(13): 1417-1431.

Peng Y, Yang ZH, Zhang H, Cui CY, Qi XB, Luo XJ, Tao X, Wu TY, Ouzhuluobu, Basang, Ciwangsangbu, Danzengduojie, Chen H, Shi H, Su B. 2011. Genetic variations in tibetan populations and high-altitude adaptation at the himalayas.Molecular Biology and Evolution,28(2): 1075-1081.

Peng Y, Cui CY, He YX, Ouzhuluobu, Zhang H, Yang DY, Zhang Q, Bianbazhuoma, Yang LX, He YB, Xiang K, Zhang XM, Bhandari S, Shi P, Yangla, Dejiquzong, Baimakangzhuo, Duojizhuoma, Pan YY, Cirenyangji, Baimayangji, Gonggalanzi, Bai CJ, Bianba, Basang, Ciwangsangbu, Xu SH, Chen H, Liu SM, Wu TY, Qi XB, Su B. 2017. Down-regulation ofEPAS1transcription and genetic adaptation of tibetans to high-altitude hypoxia.Molecular Biology and Evolution,34(4): 818-830.

Petousi N, Croft QPP, Cavalleri GL, Cheng HY, Formenti F, Ishida K, Lunn D, McCormack M, Shianna KV, Talbot NP, Ratcliffe PJ, Robbins PA. 2014.Tibetans living at sea level have a hyporesponsive hypoxia-inducible factor system and blunted physiological responses to hypoxia.Journal of Applied Physiology,116(7): 893-904.

Pickert G, Lim HY, Weigert A, H?ussler A, Myrczek T, Waldner M, Labocha S, Ferreirós N, Geisslinger G, L?tsch J, Becker C, Brüne B, Tegeder I. 2013. Inhibition of GTP cyclohydrolase attenuates tumor growth by reducing angiogenesis and M2-like polarization of tumor associated macrophages.International Journal of Cancer,132(3): 591-604.

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC. 2007. PLINK: a tool set for whole-genome association and population-based linkage analyses.The American Journal of Human Genetics,81(3): 559-575.

Qi XB, Cui CY, Peng Y, Zhang XM, Yang ZH, Zhong H, Zhang H, Xiang K, Cao XY, Wang Y, Ouzhuluobu, Basang, Ciwangsangbu, Bianba, Gonggalanzi, Wu TY, Chen H, Shi H, Su B. 2013. Genetic evidence of paleolithic colonization and neolithic expansion of modern humans on the tibetan plateau.Molecular Biology and Evolution,30(8): 1761-1778.

Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, Ackerman HC, Campbell SJ, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander ES. 2002. Detecting recent positive selection in the human genome from haplotype structure.Nature,419(6909): 832-837.

Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie XH, Byrne EH, McCarroll SA, Gaudet R, Schaffner SF, Lander ES, The International HapMap Consortium. 2007. Genome-wide detection and characterization of positive selection in human populations.Nature,449(7164): 913-918.

Schermuly RT, Ghofrani HA, Wilkins MR, Grimminger F. 2011. Mechanisms of disease: pulmonary arterial hypertension.Nature Reviews Cardiology,8(8): 443-455.

Simonson TS, Yang YZ, Huff CD, Yun HX, Qin G, Witherspoon DJ, Bai ZZ, Lorenzo FR, Xing JC, Jorde LB, Prchal JT, Ge RL. 2010. Genetic evidence for high-altitude adaptation in Tibet.Science,329(5987): 72-75.

Szpiech ZA, Hernandez RD. 2014. Selscan: an efficient multithreaded program to perform EHH-based scans for positive selection.Molecular Biology and Evolution,31(10): 2824-2827.

Tajima F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.Genetics,123(3): 585-595.

Voight BF, Kudaravalli S, Wen XQ, Pritchard JK. 2006. A map of recent positive selection in the human genome.PLoS Biology,4(3): e72.

Weir BS, Cockerham CC. 1984. Estimating F-statistics for the analysis of population structure.Evolution,38(6): 1358-1370.

Wilkins MR, Ghofrani HA, Weissmann N, Aldashev A, Zhao L. 2015. Pathophysiology and treatment of high-altitude pulmonary vascular disease.Circulation,131(6): 582-590.

Wu TY, Kayser B. 2006. High altitude adaptation in Tibetans.High Altitude Medicine & Biology,7(3): 193-208.

Xiang K, Ouzhuluobu, Peng Y, Yang ZH, Zhang XM, Cui CY, Zhang H, Li M, Zhang YF, Bianba, Gonggalanzi, Basang, Ciwangsangbu, Wu TY, Chen H, Shi H, Qi XB, Su B. 2013. Identification of a Tibetan-specific mutation in the hypoxic geneEGLN1and its contribution to high-altitude adaptation.Molecular Biology and Evolution,30(8): 1889-1898.

Xu SH, Li SL, Yang YJ, Tan JZ, Lou HY, Jin WF, Yang L, Pan XD, Wang JC, Shen YP, Wu BL, Wang HY, Jin L. 2011. A genome-wide search for signals of high-altitude adaptation in Tibetans.Molecular Biology and Evolution,28(2): 1003-1011.

Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, Xu X, Jiang H, Vinckenbosch N, Korneliussen TS, Zheng HC, Liu T, He WM, Li K, Luo RB, Nie XF, Wu HL, Zhao MR, Cao HZ, Zou J, Shan Y, Li SZ, Yang Q, Asan, Ni PX, Tian G, Xu JM, Liu X, Jiang T, Wu RH, Zhou GY, Tang MF, Qin JJ, Wang T, Feng SJ, Li GH, Huasang, Luosang JB, Wang W, Chen F, Wang YD, Zheng XG, Li Z, Bianba Z, Yang G, Wang XP, Tang SH, Gao GY, Chen Y, Luo Z, Gusang L, Cao Z, Zhang QH, Ouyang WH, Ren XL, Liang HQ, Zheng HS, Huang YB, Li JX, Bolund L, Kristiansen K, Li YR, Zhang Y, Zhang XQ, Li RQ, Li SG, Yang HM, Nielsen R, Wang J, Wang J. 2010. Sequencing of 50 human exomes reveals adaptation to high altitude.Science,329(5987): 75-78.

Zhang L, Chen WZ, Liu YJ, Hu X, Zhou K, Chen L, Peng S, Zhu H, Zou HL, Bai J, Wang ZB. 2010. Feasibility of magnetic resonance imaging-guided high intensity focused ultrasound therapy for ablating uterine fibroids in patients with bowel lies anterior to uterus.European Journal of Radiology,73(2): 396-403.

24 March 2017; Accepted: 27 April 2017

s: This study was supported by grants from the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB13010000), the National Natural Science Foundation of China (91631306 to BS, 31671329 to XQ, 31460287 to Ou., 31501013 to HZ and 31360032 to CC), the National 973 program (2012CB518202 to TW), the State Key Laboratory of Genetic Resources and Evolution (GREKF15-05, GREKF16-04), and the Zhufeng Scholar Program of Tibetan University

#Authors contributed equally to this work

*Corresponding authors, E-mail: sub@mail.kiz.ac.cn; qixuebin@mail.kiz.ac.cn DOI: 10.24272/j.issn.2095-8137.2017.037