Article Text

Bronchodilators accelerate the dynamics of muscle O2 delivery and utilisation during exercise in COPD
Free
  1. Danilo C Berton1,
  2. Priscila B Barbosa1,
  3. Luciana S Takara1,
  4. Gaspar R Chiappa1,
  5. Ana Cristina B Siqueira1,
  6. Daniela M Bravo1,
  7. Leonardo F Ferreira2,
  8. J Alberto Neder1
  1. 1Pulmonary Function and Clinical Exercise Physiology Unit (SEFICE), Division of Respiratory Diseases, Department of Medicine, Federal University of Sao Paulo (UNIFESP), São Paulo, Brazil
  2. 2Department of Physiology, University of Kentucky, Lexington, Kentucky, USA
  1. Correspondence to Professor J A Neder, Pulmonary Function and Clinical Exercise Physiology Unit (SEFICE), Respiratory Division, Department of Medicine, Federal University of São Paulo (UNIFESP), Rua Professor Francisco de Castro 54, Vila Clementino, São Paulo CEP: 04020-050, Brazil; albneder{at}pneumo.epm.br

Abstract

Background Expiratory flow limitation and lung hyperinflation promote cardiocirculatory perturbations that might impair O2 delivery to locomotor muscles in patients with chronic obstructive pulmonary disease (COPD). The hypothesis that decreases in lung hyperinflation after the inhalation of bronchodilators would improve skeletal muscle oxygenation during exercise was tested.

Methods Twelve non- or mildly hypoxaemic males (forced expiratory volume in 1 s (FEV1)=38.5±12.9% predicted; Pao2>60 mm Hg) underwent constant work rate cycle ergometer exercise tests (70–80% peak) to the limit of tolerance (Tlim) after inhaled bronchodilators (salbutamol plus ipratropium) or placebo. Muscle (de)oxygenation (∼fractional O2 extraction) was determined in the vastus lateralis by changes (Δ) in the deoxyhaemoglobin/myoglobin signal ([HHb]) from near-infrared spectroscopy, and cardiac output (QT) was monitored by impedance cardiography.

Results Bronchodilators reduced lung hyperinflation and increased Tlim compared with placebo (454±131 s vs 321±140 s, respectively; p<0.05). On-exercise kinetics of QT and pulmonary O2 uptake (V˙o2) were accelerated with active treatment; Δ[HHb] dynamics, however, were delayed by ∼78% and the signal amplitude diminished by ∼21% (p<0.01). Consequently, the ratio between V˙o2 and Δ[HHb] dynamics decreased, suggesting improved microvascular O2 delivery (τ-V˙o2/MRT-Δ[HHb]=4.48±1.57 s vs 2.08±1.15 s, p<0.05). Of note, reductions in lung hyperinflation were related to faster QT kinetics and larger decrements in τ-V˙o2/MRT-Δ[HHb] (p<0.01).

Conclusions Decreases in operating lung volumes after the inhalation of bronchodilators are associated with faster ‘central’ cardiovascular adjustments to high-intensity exercise with beneficial consequences on muscle oxygenation in patients with moderate to severe COPD.

  • Bronchodilator
  • chronic obstructive pulmonary disease
  • exercise tolerance
  • oxygen consumption
  • respiratory mechanics

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Introduction

Exercise intolerance is a multifactorial construct in patients with chronic obstructive pulmonary disease (COPD).1–3 In fact, recent studies from our laboratory suggest that expiratory flow limitation and lung hyperinflation are related to sluggish central cardiovascular adjustments at the onset of exercise which might impair the balance between O2 delivery (QO2) and uptake (V˙o2) in the exercising muscles of patients with COPD.4 5 The result, then, is faster and heightened muscle deoxygenation measured by near-infrared spectroscopy (NIRS) in patients compared with age-matched controls.4 Interestingly, these abnormalities were closely related to slower kinetics of cardiac output (QT), suggesting an important role for the central cardiovascular adjustments in setting the limits of increase in peripheral QO2 at the onset of exercise. Further experimental evidence for this contention was obtained in a subsequent study in which a strategy aimed to reduce the resistive work of breathing and dynamic hyperinflation (heliox) simultaneously accelerated the dynamics of QT and QO2 in patients with moderate to severe COPD.5 These data prompted the hypothesis that similar findings would be observed in response to the most effective and clinically feasible intervention available to alleviate lung hyperinflation and improve exercise tolerance in these patients—that is, inhaled bronchodilators.1 6 7

