Abstract
Background Cell deformability of all major blood cell types is increased in depressive disorders (DD). Furthermore, impaired glucocorticoid secretion is causally related to DD. Nevertheless, there are no longitudinal studies examining changes in glucocorticoid output and depressive symptoms regarding cell deformability in DD.
Aim To investigate, whether changes in depressive symptoms or hair glucocorticoids predict cell deformability in DD.
Methods In 136 individuals, depressive symptoms (PHQ-9) and hair glucocorticoids (cortisol and cortisone) were measured at timepoint one (T1), while one year later (T2) depressive symptoms and hair glucocorticoids were remeasured and additionally cell deformability of peripheral blood cells was assessed and DD status was determined by clinical interview.
Results Depression severity at T1 predicted higher cell deformability in monocytes and lymphocytes over the entire sample. Subjects with continuously high depressive symptoms at T1 and T2 showed elevated monocyte deformability as compared to subjects with low depressive symptoms. Depression severity at T1 of subjects with a lifetime persistent depressive disorder (PDD) was associated with elevated monocyte, neutrophil, and granulo-monocyte deformability. Depression severity at T1 of subjects with a 12-month PDD was positively associated with monocyte deformability. Furthermore, increases in glucocorticoid concentrations from T1 to T2 tended to be associated with higher immune cell deformability, while strongest associations emerged for the increase in cortisone with elevated neutrophil and granulo-monocyte deformability in the 12-month PDD group.
Conclusion Continuously elevated depressive symptomatology as well as an increase in glucocorticoid levels over one year are associated with higher immune cell deformability, particularly in PDD. These findings suggest, that persistent depressive symptomatology associated with increased glucocorticoid secretion may lead to increased immune cell deformability thereby compromising immune cell function and likely contributing to the perpetuation of PDD.
Introduction
Of all disorders, depressive disorders (DD) cause the greatest disease burden and disability worldwide (World Health Organisation (WHO), 2017). Major depressive disorder (MDD) is characterized by a period of at least two weeks of depressive mood or anhedonia in combination with other symptoms such as appetite changes, insomnia/hypersomnia, increased fatigue, or feelings of worthlessness (American Psychiatric Association, 2013a). Persistent depressive disorder (PDD; formerly dysthymia) is characterized on the basis of depressed mood over a two-year period in combination with symptoms similar to MDD, although PDD is often perceived as less severe in terms of intensity (American Psychiatric Association, 2013a). Nevertheless, both conditions cause severe functional impairment and there are major difficulties in reducing the high prevalence rates (Hasin et al., 2018; Vos et al., 2017) and improving DD treatment (Munder et al., 2019; Munkholm et al., 2019), with 33% of patients not achieving remission despite multiple treatment attempts (Guidi et al., 2011; Rush et al., 2006). Often, the poor therapeutic efficacy is attributed to the still insufficiently understood pathophysiology of DD (Biselli et al., 2019; Fiacco et al., 2019; Otte et al., 2016; Rothe et al., 2020).
Research increasingly recognizes that key pathophysiological processes occur directly at the cellular level and that DD can impair peripheral blood cell function via disturbed glucocorticoid secretion by the hypothalamus-pituitary-adrenal (HPA) axis or via inflammation (Lopez-Vilchez et al., 2016; Moylan et al., 2013; Pariante, 2017; Rodrigues et al., 2014; Walther et al., 2018; Wolkowitz et al., 2010). Blood cells fulfill multiple functions from primary immune response, over metabolite transport to overall blood flow. Real-time deformability cytometry (RT-DC) representing a new method to investigate the physical properties of blood cells such as cell mechanical features or cell size provides unprecedented insights into overall blood cell function (Otto et al., 2015a; Rosenbluth et al., 2008a). It has been highlighted that the assessment of the mechanical state of blood cells, measured by cell deformability, is appropriate to detect and classify human diseases conditions such as spherocytosis, malaria, or COVID-19 infections (Kubánková et al., 2021; Toepfner et al., 2018a) and is predictive for immune cell activation (Bashant et al., 2020, 2019).
Only recently, using RT-DC, we were able to show that in DD and especially in PDD immune cell deformability is increased as compared to healthy controls (Walther et al., 2022b). Although all major blood cells tended to be more deformable, lymphocytes, monocytes and neutrophils were particularly affected indicating immune cell mechanical changes to occur in DD. This is potentially related to disturbed HPA axis function or dysfunctional immune response. In addition, optical traps (Ekpenyong et al., 2012), atomic force microscopy (Lam et al., 2008; Rosenbluth et al., 2008b), and micropipette aspiration (Ravetto et al., 2014) have been used in proof-of-concept studies to reveal mechanical changes in immune cells under pathological conditions. Nevertheless, due to the general dominance of erythrocytes in blood and the consequent much increased statistical power to detect associations, previous studies investigated only the relationship between erythrocyte deformability and mental disorders (Jasenovec et al., 2019; Saha et al., 2019). The study by Saha et al. (2019) reported patients with chronic fatigue syndrome / myalgic encephalomyelitis to exhibit lower erythrocyte deformability than healthy controls. Further, Jasenovec et al. (2019) examined children within the autism spectrum identifying those children with more severe symptoms to show impaired erythrocyte deformability. It should be noted here, however, that the current gold standard for measuring cell deformability, namely RT-DC (Kräter et al., 2021; Nawaz et al., 2020; Otto et al., 2015a; Urbanska et al., 2020; Wu et al., 2018), was only used in our landmark study comparing individuals with DD and healthy controls (Walther et al., 2022b).
In contrast, there is a large body of literature examining the relationship between impaired glucocorticoid secretion by the HPA axis and DD (Rothe et al., 2020). Indeed, it seemed to be widely accepted that MDD is associated with increased HPA activity and increased basal cortisol concentrations in the circulation (Otte et al., 2016; Stetler and Miller, 2011). However, due to the inconclusive findings on cortisol stress reactivity in MDD, this is no longer definitive (Miller and Kirschbaum, 2019; Rothe et al., 2020; Zorn et al., 2017). Furthermore, investigating new tissue sources such as hair cortisol showing mixed associations with regard to DD highlight on the one hand the relevance of the investigation of chronic glucocorticoid output over longer time periods (Gerritsen et al., 2019; Rothe et al., 2021, 2020; Stalder et al., 2017; Steudte-schmiedgen et al., 2017; Walther et al., 2022a), but on the other hand, it emerges as increasingly important to investigate the effects of impaired chronic glucocorticoid secretion on cell function. Studies increasingly show that the extracellular microenvironment is critical to modify cellular physiology, which ultimately determines cell functionality (Kim et al., 2019). So, understanding the interactions between immune cells and DD-related alterations in glucocorticoid secretion appears to be increasingly important for understanding differences in immune cell deformability in health and disease.
In addition to increased cortisol levels and the often-identified chronic low-grade inflammation in individuals suffering from DD (Pariante, 2017), Lynall et al. (2019) further reported elevated numbers of immune cells in depressed individuals compared to controls. It is known that increased levels of glucocorticoids and catecholamines lead to increased immune cell count, as cells demarginate from the vessel walls. Interestingly, these observations were associated with cellular softening (Fay et al., 2016). It is suggested that the underlying mechanism leading to increased immune cell counts in depressed individuals is rooted in the effect of elevated glucocorticoid levels remodeling the actin cytoskeleton of blood cells and thereby softening leukocytes and enabling them to demarginate from the vessel wall (Fay et al., 2016). Not only starts the actin cytoskeleton to be remodeled with continuously elevated levels of glucocorticoids, but also are lipid metabolism and composition crucially affected resulting in increased softening and bending of blood cells (Demirkan et al., 2013; Knowles et al., 2017; Liu et al., 2016; Ravetto et al., 2014; Walther et al., 2018). This cascade of processes ultimately impairs blood cell function and may be an underlying cause of symptoms of fatigue and exhaustion in DD. The disturbed immune cell function in connection with subsequent dysregulated cytokine release could trigger sickness behavior with anhedonia and fatigue (Dantzer, 2009).
Therefore, we hypothesize that continuously high levels of depressive symptoms are associated with increased immune cell deformability. We further hypothesize that high and especially continuously high glucocorticoid levels as well as an increase in glucocorticoid secretion over one year is associated with elevated immune cell deformability in individuals with DD.
Method
Study design
This observational study used a longitudinal design with a baseline measurement and a one-year follow-up measurement. Cross-sectional results of the study were reported elsewhere (Walther et al., 2022b). Baseline examination of participants took place between October and December 2018, while follow-up examinations were conducted one year later. The study, entitled Mood-related Morpho-rheological Changes in Peripheral Blood Cells (Mood-Morph), used a participant pool of a large-scale prospective cohort study to identify eligible participants and obtain pre-assessed baseline data (Penz et al., 2018). Eligible participants were asked via e-mail whether they were willing to participate in an add-on study to the prospective longitudinal study.
For the study, participants with elevated depressive symptoms, as measured by the Patient Health Questionnaire (PHQ-9), and gender- and age-matched healthy controls were identified from the participant pool. These participants provided for the baseline measurement (T1) sociodemographic information, depressive symptomatology and further health-related data. Also, collected hair samples for glucocorticoid quantification were obtained for the study at T1. At the follow-up time point (T2), the psychometric assessment and hair sampling was repeated complemented with a clinical diagnostic interview to diagnose depression status and blood samples were obtained to measure cell deformability. The study was approved by the local ethics committee of the Dresden University of Technology (EK182042019) and all participants gave written informed consent to participate in the study.
Participants
For study inclusion, participants’ age had to fall between 18 and 68 years, and interested study participants were excluded if they had any kind of blood disorder that could affect blood cell deformability. Potential participants who were contacted to participate in the Mood-Morph study, were informed in the study invitation and informational email that certain diseases, such as spherocytosis and blood-related diseases in general, precluded participation (Toepfner et al., 2018a). Potential participants with a cold or other acute infections no longer ago than two weeks, were not included to exclude acutely infectious participants. Acute medical conditions such as blood diseases, cancer, severe cardiovascular diseases, or pregnancy led to exclusion, while in the population very prevalent conditions such as hypertension or hypothyroidism did not lead to exclusion. Differing conditions were recorded and used as covariate in the analyses as well as the specific drug categories (e.g., for hypertension—antihypertensive drugs).
In terms of depression levels measured with the PHQ-9, only subjects with a T1 PHQ-9 score >10 (high risk of depression) as well as age- and sex-matched subjects with a PHQ-9 score <5 (low risk of depression) were included. The aim of this procedure was to obtain a sample of participants with high and low depressive symptoms at baseline and to be able to identify change in depressive symptoms over one year. At T2, participants underwent a clinical diagnostic interview leading to exclusion of participants from the DD groups in case participants did not fulfill diagnostic criteria, had a positive mania screening, or suffered from another major psychological disorder. Exclusion from the healthy control group occurred if participants exhibited a major psychological disorder or psychopharmacological medication intake.
Procedures
Participants provided sociodemographic and health-related information (e.g., gender, age, medical conditions, medication) as well as depressive symptom severity and information on their general health online via the study-homepage platform. By using a personalized code, participants accessed a private study space, signed first consent forms regarding data privacy, clinical data collection and saving, and hair/blood sampling and completed subsequently their study questionnaires. For depression severity the PHQ-9 and for general health the Short Form Health Questionnaire (SF-12) was completed (Kroenke et al., 2001; Ware Jr et al., 1996). While the SF-12 measures general physical and psychological health, the PHQ-9 is a nine-item self-report questionnaire for the measurement of depression severity rating the frequency of DSM-5 diagnostic criteria for MDD (e.g., feeling down, depressed, or hopeless) during the past two weeks on a four-point Likert scale (0 = not at all, 3 = nearly every day). Scores of 10 or greater were consistently shown to identify with high sensitivity and specificity cases of MDD (He et al., 2020; Levis et al., 2020). Since the scale developers suggested that PHQ-9 scores of 5, 10, 15, and 20 represent mild, moderate, moderately severe, and severe depression (Kroenke et al., 2001), individuals with scores below 5 were considered to have minimal depressive disorder risk, whereas individuals with scores of 15 or higher were considered to have high depressive disorder risk.
To obtain capillary blood and scalp hair samples, participants arrived at the Dresden University of Technology at the Department of Biopsychology, where trained personnel performed blood and hair sampling. At T2, after the bio-sampling took place, the Composite International Diagnostic Interview (DIA-X-5/CIDI) was conducted by trained interviewers (Hoyer et al., 2020). The DIA-X-5 is a standardized clinical interview for the assessment of mental disorders. The DIA-X-5 assesses lifetime symptoms and past 12 months symptoms. To rule out bipolar disorders, a mania-screening questionnaire, consisting of the initial questions of the DIA-X-5 mania section, was conducted to detect any lifetime (hypo-) mania-symptoms. Lifetime and 12-month diagnoses were subsequently generated for PDD and MDD according to DSM-5 criteria (American Psychiatric Association, 2013b). After finishing the clinical interview, subjects received 15 € compensation for expenses. Duration of the examination at T2 was approximately one hour.
Real-time deformability cytometry (RT-DC)
A 20 µl capillary blood sample was extracted from the fingertip of participants using a safety lancet (Safety-Lancet Normal 21, Sarstedt, Nümbrecht, Germany) and harvested in a capillary (Minivette POCT, 20 µl, Sarstedt, Nümbrecht, Germany). The blood was immediately resuspended in 380 µl RT-DC measurement buffer containing 0.6% methylcellulose in phosphate buffered saline (CellCarrierB, Zellmechanik Dresden, Germany) in a microcentrifuge tube. Blood-samples were then transferred to the Department of Cellular Machines at the Biotechnology Center of the TU Dresden, where the samples were measured using an RT-DC device. Maintained at room temperature samples were measured within three hours since sampling according to a protocol published elsewhere (Toepfner et al., 2018b).
In brief, blood was flushed through a microfluidic channel constriction 20 µm x 20 µm in cross section (Flic20, Zellmechanik Dresden, Germany) by applying a constant flow rate. An image of every measured blood cell was taken by a high-speed camera and cell deformability and cell size were calculated (Otto et al., 2015b). RT-DC measurements were controlled by the acquisition software Shape-In2 (Zellmechanik Dresden, Germany). The different blood cell types were classified by utilizing artificial intelligence-based image classification as published elsewhere (Kräter et al., 2021) and mean values for cell deformability and cell size of every donor and blood cell type were extracted.
Hair sampling and liquid chromatography and mass-spectrometry (LC-MS)
Three hair strands (each containing at least 20 mg of hair) were cut as close as possible to the scalp from the posterior vertex for quantification of hair cortisol and cortisone. Hair sample processing and analysis were performed as described in our published protocols for glucocorticoids (Gao et al., 2016, 2013). Preprocessing of the samples and subsequent biochemical analysis were performed by Dresden LabService GmbH (Tatzberg 47, 01307, Dresden, Germany). Samples were cut into 3 cm segments from the scalp site representing the integrated hormone concentration over the last three months (e.g., average hair growth rate of 1 cm per month; (Wennig, 2000)) and weighed into 7.5 mg of whole, non-pulverized hair and then washed with 2.5 ml of isopropanol according to the protocol of Gao et al. (Gao et al., 2013). Biochemical analysis was performed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) as described in detail elsewhere (Gao et al., 2016). The intra- and inter-assay coefficients of variation for cortisol and cortisone are less than 8.8% and 8.5%, respectively. The lower limit of detection for cortisol and cortisone is 0.09 pg/mg.
Statistical analysis
Statistical analyses were run in R v4.2.1 (R Core Team, 2017). Data was checked for normality using Shapiro-Wilk-tests and for variance homogeneity using F tests. In case of violation of normal distribution, non-parametric tests were conducted. Differences in confounding variables (sex, age, BMI, psychopharmaceutical intake, medication category) and in PHQ-9, cortisol, cortisone and WBC count were tested between subjects suffering from depression and healthy controls using two-tailed independent t-tests (Student’s or Welch two sample t-test), Mann-Whitney U tests, univariate analysis of variance (ANOVA) and χ2-tests in dependence of normal distribution and scale of measurement. Pearson and Kendall correlation coefficients (Croux and Dehon, 2010) were tested to check for possible (non-)linearity relationships. Partial correlations for T1 and T2 data were conducted to control for gender, age, BMI and psychopharmaceutical treatment, while for change scores (T2-T1) a baseline T1 correction (e.g., PHQ-9 T1 or cortisol (or cortisone) T1) was additionally applied. All p-values were adjusted using Holm-Bonferroni correction for multiple comparisons (Holm, 1979). Cortisol and cortisone data from both measurements were log-transformed due to their non-normal distribution and tested two-tailed. Correlation analyses were performed over the entire sample size and in each subgroup with differing depressive disorders.
Results
Sample description
In total, 136 (nfemale = 100; nmale = 36) individuals, who completed the PHQ-9 and provided hair samples for glucocorticoid quantification at T1 and T2 were used for subsequent analyses investigating the longitudinal association between depressive symptoms, glucocorticoid secretion and immune cell deformability. Moreover, the selected participants (N = 136) were assigned to different diagnostic groups in dependence of their DIA-X diagnoses. Of initially 139 participants, due to incomplete PHQ-9 data of three individuals, following group size emerged from the diagnostic classification for the 136 participants: Lifetime PDD (n = 30), 12-month PDD (n = 15), lifetime MDD (n = 40), 12-month MDD (12-month MDD; n = 11) and healthy controls (HC; n = 61). Participants suffering from double depression were not treated as a separate group due to small sample size and lack of statistical power. A detailed representation of the sample flow is provided in Figure 1. Sample characteristics of the entire sample and the different diagnostic groups is presented in Table 1 and Supplementary Table 1. With respect to hair cortisol and cortisone analyses, sample size of the initially entire sample was used (N = 139). However, as hair segments could not be obtained from all subjects due to lack of sufficient hair of some participants, sample size varies within each group. Regarding cortisol analyses, one outlier was additionally excluded due to a hair cortisol concentration above 100 pg/mg. For cortisol at T1 and T2 n = 116 observations were available at each time point for analysis, while for cortisone at T1 and T2 n = 118 observations were available at each time point.
Sample Flow Leading to Relevant Case-Control Groups
Sample characteristics according to the entire sample
Control of Confounders
Participants in subgroups with differing depressive disorders showed no significant differences regarding age, sex and BMI compared to healthy controls. Expected significant group differences were observed in PHQ-9 (T1 and T2, T2-T1) scores and psychopharmaceutical treatment within all groups except PHQ-9 symptom change score (T2-T1) for the 12-month PDD group and controls. Exact group differences of the diagnostic groups are shown in Supplementary Table 1. Subgroups with differing depressive disorders and healthy controls did not show a significant difference in glucocorticoid concentrations (T1, T2, T2-T1; see Supplementary Tables 2 and 3 and Supplementary Figures 1–6). No significant correlations between depressive symptoms (PHQ-9 T1 and PHQ-9 T2) and mean cell area size over the entire sample were found in any cell type (see Supplementary Table 4 and 5). No significant differences in immune cell count were found between participants with a DD diagnosis and healthy controls (see Supplementary Table 6). Furthermore, cell deformability did not differ with respect to the six different medication groups (see Supplementary Text 1).
Association Between Depression Severity and Cell Deformability
Pearson, Kendall and Partial correlation coefficients between PHQ-9 (T1) scores and cell deformability for each cell type are depicted in Table 2. Below, only significant Pearson correlations are described. Significant correlations between PHQ-T1 and monocytes were found for the entire sample (r = .22, p = .0051) and the lifetime PDD group (r = .33, p = .0049) both corrected for confounder variables (entire sample: rp = .20, p = .01; lifetime PDD: rp = .32, p = .0085) and applied to Holm-Bonferroni correction (entire sample: p = .0245, pp = .0425; lifetime PDD: p = .0245, pp = .0425). Furthermore, within the lifetime PDD group, significant correlations between PHQ-9 scores (T1) and neutrophil deformability (r = .30, p = .0097, rp = .33, p = .0072) as well as granulo-monocytes (r = .32, p = .0057; rp = .34, p = .0057) emerged, indicating higher immune cell deformability in subjects, who suffered from PDD (see Figure 2). Pearson, Kendall and partial correlations remained significant after Holm-Bonferroni correction (p = .0485, pp = .036). No further correlations did survive correction for multiple testing. Nevertheless, with respect to the entire sample a significant correlation between PHQ-9 (T1) and lymphocyte deformability was observed (r = .16, p = .0326; rp = .17, p = .0270). Within the 12-month PDD group, significant Pearson and Kendall correlations between PHQ-9 scores (T1) and monocyte deformability were found (r = .35, p = .0306), but not when confounders were taken into account. Within the 12-month PDD group, a significant positive partial correlation with respect to erythrocyte deformability was found (r = .36, p = .0353). No stable significant correlations emerged for the lifetime MDD and 12-month MDD group. Significant results indicate a positive correlation of PHQ-9 scores at T1 with immune cell deformability in the lifetime PDD group and partially over the entire sample. Correlations between PHQ-9 scores (T2) and cell deformability for each cell type was performed over the entire sample revealing significant correlations for monocyte cell deformability (r = .17, p = .0254; rp = 16, p = .036; see Supplementary Table 7). Because of the limited variance between T1 and T2 of PHQ-9 scores, we did not focus on the association of depressive symptoms difference scores and cell deformability. Nevertheless, an overview of the correlations can be seen in the Supplementary Table 8.
Associations of Depressive Symptoms Measured by PHQ-9 (T1) and Cell Deformability According to Lifetime PDD group
Associations Between Depressive Symptoms Measured by PHQ-9 (T1) and Cell Deformability
Association Between Cortisol and Cortisone (T1 and T2) and Cell Deformability
No significant association between cortisol at T1 and immune cell deformability of any cell type was revealed over the entire sample or in any subgroup (see Supplementary Table 9). Similarly, between cortisol level T2 and immune cell deformability no significant correlation emerged (see Supplementary Table 10). Also, no significant association between cortisone level T1 and immune cell deformability emerged over the entire sample or in any subgroup (see Supplementary Table 11). A single significant Pearson correlation between cortisone T2 and lymphocyte deformability was observed in the 12-month PDD group (r = .38, p = .0419) (n = 29), which did not survive correction for multiple testing. No further significant correlations were observed between cortisone T2 and immune cell deformability for any cell type in any subgroup (see Supplementary Table 12).
Overall, no stable significant correlations between cortisol, cortisone (T1 and T2) and immune cell deformability were observed for the entire sample nor in any subgroup.
Association Between Cortisol and Cortisone Change (T2-T1) and Cell Deformability
To examine the association between cortisol or cortisone change from T1 to T2 and cell deformability independent of group affiliation, Pearson and Kendall correlations were computed. Significant correlations between cortisol and cortisone difference (T2-T1), respectively, and neutrophil deformability as well as granulo-monocyte deformability emerged for the lifetime PDD group (see Figure 3) and the 12-month PDD group. The results indicate that in those groups, subjects with an increase in cortisol or cortisone over the two assessments exhibit higher neutrophil or granulo-monocyte deformability at T2. Significant partial correlations remained when controlling for gender, age, BMI, psychopharmaceutical treatment and baseline cortisol or cortisone levels (T1) in the 12-month PDD group. However, only the association between cortisone difference (T2-T1) and neutrophil and granulo-monocyte deformability survived Holm-Bonferroni correction in the 12-month PDD group (neutrophil deformability: p = .0455, pp = .0062; granulo-monocyte deformability: p = .019, pp = .0034). No stable significant correlations for lifetime MDD and 12-month MDD were observed. An overview of correlations is presented in Table 3 and Table 4.
Associations of Cortisone Difference (T2-T1) and Cell Deformability According to 12-month PDD group
Associations Between Cortisol Difference (T2-T1) and Cell Deformability
Associations Between Cortisone Difference (T2-T1) and Cell Deformability
Cell Deformability Comparisons According to high and minimal DD risk Groups
Comparing the groups with high vs. minimal DD risk at T1 (high: PHQ-9 ≥ 15; minimal: PHQ-9 < 5), Welch two sample t-tests revealed a significant difference in monocyte deformability at T2 between the high risk group (M = .0585, SD = .0041) and the minimal risk group (M = .0549, SD = .0069; t(75.238) = −3.0642, p = .0015) indicating a moderate effect size (d = .62; see Figure 4). A significant difference in lymphocyte deformability at T2 between the high-risk group (M = .0294, SD = .0017) and minimal risk group (M = .0287, SD = .0018) was revealed by a Mann-Whitney U test (U = 677.5, p = .0337), indicating a small effect size (r = −.19). Furthermore, subjects in the consistent high risk group exhibiting at both time points PHQ-9 scores equal or above 15 (M = .0589, SD = .0035) showed higher monocyte deformability compared to the consistent minimal risk group exhibiting at both time points PHQ-9 scores lower than 5 (M = .0550, SD = .0067; t(37.138) = −2.9601, p = .0027) revealing a moderate effect size (d = .74; Figure 5). No further significant differences were found for any other cell type (see Supplementary Table 13).
Combined Boxplots for the High DD Risk Group and the Minimal DD Risk Group (T1) Regarding Blood Cell Deformability
Combined Boxplots for the Consistent High DD Risk Group and the Consistent Minimal DD Risk Group Regarding Blood Cell Deformability
Association of Cortisol and Cortisone (T1 and T2) with Immune Cell Count
Significant Pearson correlations emerged between cortisol and immune cell count in both PDD groups (see Supplementary Tables 14 and 15). Specifically, cortisol T1 was associated with the number of granulo-monocytes in the lifetime PDD (r = .29, p = .0377) and the 12-month PDD group (r = .52, p = .0035). Furthermore, in the 12-month PDD group, cortisol and cortisone at T1 were associated with neutrophil cell count (cortisol T1: r = .46, p = .0123, cortisone T1: r = .47, p = .0102), while cortisol T2 was associated with granulo-monocyte cell count (r = .3759, p = .04447) and cortisone T2 with neutrophil cell count (r = .38, .0416). Moreover, cortisone T1 and T2 were associated with the number of lymphocytes in the lifetime PDD group (T1: r = .3036, p = .0303; T2: r = .2758, p = .04996). However, only the association of cortisol and cortisone T1 with granulo-monocyte cell count remained significant in the 12-month PDD group after applying Holm-Bonferroni correction (cortisol T1: p = .0176; cortisone T1: p = .0153).
Discussion
Summary of Results
Depressive symptoms at T1 were positively associated with monocyte and lymphocyte deformability at T2 over the entire sample, while individuals with high depressive symptoms at T1 and consistently high depressive symptoms at both time points (T1 and T2) exhibited increased monocyte and lymphocyte deformability at T2 as compared to individuals with no or minimal depressive symptomatology. This pattern was most pronounced in the 12-month PDD and the lifetime PDD groups. An increase in depressive symptoms from T1 to T2 was not associated with cell deformability at T2. Basal glucocorticoid levels at T1 or T2 were further not associated with cell deformability at T2. Albeit not in the entire sample, in both the 12-month PDD and the lifetime PDD group, an increase in glucocorticoid secretion from T1 to T2 was associated with higher granulo-monocyte and neutrophil deformability. Similarly, higher glucocorticoid levels at T1 or T2 were associated with increased granulo-monocyte, neutrophil, and lymphocyte count in 12-month and lifetime PDD groups. Notably, there was no difference in cell size nor in immune cell counts between DD groups and healthy controls, nor was there an association between depressive symptoms and cell size. In addition, no difference in glucocorticoid levels was observed between DD groups and healthy controls, neither was medication intake associated with cell deformability.
Integration of Findings
The present study is the first to examine the longitudinal relationship between depressive symptoms, glucocorticoid secretion and immune cell deformability in depressed individuals and healthy controls. The finding that elevated depressive symptomatology is associated with higher immune cell deformability one year later is supported by previously reported cross-sectional findings on the positive association between depressive disorder status, depressive symptoms and immune cell deformability (Walther et al., 2022). Because these supporting findings originate from the same study, it is important that further longitudinal as well as experimental studies replicate and extend these associations.
The observed significant correlations between depressive symptoms at T1 and immune cell deformability at T2 for monocytes and lymphocytes in the entire sample and for monocytes, neutrophils, and granulo-monocytes in the lifetime PDD and healthy control subsample suggests these immune cells to be the most sensitive cells to react to depressive symptomatology with physical changes. Since these cells are rich in glucocorticoid receptors (Goecke et al., 2007; Heiske et al., 2003; Li et al., 2006; Lu et al., 2017), and HPA-axis alterations are considered a pathophysiological landmark of DD (Rothe et al., 2020; Stetler and Miller, 2011), it is reasonable to assume that changes in glucocorticoid secretion might also affect immune cell status and deformability. Indeed, it was previously shown that increased glucocorticoid levels led to increased immune cell counts due to cell demargination from the vessel walls, which was mediated by cellular softening (Fay et al., 2016). However, in the present study, no direct association was identified between glucocorticoid levels at T1 or T2 with cell deformability at T2. Nonetheless, when examining the change in glucocorticoid secretion between T1 and T2, an increase in glucocorticoid secretion was associated with increased neutrophil and granulo-monocyte deformability at T2 in the lifetime and 12-month PDD and healthy control samples. Notably, we did not observe elevated levels of circulating immune cells in individuals suffering from DD in the present study (see Supplementary Table 6), although this was previously reported by another study (Lynall et al., 2019). In line with previous reports (Gerritsen et al., 2019; Rothe et al., 2020), no direct association between hair glucocorticoid levels and depressive disorder status emerged (see Supplementary Table 2 and 3). However, higher glucocorticoid levels were positively associated with immune cell count in the 12-month and lifetime PDD groups (see Supplementary Table 14 and 15), supporting the perspective that increased glucocorticoid levels support the demargination of immune cells from vessel walls mediated by cellular softening (Fay et al., 2016).
The results suggest that changes in glucocorticoid concentrations over a one-year period may affect immune cell structure and stability. Immune cells might respond to increased glucocorticoid concentrations by reorganizing the cytoskeleton, which affects cell deformability directly as previous studies showed using dexamethasone administration (Fay et al., 2016; Ronchetti et al., 2018). Moreover, elevated glucocorticoid concentrations can cause increased permeation through the cell membrane affecting the lipidome and the cell structure leading to membrane destabilization and bending (Walther et al., 2018), and hence to increased cell deformability. Alterations in membrane-forming lipids in the central nervous system have already been linked to DD (Müller et al., 2015).
The finding that 12-month and lifetime PDD show the strongest associations with cell deformability could be explained by the persistence of disease burden. Indeed, in analyses comparing those individuals with consistently high depressive symptoms (consistent high risk) with individuals with consistently low depressive symptoms (consistent minimal risk) regardless of the DD diagnosis at T2, we found further evidence for this perspective. Individuals exhibiting severely elevated depressive symptomatology at T1 as well as consistently severely elevated depressive symptomatology at both timepoints (T1 and T2) compared to individuals with no or minimal depressive symptoms, showed increased monocyte and lymphocyte deformability at T2. The results suggest that the time factor, namely suffering from severe depressive symptoms at least over a one-year period, might be central with regard to immune cell deformability in DD. Even though significant results can only be observed for monocytes and lymphocytes in the present analysis due to the rather small group sizes, it can be seen from Figure 2 and 3 that this tends to be the case for all immune cells. Given the persistence of depressive symptoms in PDD over two years, associations with immune cell deformability are more likely to be detected than in MDD, which only requires as persistence of symptomatology for over two weeks (APA, 2013). It is also likely that in PDD, the HPA axis shows hyperactivity with increased glucocorticoid output, as there is less depressive symptomatology but over a longer period of time than in MDD. In MDD, the HPA axis appears to enter exhaustion at some point, at which point a shift from increased to decreased levels of glucocorticoids are more likely (Rothe et al., 2020; Steudte-Schmiedgen et al., 2017, 2016; Walther et al., 2022a). Thus, many reports also identify reduced hair cortisol levels or more complex patterns in MDD and especially in conjunction with experienced trauma (Cantave et al., 2022; Hinkelmann et al., 2013; Psarraki et al., 2021; Stalder et al., 2017). Therefore, future studies examining the relationship between DD, glucocorticoid output and immune cell deformability should aim to assess the time period since when individuals are suffering from depressive symptomatology as well as childhood maltreatment or trauma experience to better disentangle these seemingly complex associations.
Overall, the present study provides new insights into the pathophysiology of DD, assuming more deformed immune cells, especially monocytes, with persistent depressive symptomatology. In addition to increased glucocorticoid output, research needs to investigate whether further processes consistently linked to DD such as low-grade inflammation (Goldsmith et al., 2016; Iob et al., 2020; Pariante, 2017) or oxidative stress (Moylan et al., 2014) might be additional factors contributing to cytoskeletal alterations and membrane destabilization (Wilson & González-Billault, 2015; Wong et al., 2013) resulting in higher immune cell deformability in DD. This might be because, similar to increased glucocorticoid secretion, an inflammatory state leads to demargination of blood cells from the vessel walls into the circulation, with cell deformability being increased by the migration process (Fay et al., 2016). Additionally, it was shown that an experimentally induced immune activation with lipopolysaccharide increased monocyte deformability (Ravetto et al., 2014), highlighting chronic low-grade inflammation in depressed individuals as potential agent for increased immune cell deformability. These assumptions are supported by findings showing that individuals with increased inflammatory signaling, such as healthy individuals after inhalation of lipopolysaccharide from E. coli, or individuals suffering from acute lung injury, viral respiratory infections, or Epstein-Barr virus infections, show increased deformability of neutrophils, monocytes, or lymphocytes (Toepfner et al., 2018a). Moreover, neutrophils of individuals in the acute phase of COVID-19 infection showed higher deformability, which was still observable also seven months after the acute infection symptoms (Kubánková et al., 2021), suggesting that long-term HPA-axis or inflammatory alterations may underlie immune cell deformability in DD.
It is important to mention that medication intake was not associated with immune cell deformability. As subjects were not drug naïve, we tested the association of different medication groups with cell deformability, such as psychopharmaceutical medication, antihypertensive drugs or thyroid dysfunction medication. No association between any medication group and cell deformability was identified. Also, associations between depressive symptomatology and cell deformability cannot be traced back to cell size, as higher depression severity and cell size of immune cells were not associated with each other, nor did cell size differ between groups with DD and matched healthy controls.
Strengths and Limitations
Some limitations must be taken into account when interpreting the results. The most important limitation is the single measurement of cell deformability at T2. To investigate the predictive value of depressive symptoms, glucocorticoid concentrations as well as their changes over a one-year period on cell deformability, we measured cell deformability as an outcome measure at T2. However, analyses of changes in cell deformability in relation to changes in depressive symptoms or glucocorticoid concentrations would also be important to investigate. Furthermore, with such a study design, the overall stability of cell deformability markers could be recorded. A related point is the single assessment of the diagnostic groups by clinical interview at T2. Since we intended to capture group differences in cell deformability (Walther et al., 2022b), the clinical diagnostics as well as the cell deformability measurement had to be performed at the same time point in order to avoid possible time effects. However, in a follow-up project, clinical diagnostics will also be performed at two time points together with cell deformability measurement. Nevertheless, the present findings show for the first time that depressive symptoms are related to immune cell deformability one year later, whereas the increase in glucocorticoid secretion over one year is also associated with increased immune cell deformability. Following on from this, however, it should be noted that other important systemic markers have not been investigated in the present study. Thus, it will be important in future studies to explore potential associations between inflammatory markers such as Interleukin-6 or C-reactive protein and immune cell deformability. In addition, the present sample of 136 participants might be underpowered to detect small to moderate effects, so only relatively large effects were probably identified in the present study further supporting the strong relationship between depressive symptoms, glucocorticoid secretion and immune cell deformability. Finally, gender differences could also play a role not revealed in this study, so that larger mixed samples or samples that focus only on women, men, or gender-diverse individuals need to be further investigated.
Conclusions
Taken together, this study represents the first longitudinal investigation of the association between depressive symptoms, glucocorticoid concentrations, and peripheral immune cell deformability in individuals with DD and healthy controls. Our results suggest that higher levels of depressive symptoms are longitudinally associated with higher immune cell deformability, particularly in monocytes and lymphocytes. Moreover, we observed partial confirmation for the hypothesis that DD-related increases in glucocorticoid concentrations lead to cytoskeletal changes and cell membrane destabilization causing overall increased immune cell deformability. Increased immune cell deformability may further lead to a gradual loss of immune cell functionality, so that, for example, immune response is impaired or membrane function is compromised presenting a potentially underlying mechanism causing or maintaining DD.
Contributors
AW, JE and MK designed the study, performed experiments and data analysis. AW wrote the first draft of the manuscript and and JE wrote the statistical analysis and results section of the manuscript. MK, WG, MW, MP, NR performed experiments. CK, JG, MK, LW provided funding, infrastructure/equipment, and edited subsequent versions of the manuscript.
Declaration of interest
All authors state that they have no actual or potential conflict of interest to declare, including any financial, personal, or other relationships with other people or organisations within 3 years of beginning the submitted work that could influence or bias their work.
Data and code availability
The anonymized data and code will be made available to all interested parties upon request.
Role of the funding source
Research pool TU Dresden (F-004242-552-848-1040103) awarded to AW
Faculty of Psychology of the TU Dresden (MK201911) awarded to AW
Swiss National Science Foundation (PZPGP1_201757) awarded to AW
The funding sources are national and university funding sources and had no influence on the writing of the manuscript or the decision to submit it for publication. None of the authors received financial incentives from industrial companies to write this publication. The corresponding author had full access to all study data and had the final responsibility for the decision to submit the paper for publication.
Supplementary Material
Sample Characteristics According to Different Diagnostic Groups
Mean Group Comparisons of Cortisol Between Different Diagnostic Groups and Healthy Controls
Mean Group Comparisons of Cortisone Between Different Diagnostic Groups and Healthy Controls
Associations Between PHQ-9 (T1) and Cell Area According to the Entire Sample
Associations Between PHQ-9 (T2) and Cell Area According to the Entire Sample
Mean Difference of Immune Cell Count Between Different Diagnostic Groups and Matched Healthy Controls
Associations Between Depressive Symptoms Measured by PHQ-9 (T2) and Cell Deformation Over the Entire Sample
Associations Between Depressive Symptoms Difference Measured by PHQ-9 (T2-T1) and Cell Deformability Over the Entire Sample and According to Different Diagnostic Groups
Associations Between Cortisol T1 and Cell Deformability Over the Entire Sample and According to Different Diagnostic Groups
Associations Between Cortisol T2 and Cell Deformability Over the Entire Sample and According to Different Diagnostic Groups
Associations Between Cortisone T1 and Cell Deformability Over the Entire Sample and According to Different Diagnostic Groups
Associations Between Cortisone T2 and Cell Deformability Over the Entire Sample and According to Different Diagnostic Groups
Mean Cell Deformation Comparisons Between High and Minimal Depressive Disorder Risk Group Assessed by PHQ-9 Scores (T1)
Associations Between Cortisol T1, T2 and T2-T1 and Cell Count Over the Entire Sample and According to Different Diagnostic Groups
Associations Between Cortisone T1, T2 and T2-T1 and Cell Count Over the Entire Sample and According to Different Diagnostic Groups
Supplementary Text 1
Drug intake, measured by the subject questionnaire, was analysed by forming six groups, which are: no medication (n = 76), psychopharmaceutical medication (n = 14), antihypertensive drugs (n = 16), thyroid dysfunction medication (n = 16), others (n = 4) and an intake combination group of the preceding categories (n = 10). ANOVA revealed no significant difference between the medication intake groups and cell deformation for each cell type [red blood cells (F(5, 130) = 1.163, p = .331), monocytes (F(5, 130) = 1.112, p = .357), neutrophils (F(5, 130) = .722, p = .608), lymphocytes (F(5, 130) = .935, p = .461), granulo-monocytes (F(5, 130) = .897, p = .486), thrombocytes (F(5, 130) = 1.027, p = .405)].
Combined Boxplots for Cortisol (T1) According to Subgroups with Differing Depressive Disorders and Healthy Controls
Combined Boxplots for Cortisol (T2) According to Subgroups with Differing Depressive Disorders and Healthy Controls
Combined Boxplots for Cortisol (T2-T1) According to Subgroups with Differing Depressive Disorders and Healthy Controls
Combined Boxplots for Cortisone (T1) According to Subgroups with Differing Depressive Disorders and Healthy Controls
Combined Boxplots for Cortisone (T2) According to Subgroups with Differing Depressive Disorders and Healthy Controls
Combined Boxplots for Cortisone (T2-T1) According to Subgroups with Differing Depressive Disorders and Healthy Controls
References
Subject Area
Reviews and Context
0
Comment
0
TRIP Peer Reviews
0
Community Reviews
0
Automated Services
0
Blogs/Media
Author Videos