Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Interactions among Vascular-Tone Modulators Contribute to High Altitude Pulmonary Edema and Augmented Vasoreactivity in Highlanders

  • Zahara Ali,

    Affiliations Institute of Genomics and Integrative Biology, Delhi, India, Department of Biotechnology, University of Pune, Pune, India

  • Aastha Mishra,

    Affiliations Institute of Genomics and Integrative Biology, Delhi, India, Department of Biotechnology, University of Pune, Pune, India

  • Rahul Kumar,

    Affiliation Institute of Genomics and Integrative Biology, Delhi, India

  • Perwez Alam,

    Affiliations Institute of Genomics and Integrative Biology, Delhi, India, Department of Biotechnology, University of Pune, Pune, India

  • Priyanka Pandey,

    Affiliations Institute of Genomics and Integrative Biology, Delhi, India, Department of Biotechnology, University of Pune, Pune, India

  • Rekhbala Ram,

    Affiliation Institute of Genomics and Integrative Biology, Delhi, India

  • Tashi Thinlas,

    Affiliation Department of Medicine, SNM Hospital, Leh, Ladakh, Jammu and Kashmi, India

  • Ghulam Mohammad,

    Affiliation Department of Medicine, SNM Hospital, Leh, Ladakh, Jammu and Kashmi, India

  • M. A. Qadar Pasha

    qpasha@igib.res.in

    Affiliations Institute of Genomics and Integrative Biology, Delhi, India, Department of Biotechnology, University of Pune, Pune, India

Abstract

Background

The interactions among various biomarkers remained unexplored under the stressful environment of high-altitude. Present study evaluated interactions among biomarkers to study susceptibility for high altitude pulmonary edema (HAPE) in HAPE-patients (HAPE-p) and adaptation in highland natives (HLs); both in comparison to HAPE-free sojourners (HAPE-f).

Methodology/Principal Findings

All the subjects were recruited at 3500 m. We measured clinical parameters, biochemical levels in plasma and gene expression using RNA from blood; analyzed various correlations between and among the clinical parameters, especially arterial oxygen saturation (SaO2) and mean arterial pressure (MAP) and biochemical parameters like, asymmetric dimethylarginine (ADMA), serotonin (5-HT), 8-iso-prostaglandin F2α (8-isoPGF2α), endothelin-1 (ET-1), plasma renin activity (PRA), plasma aldosterone concentration (PAC), superoxide dismutase (SOD) and nitric oxide (NO) in HAPE-p, HAPE-f and HLs. ADMA, 5-HT, 8-isoPGF2α, ET-1 levels, and PAC were significantly higher (p<0.0001, each), whereas SOD activity and NO level were significantly lower in HAPE-p than HAPE-f (p≤0.001). Furthermore, ADMA, 5-HT, 8-isoPGF2α, NO levels and PAC were significantly higher (p<0.0001), whereas ET-1 level significantly (p<0.0001) and SOD activity non-significantly (p>0.05) lower in HLs than HAPE-f. The expression of respective genes differed in the three groups. In the correlations, SaO2 inversely correlated with ADMA, 5-HT and 8-isoPGF2α and positively with SOD in HAPE-p (p≤0.009). MAP correlated positively with 5-HT and 8-isoPGF2α in HAPE-p and HLs (p≤0.004). A strong positive correlation was observed between ADMA and 5-HT, 5-HT and 8-isoPGF2α (p≤0.001), whereas inverse correlation of SOD with ET-1 in HAPE-p and HLs (p≤0.004), with 5-HT and 8-isoPGF2α in HAPE-p (p = 0.01) and with 5-HT in HLs (p = 0.05).

Conclusions/Significance

The interactions among these markers confer enhanced vascular activity in HLs and HAPE in sojourners.

Introduction

The relevance of hypobaric hypoxia environment at high-altitude does not confine to adaptation but it also provides an excellent set up for studying pathophysiological conditions that result into various disorders [1][3]. Among these disorders high-altitude pulmonary edema (HAPE), a rare but fatal disorder, develops in susceptible sojourners upon rapid ascent to or physical exertion at altitudes above 2500 m [4]. HAPE, a non-cardiogenic pulmonary edema, is characterized by exaggerated pulmonary hypertension, which leads to vascular leakage through overperfusion and stress failure or both [4].

Different physiological pathways are activated under this environment that incorporate cellular changes like pulmonary vasoconstriction response, erythropoiesis and increased blood flow [1][2], [5], however several of the constituents of these pathways have not been evaluated. For example asymmetric dimethylarginine (ADMA), an endogenous nitric oxide synthase inhibitor, is an endothelial dysfunction and oxidative-stress marker [6][8]. It is crucially involved in various cardiovascular diseases [9] but its role in HAPE is yet to be established. Likewise, serotonin also known as 5-hydroxytryptamine, (5-HT), a mitogenic marker [10], [11] and F-iso-prostaglandins, the lipid peroxidation markers [12], [13] have not been explored in HAPE pathophysiology. Beside, the mitogenic and vasoconstrictory biomolecule, endothelin-1 (ET-1) is believed to have a role in HAPE pathogenesis [14][16]. The other constituents of vascular homeostasis cascade such as, plasma renin activity (PRA), angiotensin-1converting enzyme (ACE) and plasma aldosterone concentration (PAC) find their relevance because of their involvement in blood pressure regulation and fluid retention [17][19]. Equally imperative is the function of nitric oxide (NO) and antioxidants such as superoxide dismutase (SOD) that maintain vasodilatation and oxygenation at HA [20][23].