In the present study, therefore, we aimed to investigate the effects of acute inhalation of a combination of a short-acting β2-adrenoceptor agonist (salbutamol sulfate) and a short-acting anticholinergic (ipratropium bromide) on key determinants of QO2 and V˙o2 during high-intensity, constant work rate cycling exercise in patients with moderate to severe COPD. Confirmation of the hypothesis that bronchodilators would slow the dynamics and decrease the amplitude of leg muscle deoxygenation would provide novel evidence that amelioration of the mechanical burden of breathing might positively impact upon the provision of O2 to the working peripheral muscles during exercise in patients with COPD.

Methods

Subjects

Twelve males with moderate to severe, stable COPD (forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) <0.7 and postbronchodilator FEV1<60% predicted), and resting Pao2>60 mm Hg, volunteered to participate in the study. Patients were recruited from the COPD Outpatients Clinic of the Sao Paulo Hospital (Brazil), a University-based teaching hospital. Subjects were free of echo Doppler cardiographic evidence of severe pulmonary hypertension and left ventricular dysfunction. No patient has ever been enrolled in a pulmonary rehabilitation programme. All participants signed a written informed consent form. The study protocol was approved by the Institutional Medical Ethics Committee.

Study protocol

This was a double-blinded, placebo-controlled crossover study. After determination of peak (V˙o2) on a maximal incremental cycle ergometer exercise test, subjects were randomly assigned to inhale salbutamol sulfate 120 μg plus ipratropium bromide 20 μg per actuation, via a metered dose inhaler (MDI; Combivent, Boehringer Ingelheim, Germany) or Combivent placebo via an MDI each on a different day. At these visits, patients performed a total of four transitions from baseline pedalling to 70–80% peak work rate. The initial test was performed to the limit of tolerance (Tlim, min) 20 min after the inhalation of bronchodilators or placebo. Patients rested for 1 h before performing an additional 4 min exercise test. Thus, there were two on-exercise repetitions per day. Data from same-day transitions were averaged to reduce data noise for the kinetic analyses (see details below). Additional details are described in the supplement available online.

Measurements

Pulmonary function tests

Spirometry, single-breath lung diffusing capacity, static lung volumes by body plethysmography and arterial blood gases were measured at baseline. Recorded values were compared with those predicted for the adult Brazilian population.8 9

Exercise tests

Standard metabolic and ventilatory responses were measured breath-by-breath using a calibrated, computer-based system (CardiO2 System, Medical Graphics, St Paul, Minnesota, USA). A period of 3 min unloaded exercise preceded the imposition of the selected work rate. Serial inspiratory capacity (IC, litres) manoeuvres were performed during the test. Assuming that total lung capacity (TLC) remains constant during exercise, IC manoeuvres provide an estimate of end-expiratory lung volume.10

Skeletal muscle oxygenation

