key: cord-0055173-uvro49o4
authors: Huang, Yu-Chen; Lin, Ting-Yu; Wu, Hau-Tieng; Chang, Po-Jui; Lo, Chun-Yu; Wang, Tsai-Yu; Kuo, Chih-Hsi Scott; Lin, Shu-Min; Chung, Fu-Tsai; Lin, Horng-Chyuan; Hsieh, Meng-Heng; Lo, Yu-Lun
title: Cardiorespiratory coupling is associated with exercise capacity in patients with chronic obstructive pulmonary disease
date: 2021-01-12
journal: BMC Pulm Med
DOI: 10.1186/s12890-021-01400-1
sha: de2ad97194c7f63a1d5b35e635f24b6f75338ff7
doc_id: 55173
cord_uid: uvro49o4

BACKGROUND: The interaction between the pulmonary function and cardiovascular mechanics is a crucial issue, particularly when treating patients with chronic obstructive pulmonary disease (COPD). Synchrogram index is a new parameter that can quantify this interaction and has the potential to apply in COPD patients. Our objective in this study was to characterize cardiorespiratory interactions in terms of cardiorespiratory coupling (CRC) using the synchrogram index of the heart rate and respiratory flow signals in patients with chronic obstructive pulmonary disease. METHODS: This is a cross-sectional and preliminary data from a prospective study, which examines 55 COPD patients. K-means clustering analysis was applied to cluster COPD patients based on the synchrogram index. Linear regression and multivariable regression analysis were used to determine the correlation between the synchrogram index and the exercise capacity assessed by a six-minute walking test (6MWT). RESULTS: The 55 COPD patients were separated into a synchronized group (median 0.89 (0.64–0.97), n = 43) and a desynchronized group (median 0.23 (0.02–0.51), n = 12) based on K-means clustering analysis. Synchrogram index was correlated significantly with six minutes walking distance (r = 0.42, p = 0.001) and distance saturation product (r = 0.41, p = 0.001) assessed by 6MWT, and still was an independent variable by multivariable regression analysis. CONCLUSION: This is the first result studying the heart–lung interaction in terms of cardiorespiratory coupling in COPD patients by the synchrogram index, and COPD patients are clustered into synchronized and desynchronized groups. Cardiorespiratory coupling is associated with exercise capacity in patients with COPD.

In patients with chronic obstructive pulmonary disease (COPD), cardiovascular disease is a prevalent comorbidity and leading cause of death. The prevalence of cardiovascular disease within this population can be attributed to shared risk factors (e.g., cigarette smoking and exposure to noxious gases) as well as oxidative stress and reduced physical activity related to COPD [1] [2] [3] [4] . Deleterious pulmonary function (e.g., dynamic hyperinflation and hypoxia) can also impair cardiovascular mechanics in patients with COPD [5] [6] [7] . A more comprehensive understanding of the heart-lung interaction could be highly beneficial in efforts to treat patients with COPD [8] [9] [10] .

Heart-lung interactions can be classified according to the underlying related but different mechanisms: (1) respiratory sinus arrhythmia, (2) cardioventilatory coupling, and (3) respiratory stroke volume synchronization [11] . For example, during inspiration, central inspiratory drive [12] and negative intrathoracic pressure [13] both contribute to an increase in heart rate [14] . Negative intrathoracic pressure promotes filling of the right ventricle and impedes filling of the left ventricle [15] . A decrease in arterial blood pressure tends to increase respiratory rate and tidal volume through the baroreflex [16] .

Cardiorespiratory coupling (CRC) is an intuitive method to depict and quantify the complicated heartlung interaction by calculating the phases ratio between heartbeat and respiration. CRC in healthy subjects has been extensively studied in terms of age-related evolution [17] and its association with sleep stage transitions [18] . CRC grows in strength in the first 180 day after birth [17] , and continues to evolve with age. Association of ages and CRC synchronization is different during sleep and rest periods. In healthy adults, CRC synchronization during rest periods is not correlated with age [19, 20] ; however, the strength of CRC has been shown to decrease in elderly adults during sleep [21] . CRC has also been linked to obstructive sleep apnea (OSA) [22] and it has been suggested as a tool by which to assess OSA severity. Under the clinical observation of intimate heart-lung interaction in patients with COPD [23] , we hypothesized that CRC could provide clinically relevant information. The aim of this study was to apply the synchrogram index to evaluate CRC in patients with COPD, and to cluster patients based on their synchrogram indices.

