Association of Short-Term Heart Rate Variability With Breast Tumor Stage

Cardiac autonomic modulation, assessed by heart rate variability (HRV), is associated with tumor pathogenesis and development as well as invasion and metastasis. This study aimed to examine this association in breast cancer (BC) patients. A total of 133 patients (average age 49.2years) with BC or benign breast tumors were divided into three groups: benign group, early-stage group, and advanced-stage group. About 5-min resting ECG was collected for the analysis of linear and nonlinear HRV parameters. Multiple logistic regression models were performed to test the independent contribution of HRV to breast tumor stage. The advanced-stage group had significantly reduced HRV compared to the benign and early-stage groups. In particular, for each 1-SD increase in SD2, SD of normal-to-normal intervals, very-low frequency, total power, and low frequency, the odds of having advanced staging decreased by 69.3, 64.3, 58.3, 53.3, and 65.9%, respectively. These associations were independent of age, body mass index, mean heart rate (HR), and respiratory rate (RR). These findings suggest an association between HRV and breast tumor stage, and HRV parameters may help construct an effective early diagnostic and clinical prognostic model.

Introduction

Breast cancer (BC) is the most common cancer among women worldwide. In 2020, approximately 2.26 million new cases of BC were diagnosed. BC represents approximately 11.7% of all cancers and is the fifth leading cause of cancer deaths worldwide (680,000 deaths; Sung et al., 2021).

The vagus nerve, a major part of the parasympathetic nerve system, regulates the development and progression of cancer (De Couck and Gidron, 2013). Increased vagal nerve activity has a profound inhibitory effect on oxidative stress, DNA cell damage, inflammation, and sympathetic nervous system overreaction (Maki et al., 2006; Tsutsumi et al., 2007; Tracey, 2009; De Couck et al., 2012). Studies have demonstrated the bidirectional relationship between cancer and decreased vagal nerve activity (De Couck et al., 2018). Therefore, study of the vagus nerve could provide valuable prognostic information and guide therapy in breast cancer.

Heart rate variability (HRV) is a promising biomarker that can be used to evaluate autonomic nervous system function and it may be associated with vagal nerve function (Vanderlei et al., 2009; Karvinen et al., 2013). One group showed the clinical significance of HRV in patients with BC (Arab et al., 2016). However, few studies have explored the association of HRV with tumor-node-metastasis (TNM) in patients with BC.

In a previous study, patients with BC had lower the SD of all normal-to-normal intervals (SDNN) and the root mean square of successive interval differences (RMSSD) than women without BC, no matter how long after the surgery (Palma et al., 2016). In addition, Arab et al. (2018) found that SDNN and RMSSD negatively correlated with BC staging after analyzing time and frequency domains of HRV in patients. In patients with advanced-stage BC, low SDNN might be correlated with a poor prognosis. Kim et al. (2010) reported that an SDNN value of <21.3ms in brain metastasis predicted poor survival, while Wang et al. (2013) found that an SDNN value of <10ms in terminal-stage cancer predicted poor survival. These studies suggested a potential role of time-domain parameters of HRV as a prognostic factor. However, further research is warranted to clarify the potential role of time-domain parameters in the survival of patients with advanced-stage BC.

Previous studies indicated that higher resting high-frequency power (HF) was strongly associated with longer overall survival in patients with recurrent or metastatic BC (Giese-Davis et al., 2015). Chiang et al. (2010) found that the survival time of patients with terminal hepatocellular carcinoma was significantly related to HF. These studies suggested that, as a long-term predictor of survivors, HF may offer early estimation of clinical prognosis for cancer patients. More importantly, the results of these studies also indicated that HF strongly positively correlated with prognosis, particularly in patients with advanced-stage cancer. The vagal nerve activity might be of more importance in advanced stages.

Therefore, this study aimed to compare short-term HRV in the tumor stage of patients with BC. It was hypothesized that HRV in patients with advanced cancer would be lower than that in patients with early-stage disease. The analysis confirmed the aforementioned hypothesis and had clinical implications: HRV might be a potentially feasible tool in clinic to evaluate the prognosis of BC. More importantly, clinicians could ascertain patients at risk for disease progression through the long-term monitoring of HRV.