Skeletal muscle oxygenation profiles of the left vastus lateralis were evaluated with a commercially available NIRS system (Hamamatsu NIRO 200, Hamamatsu Photonics KK, Japan). The (deoxyhemoglobin (Hb)/myoglobin (Mb) signal ([HHb], μM/cm) during exercise has been considered a proxy of fractional O2 extraction in the microcirculation, reflecting the balance between O2 delivery and utilisation11 12—that is, Δ[HHb]V˙o2/QO2. Additional details are described in the supplement online.

Central haemodynamics

Cardiac output (QT, l/min) and stroke volume (SV, ml) were measured throughout the constant work rate test using impedance cardiography (PhysioFlow PF-0, Manatec Biomedical, Macheren, France).13 This methodology is different from previously used impedance systems as its algorithm does not require basal thoracic impedance measurement (Z0). Considering, however, the controversy on the accuracy of the system in providing reliable absolute values of SV in patients with COPD,14 15 exercise haemodynamic data were also expressed as ‘Δ’ from baseline. Additional information on response characteristics and system validation are detailed in the supplement online.

Kinetics analysis

The breath-by-breath V˙o2, Δ[HHb] and haemodynamic (QT, SV and heart rate (HR)) data from same-day exercise bouts were linearly interpolated to provide second-by-second values and, for each individual, the repetitions were time aligned to the start of exercise and ensemble averaged (SigmaPlot 10.0, Systat Software). The kinetic behaviour of haemodynamic data was estimated by calculation of the half-time of response (t1/2, s). After checking that a slow component was not discernible in the first 180 s of exercise, metabolic and NIRS data were fitted by the following monoexponential equation from the start of exercise to the third minute16:

[Y](t)=[Y](b)+Ap(1e(tTDp)τp)

where b and p refer to baseline unloaded cycling and primary component, respectively; A, TD and τ are the amplitude, time delay and time constant of the exponential response. For V˙o2 analysis we deleted the data relative to the first 20 s after exercise onset—that is, the cardiodynamic phase. Therefore, τV˙o2 represents the time course of the primary component—an estimate of the muscle V˙o2 kinetics.17 The overall kinetics of [HHb] were determined by the mean response time (MRT=τ+TD). The ratio τ-V˙o2/MRT-[HHb] was used as a semiquantitative index of microvascular QO2 kinetics, with higher values indicating slower QO2.4 5 Based on our experience with NIRS in patients with COPD,4 5 18 and previous findings with other disease populations,19 we systematically examined whether there was an ‘overshoot’ in Δ[HHb] after the initial fast response (see the Results section). In the case of a discernible ‘overshoot’ in the two interventions, we used the area under the ‘overshoot’ as an additional index of impaired microvascular QO2.12 Additional information on kinetic analyses are detailed in the supplement online.

Statistical analysis

Statistical analysis was performed using SPSS 15.0. We planned to evaluate 12 patients based on our previous results in a similar transversal and crossover study with heliox in patients with COPD.5 Paired t test was used to contrast within-subject exercise responses with the exception of symptom scores (Wilcoxon signed-rank test). A one-way repeated measures analysis of variance (ANOVA) was used to compare Δ[HHb] at quartiles of isotime—that is, the shortest Tlim between the two interventions on a given subject. Pearson product–moment correlation was used to assess the level of association between continuous variables. The level of statistical significance was set at p<0.05 for all tests.

Results

Subject characteristics and maximum exercise capacity

On average, patients had severe airflow obstruction and lung diffusing capacity for carbon monoxide (DLCO) impairment with increased static lung volumes and normal arterial blood gases at rest (table 1). All patients showed reduced maximal exercise capacity (peak V˙o2<83% predicted) (table 1)20 with increased peak V˙E/MVV ratio (>0.8). Although dyspnoea was the main limiting symptom in 8/12 patients, Borg scores for leg effort did not differ more than 2 units from the Borg dyspnoea ratings in these patients.

Table 1

Resting characteristics and responses to incremental exercise (n=12)

Effects of bronchodilators on pulmonary function and exercise tolerance

Treatment with bronchodilators increased FEV1 and IC compared with placebo (1.30±0.50 litres vs 1.17±0.56 litres and 2.28±0.83 litres vs 1.89±0.73 litres, respectively; p<0.001 for both comparisons). The criteria for a ‘positive’ bronchodilator response (FEV1 change ≥12% and 200 ml)21 were obtained in 7/12 patients, with FEV1 increases being closely related to changes in FVC.

The subjects exercised at 53±10 W. Exercise tolerance (Tlim) was improved after treatment with bronchodilators (321±140 s vs 454±131 s, for placebo and bronchodilators, respectively; p=0.01). Increases in pre-exercise IC were significantly related to Δ (bronchodilators–placebo) TLim (r=0.72 (95% CI 0.31 to 0.91); p=0.04). As previously described,13 increases in IC during exercise were accompanied by larger tidal and inspiratory reserve volumes and lower dyspnoea scores at isotime (table 2). In contrast, there were no significant effects of bronchodilators on mean Spo2 values (table 2), and improvements in Spo2 >3% at isotime were found in only 3/12 patients. Similarly, there were no significant between-interventions differences in the cardiovascular responses at the individual isotime (table 2).

Table 2

Isotime values of selected physiological and perceptual responses during high-intensity, constant work rate exercise after placebo or bronchodilators (n=12)

Kinetics response to exercise

Treatment with bronchodilators accelerated the central haemodynamic responses at the onset of exercise (t1/2QT= 75.9±10.3 s vs 58.9±18.9 s (p=0.02), t1/2HR= 78.2±13.0 s vs 62.5 ±15.5 s (p=0.03), t1/2SV=51.0±8.1 s vs 40.6±10.3 s (p=0.02) for placebo and bronchodilators, respectively). In relation to the V˙o2 response, there were no significant differences in baseline (unload exercise) and amplitude comparing bronchodilators and placebo (508±147 ml/min vs 528±117 ml/min (p=0.87) and 554±98 ml/min vs 532±108 ml/min (p=0.90)). In contrast, the dynamics of the ‘primary’ component of V˙O2 were significantly faster after bronchodilation (table 3).

Table 3

Comparative effects of inhaled placebo vs bronchodilators (BDs) on selected parameters of the kinetics of pulmonary oxygen uptake Embedded Image and deoxyhaemoglobin (Δ[HHb]) in the vastus lateralis at the onset of high-intensity, constant work rate exercise (n= 12)

Bronchodilators slowed the kinetics of leg muscle deoxygenation compared with placebo (figure 1, table 3), despite no significant effects on baseline Δ[HHb] (3±15 μM/cm vs 2±13 μM/cm; p=0.83). Therefore, τΔ[HHb] postbronchodilators was, on average, twofold the postplacebo value (table 3). Consequently, the τ-V˙o2/MRT-Δ[HHb] ratio, an index of impaired microvascular QO2 4 5 decreased after the inhalation of bronchodilators (4.48±1.57 s vs 2.08±1.15 s, for placebo and bronchodilators, respectively; p<0.001). There was a close association between reductions in t1/2QT and lower τ-V˙O2/MRT-Δ[HHb] values after active treatment (r=0.88 (95% CI=0.64 to 0.96); p=0.008). Furthermore, the area under the ‘overshoot’ of the Δ[HHb] response, a pattern of abnormality associated with impaired microvascular QO2 12 decreased (n=3; p=0.04) or disappeared (n=4) after bronchodilator treatment compared with placebo.

Figure 1

Comparative effects of inhaled placebo (A) and bronchodilators (B) on changes in deoxyhaemoglobin (Δ[HHb] by near-infrared spectroscopy) in the vastus lateralis in a patient with severe COPD. Note that bronchodilation was associated with slower Δ[HHb] adjustments. These data are consistent with improved microvascular O2 delivery to the peripheral muscles after bronchodilation. TD, time delay; τ, time constant; MRT, mean response time.

Interestingly, the decrease in lung hyperinflation (ie, higher IC) was significantly related to faster QT (figure 2A) and V˙o2 kinetics (r=−0.73 (95% CI −0.34 to −0.91); p=0.03) but slower rates of muscle deoxygenation (MRT-[HHb], r=−0.78 (95% CI −0.40 to −0.92); p=0.03). Therefore, patients with larger increases in IC after the inhalation of bronchodilators had greater decrements in τ-V˙o2/MRT-Δ[HHb] (figure 2B).

Figure 2

Increases in inspiratory capacity (IC) with inhaled bronchodilators (Δ=bronchodilator–placebo) were associated with faster cardiac output (QT) kinetics (A) and larger decrements in τ-Embedded Image/MRT-[HHb]) (B). Collectively, these data suggest that improvement in central cardiovascular adjustments after reduction in the operating lung volumes with inhaled bronchodilators (A) led to faster dynamics of microvascular O2 delivery (B). t1/2, half-time; τ, time constant; MRT, mean response time; Embedded Image, oxygen uptake; [HHb], deoxyhaemoglobin by near-infrared spectroscopy.

