Skip to main content
Download PDF

Association of variants m.T16172C and m.T16519C in whole mtDNA sequences with high altitude pulmonary edema in Han Chinese lowlanders

BMC Pulmonary Medicine volume 22, Article number: 72 (2022) Cite this article

Abstract

Background

High altitude pulmonary edema (HAPE) is a hypoxia-induced non-cardiogenic pulmonary edema that typically occurred in un-acclimatized lowlanders, which inevitably leads to life-threatening consequences. Apart from multiple factors involved, the genetic factors also play an important role in the pathogenesis of HAPE. So far, researchers have put more energy into the nuclear genome and HAPE, and ignored the relationship between the mitochondrion DNA (mtDNA) variants and HAPE susceptibility.

Methods

We recruited a total of 366 individuals including 181 HAPE patients and 185 non-HAPE populations through two times. The first time, 49 HAPE patients and 58 non-HAPE individuals were performed through whole mtDNA sequences to search the mutations and haplogroups. The second time, 132 HAPE patients and 127 non-HAPE subjects were collected to apply verifying these mutations and haplogroups of mtDNA with the routine PCR method.

Results

We analyzed and summarized the clinical characteristics and sequence data for the 49 HAPE patients and 58 non-HAPE individuals. We found that a series of routine blood indexes including systolic arterial blood pressure (SBP), heart rate (HR), white blood cell (WBC), and C-reactive protein (CRP) in the HAPE group presented higher and displayed significant differences compared with those in the non-HAPE group. Although the average numbers of variants in different region and group samples were not statistically significant (P > 0.05), the mutation densities of different regions in the internal group showed significant differences. Then we found two mutations (T16172C and T16519C) associated with the HAPE susceptibility, the T16172C mutation increased the risk of HAPE, and the T16519C mutation decreased the HAPE rating. Furthermore, the two mutations were demonstrated with 132 HAPE patients and 127 non-HAPE individuals. Unfortunately, all the haplogroups were not associated with the HAPE haplogroups.

Conclusions

We provided evidence of differences in mtDNA polymorphism frequencies between HAPE and non-HAPE Han Chinese. Genotypes of mtDNA 16172C and 16519C were correlated with HAPE susceptibility, indicating the role of the mitochondrial genome in the pathogenesis of HAPE.

Peer Review reports

Background

High altitude pulmonary edema (HAPE) is an acute idiopathic mountain disease that usually occurs quickly in lowlanders when they are exposed to an altitude exceeding 2500 m above sea level. HAPE is considered a life-threatening non-cardiogenic disease [1]. In some cases, HAPE can develop at the altitude of 1500 - 2500 m for highly susceptible populations [2]. This disease typically occurs 2 - 5 days after arriving at the altitude and is associated with insidiousness [3]. The mortality rate of untreated HAPE is up to 50%, while the mortality rate of treated HAPE is decreased to 11% [3]. Early symptoms of HAPE include non-productive cough, exertional dyspnea, chest pain, and reduced exercise tolerance. Without treatment, HAPE can progress to dyspnea at rest [3]. This disease is caused by a variety of risk factors including hypoxic ventilatory response, rapidly ascent, high altitude, tired, lacking sufficient sleeping, smoking, drinking, heredity, etc. [3,4,5,6]. In addition, the patients with basic disease such as hypertension and diabetes are easier to develop HAPE [7, 8]. Therefore, when studying the relationship between genetics and HAPE, the risk factors should be eliminated.

So far, the major pathogenic mechanism of HAPE is the imbalance and mutations of the pulmonary endothelial cells (ECs) and vascular smooth muscle cells (VSMCs), which leads to the accumulation of extravascular fluid in the lungs and then lung malfunction [9, 10]. A large number of studies have found that many pathways participate in the development of HAPE, such as renin–angiotensin–aldosterone system (RAAS) [11], hypoxia-inducible factor (HIF) pathway [12], the nitric oxide (NO) pathway [13], the mitochondrion DNA (mtDNA) [14, 15], etc. For HAPE patients, suitable treatments include oxygen suction and lowland transferring, which suggests that HAPE is very sensitive to the changes in oxygen/energy metabolism [15,16,17]. Mitochondria are major ATP energy production centers in eukaryotic cells, and their function can be affected by the mtDNA variants [18]. Therefore, the relationship between HAPE and mtDNA is important for understanding the pathogenic mechanism of HAPE.