This observational cross-sectional study was based on data obtained from a preliminary prospective study, conducted at Chang Gung Memorial Hospital (CGMH) in Linkou, Taiwan. The sample included patients who underwent regular follow-up as out-patients at CGMH between January 2019 and January 2020. Inclusion criteria included clinical diagnosis of COPD based on the Global Initiative for Obstructive Lung Disease Criteria (GOLD) [24] and post-bronchodilator FEV1 < 80% of the predicted normal value via diagnostic spirometry. Other inclusion criteria were age ≥ 40 years without known heart disease. Exclusion criteria included patients with HFlowEF [heart failure with low ejection fraction (< 40%)], known malignancy, or atrial fibrillation as well as those using oxygen or anti-arrhythmic agents for arrhythmia. All COPD patients underwent cardiac echo analysis, biochemical analysis [eosinophils, high sensitivity C-reactive protein (HS-CRP), and IgE], pulmonary function tests, chest high-resolution computed tomography (HRCT) scanning, a six-minute walking test (6MWT) and a coupling test during the first visit of enrollment. Emphysema was defined based on chest HRCT report from the radiologist and one pulmonologist [25] . Clinical profiles, a list of inhalation medicines, antipsychotic agents, result of emphysema based on HRCT and acute exacerbation history [26] were also recorded. Although there is one patient who presented high ratio of FEV 1 /FVC before exercise (0.72) and after exercise (0.73), he was not excluded as the spirometry fulfilled the GOLD guideline when he was diagnosed COPD. The remaining 55 patients with COPD [69 (51-84) years old, 54 male] (Fig. 1 ). All participants signed informed consent prior to enrollment. The study was approved by the Ethics Committee of CGMH (201702150B0).

The 6MWT was carried out on a straight aisle with smooth surface and of 30 m long. Before the exam, the patients rested in a sitting position, during which spirometry was performed to assess pulmonary function, including the flow volume and tidal volume, based on established guidelines [27] . Oxygen saturation, heart rate, arterial blood pressure, and Borg scale values were recorded to assess the degree of dyspnea. Following this preliminary assessment, the patients were instructed to walk as far as possible within a period lasting six minutes. The patients were permitted to stop and rest when they felt tired or dyspneic, and then encouraged to start again as soon as possible. The instructors avoided walking with the subjects, but rather stood within a designated area to provide encouragement with an even tone at intervals of one minute and at fifteen seconds before the end of the exam, in accordance with to the American Thoracic Society (ATS) guideline [28, 29] . Oxygen saturation and heart rate were recorded in real time while walking. At the end Fig. 1 Flow chart of the exam, walking distance, oxygen saturation, distance saturation product (i.e., the product of nadir saturation during exercise and walking), heart rate and Borg scale were recorded, and patients performed spirometry again after exercise.

Experiments were performed in a quiet room with the temperature maintained at 22-24 °C with all necessary equipment prepared beforehand, including ECG leads, pulse oximeter, a breathing tube (a disposable mouthpiece connecting with end-tidal CO 2 sensor and flow sensor), and three ACTiwave devices (CamNtech Ltd, Cambridge, UK) connected to a computer running Lab-Chart8 software. Participants were instructed to avoid inhaling short acting bronchodilators for 4 h, taking oral medicines such as beta-2 agonists, xanthene derivatives for 12 h, and consuming alcohol or caffeine-contained drink for at least eight hours prior to the test. Otherwise, participants could intake other foods before the exam. The chest skin was abraded using gel and then cleaned using alcohol to reduce electrode impedance prior to the attachment of electrocardiogram (ECG) electrodes. Prior to the examination, blood pressure, heart rate, and oxygen saturation were recorded. The subjects wore a pulse oximeter on the index finger and ECG electrodes on the chest wall. A breathing tube was inserted into the subject's mouth with his/her lips sealed around the mouthpiece and a nose clip over the nostrils [30] . Prior to the exam, the subjects were instructed to practice breathing at tidal volume for 1 min and then proceed with the exam when they felt ready. ECG and flow signals were recorded continuously for 5 min using three ACTiwave devices (CamNtech Ltd, Cambridge, UK). The recorded signals were transferred in the European Data Format to Lab-Chart 8 software (ADInstruments, Dunedin, New Zealand), and then exported to text files for analysis.