Materials and Methods

Participants and Procedures

This study followed the regulations of the National Research Ethics Committee and obtained the approval of the Clinical Medical Research Ethics Committee of the First Affiliated Hospital of Bengbu Medical College (Bengbu, Anhui, China; registration number: 2019KY031). In this study, women diagnosed with breast tumors were selected by the pathological examination from 2019 to 2020. All participants volunteered for this study and provided informed consent.

The function of the vagal nervous system was assessed by analyzing HRV using an ECG recorder (HeaLink-R211B; HeaLink Ltd., Bengbu, China). The sampling rate of the ECG signal was 400Hz. ECG data were collected at room temperature (23±1°C) and always 3days prior to radiotherapy/operation. Participants were explained the ECG collection procedure and were asked to assume the supine position and keep quiet during ECG examination. The ECG acquisition time using V5-lead was 5min.

The following conditions that are known to alter HRV were used as exclusion criteria: (1) diabetes mellitus; (2) heart diseases; (3) use of anti-arrhythmic drugs or beta-blockers; (4) pacemaker; (5) poor ECG quality; (6) ectopic beats (>10% of all beats); and (7) chemotherapy or surgery in the 3weeks before the examination. Therefore, our study analyzed the data of 133 participants.

HRV Analysis

The ECG R peaks were extracted using an algorithm based on the Pan–Tompkins algorithm (Pan and Tompkins, 1985). The technical and physiological artifacts within R-R intervals (RRI) were corrected by applying an automatic artefact correction algorithm. Subsequent HRV indices for both linear (time and frequency domain) and nonlinear methods were calculated.

Commonly used time-domain indices include SDNN and RMSSD. SDNN, a total variability index, represents the involvement of all cyclic components. RMSSD reflects parasympathetic activity (Camm et al., 1996; Vanderlei et al., 2009).

The RRI time series was converted into power spectral analysis to analyze the frequency domain, and the power spectral density was obtained using the Fast Fourier Transform algorithm. Prior to frequency-domain analysis, the RRI time series was evenly resampled at 4Hz using cubic spline interpolation. Frequency-domain parameters included total power (TP, 0–0.4Hz), high-frequency power (HF, 0.15–0.4Hz), low-frequency power (LF, 0.04–0.15Hz), very-low frequency power (VLF, 0–0.04Hz), and the ratio of LF to HF (LF/HF). LF and HF parameters were expressed in normalized units: normalized HF [HF n.u. =HF/(TP – VLF)] and normalized LF [LF n.u. =LF/(TP – VLF); Montano et al., 1994; Camm et al., 1996; Vanderlei et al., 2009].

The LF corresponds to the co-regulation of sympathetic and vagal nerve tones, HF indicates the vagal nerve tone, and the LF/HF reflects interactions of both sympathetic nervous system and parasympathetic nervous system, but they are limited to the case where the respiratory frequency is in the HF band (Hernando et al., 2016; Varon et al., 2019). The difference in respiratory rate (RR) will lead to the analysis of HRV in the standard frequency band cannot accurately estimate the activity of autonomic nervous system. For example, when the RR is higher than the upper limit of the HF band, vagus activity may be underestimated. In contrast, when the RR is within the LF band, sympathetic activity is overestimated and vagus activity is underestimated. Therefore, respiratory influences need to be separated in order to better estimate the activity of sympathetic and vagus nerves (Varon et al., 2019). An estimate of RR was calculated using an ECG-derived respiration approach (Moody et al., 1985). It is important to check whether RR is below 0.15Hz or higher than 0.4Hz in all the enrolled subjects, in order to trust in the interpretation of LF and HF related indices.

Each RRI time series included eight nonlinear HRV indices, including approximate entropy (ApEn), sample entropy (SampEn), Poincare plot: SD1, SD2, and SD2/SD1, detrended fluctuation analysis (DFA): α1 and α2, and correlation dimension (CD). The estimated ApEn and SampEn depended on three parameters: the embedding dimension m, the tolerance value r, and the data length N. The parameters are set as m=2 and r=0.2σ, where σ was the SD of each realization (Camm et al., 1996; Vanderlei et al., 2009; Voss et al., 2009; de Godoy, 2016).