MtDNA is an extranuclear double-stranded DNA found only in mitochondria. In most eukaryotes, mtDNA is a circular molecule and is inherited from the mother. The full length of mtDNA is 16,569 base pairs (bp), divided into heavy (H) brand and light (L) brand, including 37 genes encoding 2 rRNA (12S rRNA and 16S rRNA), 22 tRNA, and 13 polypeptide subunits (Complex I, III, IV, and V) that produce mitochondrial ATP and participate in the oxidative phosphorylation together with proteins coded by nuclear genome [19]. In addition to coding region, mtDNA also contains a part of the non-coding region called control region or displacement loop (D-loop), which can be divided into high variant sequencing-I (HVS-I) and high variant sequencing-II (HVS-II), and is vital to mtDNA replication and regulation [20]. Mitochondrial function can be affected by variations in mtDNA, including polymorphisms, content changes, and deletions [21]. These variants play an important role in acclimatizing or adapting to hypoxia [21]. In this study, we mainly focused on the association between mtDNA mutants and HAPE susceptibility. For example, Luo et al. used polymerase chain reaction (PCR) methods instead of whole mtDNA sequences to reveal that genotypes of mtDNA 3397G and 3552A were correlated with HAPE susceptibility in Han Chinese [22]. Some results suggested that mitochondrial haplogroups B and M7 were associated with the inadaptability of hypoxic environments, whereas haplogroups G and M9a1a1c1b may be associated with hypoxic adaptation [23]. Recently, literature presented that variants G4491A, A4944G, and A14002G associated with haplogroup M33a2′3 may be the main reasons for the susceptibility of Indian male lowlanders to HAPE, however, the study only included 15 HAPE patients and 20 non-HAPE individuals through long PCR and the results were not verified [15].

Even though these results have already established the association between mtDNA mutations and HAPE, there are still some problems worth exploring and studying, such as limited samples, verification of polymorphisms and haplogroups. In this study, we recruited two groups of people. The first group including 51 HAPE patients (2 patients were removed due to poor quality of mtDNA) and 58 non-HAPE persons was studied using whole mtDNA sequences to find the possible association between mutations and haplogroups. The second group including 132 HAPE patients and 127 non-HAPE patients was used to verify some variants of mtDNA with the routine PCR method. We found two mutations (T16172C and T16519C) associated with the HAPE susceptibility using the largest samples so far. The T16172C mutation indicated an increased susceptibility to HAPE, and the T16519C mutation could reduce the risk factor for obtaining HAPE.

Methods

Samples collections

In this study, a total of 366 subjects were recruited, including 181 HAPE patients and 185 non-HPAE subjects from September 2018 to September 2020, and divided into 2 groups. The first group contained 51 HAPE patients and 58 non-HAPE people who were low-altitude Han-Chinese and climbed rapidly to 3658 m in Lhasa, Tibet, China within 7 days (Table 1 and Additional file 1). In addition, the first group of subjects were male and met the following conditions: (1) HAPE patients were diagnosed according to the clinical standard requirements, which mainly included cough, dyspnea at rest, white or pink foamy sputum, central cyanosis, pulmonary crackles, and the presence of flake or cloud infiltrate shadows at unilateral or bilateral pulmonary hilar on chest X-ray [24, 25]. (2) They were from the General Hospital of Tibetan Military Command and had no relationship with each other. (3) They did not have a history of smoking, drinking, hypertension, diabetes, cardiopulmonary disease, acute or chronic pulmonary infection (e.g., pneumonia), and other mtDNA-related diseases. (4) The non-HAPE had no acute mountain sickness after moving to Lhasa for 7 days. The second group had 132 HAPE patients and 127 non-HAPE individuals, and also met the above screening criteria (Additional file 1).

Table 1 Basic clinical characteristics of all the subjects
Full size table

DNA extraction, sequencing, and haplogroup classification

Blood samples were collected from all the subjects and stored at − 80 ℃. Genomic DNA was extracted from peripheral blood using Magbead Blood DNA Kit (CWBio, Beijing, China). For each patient, 200 ng genomic DNA was sheared by Biorupter (Diagenode, Belgium) to acquire 150 - 200 bp fragments. The ends of the DNA fragment were repaired, and Illumina Adaptor was added (Fast Library Prep Kit, iGeneTech, Beijing, China). After the sequencing library was constructed, the whole exons were captured with AI-Mito-Cap (iGeneTech, Beijing, China) and sequenced on the Illumina platform (Illumina, San Diego, CA), with 150 base paired‐end reads. Raw reads were filtered by FastQC to remove low-quality reads. Then clean reads were mapped to the reference genome GRCh37 using Bwa. After duplications were removed, SNV and InDel were called and annotated by GATK, SamtoolS, Varscan. According to the PhyloTree Build 17 standard, using the software development tool MitoTool, the haplogroups of entire mtDNA sequences were constructed using the phylogenetic analysis method [26,27,28]. When evaluating the mtDNA variants, a series of affecters were considered, including the different variants of the particular branches of each haplogroup, the calescence time of variation distribution, and the location of protein- or RNA-encoded gene-based substitutions. All these factors were integrated into phylogenetic analysis to construct Haplogroup.

PCR conditions

We directly sequenced the PCR products to genotype the A263G, T310N, 310-311insC, T16172C, and T16519C polymorphisms in the remaining individuals. The primers (forward: 5’-CAGCCACCATGAATATTGTACG-3’, reverse: 5’-GTTAGGCTGGTGTTAGGGTTC-3’) were designed using the Primer 5.0 software (PREMIER Biosoft International, Palo Alto, CA, USA). The reactions were performed in 50 μL volume containing 50 ng template DNA, 5 mM dNTP (Takara, Dalian, China), 10 pmol forward primers, 10 pmol reverse primers, 2.5 U Taq DNA polymerase enzyme (Takara) in Taq buffer, and 75 mM MgCl2. Then PCR amplification was performed under the following conditions: a pre-denaturation cycle at 94 °C for 4 min; 45 cycles at 94 °C for 20 s, 56 °C for 30 s, 72 °C for 1 min; and a final extension cycle at 72 °C for 10 min. The obtained products were cooled to 4 °C. The PCR products were directly sequenced by Shanghai Sangon Biological Engineering Technology & Services Co., Ltd (Shanghai, China). The sequenced results were assigned with the standard mitochondrial reference genome called Revised Cambridge reference sequence (rCRS). To verify the variants A3397G and T3552A, we used the primers (forward: 5’- GTCCTACGTGATCTGAGTTCAG-3’, reverse: 5’-GCTAGGGTGACTTCATATGAG-3’).