The constituents of vascular homeostasis, thus, become pertinent at HA. Although the basic effects of most of the molecules are known, the normal physiological or pathophysiological changes cannot be credited to an individual molecule; hence, we hypothesized that cross-talk among several molecules of different pathways may lead to HA disorders and adaptation.

To test our hypothesis, we set forth objectives to measure the biochemical parameters like ADMA, 5-HT, 8-isoPGF2α, and also ET-1, PRA, PAC, SOD and NO; to quantify transcriptional activity; importantly, to evaluate the various correlations between and among these molecules and clinical parameters; and, to evaluate association of each parameter with susceptibility to HAPE and their plausible role in adaptation. The study was performed in three well-defined groups, viz. HAPE-patients (HAPE-p), HAPE-free sojourners (HAPE-f) and highland natives (HLs; residing at heights ≥3500 m).

Materials and Methods

Study subjects

Prior to participation in the study, written informed consent was obtained from each participant. The subjects were categorized into three well-defined groups: 1) HAPE-patients (HAPE-p), n = 200, are sojourners who suffered the disorder upon exposure to HA (∼3500 m); 2) HAPE-free sojourners (HAPE-f), n = 200, are the healthy subjects, who visited HA under similar conditions but did not suffer from the disorder, and 3) healthy highland natives (HLs), n = 450, are the permanent residents to HA, residing at and above an altitude of 3500 m. The HLs are of Tibeto-Burman ethnicity residing at HA for many generations, and are declared fit and healthy. All the subjects were recruited through SNM hospital, Leh (∼3500 m), India. The two human ethical committees namely; ‘the human ethical committee of Institute of Genomics and Integrative Biology, Delhi’ and ‘the human ethical committee of Sonam Norboo Memorial Hospital, Leh, Ladakh’ had approved the investigation. The experimental part was performed in our laboratory in Delhi.

Selection criteria

Diagnosis of HAPE was based on chest radiography that revealed infiltrates consistent with pulmonary edema. Pulmonary artery systolic-pressure (PASP) was measured in HAPE-p and HAPE-f using Sonos-5500 echocardiography (Hewlett-Packard, USA). The other clinical symptoms included hypoxemia, cough and dyspnea at rest, presence of pulmonary rales and cyanosis. Of note, HAPE-p acquired the disease upon first exposure to hypobaric hypoxia at 3500 m, whereas HAPE-f visited HA more than once (59%) and none had suffered from HA related disorders. The clinical parameters e.g., body mass index (BMI), blood pressure (BP), SaO2 and pulse rate (PR) were also measured. The subjects were given rest prior to BP measurement. Three measurements of BP, in supine position, using a calibrated mercury sphygmomanometer with appropriate adult cuff size were recorded. The SaO2 and PR were measured using Finger-Pulse Oximeter 503 (Criticare Systems Inc, USA). Any previous history of cardiopulmonary and other diseases was ruled out through questionnaire in the three groups. Lake Louise scoring was applied to rule out any symptoms of acute mountain sickness among HAPE-f.

Biomarker estimation

Blood samples were collected at 3500 m. Ten ml of venous blood was drawn in acid-citrate-dextrose anticoagulant, plasma was collected and stored at −80°C until analyzed. Blood in case of patients was collected once HAPE was diagnosed. The blood from the subjects of the two healthy groups was drawn following overnight fasting.

ELISA was used to measure plasma ADMA (Uscn Life Science Inc, USA), 5-HT, 8-isoPGF2α and ET-1 levels (Assay Designs, USA). PRA and PAC were measured by radioimmunoassay (Immunotech, France) on a gamma counter (Ria Calc Wizard 1470, USA). SOD was estimated by the method of Kono et al. (1978), whereas NO was estimated by an enzymatic assay kit (Cayman Chemical, USA). Barring the radioimmunoassay, all the other measurements were performed on a high-throughput SpectraMax plus384 spectrophotometer (Molecular Devices, USA). The inter- and intra-assay coefficient of variations were less than 5% for ADMA, 5-HT, SOD and NO, whereas they were less than 10% for 8-isoPGF2α, ET-1, PRA and PAC.

Gene expression

Total RNA was extracted from 2 ml whole blood by TRI reagent RT blood (Molecular Research Centre, Inc., USA). RNA quality and quantity was determined using a NanoDrop Technologies (ND-1000) spectrophotometer (Thermo Fisher Scientific, USA) and integrity was checked by running on 1.5% agarose gel. Total RNA, 1.5 µg, was used to generate cDNA by EZ-first strand cDNA synthesis kit (Biological Industries, Israel). The gene expression analysis for tryptophan hydroxylase-1 (TPH-1), gene for 5-HT, endothelin-1 (ET-1), renin (REN), cytochrome P450, family 11, subfamily B, polypeptide 2 (CYP11B2), superoxide dismutase (SOD) and nitric oxide synthase 3 (NOS3) was performed on ten samples each of HAPE-p, HAPE-f and HLs. Primers for real-time PCR were designed using the Pearl Primer software (ver. 1.1.20) and are listed in Table S1. Real-time PCR was performed in duplicate and repeated thrice for each gene and each sample on an ABI Prism 7300 Sequence Detection System (Applied Biosystems, USA) using SYBR Green PCR Master Mix (Applied Biosystems, USA). Relative transcript quantities were calculated using the ΔΔCt method with 18S rRNA as the endogenous reference.

Correlation analyses

Four types of correlations were analyzed between, 1) clinical parameters; SaO2 and MAP, 2) biomarkers; seven biomarkers namely, ADMA, 5-HT, 8-isoPGF2α, ET-1, PRA, PAC and SOD were correlated with each other using regression model, 3) next, the clinical parameters were tested against the seven biomarkers and finally, 4) the comparison was made between these biomarkers with respective gene expression.

