Abstract
Our understanding of the mechanisms of dyspnoea in fibrotic interstitial lung disease (ILD) is incomplete. The aims of this study were two-fold: 1) to determine whether dyspnoea intensity is better predicted by neural respiratory drive (NRD) or neuromechanical uncoupling (NMU) of the respiratory system in fibrotic ILD, and 2) to examine the effect of breathing 60% oxygen on NRD, NMU and dyspnoea ratings.
Fourteen patients with fibrotic ILD were included. Visit 1 comprised a familiarisation incremental cycle exercise test, Visit 2 comprised a normoxic incremental cycling test to address Aim 1, and Visits 3 and 4 consisted of constant-load cycling while breathing room air or 60% oxygen to address Aim 2. Diaphragmatic electromyography (EMGdi) was used as a surrogate of NRD. NMU was calculated as the ratio between EMGdi (%max) and tidal volume (%vital capacity).
On adjusted analysis, NMU and its constituents were all significantly associated with dyspnoea ratings during incremental cycling, with EMGdi having the strongest correlation. The between-treatment change in dyspnoea ratings during constant load cycling was only correlated with change in exercise endurance time and NMU.
Dyspnoea more strongly reflected the level of EMGdi than NMU in fibrotic ILD. However, the improvement in dyspnoea with 60% oxygen was better predicted by improvements in NMU.
Abstract
Neural respiratory drive is a strong independent predictor of dyspnoea in patients with fibrotic ILD http://ow.ly/MXjd30hcabH
Introduction
Dyspnoea is a common and debilitating symptom in patients with fibrotic interstitial lung disease (ILD). In this population, dyspnoea is an independent predictor of morbidity and mortality, and is associated with decreased health-related quality of life and reductions in functional capacity [1, 2]. Despite its clinical importance, our understanding of the mechanisms of dyspnoea in fibrotic ILD remains incomplete.
The mechanisms of dyspnoea, in both health and disease, are complex and multifactorial [3]. However, accumulating evidence suggests that neural respiratory drive (NRD) and the imbalance between NRD and the mechanical response of the respiratory system (i.e. neuromechanical uncoupling (NMU)) are key contributors [4]. A strong relationship has recently been established between NRD, estimated using crural diaphragm electromyography (EMGdi), and dyspnoea intensity ratings in patients with ILD [5]. Furthermore, a study in chronic obstructive pulmonary disease (COPD) showed that dyspnoea intensity was strongly associated with both NRD and NMU; however, the correlation between dyspnoea intensity and NRD was more significant [6]. It is unknown whether these findings are also translatable to patients with ILD.
Treatment options to reduce dyspnoea in ILD are limited; however, recent evidence suggests that breathing a hyperoxic inspirate (inspiratory oxygen fraction (FIO2) =0.60) significantly reduces dyspnoea intensity ratings and can eliminate or delay the onset of ‘unsatisfied inspiration’ during exercise [7]. To our knowledge, no study has investigated the effect of supplemental oxygen on NRD and/or NMU in ILD in order to determine which outcome is more strongly associated with reduced exertional dyspnoea ratings. Accordingly, the aims of this study were to determine whether dyspnoea intensity in fibrotic ILD is better predicted by NRD or NMU (Aim 1), and to examine the effect of 60% oxygen on EMGdi-derived measures of NRD and NMU (Aim 2).
Methods
Participants
This study included 14 fibrotic ILD patients with isolated lung involvement that simultaneously participated in a previously published study (n=20) examining the effects of hyperoxia on cycle endurance time and dyspnoea [7]. All patients had resting arterial oxygen saturation ≥92%. Patients who had recently completed pulmonary rehabilitation (≤6 weeks prior to enrolment), or with significant emphysema, pulmonary hypertension, extrapulmonary disease and/or other contraindications to exercise were excluded from participating.
Study design
This study was a component of a larger randomised, single-blind, placebo-controlled cross-over study (ClinicalTrials.gov: NCT01781793) with institutional ethical approval (H13-00059). Written informed consent was obtained from all patients prior to enrolment. Patients visited the laboratory on four occasions separated by ≥48 h between visits. Medical history, chronic activity-related dyspnoea questionnaires, pulmonary function testing, and a symptom-limited incremental cycle exercise test with dyspnoea assessment for familiarisation purposes were completed at Visit 1. During Visit 2, patients performed pulmonary function testing and a symptom-limited incremental cycle exercise test while equipped with an EMGdi catheter. Visit 2 data were used to address the primary aim. During Visits 3 and 4, patients performed pulmonary function testing and a symptom-limited constant-load cycle exercise test while breathing room air (FiO2=0.21) or hyperoxia (FiO2=0.60), in randomised order. Gases were delivered into a non-diffusing Douglas bag connected to a two-way non-rebreathing valve via large-bore tubing (Hans Rudolph, Inc., Shawnee, KS, USA). The gas delivery system was identical for both conditions, and patients were blinded to the inspired gas condition. Data from Visits 3 and 4 were used to address the secondary aim of the study.
Pulmonary function
Patients performed spirometry, plethysmography, 12-s maximal voluntary ventilation, and single-breath diffusing capacity of the lungs for carbon monoxide (DLCO), according to established guidelines [8–10]. Values were measured using a commercially available system (Vmax Encore 229, V62J Autobox; CareFusion, Yorba Linda, CA, USA) and expressed as a percentage of predicted values [11, 12].
Exercise testing protocol
Exercise was conducted using an electronically braked cycle ergometer (Ergoselect 200P; Ergoline, Bitz, Germany) at a self-selected cadence ≥60 revolutions per minute. All tests started with 6 min of rest followed by a 1-min warm-up of unloaded pedalling. Incremental tests proceeded with 15-W stepwise increases in work-rate, starting at 15 W, every 2 min until symptom limitation. Peak work-rate was defined as the highest work-rate sustained for at least 30 s. Constant-load exercise tests proceeded with an immediate increase in work-rate to 75% of the peak incremental work-rate, which was sustained until symptom limitation. Iso-time was defined as the highest equivalent submaximal time achieved on both constant-load exercise tests performed by a given patient under both conditions. No verbal encouragement was provided during constant-load exercise tests.
Symptom evaluation
Patients rated the intensity of ‘breathing discomfort’ (dyspnoea) and ‘leg discomfort’ at rest, every 2 min during exercise, and at peak exercise using the modified Borg 0–10 category-ratio scale [13]. The scale's end-points were anchored such that 0 represented ‘no breathing/leg discomfort’ and 10 represented ‘the most severe breathing/leg discomfort ever experienced or imagined’.
Cardiorespiratory responses to exercise
Standard cardiorespiratory responses were measured on a breath-by-breath basis at rest and during exercise (Vmax Encore 229; CareFusion). This metabolic cart has not been validated with an FIO2 of 0.60; therefore, differences in gas exchange parameters (i.e. oxygen consumption and carbon dioxide production) and gas exchange-derived parameters (i.e. ventilatory equivalents and the respiratory exchange ratio) were not compared between conditions (Visits 3 and 4). Flow was calibrated using a 3-L calibration syringe connected to the breathing apparatus and using the appropriate gas mixture from the Douglas bag at a wide range of flow rates. The inflection point in tidal volume relative to minute ventilation (VT/V′E) was determined for each patient by examining 30-s averaged data throughout incremental exercise (Visit 2). The time at which the corresponding VT was achieved during Visits 3 and 4, if at all, was also identified using the 30-s averaged data from those visits. Arterial oxygen saturation was estimated via finger pulse oximetry (Radical-7 Pulse CO Oximeter; Masimo Corp., Irvine, CA, USA). Heart rate was recorded using 12-lead electrocardiography.
Diaphragm electromyography
Crural diaphragm electromyography was measured using a multi-pair oesophageal electrode catheter (Guangzhou Respiratory Institute Medical Equipment and Scientific Co. Ltd, Guangzhou, China) and was used as a surrogate of NRD [14]. Diaphragm electromyography was collected, processed and analysed as previously described [15]. Briefly, the catheter consists of five consecutive recording pairs, and this was passed through the nose and positioned based on the amplitude of EMGdi signals obtained simultaneously from these recording pairs during tidal breathing. To minimise discomfort, a lidocaine hydrochloride non-aerosol spray (Lidodan; Odan Laboratories, Montreal, Canada) was administered in the nasal and pharyngeal passages. Raw EMGdi signals were sampled at 2 kHz using a PowerLab 16/35 analogue-to-digital converter (ADInstruments, Castle Hill, Australia), processed with a 60-Hz notch filter (Model 08-GL biological Amplifier, Guangzhou Respiratory Institute Medical Equipment and Scientific Co. Ltd) as well as a band-pass filter between 20 Hz and 1 kHz using LabChart Pro Version 7.3.7 software (ADInstruments), and converted to a root mean square (RMS) using a time constant of 0.1 s. The maximum RMS value for each inspired breath was manually selected between QRS complexes to avoid the influence of cardiac artefact on the EMG signal. The electrode pair with the largest amplitude for a given breath was used for the analysis. Data were reported as a percentage of maximal EMGdi, which was defined as the highest RMS value recorded during an inspiratory capacity manoeuvre at rest or during exercise. We used inspiratory capacity manoeuvres exclusively instead of using multiple maximal inspiratory manoeuvres (e.g. sniffs and maximal inspiratory pressure manoeuvres) to obtain maximal EMGdi. This was done to standardise the normalisation procedure across all subjects, and because our previous work shows that inspiratory capacity manoeuvres always result in the greatest EMGdi activity [16]. NMU was defined as the ratio of EMGdi expressed as a percentage of maximum to VT expressed as a percentage of vital capacity (VC) (i.e. EMGdi (%max):VT (%VC)).
Statistical analyses
Exercise variables were measured continuously, with mean values determined over a 30-s epoch linked to the time of dyspnoea assessment. We determined the association of dyspnoea intensity ratings with NMU and its individual components (i.e. EMGdi (%max) and VT (%VC)), using an unadjusted repeated-measures linear regression. All of these predictors were then applied to a multivariate model, also using repeated-measures analysis, to determine independent predictors. Variables included in the multivariate model were assessed for multicollinearity. Paired t-tests were used to compare outcomes between the room air and hyperoxic conditions at iso-time and peak exercise. Spearman's correlation coefficients were used to make intra-subject, between-treatment comparisons of dyspnoea intensity ratings with select physiological variables during constant-load exercise at iso-time. Data are presented as mean±sd or median (interquartile range) where appropriate. Statistical significance was set at p<0.05. Data were analysed using Stata 11.2 software (StataCorp LP, College Station, TX, USA).
Results
Participants
Patient characteristics and peak exercise data are summarised in table 1. On average, patients had moderate physiological impairment, as demonstrated by reduced forced vital capacity and diffusing capacity of the lungs for carbon monoxide, and reduced functional capacity, as evidenced by reduced peak aerobic capacity and work-rate relative to predicted values [11, 12, 17]. The VT/V′E inflection point occurred at 76±28 W during incremental exercise. None of the patients qualified for supplemental oxygen, at rest or with ambulation, according to regional criteria [18]. No cardiac abnormalities were noted during exercise testing.
Patient characteristics and maximal exercise data
Associations with dyspnoea
During the incremental cycle exercise test on Visit 2, 97 observations were obtained. One patient was excluded from this analysis (Aim 1) owing to lack of EMGdi data on Visit 2; the patient refused the catheter. Dyspnoea intensity ratings were significantly associated with NMU, EMGdi (%max) and VT (%VC) (all p<0.001, fig. 1). With adjusted analysis, EMGdi (%max), NMU and VT (%VC) remained significantly associated with dyspnoea intensity ratings (p<0.001, p=0.006 and p=0.02, respectively), with EMGdi having the strongest correlation (fig. 1).
Trend lines for individual patients demonstrating the relationship between dyspnoea intensity and a) electromyography of the diaphragm (EMGdi), b) neuromechanical uncoupling, and c) tidal volume (VT) expressed as a percentage of vital capacity (VC). Reported r-values are expressed as median (interquartile range) from Spearman's correlations within each patient.
Effects of hyperoxia
The sensory and physiological effects of hyperoxia during constant-load cycle exercise are summarised in table 2. Three patients were excluded from this analysis (Aim 2); these patients refused the catheter for at least one of these visits and/or did not complete both visits. The average time between constant-load cycle exercise tests was 6±2 days. The mean work-rate for the constant-load cycle exercise tests was 71±27 W, and the median iso-time was 6.0 (4.0–10.0) min. Constant-load cycle exercise endurance time increased significantly with hyperoxia compared to that with room air (24.2±11.1 versus 9.6±8.0 min, respectively; p<0.001).
Selected sensory and physiological parameters at rest, iso-time and peak exercise for constant-load cycle exercise tests with room air and hyperoxia
Dyspnoea intensity ratings were similar between room air and hyperoxia at peak exercise (p=0.38), but significantly reduced at iso-time with hyperoxia versus room air (1.8±1.9 versus 3.8±3.0 Borg 0–10 units, respectively; p=0.005; table 2). EMGdi was significantly reduced during constant-load cycle exercise at iso-time (p<0.001) and peak exercise (p=0.005) with hyperoxia versus room air (table 2). NMU was significantly improved at iso-time with hyperoxia versus room air (p=0.003; table 2). The between-treatment change in iso-time dyspnoea ratings during constant-load cycle exercise was only correlated with the change in exercise endurance time (r=−0.67, p=0.02) and NMU (r=0.63, p=0.04). Additional Spearman's correlation coefficients are shown in table 3.
Spearman's correlations for selected parameters with the change in dyspnoea intensity ratings at iso-time with hyperoxia versus room air
During constant-load cycle exercise with room air, all 11 patients included in the analysis for Aim 2 achieved or surpassed a VT that corresponded to the VT/V′E inflection point identified on the incremental test. Two of these patients did not achieve this VT with hyperoxia, and the time at which this VT was achieved was delayed with hyperoxia versus room air in the remaining nine patients (5.5±5.1 versus 1.9±0.7 min, respectively; p=0.06). Measures of dyspnoea, EMGdi (%max) and NMU were not different at this measurement time between treatments (p=0.12, p=0.12 and p=0.09, respectively).
Discussion
We sought to determine whether dyspnoea intensity in fibrotic ILD is primarily determined by NRD or NMU, and to examine the impact of 60% oxygen on NRD and NMU. Our main findings were as follows: 1) dyspnoea intensity ratings were significantly correlated with both EMGdi and NMU during incremental cycle exercise, but with a stronger correlation between dyspnoea and EMGdi; and 2) the improvements in dyspnoea intensity ratings at iso-time during constant-load cycle exercise with hyperoxia were significantly correlated with the change in NMU, but not with the change in EMGdi.
It has been proposed that increased dyspnoea is a reflection of increased levels of NRD in both health and disease [4]. There is also evidence that central and peripheral afferent feedback pathways may neuromodulate the intensity of perceived dyspnoea through direct effects on NRD [4]. Accordingly, increased dyspnoea intensity ratings in patients with cardiopulmonary disease likely correspond to the awareness of an increased NRD needed to meet an increased ventilatory demand in the context of increased mechanical loading and/or reduced capacity of the respiratory muscles. For example, patients with ILD demonstrate reduced compliance and increased static recoil pressure of the lung due to fibrosis. The static pressure-volume curve of the lungs is therefore shifted downwards and to the right [19, 20]. The pressure-volume relationship of the entire respiratory system is therefore contracted along its volume axis [21], whereby resting inspiratory capacity and inspiratory reserve volume are generally diminished. This reduction in baseline reserve translates to less potential for VT expansion during exercise [22, 23]. The ILD patients in the present study achieved critical constraints on VT expansion at markedly lower absolute work-rates than what has been previously documented in healthy men and women of similar ages [24]. Thus, ILD patients are likely to have increased elastic loading of the respiratory muscles [21], which necessitates greater respiratory muscle force for inspiration [25], and therefore an increased NRD [5], in order to meet ventilatory requirements.
Faisal et al. [5] showed that the relationship between EMGdi and dyspnoea intensity is preserved in patients with ILD and COPD, as well as healthy controls. This supports the hypothesis that NRD is an important component of the physiological basis of increased dyspnoea intensity, as disease-specific differences in respiratory function and/or capacity, and therefore afferent feedback pathways, do not influence this association [4]. Indeed, another recent study in COPD patients [6] examined the relationships of dyspnoea intensity with NRD and NMU, concluding that dyspnoea intensity is closely related to EMGdi and NMU, but that the relationship between dyspnoea intensity ratings and EMGdi was comparatively stronger. This is consistent with our findings, suggesting that increased NRD largely contributes to increased dyspnoea intensity in patients with fibrotic ILD, at least during normoxic incremental exercise.
We recently demonstrated significant reductions in iso-time dyspnoea intensity ratings during constant-load exercise with hyperoxia versus room air [7]. However, no signification correlations were found between the reduction in dyspnoea intensity ratings and the between-treatment changes in standard ventilatory measures during exercise. We reasoned that a reduction in dyspnoea intensity ratings with hyperoxia may therefore be associated with improvements in neurophysiological parameters such as NRD and/or NMU. However, despite EMGdi being more significantly correlated with dyspnoea intensity ratings during normoxic incremental exercise, the improvement in dyspnoea intensity ratings with hyperoxia were only significantly correlated with an improvement in NMU and exercise performance (i.e. endurance time). This suggests that hyperoxia reduces dyspnoea via multiple mechanisms, with a combined beneficial effect on both NRD and respiratory mechanics.
Previous studies have shown that hyperoxia improves oxygen delivery to the skeletal muscles by widening the arteriovenous oxygen content difference, thereby increasing the anaerobic threshold and delaying the onset of metabolic acidosis and skeletal muscle fatigue [26, 27]. The subsequent reduction in V′E at iso-time and peak exercise in the present study may have reflected, at least in part, both a reduction in EMGdi as well as VT, with attendant improvements in perceived breathlessness. This is further supported by the delayed onset of the VT/V′E inflection with hyperoxia versus room air. NMU was defined in this study as the ratio of EMGdi (%max) to VT (%VC), and was therefore a direct representation of the interaction between EMGdi and VT. In this context, had only EMGdi been reduced and VT held constant with hyperoxia, NMU would have decreased (i.e. improved). Alternatively, had only VT been reduced and EMGdi held constant with hyperoxia, NMU would have increased (i.e. worsened). We observed reductions in both EMGdi and VT with hyperoxia versus room air, with a relatively larger reduction in EMGdi (26±11 versus 12±6%, respectively). Accordingly, our study suggests that hyperoxia reduces dyspnoea mainly via an improvement in NMU, which reflects an optimised relationship between EMGdi and VT. The interaction between the improvement in NRD and the change in VT seems to have a more important effect than the improvement in NRD alone.
Limitations
There are some limitations to this study that must be acknowledged. First, we used EMGdi as a surrogate of NRD, and the limitations of this method have been well described [5]. Second, while the diaphragm is the primary inspiratory muscle, we cannot ignore the potential contribution of extra-diaphragmatic inspiratory muscles, such as the sternocleidomastoid, scalene and intercostal muscles. Changes in the ventilatory response with hyperoxia could elicit changes in respiratory muscle recruitment patterns and their relative contributions, which could influence measures of NMU. Third, owing to safety reasons, it was not feasible to blind all study personnel to the study condition. However, we do not believe this impacted our results as all test procedures were rigorously standardised, including an absence of verbal encouragement for all constant-load exercise tests. Fourth, we acknowledge that this invasive study had a small sample size, and therefore may not be representative of all patients with fibrotic ILD. Finally, we included a variety of ILD subtypes; however, all patients had isolated lung involvement. Owing to the similar physiology of these diseases, we expect the physiological findings to be similar across these subtypes.
Conclusion
Our data suggest that the intensity of dyspnoea more likely reflects increased NRD, as estimated by EMGdi, rather than NMU, in patients with ILD during normoxic incremental exercise. However, attenuation of dyspnoea intensity ratings with hyperoxia was correlated with improvements in NMU, but not EMGdi. This suggests that the combined benefits of hyperoxia on the coupling between NRD and the mechanical output of the respiratory system are more important in relieving dyspnoea than the effects on NRD alone. Future research is needed to better isolate the independent effects of EMGdi and the mechanical output of the respiratory system (i.e. VT) on dyspnoea during exercise in patients with ILD. This knowledge may aid in the development of symptom-based treatment options in this population.
Disclosures
Supplementary Material
C.J. Ryerson ERJ-01726-2017_Ryerson
P.G. Camp ERJ-01726-2017_Camp
J.A. Guenette ERJ-01726-2017_Guenette
Acknowledgements
The authors would like to thank all of the patients for their enthusiastic participation in this study. All authors played a role in the content and writing of the manuscript. M.R. Schaeffer, C.J. Ryerson, A.W. Sheel, N. Khalil, P.G. Camp and J.A. Guenette had input into the study design and conduct of the study; M.R. Schaeffer, C.J. Ryerson, A.H. Ramsook, Y. Molgat-Seon, S.S. Wilkie, R.A. Mitchell and S.S. Dhillon collected the data; and M.R. Schaeffer, C.J. Ryerson, A.H. Ramsook, S.S. Dhillon, R.A. Mitchell and J.A. Guenette performed data analysis.
Footnotes
Support statement: This study was funded by the British Columbia Lung Association (BCLA), an Emerging Research Leaders Initiative Grant from the Canadian Respiratory Research Network, and seed funding from the Dept of Physical Therapy at the University of British Columbia (Vancouver, BC, Canada). M.R. Schaeffer was supported by fellowships from the University of British Columbia and BCLA. C.J. Ryerson and P.G. Camp were supported by Scholar Awards from the Michael Smith Foundation for Health Research (MSFHR). A.H. Ramsook was supported by the University of British Columbia 4-Year Fellowship (UBC 4YF) and an Alexander Graham Bell Canada Graduate Scholarship from the Natural Sciences and Engineering Research Council of Canada (NSERC). Y. Molgat-Seon was supported by the UBC 4YF and a Postgraduate Scholarship from the NSERC. J.A. Guenette was supported by a Scholar Award from the MSFHR, a Clinical Rehabilitation New Investigator Award from the Canadian Institutes of Health Research, and a New Investigator Award from the Providence Health Care Research Institute and St Paul's Hospital Foundation. The funders had no role in the study design, data collection and analysis, or preparation of the manuscript. Funding information for this article has been deposited with the Crossref Funder Registry.
Conflict of interest: Disclosures can be found alongside this article at erj.ersjournals.com
- Received August 23, 2017.
- Accepted October 14, 2017.
- Copyright ©ERS 2018