R peaks were detected using a standard R peak detection algorithm from the ECG signal (Fig. 2a5, b5) . The time differences between consecutive R peaks were calculated, and then converted into an instantaneous heart rate (IHR) time series using the standard interpolation algorithm (Fig. 2a4, b4 ) [31] . The phase of the respiratory signal (denoted as) was extracted using the synchrosqueezing transform (SST) (Fig. 2a2, b2 ) [32] . The phase of IHR (denoted as) was extracted by the same method (Fig. 2a3, b3) . After obtaining the phases of the IHR and the respiratory signal, the synchrogram was used to quantify the cardiorespiratory coupling [33, 34] . The output is the synchrogram index, which is a non-unit quantity between 0 and 1. When the cardiorespiratory coupling is strong, the synchrogram index is close to 1; Otherwise, it is close to 0.

The synchrogram is a signal processing tool used to depict coupling between two oscillatory signals. In the current study, we first obtained the timestamps [ Fig. 2 (dashed line between a2 and a3, b2 and b3)], where the IHR phase attained modulo. We then measured the respiratory phase at as follows: that is, we evaluate the phase of the respiratory signal at [ Fig. 2 (circle points at a2 and b2)] Finally, plot against. When the cardiorespiratory coupling is strong, the phase of the respiratory signal at would be fixed for all, and hence we obtain a horizontal stripe in the plot (Fig. 2a6) ; Otherwise, we obtained scattered points in the plot (Fig. 2b6) . The synchrogram index [34] is aiming to quantify if the plot is scattered or fixed along a horizontal line. It is defined by, where M is the number of detected cycles in the IHR.

All results are presented as median (range) or mean ± standard deviation. The nonparametric exact two-tailed Mann-Whitney U test was used to determine the statistical significance between two groups of continuous variables, and Fisher's exact tests were used for categorical variables. Pearson's correlation coefficient was used to examine the association between 6MWD, distance saturation product (DSP), and other clinical parameters, including the synchrogram index. To avoid variable selection caused by spurious correlations, multivariable regression analysis was based on variables that presented a significant linear relationship (defined as p ≤ 0.05) with 6MWD and DSP. Multivariable regression analysis was performed using the lm function of the R statistical software package. K-means clustering was applied to cluster COPD patients based on their synchrogram indices. Silhouette analysis was performed to select optimal cluster numbers. All reported P values were two-sided, with P < 0.05 considered statically significant. Signals were analyzed using MATLAB. All data were analyzed using R version 3.5.2 (R foundation for statistical computing).

Among 55 COPD patients, 54 (98.2%) were male, 49 (89.1%) had smoking history, 36 (65.5%) were in allergic status, 33 (60%) were confirmed with emphysema from chest HRCT, and only 1 (1.8%) patient fulfilled the criteria of Asthma-COPD overlap (ACO) [35, 36] . The median synchrogram indices in the COPD group was 0.87 and the distribution was skewed (range: 0.02-0.97). The median BMI was 24.7 (range 16.7-32.1), the modified medical research council (mMRC) was 1 (range 0-4), and the COPD assessment test (CAT) was 10 (range 2-29). The median ejection fraction was 65.5% (range 52-90), suggesting that there was no patient with heart failure with mid-range ejection fraction (HFmrEF) in this study. However, there were 24 patients (44.3%) presented diastolic dysfunction. The median left atrial size was 34 mm (range 23-46) and E/e′ ratio (the ratio of the trans mitral early peak velocity over early diastolic mitral annulus velocity) was 8.9 (4.5-20.0). In addition, the median eosinophil count was 129 (range 0-615.6), IgE level was 59.7 (2-1652) and 19 (34.5%) patients had a history of acute exacerbation one year prior to enrollment in the study. Most patients used combination therapy of long-acting 2 agonist (LABA) with long-acting muscarinic antagonist (LAMA) (20 (36.4%)) and triple therapy of LABA with LAMA and ICS (26 (47.3%)) ( Table 1) .

