Abstract
We used meta-analysis to measure the effect of high-altitude climate therapy (HACT) on lung function outcomes in asthma, and systematically searched PubMed, Embase and www.elibrary.ru for publications appearing from 1970 to mid-2015. We included studies carried out with children or adults with an exposure of up to 12 weeks at an altitude of ≥1500 m above sea level. Changes in forced expiratory volume in 1 s (FEV1), FEV1/vital capacity ratio or peak expiratory flow rate as the HACT intervention outcomes were analysed. We included data for 907 participants (age range 4–58 years) from 21 studies, altogether including 28 substrata based on asthma type or severity. Only three of 21 included studies had high quality, whereas 93% of substudies reported lung function improvement with an overall pooled standardised mean difference (SMD) of 0.53 (95% CI 0.43–0.62). The measured effect of HACT was greater in adults (SMD 0.75, 95% CI 0.63–0.88, n=14) than in children (SMD 0.24, 95% CI 0.09–0.38, n=14). Studies at altitude >2000 m above sea level yielded the same effect as those at lower altitude. Based on a cut-point of a 0.50 change in SMD to define a meaningful clinical difference, HACT appears to have efficacy as an intervention. This extent of benefit appears to be limited to adults with asthma.
Abstract
High-altitude alpine therapy may be an effective intervention to improve lung function in patients with asthma http://ow.ly/u3i23008vU5
Introduction
Environmental modification remains a key component of asthma management, ranked together with appropriate maintenance pharmacologic therapy in stepwise treatment protocols [1]. Environmental interventions often begin with household allergen reduction (e.g. dust mite controls) or irritant removal (e.g. avoidance of second-hand cigarette smoke). More advanced environmental modification strategies also have an established role in the care of children and adults with asthma. Examples include speleotherapy, hypobaria and, most notably, high-altitude climate therapy (HACT, also known as “alpine” therapy). Historically, this latter treatment is the most widely used single modality of such environmental modulation for asthma [2]. Due to inconvenience, expense and questionable efficacy, none of these environmental interventions is currently accepted widely, including HACT. Nevertheless, there are biological mechanisms through which a beneficial effect of HACT is plausible, including reduced antigen exposure [3, 4], immune modulation [5, 6] and lowered air viscosity possibly ameliorating airway resistance.
Multiple investigators have studied the effects of HACT at various treatment sites around the world, characteristically studies that were not blinded or properly randomised and that provided single-arm, post-treatment to pre-treatment comparisons. Elements of the HACT literature are poorly accessible for critical review, however; in particular due to non-English language publication. These accumulated data present a rich opportunity for a systematic assessment of the presence or absence of treatment efficacy for HACT. We wished to address the existing knowledge gap in the treatment effect of HACT for asthma by subjecting data from these studies to a formal meta-analysis, including publications that had not appeared in the English language literature.
Methods
Data sources and searches
Two researchers (D.V. and A.K.) independently performed the initial literature search and abstract screening. For differences in initial assessment necessitating adjudication, we strove for consensus among all four co-authors. The search was completed in August 2015 and utilised three databases: PubMed, Embase and a major Russian resource, www.elibrary.ru. The search was limited to publications appearing from 1970 onwards. We used the following search term strategies: “asthma” AND “altitude” in Embase; asthma(Title/Abstract) AND altitude(MeSH Terms “altitude” including subheadings: adverse effects, epidemiology, history, instrumentation, pharmacology, physiology, statistics and numerical data, therapeutic use, and therapy” OR “altitude sickness” OR “pulmonary edema of mountaineers”) in PubMed; and “asthma” AND “high-altitude” OR “high-altitude climate therapy” OR “alpine therapy” (in Russian) for www.elibrary.ru. The results of the search algorithm using the PRISMA protocol, and consistent with its checklist and steps of study identification and inclusion, are shown in figure 1 [7].
Search algorithm. PFT: pulmonary function test.
Study selection, data extraction and assessment of risk of bias
PubMed, Embase and www.elibrary.ru were searched, returning 342 cited publications that appeared from 1970 onwards. We considered eligible studies that represented original research on human asthma treatment involving exposure to conditions of altitude ≥1500 m above sea level (MASL). There was no age restriction: studies both of adult and childhood asthma were considered eligible. We also searched the references cited by papers in a systematic review [8]. Of 342 citations identified initially, 307 (90%) did not meet the eligibility criteria at the initial screening stage. The remaining 35 studies underwent further review. Critical review of complete texts of the 35 eligible studies identified two publications [9, 10] that were excluded because they were based on data ultimately subsumed in two other publications [11, 12]. We also found that a study set of one publication [13] was included within a larger cohort [14], so only the latter was included. This produced 32 eligible papers for this review. At this stage of review, we further excluded initially eligible studies that did not include at least one of the following measures derived from pulmonary function testing: forced expiratory flow in 1 s (FEV1), the ratio of FEV1 to forced vital capacity (FVC) or vital capacity (VC), or peak expiratory flow rate (PEFR). There were nine studies excluded due to a lack of such pulmonary function data. We additionally excluded two other studies that only examined longer-term effects of HACT (>12 weeks duration at altitude). After all exclusions, 21 publications remained in the final meta-analysis.
Quality assessment using available scales was hard to accomplish because of significant heterogeneity arising from absence of uniform criteria for performing such studies across different locations and time. Therefore, we set our own criteria to assess quality of studies, and minimum quality parameters for study inclusion were pulmonary function measurements before and during altitude, no evidence of dropout rates >60%, no other major reported change in therapy, and altitude of at least 1500 MASL.
Data synthesis and analysis
The end-point for analysis was one of the measures of pulmonary function required for eligibility. For studies that reported more than one pulmonary function end-point, we preferentially selected the FEV1 % pred for analysis. If the raw FEV1 only was available rather than as a percent predicted, we extracted that outcome. If these were not provided, we selected the PEFR (actual or as a percent predicted). The FEV1/FVC (or VC) ratio was used only if neither of the other two measures was available. The pool of identified studies included both randomised controlled studies and single-arm studies without a control group. For consistency, we extracted the pulmonary function end-point of interest before and after altitude exposure and did not utilise data from control study arms (i.e. study of persons without asthma at altitude). We calculated the mean difference in end-point between the baseline pre-treatment value and after the high-altitude intervention. The pre-treatment value could represent measurement at low altitude or immediately upon arrival at high altitude. Several studies stratified patients based on asthma severity (mild, moderate or severe) or aetiology (atopic and nonatopic asthma), reporting grouped pulmonary function data within such strata. In those cases, because they represented independent samples (an assumption of the meta-analysis model), we treated each stratum as a separate substudy. Taking such strata into account, among the 21 studies included we ultimately analysed 28 end-points.
The effect size was calculated as the mean change in the selected pulmonary function end-point subtracting the baseline value from the post-treatment value. Thus, a positive mean value represented a treatment-associated improvement in lung function. In addition to studies reporting effect sizes by varying scales (i.e. using differing lung function indices), we calculated the standardised mean differences (SMDs). A few studies reported the standard error rather than the standard deviation. For those studies, we derived the standard deviation from the reported standard error using the formula: sd=se×n1/2. For a single study reporting the median value of pulmonary function end-point together with its 95% confidence interval, we calculated the standard deviation assuming the mean to be the mid-point of the 95% confidence interval, using the formula: sd=n1/2×((95% upper bound–mean)/1.96). We used Stata 14.1 software to run a meta-analysis function to calculate SMDs. We recalculated SMDs in stratified analyses based on age, altitude and duration of study. In a sensitivity analysis, we calculated the nonstandardised outcome for the pool of studies reporting mean differences in FEV1 % pred. Further sensitivity analyses recalculated the pooled SMD after selected exclusions on studies to assess their impact on the overall results. Publication bias was assessed using Egger's and Begg's tests [15, 16]. The degree of heterogeneity due to between-study variability was estimated using the Chi-squared test with the relevant p-value. For interpretation purposes, heterogeneity was also supplemented with the I2 metric [17]. We also report the tau value (τ2) in absolute units as an estimate of heterogeneity among the studies.
Results
Table 1 summarises the 21 studies included in the analysis. These subsumed, in total, 28 strata treated as independent substudies for the purposes of the meta-analysis. All but four studies reported the FEV1 % pred, allowing us to analyse this in the majority of cases (17 of 21 included studies; 21 of 28 substudy strata). Data for raw FEV1 and for PEFR were available from one study each, and the FEV1/VC ratio from two others (altogether end-points other than FEV1 % pred from seven substudy strata). In total, data for 907 participants were included, of whom 568 (63%) were adults. 13 studies were conducted in Western Europe (Italy and Switzerland), four in Eastern Europe, three in Nepal and one in South America. The publication dates for the studies ranged from 1973 to 2013; 12 (57%) were published before 2000. In all but one study [11], pre-altitude asthma pharmacotherapy was continued when exposure to altitude was initiated. Such baseline therapy was heterogeneous, varying from patient to patient and over time with changing standards of care. Short-acting β-agonists were typical across studies as rescue medication.
Characteristics of the 22 studies including 29 substrata included in the meta-analysis
As shown in figure 2, the SMD post-treatment ranged from a fall in lung function of −1.79 (95% CI −3.31–−0.27) to an increase of 2.96 (95% CI 2.36–3.56) standardised units. Heterogeneity across studies was high (τ2=1.58, I2=85%, p<0.001). Despite this statistical heterogeneity, the studies were mostly consistent in their direction of effect. In total, 26 of 28 strata (93%) observed a positive response (SMD >0), while a lung function reduction was observed in only two strata. Of note, in both of these, which accounted for only 3.2% of the overall study weight, the exposure was to altitudes >4000 MASL [31, 33]. Of the 26 positive studies/substudies, 13 yielded an SMD >0 but <0.50, while 13 were at or above that value, which is an accepted threshold for a meaningful clinical difference (MCD) [37, 38]. Of these, eight yielded an SMD of ≥0.50 or more but ≤1.0, while for six, the SMD was >1.0. The studies/substudies yielding a positive treatment effect size of ≥0.5 SMD subsumed 403 of the 907 (44%) subjects in the total study pool. There was a trend to more attenuated effects among the publications with a more recent year of publication (figure 2).
Pooled meta-analysis of included studies. Diamonds indicate the standardised mean difference (SMD) and shaded boxes represent the individual study weights.
In order to investigate further potential differences in HACT efficacy by intervention and by study population, we carried out stratified analyses by study site altitude, age group (childhood versus adult) and treatment duration. These stratified analyses are shown in table 2. Six studies/substudies carried out at an altitude >2000 MASL reported a lung function improvement (SMD 0.53, 95% CI 0.32–0.74) that was similar to that seen in the 22 studies/substudies carried out at an altitude of ≤2000 MASL (SMD 0.53, 95% CI 0.42–0.64)). Stratification by age (children versus adults) yielded a positive effect in both groups, although the estimated SMD was greater in adults (SMD 0.75, 95% CI 0.63–0.88) compared with that among children (SMD 0.24, 95% CI 0.09–0.38). Treatment for >4 weeks was associated with a marginally greater effect size than more brief treatment (SMD 0.63, 95% CI 0.50–0.77 versus SMD 0.42, 95% CI 0.28–0.55, respectively). In the single study in which baseline medicines were discontinued, the observed SMD was small and not statistically significant (SMD 0.10, 95% CI −0.73–0.94).
Re-analysis stratified by altitude, age and duration of high-altitude climatic therapy
We performed a sensitivity analysis excluding studies reporting outcome measures other than FEV1 % pred (n=7 studies/substudies excluded). After this exclusion, the overall SMD was somewhat lower (SMD 0.34, 95% CI 0.24–0.45). As these studies reported effect using the same scale, we also calculated the mean difference in this stratum without standardisation. At baseline, FEV1 % pred ranged from 64.5% to 105.8% (median 90.2%). The overall FEV1 % pred increase was 5.79% (95% CI 4.36–7.23) with moderate heterogeneity (τ2=0.29, I2=58%, p<0.001). A forest plot of this analysis is presented in figure 3. In a 170 cm 30-year-old adult male, this would equate to 229 mL (assuming FEV1 100% predicted to be 3.95 L) [36]. We also carried out another sensitivity analysis by re-estimating the effect of HACT excluding each individual study one by one. This re-analysis showed that the study by Dubiley et al. [18] made the biggest contribution to the overall SMD: when excluded (all three substudies from that publication combined), the overall SMD fell to 0.35, but was still statistically significant (95% CI 0.25–0.46) (table 3). In the opposite direction, the study of Verkleij et al. [14] brought the pooled SMD down to the greatest degree: when excluded, the overall SMD increased to 0.60 (95% CI 0.49–0.70). When excluded, the only study that had reported lung function as a median rather than as a mean value [30] did not substantively change the observed overall SMD (0.55, 95% CI 0.45–0.65). We also wished to calculate SMD for studies performed at altitudes of >2500 MASL, which is an accepted threshold for high altitude. Six remaining substudies from altitudes of 3200–6965 MASL yielded SMD 0.53 (95% CI 0.32–0.74). Finally, when studies that used control group designs were excluded in order to determine if these had impacted the findings even though the control arm data had not been included [26, 31, 33], the pooled SMD slightly increased to 0.59 (95% CI 0.48–0.69). There was no strong evidence of systematic publication bias (figure 4). In addition to the visual pattern displayed in the funnel plot, Begg's test (p=0.47) and Egger's test (p=0.86) showed no clear evidence of publication bias.
Pooled meta-analysis of studies reporting forced expiratory volume in 1 s % pred. Diamonds indicate the weighted mean difference (WMD) and the shaded boxes represent the individual weight of each study.
Sensitivity analysis
Funnel plot (with pseudo 95% CI) of 28 substrata analysed.
Discussion
This comprehensive quantitative assessment of the effect of HACT on asthma using a standard meta-analytic approach found a statistically significant improvement in lung function following this intervention. Moreover, the positive effect is present among both adults and children with asthma. Particularly noteworthy, in aggregate the magnitude of this effect meets a standard criterion for an MCD, although this benchmark was achieved only among adults. This MCD metric is defined a ≥0.50 SMD effect size [37, 38]. Indeed, this magnitude of effect is very close to that observed in asthma in clinical trials of moderate doses of inhaled beclomethasone (SMD 0.49, 95% CI 0.30–0.69) [39]. Using the metric of FEV1 absolute change (standardised to an adult male) available for 21 of the included substudies, HACT appears to be almost as potent compared with the maximum estimated pooled effect of inhaled budesonide for asthma in clinical trials (229 versus 290 mL) [40].
Our meta-analysis has a number of strengths. We used a search strategy targeting inclusion of studies published in both English and languages other than English, spanning a wide time range for inclusion. We were conservative in our approach by relying on pulmonary function-defined outcomes as the objective measure of clinical improvement in asthma. By limiting the pool of studies to those ≥1500 MASL we also insured that the HACT intervention was at least moderate in its environmental attributes. We acknowledge, however, that some of the studies carried out at the higher altitudes included were not specifically designed to evaluate the effect of HACT, but rather designed as a respite sojourn for disease sufferers. Our retention of study subsets in the analysis rather than collapsing these into single values for each included publication should be considered a study strength rather than a limitation, because this conservative approach enabled estimation of the effect of HACT in different asthma subgroups.
This meta-analysis approach has a number of inherent limitations. We only considered pulmonary function as a measure of clinical effect, although other clinical and functional criteria are used widely to assess asthma. Bronchial hyperreactivity is one of these. Although this was reported in a limited number of the studies that we accessed, there were too few to contribute to an effective meta-analysis strategy. In many but not all of these studies, nonspecific hyperresponsiveness decreased after the intervention, consistent with a salutary effect of HACT on asthma. Another potential limitation in our study is the use of FEV1 as an end-point in children because it may be a poor measure of asthma control, especially in those of younger age [41, 42]. Nevertheless, HACT has been used as ancillary treatment for asthma in the paediatric population, as the number of publications included in this analysis indicate. Furthermore, many of the nonadults included were, in fact, adolescents, thus limiting the likely impact on our findings of using the FEV1 as an asthma measure in nonadults. The studies we included were limited to those conducted at altitude >1500 MASL. We have elected to exclude studies <1500 MASL to reduce heterogeneity identified in the previous qualitative review [8].
Critically, we did not assess the potential long-term benefit of HACT after return to lower altitude or in combination with other treatments (e.g. standard stepwise pharmacologic therapy). In many of the studies included in our analysis the participants were on such therapies. Further, we did not have access to the direct and indirect economic costs of HACT, such that this could be compared to other interventions. We also did not study the effects on asthma of dwelling longer-term at high altitude. Another limitation is that many of the studies that we identified were relatively small in size. This, however, is an argument in favour of a meta-analysis strategy, since statistical power is greater and more subtle associations can be identified when results from multiple small studies can be combined in meta-analysis. Although small studies tend to be impacted by publication bias to a greater degree than larger studies [15, 43], we did not find evidence of that in this analysis. Selection bias may be present across all of the selected publications because HACT is often used as a supplementary treatment in relatively stable asthma, whereas patients with severe or difficult to control asthma may be excluded from such investigations. There also may be information bias (disease misclassification) because publications are included that predate current asthma diagnosis and treatment standardisation. Finally, unmeasured confounding cannot be excluded, including unaccounted for climatic or other geographic parameters.
A recent systematic review assessed the effects of HACT on lung function and airway hyperreactivity using a systematic review but without a statistical meta-analytic component [8]. That review concluded that there was insufficient evidence to establish efficacy, invoking heterogeneity of studies and small sample sizes. Of note, that review was limited to studies at altitudes of 1500–2500 MASL published in English and German only, and included longer duration studies that may have added additional heterogeneity. We also found statistically significant heterogeneity in our meta-analysis. Given the variety of different study designs, populations and research methods used, statistical heterogeneity is expected. In contradistinction to this statistical phenomenon, it is important to note that all of the 26 studies carried out at an altitude of <4000 MASL consistently observed lung function improvement, a homogeneous direction of effect. The findings of another relatively recent systematic analysis of HACT as an intervention to treat atopic dermatitis, a condition that has some shared attributes with asthma, are also relevant [44]. That review concluded that there was a positive HACT-associated treatment effect. In a broader sense, it also is important to consider whether a meta-analytic approach is appropriate to the data at hand, given that HACT is an environmental intervention that potentially is heterogeneous, and varies by elevation and in its effects among different asthma subgroups. Our results, including the overall findings as illustrated in the forest and funnel plots and in the sensitivity analyses that we performed, indicate that the estimates yielded do provide a useful assessment of the effect of HACT on lung function in asthma.
We found that HACT may be an effective intervention to improve lung function in asthma. Nonetheless, a number of questions remain to be addressed including optimal altitude and duration of treatment, and most importantly, strategies to extend a retained benefit upon return to the home environment and subgroups (e.g. adults rather than children) most likely to benefit from this intervention. Even when geographic conditions allow (i.e. available mountainous areas), HACT is a relatively resource-intensive intervention that may limit its application and be subject to inequalities of access. This clear direction and magnitude of effect shown in this meta-analysis may help to stimulate further inquiry into these outstanding research topics.
Acknowledgements
Authors would like to thank the authors who have provided additional data on their studies in reply to enquiries (Oktam I. Bobokhodzhaev (Ministry of Health of Tajikistan, Dushanbe, Tajikistan) and Anzor T. Kokov (“Elbrus” Sanatorium, Nalchik, Russia)). We also thank Michaela Hornakova (University Library, Bratislava, Slovakia), Professor Sergey L. Scherbin (Kabardino-Balkarian State University, Nalchik, Russia) and Valery S. Shirinsky (Scientific Institute of Fundamental and Clinical Immunology, Novosibirsk, Russia) for their help in finding original references in the libraries.
Footnotes
Conflict of interest: None declared.
- Copyright ©ERS 2016
This article is open access and distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0.