Statistical analyses

All the values are presented as mean ± SD. Data were analyzed using a statistical package for social sciences for windows (SPSS ver. 15.0) and EPIINFO ver. 6.0 software. Unpaired Student's t-test with two tailed values was performed to compare the differences between the two groups. Log transformed 5-HT and 8-isoPGF2α levels were used for the analysis. Difference in 5-HT and 8-isoPGF2α levels between the three groups was analyzed by one-way analysis of variance (ANOVA). The regression model was used to witness the overall correlation among biomarkers and clinical parameters. The p value was calculated by analysis of covariance method (ANCOVA) after adjustment with the covariates age, gender and BMI followed by Bonferroni's correction test. A p value of <0.05 was considered statistically significant.

Results

Clinical characteristics

Table 1 portrays the clinical characteristics of HAPE-p versus HAPE-f and HLs versus HAPE-f. The systolic blood pressure (SBP), mean arterial pressure (MAP), respiratory rate (RR), pulse rate (PR), SaO2 and PASP differed between HAPE-p and HAPE-f (p<0.05), and SBP, DBP, MAP and SaO2 between HLs and HAPE-f (p<0.05). Of note, SaO2 at 69±13.5% was the lowest in HAPE-p (p<0.0001); the decrease in SaO2 in patients was almost 20% of the average SaO2 of HAPE-f. Interestingly, SaO2 at 87±5.6% was also lower in HLs as against the 91±3.9% in HAPE-f (p<0.0001). PASP was 1.73 times higher in HAPE-p than HAPE-f (p<0.0001). Table 2 represents the baseline characteristics of HAPE-p versus HAPE-f and HLs versus HAPE-f. The height, weight, and BMI differed between the respective groups (p<0.05).

thumbnail
Table 1. Clinical characteristics of the HAPE-p, HAPE-f and HLs.

https://doi.org/10.1371/journal.pone.0044049.t001

thumbnail
Table 2. Baseline characteristics of the HAPE-p, HAPE-f and HLs.

https://doi.org/10.1371/journal.pone.0044049.t002

Biomarker status

Table 3 represents the status of biomarkers in HAPE-p against HAPE-f and HLs against HAPE-f. Significantly higher levels of ADMA, 5-HT, 8-isoPGF2α, ET-1, PAC (p<0.0001, each) and PRA (p = 0.02); whereas, significantly lower SOD activity and NO level were obtained in HAPE-p than HAPE-f (p≤0.001). In terms of fold difference, the levels of ADMA, 5-HT, 8-isoPGF2α, ET-1, PRA and PAC were higher by 4.19, 1.5, 1.5, 2.3, 2.3 and 2.0 times, while SOD activity and NO level were lower by 1.7 and 1.27 times, respectively, in HAPE-p.

thumbnail
Table 3. Plasma levels of biomarkers in HAPE-p and HAPE–f and HLs.

https://doi.org/10.1371/journal.pone.0044049.t003

As can be seen from table 3, significantly higher levels of ADMA, 5-HT, 8-isoPGF2α, PAC, NO (p<0.0001, each) and PRA (p = 0.04) were obtained in HLs than HAPE-f; contrary to this, significantly lower ET-1 level and lower SOD activity were obtained in HLs (p<0.0001 and p = 0.40, respectively). In terms of fold difference, the levels of ADMA, 5-HT, 8-isoPGF2α, PRA, PAC and NO were higher by 1.72, 1.2, 1.3, 1.4, 2.0 and 1.65 times, respectively, when compared with HAPE-f, whereas ET-1 level and SOD activity were lower by 2.1 and 1.2 times, respectively, in HLs than HAPE-f.

Gene expression

Figure 1 represents the expression pattern of TPH-1, ET-1, REN, CYP11B2, SOD and NOS3. In HAPE-p, the expression of TPH-1, ET-1, CYP11B2 and REN was up-regulated by 2.0, 2.45, 3.06 and 1.84 fold (p = 0.01, p = 0.015, p = 0.032 and p>0.05, respectively), whereas the expression of SOD and NOS3 was down-regulated by 2.46 and 1.55 fold, respectively, against HAPE-f (p = 0.01, p>0.05). In HLs, the expression of TPH-1, CYP11B2, NOS3 and REN was up-regulated by 3.0, 2.76, 5.2 and 2.05 fold (p = 0.001, p = 0.023, p = 0.001 and p>0.05), against HAPE-f, respectively, whereas ET-1 and SOD was down-regulated by 1.24 and 1.07 fold (p>0.05), respectively, against HAPE-f.

thumbnail
Figure 1. Relative gene expression by real-time PCR expressed as fold change in HAPE-p and HLs.

Bars (mean ± SD) above the base line show upregulation and bars below the base line show downregulation. The differential expression (+/−) varied between 1.01–5.2. TPH-1, tryptophan hydroxylase-1; ET-1, endothelin-1; REN, renin; CYP11B2, cytochrome P450, family 11, subfamily B, polypeptide 2; SOD, superoxide dismutase; NOS3, nitric oxide synthase 3; HAPE-p, high altitude pulmonary edema patients; HLs, healthy highland natives; HAPE-f, High altitude pulmonary edema free sojourners. The comparisons of HAPE-p and HLs were made against HAPE-f.

https://doi.org/10.1371/journal.pone.0044049.g001

Correlations, the extent of interactions