Ethical statement

All the samples of this study collected had written informed content from patients. All clinical and experimental steps were approved by the Ethical Committee of the Seventh Medical Center of PLA General Hospital and the Ethical Committee of the General Hospital of Tibetan Military Command. Meanwhile, these contents were conducted in accordance with the Helsinki Declaration guideline.

Statistical analyses

All the statistical analyses of data were performed using SPSS 17.0. Pearson χ2 test was used to assess the significance of differences in haplogroup and SNP frequencies between HAPE and non-HAPE groups. Quantitative data were expressed as mean ± standard deviation (mean ± SD) and the unpaired student’s t-tests were used to compare the HAPE group and non-HAPE group. The P-values, odds ratios (ORs), and 95% confidence intervals (95% CI) were calculated. A two-sided P-value < 0.05 was considered statistically significant. All P-values were adjusted by Bonferroni correction.

Results

Clinical characteristics of all the individuals

Firstly, in the first group, 2 HAPE patients were removed because of the poor-quality DNA levels. Then, a total of 49 HAPE patients and 58 non-HAPE individuals were further analyzed. They had same screening conditions and suffered similar environments. Their detailed clinical information is listed in Table 1 and Additional file 1.

As shown in Table 1, the age distribution was consistent between the HAPE patients and non-HAPE individuals. Except for red blood cell (RBC), all the routine blood indexes such as systolic arterial blood pressure (SBP), white blood cell (WBC), and C-reactive protein (CRP) in HAPE patients were significantly different from non-HAPE (P < 0.05), which was probably related to the adaption of high environments [1, 29, 30].

Data summary of sequencing and quality control

The mtDNA of the first group samples were isolated and sequenced using the Illumina platform. Nearly 92% of total data generated a Q > 30 Phred score. In addition, after removing the duplications, the mean sequence depth of 8400 × was observed, which accounted for 24.6% on average. The whole mitochondrial genome region was covered for all the samples. All the variants of HAPE patients and non-HAPE individuals are listed in Additional file 2 and Additional file 3. After processing these variants, the detailed and summarized results are shown in Figs. 1 and 2, respectively.

Fig. 1
figure 1

Distribution of variants in HAPE and non-HAPE groups. a Number of variants per sample. b Number of variants in coding regions and non-coding regions in per sample. c Number and proportion of mutations in different regions. Left, HAPE group; Right, non-HAPE group. P-value was obtained in t-test; *, P-value < 0.05; **, P-value < 0.01

Full size image
Fig. 2
figure 2

Gene distributions in mitochondrial genome of HAPE and non-HAPE groups. a Mutation densities (mutation number in a region/region bp length) in coding regions and non-coding regions of HAPE or non-HAPE group. b Number and proportion of frameshift variants, missense variants, stop gained, stop lost, synonymous variants and upstream gene variants. c Mutation rates in HAPE or on-HAPE individuals. P-value was obtained in t-test; *, P-value < 0.05; **, P-value < 0.01

Full size image

Compared to the non-HAPE group, the mean numbers of variants in HAPE patients were higher; however, the difference between non-HAPE and HAPE groups was not significant (P = 0.799) (Fig. 1a). In HAPE group, a total of 1945 mutations were distributed in non-coding region (602, 30.95%), tRNA region (20, 1.03%), rRNA region (255, 13.11%), complex I (527, 27.10%), complex III (237, 12.19%), complex IV (173, 8.89%), and complex V (131, 6.74%) (Fig. 1b). There were 2284 mutations in non-HAPE individuals, which were located in non-coding region (767, 33.58%), tRNA region (18, 0.79%), rRNA region (299, 13.09%), complex I (581, 25.44%), complex III (260, 11.38%), complex IV (214, 9.37%), and complex V (145, 6.35%) (Fig. 1b). Unfortunately, the two randomly corresponding regions did not show significant differences (Additional file 4). Then, we analyzed the numbers of variants of coding and non-coding regions in HAPE and non-HAPE groups. The results indicated that the differences between HAPE and non-HAPE groups were not significant (Fig. 1c). However, the average variants between the coding region and the non-coding region in HAPE patients were significantly different, which may be related to the length of coding regions and non-coding regions in the sequence. Similar findings were found in the non-HAPE group (Fig. 1c).