Muscle oxygenation and leg effort scores at steady-state exercise

The amplitude of Δ[HHb] after the initial ‘fast’ response was significantly reduced with bronchodilators compared with placebo (185±115 μM/cm vs 153±108 μM/cm; p=0.009), remaining lower throughout exercise (figure 3). Δ[HHb] at isotime was significantly related to larger postbronchodilator increases in IC (r=−0.82 (95% CI −0.62 to −0.94); p=0.01) and improvements in Tlim (r=−0.80 (95% CI=−0.54 to −0.92); p=0.01). Leg effort scores at isotime were also lower after active treatment compared with placebo (table 2) and they were marginally related to changes in Tlim (r=−0.53; p=0.08).

Figure 3

Effects of inhaled bronchodilators (filled symbols) and placebo (open symbols) on skeletal muscle deoxyhaemoglobin (Δ[HHb], % of maximal value obtained on a previous leg occlusion) in the vastus lateralis at standardised time points during high-intensity exercise. Isotime, the shortest exercise test between the two interventions on a given subject. Values are mean and standard error.

Discussion

This seems to constitute the first study to investigate the effects of inhaled bronchodilators on the determinants of (estimated) microvascular O2 delivery (QO2) and utilisation (phase II pulmonary V˙o2 kinetics) during exercise in patients with COPD. The main novel findings of the present study are that treatment with inhaled bronchodilators to reduce airflow obstruction and lung hyperinflation were associated with: (1) faster on-exercise dynamics of QT and V˙o2; (2) slower kinetics and decreased amplitude of muscle deoxygenation, estimated by Δ[HHb]; and (3) lower leg effort scores. The relationship between lower operating lung volumes after bronchodilators (ie, larger IC) with faster ‘central’ haemodynamic adjustments to exercise (QT) and improved muscle oxygenation (τ-V˙o2/MRT-Δ[HHb]) suggest a mechanistic link between dynamic hyperinflation and impaired QO2 in these patients. Our data, therefore, provide novel evidence that by alleviating the mechanical burden of breathing in patients with COPD pharmacological bronchodilation promotes an enhancement of O2 delivery to the working peripheral muscles.