Figure 2 and Figure S1, S2, S3, S4 represent various correlation analyses in the three studied groups.

thumbnail
Figure 2. Clinical correlations in the three groups i.e. HAPE-p, HAPE-f and HLs.

An inverse correlation was obtained between the clinical parameters viz MAP and SaO2 in the three groups; MAP, mmHg; SaO2, %.

https://doi.org/10.1371/journal.pone.0044049.g002

Correlations between clinical parameters

A strong inverse correlation was observed between SaO2 and MAP in HAPE-p (p = 0.0006), HAPE-f (p≤0.0001) and HLs (p≤0.0001) as can be seen from figure 2.

Correlations between Biomarkers

Here, Figure S1, represents nine correlation comparisons between the biomarkers except NO, which has been presented in Figure S2. A strong positive correlation was observed between ADMA and 5-HT in HAPE-p and HAPE-f (Figure S1 a), ADMA also showed an inverse correlation with SOD in HAPE-p and HAPE-f (p≤0.0001; Figure S1b). ADMA also showed a marginal correlation with 8-isoPGF2α only in HLs (p = 0.04; Figure S1c). 5-HT showed a strong positive correlation with 8-isoPGF2α in the three groups (p≤0.001; Figure S1d) and with ET-1 in HLs (p = 0.03, Figure S1e); while, it inversely correlated with SOD in HLs (p = 0.05; Figure S1f). Although, 8-isoPGF2α was positively correlated with 5-HT (p≤0.001; Figure S1d), it did not correlate with ET-1 in the three groups (p>0.05; Figure S1g), whereas, it inversely correlated with SOD in HAPE-p and HAPE-f (p = 0.01; Figure S1h). ET-1 showed significant inverse correlation with SOD in HAPE-p and HLs (p≤0.004; Figure S1i). In addition, NO was inversely correlated with ADMA, 5-HT and 8-isoPGF2α in the three groups (p≤0.031, p≤0.06, and p≤0.022, respectively, Figure S2). The correlation outcome on renin and aldosterone was not encouraging hence it is not included here.

Correlations between clinical parameters and Biomarkers

Here, seven correlation comparisons are provided, between clinical parameters and biomarkers in Figure S3, and between clinical parameters and NO in Figure S4. SaO2 showed inverse correlation with ADMA in the three groups (p≤0.0084; Figure S3a), with 5-HT and 8-isoPGF2α in HAPE-p (p = 0.0001; Figure S3b; p = 0.0003; Figure S3c, respectively); whereas, it was positively correlated with SOD in HAPE-p (p = 0.009; Figure S3d), and with NO in the three groups (p≤0.022; Figure S4a). MAP showed significant positive correlation with ADMA in HAPE-p (p<0.05; Figure S3e), with 5-HT in the three studied groups (p = 0.0001; Figure S3f), with 8-isoPGF2α in HAPE-p and HLs (p≤0.004; Figure S3g) and inversely with NO in the three groups (p≤0.02; Figure S4b).

Correlation between Biomarkers and Gene expression

The correlation comparison between biomarkers and the respective gene expression profile was evaluated. Plasma 5-HT concentration and TPH-1 gene expression were upregulated by 1.5 and 2.0 fold, respectively, in HAPE-p, and 1.2 and 3.0 fold, respectively, in HLs. Plasma ET-1 concentration and ET-1 gene expression in HAPE-p were upregulated by 2.3 and 2.45 fold, respectively, whereas in HLs both were downregulated by 2.1 and 1.24 fold, respectively. Plasma renin and aldosterone concentration and gene expression in HAPE-p were upregulated by 2.3, 1.84 and 2.0, 3.06 fold, respectively; also in HLs, the two were upregulated by 1.4, 2.05 and 2.0, 2.76 fold, respectively. Contrary to these findings, plasma SOD activity and SOD gene expression were downregulated by 1.7 and 2.46 fold, respectively, in HAPE-p; it followed with 1.2 and 1.07 fold downregulation, respectively, in HLs. A decrease of 1.27 and 1.55 fold in plasma NO and NOS3 gene expression was observed in HAPE-p, whereas the fold changes were elevated by 1.65 and 5.2 fold in plasma and gene expression, respectively, in HLs. This comparison compliments the biolevels and its respective gene expression.

The multiple correlations

The multiple regression analysis provided the depiction of interaction among these biomarkers in the etiology of HAPE and adaptation (Figure S5). Positive correlations were observed among the vasoconstrictors like ADMA, 5-HT, 8-isoPGF2α and ET-1, (p between 0.03–0.000) and inverse correlations of these vasoconstrictors were observed with SOD and NO (p value ranges from 0.066–0.000) in the three groups. Among the clinical parameters, MAP showed positive correlation with ADMA, 5-HT, 8-isoPGF2α and ET-1 (p>0.05) and inverse correlation with SOD and NO (p value ranges from 0.37–0.018). SaO2 correlated inversely with ADMA, 5-HT, 8-isoPGF2α and ET-1 (p value ranges from 0.05–0.001) and positively with SOD and NO (p value ranges from 0.4–0.02). The schematic interactions of the studied molecules are profiled in Figure 3. The interactions between these molecules clearly suggested impairment of the endothelial function and vascular homeostasis (Figure 3) thereby leading to HAPE.

thumbnail
Figure 3. Schematic presentation of the studied biomarkers in a physiological function under hypobaric hypoxia.