To further explore the variants’ distributions characteristics of the mitochondrial genome in HAPE and non-HAPE patients, we used the mutation densities (mutation number in a region/region bp length) to describe the variants distributions in coding regions and non-coding regions. The results indicated that the mutation densities in coding regions or non-coding regions were consistent between the HAPE and non-HAPE groups (Fig. 2a). However, within the same group (HAPE or non-HAPE group), the mutation densities in the non-coding regions were higher than those in the coding regions (Fig. 2a). This suggested the sequences of non-coding regions were inclined to mutations. In all the variants, the top three molecular consequences in both the HAPE group and non-HAPE group were 44.06% upstream gene, 31.36% synonymous, and 23.55% missense variants (Fig. 2b), and there was no significant difference between the HAPE group and non-HAPE group (Additional file 5). Interestingly, a few positions in the mitochondrial genome had very high mutation rates (variants number/samples number), such as 73rd position (1.0), 2706th position (1.0), 8860th position (1.0) in HAPE patients (Fig. 2c and Additional file 6). Meanwhile, these positions in the non-HAPE group also showed similar results and the difference between HAPE and non-HAPE groups was not significant (Fig. 2c and Additional file 6). The results suggested that the mitochondrial genome had some high-frequency mutation sites that were not associated with HAPE.

Association of mtDNA variants with HAPE susceptibility

In this study, a total of 1945 mutations and 2284 mutations were observed in HAPE and non-HAPE groups, respectively, which were located in 483 variants’ sites and 530 variants’ sites (combined with 742 variants’ sites), respectively (Additional file 6). There were 18 different and typical mutations in Table 2, including 5 variants sites (A263G, T310N, (310–311) insC, T16172C, and T16519C) with significant differences between both groups, 6 variants sites (A8459G, C8684T, C14067T, T14470C, and A16164G) only found in HAPE group, 3 variants sites (C194T, C6960T, and G16390A) only found in the non-HAPE group, and 4 variants sites (T11944C, A15235G, C16291T, and A16316G) had high odds ratio (OR) values (Table 2). After Bonferroni correction, only 4 mutations (T310N, (310–311) insC, T16172C, and T16519C) were statistically significant between HAPE and non-HAPE groups (P < 0.05), and they were located in the D-loop of mtDNA non-coding regions. Since the 310 site was short repeats sequences C bases, this site was removed in further analysis. Our results showed that the frequencies of T16172C in the HAPE group (20.41%) were higher than those in the non-HAPE group (5.17%, p-adjust = 0.035, OR = 4.701, and 95%CI = 1.214–18.204). In contrast, the frequencies of T16519C in HAPE patients (34.69%) were significantly lower than those in non-HAPE individuals (67.24%, p-adjust = 0.002, OR = 0.259, and 95%CI = 0.116–0.578) (Table 2). The two mutations were previously reported in other studies, but not in high-altitude diseases [31,32,33]. Then, we also statistically analyzed the number of T16172C mutation and T16519C mutation existing in a certain body at the same time. there was no difference between HAPE patients and non-HAPE individuals (Additional file 7).

Table 2 Comparison of variant frequencies between HAPE and non-HAPE groups
Full size table

Haplogroup comparison between HAPE and non-HAPE patients

A total of 107 complete mtDNA sequences of the first group’s patients were analyzed using MitoTool software according to the PhyloTree Build 17 criteria (GenBank J01415.2) [26,27,28]. We found that all the mtDNA sequences were mapped to the single initial haplogroup L3, which originated from Africa and was composed of two macrohaplogroups, M and N (Table 3). All the subjects belonged to 16 haplogroups. No differences were observed in the haplogroups between the HAPE and non-HAPE groups (P > 0.05, Table 3). These haplogroups were further divided into the next-level haplogroups, which were not statistically significant (P > 0.05, Additional file 8).

Table 3 Distribution of haplogroups in HAPE and non-HAPE groups
Full size table

Validation of two mutations (T16172C and T16519C)

In order to validate the association between two mutations (T16172C and T16519C) and HAPE susceptibility, the subjects of the second group were collected (Additional file 1). Their basic clinical characteristics were listed in Additional file 9. The PCR results of these variants are shown in Table 4. The frequencies of T16172C in HAPE patients (44, 33.33%) were higher than those in non-HAPE individuals (27, 21.26%, p-adjust = 0.042, OR = 1.852, and 95%CI = 1.060–3.236). The frequencies of T16519C in HAPE group (52, 39.39%) were lower than the non-HAPE group (67, 52.76%, p-adjust = 0.042, OR = 0.582, and 95%CI = 0.355–0.953). The above results indicated that the T16172C mutation increased the HAPE susceptibility and the T16519C mutantation decreased the impact from HAPE. Meanwhile, we also analyzed the number of 16172C-16519C mutation between HAPE patients and non-HAPE individuals in the second group and all the subjects. Similar to the result of the first group, there was no difference (Additional file 7).

Table 4 Validation of association between mt16172C, mt16519C genotypes and HAPE susceptibility
Full size table

Discussion

HAPE is a severe acute disease caused by high-altitude hypoxia. Several studies have found that the mutations of mtDNA associated with the HAPE susceptibility; however, these results can’t be verified by other samples [22, 34]. In this study, we found two mutations (T16172C, T16519C) located in the D-loop of the non-coding region were associated with HAPE susceptibility in Han Chinese, which were validated with PCR method in the other group individuals. The frequencies of T16172C mutation in the HAPE group were higher than those in the non-HAPE group. In contrast, the 16519C mutations tended to reduce the HAPE susceptibility. Two mutations were also reported previously in the other diseases, but not in high-altitude diseases [31,32,33].