All the above processing steps were performed using the Kubios HRV Premium software (version 3.1.0, Kubios Oy, Kuopio, Finland).¹

Breast Tumor Groups

Patients with breast tumors were divided into three groups: benign group, early-stage group, and advanced-stage group. Patients with benign breast tumors (i.e., benign epithelial proliferations, intraductal papilloma, phyllodes tumor, breast hyperplasia, and fibroadenomas) were selected as the controls. According to the National Comprehensive Cancer Network Clinical Practice Guidelines TNM staging version 3.2020 (Gradishar et al., 2020), the remaining participants were divided into early-stage and advanced-stage groups. The early-stage group consisted of T1–2, N0–1, and M0, and T3N0M0 cancers, while the advanced-stage group consisted of T0–4, any N, and M0–1 cancers. The non-advanced-stage group included controls and early-stage patients.

Statistical Analysis

Descriptive statistical data were expressed as mean (SD), median (Q1, Q3), or percentage. The Shapiro–Wilk test was used to test the normality of HRV indices. A chi-square test was used to analyze the difference between the two cancer groups. Dependent variable analyses for linear and nonlinear HRV parameters were separately conducted using parametric and nonparametric tests. One-way ANOVA was used to calculate normal data and Fisher’s least significant difference (LSD) was used for multiple comparisons between groups. The Kruskal–Wallis test was used to analyze non-normal data. Finally, separate multiple logistic regression models were performed with breast tumor stage as an outcome and with each significant HRV parameter set as a predictor while adjusted for age, body mass index (BMI), mean heart rate (HR), and RR. SPSS Statistics 25.0 (IBM Corp., Chicago, Illinois, United States) was used, and a value of p<0.05 was considered statistically significant.

Results

Table 1 presents the demographics and HRV indices of patients with breast tumors. In TNM staging, the early-stage group mainly comprised patients with stages T1–2, N0–1, and M0 cancers, while the advanced-stage group commonly comprised patients with stages T1–3, N2–3, and M0 cancers. Twelve patients in the advanced-stage group had distant metastases to the brain, bones, and lungs. Among patients with stages I–IV, the early-stage group mostly comprised stages Ia, IIa, and IIb; however, the advanced-stage group mostly comprised stage IIIa, IIIc, and IV. Invasive carcinoma with no special type is mostly frequent in patients with BC. The early-stage group had a more noninvasive type of BC compared to the advanced-stage group. The results of the molecular typing revealed no significant differences between the groups (Table 2).

TABLE 1

Table 1. Demographics and heart rate variability (HRV) of breast tumor patients.

TABLE 2

Table 2. Clinical characteristics of breast cancer (BC) patients.

Differences in SDNN (p<0.001), RMSSD (p=0.006), VLF (p<0.001), LF (p<0.001), HF (p=0.005), TP (p<0.001), SD1 (p=0.007), SD2 (p<0.001), and CD (p<0.001) between groups were shown using one-way ANOVA and Kruskal–Wallis test. Furthermore, SDNN, RMSSD, VLF, LF, HF, TP, SD1, SD2, and CD were significantly decreased in the advanced-stage group compared to the corresponding values in the benign and early-stage groups. However, differences in HRV indices between the benign and early-stage groups were not significantly different. Moreover, there were no statistically significant differences in the LF n.u, HF n.u, LF/HF, SD2/SD1, ApEn, SampEn, α₁, and α₂ among the groups (Table 3).

TABLE 3

Table 3. Comparison of HRV parameters among the benign, early-stage, and advanced-stage groups.

With SDNN=20ms as the cutoff value, the subjects were divided into two subgroups: a low-SDNN subgroup (SDNN<20ms; n=31) and a high-SDNN subgroup (SDNN>20ms; n=102). There was no statistically significant difference in age between the non-advanced and advanced groups (low-SDNN subgroup, p=0.311; high-SDNN subgroup, p=0.218; Figure 1).

