Does psychological treatment of major depression reduce cardiac risk biomarkers? An exploratory randomized controlled trial

Background Depression is associated with an increased risk for cardiovascular disease (CVD). Biological cardiac risk factors are already elevated in depressed patients without existing CVD. The purpose of this exploratory trial was to examine whether treating Major Depression (MD) with cognitive behavioral therapy (CBT) is associated with improvements in cardiac risk biomarkers and whether depressive symptom severity at baseline moderates treatment effects. Methods Eighty antidepressant-free patients with MD were randomly assigned to CBT or waiting list (WL). Biological outcomes included long-term recordings (24-h, daytime, nighttime) of heart rate, heart rate variability (HRV), and blood pressure, as well as inflammatory markers such as C-reactive protein (CRP), interleukin (IL)-6, and tumor necrosis factor (TNF)-α. A sample of 40 age- and sex-matched non-clinical controls was also involved to verify biological alterations in MD at study entry. Results Compared to WL, CBT was associated with a significant increase in overall HRV, as indexed by the 24-h and daytime HRV triangular index, as well as trend improvements in 24-h low-frequency HRV and daytime systolic blood pressure. Self-rated depressive symptom severity moderated (or tended to moderate) improvements in CBT for 24-h and daytime heart rate and several indices of HRV (especially daytime measures). Inflammatory treatment effects were not observed. Conclusions CBT increased overall HRV in patients with MD. Initially more depressed patients showed the most pronounced cardiovascular improvements through CBT. These exploratory findings may provide new insights into the biological effects of psychological treatment against depression and must be confirmed through future research.