Figures 2a and 2b illustrate CRC analysis based on the synchrogram of IHR and respiratory flow signals. Since there are no definitions of good or poor synchronization, we applied K-means [37] . Two clusters were identified, i.e., synchronized group (n = 43) and desynchronized group (n = 12) according to the optimal cluster number based on the silhouette analysis. The median synchrogram index values in these two groups were as follows: synchronized group (0.89; 0.64-0.97) and desynchronized group (0.23; 0.02-0.51) (Fig. 3b) 

In the 6MWT, patients in the synchronized group achieved longer walking distances than did their counterparts in the desynchronized group (468 (328-624) vs. 408 (182-517), unit = m, p = 0.009) and a higher distance saturation product (DSP) (421.2 (255.6-536.6), vs 373.2 (149.2-464.6), unit = m%, p = 0.02) ( Table 2) . Table 3 lists the correlation between distance and DSP based on clinical parameters recorded during the 6MWT. Age, synchrogram index, CAT, mMRC, and emphysema were all correlated with 6MWD and DSP. The synchrogram index was significantly correlated with distance (r = 0.42, p = 0.001) (Fig. 3c) and DSP (r = 0.41, p = 0.001) (Fig. 3d) . In the multivariable regression model, age, mMRC, and synchrogram index were independent variables predictive of distance. Age, synchrogram index, mMRC, emphysema were independent variables predictive of DSP (Tables 4, 5). Distance

Data was presented as mean ± SD or median (range) or number (%). Presence of emphysema was defined as presence of emphysema from high resolution computed tomography. Significantly difference between patients in the synchronized group and desynchronized group was defined as P < 0.05 COPD chronic obstructive pulmonary disease, BMI body mass index, HS-CRP high sensitivity C reactive protein, CAT chronic obstructive pulmonary disease assessment test, mMRC modified medical research council, ACO asthma and COPD overlap, DM diabetes mellitus, E/e′ ratio the ratio of the transmitral early peak velocity over early diastolic mitral annulus velocity, EF ejection fraction, LA left atrial, LABA Long acting beta agonists, LAMA long acting antimuscarinic agents, ICS inhaled corticosteroids, Triple LABA + LAMA + ICS, OCS oral corticosteroids, AE acute exacerbation 

This is the first study to study CRC of patients with COPD by clustering them into synchronized or desynchronized groups. Patients in the synchronized group had higher 6MWD and DSP compared with those in the desynchronized group. Our results identified the synchrogram index as a novel independent variable by which to predict DSP and 6MWD, which is a well-established predictor of mortality [38] and acute exacerbation [39] in patients with COPD. The synchrogram index depicts the heart-lung interaction; therefore, its relationship with the 6MWD suggests that it could potentially provide clinically useful information from a dimension other than 6MWD. Confirming the clinical applicability of this index to COPD patients (e.g., predicting mortality or acute exacerbation) will require following up patients for an extended period of time. The fact that CRC can be easily obtained using widely available non-invasive equipment means that it is applicable to a variety of healthcare environments, such as long-term homecare monitoring with the assistance of mobile technologies.

Researchers have previously demonstrated that 6MWD is an important predictor of survival in COPD patients [40, 41] . The limited walking distance demonstrated by COPD patients can be attributed to age [42] , impaired heart function with low ventricular ejection fraction (LVEF < 50%) [43] , and impaired respiratory function including desaturation [44] , emphysema severity [45] , dyspnea scores [46] , inspiratory capacity, and dynamic hyperinflation [47] . Several comorbid conditions, such as skeletal muscle dysfunction, impaired autonomic regulation, and nutritional factors, also contribute to exercise intolerance in patients with COPD [48, 49] .

