Participants
Patients with COPD were diagnosed and referred by a physician from the Department of Thoracic Medicine of Kaohsiung Municipal Siaogang Hospital, Kaohsiung City, Taiwan. The inclusion criteria based on the Global Initiative for Chronic Obstructive Lung Disease (2024) were: (1) Exposure to risk factors for COPD. (2) Presence of respiratory symptoms such as difficulty breathing, coughing, or sputum production. (3) A forced expiratory volume in the first second (FEV1) to forced vital capacity (FVC) ratio (FEV1/FVC ratio) < 70%. (4) Age between 40 and 80 years.
The exclusion criteria were as follows: (1) Presence of other lung diseases (e.g., asthma or lung cancer), heart disease (e.g., arrhythmia or pacemaker), cognitive impairment, seizures, or severe mental disorders. (2) Employment with shifts. (3) Unstable vital signs at rest, defined as systolic blood pressure ≥ 150 mmHg, heart rate ≥ 110 bpm, and oxygen saturation (SpO2) < 90%.
This study included 65 patients with COPD who completed the pre-test and were assigned to the HRVB (n = 30) and control groups (n = 35), with participants matched by age and sex. A total of 53 participants completed the post-test (26 in the HRVB group and 27 in the control group). A total of 12 participants were excluded from data analysis because of electrocardiogram (ECG) artifacts (n = 8) or dropout (n = 4). Detailed participant information has been published in our previous study (please see Wu et al., 2024).
The institutional review board of Kaohsiung Medical University Hospital, Kaohsiung City, Taiwan, approved this study (Approval No.: KMUHIRB-F(II)-20200060). Informed consent was obtained from all participants before the study began. After completing all study procedures, the participants received a gift card (approximately USD 20) for six visits in the HRVB group and a gift card (approximately USD 7) for the pre-test and post-test measurements in the control group.
Measurements
Demographic data (age and sex), smoking years and status, Global Initiative for Chronic Obstructive Lung Disease (GOLD) ABE Assessment Tool, COPD Self-Efficacy Scale (CSES), St. George's Respiratory Questionnaire (SGRQ), Beck Depression Inventory-II (BDI-II), and Beck Anxiety Inventory (BAI) were collected at both pre-test and post-test. After completed the psychological questionnaires, the 6MWT was conducted under baseline, walking, and recovery stages, and Borg Rate of Perceived Exertion Scale (Borg Scale) was measured under these three stages.
(1) GOLD ABE Assessment Tool was classified participants into three subgroups using the Modified Medical Research Council Dyspnea Scale (mMRC) and exacerbation history. Group A: mMRC score < 2, and 0 or 1 moderate exacerbation not leading to hospitalization in the past year. Group B: mMRC score ≥ 2, and 0 or 1 of moderate exacerbations, and not leading to hospitalization in the past year. Group E: mMRC score ≥ 2, and moderate exacerbations or ≥ 1 leading to hospitalization in the past year (Global Initiative for Chronic Obstructive Lung Disease, 2024).
(2) Chronic Obstructive Pulmonary Disease Self-Efficacy Scale was developed by Wigal et al. (1991) and was used to assess self-efficacy in patients with COPD. The CSES comprises 34 items rated on a Likert scale and includes five subscales: 12-item for negative affect, 8-item for intense emotional arousal, five-item for physical exertion, 6-item for weather/environment, and three-item for behavioral risk factors. The total score ranged from 34 to 170, with higher scores indicating greater self-efficacy.
Chiang et al. (2011) translated the Chinese version of CSES and used in Hong Kong. Its internal consistency reliability (Cronbach’s alpha) was .95, with a two-week rest-retest reliability of .91 (p < .001). In a sample of 30 COPD patients, the criterion-related validity between the CSES and the SGRQ was − .66 (p < .001). In this study, the CSES was translated into traditional Chinese and examined for internal consistency and reliability (Cronbach’s alpha) in 53 patients with COPD. Cronbach’s alpha was .98 for total scores, .95 for negative affect, .93 for intense emotional arousal, .90 for physical exertion, .88 for weather/environment, and .87 for behavioral risk factors. The six-week test-retest reliabilities were .93 (p < .001) for the total score, .91 (p < .001) for negative affect, .90 (p < .001) for intense emotional arousal, .94 (p < .001) for physical exertion, .92 (p < .001) for weather/environment, and .82 (p < .001) for behavioral risk factors. The criterion-related validity between CSES and SGRQ was r = − .30 (p = .031). These results indicate strong psychometric properties of the CSES.
(3) St. George’s Respiratory Questionnaire was developed by Jones et al. (1992), and was used to assess the impact of COPD-related dyspnea on daily life. The SGRQ contains 50-item on a Likert scale and three subscales: eight-item assessing respiratory symptoms, 16-item measuring activity limitations in daily life, and 26-item evaluating the impact of illness. The total SGRQ score ranges from 0–100%, with higher scores indicating a poor quality of life and lower scores indicating a better quality of life. A score higher than 25% is considered the cutoff for poor quality of life (Global Initiative for Chronic Obstructive Lung Disease, 2024). The Chinese version of the SGRQ was translated by David et al. (2004) for use in Hong Kong. The two-week rest-retest reliability ranged from .70 to .86, and the internal consistency reliability ranged from .74 to .98. The SGRQ and 36-item Short Form Health Survey (SF-36) validity ranged from .33 to .92.
In this study, the SGRQ was translated into traditional Chinese, and its internal consistency reliability (Cronbach’s alpha) was assessed in 53 patients with COPD. Cronbach’s alpha was .91 for the total score, .76 for symptoms, .87 for activity, and .85 for impact. The six-week rest-retest reliability (n = 27) was .88 (p < .001) for the total score, .88 (p < .001) for symptoms, .92 (p < .001) for activity, and .83 (p < .001) for impact. The validity of the SGRQ and BODE index in this study was r = .51 (p < .001). These results confirmed the strong psychometric properties of the SGRQ.
(4) Beck Depression Inventory-II comprises 21 items on a Likert scale designed to evaluate depressive symptoms, including somatic and cognitive depression subscales. The total BDI-II score ranges from 0 to 63, with higher BDI-II scores indicating more severe depressive symptoms. A total score below 13 indicated normal depression, 14–19 indicated mild depression, 20–28 indicated moderate depression, and ≥ 29 indicated severe depression (Beck et al., 1996). The Chinese version of the BDI-II, translated by Chen (2000), has demonstrated strong psychometric properties, including a Cronbach’s alpha of 0.94, split-half reliability of 0.91, and a correlation of 0.69 with the Chinese Health Questionnaire of 0.69 (Lu et al., 2002).
(5) Beck Anxiety Inventory comprises 21-item on a Likert scale to evaluate anxiety symptoms. It includes subscales assessing somatic and cognitive anxieties. The total BAI score ranges from 0 to 63, with higher scores indicating more severe anxiety symptoms. A total score below 7 indicated normal anxiety, 8–15 indicated mild anxiety, 16–25 indicated moderate anxiety, and ≥ 26 indicated severe anxiety (Beck et al., 1988). The Chinese version of BAI, translated by Lin (2000), has shown good psychometric properties, with a Cronbach’s alpha of .95, split-half reliability of .91, and a correlation of .72 with the Hamilton Anxiety Scale of 0.72 (Che et al., 2006).
(6) Six-Minute Walking Test was used to measure the exercise capacity for patients with COPD. Based on the guidelines and manual of the American Thoracic Society (ATS Committee on Proficiency Standards for Clinical Pulmonary Function Laboratories, 2002), participants were asked to walk at a fast pace for 6 minutes along a flat corridor. This study based on the Vagal tank theory (Laborde et al., 2018) to examine the HRV changes from resting baseline to walking (HRV reactivity), and from walking to recovery (HRV recovery). Therefore, ECG data was collected under resting baseline, walking, and recovery using the ProComp Infiniti version 6.0 with ECG sensor with 2048 Hz/s sampling rate (Thought Technology, Montreal, QC, Canada) at baseline and recovery stages. The Zephyr BioHarness™ with 250 Hz/s sampling rate (Zephyr Technology Corp., Annapolis, Maryland, USA) attached to chest strap and transfer to smartphone via Bluetooth during walking stage.
(7) The Borg Rate of Perceived Exertion Scale was developed by Dr. Gunnar Borg and is used to rate how hard individuals feel they are exerting themselves during baseline, walking, and recovery stages of the 6MWT (Borg, 1982). The Borg Scale typically ranges from 0 (nothing at all) to 10 (maximum exertion).
Heart Rate Variability Biofeedback Training Equipment
During the HRVB training, a ProComp Infiniti version 6.0 (Thought Technology, Montreal, QC, Canada) with a lead II ECG sensor (sampling rate: 2048 Hz/s) and a respiration sensor (sampling rate: 256 Hz/s) were placed on the participants’ chests to collect raw ECG and respiration signals. During the home practice, a wearable Eureka device with a photoplethysmography (PPG) sensor (Finesse Lifecare Co., Ltd., Taiwan) was placed on participants’ earlobes and connected to a smartphone via Bluetooth. The accompanying breathing app guided patients to slow their breathing rate and provided feedback after the home training.
Research Design and Experimental Procedure
A two-factor mixed design encompassing a between-participants design (comparing the HRVB and control groups) and a within-subjects design (pre-test and post-test). Patients completed the psychological questionnaires, and a research assistant ensured that patients’ vital signs were stable (SBP < 150 mmHg, heart rate < 110 bpm, and SpO2 ≥ 90%). Raw ECG signals and breathing rates were recorded using a lead II ECG and respiration sensor under the resting baseline at pre-test and post-test, during the 6MWT, and during the HRVB training period. This study primarily focused on the effects of the HRVB on psychological characteristics.
Heart Rate Variability Biofeedback Protocol
Heart Rate Variability Biofeedback Protocol
This study implemented Lehrer's HRVB protocol (Lehrer et al., 2000, 2013), with modifications based on our prior clinical studies: (1) To accommodate the aging patients and their dyspnea symptoms, breathing rates were gradually reduced until reaching each individual’s resonance frequency, as some patients found slower breathing intolerable during the initial session. (2) In the first session, patients were introduced to self-guided muscle relaxation through an audiotape and were guided to breathe in a relaxed manner. (3) Based on our previous studies, the HRVB training was structured into six weekly sessions over six consecutive weeks (Lin et al., 2015).
The patients’ optimal resonance frequency was determined during the first session by assessing the highest low frequency (LF) power across five breathing rates (6.5, 6, 5.5, 5, and 4.5 breaths per minute), followed by self-guided relaxation training. In sessions 2–6, the patients were trained in the abdominal and pursed lip-breathing techniques, gradually adjusting their breathing rates from baseline to achieve their target resonance frequency. In the third session, the Eureka (Finesse Technology Ltd., Zhubei, Hsinchu, Taiwan) were instructed to engage in home practice sessions twice daily for 5 min, 5 days a week. After each home practice session, the patients emailed their raw data to the researchers and received feedback in the following session. Psychoeducation on breathing techniques was provided through a brochure as part of the HRVB training.
Data reduction and statistical analysis
During data processing, a research assistant visually inspected R spike from ECG raw signals and peak wave from PPG signals, removed arrhythmia and movement artifacts, and then converted into HRV indices using Kubios-Premium version 3.5 (Kubios Co., Kuopio, Finland). ECG artifacts or arrhythmias that occurred in more than 20% of a 5-minute period were excluded from the statistical analysis (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). The HRV indices included the time and frequency domains of HRV: (1) Standard deviation of the normal-to-normal intervals (SDNN), representing total HRV; (2) Root mean square of successive differences (RMSSD), reflecting parasympathetic activation; (3) low-frequency power (LF, 0.04–0.15 Hz, influenced by both sympathetic and parasympathetic activation or baroreceptor gain); (4) high-frequency power (HF, 0.15–0.4 Hz, indicative of parasympathetic activation; (5) low-frequency/high-frequency ratio (LF/HF ratio, reflecting sympathetic activation); (6) total power (TP, representing total HRV). To meet the homogeneity assumption of HRV, the HRV indices were transformed using the natural logarithms of lnSDNN, lnRMSSD, lnLF, lnHF, and lnTP.
To assess the accuracy and consistency from three devices, we collected simultaneously from 15 healthy students using the ProComp Infiniti, Zephyr BioHarness™ and Eureka for five-minute resting baseline. We explore interbeat interval (IBI) with 256 Hz/s sampling rate from the ProComp Infiniti, IBI with 250 Hz/s sampling rate from the Zephyr BioHarness™, and peak to peak interval with 256 Hz/s sampling rate from the Eureka, and then analyzed using Kubios-Premium version 3.5 (Kubios Co., Kuopio, Finland) to derive HRV indices.
Due to Food and Drug Administration approval of ProComp Infiniti, this study analyzed the correlations between ProComp Infiniti and Zephyr BioHarness™ to confirm the, finding r = 1.000 (p < .001) for mean IBI, low frequency (LF), high frequency (HF), and low frequency/high frequency ratio (LF/HF ratio). High correlations were observed between ProComp Infiniti and Eureka: r = .998 (p < .001) for mean IBI, r = .990 (p < .001) for SDNN, r = .943 (p < .001) for the root mean square of successive heartbeat interval differences, r = .999 (p < .001) for LF, r = .998 (p < .001) for HF, and r = .992 (p < .001) for LF/HF ratio. These findings indicate a high level of consistency between the ProComp Infiniti, Zephyr BioHarness™, and Eureka devices.
All statistical analyses were conducted using SPSS Windows software version 24.0; IBM Corp., Armonk, NY). A t-test and chi-square test were employed to compare the pre-test demographic data and psychological questionnaires between the HRVB and control groups. Two-way analysis of variance (ANOVA) was used to examine the interaction effects of Group (HRVB vs. control) × Time (pre-test and post-test) on the psychological questionnaires (CSES, SGRQ, BDI-II, BAI, and HRV indices). When a significant Group × Time interaction effect was found, a simple main effects analysis with a Bonferroni post hoc comparison was conducted to examine the differences between groups or across time points in the psychological questionnaires. If the Group × Time interaction effect was not significant, a main effect analysis was conducted to examine the differences between groups or across time points in the psychological questionnaires. Additionally, the Mauchly sphericity test was applied, and the Greenhouse-Geisser correction was used to adjust for any violation of sphericity in the repeated-measures ANOVA.
Partial eta squared (np2) was used to calculate the effect size. Effect size (np2) were interpreted as follows: 0.01 ≤ np2 < 0.58 for a small effect size, .058 ≤ np2 < 0.138 for a medium effect size, and np2 ≥ 0.138 for a large effect size (Cohen, 1988).