Previous study had found that the 16172 (rs2853817) and 16519 (rs3937033) mutations showed preliminary association with mitochondrial function without changing primary sequences of protein [35]. The site of T16172C mutation was located in the hypervariable region of mtDNA. So far, little is known about the T16172C mutation in other diseases including mountain sickness, except for one paper that demonstrated the relationship between the T16172C and ataxia telangiectasia for twins [33]. The site of mtDNA 16519 has been reported to have the highest worldwide mutation rates [35], and has an association with some diseases. Navaglia et al. found that the T allele of mtDNA 16519 SNP was correlated with shorter life expectancy in pancreatic cancer [31]. The relationship between 16519 site and other cancers such as breast cancer and familial nasopharyngeal carcinoma has also been demonstrated [36, 37]. Liao et al. reported that the T16519C mutation showed a susceptible tendency to type-2 diabetes mellitus (T2DM) in the Chinese Han population [38]. Significant differences were also found in the change rate of VO (2 max) and citrate synthase activity as a result of training between the two groups at 16519 site [39]. Besides, this mutation might cause other common disorders such as migraine headaches or psychiatric disorders such as schizophrenia [32, 40, 41]. So far, the T16172C and T16519C variants in mtDNA were first reported to be associated with the HAPE of mountain sickness. Although both the 16172 and 16519 variants are located in the D-loop region of mtDNA, they may act by changing the transcription levels of mitochondrial proteins, which are related to oxidative phosphorylation, and the changes in this process may lead to β-cell failure or ATP production lack [42, 43]. These phenomena need to be further verified in future research.

About ten years ago, Luo et al. [22] found that genotypes of mtDNA 3397G and 3552A were correlated with HAPE susceptibility, which were not verified in this paper (Additional file 10). As we all know, HAPE is caused by caused by multiple factors including environmental and genetic reasons, such as hypoxic ventilatory response, rapid ascent, high altitude, tiredness, lacking sufficient sleeping, smoking, drinking, heredity, etc. [3,4,5,6]. In addition, the patients with basic diseases such as hypertension and diabetes are easier to develop HAPE [7, 8]. Thus, the different clinical samples could lead to different results.

Besides variant genotypes of mtDNA, the mitochondrial haplogroup classification in HAPE and non-HAPE individuals was performed using MitoTool. mtDNA is a well-known genetic marker due to the characteristics of high mutation rate, maternal inheritance, high copy number, and lack of recombination, which makes it different from the nuclear genome [34]. On this basis, mtDNA haplogroups such as the non-recombining region of Y-chromosome are usually considered as the main measures for investigating human evolution and origin [44, 45]. Mitochondrial DNA haplogroups are a characteristic cluster of tightly linked mtDNA polymorphism that forms continent-specific genotypes [44]. Studies have also shown that certain mitochondrial haplogroups are predisposed to high-altitude diseases while others can prevent mountain sickness. For example, the studies on the mitochondrial genome had revealed that there was a significant difference in the frequency of the mitochondrial 3010G-3970C haplotype between high- and low-altitude populations, which was believed to be associated with the improved adaptability of the Tibetan population to the low-oxygen environment [46]. Other studies had also shown that the mitochondrial haplogroups were associated with high-altitude adaptation, such as haplogroups B, M7, M9a1a1c1b, B4b, D4, etc. [23, 34, 47]. These mtDNA haplogroups were also observed in our samples; however, we found that they were not significantly different between HAPE and non-HAPE groups (Additional file 8). This phenomenon was inconsistent with the results in the previous study that mtDNA haplogroup M33a2′3 may be linked with the HAPE susceptibility of 15 HAPE patients and 20 non-HAPE individuals who were male and came from Indian [15].

Although we have found two mutations (T16172C, T16519C) associated with the HAPE susceptibility, this study still has limitations, including relatively small sample size and lack of functional experiments to identify the effects of the genetic variants on gene structure/function. Therefore, our findings need to be validated in a larger population and the functional significance of the important genetic variants should be further investigated.

Conclusions

In this study, we established the associations between mtDNA mutations and HAPE susceptibility. The T16172C mutation increased the risk of HAPE, while the variant of T16519C tended to be resistant to HAPE. Then, these two mutations were verified with the other group of patients. We also found haplogroup of mtDNA may not play a role in the HAPE susceptibility. Meanwhile, we found that a few routine blood indexes, such as SBP, HR, WBC, and CRP, showed significant differences between the HAPE group and the non-HAPE group. So, these findings would enrich the association between mtDNA and HAPE and can provide a new perspective to reveal the genetic mechanism of HAPE.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Abbreviations

bp:

Base pairs

CI:

Confidence intervals

CRP:

C-reactive protein

D-loop:

Displacement loop

ECs:

Endothelial cells

H:

Heavy

HAPE:

High altitude pulmonary edema

HIF:

Hypoxia-inducible factor

HR:

Heart rate

HVS-I:

High variant sequencing-I

HVS-II:

High variant sequencing-II

L:

Light

mtDNA:

Mitochondrion DNA

NO:

Nitric oxide

ORs:

Odds ratios

PCR:

Polymerase chain reaction

RAAS:

Renin angiotensin aldosterone system

RBC:

Red blood cell

rCRS:

Revised Cambridge reference sequence

SBP:

Systolic arterial blood pressure

T2DM:

Type-2 diabetes mellitus

VSMCs:

Vascular smooth muscle cells

WBC:

White blood cell

References

  1. Swenson ER, Bartsch P. High-altitude pulmonary edema. Compr Physiol. 2012;2(4):2753–73.

    Article  PubMed  Google Scholar 

  2. Burtscher M, Mairer K, Wille M, Gatterer H, Ruedl G, Faulhaber M, Sumann G. Short-term exposure to hypoxia for work and leisure activities in health and disease: which level of hypoxia is safe? Sleep Breath. 2012;16(2):435–42.

    Article  PubMed  Google Scholar 

  3. Jensen JD, Vincent AL: High altitude pulmonary edema (HAPE). In: StatPearls. edn. Treasure Island (FL); 2020.

  4. Honigman B, Theis MK, Koziol-McLain J, Roach R, Yip R, Houston C, Moore LG, Pearce P. Acute mountain sickness in a general tourist population at moderate altitudes. Ann Intern Med. 1993;118(8):587–92.

    Article  CAS  PubMed  Google Scholar 

  5. Basnyat B, Murdoch DR. High-altitude illness. Lancet. 2003;361(9373):1967–74.

    Article  PubMed  Google Scholar 

  6. Thapa SS, Basnyat B. Higher ascent, trouble breathing: High altitude pulmonary edema (HAPE). Pan Afr Med J. 2018;30:43.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Bartsch P, Maggiorini M, Ritter M, Noti C, Vock P, Oelz O. Prevention of high-altitude pulmonary edema by nifedipine. N Engl J Med. 1991;325(18):1284–9.

    Article  CAS  PubMed  Google Scholar 

  8. Paralikar SJ. High altitude pulmonary edema-clinical features, pathophysiology, prevention and treatment. Indian J Occup Environ Med. 2012;16(2):59–62.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Dunham-Snary KJ, Wu D, Sykes EA, Thakrar A, Parlow LRG, Mewburn JD, Parlow JL, Archer SL. Hypoxic pulmonary vasoconstriction: from molecular mechanisms to medicine. Chest. 2017;151(1):181–92.

    Article  PubMed  Google Scholar 

  10. Lumb AB, Slinger P. Hypoxic pulmonary vasoconstriction: physiology and anesthetic implications. Anesthesiology. 2015;122(4):932–46.

    Article  CAS  PubMed  Google Scholar 

  11. Qi Y, Niu W, Zhu T, Zhou W, Qiu C. Synergistic effect of the genetic polymorphisms of the renin-angiotensin-aldosterone system on high-altitude pulmonary edema: a study from Qinghai-Tibet altitude. Eur J Epidemiol. 2008;23(2):143–52.

    Article  CAS  PubMed  Google Scholar 

  12. Soree P, Gupta RK, Singh K, Desiraju K, Agrawal A, Vats P, Bharadwaj A, Baburaj TP, Chaudhary P, Singh VK, et al. Raised HIF1alpha during normoxia in high altitude pulmonary edema susceptible non-mountaineers. Sci Rep. 2016;6:26468.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Beall CM, Laskowski D, Erzurum SC. Nitric oxide in adaptation to altitude. Free Radic Biol Med. 2012;52(7):1123–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Gu M, Dong X, Shi L, Shi L, Lin K, Huang X, Chu J. Differences in mtDNA whole sequence between Tibetan and Han populations suggesting adaptive selection to high altitude. Gene. 2012;496(1):37–44.

    Article  CAS  PubMed  Google Scholar 

  15. Sharma S, Singh S, Gupta RK, Ganju L, Singh SB, Kumar B, Singh Y. Mitochondrial DNA sequencing reveals association of variants and haplogroup M33a2’3 with High altitude pulmonary edema susceptibility in Indian male lowlanders. Sci Rep. 2019;9(1):10975.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Fuhrmann DC, Brune B. Mitochondrial composition and function under the control of hypoxia. Redox Biol. 2017;12:208–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Kang L, Zheng HX, Zhang M, Yan S, Li L, Liu L, Liu K, Hu K, Chen F, Ma L, et al. MtDNA analysis reveals enriched pathogenic mutations in Tibetan highlanders. Sci Rep. 2016;6:31083.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Cawthon RM. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucl Acids Res. 2009;37(3):e21.

    Article  PubMed  PubMed Central  Google Scholar 

  19. St John JC, Sakkas D, Barratt CL. A role for mitochondrial DNA and sperm survival. J Androl. 2000;21(2):189–99.

    CAS  PubMed  Google Scholar 

  20. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999;23(2):147.

    Article  CAS  PubMed  Google Scholar 

  21. Luo Y, Yang X, Gao Y. Mitochondrial DNA response to high altitude: a new perspective on high-altitude adaptation. Mitochondrial DNA. 2013;24(4):313–9.

    Article  CAS  PubMed  Google Scholar 

  22. Luo Y, Gao W, Chen Y, Liu F, Gao Y. Rare mitochondrial DNA polymorphisms are associated with high altitude pulmonary edema (HAPE) susceptibility in Han Chinese. Wilderness Environ Med. 2012;23(2):128–32.

    Article  PubMed  Google Scholar 

  23. Li Q, Lin K, Sun H, Liu S, Huang K, Huang X, Chu J, Yang Z. Mitochondrial haplogroup M9a1a1c1b is associated with hypoxic adaptation in the Tibetans. J Hum Genet. 2016;61(12):1021–6.

    Article  CAS  PubMed  Google Scholar 

  24. Leon-Velarde F, Maggiorini M, Reeves JT, Aldashev A, Asmus I, Bernardi L, Ge RL, Hackett P, Kobayashi T, Moore LG, et al. Consensus statement on chronic and subacute high altitude diseases. High Alt Med Biol. 2005;6(2):147–57.

    Article  PubMed  Google Scholar 

  25. West JB. English translation of "Nomenclature, classification, and diagnostic criteria of high altitude disease in China". High Alt Med Biol. 2010;11(2):169–72.

    Article  PubMed  Google Scholar 

  26. Fan L, Yao YG. MitoTool: a web server for the analysis and retrieval of human mitochondrial DNA sequence variations. Mitochondrion. 2011;11(2):351–6.

    Article  CAS  PubMed  Google Scholar 

  27. Fan L, Yao YG. An update to MitoTool: using a new scoring system for faster mtDNA haplogroup determination. Mitochondrion. 2013;13(4):360–3.

    Article  CAS  PubMed  Google Scholar 

  28. van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009;30(2):E386-394.

    Article  PubMed  Google Scholar 

  29. Strapazzon G, Vezzaro R, Hofer G, Dal Cappello T, Procter E, Balkenhol K, Platzgummer S, Brugger H. Factors associated with B-lines after exposure to hypobaric hypoxia. Eur Heart J Cardiovasc Imaging. 2015;16(11):1241–6.

    Article  PubMed  Google Scholar 

  30. Richalet JP, Larmignat P, Poitrine E, Letournel M, Canoui-Poitrine F. Physiological risk factors for severe high-altitude illness: a prospective cohort study. Am J Respir Crit Care Med. 2012;185(2):192–8.

    Article  PubMed  Google Scholar 

  31. Navaglia F, Basso D, Fogar P, Sperti C, Greco E, Zambon CF, Stranges A, Falda A, Pizzi S, Parenti A, et al. Mitochondrial DNA D-loop in pancreatic cancer: somatic mutations are epiphenomena while the germline 16519 T variant worsens metabolism and outcome. Am J Clin Pathol. 2006;126(4):593–601.

    Article  CAS  PubMed  Google Scholar 

  32. Sequeira A, Martin MV, Rollins B, Moon EA, Bunney WE, Macciardi F, Lupoli S, Smith EN, Kelsoe J, Magnan CN, et al. Mitochondrial mutations and polymorphisms in psychiatric disorders. Front Genet. 2012;3:103.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Bouhlal Y, Martinez S, Gong H, Dumas K, Shieh JT. Twin mitochondrial sequence analysis. Mol Genet Genomic Med. 2013;1(3):174–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Luo YJ, Gao WX, Li SZ, Huang XW, Chen Y, Liu FY, Huang QY, Gao YQ. Mitochondrial haplogroup D4 confers resistance and haplogroup B is a genetic risk factor for high-altitude pulmonary edema among Han Chinese. Genet Mol Res. 2012;11(4):3658–67.

    Article  CAS  PubMed  Google Scholar 

  35. Soares P, Ermini L, Thomson N, Mormina M, Rito T, Rohl A, Salas A, Oppenheimer S, Macaulay V, Richards MB. Correcting for purifying selection: an improved human mitochondrial molecular clock. Am J Hum Genet. 2009;84(6):740–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Bai RK, Leal SM, Covarrubias D, Liu A, Wong LJ. Mitochondrial genetic background modifies breast cancer risk. Cancer Res. 2007;67(10):4687–94.

    Article  CAS  PubMed  Google Scholar 

  37. Peng Z, Xie C, Wan Q, Zhang L, Li W, Wu S. Sequence variations of mitochondrial DNA D-loop region are associated with familial nasopharyngeal carcinoma. Mitochondrion. 2011;11(2):327–33.

    Article  CAS  PubMed  Google Scholar 

  38. Liao WQ, Pang Y, Yu CA, Wen JY, Zhang YG, Li XH. Novel mutations of mitochondrial DNA associated with type 2 diabetes in Chinese Han population. Tohoku J Exp Med. 2008;215(4):377–84.

    Article  CAS  PubMed  Google Scholar 

  39. Murakami H, Ota A, Simojo H, Okada M, Ajisaka R, Kuno S. Polymorphisms in control region of mtDNA relates to individual differences in endurance capacity or trainability. Jpn J Physiol. 2002;52(3):247–56.

    Article  CAS  PubMed  Google Scholar 

  40. Zaki EA, Freilinger T, Klopstock T, Baldwin EE, Heisner KR, Adams K, Dichgans M, Wagler S, Boles RG. Two common mitochondrial DNA polymorphisms are highly associated with migraine headache and cyclic vomiting syndrome. Cephalalgia. 2009;29(7):719–28.

    Article  CAS  PubMed  Google Scholar 

  41. Hjelm BE, Rollins B, Mamdani F, Lauterborn JC, Kirov G, Lynch G, Gall CM, Sequeira A, Vawter MP. Evidence of mitochondrial dysfunction within the complex genetic etiology of schizophrenia. Mol Neuropsychiatry. 2015;1(4):201–19.

    PubMed  PubMed Central  Google Scholar 

  42. Maassen JA, Janssen GM, t Hart LM. Molecular mechanisms of mitochondrial diabetes (MIDD). Ann Med. 2005;37(3):213–21.

    Article  CAS  PubMed  Google Scholar 

  43. Simmons RA, Suponitsky-Kroyter I, Selak MA. Progressive accumulation of mitochondrial DNA mutations and decline in mitochondrial function lead to beta-cell failure. J Biol Chem. 2005;280(31):28785–91.

    Article  CAS  PubMed  Google Scholar 

  44. Ayala FJ. The myth of Eve: molecular biology and human origins. Science. 1995;270(5244):1930–6.

    Article  CAS  PubMed  Google Scholar 

  45. Underhill PA, Kivisild T. Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet. 2007;41:539–64.

    Article  CAS  PubMed  Google Scholar 

  46. Luo Y, Gao W, Liu F, Gao Y. Mitochondrial nt3010G-nt3970C haplotype is implicated in high-altitude adaptation of Tibetans. Mitochondrial DNA. 2011;22(5–6):181–90.

    Article  CAS  PubMed  Google Scholar 

  47. Li FX, Ji FY, Zheng SZ, Yao W, Xiao ZL, Qian GS. MtDNA haplogroups M7 and B in southwestern Han Chinese at risk for acute mountain sickness. Mitochondrion. 2011;11(4):553–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would thank the patients for their participation in this study. We also thank the teachers for providing critical advice in data analysis. This project has received funding from Military Logistics Research Projects (CLJ16J002).