FIGURE 1

Figure 1. Subgroup analysis of age between non-advanced and advanced groups.

To test the independent contribution of HRV to breast tumor stage, we performed multiple logistic regression models after adjusting for age, BMI, mean HR, and RR (Table 4). The associations of tumor stage with SD2, SDNN, VLF, TP, and LF were significant in logistic regression analysis. Specifically, for each 1-SD increase in SD2, SDNN, VLF, TP, and LF, the odds of having advanced staging decreased by 69.3% [odds ratio (OR): 0.307, 95% CI: (0.155, 0.606)], 64.3% [OR: 0.357, 95% CI: (0.183, 0.694)], 58.3% [OR: 0.417, 95% CI: (0.216, 0.807)], 53.3% [OR: 0.467, 95% CI: (0.239, 0.910)], and 65.9% [OR: 0.341, 95% CI: (0.130, 0.898)], respectively. The associations of tumor stage with CD, RMSSD, SD1, and HF were not significant in the logistic regression models.

TABLE 4

Table 4. Results from Logistic regression models (adjusted for age, BMI, Mean HR, and RR).

Discussion

This study aimed to compare the HRV of patients with breast tumors (benign tumors, early-stage BC, and advanced-stage BC) and evaluate the feasibility of HRV as a tool for the early diagnosis and prognosis of BC patients. Our results revealed that patients with advanced-stage BC had lower HRV than those with benign tumors and early-stage BC. However, no statistically significant difference was observed in the HRV indices between the benign and early-stage groups. After adjusting for age, BMI, Mean HR, and RR, our results showed that SD2, SDNN, VLF, TP, and LF were associated with tumor stage.

The vagus nerve also called the wandering nerve, works via many neurotransmitters and plays an important role in multiple systems, such as cardiovascular, neuroendocrine, and immunological (Tracey, 2009). Studies have showed that the vagal nerve system transmits a variety of signals to the brain in order to restore the body to a steady state (Tracey, 2009; Ohira et al., 2013). Moreover, published preliminary studies have also demonstrated that vagal nerve tone is vital in the prognosis of cancer. A higher vagal nerve tone may protect cancer patients by reducing inflammation (De Couck et al., 2016). Measuring the HRV is a noninvasive approach of measuring vagal nerve activity (Vanderlei et al., 2009). HRV is known to be associated with vagal nerve tone (r=0.88; Kuo et al., 2005). In recent years, various linear and nonlinear methods have been applied to analyze the time series of heartbeat cardiac intervals, reflecting the physiological and pathological information contained in the HRV signal from different angles.

Few studies have explored the association of HRV with TNM in patients with breast cancer. In one study, Mouton et al. (2012) examined the data of 72 patients with colorectal cancer and found that the baseline HRV could predict the carcinoembryonic antigen levels at 12months. Moreover, they found that SDNN<20ms was associated with significantly higher CEA at 12months. Arab et al. (2016) showed the clinical significance of HRV in patients with breast cancer.

Several studies indicated that higher resting HF was strongly associated with longer overall survival in patients with recurrent or metastatic BC (Giese-Davis et al., 2015). Chiang et al. (2010) found that the survival time of patients with terminal hepatocellular carcinoma was significantly related to HF. More importantly, the results of these studies also indicated that HF strongly positively correlated with prognosis, particularly in patients with advanced-stage cancer. However, despite the interesting results of our study, we did not find any statistically significant differences in the HF among the three cancer stage groups. The difference in mean HR and RR may cause HF to not accurately estimate the cardiac autonomic regulation activity. In our study, the difference was not statistically significant in mean HR and RR between the groups, and the RR of all the enrolled subjects were in the HF band. Therefore, the interpretation of spectral HRV indices as autonomic cardiac regulation markers in our study was more accurate and reliable. The findings of the present study could be verified through a prospective study performed over a longer follow-up period.