Interactions of several of these biomarkers that translate into the incidence of HAPE are also interpretable through the renin–angiotensin–aldosterone system, Kinin–kallikrein system, and the pathways of ET-1, 5-HT, NO signaling and oxidative-stress. ADMA, asymmetric dimethylarginine; ATII, angiotensin II; AT-1R, angiotensin-II type I receptor; AT-2R, angiotensin-II type II receptor; AngI, angiotensin I; AngII, angiotensin II; ACE, angiotensin-I converting enzyme; AGT, angiotensinogen; 8-isoPGF2α, 8-iso-prostaglandinF2α; ROS, reactive oxygen species; SOD, superoxide dismutase; O2˙-, superoxide anion; ONOO, peroxynitrite; 5-HT, serotonin; ET-1, endothelin-1; ET-A and ET-B, endothelin receptors A and B; NO, nitric oxide; NOS3, endothelial nitric oxide synthase; VSMC, vascular smooth muscle cell; H2O2, hydrogen peroxide; H2O, water; EFA, essential fatty acids.

https://doi.org/10.1371/journal.pone.0044049.g003

Discussion

The highlight of our study is the investigation of several biomarkers with respect to HAPE and adaptation; few biomarkers were investigated until now.

It was encouraging to note that these molecules provided a glimpse of individual and interactive-mode function and the phenotypic outcome. ADMA, one of the endothelial dysfunction and oxidative-stress markers was significantly elevated in HAPE. It is synthesized during the methylation of arginine residues in a protein, and is released during proteolysis. ADMA plays three crucial roles; it competitively inhibits the cellular L-arginine uptake [24], [25], thereby inhibiting endogenous NOS3, enhances ACE activity [26] and elevates reactive oxygen species [8], [26]. We figured out that all the three functions were involved in HAPE, because of the fact that NO levels were found consistently lower, whereas ACE activity and aldosterone levels and the oxidative stress markers were higher in the present and previous studies [14], [27]. It is known that in renin-angiotensin aldosterone system, angiotensin II produced by ACE stimulates aldosterone synthase to secrete aldosterone that regulates water and salt retention [28]; on the other hand, in kallikren-kinin system, ACE proteolyses bradykinin deterring NO production, the net result of this cascade being increased vascular resistance and high BP [29], [30]. Likewise, our observation of lower activity of SOD in the present study and higher level of 8-isoPGF2α [31] were reminiscent of dysfunction of the oxidative-stress pathway and its causative role in HAPE. In fact, the expression level of two related genes was in accordance to these markers such as higher expression of pro-oxidant gene, cytochrome b-245, alpha polypeptide and lower expression of another antioxidant gene, glutathione S-transferase pi 1 [31]. The formation of abnormal phospholipids [32] such as in the present situation apparently exerted profound effect on membrane fluidity and integrity, a sign of oxidant injury. The positive correlation of ADMA with 5-HT and ET-1 pointed to exaggerated vasoconstriction.

Though, these effects were more pronounced in HAPE, similar trends were noticed in natives, pointing to the fact that they were also under continuous stress but were more or less adapted to such a stressful environmental condition. It, importantly, points to the fact that chronic exposure results into elevated vascular reactivity but not as severely as in acute exposure in susceptible sojourners.

Another important finding was the elevated levels of 5-HT that underline its causative role in the pathogenesis of HAPE. 5-HT exerts mitogenic and co-mitogenic effects on pulmonary artery smooth muscle cells [10], [33]. The constricting action of 5-HT on smooth muscle cells is mainly mediated through G-protein coupled receptors, 5-HT1B and 5-HT2A, and internalization by its transporter, which eventually potentiates PH and associated diseases [10], [34][35]. 5-HT is also involved in ROS generation through the interaction with NADPH oxidase, and proliferation through the interaction with angiotensin II [36], [37]. Although investigations have delineated the potential role of ADMA and 5-HT in the pathophysiology of PH [9][11], [38], [39], knowledge with respect to HAPE was acquired only through present study. Among the other physiological systems, ET-1 is a known contributor to the pathogenesis of HAPE [14][16]. Berger and colleagues (2009) had reported its involvement in PASP through hypoxia-induced pulmonary vasoconstriction [40], which is concurrent to our findings of elevated ET-1 and PASP in HAPE.

Equally pertinent would be to describe here the role of NO, which being a critical mediator of vasodilation remains cardinal to these physiological inferences. However, under the stressful environment of HA, several of the mediators of vasoconstriction interfere with the production of NO. True to this statement, we could observe inverse interaction of NO with vasoconstrictors; the net result was upsetting the normal physiological function.

It would further be pertinent to highlight the correlations between SaO2 and stress markers; the former is an index of oxygen affinity in blood at HA. An inverse correlation of SaO2 with MAP, ADMA, 5-HT and 8-isoPGF2α was an indication that the homeostasis of body had shifted towards disease susceptibility, which was observed in the patients. Higher oxygen availability in blood and controlled production of stress markers may prove fruitful in defining physical activity and even high endurance. This study, thus, imparts imperative knowledge about the SaO2 phenotype in relation to health and adaptation.

Finally, our evaluation of biomarkers in native population has been equally motivating and interesting. Although, Tibetans are one of the best-adapted populations, however, until present study literature was scant on the role of these biomarkers. It is of consequence to add that HLs revealed higher ADMA, 5-HT, 8-isoPGF2α, PRA and PAC and among the clinical parameters, relatively higher DBP, SBP and MAP. Yet, by virtue of their long term residence, any kind of negative consequences may have been counterbalanced by factors that maintained vascular tone at HA, such as reduced plasma ET-1 and higher NO and SaO2 levels as was observed in the present study, and as was reported earlier [41][43]. It is likely that the HLs have evolved few assorted mechanisms in response to hypoxia. Thus, a normal or enhanced vagal tone and preserved vasomotion in HLs are probably evidence of adaptation.