Funding

This project has received funding from Military Logistics Research Projects (CLJ16J002).

Author information

Author notes

  1. Yan Wang, Xuewen Huang and Fujun Peng have contributed equally to this work

Authors and Affiliations

  1. Clinical Biobank Center, Medical Innovation Research Division of Chinese, PLA General Hospital, No. 28 Fu Xin Road, Hai Dian District, Beijing, 100853, China

    Yan Wang

  2. BaYi Children’s Hospital, The Seventh Medical Center of PLA General Hospital, No.5 Nan Men Cang, Dong Cheng District, Beijing, 100700, China

    Yan Wang, Huiling Han, Yanan Gu, Xin Liu & Zhichun Feng

  3. The Mountain Sickness Prevention Research Center of the General Hospital of Tibet Military Command, Tibet, China

    Xuewen Huang

  4. School of Basic Medical Sciences, Weifang Medical University, Weifang, Shandong, China

    Fujun Peng

Authors
  1. Yan Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xuewen Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fujun Peng
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Huiling Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yanan Gu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xin Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Zhichun Feng
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YW and ZCF conceived the project, supervised the study, and provided guidance in bioinformatics analysis. YW and FJP analyzed the data and wrote the manuscript drafting. XWH and HLH collected the samples, selected the samples, and analyzed the clinical information. YNG and XL performed PCR validation for mutations. YW, XWH, and ZCF revised the manuscript for scholarly publication. All authors reviewed and approved the final manuscript as submitted.

