Abstract
Background:
Poor healing is an important contributing factor to amputation among patients with diabetic foot ulcers (DFUs). Physiological stress may slow wound healing and increase susceptibility to infection.
Objectives:
The objective was to examine the association between heart rate variability (HRV) as an indicator of physiological stress response and healing speed (HealSpeed) among outpatients with active DFUs.
Design and Methods:
Ambulatory patients with diabetes with DFUs (n = 25, age: 59.3 ± 8.3 years) were recruited. HRV during pre–wound dressing was measured using a wearable sensor attached to participants’ chest. HRVs were quantified in both time and frequency domains to assess physiological stress response and vagal tone (relaxation). Change in wound size between two consecutive visits was used to estimate HealSpeed. Participants were then categorized into slow healing and fast healing groups. Between the two groups, comparisons were performed for demographic, clinical, and HRV derived parameters. Associations between different descriptors of HRV and HealSpeed were also assessed.
Results:
HealSpeed was significantly correlated with both vagal tone (
Conclusion:
This study confirms an association between stress/vagal tone and wound healing in patients with DFUs. In particular, it highlights the importance of vagal tone (relaxation) in expediting wound healing. It also demonstrates the feasibility of assessing physiological stress responses using wearable technology in outpatient clinic during routine clinic visits.
Diabetic foot ulcers (DFUs), a leading cause of approximately 80% of lower-limb amputations, 1 develop in at least 25% of patients with diabetes 2 largely due to diabetic peripheral neuropathy. 3 In the US adult population reported estimates of the prevalence of diagnosed diabetes have increased from 5.9% in 2000 4 to 6.4% (~28.5 million) in 2010. 5 Whereas in 2005 a limb was lost to diabetes every 30 seconds, the frequency of limb loss has increased to every 20 seconds in 2011.5,6 Chronic DFUs also dramatically compromise physical activity pattern which is significantly associated with falls and can lead to serious complications, such as fractures, hospital admissions and poor quality of life.7,8 The economic costs associated with diabetic foot care, including amputation care, represent the single largest category of excess medical expenses related to diabetes. In the United States, the total annual cost for diabetic limb complications is estimated to be $17 billion, which exceeds the annual costs of expenses related to breast cancer and to colorectal cancer. 9
To avoid the adverse outcomes associated with DFUs and reduce the risk of amputation, it is necessary to expedite the wound healing process. Wound healing is a dynamic process, 10 and multiple factors including local (e.g. oxygenation, infection, etc) and systematic (e.g. demographics, nutrition, medication, stress, etc). Both local and systematic factors can impair the process of wound healing.11-14
Physiological stress is a systematic response to a stressor that facilitates adaptation to meet the challenge. The autonomic responses involved in modulating physiological stress include activation of the sympathetic and parasympathetic nervous systems, which work in tandem to keep the body in a state of homeostasis. During a stressful event, the sympathetic nervous system predominates, resulting in the fight-or-flight response. Because the body cannot maintain this state for extended periods of time, the parasympathetic system returns the body’s physiological conditions to the normal rest-and-digest state. Although the sympathetic physiological response is essential to protect the body and adapt to the stressors, prolonged exposure to stressors, referred to as episodic acute stress, can have an adverse impact on both psychological and physiological health and may affect the wound healing process. 15
Wound healing is a complex and fragile process, and stress may interrupt or lead to the formation of nonhealing chronic wounds in DFU patients. Stress disrupts the wound healing process primarily by mediating the hypothalamic-pituitary-adrenal and sympathetic-adrenal medullary axes as well as the psychological response, inducing unhealthy behaviors such as sedentary lifestyle, smoking, and so on. 11 Several studies confirmed the association between poor wound healing and stress by quantifying the presence of physiological stress.11-14
Several approaches have been proposed to quantify the presence of stress. These approaches include use of both subjective and objective tools. Questionnaires are a subjective method used to measure the patient’s perceived stress or mental strain. Objective methods record measurements induced stress via vital signs or biomarkers, such as cortisol sampling, heart rate, heart rate variability (HRV), blood pressure, sweating, and galvanic skin response.13,16-22
The effect of physiological stress on the wound healing process has not been explored specifically within a population with diabetes at high risk of limb amputation. However, studies have demonstrated that management of stress has a positive effect in diabetes patients. 23 Studies have also measured pain and stress during dressing changes in patients with a variety of chronic wounds using objective measures such as heart rate, blood pressure, and respiration rate. It has been suggested that wound dressing is a painful and stressful condition. Applying interdisciplinary approaches may help to reduce adverse outcomes such as pain and lower limb amputations in DFUs. 9
Advances in technology allow assessing phsyiological stress response by wearable wireless sensors. These sensors have few important advantages over wired alternatives including ease of use, comfort, enhanced freedom of mobility, and long term monitoring irrespective of setting (e.g. in home, in clinic, or in research laboratory).24-28 Our previous research assessed the level of physiological stress response in patients with DFUs during outpatient clinical visits and found the existence of high stress during wound dressing and wound care (e.g debridement, dressing change, etc). 29 In this exploratory study of patients with DFUs, we examined the physiological stress response in an outpatient wound care setting to determine the association between wound healing speed and phsyiological stress responses. This was done through quantification of HRV in time and frequency domains, which allow quantification of sympathetic regulation (indicator of stress) and vagal tone activity (indicator of relaxation).
Method
Participants’ Recruitment and Instruments
A convenience sample of ambulatory patients with non-infected and non-ischemic DFUs was recruited from The Diabetic Foot and Wound Center at the Hamad Medical Corporation, Doha, Qatar. Informed consent in agreement with the procedure of Hamad Medical Corporation was collected from each participant. The existence of neuropathy was confirmed via medical records and by Semmes-Weinstein Monofilament testing. Real-time ECG and heart rate data were collected using a body-worn device, BioHarness3 (Zephyr Technology Corp, Annapolis, MD) using the procedure described in our previous study. 29 Briefly, the device was attached to the chest using two electrocardiography (ECG) patches, and it recorded uni-channel ECG (250 Hz), respiration speed, accelerations and approximate core body temperature. ECG data were collected for approximately 10-15 minutes during the waiting period before wound dressing. The time was logged by a trained research assistant to synchronize with the bio-patch clock. The sensor allows for recording ECG during the entire duration of outpatient clinic visit—from the waiting room period (pre–wound dressing change [preWD]) to the completion of the wound care process, including debridement and wound dressing. The ECG measurements were recorded only during the first visit, but wound healing outcomes were assessed at baseline and the next clinic visits. Specifically, the wound area at the first and the second clinical visit were measured (21 ± 4 days apart) and the percentage of change was used to determine the wound healing speed.
Clinical Assessment
Along with demographic information, several questionnaires were used to assess participants’ characteristics including self-reported Mobility-Tiredness (Mob-T) Scale, 30 Barthel index, 31 depression scale (CES-D), 32 perceived stress scale, 33 Falls Self-Efficacy Scale (FES-I), 34 Short Health Survey (SF12), 35 and 0-10 Numeric Pain Rating Scale. 36 In addition, we collected the glycated hemoglobin or HbA1c level.
HRV Analysis and Quantification of Physiological Stress Response
The recorded sensor ECG data were transferred to a computer and fed into an open source software, Kubios, for HRV analysis. 37 The R-waves were extracted from the ECG signal and manually inspected for correction to estimate R-wave-to R-wave (NN) intervals. The time and frequency domain parameters were obtained for each nonoverlapping five-minute windows following the recommendations of the Task Force of the European Society of Cardiology and the North American Society of Pacing Electrophysiology. 38 Average NN intervals and average heart rate (HR-wave) were reported. The standard deviation of NN intervals (SDNN), which represents the variations in the R-to-R waves, was calculated with lower SDNN indicating the presence of physiological stress. Root mean squared of successive NN intervals (RMSSD), another time domain parameter commonly used as an index of vagally (vagal nerve) mediated cardiac control (relaxation), was calculated. 38 To measure frequency domain HRV parameters, we used the recommendation provided by the European Society of Cardiology guideline. 38 In summary, the power spectrum density (PSD) of time series representing R-wave-to-R-wave intervals was calculated. The frequency band of the PSD was divided into four bands: very low frequency (VLF, 0-0.04 Hz), low frequency (LF, 0.04-0.15 Hz), and high frequency (HF, 0.15-0.4 Hz). Then, the area under the curve for each frequency band, was calculated to estimate the energy localized at each frequency band. The normalized energy value at each frequency band was calculated using the following formulas: normalized LF (nLF) = LF/(LF+HF) × 100, normalized HF (nHF) = HF/(LF+HF)x100, and LF-to-HF ratio (LF/HF) = LF/HF. 38 The nLF represents the sympathetic regulation (stress indicator). The nHF represents the parasympathetic regulations (relaxation indicator) of the heart. LF/HF ratio, which indicates the balance between the sympathetic and parasympathetic activity of heart, was also estimated with higher values representing a more stressful condition. 38
Statistical Analysis
Correlations between various measured parameters and HRV were determined using Pearson correlation coefficients (
Based on the distribution of the “speed of wound healing” variable, the median split (–46.7%; the negative value indicates reduction in wound size in the follow-up visit)
13
was chosen as a cutoff for categorizing participants into two groups: fast healing and slow healing. Independent
Results
A total of 25 eligible participants recruited and completed baseline and follow-up visits. However, the data from 19 subjects (age = 59.3 ± 8.3 years, 83% male), who had valid recorded data included in the final analysis. The reasons for excluding six subjects were, observed arrhythmias and/or low signal quality (e.g. sensor frequent disconnections, poor sensor attachment, etc), which make difficult to accurantly estimate HRV (retention = 76%). Ten participants classified as slow-healers and 9 participants classified as fast-healers. Table 1 summarizes demographic and clinical data for each group. The speed of wound healing was 80% higher in the fast-healing group compared to the slow-healers (
Univariate Analysis of Descriptive Variables and Heart Rate Variability (HRV) Between Slow Healing and Fast Healing Groups (n = 19).
The results were adjusted by the wound size at the baseline.
When the correlation between speed of wound healing and HRV drived parameters was examined irrespective of group assignment, a strong correlation was observed between wound healing speed and nHF (indication of vagal tone or relaxation,

The association between the speed of wound healing and heart rate variability (HRV). Column A: the nHF, which represents the vagal tone (relaxation) is higher in individuals with better healing pace. Column B: the LF/HF ratio, which represents the sympathetic (stress) is lower in the individuals with better healing speed. The y-axis represented the speed of wound healing in percentage ( =
Among the HRV parameters, nLF (indicator of physiological stress response) was significantly lower on average by 41.7% in the fast healers than the slow healers (

Comparing the speed of wound healing in the two groups: patients with fast healing speed and patients with slow healing speed. The univariate analysis shows that several heart rate variability parameters are independent predictors between the two groups such as normalized low frequency (nLF), normalized high frequency (nHF), and the ratio of low frequency to high frequency. *Difference is statistically significant (
The predictive powers of demographic, clinical, and HRV parameters were assessed using discriminant analysis. Model 1 included demographics as well as patient’s level of HbA1c and years of diabetes, as potential independent predictors of wound healing rate, as suggested in the previous study
43
(Figure 3). However, test of equality found both parameters to be nonsignificant (

Receiver operating curve (ROC) of different models for predicting slow wound healing speed. Model 1. Demographic information of participants, ie, years in diabetes and HbA1c level. Model 2. The HRV parameters from the pre–wound dressing (preWD) interval.
Comparison of Logistic Models.
Demographic parameters, model using years in diabetes and A1c as input parameters; HRV parameters, model using HRV parameters.
Model 2 included nLF, nHF, and LF/HF. Test of equality found all three parameters to be significant (
Discussion
The study assessed the association between physiological stress responses and wound healing outcomes among patients with diabetes who had plantar ulcers. While prior studies have demonstrated a negative association between psychological stress and acute wound healing processes 44 or chronic wound healing process, 45 little is known about the association between physiological stress response and chronic wound healing in patients with DFUs. This study proposed heart rate variability (HRV) measurement to objectively quantify physiological stress responses and examined its association with wound healing outcomes in diabetes patients with DFUs.
We did not observe an association between subjective stress assessment (eg, questionnaires) and wound healing speed. In contrast, objective stress assessment of physiological stress showed significant correlation with wound healing speed. Subjective tools like questionnaires assess patients’ emotional and psychological responses. 20 Thus, questionnaires can sometimes cause patients to report high levels of stress to meet the expectations of the health care assistants. 20 Similarly, patients may be accustomed to answering questions in a certain way, such as selecting the middle options on a rating scale rather than the extreme responses.20,46 Miller et al. 47 examined more than 300 empirical articles describing a relationship between psychological stress and parameters of the immune system in human participants and revealed that while chronic stressors are associated with suppression of both cellular and humoral measures, subjective reports of stress generally did not associate with immune changes, highlighting limitation of subjective assessment of stress for chronic wound healing studies. A newly published study confirmed a lack of efficacy in subjective stress assessments to indicate the level of physiological stress in the participants. 48 Thus, it has been reported that no correlation existed between subjective and objective metrics of stress; therefore, the self-stress evaluation questionnaires may not adequately capture the existence of physiological stress.
Several studies indicated that wound pain is a common cause of distress in patients with acute or chronic wounds.20,49 In patients with diabetes who had foot ulcers, wounds developed because of loss of protective sensation (diabetic peripheral neuropathy). Therefore, the previous conceptual model for the chronic wound which claimed wound pain as a primary source of discomfort may not be applicable in DFU patients. In this study, 67% of participants had mild pain (pain score ≤ 3.4) and no significant association between wound healing and pain observed. However, there is a need for further investment in DFU patients to build a better conceptual model which represents the underlying physiological process related to wound healing process and pain. One of the physiological factors is poor sleep, which has been shown to be underdiagnosed in the patient with chronic wounds. 50 Poor sleep reduces the vagal tone/relaxation and increases the sympathetic/stress activation. Such a process may cause vasoconstriction which decreases perfusion in the wound area 51 or may cause higher cortisol level in blood which increases blood sugar (HbA1c) and slows the wound healing process. 51 In addition, sleep apnea, with its high prevalence in patients with diabetes, 50 reduces the saturated oxygen level in the blood 50 and speed of wound healing. Further study is needed to examine whether management of sleep disorders/stress may assist in reducing stress/sleep quality and improving wound healing outcomes.
A negative correlation (
HRV parameters extracted form frequency domains (ie, nLF, nHl, and LF/HF) have achieved to significantly discriminate between slow or fast healers. While a trend was also observed for HRV parameters extracted from time-domain (i.e., SDDNN and RMSSD), these parameters didn’t achieve statistical significant level in our sample to discriminate between groups. This may suggest that HRV frequency derived parameters may be more sensitive to track stress/relaxation in people with diabetes and foot ulcers during short term monitoring (ie, pre–wound dressing interval). This finding is aligned with a previous study demonstrating that while there is an agreement between short-term (30 minutes) and long-term (24 hours) frequency domain HRV parameters, this agreement is diminished when time domain HRV parameters were examined. 56 This may suggest that for an accurate information extracted from time domain HRV parameters, a longer recording (ie, minimum 2 hours 57 ) may be needed.
This study suggests that frequency domain HRV parameters derived from a single chest worn sensor module yields a robust prediction model superior than a model based on demographic or clinical parameters (AUC > 0.9 for the model based on HRV vs AUC < 0.50 for a model based on demographic and clinical data; Figure 3). This finding is in line with a previous study 43 demonstrating that nonmodifiable demographic (ie, sex, age, years in diabetes) or clinical (baseline wound size) factors are unable to predict the outcome of wound healing in a single visit.
Limitations
In this exploratory study, we faced a few limitations. The sample size in this study was small and may not have sufficient power to confirm significant statistical correlation with parameters of interest. The results should be confirmed in higher powered study. The time between the second visit and the first visit was not strictly controlled, which may affect the estimation of wound healing speed in our sample. We focused only on the pre–wound dressing period when the patient was in the waiting room and prior visiting his or her doctor, assuming that this time may better represent overall physiological stress in the subject. Another study is warranted to confirm whether the pre–wound dressing period could represent daily physiological stress response in diabetes patients with DFU. Despite these limitations, the results are encouraging and warrant future longitudinal studies with larger sample size to confirm the results.
Conclusions
The current study highlights the associations between physiological stress and speed of wound healing in patients with DFUs. The results also emphasize the potential of HRV and in particular the metrics representing the ability of relaxation to be a significant predictor of wound healing outcomes. The results of this exploratory study warrant further investigation in a large sample size. Our approach to measuring the association between physiological stress and speed of wound healing may open new avenues to better target slow healers and implement stress management strategies to speed up wound healing and prevent limb amputation in patients with diabetes.
Footnotes
Acknowledgements
We wish to thank Dr Saman Parvaneh, Ms Yuxi Xia, Mr Robert Menzies, and Mrs Samira Abdulla, who contributed in part to data analysis, data collection, or patient recruitment. The authors would like to thank Kimberly Macellaro, PhD, a member of the Baylor College of Medicine Michael E. DeBakey Department of Surgery Research Core Team, for her editorial assistance during the preparation of this manuscript.
Abbreviations
ANOVA, analysis of variance; AUC, area under the curve; BMI, body mass index; CES-D, Center for Epidemiologic Studies Depression; CI, confidence interval; DFU, diabetic foot ulcer; ECG, electrocardiography; FES-I, Falls Self-Efficacy Scale; HRV, heart rate variability; IRB, institutional review board; LF/HF, ratio of low frequency components to high frequency components of normalized R-to-R waves; nHF, normalized high frequency components of normalized R-to-R waves; nLF, normalized low frequency components of normalized R-to-R waves; RMSSD, root mean squared of successive normalized R-to-R wave; ROC, receiver operating characteristic; SD, standard deviation; SDNN, standard deviation of normalized R-to-R wave; SF12, Short Health Survey; SF12-MCS, Mental Health Composite Scale score of SF12; SF12-PCS, Physical Health Composite Scale score of SF12.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The project described was supported in part by Award Number NPRP 4-1026-3-277 from the Qatar National Research Foundation. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Qatar National Research Foundation.