Similar to the above results, the median HF at baseline was 122 (IQR, 55, 255) ms² in our study, and the median HF values in benign, early-stage, and advanced-stage BC were 151 (78, 300), 151 (55, 270), and 60 (23, 136), respectively, indicating that a higher HF may be associated with benign or early-stage BC, and consequently a better survival rate. However, unlike Mouton et al. (2012), we used a cutoff of SDNN=20ms, and found that the SDNN showed a significant inverse association with tumor stage. A lower SDNN was found to be associated with more advanced-stage BC.

The interpretation of LF/HF is controversial. Malliani et al. (1991) showed that the LF/HF could reflect the balance between the vagus nerve and the sympathetic nerve. But LF/HF has been largely criticized as a marker of sympatho-vagal balance (Billman, 2011, 2013). First, LF power is not a pure index of sympathetic nerve, it may also be affected by vagus nerve and other unspecified factors. Second, sympathetic and vagus nerve can be simultaneously active, and their interaction is complex and nonlinear. Third, respiratory parameters and mechanical factors will also cause uncertainty in the contribution of sympathetic and vagus nerve to LF/HF. Finally, HR can affect LF/HF independently of cardiac autonomic nerve activity (Billman, 2013). In our study, no statistically significant difference was observed in LF/HF among the three cancer stage groups. This could be because of the complex physiological basis of LF/HF and the other unidentified factors.

Linear methods cannot be used to describe properly the complex nonlinear behavior, which is predominant in human systems. Therefore, it is necessary to search the novel indexes to reflect the correlation and the complexity characteristics of the HRV signal. The characteristic of HRV nonlinear analysis can better express the irregularity, complexity, and other dynamic characteristics of heartbeat fluctuations. Across various studies in the field of cardiovascular disorders, nonlinear dynamical HRV analysis is significantly superior to linear time-domain and frequency-domain methods (Mäkikallio et al., 2001; Stein et al., 2005; Voss et al., 2009; de Godoy, 2016). Some preliminary studies also explored the correlation between cancer and several nonlinear heartbeat dynamics measurements. Bettermann et al. (2001) showed that the variability, complexity, or rhythmicity of HRV in patients with BC was lower than that in patients with diabetes and age-matched healthy women. In particular, while comparing patients with and without BC metastasis, patients with metastasis had lower ApEn compared with those without metastasis. Shi et al. (2019) explored the perturbations of HRV nonlinear dynamical patterns to predict the increase in the severity of gastric cancer and found that nonlinear HRV parameters were the markers of autonomic nervous function to tumor progression.

Although, analysis of HRV by methods based on nonlinear dynamics do not reflect vagal or sympathetic regulation, we found significant correlations between time- and frequency domain indices and some of the nonlinear HRV parameters in patients with BC. For example, SampEn, α₁, and CD correlated with RMSSD [SampEn (r=0.546, p<0.001), α₁ (r=−0.564, p<0.001), and CD (r=0.800, p<0.001)], and also correlated with HF [SampEn (r=0.521, p<0.001), α₁ (r=−0.530, p<0.001), and CD (r=0.792, p<0.001)]. This illustrates that the nonlinear parameters also contain the component of time-frequency domain index. The nonlinear analysis method is still in the preliminary exploration stage, and the exact physiological and pathological background has not been fully clarified. The findings of the current study might provide new evidence on the role of nonlinear HRV in cancer. Further studies are needed to clarify the correlation of nonlinear HRV as a long-term predictor of survival.

This study had some limitations. First, its cross-sectional study design was a major notable limitation. The correlation between HRV parameters and outcomes could not be inferred. Second, a comprehensive understanding of the connections of nonlinear HRV with BC prognosis is currently lacking. Third, more background variables, such as physical activity, stress levels, use of medications, and other relevant medical variables, could not be included. To address these limitations, studies with larger sample sizes, more detailed background variables, and a prospective design should be conducted to clarify the correlation of linear and nonlinear HRV parameters with BC prognosis.

Conclusion

