Early bereavement is associated with increased cardiovascular events. The mechanism, however, has not been well studied. We assessed whether bereavement is associated with an increased heart rate (HR) and decreased heart rate variability that might contribute to increased cardiovascular risk. A total of 78 bereaved spouses and parents (55 women and 23 men; aged 34 to 87 years, mean 65) were studied with 24-hour Holter monitoring within 2 weeks of bereavement (acute) and at 6 months. Their findings were compared to those from a nonbereaved reference group (52 women and 27 men) aged 33 to 91 years (mean 63.6). All participants were in sinus rhythm. We assessed the mean HR, atrial and ventricular arrhythmias, and both time and frequency domain heart rate variability measures. Acute bereavement was associated with increased 24-hour HR (mean ± SE, 75.1 ± 1.1 vs 70.7 ± 1.0; p = 0.004) and reduced heart rate variability, as indicated by lower standard deviation of the NN intervals index (median 45.4 vs 49.9, p = 0.017), total power (7.78 ± 0.10 vs 8.02 ± 0.09, p = 0.03), very low frequency (7.23 ± 0.09 vs 7.44, p = 0.046) and low frequency (5.76 ± 0.12 vs 6.16 ± 0.09, p = 0.01). At 6 months, the bereaved had a significantly lower HR (p = 0.001) and increased standard deviation of the NN intervals index (p = 0.02), square root of the mean square of differences of successive intervals (p = 0.045), number of interval differences of successive NN intervals >50 ms divided by the number of NN intervals (p = 0.039), low-frequency power (p = 0.02), and high frequency (p = 0.002) compared to the initial acute levels. In conclusion, the present study, the first to report 24-hour HR monitoring in the early weeks of bereavement, has demonstrated increased HR and altered autonomic function that might contribute to the increased cardiovascular events in early bereavement.
Bereavement is associated with increased morbidity and mortality, particularly in the surviving spouse and parents. Cardiovascular disease (CVD) accounts for 20% to 53% of excess deaths during spousal bereavement, with the greatest risk in the immediate weeks after the bereavement but remaining elevated for the first 6 months. The acute effect is exemplified by a recent analysis in which bereavement increased the relative risk of acute myocardial infarction by 21-fold in the first 24 hours and 4.8-fold at 1 month. The autonomic changes, elevated heart rate (HR), and arrhythmia are potential contributors to the increased CVD risk. An increased HR predicts arrhythmias, sudden cardiac death, and CVD mortality, even after controlling for traditional risk factors. Lower heart rate variability (HRV) is also predictive of increased CVD risk, with increased sympathetic activity associated with increased inflammation, blood pressure, and platelet reactivity. To date, the effects of bereavement on autonomic function have not been well studied. We evaluated the 24-hour HR and HRV measurements in the first 2 weeks of bereavement and compared these findings to those 6 months later, a period of lower CVD risk. We also recruited a sample of nonbereaved controls as a reference group.
Methods
The present study was a prospective, controlled cohort study. From February 2005 through June 2008, we recruited bereaved spouses or partners or parents (74 spouses and 6 parents) of deceased patients from critical care units in 5 urban hospitals. The bereaved subjects were recruited and studied within 2 weeks of bereavement and at 6 months. Concurrently, we enrolled 80 nonbereaved relatives of patients using the same hospitals as the bereaved. Of the participants, 3 (2 bereaved and 1 nonbereaved) were excluded from the present analysis because of underlying atrial fibrillation. Therefore, 78 bereaved and 79 nonbereaved participants were included in the final analysis. The institutional ethics committee approved the study protocol.
The bereaved spouses or parents who met the eligibility criteria at the participating institutions were contacted by a social worker or study investigator at the hospital or by telephone within 72 hours after the bereavement. Interested subjects were assessed at home within 2 weeks of their bereavement (average 11 days, range 3 to 15). Of those who met the inclusion criteria, 60% agreed to participate; 3 bereaved (4%) withdrew from the study before the 6-month assessment and 1 died before the final assessment from multiorgan failure. The control group consisted of family members of patients, who were recruited from clinical areas (spouses or parents) such as day surgery, orthopedic, or endoscopy clinics of the same hospitals as the bereaved. The initial evaluation of nonbereaved participants was approximately 1 month after hospital discharge of their hospitalized family member.
The bereaved and nonbereaved participants were excluded if they reported severe illness (ie, respiratory, heart, liver, renal failure), malignancy, history of coagulopathy, cognitive impairment, psychosis, or immunosuppressive illness, were a nursing home resident, were unable to speak or read English, or had baseline rhythm of atrial fibrillation. We excluded nonbereaved subjects who had experienced the death of a close family member or friend in the previous 2 years. Only 1 bereaved participant had experienced a second close bereavement in the previous 2 years.
The symptoms of depression were measured using the Center for Epidemiologic Studies Depression (CES-D) questionnaire, and the symptoms of anxiety and anger were measured using the Spielberger State Anxiety and Anger scales. Social support was also measured using a validated questionnaire (Social Support Questionnaire-6), and the morning blood cortisol level was collected and analyzed, as described previously. The participants were instructed to undertake their usual daily activity and completed an activity diary every 30 minutes throughout the waking hours of the monitoring period. The activity level was documented in the diary on a 1 to 10 Likert scale with greater levels, indicating more activity.
Five-lead Medilog AR4 Holter monitors (Schiller, Australia) were used for 24-hour monitoring and analysis conducted using Medilog Prima 3 software (Schiller, Australia). All RR interval time series were first edited automatically, after which careful manual editing was performed by visual inspection. All abnormal beats were manually checked and arrhythmia strips classified by a cardiologist who was unaware of the study group. After editing of the RR intervals, the RR interval spectrum was computed for the entire recording interval. A Fast Fourier transform algorithm was used to estimate the power spectrum densities of HR variability. The frequency domain measures of RR interval variability were computed by integrating the point power spectrum over the frequency intervals. The power spectra were quantified by measuring the areas in the following frequency bands: (1) overall band (0.00 to 0.40 Hz range), (2) ultra-low-frequency power (<0.0033 Hz), (3) very-low-frequency power (0.0033 to <0.04 Hz), (4) low-frequency (LF) power (0.04 to 0.15 Hz), and (5) high-frequency (HF) power (0.15 to <0.40 Hz).
The time domain measures included the standard deviation of the NN intervals index (SDNNi), an estimate of the changes in HR due to cycles <5 minutes and sensitive to short-term changes in HRV; the square root of the mean squared differences of successive intervals (rMMSD), and number of interval differences of successive NN intervals >50 ms divided by the number of NN intervals (pNN50), which similar to the rMMSD, is associated with parasympathetic modulation of the HR.
Differences in the sociodemographic data were analyzed using Student’s t test for interval data and the chi-square test for categorical data. The primary analysis comparing the bereaved to nonbereaved at the acute assessment used the t test or Mann–Whitney U test, as appropriate. The longitudinal changes from acute to 6 months were evaluated using paired analysis (paired t test or McNemar test, as appropriate). The skewed HRV power parameters were logarithmically transformed before statistical testing. To explore the relations between HR and HRV variables with the measure of psychological stress and cortisol, Spearman’s rank correlation and Pearson’s correlation tests were used. The data were analyzed using SPSS for Windows, version 17.0 (SPSS, Chicago, Illinois).
Results
The bereaved and nonbereaved had similar sociodemographic characteristics, including the level of social support ( Table 1 ). Most (n = 72) were spouses or partners of the deceased, and 6 were parents. For bereaved spouses, the number of years as a couple ranged from 3 to 62 years (mean 38). The cause of death for the deceased family members of bereaved participants was CVD in 21 (cardiac arrest in 17, heart failure in 2, aortic aneurysm in 2), respiratory failure in 17, cerebral hemorrhage in 13, multiorgan failure in 9, cancer in 9, liver disease in 4, trauma in 2, suicide in 2, and meningitis in 1.
Variable | Bereaved | Nonbereaved | p Value |
---|---|---|---|
(n = 78) | (n = 79) | ||
Age (years) | |||
Mean | 65.0 | 63.5 | 0.36 |
Range | 34–87 | 33–91 | |
Women | 52 (69%) | 52 (66%) | 0.43 |
Living alone (at assessment) | 45 (58%) | 6 (8%) | <0.001 |
Education (highest level) | |||
High school | 31 (40%) | 28 (36%) | |
Technical | 26 (33%) | 20 (25%) | |
University | 18 (23%) | 28 (36%) | |
Did not state | 3 (4%) | 3 (4%) | 0.24 |
Household income | |||
<$20,000 | 16 (21%) | 12 (15%) | |
$20,001–$60,000 | 20 (26%) | 26 (34%) | |
>$60,000 | 24 (30%) | 28 (36%) | |
Did not state | 18 (23%) | 13 (15%) | 0.43 |
Body mass index (kg/m 2 ) | 26.1 ± 5.0 | 27.1 ± 5.0 | 0.14 |
Waist circumference | 93.3 ± 12.7 | 94.9 ± 15.8 | 0.48 |
History of myocardial infarction or angina | 13 (17%) | 4 (5%) | 0.11 |
Diabetes mellitus | 4 (5%) | 6 (8%) | 0.22 |
Stroke | 5 (6%) | 4 (5%) | 0.59 |
Peripheral vascular disease | 3 (4%) | 1 (1%) | 0.22 |
Hypertension (by history) | 32 (41%) | 30 (38%) | 0.11 |
Heart rate-lowering drugs | 8 (11%) | 3 (4%) | 0.10 |
Depression symptoms | 25.9 ± 11.8 | 6.15 ± 5.1 | <0.001 |
State anxiety | 46.4 ± 13.9 | 28.9 ± 8.0 | <0.001 |
Symptoms of anger | |||
Median | 16.0 | 15.0 | <0.001 |
Interquartile range | 16–37 | 15–15 ⁎ | |
Social support | |||
Availability (possible range 0–9) | 4.7 ± 2.5 | 4.4 ± 8.0 | 0.43 |
Satisfaction (possible range 1–6) | 5.2 ± 1.1 | 5.2 ± 1.1 | 0.87 |
⁎ Nonbereaved (88%) reported a score of 15; therefore, no interquartile range.
The level of physical activity, as recorded during 24-hour monitoring, was similar between the bereaved and nonbereaved at the initial assessment (mean ± SE 5.1 ± 0.2 vs 4.9 ± 0.2, p = 0.23; possible range 1 to 10). Among the bereaved, no difference was found in the activity level from acute assessment to 6 months (5.1 ± 0.2 vs 5.3 ± 0.2; p = 0.33). HR-lowering medication was used by 8 of the bereaved (β blockers in 3 and calcium channel blockers in 5) at baseline and by 3 of the nonbereaved controls (β blocker in 1 and calcium channel blockers in 2; p = 0.16).
The symptoms of depression, anxiety, and anger were all elevated acutely in the bereaved compared to the nonbereaved and reduced in the bereaved at 6 months compared to at the acute assessment ( Table 1 ). Additionally, the morning cortisol levels were greater in the bereaved participants than in the nonbereaved acutely (median 312 mmol/L, interquartile range 241 to 414; vs 271 mmol/L, interquartile range 222 to 353; p = 0.01). The cortisol levels did not change from the acute to the 6-month assessment in the bereaved (median 326 mmol/L, interquartile range 237 to 400; p = 0.47) or the nonbereaved (median 279 mmol/L, interquartile range 229 to 361; p = 0.39).
The results of 24-hour HR analysis are listed in Table 2 . At the initial acute assessment, the bereaved had a greater 24-hour average and minimum HR and more atrial premature beats and episodes of nonsustained (>3 consecutive beats) supraventricular tachycardia (SVT) compared to the nonbereaved ( Table 2 ). Similar between-group differences were observed after excluding the 8 bereaved and 3 nonbereaved participants taking HR-lowering medications (data not shown). At 6 months, the average and minimal HR and episodes of SVT were all significantly lower in the bereaved participants compared to their initial acute assessment ( Table 2 ). However, there were no differences in HR or arrhythmia in the nonbereaved from the acute assessment to the 6-month assessment ( Table 2 ). There were too few ventricular premature complexes or episodes of ventricular arrhythmias in either bereaved or nonbereaved to enable statistical comparison.
Variable | Bereaved | Nonbereaved | p Value | ||||
---|---|---|---|---|---|---|---|
Acute | 6-mo | Acute | 6-mo | B vs NB Acute | B Acute to 6 mo | NB Acute to 6 mo | |
Heart rate | |||||||
Average | 75.1 ± 9.8 | 71.0 ± 9.7 | 70.7 ± 8.9 | 70.3 ± 8.5 | 0.004 | 0.001 | 0.10 |
Minimum | 50.1 ± 6.5 | 48.8 ± 6.5 | 48.0 ± 7.1 | 48.9 ± 6.0 | 0.039 | 0.06 | 0.65 |
Maximum | 135 ± 19 | 131 ± 20 | 132 ± 20 | 126 ± 18 | 0.39 | 0.05 | 0.006 |
APBs | |||||||
Median | 75 | 59 | 45.0 | 43.0 | 0.05 | 0.42 | 0.84 |
Interquartile range | 32–158 | 24–199 | 18–111 | 11–156 | |||
Episodes of nonsustained SVT | 2.2 ± 6.3 | 0.9 ± 13.8 | 1.2 ± 4.4 | 0.7 ± 1.3 | 0.002 | 0.002 | 0.61 |
VPBs | |||||||
Median | 13 | 12 | 7 | 5 | 0.74 | 0.39 | 0.27 |
Interquartile range | 2–131 | 2–107 | 1–221 | 1–83 | |||
Episodes of VT (>3 beats) | 1 subject with 1 episode | 2 subjects with 1 episode each | 2 subjects with 1 episode each | 2 subjects with 1 episode each | NA | NA | NA |
At the initial acute assessment, the bereaved had lower SDNNi, total power, very-low-frequency, and LF power than the nonbereaved ( Table 3 ). In the bereaved, the SDNNi, rMMSD, pNN50, and LF and HF power parameters were all greater at 6 months compared to the initial acute assessment. No significant differences were seen in any HRV measures from the initial assessment to the 6-month assessment in the nonbereaved participants ( Table 3 ).
Variable | Bereaved | Nonbereaved | p Value | ||||
---|---|---|---|---|---|---|---|
Acute | 6 mo | Acute | 6 mo | B vs NB, Acute | B, Acute to 6 mo | NB, Acute to 6 mo | |
SDNNi (ms) | 45.4 (35.8–61.7) | 46.0 (40.8–69.4) | 49.9 (42.8–59.3) | 49.8 (42.3–59.4) | 0.02 | 0.02 | 0.38 |
RMSSD (ms) | 30.1 (18.8–44.0) | 36.5 (23.9–53.5) | 26.8 (22.2–39.8) | 30.4 (19.9–38.8) | 0.15 | 0.04 | 0.74 |
PNN50 (ms) | 5.6 (1.6–12.0) | 7.2 (3.17–14.2) | 5.2 (2.4–9.9) | 5.5 (2.4–10.2) | 0.49 | 0.04 | 0.54 |
Total power (ms ) | 2,225 (1378–3578) | 2,700 (1892–6171) | 2,935 (1980–4063) | 2,891 (1950–3843) | 0.03 | 0.08 | 0.23 |
Ultra-low-frequency power (ms ) | 294 (169–524) | 320 (189–947) | 404 (212–608) | 262 (238–485) | 0.07 | 0.62 | 0.13 |
Very-low-frequency power (ms ) | 1347 (853–2,104) | 1489 (1,110–3,373) | 1677 (1,185–2,421) | 1705 (1,176–2,353) | 0.046 | 0.06 | 0.24 |
Low-frequency power (ms ) | 386 (128–582) | 521 (298–1,041) | 521 (295–753) | 544 (303–804) | 0.01 | 0.02 | 0.42 |
High-frequency power (ms ) | 152 (109–258) | 257 (138–573) | 177 (137–296) | 177 (130–326) | 0.06 | 0.002 | 0.89 |
LF/HF power ratio | 2.1 (1.6–2.9) | 2.1 (1.5–2.9) | 2.4 (1.8–3.9) | 2.3 (1.7–3.2) | 0.076 | 0.487 | 0.924 |