Our study identified a correlation between 6MWD and age, mMRC, and synchrogram index. Since the severity of emphysema, desaturation, and diastolic heart failure do not show a significant contribution in our patients, we suggest that impaired pulmonary and heart function are not directly related to 6MWD in this study. Rather, we should consider nutritional status [50] , oxygen utilization by peripheral muscles, and/or negative a b c Fig. 3 a Distribution of synchrogram index in synchronized group and desynchronized group. b Scatterplot of synchrogram index against distance (m) from six-minute walking test of all patients. c Scatterplot of synchrogram index against distance saturation product (m%) from six-minute walking test of all patient cardiorespiratory-muscle interactions [51] . Previous studies posited that integrated cardiopulmonary function and muscle condition could reflect 6MWD [49] . These assertions are in line with our findings indicating that the synchrogram index (a quantification of heart-lung interactions), is an independent factor contributing to 6MWD. These findings warrant further investigation into the relationships among oxygen utilization by peripheral muscle, cardiopulmonary-muscle interactions, and muscle strength. DSP is a reliable factor to predict mortality among patients with bronchiectasis [52] , interstitial lung disease [53] , and COPD [54, 55] . In this study, patients in the desynchronized group present a lower DSP, implying an elevated likelihood of poor outcomes but need adequate follow-up duration to confirm. To our knowledge, this is the first study to evaluate factors that associated with DSP in COPD patients. Age, mMRC, synchrogam index and emphysema are independent variables to predict DSP. Emphysema is an independent factor in determining DSP but not 6MWD in this study, which may be related to the correlation of emphysema among desaturation during exercise [56] and its contribution to the desaturation component of DSP.

A strong heart-lung interaction may improve ventilation and perfusion matching, resulting in a better oxygen transport [57] . However, we did not observe any discrepancy between the synchronized and desynchronized groups in terms of saturation. This may be explained by the fact that we excluded patients who were using oxygen daily and by the reason that there was similar proportion of emphysema. Note that there may be a link between desaturation and coupling in those patients. In order to evaluate this relationship, it is necessary to explore COPD patients with chronic hypoxemic failure in the next program.

This study faced a few limitations. First, despite measuring and quantifying the coupling between respiration flow signals and IHR, we cannot conclude causality. Second, strict inclusion criteria prevented us from analyzing patients who were using oxygen on a daily basis, with the result that the study population was small, particularly in the desynchronized subgroup. Third, most of the patients in this study were male, with low CAT scores, and of East Asian descent, such that our findings are not necessarily generalizable to all COPD patients. Finally, this is a cross-sectional and preliminary data of a prospective-designed study. Due to the insufficient follow-up time, we cannot evaluate mortality outcomes and cardiac vascular events. We will continue monitoring the subjects in this study in order to observe the clinical impact of synchronization in heart-lung interactions.

This study first conducted the CRC analysis to describe heart-lung interactions of COPD patients. Asides from age and mMRC, synchrogram index is an independent variable that could predict 6MWD and DSP. 

Cardiovascular Comorbidities in Chronic Obstructive Pulmonary Disease (COPD)-current considerations for clinical practice

Defining the relationship between COPD and CVD: what are the implications for clinical practice? TherAd-vRespir Dis

The cardiopulmonary continuum systemic inflammation as 'common soil'of heart and lung disease

Chronic obstructive pulmonary disease and cardiac diseases. An urgent need for integrated care

Reduced intrathoracic blood volume and left and right ventricular dimensions in patients with severe emphysema: an MRI study

Influence of respiration on venous return in pulmonary emphysema

Effective bronchoscopic lung volume reduction accelerates exercise oxygen uptake kinetics in emphysema

Complex chronic comorbidities of COPD

COPD-related morbidity and mortality after smoking cessation: status of the evidence

Systemic manifestations of COPD

Cardiorespiratory interactions in humans and animals: rhythms for life

Pontomedullary transection attenuates central respiratory modulation of sympathetic discharge, heart rate and the baroreceptor reflex in the in situ rat preparation

Proposed respiratory 'gating' mechanism for cardiac slowing

Central regulation of heart rate and the appearance of respiratory sinus arrhythmia: new insights from mathematical modeling

Respiratory variations in pulmonary and systemic blood flow in healthy humans

Interactions between respiration and circulation

Quantitative analysis of cardiorespiratory synchronization in infants

Association between phase coupling of respiratory sinus arrhythmia and slow wave brain activity during sleep

Evolution of cardiorespiratory interactions with age

Fractal scale-invariant and nonlinear properties of cardiac dynamics remain stable with advanced age: a new mechanistic picture of cardiac control in healthy elderly