This novel study investigated linear and nonlinear HRV parameters in breast tumor groups. It found that the HRV was related to BC staging, indicating a correlation between tumor and HRV. The results of our study showed that patients with advanced-stage BC had lower HRV and might have a poor prognosis, and demonstrated that nonlinear HRV parameters might predict tumor staging in patients with breast tumors. Nonlinear approaches are of great significance in coping well with the nonstationary and nonlinear nature of heartbeat fluctuations. It is suggested that the combined measurement of linear and nonlinear HRV parameters may benefit future investigations. Researchers should identify a comprehensive biomarker for predicting BC prognosis by leveraging existing linear methods and nonlinear indicators. In addition to evaluating tumor stage, vagal nerve activity should be considered to estimate the prognosis of a cancer patient. Vagal nerve activity can be easily assessed and has the potential to provide healthcare professionals with incremental information based on the treatment plans. Future research should investigate the therapeutic potential of vagal nerve activation in cancer treatments through different supportive therapies such as relaxation, exercise interventions, and Traditional Chinese Medicine treatments (Niederer et al., 2013).

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics Statement

This study followed the regulations of the National Research Ethics Committee and obtained the approval of the Clinical Medical Research Ethics Committee of the First Affiliated Hospital of Bengbu Medical College, China. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

BS: conceptualization, resource allocation, and review and editing of the manuscript. SW and MC: data collection, interpretation of the results, and manuscript preparation. JW: data analysis. YZ: supervision and resource allocation. All authors contributed to the article and approved the submitted version.

Funding

This study was funded by the “512” Outstanding Talents Fostering Project of Bengbu Medical College under Grant BY51201312, the Scientific Research Innovation Project of Bengbu Medical College under Grant BYKC201905, and the Students Scientific Research Innovation Project of Bengbu Medical College under Grant byycx20103.

Conflict of Interest

A direct family member of BS owns stock in HeaLink Ltd., Bengbu, China.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

1. ^https://www.kubios.com

References

Arab, C., Dias, D. P. M., Barbosa, R. T. A., Carvalho, T. D., Valenti, V. E., Crocetta, T. B., et al. (2016). Heart rate variability measure in breast cancer patients and survivors: a systematic review. Psychoneuroendocrinology 68, 57–68. doi: 10.1016/j.psyneuen.2016.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Arab, C., Vanderlei, L. C. M., da Silva Paiva, L., Fulghum, K. L., Fristachi, C. E., Nazario, A. C. P., et al. (2018). Cardiac autonomic modulation impairments in advanced breast cancer patients. Clin. Res. Cardiol. 107, 924–936. doi: 10.1007/s00392-018-1264-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Bettermann, H., Kröz, M., Girke, M., and Heckmann, C. (2001). Heart rate dynamics and cardiorespiratory coordination in diabetic and breast cancer patients. Clin. Physiol. 21, 411–420. doi: 10.1046/j.1365-2281.2001.00342.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Billman, G. E. (2011). Heart rate variability–a historical perspective. Front. Physiol. 2:86. doi: 10.3389/fphys.2011.00086

PubMed Abstract | CrossRef Full Text | Google Scholar

Billman, G. E. (2013). The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front. Physiol. 4:26. doi: 10.3389/fphys.2013.00026

PubMed Abstract | CrossRef Full Text | Google Scholar

Camm, A. J., Malik, M., Bigger, J. T., Breithardt, G., Cerutti, S., Cohen, R. J., et al. (1996). 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. Circulation 93, 1043–1065. doi: 10.1161/01.CIR.93.5.1043

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiang, J. K., Koo, M., Kuo, T. B. J., and Fu, C. H. (2010). Association between cardiovascular autonomic functions and time to death in patients with terminal hepatocellular carcinoma. J. Pain Symptom Manag. 39, 673–679. doi: 10.1016/j.jpainsymman.2009.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

De Couck, M., Caers, R., Spiegel, D., and Gidron, Y. (2018). The role of the vagus nerve in cancer prognosis: a systematic and a comprehensive review. J. Oncol. 2018, 1–11. doi: 10.1155/2018/1236787

PubMed Abstract | CrossRef Full Text | Google Scholar

De Couck, M., and Gidron, Y. (2013). Norms of vagal nerve activity, indexed by heart rate variability, in cancer patients. Cancer Epidemiol. 37, 737–741. doi: 10.1016/j.canep.2013.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