It may be of relevance to add on genetic aspects in relation to biomarkers in HAPE and adaptation. Last two decades focused on the candidate markers of vascular homeostasis pathway and oxidative stress. The notable genetic variants from among the several genes studied were those of angiotensin-I converting enzyme (ACE); angiotensinogen; aldosterone synthase; nitric oxide synthase 3 (NOS3); endothelin-1 (ET-1); β-adrenergic receptor type 2; cytochrome b-245 alpha polypeptide (CYBA) and glutathione S-transferase pi 1 (GSTP1). Several of the variants were reported to associate with adaptation and disorder. Interestingly, the variants were associated with the biolevels. Such as the ACE D allele of I/D polymorphism was associated with higher ACE activity and lower SaO2 level in patients or the I allele with lower activity and higher SaO2 level in healthy subjects [14], [43]; ET-1 longer repeats of (CT)n-(CA)n and NOS3 4a allele of 4b/4a with lower levels of ET-1 and NO, respectively [14], [42]. CYBA A and T alleles of −930A/G and H72Y (C/T), GSTP1 A and C alleles of I105V (A/G) and A114V (C/T) were associated with lower levels of 8-iso PGF2α [31]. A reverse trend was observed with the other alleles on the respective same loci. In recent years the oxygen sensing pathway has been reported to associate with adaptation in highlanders; the genome-wide association studies have revealed the involvement of the genes EGL nine homolog 1 and endothelial PAS domain protein 1 [44][46] in adaptation.

The prominent findings of the present study are the emergence of ADMA, 5-HT and 8-isoPGF2α as novel markers having promising role in the regulation of physiological functions at HA; as a result increased oxidative and pulmonary stress could be predicted. Furthermore it was also evident from our results that these biomarkers bonded together to produce the phenotypic outcome. The relative involvement, however, has to be evaluated through continued efforts. The added advantage of this study was that it was carried out in a fairly larger sample size. Among the limitations, PASP could not be measured in highlanders, as these being healthy native individuals did not consent to additional tests unless significantly required. PASP was measured in HAPE-p and HAPE-f using echocardiography and not through right heart catheterization, as the latter is mandatory for severe heart and pulmonary ailments like IPAH, FPAH etc and is a risky procedure thus is not permitted in SNM hospital. The two baseline characteristics, age and BMI need mention here. The HLs were older in age than the two other groups; however this age difference apparently did not influence the levels, although a younger age group would have validated this. The BMI in this context was in support of the above notion. Despite the age difference, the BMI was within the normal range as per national guidelines. Moreover, the two controls HAPE-p and HLs although revealed statistical difference but the mean BMI value was 19 and 20 kg/m2. Also BMI is not a phenotype of importance in HAPE hence in the present study BMI and age were used as confounders. Hence it seems these two parameters did not influence the biolevels. Nevertheless these parameters would be scrutinized in further studies.

Supporting Information

Figure S1.

(a–i): Correlation analyses among biochemical parameters in the three groups i.e. HAPE-p, HAPE-f and HLs. Proportional or inverse correlations were obtained between a) ADMA, nmol/mL and 5-HT, ng/mL; b) ADMA, nmol/mL and SOD, U/mL; c) ADMA, nmol/mL and 8-isoPGF2α, pg/mL; d) SOD, U/mL and ET-1, pg/mL; e) 5-HT, ng/mL and ET-1, pg/mL; f) and 8-isoPGF2α, pg/mL and ET-1, pg/mL; g) 8-isoPGF2α, pg/mL and 5-HT, ng/mL; h) SOD, U/mL and 5-HT, ng/mL; i) SOD, U/mL and 8-isoPGF2α, pg/mL.

https://doi.org/10.1371/journal.pone.0044049.s001

(TIF)

Figure S2.

(a–e): Correlation analyses among biochemical parameters with NO in the three groups i.e. HAPE-p, HAPE-f and HLs. Proportional or inverse correlations were obtained between a) NO, µmol/mL and ADMA, nmol/mL; b) NO, µmol/mL and 5-HT, ng/mL; c) NO, µmol/mL and 8-isoPGF2α, pg/mL; d) NO, µmol/mL and SOD, U/mL; e) NO, µmol/mL and ET-1, pg/mL.

https://doi.org/10.1371/journal.pone.0044049.s002

(TIF)

Figure S3.

(a–g): Correlation analyses between biochemical parameters and clinical parameters in the three groups i.e. HAPE-p, HAPE-f and HLs. Proportional or inverse correlations were obtained between a) ADMA, nmol/mL and SaO2, %; b) 5-HT, ng/mL and SaO2,%; c) 8-isoPGF2α, pg/mL and SaO2, %; d) SOD, U/mL and SaO2, %; e) ADMA, nmol/mL and MAP, mmHg; f) 8-isoPGF2α, pg/mL and MAP, mmHg; g) 5-HT, ng/mL and MAP, mmHg.

https://doi.org/10.1371/journal.pone.0044049.s003

(TIF)

Figure S4.

(a–b): Correlation analyses between NO and clinical parameters in the three groups i.e. HAPE-p, HAPE-f and HLs. Proportional or inverse correlations were obtained between a) NO, µmol/mL and SaO2, %; b) NO, µmol/mL and MAP, mmHg.

https://doi.org/10.1371/journal.pone.0044049.s004

(TIF)

Figure S5.