Phase transitions in physiologic coupling

Cardiorespiratory phase-coupling is reduced in patients with obstructive sleep apnea

Heart-lung interaction in a model of COPD: importance of lung volume and direct ventricular interaction

Global initiative for the diagnosis, management, and prevention of chronic obstructive lung disease: the 2020 GOLD Science Committee Report on COVID-19 & COPD

Characterisation of phenotypes based on severity of emphysema in chronic obstructive pulmonary disease

Toward a consensus definition for COPD exacerbations

Polish Respiratory Society guidelines for the methodology and interpretation of the 6 minute walk test (6MWT)

Standardization of spirometry: 1994 update

ATS committee on proficiency standards for clinical pulmonary function laboratories

Standardisation of spirometry

Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology

Synchrosqueezed wavelet transforms: an empirical mode decomposition-like tool

Heartbeat synchronized with ventilation

Automated quantification of the synchrogram by recurrence plot analysis

What is asthma− COPD overlap syndrome? Towards a consensus definition from a round table discussion

Clinical approach to the therapy of asthma-COPD overlap

AK-means clustering algorithm

Distance and oxygen desaturation during the 6-min walk test as predictors of long-term mortality in patients with COPD

Predicting outcomes from 6-minute walk distance in chronic obstructive pulmonary disease

Six-minute-walk test in chronic obstructive pulmonary disease: minimal clinically important difference for death or hospitalization

The 6-min walk distance: change over time and value as a predictor of survival in severe COPD

Chronic obstructive pulmonary disease (COPD) in elderly subjects: impact on functional status and quality of life

convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year • At BMC, research is always in progress. Learn more biomedcentral.com/submissions Ready to submit your research Ready to submit your research ? Choose BMC

What is the impact of impaired left ventricular ejection fraction in COPD after adjusting for confounders?

Determinants of poor 6-min walking distance in patients with COPD: the ECLIPSE cohort

Emphysema and DLCO predict a clinically important difference for 6MWD decline in COPD

Value of CAT and mMRC in clinical evaluation of chronic obstructive pulmonary disease

Inspiratory capacity, dynamic hyperinflation, breathlessness, and exercise performance during the 6-minutewalk test in chronic obstructive pulmonary disease

The body-mass index, airflow obstruction, dyspnea, and exercise capacity index in chronic obstructive pulmonary disease

Exercise response to oxygen supplementation is not associated with survival in hypoxemic patients with obstructive lung disease

The 6-minute walk test as a new outcome measure in Duchenne muscular dystrophy

Exercise intolerance in comorbid COPD and heart failure: the role of impaired aerobic function

Distance-saturation product of the 6-minute walk test predicts mortality of patients with non-cystic fibrosis bronchiectasis

The distance-saturation product predicts mortality in idiopathic pulmonary fibrosis

Prognostic value of variables derived from the six-minute walk test in patients with COPD: results from the ECLIPSE study

Prognostic ability of the distance-saturation product in the 6-minute walk test in patients with chronic obstructive pulmonary disease

Determinants of exerciseinduced oxygen desaturation including pulmonary emphysema in COPD: results from the ECLIPSE study

Respiratory sinus arrhythmia: why does the heartbeat synchronize with respiratory rhythm?

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

This result has been previously announced as a poster in 2020 Annual Congress of Taiwan Society of Pulmonary and Critical Care Medicine, Taiwan. We thank all the investigators and members of the Department of Thoracic Medicine for their efforts.

Authors' contributions TYL, HTW, PJC, CYL and YLL conceived and designed the analysis. SML, MHH, FTC, TYW, CSK and HCL contributed to the clinical and laboratory work for the study. YLL, HTW, and YCH designed and performed the statistical analyses. All authors were involved in data analysis, data interpretation, and preparation of the final manuscript. All the authors read and approved the final manuscript.

This work was supported by a grant from the Taiwan Ministry of Science and Technology, 107-2314-B-182A-145-MY2. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

The data sets analyzed during the current study are available from the corresponding author upon reasonable request.

The study was approved by the Ethics Committee of Chang Gung Memorial Hospital, Linkou, Taiwan (201702150B0). Patients provided signed informed consent prior to screening.

Not applicable.