Interpreting the kinetics of muscle oxygenation during exercise

NIRS is a non-invasive technique that allows continuous measurement of key determinants of peripheral muscle oxygenation during exercise. The time course of Δ[HHb] reflects the local balance between QO2 and V˙o2 in the area under investigation.11 12 Previous studies have shown that Δ[HHb] provides a surrogate of fractional O2 extraction,22 23 and its measurement in parallel with (muscle) phase II of the kinetics of pulmonary V˙o2 can be used to make useful inferences regarding microvascular QO2 via the Fick principle.4 5 11 12 22 23 Therefore, the muscle haemodynamic response at the onset of exercise is characterised by an initial ‘phase I’ response centrally related to the mechanical effect of muscle contraction (ie, muscle pump) and a slower ‘phase II’ matched with metabolic demand. Conversely, muscle (V˙o2) response seems to be well characterised by a single monoexponential function from the start of exercise.24 Therefore, the resulting response of fractional O2 extraction (∼[HHb]) encompasses an early delay where O2 delivery meets (or exceeds) demand (TD) and a rapidly increasing response when the dynamics of capillary flow blood are slower than the rate of change in V˙o2.11 12 Consequently, the overall [HHb] kinetics are expected to be inversely related to the dynamics of QO2.4 5 11 12 22 23

In some specific circumstances, this pattern of response is further complicated by a transient ‘overshoot’ in Δ[HHb],25–27 indicating an increase in fractional O2 extraction above the steady-state level due to impaired phase II blood flow relative to muscle V˙o2.12 Consistent with these assumptions, previous studies found an ‘overshoot’ in Δ[HHb] in patients with type 2 diabetes,25 chronic heart failure19 and COPD.4 The effects of bronchodilators on these complex inter-relationships, however, have not been previously investigated in patients with COPD.

Effects of bronchodilators on leg O2 delivery

The present study provides novel evidence that alleviation of airflow obstruction after pharmacological bronchodilation, with consequent reductions in lung hyperinflation (table 2), resulted in slower kinetics and lower amplitude of muscle deoxygenation (table 3 and figure 3, respectively). In addition, bronchodilators lessened or eliminated the ‘overshoots’ in Δ[HHb], further indicating that muscle oxygenation improved after active treatment. Although the precise mechanism(s) leading to enhanced leg oxygenation could not be determined in this non-invasive study, it is noteworthy that active treatment accelerated the kinetics of the ‘central’ haemodynamic adjustments to exercise. In this context, a decrease in the operating lung volumes might have been associated with increased reliance of parasympathetic control of HR,26 an effect that could be related to less airways compression27 and lower pulmonary arterial pressure.28 Reduction in air trapping and expiratory positive pleural pressure may have increased venous return, decreased pulmonary vascular resistance and lessened the compression of the left ventricle by the overfilled right ventricle.29 As there was a close association between faster QT kinetics and reductions in τ-(V˙o2)/MRT-Δ[HHb], an index of impaired microvascular QO2 4 5 it is likely that the slower changes in muscle deoxygenation were driven by quicker ‘central’ haemodynamic responses at the onset of exercise. Consistent with this notion, Laveneziana and colleagues recently found significant correlations between faster HR dynamics and improved V˙o2 kinetics after inhaled bronchodilators in patients with moderate COPD (FEV1=57±19% predicted).30