Corresponding authors

Correspondence to Yan Wang or Zhichun Feng.

Ethics declarations

Ethics approval and consent to participate

All clinical and experimental steps were approved by the Ethical Committee of the Seventh Medical Center of PLA General Hospital and the Ethical Committee of the General Hospital of Tibetan Military Command. All samples were obtained with written informed consent from populations including HAPE and non-HAPE persons.

Consent for publication

Consent for publication of all the individuals was obtained for the HAPE patients and non-HAPE persons including clinical data and mtDNA sequences results. And all authors have read and approved the submission of the manuscript.

Competing interests

The authors declare no competing financial interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Basic clinical information of all the patients.

Additional file 2

: The mtDNA sequences results of all the HAPE patients.

Additional file 3

: The mtDNA sequences results of all the non-HAPE persons.

Additional file 4

: Number and proportion of mutations in different regions of mtDNA.

Additional file 5

: Number and proportion of molecular consequence in mtDNA.

Additional file 6

: Comparison of variant frequencies between HAPE and non-HAPE groups.

Additional file 7

: Validation of association between 16172C-16519C genotypes and HAPE susceptibility.

Additional file 8

: Distribution of haplogroups in HAPE and non-HAPE groups.

Additional file 9

: Basic clinical characteristics of all the validating subjects.

Additional file 10

: Validation of association between mt3397G, mt3552A genotypes and HAPE susceptibility.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Huang, X., Peng, F. et al. Association of variants m.T16172C and m.T16519C in whole mtDNA sequences with high altitude pulmonary edema in Han Chinese lowlanders. BMC Pulm Med 22, 72 (2022). https://doi.org/10.1186/s12890-021-01791-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-021-01791-1

Keywords

BMC Pulmonary Medicine

ISSN: 1471-2466

Contact us