De Couck, M., Maréchal, R., Moorthamers, S., Laethem, J.-L. V., and Gidron, Y. (2016). Vagal nerve activity predicts overall survival in metastatic pancreatic cancer, mediated by inflammation. Cancer Epidemiol. 40, 47–51. doi: 10.1016/j.canep.2015.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

De Couck, M., Mravec, B., and Gidron, Y. (2012). You may need the vagus nerve to understand pathophysiology and to treat diseases. Clin. Sci. 122, 323–328. doi: 10.1042/cs20110299

PubMed Abstract | CrossRef Full Text | Google Scholar

de Godoy, M. F. (2016). Nonlinear analysis of heart rate variability: a comprehensive review. J. Cardiol. Ther. 3, 528–533. doi: 10.17554/j.issn.2309-6861.2016.03.101-4

CrossRef Full Text | Google Scholar

Giese-Davis, J., Wilhelm, F. H., Tamagawa, R., Palesh, O., Neri, E., Taylor, C. B., et al. (2015). Higher vagal activity as related to survival in patients with advanced breast cancer. Psychosom. Med. 77, 346–355. doi: 10.1097/PSY.0000000000000167

PubMed Abstract | CrossRef Full Text | Google Scholar

Gradishar, W. J., Anderson, B. O., Abraham, J., Aft, R., Agnese, D., Allison, K. H., et al. (2020). Breast cancer, version 3.2020, NCCN clinical practice guidelines in oncology. J. Natl. Compr. Cancer Netw. 18, 452–478. doi: 10.6004/jnccn.2020.0016

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernando, A., Lazaro, J., Gil, E., Arza, A., Garzon, J. M., Lopez-Anton, R., et al. (2016). Inclusion of respiratory frequency information in heart rate variability analysis for stress assessment. IEEE J. Biomed. Health Inform. 20, 1016–1025. doi: 10.1109/JBHI.2016.2553578

PubMed Abstract | CrossRef Full Text | Google Scholar

Karvinen, K. H., Murray, N. P., Arastu, H., and Allison, R. R. (2013). Stress reactivity, health behaviors, and compliance to medical care in breast cancer survivors. Oncol. Nurs. Forum 40, 149–156. doi: 10.1188/13.ONF.149-156

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D. H., Kim, J. A., Choi, Y. S., Kim, S. H., Lee, J. Y., and Kim, Y. E. (2010). Heart rate variability and length of survival in hospice cancer patients. J. Korean Med. Sci. 25, 1140–1145. doi: 10.3346/jkms.2010.25.8.1140

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuo, T. B., Lai, C. J., Huang, Y. T., and Yang, C. C. (2005). Regression analysis between heart rate variability and baroreflex-related vagus nerve activity in rats. J. Cardiovasc. Electrophysiol. 16, 864–869. doi: 10.1111/j.1540-8167.2005.40656.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Maki, A., Kono, H., Gupta, M., Asakawa, M., Suzuki, T., Matsuda, M., et al. (2006). Predictive power of biomarkers of oxidative stress and inflammation in patients with hepatitis c virus-associated hepatocellular carcinoma. Ann. Surg. Oncol. 14, 1182–1190. doi: 10.1245/s10434-006-9049-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Mäkikallio, T. H., Huikuri, H. V., Mäkikallio, A., Sourander, L. B., Mitrani, R. D., Castellanos, A., et al. (2001). Prediction of sudden cardiac death by fractal analysis of heart rate variability in elderly subjects. J. Am. Coll. Cardiol. 37, 1395–1402. doi: 10.1016/s0735-1097(01)01171-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Malliani, A., Pagani, M., Lombardi, F., and Cerutti, S. (1991). Cardiovascular neural regulation explored in the frequency domain. Circulation 84, 482–492. doi: 10.1161/01.CIR.84.2.482

PubMed Abstract | CrossRef Full Text | Google Scholar

Montano, N., Ruscone, T. G., Porta, A., Lombardi, F., Pagani, M., and Malliani, A. (1994). Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. Circulation 90, 1826–1831. doi: 10.1161/01.CIR.90.4.1826

PubMed Abstract | CrossRef Full Text | Google Scholar