(a–c): Multiple regression analysis in the three groups i.e. HAPE-p, HAPE-f and HLs. The boxes with shades of red depicted positive correlations and boxes with shades of green depicted inverse correlations. Each box contains the coefficient of correlation along with p value. The p value was calculated placing one biomarker against all the studied biomarkers.

https://doi.org/10.1371/journal.pone.0044049.s005

(TIF)

Table S1.

Real-Time PCR conditions for TPH-1 , ET-1 , REN , CYP11B2 , SOD and NOS3 .

https://doi.org/10.1371/journal.pone.0044049.s006

(DOC)

Acknowledgments

We acknowledge the support and encouragement of director, IGIB. The technical support of the staff at Sonam Norboo Memorial Hospital, Leh and all the volunteers is appreciated.

Author Contributions

Conceived and designed the experiments: MAQP. Performed the experiments: ZA AM PP RR TT GM. Analyzed the data: RK PA. Contributed reagents/materials/analysis tools: ZA RK PA. Wrote the paper: ZA MAQP.

References

  1. 1. Hackett PH, Roach RC (2001) High-Altitude Illness. N Engl J Med 345: 107–114.
  2. 2. Beall CM (2007) Detecting natural selection in high-altitude human populations. Respir Physiol Neurobiol 158: 161–171.
  3. 3. Moore LG, Niermeyer S, Zamudio S (1998) Human adaptation to high altitude: regional and life-cycle perspectives. Am J Phys Anthropol 41: 25–64.
  4. 4. Basnyat B, Murdoch DR (2003) High-altitude illness. Lancet 361: 1967–1974.
  5. 5. León-Velarde F, Villafuerte FC, Richalet JP (2010) Chronic mountain sickness and the heart. Prog Cardiovasc Dis 52: 540–549.
  6. 6. Skoro-Sajer N, Mittermayer F, Panzenboeck A, Bonderman D, Sadushi R, et al. (2007) Asymmetric dimethylarginine is increased in chronic thromboembolic pulmonary hypertension. Am J Respir Crit Care Med 176: 1154–1160.
  7. 7. Zakrzewicz D, Eickelberg O (2009) From arginine methylation to ADMA: a novel mechanism with therapeutic potential in chronic lung diseases. B M C Pulm Med 9: 5.
  8. 8. Arrigoni FI, Vallance P, Haworth SG Leiper JM (2003) Metabolism of asymmetric dimethylarginines is regulated in the lung developmentally and with pulmonary hypertension induced by hypobaric hypoxia. Circulation 107: 1195–1201.
  9. 9. Vallance P, Leiper J (2004) Cardiovascular biology of the asymmetric dimethylarginine:dimethylarginine dimethylaminohydrolase pathway. Arterioscler Thromb Vasc Biol 24: 1023–1030.
  10. 10. Eddahibi S, Fabre V, Boni C, Martres MP, Raffestin B, et al. (1999) Induction of serotonin transporter by hypoxia in pulmonary vascular smooth muscle cells. Relationship with the mitogenic action of serotonin. Circ Res 84: 329–336.
  11. 11. Herve P, Launay JM, Scrobohaci ML, Brenot F, Simonneau G, et al. (1995) Increased plasma serotonin in primary pulmonary hypertension. Am J Med 99: 249–254.
  12. 12. Montuschi P, Barnes PJ, Roberts II LJ (2004) Isoprostanes: markers and mediators of oxidative stress. FASEB J 18: 1791–1800.
  13. 13. Morrow JD, Roberts LJ (2002) The isoprostanes: their role as an index of oxidative stress status in human pulmonary disease. Am J Respir Crit Care Med 166: S25–S30.
  14. 14. Charu R, Stobdan T, Ram RB, Khan AP, Qadar Pasha MA (2006) Susceptibility to high-altitude pulmonary oedema: role of ACE and ET-1 polymorphisms. Thorax 61: 1011–1012.
  15. 15. Sartori C, Vollenweider L, Löffler BM, Delabays A, Nicod P, et al. (1999) Exaggerated endothelin release in high-altitude pulmonary edema. Circulation 99: 2665–2668.
  16. 16. Droma Y, Hayano T, Takabayashi Y, Koizumi T, Kubo K, et al. (1996) Endothelin-1 and interleukin-8 in high altitude pulmonary oedema. Eur Respir J 9: 1947–1949.
  17. 17. MacGregor GA, Markandu ND, Roulston JE, Jones JC, Morton JJ (1981) Maintenance of blood pressure by the renin-angiotensin system in normal man. Nature 291: 329–331.
  18. 18. Rock PB, Kraemer WJ, Fulco CS, Trad LA, Malconian MK, et al. (1993) Effects of altitude acclimatization on fluid regulatory hormone response to submaximal exercise. J Appl Physiol 75: 1208–1215.
  19. 19. Milledge JS, Catley DM (1982) Renin, aldosterone, and converting enzyme during exercise and acute hypoxia in humans. J Appl Physiol 52: 320–323.
  20. 20. Busch T, Bärtsch P, Pappert D, Grünig E, Elser H, et al. (2001) Hypoxia decreases exhaled nitric oxide in mountaineers susceptible to high altitude pulmonary edema. Am J Respir Crit Care Med 163: 368–373.
  21. 21. Beall CM, Laskowski D, Strohl KP, Soria R, Villena M, et al. (2001) Pulmonary nitric oxide in mountain dwellers. Nature 414: 411–412.
  22. 22. Nakanishi K, Tajima F, Nakamura A, Yagura S, Ookawara T, et al. (1995) Antioxidant system in hypobaric-hypoxia. J Physiol 489: 869–876.
  23. 23. Dosek A, Ohno H, Acs Z, Taylor AW, Radak Z (2007) High altitude and oxidative stress. Respir Physiol Neurobiol 158: 128–131.
  24. 24. Moncada S, Higgs A (1993) The L-arginine-nitric oxide pathway. N Engl J Med 329: 2002–2012.
  25. 25. Cooke JP, Dzau VJ (1997) Derangements of the nitric oxide synthase pathway, L-arginine, and cardiovascular diseases. Circulation 96: 379–382.
  26. 26. Suda O, Tsutsui M, Morishita T, Tasaki H, Ueno S, et al. (2004) Asymmetric dimethylarginine produces vascular lesions in endothelial nitric oxide synthase-deficient mice: involvement of renin-angiotensin system and oxidative stress. Arterioscler Thromb Vasc Biol 24: 1682–1688.
  27. 27. Qadar Pasha MA, Khan AP, Kumar R, Grover SK, Ram RB, et al. (2001) Angiotensin converting enzyme insertion allele in relation to high altitude adaptation. Ann Hum Genet 65: 531–536.
  28. 28. Struthers AD, MacDonald TM (2004) Review of aldosterone-and angiotensin II-induced target organ damage and prevention. Cardiovasc Res 61: 663–670.
  29. 29. Lemarié CA, Paradis P, Schiffrin EL (2008) New insights on signaling cascades induced by cross-talk between angiotensin II and aldosterone. J Mol Med 86: 673–678.
  30. 30. Linz W, Wohlfart P, Schölkens BA, Malinski T, Wiemer G (1999) Interactions among ACE, kinins and NO. Cardiovasc Res. 43: 549–561.
  31. 31. Mishra A, Ali Z, Vibhuti A, Kumar R, Alam P, et al. (2011) CYBA and GSTP1 variants associate with oxidative stress under hypobaric hypoxia as observed in high-altitude pulmonary edema. Clin Sci 122: 299–309.
  32. 32. Jefferson JA, Simoni J, Escudero E, Hurtado ME, Swenson ER, et al. (2004) Increased oxidative stress following acute and chronic high altitude exposure. High Alt Med Biol 5: 61–69.
  33. 33. Marcos E, Fadel E, Sanchez O, Humbert M, Dartevelle P, et al. (2004) Serotonin-induced smooth muscle hyperplasia in various forms of human pulmonary hypertension. Circ Res 94: 1263–1270.
  34. 34. Königshoff M, Dumitrascu R, Udalov S, Amarie OV, Reiter R, et al. (2010) Increased expression of 5-hydroxytryptamine2A/B receptors in idiopathic pulmonary fibrosis: a rationale for therapeutic intervention. Thorax 65: 949–955.
  35. 35. MacLean MR (2007) Pulmonary hypertension and the serotonin hypothesis: where are we now? Int J Clin Pract Suppl 61: 27–31; 2007.
  36. 36. Lee SL, Wang WW, Fanburg BL (1998) Superoxide as an intermediate signal for serotonin-induced mitogenesis. Free Radic Biol Med 24: 855–858.
  37. 37. Morrell NW, Adnot S, Archer SL, Dupuis J, Jones PL, et al. (2009) Cellular and molecular basis of pulmonary arterial hypertension. J Am Coll Cardiol 54: S20–31.
  38. 38. Gorenflo M, Zheng C, Werle E, Fiehn W, Ulmer HE (2001) Plasma levels of asymmetrical dimethyl-L-arginine in patients with congenital heart disease and pulmonary hypertension. J Cardiovasc Pharmacol 37: 489–492.
  39. 39. Kielstein JT, Bode-Böger SM, Hesse G, Martens-Lobenhoffer J, Takacs A, et al. (2005) Asymmetrical dimethylarginine in idiopathic pulmonary arterial hypertension. Arterioscler Thromb Vasc Biol 25: 1414–1418.
  40. 40. Berger MM, Dehnert C, Bailey DM, Luks AM, Menold E, et al. (2009) Transpulmonary plasma ET-1 and nitrite differences in high altitude pulmonary hypertension. High Alt Med Biol 10: 17–24.
  41. 41. Rajput C, Najib S, Norboo T, Afrin F, Qadar Pasha MA (2006) Endothelin-1 gene variants and levels associate with adaptation to hypobaric hypoxia in high-altitude natives. Biochem Biophys Res Commun 341: 1218–1224.
  42. 42. Ahsan A, Charu R, Norboo T, Afrin F, Baig MA, et al. (2004) eNOS allelic variants at the same locus associate with HAPE and adaptation. Thorax 59: 1000–1002.
  43. 43. Qadar Pasha MA, Kocherlakota KS, Khan AP, Norboo T, Grover SK, et al. (2003) Arterial oxygen saturation under hypoxic environment of high-altitude associates with routine physical activities of natives. Current Science 85: 502–506.
  44. 44. Aggarwal S, Negi S, Jha P, Singh PK, Stobdan T, et al. (2010) EGLN1 involvement in high-altitude adaptation revealed through genetic analysis of extreme constitution types defined in Ayurveda. Proc Natl Acad Sci U S A 107: 18961–18966.
  45. 45. Simonson TS, Yang Y, Huff CD, Yun H, Qin G, et al. (2010) Genetic evidence for high-altitude adaptation in Tibet. Science 329: 72–75.
  46. 46. Beall CM, Cavalleri GL, Deng L, Elston RC, Gao Y, et al. (2010) Natural selection on EPAS1 (HIF2alpha) associated with low hemoglobin concentration in Tibetan highlanders. Proc Natl Acad Sci U S A 22; 107: 11459–11464.