Although considered the gold-standard for the psychological treatment of MD, one important and understudied question is, however, whether CBT reduces cardiac risk biomarkers in populations with MD without CVD. In terms of HRV, two randomized controlled trials (RCT) indicated that short CBT interventions (i.e. 4-6 weeks) combined with HVR biofeedback (Caldwell & Steffen, 2018) or breathing exercises (Chien, Chung, Yeh, & Lee, 2015) may improve specific measures of HRV. Both studies have important limitations, such as a small sample size of younger females (Caldwell & Steffen, 2018), a short treatment duration (Caldwell & Steffen, 2018;Chien et al., 2015), the inclusion of patients with antidepressant medication (Chien et al., 2015), as well as the use of short-term measures of HRV (Caldwell & Steffen, 2018;Chien et al., 2015) instead of long-term HRV recordings which are more stable, less affected by placebo effects and more predictive for cardiovascular events in populations without known CVD (Hillebrand et al., 2013; Task Force of The European Society of Cardiology and The North American & Society of Pacing and Electrophysiology, 1996). In terms of inflammation, one previous RCT suggests that CBT lowers IL-6 and TNF-α in a sample of young adults (Moreira et al., 2015), while two trials did not observe an overall effect of CBT on inflammatory markers (Euteneuer et al., 2017;Taylor et al., 2009).
The present RCT examined whether established CBT for depression is accompanied by reductions in cardiac risk biomarkers in patients with MD without CVD. We aimed to extend previous research in several regards. First, we exploratory studied biological effects of CBT for a wide range of measures with predictive value for CVD, such as heart rate (Greenland et al., 1999;Kannel, Kannel, Paffenbarger, & Cupples, 1987), indices of HRV (Hillebrand et al., 2013), blood pressure (Lewington, Clarke, Qizilbash, Peto, & Rory, 2002), as well as inflammatory markers such as CRP, IL-6 and TNF-α (Avan et al., 2018;Pearson, 2003;Subirana et al., 2018). Second, this RCT only included patients without antidepressant medication. This point is important because antidepressants, in particular tricyclic antidepressants, may reduce HRV (Kemp et al., 2016(Kemp et al., , 2010Kemp, Quintana, & Malhi, 2011;Yeh et al., 2016). Third, we used 24-h long-term recordings for cardiac measures instead of short-term recordings. Long-term recordings are more predictive for CVD (Hillebrand et al., 2013) and provide the possibility to analyze geometric indices of HRV such as the triangular index, which has been suggested for prospective research because of the highest reliability among all HRV indices (Ziegler, Piolot, Strassburger, Lambeck, & Dannehl, 1999). In addition, long-term recordings enable the possibility to consider circadian aspects of cardiac measures (Bilan, Witczak, Palusiński, Myśliński, & Hanzlik, 2005;Bonnemeier et al., 2003;Carney et al., 2000;Li et al., 2011;Ohkubo et al., 2002;Parati et al., 2014). Based on a pre-post study suggesting that CBT improved heart rate and daytime HRV in severely depressed patients but not in mildly depressed patients with coronary heart disease (Carney et al., 2000), we examined whether the intensity of depressive symptoms at study entry moderates potential effects on biological outcomes. This trial should be considered exploratory, because a wide range of biological markers were assessed, aiming to identify potential biological responsiveness for future research.

Materials and methods
This RCT was conducted from October 2015 to October 2019 in accordance with the World Medical Association Declaration of Helsinki and the ethical guidelines of the German Society for Psychology (see online Supplementary material for CONSORT 2010 checklist). The approval was given by the Institutional Review Board of the Department of Psychology at the University of Marburg (Approval number: 2014-26k). The study was funded by the German Research Foundation (DFG EU 154/2-1), and study registration took place (www.clinicaltrial.gov NCT 02787148).

Participants
Eighty antidepressant-free patients with MD according to DSM-IV were (1:1) urn randomized to 14 weeks of CBT or a waitlist control condition (see online Supplementary material for recruiting, setting and intervention). A sample of 40 (1:2) ageand sex-matched non-clinical controls from the same community was also involved to verify potential baseline alterations in biological markers in MD.

Inflammatory markers
Non-fasting blood samples were collected in EDTA-treated tubes (S-Monovette, Sarstedt, Nümbrecht, Germany) between 8:00 AM and 11:00 AM. Before each blood sampling, participants were queried about acute infections during the last 14 days, chronic infections, and illness. Plasma for CRP and cytokine measurements were separated by centrifugation at 2000 × g for 10 min at 4°C, and plasma was stored at −80°C until analysis. CRP and cytokines were analyzed using the V-PLEX Human CRP Kit and MSD cytokine assays (Meso Scale Discovery, Rockville, USA) according to the manufacturer's instructions. The sensitivity of the assays was 1.33 pg/ml for CRP, 0.04 pg/ml for TNF-α, and 0.06 pg/ml for IL-6.

Ambulatory cardiovascular monitoring
Twenty-four-hour oscillometric blood pressure and electrocardiogram (ECG) monitoring for heart rate and HRV analyses was conducted during weekdays using a hybrid blood pressure/ECG monitor (card(X)plore, Meditech Ltd., Hungary). Participants were instructed to follow their typical daily activities but refrain from intense physical activity. Monitoring started between 8 AM and 11 AM. Considering circadian variations in heart rate, blood pressure, and HRV (Bilan et al., 2005;Bonnemeier et al., 2003;Li et al., 2011), these measures were calculated for the total 24-h recording and additionally for day-and nighttime. Using a fixed-narrow time interval approach, measures for dayand nighttime were calculated from standardized periods (i.e. daytime: 9 AM to 9 PM; nighttime: 1 AM to 6 AM) in which the retiring and rising periods (which are subject to considerable variation) were eliminated to increase reliability and comparability (Fagard, Brguljan, Thijs, & Staessen, 1996;O'Brien et al., 2013). CardioVisions 1.24 (Meditech Ltd., Hungary) was used to obtain parameters for heart rate, HRV, and blood pressure. Heart rate was recorded continuously by a 3-channel Holter ECG at a sampling rate of 600 Hz using seven Ag/AgCl electrodes (Kendall H92SG, Cardinal Health Germany 507 GmbH, Germany) placed on the thoracic region. First, similar QRS complexes were automatically grouped and defined as normal or artifact beats. Next, QRS complexes were visually screened for incorrect beat detections and, in case of an incorrect reading, marked as artifact beat. All artifacts were excluded from the data. A minimum of 18 h of analyzable data and less than 5% artifacts of all analyzable beats were required for a recording to be accepted for further analyses (Crawford et al., 1999).
As recommended by common guidelines (Task Force of The European Society of Cardiology and The North American & Society of Pacing and Electrophysiology, 1996), the following timedomain HRV measures were assessed: the HRV triangular index and the standard deviation of all NN intervals (SDNN) for overall HRV, as well as the square root of the mean of the sum of the squares of differences between adjacent NN intervals (RMSSD) for estimating vagally mediated changes reflected in HRV. In terms of frequency-domain HRV, low-frequency (LF)-HRV (0.04-0.15 Hz) and high-frequency (HF)-HRV (0.15-0.4 Hz) were analyzed. HF-HRV reflects vagal modulation of HR. LF-HRV is more complex and may include both sympathetic and parasympathetic influences (Shaffer & Ginsberg, 2017; Task Force of The European Society of Cardiology and The North American & Society of Pacing and Electrophysiology, 1996). To determine LF and HF parameters, the software utilizes a Fast Fourier Transform (FFT) algorithm considering 4-min Hanning-windowed samples.
For ambulatory blood pressure monitoring (ABPM), the device was programmed to obtain blood pressure readings every 15 min during daytime (7 AM to 10 PM) and every 30 min during the night-time (10 PM to 7 AM). The inflatable cuff was placed on the participant's right arm. After the measurement, ABPM readings were manually screened for artifacts. As recommended, only strongly incorrect readings were deleted from recordings (O'Brien et al., 2003). Systolic blood pressure (SBP) values lower than 60 or higher than 260, and diastolic blood pressure (DBP) values lower than 40 and higher than 150 were defined as artificial readings and excluded. Daytime, nighttime and 24-h calculations based on weighted average blood pressure. Nocturnal SBP and DBP dipping was determined by subtracting the nighttime average blood pressure from the daytime average blood pressure. Thus, lower values on these different measures reflect lower nocturnal blood pressure dipping. This procedure yields adequate reliability for dipping values (Dimsdale et al., 2000).

Data analysis
Baseline measures are reported as means with SDs for continuous variables and as numbers with percentages for categorical variables. Data distributions were inspected, and screening for extreme outliers was conducted (i.e. values more than three interquartile ranges above quartile 3). Potential differences in study variables between patients with MD and non-clinical controls were examined using Welch's t tests and χ 2 tests. Constrained longitudinal data analysis (cLDA) with linear mixed models (Coffman, Edelman, & Woolson, 2016;Fitzmaurice, Laird, & Ware, 2011;Liang & Zeger, 2000;Liu, Lu, Mogg, Mallick, & Mehrotra, 2009) was applied to examine differences in changes from baseline to the end of treatment between both treatment groups. Constrained linear mixed models incorporated baseline adjustment (i.e. baseline means are constrained to be equal between the randomized groups;) with time (baseline, end of treatment) and time × treatment interactions as fixed factors and with a random intercept at patient level. Longitudinal data were analyzed on an intention-to-treat base using maximum likelihood estimation to account for missing data and dropouts, respectively. To examine moderating effects of depressive symptom severity, time × treatment × baseline severity interactions, as well as corresponding lower-order terms, were added to constrained linear mixed models. Following common suggestions for exploratory research (Althouse, 2016;Bender & Lange, 2001;Rothman, 1990;Rubin, 2017), we report and interpret results without adjusting for multiple testing. All p-values are two-tailed and analyses were carried out with SPSS version 20.0 for Windows (Chicago, SPSS, Inc.) and Mplus7 (Muthén & Muthén, 1998.
Descriptive statistics for all groups and comparisons between patients with MD and non-clinical controls are presented in Table 1 (see online Supplementary material Table S1 for correlations between biological variables at study entry and Table S2 for correlations between changes in biological variables for each treatment group). Compared to non-clinical controls, patients with MD exhibited significantly higher levels in all inflammatory markers (i.e. CRP, IL-6 and TNF-α). Patients with MD showed significantly higher levels of 24-h and daytime heart rate, relative to non-clinical controls. With respect to HRV, patients had significantly lower levels of 24-h HRV triangular index, daytime HRV triangular index, daytime HF-HRV, daytime LF-HRV and daytime RMSSD, and a trend of reduced daytime SDNN. For blood pressure, higher levels of 24-SBP and daytime DBP, and a trend of higher daytime SBP and 24-h SBP, were observed in patients with MD compared to controls. There were no significant differences between groups for any other biological measures. Table 2 reports test statistics for differences in changes from baseline to the end of treatment between both treatment groups.

Main treatment effects
Compared to patients in the WL group, patients who received CBT had stronger reductions in self-rated depressive symptoms (d = 0.52) and in clinician-rated depressive symptoms (d = 0.61). Changes in inflammatory markers and measures of heart rate did not significantly differ between patients in the CBT group and patients in the WL group. In terms of HRV measures which were reduced in depressed patients at baseline, the CBT group showed a significant increase in 24-h HRV triangular index (d = 0.48) and in daytime HRV triangular index (d = 0.43) from baseline to the end of treatment, compared to the WL group (shown in Fig. 2). Although 24-h LF-HRV was not significantly reduced in depressed patients at baseline, there was a trend increase in 24-h LF-HRV (d = 0.35) in the CBT group, compared to the WL group (see online Supplementary material Fig. S1). With respect to blood pressure measures which were increased (or marginally increased) in depressed patients at baseline, the CBT group showed a trend decrease in daytime SBP (d = 0.36), compared to patients in the WL group (shown in Fig. 2). There was no evidence for the main treatment effect in other variables of HRV or blood pressure. Compared to intention-to-treat analyses, analyses per protocol (results not shown) did not provide any substantial changes in the estimates and showed the same pattern of significance for treatment effects with one exception: the treatment effect for 24-h LF-HRV became significant in the per-protocol analysis (Intention-to-treat: p = 0.051; per-protocol: p = 0.0499). Table 3 shows results for moderation analyses with baseline selfrated (i.e. BDI-II) and clinician-rated (i.e. MADRS) depressive

3738
Frank Euteneuer et al. symptom severity as a predictor for changes in biological measures from baseline to the end of treatment. In cases where there were any significant or trend moderating effects ( p < 0.10), differences between changes were plotted from lower (25 th percentile) to higher (75 th percentile) levels of the moderator to explore the nature of the moderation. Figure 3 plots moderating effects for biological measures which were altered in patients with MD at baseline (see online Supplementary material Fig. S2 for all other moderating effects). Self-rated (but not clinician-rated) baseline depressive symptom severity significantly moderated differences between changes in 24-h heart rate. In patients with more severe depressive symptoms at study entry, there was a stronger decrease in 24-h heart rate from baseline to the end of treatment in the CBT group, compared to the WL group (see Fig. 3). A significant treatment effect on 24-h heart rate in favor of CBT over WL emerged with a baseline severity score of BDI-II = 30 (d = 0.39). At the 75 th percentile (BDI = 35), the effect on 24-h heart rate was moderate (d = 0.57).
Although not altered at study entry, there was also a trend for this moderating effect in terms of changes in nighttime heart rate (see online Supplementary material Fig. S2).
In terms of HRV measures which were impaired (or marginally impaired) in depressed patients at baseline, higher self-rated depressive symptom severity at baseline was associated with significant improvements in daytime HRV triangular index, daytime HF-HRV, and daytime LF-HRV from baseline to the end of treatment in the CBT compared to the WL group (see Fig. 3 Fig. S2). Depressive symptom severity did not moderate the impact of treatment on inflammatory markers. Compared to intention-to-treat analyses, analyses per protocol (results not shown) did not provide any substantial changes in the estimates and showed the same pattern of significance for moderated treatment effects with one exception: the moderated treatment effect (i.e. self-rated depressive symptom severity) for 24-h HRV triangular index became significant in the per-protocol analysis (Intention-to-treat: p = 0.050; per-protocol: p = 0.049).

Correlations between changes in depressive symptom severity and biological outcomes
In cases where cLDA resulted in any significant or trend effects ( p < 0.10), we performed post-hoc intention-to-treat correlational analyses with maximum likelihood estimation to examine whether a reduction in depressive symptom severity correlated with corresponding improvements in biological outcomes in the CBT group. A decline in self-rated depressive symptom severity from baseline to the end of CBT significantly correlated with a decrease in 24-h heart rate (r = 0.35; p = 0.010) and tended to correlate with a decrease in nighttime heart rate (r = 0.25; p = 0.076). In addition, decrease in self-rated depressive symptom severity significantly correlated with improvement in 24-h HF-HRV (r = −0.43; p < 0.001), daytime HF-HRV (r = −0.35; p = 0.004), 24-h RMSSD (r = −0.49; p < 0.001) and daytime RMSSD (r = −0.41;   p = 0.001) but not with changes in measures of SDNN, LF-HRV and HRV triangular index, as well as daytime SBP (all ps ⩾ 0.1).
Correlations between changes in clinician-rated depressive symptom severity and changes in cardiac measures were non-significant (all ps ⩾ 0.1).

A posterori multiplicity adjustments
As outlined in the methods section, all results are not adjusted for multiple testing and should be considered exploratory (Althouse, 2016;Bender & Lange, 2001;Rothman, 1990;Rubin, 2017). However, because there may not always be consensus on the need for multiplicity adjustment in exploratory studies, we decided to provide separate post-hoc multiplicity adjustments with a false discovery rate (FDR) of 10% in cases where at least one significant treatment effect was observed within a reasonable family of potentially interrelated HR and HRV outcomes. More specifically, we separately adjusted within the domain of 24-h HR and HRV measures, within the domain of daytime HR and HRV measures, as well as within the domain of nighttime HR and HRV measures via the less-conservative two-stage step-up method of Benjamini, Krieger, and Yekutieli (2006) using Prism 9 (GraphPad) for Windows. As shown in Table 2, the main treatments effects on 24-h and daytime HRV triangular index could not be observed after FDR adjustments. As seen in Table 3, all originally significant moderated treatment effects would also be prominent after FDR correction. Importantly, multiplicity adjustments were conducted a posteriori and thus the sample size in this study may not account for these adjustments.

Discussion
This exploratory RCT examined whether CBT is accompanied by reductions in cardiac risk biomarkers in antidepressant-free patients with MD. Biological outcomes included a wide range of measures with potential predictive value for CVD such as heart rate (Kannel et al., 1987), indices of HRV (Hillebrand et al., 2013), blood pressure (Lewington et al., 2002), as well as inflammatory markers (i.e. CRP, IL-6, and TNF-α) (Avan et al., 2018;Pearson, 2003;Subirana et al., 2018). Compared to WL, CBT was associated with a significant increase in 24-h and daytime HRV triangular index and a borderline significant increase in 24-h LF-HRV, as well as a trend decrease in daytime SBP. Self-rated depressive symptom severity at study entry moderated (or tended to moderate) improvements through CBT for 24-h and daytime heart rate, and for several indices of HRV, particularly those daytime measures which were altered in depressed patients at study entry. With increasing symptom severity, CBT increased or tended to increase these indices of HRV and reduced heart rate, compared to the control condition. Inflammatory effects were not observed in this trial. At study entry, the 24-h and daytime HRV triangular index, a robust geometric measure of overall HRV (Stapelberg, Neumann, Shum, McConnell, & Hamilton-Craig, 2018;Vila, Lado, & Cuesta-Morales, 2019;Ziegler et al., 1999), was reduced in patients with MD. Compared to the control condition, CBT increased the 24-h and daytime HRV triangular index. Several studies suggest that overall HRV predict cardiovascular health outcomes in populations without known CVD at study entry, as well as in clinical populations (Fang, Wu, & Tsai, 2020a, 2020bHämmerle et al., 2020;Hillebrand et al., 2013). Although our finding needs confirmatory replication, an improvement of overall HRV during CBT may suggest that this common form of psychotherapy does not only improve depressive symptom severity but also cardiovascular health. Given the increased incidence of CVD in patients with MD, it is of clinical relevance whether common treatments affect cardiovascular health (Gan et al., 2014;Harshfield et al., 2020;Li et al., 2019;Wu & Kling, 2016). Our finding may thus inspire future research to include the HRV triangular index when evaluating biological effects of psychotherapy in MD, and also to examine whether potential improvements in overall HRV during treatment translate into reduced risk for CVD.
Although CBT numerically increased all other indices of HRV, differences in changes between CBT and WL were not significant. One explanation for this finding might be, that the triangular index is more reliable than non-geometric measures of HRV in terms of high intraindividual reproducibility and thus more appropriate for prospective research (Stapelberg et al., 2018;Vila et al., 2019;Ziegler et al., 1999). The major disadvantage of the triangular index is that time-consuming long-term recordings should be used to ensure valid calculation (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, 1996). However, because of its robustness and responsiveness in this trial, the triangular index might be a good measure to evaluate changes in overall HRV through psychological treatment.
The observation that self-rated depressive symptom severity, assessed by the BDI, moderated treatment effects on heart rate and daytime HRV is in line with findings from an earlier pre-post study in depressed patients with coronary heart disease (Carney  Table 2 for test statistics). HRV, heart rate variability. +p < 0.10 *p < 0.05.

3744
Frank Euteneuer et al. et al., 2000). In this previous work, patients were classified according to their BDI score as mildly depressed or moderately to severely depressed. The main result was that heart rate and daytime RMSSD improved significantly in the moderately to severely depressed patients but remained unchanged in the mildly depressed patients (Carney et al., 2000). This earlier study together with our present RCT suggests that potential changes in heart rate and especially in daytime measures of HRV are primarily observable in patients with more severe depressive symptoms. Patients with more severe depressive symptoms thus seem BDI-II, Beck Depression Inventory-II; CBT, Cognitive behavioral therapy; CRP, C-reactive protein; HF-HRV, high-frequency heart rate variability; IL-6, interleukin 6; LF-HRV, low-frequency heart rate variability; MADRS, Montgomery-Asberg Depression Rating Scale; RMSSD, square root of the mean of the sum of the squares of differences between adjacent NN intervals; SDNN, a standard deviation of all NN intervals; SQRT, square root transformed; TNF-α, tumor necrosis factor α; WL, waitlist. a A finding is observed after adjustment with a false discovery rate of 10% (see methods section for details on false discovery adjustment). b Models include corresponding 24-h blood pressure as a time-varying covariate. Fig. 3. Baseline self-rated depressive symptom severity (i.e. BDI-II) as moderator of differences in changes in cardiac measures from baseline to the end of treatment.
Note. Estimated marginal means (standard errors) from constrained linear mixed models are plotted from lower (25 th percentile) to higher (75 th percentile) levels of the moderator (see Table 3 for test statistics). HRV, heart rate variability; HF-HRV, high-frequency heart rate variability; LF-HRV, low-frequency heart rate variability. +p < 0.10 *p < 0.05 **p < 0.01.
to improve their cardiovascular risk profile the most. This finding is of particular interest since also the risk for CVD increases with increasing self-rated depressive symptom severity (Harshfield et al., 2020) and an improvement of HRV during treatment might thus be of stronger prognostic relevance in more severe depression. Higher depression scores have greater potential for change than lower scores. If the degree of change in depressive symptoms relates to changes in heart rate and HRV, one would expect that a reduction in depressive symptom severity relates to improvements in heart rate and HRV during treatment. Indeed, the present study found that a reduction in self-rated depressive symptom severity through CBT significantly correlated with improvements in HF-HRV, RMSSD and heart rate. Importantly, clinician-rated depressive symptoms, as assessed by the MADRS, did neither moderate biological treatment effects nor correlate with changes in a cardiac measure in CBT. One explanation for these diverging findings may be that the BDI, in contrast to the MADRS, has been conceptualized to reflect outcomes of psychotherapy (Demyttenaere & De Fruyt, 2003) and may thus stronger relate to the biological effects of CBT. Moreover, expert ratings are prone to over-estimate clinical improvements (Cuijpers, Li, Hofmann, & Andersson, 2010;Rief et al., 2009). For future psychotherapeutic trials, we suggest considering moderating effects of patients' self-rated symptom severity when studying changes in heart rate and HRV. In the present study, patients exhibited increased levels of inflammatory levels at study entry. There was no evidence for reductions in these markers through CBT. A previous RCT with partially medicated patients with depression and elevated cardiovascular risk did also not observe an effect of CBT on CRP (Taylor et al., 2009). Another RCT with partially medicated patients found that CBT reduced CRP (but not IL-6) only in patients with increased baseline CRP levels, and only in combination with physical exercise (Euteneuer et al., 2017). On the other hand, reductions in IL-6 and TNF-α during CBT (but not during Narrative Cognitive Therapy) have been observed in a study with young patients with MD (18 to 29 years old), but this RCT did only analyze within-group changes and provided no intention-to-treat analyses for differences in changes between both treatment groups (Moreira et al., 2015). Based on all these findings, there may be little or no evidence for inflammatory effects of standard CBT in the context of depression, although some findings in other populations suggest that CBT may improve inflammation, for example, in insomnia (Irwin et al., 2015) or fibromyalgia (Zabihiyeganeh et al., 2019). Even if inflammation may neither be necessary nor sufficient to induce or sustain depression in general (Kiecolt-Glaser et al., 2015), a potential insufficiency of CBT to improve low-grade inflammation strengthens the importance of research efforts to develop more specific anti-inflammatory treatments for subgroups of patients (Jones et al., 2020).
Although findings for associations between depression and blood pressure alterations are mixed (Gould & Beaudreau, 2013;Hildrum et al., 2007;Rogeness, Cepeda, Macedo, Fisher, & Harris, 1990), depression has been shown to predict the development of hypertension (Ginty et al., 2013;Meng et al., 2012). In addition, research using long-term assessment of blood pressure (i.e. 7 days) suggests increased SBP and DBP, as well as reduced nocturnal blood pressure dipping in depressed individuals (Shinagawa et al., 2002). Compared to non-clinical controls, patients in the present study showed significantly increased 24-h SBP and daytime DBP, as well as marginally higher daytime SBP and 24-h DBP. Although these alterations numerically improved in the CBT group compared to the WL group, there was only a trend effect for daytime SBP. To the best of our knowledge, there are no previous RCTs that evaluate CBT effects on ambulatory blood pressure in patients with MD. So, this finding is of potential interest but needs further investigation.
Notwithstanding the strengths of the present study, such as the inclusion of a non-clinical control group, the randomized controlled design, the comprehensive assessment of potential cardiac risk biomarkers, the inclusion of patients without antidepressant medication, and the use of 24-h long-term recordings for cardiac measures instead of short-term recordings, limitations need to be reflected. First, given the exploratory nature of this trial, further research is necessary to confirm the observed findings. In particular, treatment effects for the HRV triangular index may have occurred by chance. Second, our sample consists of outpatients with MD who were eligible for psychological treatment. Thus, findings may not generalize to other samples of depressed patients (e.g. MD patients with psychotic features). Third, we did not assess patients`experience with CBT, a factor which might of relevance for treatment responsiveness. Fourth, although studies on depression and dysregulation of the hypothalamic-pituitary-adrenal (HPA) stress axis are less consistent (Dedovic & Ngiam, 2015;Psarraki et al., 2021;Stetler & Miller, 2011), there is meanwhile evidence for a prospective association between higher levels of morning plasma cortisol and incident CVD (Crawford et al., 2019). It might thus have been useful to include HPA-axis biomarkers in this study. Moreover, this RCT was not designed to examine specific mechanisms (as described in the introduction section) and the relative contribution of specific components of CBT to potential biological effects, which may be of interest for future dismantling studies. Finally, despite longitudinal evidence linking heart rate and HRV to cardiovascular risk, the clinical implications of the present study are not clear. The most important question is whether improving biological cardiac risk factors in patients with MD will translate into reduced risk for CVD. Given the lack of studies which focus on such long-term effects of psychological treatment, it is not possible to answer this question.
To conclude, this study found that CBT for MD is accompanied by an increase in overall HRV, as indexed by the triangular index. If this exploratory result is replicated in confirmatory studies, this could imply that the cardiovascular risk profile is reduced after successful CBT. Moreover, self-rated depressive symptom severity was identified as a potential moderator for improvements in several indices of heart rate and HRV. The findings provide new insights into biological effects of psychological treatment against depression.
Supplementary material. The supplementary material for this article can be found at https://doi.org/10.1017/S0033291722000447.