AbstractRespiratory muscle dysfunction after open-heart surgery may influence the cardiopulmonary interactions. The purpose of this study was to examine the correlation between change in the neural respiratory drive (NRD) and change in heart rate variability (HRV) in patients submitted to open-heart surgery. An observational cross-sectional study was conducted among 32 participants. NRD was assessed via a surface electromyogram of the parasternal intercostal muscle (sEMGpara). Polar heart rate monitor was used to measure HRV during the deep breathing maneuver. Evaluations were performed on the day of admission and discharge. There were statistically significant differences in NRD and HRV indices between admission and discharge periods (P<0.05). The difference in peak root mean square of sEMGpara recorded during resting (ΔRMS sEMGpara tidal), during maximal inspiratory maneuver (ΔsEMGpara max), and its normalized values (ΔRMS sEMGpara%max) were significantly correlated with the difference in total power (ΔTotal power), mean of heart rate (ΔMeanHR), and mean of R to R intervals (ΔMeanRR) (r=−0.844, P=0.004, r=−0.835, P=0.005, and r=0.643, P=0.043, respectively). It can be concluded that NRD correlated well with HRV in patients who had undergone open-heart surgery.
INTRODUCTIONOpen-heart surgery, which is performed through a median sternotomy can cause postoperative pulmonary changes by numerous factors including surgical manipulation, anesthetic agents, cardiopulmonary bypass, chest drain, thoracotomy pain, and immobilization (Weissman, 2004; Wynne and Botti, 2004). These changes are associated with respiratory muscle dysfunction and restrictive pulmonary pattern postoperatively (El-Sobkey and Gomaa, 2011; Morsch et al., 2009).
Neural respiratory drive (NRD) is a noninvasive method to measure the neural output of the brainstem respiratory center indirectly by quantifying the electromyogram (EMG) of the parasternal intercostal muscles via surface electrodes (sEMGpara) (Reilly et al., 2011; Reilly et al., 2013). According to previous studies, the sEMGpara has been observed and indicated the load on respiratory muscles in cystic fibrosis (Reilly et al., 2011; Reilly et al., 2012) and chronic obstructive pulmonary disease patients (Jolley et al., 2009; Suh et al., 2015).
Heart rate variability (HRV) is a noninvasive marker used to investigate cardiac autonomic modulations that control the oscillation of instantaneous heart rates (McCraty and Shaffer, 2015). The fluctuation of cardiac rhythm has synchronism with the respiratory cycle, which is known as respiratory sinus arrhythmia (RSA). This phenomenon reflects the parasympathetic integrity markers over the cardiac sinus node. Heart rate was increased during inhalation due to parasympathetic withdrawal to the sinus node and reversed it during exhalation (Grossman and Taylor, 2007). The attenuation of cardiac vagal tone was found in patients who had undergone coronary artery bypass graft (Soares et al., 2005) and cardiac valve surgery (Lakusic et al., 2008).
Although the impact of cardiac surgery on cardiopulmonary functions has been demonstrated, respiratory load and its association with cardiac autonomic function have not been elucidated. Therefore, this study attempted to clarify the correlation between change in NRD and change in HRV among patients submitted to open-heart surgery.
MATERIALS AND METHODSStudy design and participantsThis observational cross-sectional study was accomplished in 32 participants who were undergoing open-heart surgery at Thammasat University Hospital, Thailand within the age of 35–60 years. The participants who had chronic heart failure, cardiac arrhythmia, myocardial infarction, unstable angina, implanted a cardiac pacemaker, uncontrolled diabetes mellitus, uncontrolled blood pressure, pulmonary or neurological diseases, experienced previous cardiac surgery and using mechanical ventilator more than 24 hr after surgery were excluded. This study was approved by the human ethics committee of Thammasat University (152/2560) and all participants have signed written informed consent.
NRD measurementThe sEMGpara was recorded at rest with normal tidal breath using bipolar surface electrodes (Kendall Arbo, Tyco healthcare, Neustadt, Germany) attached at the second intercostal space of both sides, lateral to the sternum for 3 cm. The amplifier and bandpass filter, which were represented as a frequency of EMG, was set at 1 kHz and between 10 Hz–2 kHz using a wireless EMG system (TeleMyo 2400T G2, Noraxon USA Inc., Scottsdale, AZ, USA), respectively. Data were analyzed using MyoResearch XP (ver. 1.07.25; Noraxon USA Inc.). The analog signal was converted to a digital signal at 10 kHz, and peak root means square per breath was calculated and averaged over 1 min (RMS sEMGpara tidal). EMG recordings at rest were normalized to the EMG signal obtained during a maximal static inspiratory pressure (RMS sEMGpara max) and the numerically largest EMG signal from 5 times repeated maneuver was used for normalization (sEMGpara% max).
HRV measurementRR intervals were recorded using a V800 Polar heart rate monitor (Polar Electro Ltd., Kempele, Finland) at a sampling rate of 1,000 Hz. HRV was recorded during deep breathing maneuver for 4 min in a supine position. Participants were instructed to perform a series of deep and slow inhalation and exhalation to provide a maximal pulmonary volume that varied from total lung capacity to residual volume. Each breathing cycle was performed for 10 sec, divided into 5 sec for inhalation and 5 sec for exhalation, which provides a maximal RSA response (Hayano et al., 1994; Song and Lehrer, 2003). Participants controlled their respiratory cycles via a pointer clock on the computer screen and received verbal feedback from the researcher.
HRV was analyzed in a linear method (time and frequency domains) by Kubios HRV software version 3.0.2 (Biosignal Analysis and Medical Imaging Group, University of Eastern Finland, Kuopio, Finland) following to the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology (1996).
Time-domain analysis was calculated from the mean of the longest RR intervals obtained during the expiratory phase divided by the mean of the shortest RR intervals obtained during the inspiratory phase (expiratory/inspiratory ratio, E/I), the difference between the mean of the highest heart rate obtained during the inspiratory phase and the mean of the lowest heart rate obtained during the expiratory phase (inspiratory-expiratory difference; ΔIE), the mean of RR intervals (MeanRR), the standard deviation of all normal RR intervals (SDNN), and the root mean square of the successive difference (RMSSD). E/I ratio and ΔIE has been indicated cardiac sympathovagal balance, MeanRR and SDNN have been indicated cardiac autonomic modulation, and RMSSD has been indicated cardiac vagal modulation, subsequently.
The fast Fourier transform (FFT) on the time series was utilized as a frequency domain. This FFT algorithm was applied to determine power spectrum density which consisted of low-frequency power (LF: 0.04–0.15 Hz) and high-frequency power (HF: 0.15–0.4 Hz). Spectral components were obtained in normalized units (nu). LF power has been predominantly indicated a sympathetic tone, HF power has been indicated a parasympathetic tone, while the LF/HF ratio has been indicated cardiac sympathovagal balance, respectively.
NRD and HRV indices were always evaluated in the afternoon on the day of admission and discharge. All participants were received conventional physiotherapy intervention following the standard cardiac rehabilitation protocol of the hospital including breathing exercises, airway clearance techniques, early mobilization, exercise training, and ambulation training during the postoperative period.
Statistical analysisShapiro–Wilk test was used to determine data distribution. Paired t-test and Pearson correlation were used for data fitting normal distribution, while the Wilcoxon signed-rank test and Spearman’s rank correlation were used for data not consistent with a normal distribution. The statistical significance was considered as P<0.05 using the IBM SPSS Statistics ver. 23.0 (IBM Co., Armonk, NY, USA).
RESULTSA total of 32 participants, 62.5% male (52.25±5.25 years) were included in the study (Table 1). The comparison of NRD and HRV indices between admission and discharge period are summarized in Table 2. For NRD parameters, RMS sEMGpara tidal and sEMGpara% max were significantly increased, while RMS sEMGpara max was significantly decreased in the postoperative period (all P<0.05). For HRV indices, MeanRR, SDNN, RMSSD, total power, and HFnu were significantly reduced, whereas ΔIE, E/I ratio, mean of heart rate (MeanHR), and LFnu were significantly increased at the discharge period (all P<0.05). There were no significant differences between time points in the LF/HF ratio. Correlations between change (Δ) in NRD and HRV indices are given in Table 3. There were significant negative correlations between ΔRMS sEMGpara tidal and ΔTotal power (P<0.01), and between ΔRMS sEMGpara max and ΔMeanHR (P<0.01). In contrast, ΔsEMGpara%max was a significantly positive correlation with ΔMeanRR (P<0.05).
DISCUSSIONThe finding of the present study was showed the correlation between the alteration of both NRD and HRV indices indicates that increased a load of respiratory muscle after open-heart surgery was associated with the attenuated of cardiac vagal control.
This study demonstrated postoperative hyperactivation of NRD which was detected by measuring neural control of parasternal intercostal muscle via surface EMG. Although the study of sEMGpara in patients who had undergone cardiac surgery is still lacking, the purpose mechanisms of this phenomenon may be involved in respiratory muscle workload and could be explained by physiological responses. The phasic mode of parasternal intercostal muscle plays a role in the thoracic expansion and increase lung volume during inspiration (De Troyer et al., 2005). It was possible that a load of parasternal intercostal muscle was increased due to respiratory muscle dysfunction and restrictive pulmonary pattern postoperatively (El-Sobkey and Gomaa, 2011; Morsch et al., 2009). The worsening of these respiratory functions is related to the surgical manipulation, anesthetic agents, thermal damage, cardiopulmonary bypass, chest drain, thoracotomy pain, and immobilization that lead to altering the respiratory mechanic and pulmonary compliance (Weissman, 2004; Wynne and Botti, 2004).
HRV indices obtained during deep breathing maneuver provided a marker of vagal integrity over the cardiac sinus node (McCraty and Shaffer, 2015). This study revealed a decrease in postoperative cardiac vagal control and reaffirmed with the previous study that showed a significant HRV reduction after coronary artery bypass graft (Soares et al., 2005) and cardiac valve surgery (Lakusic et al., 2008). Sinus arrhythmia is modulated by the synchronicity between the cardiovascular and respiratory system (Grossman and Taylor, 2007). In this context, alteration of respiratory variables can alter this phenomenon. Because the breathing frequency was controlled during deep breathing maneuver, it possible that lower HRV values are a result of the reduction of tidal volume. Altered pulmonary compliance may be the primary mechanism to explain this phenomenon. As a result of the postoperative restrictive pulmonary pattern, the vital capacity was decreased, thus reducing the range of the tidal volume displacement (El-Sobkey and Gomaa, 2011). So, although the patients were instructed to take deep and slow breathing during deep breathing maneuver, the tidal volume assembled could have been minimal due to the decreased lung capacity.
The correlation between alteration of NRD and HRV indices as a result of this study indicated that modifications of the respiratory muscle workload could profoundly influence the RSA magnitude and HRV behavior; slow and deep breathing will amplify RSA magnitude and HRV behavior, whereas fast and shallow breathing may contribute to reduced RSA magnitude and HRV behavior (Grossman and Taylor, 2007). Assuming that respiratory muscle weakness with restrictive pulmonary pattern leads to shallow breathing and considering that the neural output of the brainstem respiratory center drive to the respiratory muscles was shown to be increased, the baroreceptor and pulmonary stretch receptor may be early activated and may consequently be responsible for the fast central response in cardiac autonomic modulation (Yasuma and Hayano, 2004).
This study has some limitations which have to be pointed out. First, this study was not measured and controlled tidal volume and end-tidal of carbon dioxide during deep breathing protocol, but we instructed the participants to breathe as deep and slow as they can follow the pointer clock to provide a maximum lung volume that varied from total lung capacity to the residual volume also with monitored pulse oxygen saturation. Second, the RMS sEMGpara max was conducted in one maximal volitional maneuver, so normalized to the EMG signal at rest may not be the most value. Furthermore, all participants were received beta-blockers, which could influence on HRV, however, this study intended to investigate the real-life situations.
In conclusion, this study revealed that the patients who had undergone open-heart surgery were increased NRD that measured by sEMGpara and it was correlated to altering cardiac autonomic function, which was characterized by a decrease in HRV. This finding could provide the rehabilitation targets to restore respiratory workload and improve cardiac vagal tone.
ACKNOWLEDGMENTSWe would like to thank all participants and medical staffs from Thammasat University Hospital, Thailand, especially Assistant Professor Opas Sutdhabudha and Miss Chitima Kulchanarat for their involvement. This study was financially supported by a research grant from Chiangrai College and Faculty of Allied Health Sciences, Thammasat University (grant number 4/2560).
REFERENCESDe Troyer A, Kirkwood PA, Wilson TA. Respiratory action of the intercostal muscles. Physiol Rev. 2005;85:717–756.
El-Sobkey SB, Gomaa M. Assessment of pulmonary function tests in cardiac patients. J Saudi Heart Assoc. 2011;23:81–86.
Grossman P, Taylor EW. Toward understanding respiratory sinus arrhythmia: relations to cardiac vagal tone, evolution, and biobehavioral functions. Biol Psychol. 2007;74:263–285.
Hayano J, Mukai S, Sakakibara M, Okada A, Takata K, Fujinami T. Effects of respiratory interval on vagal modulation of heart rate. Am J Physiol. 1994;267:1 Pt 2. H33–40.
Jolley CJ, Luo YM, Steier J, Reilly C, Seymour J, Lunt A, Ward K, Rafferty GF, Polkey MI, Moxham J. Neural respiratory drive in healthy subjects and in COPD. Eur Respir J. 2009;33:289–297.
Lakusic N, Slivnjak V, Baborski F, Sonicki Z. Heart rate variability in patients after cardiac valve surgery. Cent Eur J Med. 2008;3:65–70.
McCraty R, Shaffer F. Heart rate variability: New perspectives on physiological mechanisms, assessment of self-regulatory capacity, and health risk. Glob Adv Health Med. 2015;4:46–61.
Morsch KT, Leguisamo CP, Camargo MD, Coronel CC, Mattos W, Ortiz LD, Lima GG. Ventilatory profile of patients undergoing CABG surgery. Rev Bras Cir Cardiovasc. 2009;24:180–187.
Reilly CC, Jolley CJ, Elston C, Moxham J, Rafferty GF. Measurement of parasternal intercostal electromyogram during an infective exacerbation in patients with cystic fibrosis. Eur Respir J. 2012;40:977–981.
Reilly CC, Jolley CJ, Ward K, MacBean V, Moxham J, Rafferty GF. Neural respiratory drive measured during inspiratory threshold loading and acute hypercapnia in healthy individuals. Exp Physiol. 2013;98:1190–1198.
Reilly CC, Ward K, Jolley CJ, Lunt AC, Steier J, Elston C, Polkey MI, Rafferty GF, Moxham J. Neural respiratory drive, pulmonary mechanics and breathlessness in patients with cystic fibrosis. Thorax. 2011;66:240–246.
Soares PP, Moreno AM, Cravo SL, Nóbrega AC. Coronary artery bypass surgery and longitudinal evaluation of the autonomic cardiovascular function. Crit Care. 2005;9:R124–131.
Song HS, Lehrer PM. The effects of specific respiratory rates on heart rate and heart rate variability. Appl Psychophysiol Biofeedback. 2003;28:13–23.
Suh ES, Mandal S, Harding R, Ramsay M, Kamalanathan M, Henderson K, O’Kane K, Douiri A, Hopkinson NS, Polkey MI, Rafferty G, Murphy PB, Moxham J, Hart N. Neural respiratory drive predicts clinical deterioration and safe discharge in exacerbations of COPD. Thorax. 2015;70:1123–1130.
Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation. 1996;93:1043–1065.
Weissman C. Pulmonary complications after cardiac surgery. Semin Cardiothorac Vasc Anesth. 2004;8:185–211.
Table 1Table 2
sEMGpara, surface parasternal electromyogram; RMS, root mean square; RMS sEMGpara tidal, peak root mean square of sEMGpara recorded during resting; RMS sEMGpara max, peak root mean square of sEMGpara recorded during maximal inspiratory maneuver; sEMGpara% max, resting parasternal EMG activity normalized to the maximal EMG evoked during maximal inspiratory maneuver; ΔIE, inspiratory-expiratory differences; E/I ratio, expiratory/inspiratory ratio; MeanRR, mean of R to R intervals for normal beats; MeanHR, mean of heart rate; SDNN, standard deviation of all R to R intervals; RMSSD, square root of the mean of the sum of the squares of differences between adjacent R to R intervals; LFnu, low frequency power in normalized units; HFnu, high frequency power in normalized units; LF/HF, the ratio between low and high frequency power. Table 3
sEMGpara, surface parasternal electromyogram; RMS, root mean square; ΔsEMGpara tidal, delta change of peak root mean square of sEMGpara recorded during resting; ΔsEMGpara max, delta change of peak root mean square of sEMGpara recorded during maximal inspiratory maneuver; ΔsEMGpara% max, delta change of resting parasternal EMG activity normalized to the maximal EMG evoked during maximal inspiratory maneuver; ΔIE difference, delta change of inspiratory-expiratory differences; ΔE/I ratio, delta change of expiratory/inspiratory ratio; ΔMeanRR, delta change of mean of R to R intervals for normal beats; ΔMeanHR, delta change of mean of heart rate; ΔSDNN, delta change of standard deviation of all R to R intervals; ΔRMSSD, delta change of square root of the mean of the sum of the squares of differences between adjacent R to R intervals; ΔLFnu, delta change of low frequency power in normalized units; ΔHFnu, delta change of high frequency power in normalized units; ΔLF/HF, delta change of the ratio between low and high frequency power. |
|