Other potential mechanisms for the faster on-exercise dynamics of QT and improved muscle oxygenation are related to the pharmacological properties of the bronchodilators. Although the limited available evidence suggests that inhaled bronchodilators have few haemodynamic effects during exercise in patients with COPD,6 28 their role in modulating skeletal muscle blood flow at the exercise transient in animals is still controversial.31 Additional studies, therefore, are warranted to address this specific issue in patients with COPD.

Interestingly, Δ[HHb] remained lower during steady-state exercise after the inhalation of bronchodilators (figure 3) despite no significant increases in QT or Spo2 (table 2). This finding raises the intriguing hypothesis that alleviation of the respiratory muscle work with bronchodilators may have redirected a fraction of the QT from these muscles back to the appendicular muscles. In fact, we previously reported that respiratory muscle unloading by proportional assisted ventilation increased microvascular QO2 at the same QT, suggesting improved perfusion due to blood flow redistribution from respiratory to working peripheral muscles.18 Further studies are needed to apportion the role of lower operating lung volumes and improved muscle oxygenation in enhancing exercise tolerance in patients with COPD.

Study limitations

Constant work rate exercise tests at 70–80% peak have been recognised as the most sensitive testing modality to unravel the functional benefits of bronchodilators.32 However, our findings may not be relevant to milder exercise intensities where lower ventilatory requirements may decrease, or even abolish, the beneficial effects of bronchodilators on muscle deoxygenation. We also recognise that the use of cycling instead of walking may have contributed to magnify the positive effects of bronchodilation on exercise tolerance, as contractile fatigue33 and the sense of leg effort34 are more relevant contributors to decrease exercise capacity in the former modality. It should also be acknowledged that our non-hypercapnic patients with advanced COPD were severely hyperinflated and they presented with substantial reductions in DLCO (table 1)—that is, a disease profile more compatible with the emphysematous subtype. Consequently, our results should be cautiously extrapolated to other COPD phenotypes. Finally, most of the correlation coefficients had wide confidence intervals, a likely consequence of the small sample size. Therefore, it is advisable that these specific results be confirmed in larger studies.

Conclusions

The present study showed that administration of inhaled short-acting bronchodilators to reduce lung hyperinflation slowed the dynamics and reduced the amplitude of muscle deoxygenation during high-intensity cycling exercise in patients with moderate to severe COPD. Considering that these findings were associated with faster on-exercise kinetics of QT and the ‘metabolic’ phase of pulmonary V˙o2, our data are consistent with a speeding and heightening of O2 delivery to skeletal muscles of patients with COPD treated with inhaled bronchodilators.

Acknowledgments

The authors would like to thank all colleagues from the Pulmonary Function and Clinical Exercise Physiology Unit (Division of Respiratory Diseases, Federal University of Sao Paulo (UNIFESP), Brazil) for their friendly collaboration. More importantly, however, they are indebted to the patients for their effort and enthusiastic cooperation throughout the study.

References

Supplementary materials

  • Web Only Data thx.2009.120857

    Files in this Data Supplement:

Footnotes

  • Linked articles 135129

  • Funding DCB was suppported by a Doctoral Research Fellowship from Coordenadoria de Aperfeiçoamento do Pessoal de Nível Superior (CAPES), Brazil. GRC was suppported by a Postdoctoral Research Fellowship from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil. JAN is an Investigator of the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Brazil. The study was supported by a Research Grant from FAPESP, Brazil.

  • Competing interests None.

  • Ethics approval This study was conducted with the approval of the Institutional Medical Ethics Committee.

  • Provenance and peer review Not commissioned; externally peer reviewed.

Linked Articles