Introduction
The leading cause of mortality in women is cardiovascular disease, accounting for a quarter of deaths in developed as well as low- to middle-income societies. Reference Vogel, Acevedo and Appelman1 It is known that women with a history of pregnancy complications, such as preeclampsia, have more than double the risk of future cardiovascular disease compared to women with healthy pregnancies. Reference Ray, Vermeulen, Schull and Redelmeier2,Reference Çağlar, Ozde and Bostancı3 Therefore, diagnosis of complicated pregnancy offers a window of opportunity to identify and treat women at risk of future cardiovascular complications, with the hope of ameliorating the female mortality rate.
It is also established that pregnancy at high altitude is an independent risk factor for preeclampsia and increased maternal mortality. Reference Keyes, Armaza, Niermeyer, Vargas, Young and Moore4,Reference Grant, Giussani and Aiken5 This is a particular problem for Bolivia, where a quarter of its 11.6 million population lives in the cities of La Paz and El Alto above 3400 m or 11,000 ft. Together with other populations residing at altitude, this makes Bolivian women part of the largest single human group at risk for preeclampsia during pregnancy with an increased future risk of cardiovascular mortality. Preeclampsia in Bolivia, with an incidence as high as 14% in the city of La Paz, is the second greatest cause of maternal death. Reference Gutiérrez, Ochoa and Castillo6,Reference De la Gálvez7 It is estimated that 38% of maternal mortality in Bolivia stems from maternal deaths in La Paz and El Alto. Reference Gutiérrez, Ochoa and Castillo6,Reference De la Gálvez7 This comprises a maternal mortality rate in Bolivia of 346/100,000 live births, compared to just 7/100,000 in the UK. Reference Kassebaum, Bertozzi-Villa and Coggeshall8 Unsurprisingly rates for maternal mortality in Bolivia are the second highest in the Western hemisphere. Reference Kassebaum, Bertozzi-Villa and Coggeshall8
While a link between preeclampsia and an increased risk of pulmonary arterial hypertension is established in sea-level pregnancies not only during pregnancy but also long after birth, Reference Ray, Vermeulen, Schull and Redelmeier2,Reference Çağlar, Ozde and Bostancı3 whether women who develop preeclampsia during highland pregnancy are also more likely to develop pulmonary arterial hypertension before and after pregnancy has not been studied. This is important as both preeclampsia and pulmonary hypertension independently contribute to significant maternal mortality, Reference Ray, Vermeulen, Schull and Redelmeier2,Reference Martin and Edwards9 and the risk of both conditions are elevated at high altitude. Reference Keyes, Armaza, Niermeyer, Vargas, Young and Moore4,Reference Grant, Giussani and Aiken5 Therefore, highland women may be at much greater risk of preeclampsia and pulmonary hypertension, and both conditions may contribute to their much greater long-term cardiovascular risk and mortality. To address this, we carried out a prospective cohort study in Bolivia to test the hypothesis that during highland pregnancy, women who develop preeclampsia are at greater risk of pulmonary arterial hypertension compared to women undergoing healthy pregnancy at high altitude. The hypothesis was tested by determining, in the third-trimester, indices of right heart and pulmonary vascular function via echocardiography in Bolivian women undergoing healthy pregnancy and pregnancy diagnosed with preeclampsia. To determine if any differences persisted, assessments were repeated ca. 6 weeks after birth.
Methods
Data, materials and code disclosure statement
The data that support the findings of this study are available herein. Additional relevant data can be requested from the corresponding author.
Study cohort
Using a multi-institutional approach, women undergoing healthy pregnancies or pregnancies complicated by preeclampsia were recruited following their written consent from 5 hospitals with maternity units in La Paz. The institutions involved were the Hospital Materno-Infantil (Caja Nacional de la Salud), the Hospital de la Mujer, Hospital Nuestra Señora de La Paz, the Hospital Cotahuma, and the Centro de Salud Materno Infantil Tembladerani. Preeclampsia was diagnosed by obstetricians according to guidelines set by the Royal College of Obstetrics & Gynaecology in the UK and the American College of Obstetrics and Gynecology (Hypertension in Pregnancy Task Force, 2013). 10 These guidelines diagnose preeclampsia as new-onset hypertension (defined as a systolic blood pressure ≥ 140 mmHg and/or a diastolic blood pressure ≥ 90 mmHg), most often after 20 weeks of gestation and frequently near term. Abnormal blood pressure readings were usually accompanied by other symptoms including proteinuria, a low platelet count, impaired liver function, signs of kidney problems other than protein in the urine, pulmonary oedema and new-onset headaches or visual disturbances.
The inclusion criteria were that women should be highland residents for at least two generations, non-smoking mothers of normal BMI and no history of health complications carrying a singleton pregnancy and a maternal age between 18–45 years. The study recruited 208 pregnant women who satisfied the inclusion criteria at the time of enrolment. This initial cohort was divided into 110 control (C) and 98 preeclamptic (PE) women. Several women were lost during follow-up for reasons including change of mind, incomplete assessments, and infant death, to yield a final cohort of C (n = 70) and PE (n = 70) women.
Cardiovascular assessment
Maternal B-mode, M-mode, spectral pulsed-wave (PW) Doppler and spectral tissue Doppler imaging (PW-TDI) were applied at 35 (4.3) weeks of pregnancy and 6 (2.2) weeks after birth (median, IQR). One investigator (C.E.S.) performed all pre- and post-natal maternal ultrasound examinations using the same machine (VIVID, Model 2007, Digital Graphic Printer Up- D897). M-mode ultrasound was used for the assessment of left ventricular geometry and for right ventricular (RV) longitudinal contractility by measurement of the tricuspid annular plane systolic excursion (TAPSE). Reference Lang, Badano and Mor-Avi11 B-mode imaging was performed for measurements of pulmonary artery size, RV and right atrial (RA) dimensions in systole and diastole. RV relative wall thickness was estimated as twice the RV free wall thickness divided by RV ventricular end-diastolic diameter. The PW Doppler technique was used to obtain Doppler signals from the RV outflow tract (RVOT) with measurements of peak velocity (Vmax), peak gradient and the velocity-time integral (VTI). The tricuspid regurgitation velocity (TR Vmax) and the RV-RA peak gradient were obtained for the calculation of the pulmonary artery pressure (PAP) by application of the modified Bernoulli equation 4x[TR Vmax]2+expected RA pressure. Reference Rudski, Lai and Afilalo12 Assessment of changes in vena caval collapse did not show any differences between groups. Therefore, a fixed RA pressure of 10 mmHg was used for the estimation of PAP in all women. Pulmonary vascular resistance was then estimated as 10 × (TR Vmax/VTIRVOT) + 0.16. Reference Abbas, Fortuin, Schiller, Appleton, Moreno and Lester13,Reference Augustine, Coates-Bradshaw and Willis14 The PW-TDI technique was applied to derive RV myocardial time-intervals: isovolumetric contraction time (IVCT’), isovolumetric relaxation time (IVRT’) and ejection time (FT’). These time-interval indices were then used to estimate the RV myocardial performance index (MPI’), calculated as MPI’=[(IVCT’+IVRT’)]/ET’. Reference Lang, Badano and Mor-Avi11 All echocardiographic measurements were performed following the American Society of Echocardiography guidelines for echocardiographic assessment of the right heart in adults and recommendations from the British Society of Echocardiography for the assessment of pulmonary hypertension. Reference Lang, Badano and Mor-Avi11,Reference Rudski, Lai and Afilalo12,Reference Augustine, Coates-Bradshaw and Willis14
Maternal oxygen saturation and spirometry
At the same gestational age and post-partum period than echocardiography was performed, maternal oxygen peripheral saturation was measured by pulse oximetry (Criticare System INC, model 504-504P). Maternal respiratory function was determined by standard spirometry (Dymamic Spirometry, Microloop, Sensor Medics) following guidelines set by the American Thoracic Society and the European Respiratory Society. Reference Macintyre, Crapo and Viegi15,Reference Miller, Hankinson and Brusasco16 The forced vital capacity (FVC), the forced expiratory volume in one second (FEV1), the FEV1/FVC ratio (the bronchial permeability index or the amount of air one can forcefully exhale from the lungs), and the forced expiratory flow at 50% (FEF50) and 75% (FEF75) of vital capacity were determined.
Statistical analyses
Wherever possible, scientists measuring outcomes were blinded to the maternal preeclampsia status by coding the patients. Data are presented as the mean ± SEM or SD. Distribution was verified with the Shapiro-Wilk test, and non-parametric data were log transformed to approximate a normal distribution where necessary. Data were analysed by Two-Way ANOVA with the Tukey post hoc test. For comparison of the percentage of cases with a TR Vmax > 2.8 m/s between Control and PE women, the Chi-Squared test and the Fisher’s Exact test were used. For all comparisons, statistical significance was accepted when P < 0.05.
Results
Maternal characteristics
Values for maternal age, height, weight, body mass index (BMI) and body surface area (BSA) were not different for control and PE women (Table 1). Similarly, gravidity and primiparity were similar in control and PE women (Table 1). Post-partum, all women lost weight, BMI and BSA (Table 1).
Data presented as mean ± SD or median (IQR). *, significant difference (P < 0.01) Preeclampsia vs. Controls before and after birth; †, significant differences within Preeclampsia group; #, significant differences within Control group. Abbreviations: BMI, body mass index; BSA, body surface area; EDD, end-diastolic dimension; EDL, end-diastolic length; EF, ejection fraction by Teicholtz method; ESA, end-systolic area; ESD, end-systolic dimension; LV, right ventricle; MPI’, myocardial performance index by pulsed wave tissue Doppler; PA, pulmonary artery; PAP, pulmonary artery pressure; PG, peak gradient; PLAX, parasternal long axis view; RA, right atrium; RV, right ventricle; PVR, pulmonary vascular resistance; SI, sphericity index; SF, shortening fraction; TAPSE, tricuspid annular plane systolic excursion; TR, tricuspid regurgitation; Vmax, maximal systolic velocities; VTI, velocity time integral; 4CHV, four-chamber view. RV-RA pressure gradient calculated by application of modified Bernoulli equation: 4 x (TR Vmax)2.
Infant characteristics
Relative to babies from control mothers, babies from PE mothers were born with a lower birth weight and shorter body length, by a greater incidence of Caesarean section, at an earlier stage of gestation (Table 1). The infant sex ratio was similar between groups (Table 1).
Cardiovascular assessment
Right ventricular and pulmonary geometry . Relative to controls, PE women had significantly greater values for RV end-systolic and end-diastolic diameter, and for RV end-systolic and end-diastolic length (Fig 1 A-D and Table 1). Relative to controls, PE women had significantly lower values for RV relative wall thickness (Fig 1 E and Table 1). Relative to controls, PE women had significantly greater values for RA end-systolic diameter and area (Fig 1 F & G and Table 1). Relative to controls, PE women had significantly greater values for pulmonary artery dimension (Fig 1 H and Table 1). Combined, the data show that relative to controls, PE women had enlarged RV and RA chamber sizes and pulmonary artery dimensions, and that these differences remained ca. 6 weeks after birth (Fig 1 A-H and Table 1).
Right ventricular and pulmonary artery function . Relative to controls, PE women had significantly greater values for the TAPSE and the maximal TR Vmax (Fig 2 A & B and Table 1). Combined, these data show that relative to controls, PE women had significantly increased RV global longitudinal contractility and blood flow to the lungs (Fig 2 A & B). Relative to controls, PE women also had significantly greater values for the RV to RA pressure gradient (RV-RA PG), the pulmonary valve velocity time integral (PV VTI), estimated mean PAP, estimated systolic PAP, and estimated mean pulmonary vascular resistance (PVR; Fig 2 C-F and Table 1). Most differences remained significant ca. 6 weeks after birth (Fig 2 A-F and Table 1). Values for TAPSE, TR Vmax, RA-RV PG and mean and systolic PAP were significantly lower in PE women after compared to before birth. However, these variables remained significantly elevated after birth relative to control women (Fig. 2 and Table 1). Mean estimated PAP values in control women were 35.8 ± 0.4 mmHg before birth and 34.9 ± 0.4 after birth. Mean estimated PAP values in PE women were 45.2 ± 0.5 mmHg before birth and 40.6 ± 0.9 after birth. The British Society of Echocardiography has provided guidelines for the estimation of pulmonary hypertension based on a flow chart that assesses the probability of pulmonary hypertension using values for the maximum TR Vmax and other echocardiographic indices of pulmonary hypertension. Values for the TR Vmax > 2.8 m/s coupled with two other positive echocardiography indices are classified as an intermediate to high probability of significant pulmonary hypertension in affected patients. Reference Augustine, Coates-Bradshaw and Willis14 In our Bolivian study, the percentage of cases with a TR Vmax > 2.8 m/s with two other echocardiographic indices of pulmonary hypertension was greatly elevated in PE relative to Control pregnant women at high altitude before (80% vs. 6%) and ca. 6 weeks after (54% vs. 4%) birth (both P < 0.05).
Respiratory assessment
Relative to controls, PE women had significantly lower values for peripheral O2 saturation (Fig 3 A). Relative to controls, PE women had significantly lower values for the FVC, the FEV1, the bronchial permeability index (FEV1/FVC), and the FEF50 and FEF75 (Fig 3 B-F). These differences, except FVC and FEF50, remained significant post-partum (Fig 3 A-F). In contrast, basal respiratory rate was similar between groups pre- and post-partum (Table 1). Values for peripheral O2 saturation were lower in control and PE women after compared to before birth. Values for FEV1/FVC and FEF50 were lower in control women after compared to before birth (Fig. 3).
Discussion
In this prospective cohort clinical study in the high-altitude city of La Paz, Bolivia, the data show that relative to women undergoing healthy pregnancy at high altitude, women diagnosed with preeclampsia at high altitude had greater values for RV end-systolic and end-diastolic diameter and length. Women diagnosed with preeclampsia at high altitude also showed greater pulmonary artery dimension and maximal TR velocity, increased RV contractility and performance index. Women diagnosed with preeclampsia at high altitude also showed elevated estimated values for PAP and pulmonary vascular resistance. All outcomes measured in women diagnosed with preeclampsia at high altitude are consistent with pulmonary hypertension. Relative to women undergoing highland healthy pregnancy, women diagnosed with preeclampsia at high altitude also had lower values for peripheral oxygen saturation, forced expiratory flow and a reduced bronchial permeability index. All cardiorespiratory differences were independent of body weight, and most differences remained ca. 6 weeks after birth.
The present study was a multi-institutional clinical investigation, whose study population was pregnant women undergoing pregnancy at the high-altitude in the city of La Paz, Bolivia (3,625m; 11,893 ft). Patients were transferred from five different first, second and third-level health care centres with either a diagnosis of preeclampsia or normal pregnancy to undergo cardiorespiratory assessment at the Instituto Boliviano de Biologia de Altura (IBBA, 3,650m). All patients wished to participate in the study on a voluntary basis. All patients were born and had been permanent residents in La Paz for at least two generations.
The annular plane systolic excursion (TAPSE) is an established measure of global RV function, which describes load-dependent longitudinal contractility. Reference Alam, j, Andersson, Samad and Nordlander17,Reference Ueti, Camargo, Ueti Ade, de Lima-Filho and Nogueira18 The TR Vmax is the most widely used variable by transthoracic echocardiography in patients with suspected pulmonary hypertension. Reference Marra, Naeije, Ferrara, Vriz, Stanziola, D'Alto, D'Andrea, Carannante, Eichstaedt, Cittadini, Benjamin, Grünig and Bossone19 In presence of significant TR and normal contractility, TAPSE is increased. Reference Ueti, Camargo, Ueti Ade, de Lima-Filho and Nogueira18,Reference Marra, Naeije, Ferrara, Vriz, Stanziola, D'Alto, D'Andrea, Carannante, Eichstaedt, Cittadini, Benjamin, Grünig and Bossone19 An increase in the right ventricular-right atrial (RV-RA) pressure gradient is associated with pulmonary hypertension. Reference Casazza, Bongarzoni and Forgione20 Pulmonary arterial hypertension results in increased RV afterload and increased tricuspid valve regurgitation. Reference Farber and Loscalzo21 Therefore, increased values for the pulmonary artery dimension, TR Vmax, TAPSE, together with enlarged RV and RA chamber sizes are all strongly indicative of pulmonary hypertension. The latter was further supported by calculation of a significant increase in the estimated PAP and in the pulmonary vascular resistance in Bolivian women diagnosed with preeclampsia at high altitude relative to those undergoing healthy pregnancies at high altitude. Values for the estimated mean PAP in women diagnosed with preeclampsia in our study were much greater than the elevated normal threshold for those living at high altitude as proposed by Leon Velarde et al. Reference León-Velarde, Maggiorini and Reeves22 In addition, according to Augustine et al. Reference Augustine, Coates-Bradshaw and Willis14 and the flow chart provided in the British Society of Echocardiography guideline for the evaluation of pulmonary hypertension, calculation of the percentage of women with TR Vmax > 2.8 m/s coupled with two other echocardiography categories, also indicated a significantly greater probability of pulmonary hypertension in Bolivian women diagnosed with preeclampsia at high altitude relative to those undergoing healthy pregnancy at high altitude.
Measurements of forced expiratory volumes (FEV) and flow (FEF) are used in the diagnosis of obstructive pulmonary function. Reference Langan and Goodbred23 The bronchial permeability index or FEV1/FVC ratio indicates how much air one can forcefully exhale. The FEV1 measures how much air one can exhale in one second, and the FVC measures the total amount of air one can exhale forcefully in one breath. A reduced value for FEV1 and for the FEV1/FVC ratio is indicative of increased airway resistance to expiratory flow. Reference Langan and Goodbred23 Several studies have also confirmed increased peripheral airway obstruction in patients diagnosed with pulmonary hypertension. Reference Jing, Xu and Badesch24 Therefore, data in the present study also show that relative to women undergoing healthy highland pregnancy, women diagnosed with preeclampsia during pregnancy at high altitude are not only at greater risk of pulmonary hypertension but also of impaired peripheral airway function, leading to compromised oxygenation. The latter was again confirmed by measurement of significantly ílower values for the peripheral O2 saturation in Bolivian women diagnosed with preeclampsia relative to those undergoing healthy pregnancies at high altitude. This is striking considering that some of the inclusion criteria for patients in this study were that women should be non-smoking mothers of normal BMI and no history of cardiovascular or respiratory health complications. Although many differences in respiratory function were statistically significant, absolute differences between groups were generally modest, and the long-term implications of these differences are unclear.
Previous studies in the USA, Saudi Arabia and Bolivia have revealed that the prevalence of preeclampsia is significantly increased during pregnancy at high altitude. Reference Keyes, Armaza, Niermeyer, Vargas, Young and Moore4,Reference Mahfouz and al-Erian25,Reference Palmer, Moore, Young, Cregger, Berman and Zamudio26 . Two studies of preeclampsia in Bolivia have hypothesised that vascular toxic factors may be released into the maternal circulation by the diseased placenta, and that these factors can pass the placental barrier and affect the circulation of the offspring, increasing their cardiovascular risk in later life. Reference Jayet, Rimoldi and Stuber27,Reference Heath-Freudenthal, Toledo-Jaldin and von Alvensleben28 Augmented oxidative stress and an angiogenic imbalance have been linked with preeclampsia during pregnancy at high altitude and pulmonary abnormalities in the newborn and young children. Reference Jayet, Rimoldi and Stuber27,Reference Heath-Freudenthal, Toledo-Jaldin and von Alvensleben28 Pre-clinical studies in pregnant sheep have reported that pregnancy at high altitude also increases uterine vascular resistance and promotes maternal hypertension. Reference Hu, Dasgupta, Xiao, Huang, Yang and Zhang29 Studies in sheep, mice, rat and guinea pigs have also confirmed that it is the hypoxia of high-altitude pregnancy that promotes the adverse effects on the placenta and mother, since pregnant animals undergoing isobaric hypoxic pregnancy also show symptoms associated with preeclampsia. Reference Zhou, Xiao and Hu30–Reference Tong, Allison and Brain34 One study in sheep undergoing isobaric hypoxic pregnancy reported evidence of increased placental vascular resistance, placental oxidative stress, placental activation of the unfolded protein response and an angiogenic imbalance in the maternal circulation. Reference Tong, Allison and Brain34,Reference Giussani35 The present study extends these ideas and reports the novel link between greater risk of preeclampsia and greater risk of concurrent as well as future pulmonary hypertension in the mother at high altitude. Similar placental and maternal vascular toxic factors may provide a mechanism that underlies this link. However, this needs to be confirmed in future investigations. A study of women with severe preeclampsia at sea level in Baltimore, Maryland, USA, Reference Vaught, Kovell and Szymanski36 also showed higher estimated PAP compared to corresponding controls, consistent with the present finding of higher estimated PAP in women with preeclampsia in high altitude Bolivia. Interestingly, both control and pre-eclampsia groups in the present study in Bolivia had higher values of estimated PAP and cardiac output than preeclamptic women in the Baltimore study, which may be related to living at altitude.
The present study has several strengths and limitations. Strengths include that this was a prospective clinical study of cardiovascular as well as respiratory outcomes with follow-up after pregnancy. The follow-up data showed that significant differences between control and preeclamptic women remained at ca. 6 weeks after birth. The study was also multi-institutional and conducted in highland women of homogenous mestizo ethnic background who were residents at ∼ 3625 m or 11,893 ft in La Paz for at least two generations derived from five maternity hospitals. Since an increase in body weight is related to an increase in cardiovascular risk, including pulmonary hypertension, matching control and preeclamptic women with similar values for body weight and body mass surface area in addition to gestational age at the prenatal ultrasound examination was also important. Another strength was that estimated PAP as well as pulmonary vascular resistance were calculated. Finally, the echocardiography measurements were performed by a single expert (C.E.S.) using the same machine and probe, and every outcome was the average of several measurements to satisfy intra-observer variability. The limitations of the study are just as important. They include the lack of pre-pregnancy measurements and longer-term follow-up with cardiorespiratory measurements past 6 weeks. The study also did not collect maternal plasma samples, which could have been used to give greater insight into underlying mechanisms.
In conclusion, the data in the present study support the novel hypothesis tested that during highland pregnancy, women who develop preeclampsia are at greater risk of pulmonary hypertension compared to women undergoing healthy highland pregnancy. Such information is important for exploring whether there are common underlying pathways between preeclampsia and future cardiovascular risk, and for designing public health prevention initiatives in women with an increased cardiovascular risk, particularly in those who are residents at high altitude. Since preeclampsia and maternal pulmonary hypertension are independent risk factors for maternal mortality, women undergoing pregnancy at high altitude should be counselled regarding the risks, and a multi-professional approach with expert care in pulmonary hypertension centres should be adopted to improve outlook and reduce risks. Such an approach will hopefully lead to a decrease in the alarming rates of maternal mortality in high-altitude countries, particularly in Bolivia.
Data availability statement
All data supporting the results are presented in the manuscript.
Acknowledgements
We would like to thank all the Bolivian women who volunteered to take part in this study. We also acknowledge the invaluable help and work of Dr Patricio Andrade, Dr Nelson Patiño, Dr Jhonny Gonzales, Dr Patricio Gutierrez, Dr Luis Niño de Guzman, Dr Christian Aguilar Ledezma, Dr Carlos Loayza, Dr Jaime Huaycho, Dr Daniel Cardenas, Dra. Nancy Villena, Dra. Mercedes Villena, Dra. Virignia Condori, Sra. Ivonne Contreras Tapia, Sra, Ana Maria Alarcon, Sra. Cristina Gonzales Davila and Sr. Daniel Mamani Baltazar. We also would like to acknowledge the Centro de Salud Integral Tembladerani, the Red de Salud Sur Oeste Numero 1, the Hospital de la Mujer, Hospital, the Hospital Materno Infantil CNS, the Hospital de Cotahuma and the Hospital Nuestra Señora de La Paz.
Author contribution
Conceptualisation, C.E.S, O.V.P., D.A.G.; Methodology, C.E.S., O.V.P., C.M., M.G., V.E., S.M., R.R., R.V., Y.P., J.M., L.T., A.B., J.J., J.P., C.R., M.A., L.R., D.A.G. Writing – Original Draft, C.E.S, O.V.P., D.A.G. Writing – Review & Editing, C.E.S, O.V.P., D.A.G.; Methodology, C.E.S., O.V.P., C.M., M.G., V.E., S.M., R.R., R.V., Y.P., J.M., L.T., A.B., J.J., J.P., C.R., M.A., L.R., D.A.G. Visualisation, C.E.S, O.V.P., D.A.G.; Methodology, C.E.S., O.V.P., C.M., M.G., V.E., S.M., R.R., R.V., Y.P., J.M., L.T., A.B., J.J., J.P., C.R., M.A., L.R., D.A.G. Project Administration, C.E.S, O.V.P., D.A.G.; Methodology, C.E.S., O.V.P., C.M., M.G., V.E., S.M., R.R., R.V., Y.P., J.M., L.T., A.B., J.J., J.P., C.R., M.A., L.R., D.A.G. Funding Acquisition, C.E.S, D.A.G.
Funding
This work was supported by Embajada Francesa and La Cooperación Francesa in Bolivia (CIF: UAP/FRA/1272/2015; CES) and the British Heart Foundation (FS/02/072/13469 and RG/06/006/22028; DAG).
Competing interests
The authors declare no competing interests.
Ethical approval
The study was carried out according to standards set by the Declaration of Helsinki and all procedures were approved by the Local Ethics Committee of the Bolivian Institute for High Altitude Biology (Consejo Técnico, Instituto Boliviano de Biología de Altura, Universidad Mayor de San Andrés, La Paz, Bolivia).