- Open Access
EPAS1 and VEGFA gene variants are related to the symptoms of acute mountain sickness in Chinese Han population: a cross-sectional study
Military Medical Research volume 7, Article number: 35 (2020)
More people ascend to high altitude (HA) for various activities, and some individuals are susceptible to HA illness after rapidly ascending from plains. Acute mountain sickness (AMS) is a general complaint that affects activities of daily living at HA. Although genomic association analyses suggest that single nucleotide polymorphisms (SNPs) are involved in the genesis of AMS, no major gene variants associated with AMS-related symptoms have been identified.
In this cross-sectional study, 604 young, healthy Chinese Han men were recruited in June and July of 2012 in Chengdu, and rapidly taken to above 3700 m by plane. Basic demographic parameters were collected at sea level, and heart rate, pulse oxygen saturation (SpO2), systolic and diastolic blood pressure and AMS-related symptoms were determined within 18–24 h after arriving in Lhasa. AMS patients were identified according to the latest Lake Louise scoring system (LLSS). Potential associations between variant genotypes and AMS/AMS-related symptoms were identified by logistic regression after adjusting for potential confounders (age, body mass index and smoking status).
In total, 320 subjects (53.0%) were diagnosed with AMS, with no cases of high-altitude pulmonary edema or high-altitude cerebral edema. SpO2 was significantly lower in the AMS group than that in the non-AMS group (P = 0.003). Four SNPs in hypoxia-inducible factor-related genes were found to be associated with AMS before multiple hypothesis testing correction. The rs6756667 (EPAS1) was associated with mild gastrointestinal symptoms (P = 0.013), while rs3025039 (VEGFA) was related to mild headache (P = 0.0007). The combination of rs6756667 GG and rs3025039 CT/TT further increased the risk of developing AMS (OR = 2.70, P < 0.001).
Under the latest LLSS, we find that EPAS1 and VEGFA gene variants are related to AMS susceptibility through different AMS-related symptoms in the Chinese Han population; this tool might be useful for screening susceptible populations and predicting clinical symptoms leading to AMS before an individual reaches HA.
Chinese Clinical Trial Registration, ChiCTR-RCS-12002232. Registered 31 May 2012.
Every year, millions of people ascend to high altitude (HA) above 2500 m for travel, work or military activities. However, some individuals are susceptible to HA illness or diseases, including acute mountain sickness (AMS), high-altitude pulmonary edema (HAPE) and high-altitude cerebral edema (HACE) . AMS, a hypoxic syndrome, has characteristic symptoms, including headache, anorexia, dizziness, nausea, vomiting, fatigue, and sleep disturbance . The incidence of AMS ranges from 10 to 85%, depending on the altitude reached, the ascension rate and individual susceptibility . Most symptoms are self-limiting, and affected individuals recover after days of rest and acclimatization. Under severe conditions, some patients may develop life-threatening HAPE or HACE . AMS has become a general complaint that limits activities of daily living and work capabilities at HA. Therefore, it is important to understand the parameters of susceptibility to AMS.
To date, the pathophysiology of AMS remains largely unknown. Although the primary etiology is believed to be hypobaric hypoxia, populations in highland regions, such as the Ethiopian Highlands, the Andes and the Tibetan Plateau, have a much lower risk of suffering from any mountain illness . Interestingly, even under the same acute hypoxic conditions, certain lowland individuals are more susceptible to AMS than others. It is generally assumed that genetic background plays a pivotal role in the differences in susceptibility between lowlanders and highlanders . Many genome-wide studies have revealed genetic variations in native highlanders related to adaptation at HA [5, 7,8,9,10]. Some important genes in the hypoxia-inducible factor (HIF) pathway are closely related to hypoxia adaptation in Tibetan and Andean populations . Thus, genes involved in the HIF signaling pathway may be the basis for adaptive changes in plateau populations.
HIFs are transcription factors that directly or indirectly regulate hundreds of genes involved in angiogenesis, cell growth, apoptosis, energy metabolism and vasomotor regulation [11, 12]. HIF1α and HIF2α (also known as endothelial PAS domain-containing protein 1, EPAS1) are basic helix-loop-helix PAS domain transcription factors that are expressed in all metazoan organisms and are composed of α and β subunits . Under hypoxic conditions, HIF1α regulates the transcription of hundreds of genes in a cell type-specific manner . Hypoxia-sensitive gene pathways play key regulatory roles in red blood cell production (erythropoiesis), development, energy metabolism, angiogenesis, iron metabolism, heart and lung regulation and tumorigenesis in mammals . HIF-1α accumulation and the subsequent upregulation of vascular endothelial growth factor (VEGF) potentially contribute to basement membrane damage, and exaggerated edema might be involved in the pathophysiology of AMS . VEGF is a major angiogenic factor that plays an important role in regulating endothelial cell proliferation and vascular permeability, which has been shown to be markedly upregulated under hypoxic conditions. Two oxygen sensors, prolyl hydroxylase domain protein 2 (PHD2), also known as Egl nine homolog 1 (EGLN1), and factor inhibiting HIF-1α (FIH-1), also known as hypoxia-inducible factor 1-α inhibitor (HIF-1AN), play pivotal roles in regulating the HIF-1α pathway.
However, the role of polymorphisms in HIF pathway genes in the response to acute hypoxia, which may lead to AMS, has not been clarified. In a preliminary genome-wide association study of the Nepalese population, four single nucleotide polymorphisms (SNPs) in the FAM149A gene were found to be associated with AMS, but these results were not verified in two validation cohorts . Moreover, few studies have examined the genetic component of susceptibility to AMS in lowlanders, such as the Chinese Han population, following acute exposure to HA. Our group showed that SNPs in EPAS1 and the 5′-untranslated region (UTR) of EGLN1 may be associated with a high risk of AMS [17, 18]. Ding et al.  reported a possible relationship between VEGFA polymorphisms and AMS in a relatively small Chinese Han population. In another study, VEGFA SNPs in transcription factor binding sites were analyzed and compared between relatively small Han Chinese and Tibetan populations (< 100), but no sites were associated with HA sickness, although some sites were associated with arterial oxygen saturation (SaO2) . Since most of the Chinese Han population resides in the lowland region, they may be an appropriate population for investigating the association between genetic factors and susceptibility or resistance to AMS. Elucidating these associations may help millions of lowlanders predict their response to travel to plateaus and provide novel genetic targets for potential clinical intervention. Furthermore, no previous studies have reported genetic associations with AMS-related symptoms.
The AMS diagnosis system was recently revised by the Lake Louise AMS Score Consensus Committee in 2018 . Based on two independent studies [22, 23] and suggestions from the consensus committee, disrupted sleep was excluded as a necessary symptom for the diagnosis of AMS . It is unknown whether this exclusion will change our previous genetic findings regarding AMS. Thus, based on our field study of a relatively large population in 2012, we reevaluated 30 SNPs in 7 genes (EGLN1, EGLN3, EPAS1, HIF1A, HIF1AN, PPARA, and VEGFA) closely related to the HIF pathway according to the new 2018 scoring instructions and analyzed symptoms associated with the onset of AMS. We aimed to clarify the relationship between the genetic background of the Han population and the pathogenesis of AMS as well as related symptoms.
A total of 604 unrelated Chinese Han males aged 18–45 years were recruited in June and July of 2012 in Chengdu (500 m above sea level). All subjects were military officers and soldiers who travelled from Chengdu to Lhasa (3700 m above sea level) by plane in 2 h. Health examinations were performed at Xinqiao Hospital (Chongqing, China) before the trip. Participants were enrolled only if they met the inclusion criteria of residing on the plain for generations without plateau exposure in the past 6 months. Subjects who suffered from diseases with similar clinical features to AMS, cardiovascular or respiratory system diseases, neurological diseases, cerebral vascular diseases, gastrointestinal system diseases or cancer, were excluded. Ethical approval was obtained from the Ethics Committee of Xinqiao Hospital, Army Medical University (number: 2012014; approved May 9, 2012). Each participant provided written informed consent for the collection of basic information, physical parameters and blood samples for further analyses. This study was registrated at Chinese Clinical Trial Registration (ChiCTR-RCS-12002232; Registered 31 May 2012).
Data collection and AMS evaluation
Structured case report questionnaires were designed for participants to record their demographic information, including age, height, weight, body mass index (BMI), and smoking status, which was collected 7 days before travel. Physiological parameters, including heart rate (HR), pulse oxygen saturation (SpO2) and blood pressure (BP), and AMS-related symptoms, including headache, dizziness/light-headedness, gastrointestinal (GI) symptoms and fatigue, were collected 18–24 h after arriving in Lhasa. Participants rested for at least 10 min before physiological parameters were collected. A mercury sphygmomanometer was used to record HR and BP in the upper arm at heart level, and a Finger-Pulse Oximeter 503 (Criticare Systems, Inc., Waukesha, WI, USA) was used to record SpO2. AMS-related symptoms were recorded using the latest consensus of the 2018 Lake Louise scoring system (LLSS) , and subjects with a LLSS score ≥ 3 points and a headache were assigned to the AMS group, while those with a score ≥ 3 but no headache or a score < 3 were assigned to the non-AMS group.
Selection of SNPs
Previous genome-wide association studies suggested that mutations from multiple genetic regions were associated with HA adaption between highlanders and lowlanders. We performed a search on the National Center for Biotechnology Information (NCBI) SNP database (http://www.ncbi.nlm.nih.gov/projects/SNP/) for basic information and the minor allele frequency (MAF) of selected SNPs in Chinese general population. Linkage disequilibrium (LD) analyses were performed and linkage maps were drawn using Haploview. As shown in Additional file 1, the correlation intensity between SNPs was evaluated by the values of LD coefficient (D’), logarithm of odds (LOD) and r2. Based on LD analyses, tag SNPs which represented neighboring SNPs in the genomic region were retained. In EGLN1, we selected rs2066140, rs508618, rs1538667, rs1339891 and rs2275279 which located in introns, rs2153364, rs1361384 and rs1339894 at 5′-UTR, rs12757362 and at the upstream of the gene [8, 10, 12, 18]. In EPAS1, we selected rs13419896, rs4953354 and rs6756667 in introns and rs1868092 at the downstream of the gene [7, 12, 19]. In VEGFA, we selected rs3025039 and rs10434 at 3′-UTR and rs1413711 in introns . In PPARA, we selected rs4253623, rs135538, rs4253681 and rs4253747 in introns, rs7292407 which located 90 kbp from the transcriptional initiation site of gene [7, 8]. In EGLN3, we selected rs1680710 and rs11156819 at 3′-UTR . In HIF1A, we selected rs2301104, rs12434438, rs2301112, rs2301113 and rs11549467 in introns, while in HIF1AN, we selected rs2295778 in introns, rs10883512 at 3′-UTR [9, 18]. In total, 30 tag SNPs on 7 HIF-related genes were included in this study.
DNA extraction and genotyping
For DNA extraction, we collected approximately 5 ml of venous blood from each participant into EDTA tubes for anti-coagulation purposes. Genomic DNA was extracted from whole blood using the PAXgene Blood DNA kit (Qiagen, Hilden, Germany) and was stored at − 20 °C. Sequenom® Assay Design software (version 3.1, Sequenom Inc., San Diego, CA, USA) was used to design the polymerase chain reaction (PCR) primers listed in Additional file 2. After PCR amplification and purification, the DNA sequence containing the target SNP was subjected to single-base extension in a system with ddNTPs using iPLEX Gold technology. Then, the products were analyzed by matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) (Sequenom Inc., San Diego, CA, USA), and polymorphisms were genotyped on the basis of different detection peaks. All genotyping was performed in a blinded fashion, and 10% of the DNA samples were subjected to repeat detection.
The statistical analysis was performed with SPSS software (version 24.0; SPSS, Inc., Chicago, IL, USA) and SNPStats online software (http://bioinfo.iconcologia.net/SNPstats). Continuous data were tested for normality using the Kolmogorov-Smirnov test. Variables with a normal distribution are presented as the mean ± standard deviation, and comparisons between groups were evaluated with the two-sided Student’s t-test. Variables that did not show a normal distribution are presented as medians and quartiles, and comparisons between groups were evaluated with the Mann-Whitney U test. Categorical variables are expressed as cases and percentages. The chi-square test was used to compare categorical variables, Hardy-Weinberg equilibrium (HWE) and allelic frequencies of SNPs between groups. Binary logistic regression was used to analyze associations between genotypes and AMS/AMS-related symptoms, which were further adjusted for potential confounders, including age, BMI and smoking status. Linkage disequilibrium analyses were performed with Haploview 4.2 software. False discovery rate (FDR) created by Benjamini and Hochberg was used for multiple hypothesis testing correction . A two-sided P value < 0.05 indicated a statistically significant difference.
Table 1 shows the basic characteristics of the subjects; 53.0% [n = 320; 23 (20, 25) years old] were diagnosed with AMS, and 47.0% [n = 284; 22 (20, 25) years old] were non-AMS. All participants were healthy Chinese Han soldiers who had lived in the plains for generations. No significant differences between groups were found in age, height, weight, body mass index (BMI), smoking status, HR, systolic BP (SBP) or diastolic BP (DBP). AMS patients had significantly lower SpO2 than non-AMS subjects [89% (86, 90) vs. 89% (87, 91), P = 0.003]. Although the mean values were close, the distribution frequency of SpO2 < 90% was higher in AMS than that in non-AMS group (66.2% vs. 56.0%, Fig. 1).
Associations between SNPs and AMS
Basic information of 30 SNPs in the EGLN1, EPAS1, VEGFA, PPARA, HIF1A, HIF1AN and EGLN3 genes among subjects is presented in Additional file 3. Of these SNPs, rs1538667, rs1361384, rs1339894, rs12757362, rs1680710 and rs11549467 exhibited MAF less than 0.05, and the distribution of the rs1413711 G/A allele was not in HWE. Therefore, the above 7 SNPs were excluded, and further association analyses were conducted between the other 23 SNPs and AMS as well as AMS-related symptoms.
The exact mode of genetic inheritance was not clear, so we applied codominant, dominant and recessive models for genetic analyses. Ultimately, 4 SNPs (EGLN1-rs2153364, EPAS1-rs6756667, VEGFA-rs3025039 and PPARA-rs7292407) were found to be associated with the risk of developing AMS, even after adjusted by age, BMI and smoking status, as well as the results after FDR correction. The negative results of the other 19 SNPs were shown in Additional file 4.
The distributions of the EGLN1-rs2153364 A allele (AMS, 49.3%; non-AMS, 55.4%) and G allele (AMS, 50.7%; non-AMS, 44.6%) were significantly different between the AMS group and the non-AMS group (P = 0.044). Under the dominant model, the AG/GG genotype was associated with significantly increased AMS risk (OR = 1.56; 95% CI 1.07–2.28; P = 0.021) after adjusting for age, BMI and smoking status. However, the result was not significant after FDR correction (Q = 0.242). No significant results were found under the codominant or recessive model (Table 2).
The distributions of the EPAS1-rs6756667 G allele (AMS, 91.1%; non-AMS, 85.0%) and A allele (AMS, 8.9%; non-AMS, 15.0%) were significantly different between groups (P = 0.001). Compared with individuals carrying the G allele, those carrying the A allele had a significantly decreased risk of AMS. The genotype frequencies of GG, GA and AA were 82.2% (n = 263), 17.8% (n = 57) and 0.0% (n = 0), respectively, in the AMS group and 72.2% (n = 205), 25.7% (n = 73) and 2.1% (n = 6), respectively, in the non-AMS group (P = 0.0005) and retained significance after FDR (Q = 0.012). Under the dominant model, the GA/AA genotype was associated with decreased AMS risk (OR = 0.56; 95% CI 0.38–0.82; P = 0.003) after adjusted for age, BMI and smoking status. FDR testing turned out to be insignificant (Q = 0.069) (Table 3).
The frequencies of the VEGFA-rs3025039 C allele (AMS, 82.2%; non-AMS, 87.5%) and T allele (AMS, 17.8%; non-AMS, 12.5%) were significantly different (P = 0.011). The frequencies of the CC, CT and TT alleles were 67.3% (n = 214), 29.9% (n = 95) and 2.8% (n = 9), respectively, in the AMS group and 75.7% (n = 215), 23.6% (n = 67) and 0.7% (n = 2), respectively, in the non-AMS group. Under the dominant model, the CT/TT genotype was associated with increased AMS risk (OR = 1.51; 95% CI 1.05–2.16; P = 0.024) after adjusting for age, BMI and smoking status, but insignificant after FDR (Q = 0.184). No significant results were found under the recessive model (Table 4).
As for PPARA-rs7292407, the allelic frequency was significantly different between the two groups (P = 0.014). Compared with individuals carrying the C allele, those carrying the A allele showed a significantly decreased risk of developing AMS. Under the dominant model, the AC/AA genotype was associated with decreased AMS risk (OR = 0.68; 95% CI 0.47–0.99; P = 0.046) after adjusting, but insignificant after FDR (Q = 0.265). And no significant results were obtained with the codominant or recessive model (Table 5).
Associations between SNPs and AMS-related symptoms
Binary logistic regression analyses were conducted, and the dominant model was chosen to show correlations between SNPs and AMS-related symptoms (Table 6). Under the dominant model, subjects with the EPAS1-rs6756667 GA/AA genotype were less likely to suffer from AMS-related GI symptoms (OR = 0.53; 95% CI 0.32–0.89; P = 0.013), and the result was marginally significant after FDR correction (Q = 0.052). Meanwhile, VEGFA-rs3025039 was associated with AMS-related headaches. Under the dominant model, the rs3025039 CT/TT genotype showed a 2.08-fold increased risk of AMS-related headaches (P = 0.0007), and the results remained significant after FDR correction (Q = 0.028). No associations were found for EGLN1-rs2153364 and PPARA-rs7292407 with any AMS-related symptoms even before FDR testing. Based upon our findings with AMS-related symptoms, only EPAS1-rs6756667 and VEGFA-rs3025039 were applied for the following analysis. Further information, including the distribution of SNP alleles/genotypes and results under the recessive model, was presented in Additional files 5-8.
Associations between SNPs and the intensity of AMS/AMS-related symptoms
Association analyses between EPAS1-rs6756667, VEGFA-rs3025039 and the intensity of AMS/AMS-related symptoms were performed using the chi-square test. The rs6756667 GG genotype was associated with a significantly increased risk of AMS (P = 0.003), especially mild AMS (P = 0.018). Moreover, the rs6756667 GG genotype showed a significantly increased risk of mild GI symptoms (P = 0.040, Fig. 2a). The rs3025039 CT/TT genotype was associated with a significantly increased risk of AMS (P = 0.017), especially mild AMS (P = 0.019). In addition, the rs3025039 CT/TT genotype exhibited a significantly increased risk of mild headache (P = 0.003, Fig. 2b). However, no significant associations were found between SNPs and moderate-severe AMS/AMS-related symptoms (Fig. 2c-d).
The combined effects of EPAS1-rs6756667 and VEGFA-rs3025039
To explore the combined effects of EPAS1 and VEGFA, we divided the study subjects into 4 subgroups according to the rs6756667 GG and rs3025039 CT/TT status (Table 7). Subjects carrying both the rs6756667 GG and rs3025039 CC genotypes were more likely to develop AMS than those in subgroup 1 (OR = 2.06, 95% CI 1.30–3.27, P = 0.002). Subjects carrying both the rs6756667 GA/AA and rs3025039 CT/TT genotypes had an increased risk of AMS compared with those in subgroup 1 (OR = 2.35, 95% CI 1.07–5.17, P = 0.033). Moreover, compared with non-AMS subjects, AMS patients were more often carriers of both the rs6756667 GG genotype and the rs3025039 CT/TT genotype: 25.6% vs. 19.0% (OR = 2.70; 95% CI 1.59–4.61; P < 0.001). All results remained significant after adjusting for age, BMI and smoking status.
In our present cross-sectional study, the contributions of 23 SNPs in 7 genes involved in sensing low oxygen levels to the risk of AMS were assessed in 604 young, healthy Han Chinese men who rapidly ascended to HA. Of them, we found that 4 SNPs (rs6756667 in EPAS1, rs3025039 in VEGFA, rs7292407 in PPARA and rs2153364 in EGLN1) were associated with the risk of developing AMS even though the results turned insignificant after FDR testing. Further stratified analysis revealed that AMS patients with the rs6756667 GG genotype had a higher risk of developing mild AMS and mild GI symptoms than those with other genotypes. In contrast, the rs3025039 CC genotype was associated with a lower risk of mild AMS and mild headache. These data indicate that EPAS1 and VEGFA SNPs play important but different roles in the physiological effects of AMS in the Chinese Han population. To the best of our knowledge, this is the first epidemiological study to demonstrate such an association with AMS-related symptoms. The results of the current study provide novel information about EPAS1 and VEGFA under conditions of acute hypobaric hypoxia, perhaps revealing a role for these genes in AMS; however, further study is required to confirm this hypothesis.
We showed 4 SNPs (rs6756667 in EPAS1, rs3025039 in VEGFA, rs7292407 in PPARA and rs2153364 in EGLN1) were associated with the risk of developing AMS. Comparing these 4 SNPs under dominant model, EPAS1-rs6756667 showed strongest association with the occurrence of AMS, while the other 3 SNPs showed relatively weak association (P = 0.003 vs. P > 0.01), which implicating at least in our tested gene SNPs, EPAS1 variants might play more important role than other genes on HIF pathway. Furthermore, following the revised diagnosis criteria of LLSS, rs6756667 in EPAS1 also showed positive association with AMS even without the influence of sleeping disturbance . Although after FDR correction, only rs6756667 in EPAS1 retained significance under codominant model. Our results need to be further validated in a larger sample size.
Although our results showed that rs6756667 in EPAS1 correlates significantly with the development of GI symptoms, leading to a higher risk of AMS, the relationship may be difficult to interpret. EPAS1 has 16 exons extending over 90 kb on chromosome 2, with a large first intron (50 kb) . In our present study, all the selected EPAS1 SNPs, rs13419896, rs4953354, rs6756667 and rs1868092, are located in EPAS1 introns, which may affect binding to various transcription factors. The exonic SNPs were not selected because previous studies showed that the EPAS1 SNPs with the greatest differences in frequency between Tibetans and non-Tibetan Chinese individuals are mainly located within introns [8, 27,28,29]. Hypoxia alters genome-wide gene expression at both the transcription and translation levels. Ribosome density at the 5′-UTR significantly increases under hypoxic conditions and is positively correlated with the presence of upstream open reading frames in the 5′-UTR of mRNAs . Interestingly, in a study of native Tibetans residing at different HAs, the rs6756667 A allele was dominant. As the altitude increased, there was a linear decreasing trend in GG and AG carriers and an increasing trend in AA carriers . Studies of Tibetans revealed that this population comprises a mixture of ancestral populations related to the Sherpa and Han Chinese . Thus, it is reasonable to speculate that the different genetic backgrounds between lowlanders, such as Han Chinese, and highlanders, such as Tibetans, as related to EPAS1 rs6756667 might be attributed to adaptation at HA. The rs6756667 A allele may be an allele that increases fitness to HA-related hypoxia. EPAS1 was initially detected in endothelial cells, but expression has also been found in the lung, placenta, kidney, heart, liver, small intestine and many other organs associated with the regulation of oxygen metabolism [33, 34]. EPAS1 was shown in genome-wide association studies to be associated with lower hemoglobin concentration and hypoxia adaptation in Tibetans [27, 35]. The adaptability of the Tibetan and plateau Han populations is mainly due to the strong affinity of hemoglobin for oxygen, which provides sufficient oxygen to tissues and organs . EPAS1 was also shown by whole-exome sequencing to be highly significantly associated with HA-related pulmonary hypertension in Angus cattle . A recent study showed that HIF-1α and EPAS1/HIF-2α are upregulated in intestinal tissues from patients with GI vascular malformations . VEGF is a direct target of EPAS1, and EPAS1 enhances vessel formation by endothelial cells through NOTCH1, Ang2 and DLL4 . Iron and glucose homeostasis in the gut and liver may also be involved in regulating the GI system through the EPAS1 gene under acute hypoxia [39, 40]. However, how EPAS1, and the intronic SNP in particular, affects GI symptoms remains largely unknown.
We also showed that VEGFA rs3025039 is associated with the risk of developing AMS, which is consistent with previous findings [19, 41]. Furthermore, we revealed that VEGFA rs3025039 is associated with headache, the most important symptom for the diagnosis of AMS. The SNP rs3025039 is a C/T single-nucleotide variation denoted as 936 C/T at the beginning of chromosome 6; the affected gene has a 14-kb coding region with 8 exons and 7 introns . Unlike rs6756667 in EPAS1, relationships with the SNP rs3025039 and other VEGFA gene variants have been widely reported for HA illness and for cancer, heart disease, diabetes, neurological diseases and many other types of diseases. It is reasonable to assume that VEGFA gene variants may contribute to variations in VEGF serum levels and to altered susceptibility to diseases related to altered angiogenesis. In Bahraini  and southern Italian subjects , the TT genotype of rs3025039 was associated with lower VEGF levels. However, we could not find any reports on this topic in the Chinese Han population at sea level.
It remains to be determined whether plasma VEGF contributes to the pathogenesis of AMS. VEGF is a selective endothelial cell mitogen that induces vascular permeability and may be involved in increasing blood-brain barrier permeability in rather severe illness. The plasma levels of both VEGF and its soluble receptor (sFlt-1) are significantly increased in HAPE patients compared to non-HAPE controls . Evidence regarding VEGF from AMS studies remains controversial. Absolute values or changes in VEGF were not correlated with AMS scores [46, 47]. On the other hand, Ding et al. reported that among Chinese Han subjects (n = 40), the AMS group had a higher VEGF plasma concentration than the control group at both sea level and HA . In indirect agreement with the results of this previous study, we showed that individuals carrying rs3025039 CT/TT exhibited a higher frequency of developing AMS than those carrying the wild-type CC allele. We speculated that the same SNPs may not have the same contributions to protein production in different races and regions; in support this hypothesis, rs3025039 has different effects on VEGF plasma levels in populations from three villages in Southern Italy . Further investigation is necessary to elucidate the relationship between VEGFA variants and VEGF serum levels and the underlying mechanisms in the Chinese Han population.
This is the first study to show that VEGFA rs3025039 is associated with headache in AMS. However, no direct evidence has been reported regarding the association of VEGF plasma levels with HA-related headache. In a rat model, recurrent headache induced by episodic inflammatory soup blood-brain barrier disruption was accompanied by increased VEGF expression . Recently, Zhao et al.  reported that the rs3025039 TC/TT genotypes significantly increased the risk of poor recovery from ischemia and stroke in Chinese patients, which may indirectly explain why rs3025039 CT/TT carriers had a higher risk of headache and AMS in our study. Although we did not measure serum VEGF levels in our subjects, our results support the hypothesis that VEGF is an important component of the pathogenesis of AMS. The significant association of SNP rs3025039 with AMS needs to be validated in further studies with larger populations.
The present study does have certain limitations. First, not all tag SNPs in the HIF pathway were selected. Second, the number of subjects carrying the EPAS1 rs6756667 AA genotype in this study was quite small (n = 6), and statistical errors may exist; more subjects must be included in future studies to verify the association between SNPs and AMS. Furthermore, our study was limited to young Chinese Han males, and therefore, the results may not be applicable to all Han populations, including individuals of other genders and ages.
Despite the limited power of clinical association studies, the results using the novel LLSS standard have shown for the first time that variants in EPAS1 and VEGFA are related to the susceptibility to AMS in the Chinese Han population. More importantly, EPAS1 and VEGFA variants may affect the risk of AMS through different AMS-related symptoms. Our results provide novel evidence regarding the involvement of genetic factors in the development of AMS. These gene variants have potential utility in screening susceptible populations and predicting the clinical symptoms that lead to AMS, which will benefit soldiers and travelers who rapidly ascend to plateaus.
Availability of data and materials
All data generated or analyzes during this study are included in this published article and its supplementary information files.
Acute mountain sickness
Body mass index
linkage disequilibrium coefficient
Endothelial PAS domain-containing protein 1
Egl nine homolog 1
Factor inhibiting HIF-1α
High-altitude cerebral edema
High-altitude pulmonary edema
Hypoxia-inducible factor 1-α inhibitor
Lake Louise scoring system
Logarithm of odds
- OR :
Prolyl hydroxylase domain protein 2
- SaO2 :
Arterial oxygen saturation
Single nucleotide polymorphism
- SpO2 :
Pulse oxygen saturation
Vascular endothelial growth factor
Hartman-Ksycinska A, Kluz-Zawadzka J, Lewandowski B. High altitude illness. Przegl Epidemiol. 2016;70(3):490–9.
Jin J. Acute Mountain sickness. JAMA. 2017;318(18):1840.
Bartsch P, Swenson ER. Clinical practice: acute high-altitude illnesses. N Engl J Med. 2013;368(24):2294–302.
Luks AM, Swenson ER, Bartsch P. Acute high-altitude sickness. Eur Respir Rev. 2017;26(143):160096.
Simonson TS. Altitude adaptation: a glimpse through various lenses. High Alt Med Biol. 2015;16(2):125–37.
MacInnis MJ, Koehle MS. Evidence for and against genetic predispositions to acute and chronic altitude illnesses. High Alt Med Biol. 2016;17(4):281–93.
Simonson TS, McClain DA, Jorde LB, Prchal JT. Genetic determinants of Tibetan high-altitude adaptation. Hum Genet. 2012;131(4):527–33.
Simonson TS, Yang Y, Huff CD, Yun H, Qin G, Witherspoon DJ, et al. Genetic evidence for high-altitude adaptation in Tibet. Science. 2010;329(5987):72–5.
Bigham AW. Genetics of human origin and evolution: high-altitude adaptations. Curr Opin Genet Dev. 2016;41:8–13.
Lorenzo FR, Huff C, Myllymaki M, Olenchock B, Swierczek S, Tashi T, et al. A genetic mechanism for Tibetan high-altitude adaptation. Nat Genet. 2014;46(9):951–6.
Zimna A, Kurpisz M. Hypoxia-inducible Factor-1 in physiological and pathophysiological angiogenesis: applications and therapies. Biomed Res Int. 2015;2015:549412.
Bigham AW, Lee FS. Human high-altitude adaptation: forward genetics meets the HIF pathway. Genes Dev. 2014;28(20):2189–204.
Kumar H, Choi DK. Hypoxia inducible factor pathway and physiological adaptation: a cell survival pathway? Mediat Inflamm. 2015;2015:584758.
Dengler VL, Galbraith M, Espinosa JM. Transcriptional regulation by hypoxia inducible factors. Crit Rev Biochem Mol Biol. 2014;49(1):1–15.
Wilson MH, Newman S, Imray CH. The cerebral effects of ascent to high altitudes. Lancet Neurol. 2009;8(2):175–91.
MacInnis MJ, Widmer N, Timulsina U, Subedi A, Siwakoti A, Pandit BP, et al. A preliminary genome-wide association study of Acute Mountain sickness susceptibility in a Group of Nepalese Pilgrims Ascending to 4380 m. High Alt Med Biol. 2015;16(4):290–7.
Zhang E, Zhang J, Jin J, Qin J, Li H, Huang L. Variants of the low oxygen sensors EGLN1 and HIF-1AN associated with acute mountain sickness. Int J Mol Sci. 2014;15(12):21777–87.
Guo LI, Zhang J, Jin J, Gao X, Yu J, Geng Q, et al. Genetic variants of endothelial PAS domain protein 1 are associated with susceptibility to acute mountain sickness in individuals unaccustomed to high altitude: a nested case-control study. Exp Ther Med. 2015;10(3):907–14.
Ding H, Liu Q, Hua M, Ding M, Du H, Zhang W, et al. Associations between vascular endothelial growth factor gene polymorphisms and susceptibility to acute mountain sickness. J Int Med Res. 2012;40(6):2135–44.
Buroker NE, Ning XH, Zhou ZN, Li K, Cen WJ, Wu XF, et al. VEGFA SNPs and transcriptional factor binding sites associated with high altitude sickness in Han and Tibetan Chinese at the Qinghai-Tibetan plateau. J Physiol Sci. 2013;63(3):183–93.
Roach RC, Hackett PH, Oelz O, Bartsch P, Luks AM, MacInnis MJ, et al. The 2018 Lake Louise Acute Mountain sickness score. High Alt Med Biol. 2018;19(1):4–6.
Macinnis MJ, Lanting SC, Rupert JL, Koehle MS. Is poor sleep quality at high altitude separate from acute mountain sickness? Factor structure and internal consistency of the Lake Louise score questionnaire. High Alt Med Biol. 2013;14(4):334–7.
Hall DP, MacCormick IJ, Phythian-Adams AT, Rzechorzek NM, Hope-Jones D, Cosens S, et al. Network analysis reveals distinct clinical syndromes underlying acute mountain sickness. PLoS One. 2014;9(1):e81229.
Ge RL, Simonson TS, Gordeuk V, Prchal JT, McClain DA. Metabolic aspects of high-altitude adaptation in Tibetans. Exp Physiol. 2015;100(11):1247–55.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001;125(2001):279–84.
Tian H, McKnight SL, Russell DW. Endothelial PAS domain protein 1 (EPAS1), a transcription factor selectively expressed in endothelial cells. Genes Dev. 1997;11(1):72–82.
Beall CM, Cavalleri GL, Deng L, Elston RC, Gao Y, Knight J, et al. Natural selection on EPAS1 (HIF2alpha) associated with low hemoglobin concentration in Tibetan highlanders. Proc Natl Acad Sci U S A. 2010;107(25):11459–64.
Wang B, Zhang YB, Zhang F, Lin H, Wang X, Wan N, et al. On the origin of Tibetans and their genetic basis in adapting high-altitude environments. PLoS One. 2011;6(2):e17002.
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZX, Pool JE, et al. Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010;329(5987):75–8.
Jiang ZW, Yang JQ, Dai AM, Wang YM, Li W, Xie Z. Ribosome profiling reveals translational regulation of mammalian cells in response to hypoxic stress. BMC Genomics. 2017;18(1):638.
Basang Z, Wang B, Li L, Yang L, Liu L, Cui C, et al. HIF2A variants were associated with different levels of high-altitude hypoxia among native Tibetans. PLoS One. 2015;10(9):e0137956.
Jeong C, Alkorta-Aranburu G, Basnyat B, Neupane M, Witonsky DB, Pritchard JK, et al. Admixture facilitates genetic adaptations to high altitude in Tibet. Nat Commun. 2014;5:3281.
Patel SA, Simon MC. Biology of hypoxia-inducible factor-2alpha in development and disease. Cell Death Differ. 2008;15(4):628–34.
Putra AC, Eguchi H, Lee KL, Yamane Y, Gustine E, Isobe T, et al. The a allele at rs13419896 of EPAS1 is associated with enhanced expression and poor prognosis for non-small cell lung Cancer. PLoS One. 2015;0(8):e0134496.
Tashi T, Scott Reading N, Wuren T, Zhang X, Moore LG, Hu H, et al. Gain-of-function EGLN1 prolyl hydroxylase (PHD2 D4E:C127S) in combination with EPAS1 (HIF-2alpha) polymorphism lowers hemoglobin concentration in Tibetan highlanders. J Mol Med (Berl). 2017;95(6):665–70.
Li C, Li X, Liu J, Fan X, You G, Zhao L, et al. Investigation of the differences between the Tibetan and Han populations in the hemoglobin-oxygen affinity of red blood cells and in the adaptation to high-altitude environments. Hematology. 2018;23(5):309–13.
Newman JH, Holt TN, Cogan JD, Womack B, Phillips JA, Li C, et al. Increased prevalence of EPAS1 variant in cattle with high-altitude pulmonary hypertension. Nat Commun. 2015;6:6863.
Feng N, Chen H, Fu S, Bian Z, Lin X, Yang L, et al. HIF-1alpha and HIF-2alpha induced angiogenesis in gastrointestinal vascular malformation and reversed by thalidomide. Sci Rep. 2016;6:27280.
Mastrogiannaki M, Matak P, Peyssonnaux C. The gut in iron homeostasis: role of HIF-2 under normal and pathological conditions. Blood. 2013;122(6):885–92.
Ramakrishnan SK, Shah YM. A central role for hypoxia-inducible factor (HIF)-2alpha in hepatic glucose homeostasis. Nutr Healthy Aging. 2017;4(3):207–16.
Buroker NE, Ning XH, Zhou ZN, Li K, Cen WJ, Wu XF, et al. AKT3, ANGPTL4, eNOS3, and VEGFA associations with high altitude sickness in Han and Tibetan Chinese at the Qinghai-Tibetan plateau. Int J Hematol. 2012;96(2):200–13.
Renner W, Kotschan S, Hoffmann C, Obermayer-Pietsch B, Pilger E. A common 936 C/T mutation in the gene for vascular endothelial growth factor is associated with vascular endothelial growth factor plasma levels. J Vasc Res. 2000;37(6):443–8.
Al-Habboubi HH, Sater MS, Almawi AW, Al-Khateeb GM, Almawi WY. Contribution of VEGF polymorphisms to variation in VEGF serum levels in a healthy population. Eur Cytokine Netw. 2011;22(3):154–8.
Ruggiero D, Dalmasso C, Nutile T, Sorice R, Dionisi L, Aversano M, et al. Genetics of VEGF serum variation in human isolated populations of cilento: importance of VEGF polymorphisms. PLoS One. 2011;6(2):e16982.
Zhang S, Liu J, Jiang D, Wuren T, Ma S, Du Y, et al. The plasma level changes of VEGF and soluble VEGF receptor-1 are associated with high-altitude pulmonary edema. J Med Invest. 2018;65(1.2):64–8.
Schommer K, Wiesegart N, Dehnert C, Mairbaurl H, Bartsch P. No correlation between plasma levels of vascular endothelial growth factor or its soluble receptor and acute mountain sickness. High Alt Med Biol. 2011;12(4):323–7.
Nilles E, Sayward H, D'Onofrio G. Vascular endothelial growth factor and acute mountain sickness. J Emerg Trauma Shock. 2009;2(1):6–9.
Mi X, Ran L, Chen L, Qin G. Recurrent headache increases blood-brain barrier permeability and VEGF expression in rats. Pain Physician. 2018;21(6):E633–42.
Zhao J, Bai Y, Jin L, Weng Y, Wang Y, Wu H, et al. A functional variant in the 3′-UTR of VEGF predicts the 90-day outcome of ischemic stroke in Chinese patients. PLoS One. 2017;12(2):e0172709.
Authors are grateful to all subjects for participating in this study. Authors thank Can Chen, Meng-Yang Deng, Hua-Li Kang for laboratory technical support. Authors are also grateful for comments from reviewers which are valuable for improving manuscript.
This project was support by the National Natural Science Foundation of China (Grant No. 81730054, Grant No.81873519), the Ministry of Health of P.R. China (Grant No. 201002012) and Research Project of PLA (Grant No. BLJ18J007).
Ethics approval and consent to participate
Ethical approval was obtained from the Ethics Committee of Xinqiao Hospital, Army Medical University (number: 2012014; approved May 9, 2012). Each participant provided written informed consent, and our study conformed to the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Haplotype block maps for SNPs in EGLN1, HIF1A, HIF1AN, PPARA, and VEGFA.
Table S1. Primer sequences for MALDI-TOF MS.
Table S2. Basic SNP information among subjects.
Negative associations between SNP sites and AMS.
Associations between SNPs and AMS related-headache.
Associations between SNPs and AMS related-dizziness/light-headedness.
Associations between SNPs and AMS related-gastrointestinal symptoms.
Associations between SNPs and AMS related-fatigue.
About this article
Cite this article
Zhang, JH., Shen, Y., Liu, C. et al. EPAS1 and VEGFA gene variants are related to the symptoms of acute mountain sickness in Chinese Han population: a cross-sectional study. Military Med Res 7, 35 (2020). https://doi.org/10.1186/s40779-020-00264-6
- Acute mountain sickness
- Single nucleotide polymorphism
- Endothelial PAS domain protein 1 (EPAS1)
- Vascular endothelial growth factor (VEGFA)