Moody, G., Mark, R., Zoccola, A., and Mantero, S. (1985). Derivation of respiratory signals from multi-lead ECGs. Comput. Cardiol. 12, 113–116.

Google Scholar

Mouton, C., Ronson, A., Razavi, D., Delhaye, F., Kupper, N., Paesmans, M., et al. (2012). The relationship between heart rate variability and time-course of carcinoembryonic antigen in colorectal cancer. Auton. Neurosci. 166, 96–99. doi: 10.1016/j.autneu.2011.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Niederer, D., Vogt, L., Thiel, C., Schmidt, K., Bernhörster, M., Lungwitz, A., et al. (2013). Exercise effects on HRV in cancer patients. Int. J. Sports Med. 34, 68–73. doi: 10.1055/s-0032-1314816

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohira, H., Matsunaga, M., Osumi, T., Fukuyama, S., Shinoda, J., Yamada, J., et al. (2013). Vagal nerve activity as a moderator of brain–immune relationships. J. Neuroimmunol. 260, 28–36. doi: 10.1016/j.jneuroim.2013.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Palma, M. R., Vanderlei, L. C. M., Ribeiro, F. E., Mantovani, A. M., Christofaro, D. G. D., and Fregonesi, C. E. P. T. (2016). The relationship between post-operative time and cardiac autonomic modulation in breast cancer survivors. Int. J. Cardiol. 224, 360–365. doi: 10.1016/j.ijcard.2016.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, J., and Tompkins, W. J. (1985). A real-time QRS detection algorithm. IEEE Trans. Biomed. Eng. 32, 230–236. doi: 10.1109/TBME.1985.325532

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, B., Wang, L., Yan, C., Chen, D., Liu, M., and Li, P. (2019). Nonlinear heart rate variability biomarkers for gastric cancer severity: a pilot study. Sci. Rep. 9:13833. doi: 10.1038/s41598-019-50358-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, P. K., Domitrovich, P. P., Huikuri, H. V., and Kleiger, R. E. (2005). Traditional and nonlinear heart rate variability are each independently associated with mortality after myocardial infarction. J. Cardiovasc. Electrophysiol. 16, 13–20. doi: 10.1046/j.1540-8167.2005.04358.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

Tracey, K. J. (2009). Reflex control of immunity. Nat. Rev. Immunol. 9, 418–428. doi: 10.1038/nri2566

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsutsumi, T., Ide, T., Yamato, M., Kudou, W., Andou, M., Hirooka, Y., et al. (2007). Modulation of the myocardial redox state by vagal nerve stimulation after experimental myocardial infarction. Cardiovasc. Res. 77, 713–721. doi: 10.1093/cvr/cvm092

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanderlei, L. C. M., Pastre, C. M., Hoshi, R. A., Carvalho, T. D., and Godoy, M. F. (2009). Basic notions of heart rate variability and its clinical applicability. Rev. Bras. Cir. Cardiovasc. 24, 205–217. doi: 10.1590/S0102-76382009000200018

PubMed Abstract | CrossRef Full Text | Google Scholar

Varon, C., Lazaro, J., Bolea, J., Hernando, A., Aguilo, J., Gil, E., et al. (2019). Unconstrained estimation of HRV indices after removing respiratory influences from heart rate. IEEE J. Biomed. Health Inform. 23, 2386–2397. doi: 10.1109/JBHI.2018.2884644

PubMed Abstract | CrossRef Full Text | Google Scholar

Voss, A., Schulz, S., Schroeder, R., Baumert, M., and Caminal, P. (2009). Methods derived from nonlinear dynamics for analysing heart rate variability. Philos. Trans. A Math. Phys. Eng. Sci. 367, 277–296. doi: 10.1098/rsta.2008.0232

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. M., Wu, H. T., Huang, E. Y., Kou, Y. R., and Hseu, S. S. (2013). Heart rate variability is associated with survival in patients with brain metastasis: a preliminary report. Biomed. Res. Int. 2013, 1–6. doi: 10.1155/2013/503421

PubMed Abstract | CrossRef Full Text | Google Scholar