Depression is a major mental health concern that is becoming a more widespread global issue. In 2015, the World Health Organization estimated that approximately 4.4% of the world’s population, or approximately 322 million people, were affected by depression1. Additionally, this report found that a similar percentage of the Thai population, or approximately 2.8 million people, suffers from depression1.

There is now evidence that major depressive disorder (MDD) and bipolar disorder (BD) are associated with immune-inflammatory pathway activation. Studies have demonstrated that inflammatory markers, such as cytokines, chemokines, complement factors, and acute phase reactants, are elevated in depressed individuals2,3,4,5. Both MDD and BD are characterized by activation of the immune-inflammatory response system (IRS) as indicated by elevated levels of pro-inflammatory cytokines, such as interleukin (IL)-1β, IL-6, IL-8 (CXCL8), IL-12, interferon (IFN)-γ, and tumor necrosis factor (TNF)-α6, and the compensatory immune-regulatory system (CIRS), which prevents hyperinflammation by downregulating the IRS and producing anti-inflammatory cytokines, such as IL-4, IL-10 and transforming growth factor (TGF) -β6. Under physiological conditions, IRS and CIRS interact in a balanced manner6. During acute mood episodes, such as depressive episodes and mania, this balance is disrupted, resulting in an increase in the production of IRS and CIRS cytokines and a new IRS/CIRS setpoint6,7,8. This is significant because certain IRS and CIRS cytokines have neurotoxic effects that may cause functional injury to neuronal brain circuits, particularly neuronal and astroglial projections9.

As is the case with IRS/CIRS cytokines, there is evidence that imbalances in the number of T helper (CD4+) and T cytotoxic (CD8+) T cells and T regulatory (Treg) cells may contribute to the development of severe mood disorders10,11,12. Flow cytometry-based immunophenotyping is a method for identifying activated T cells that contribute to IRS and CIRS profiles by producing pro- and anti-inflammatory cytokines. This technique enables us to distinguish between various categories of activated T effector (Teff) and Treg cells. Activated CD4+ and CD8+ cells that express activation markers such as CD69+, CD71+, CD40L+ (CD154+), and HLA-DR+ contribute to the IRS10,13,14 by producing pro-inflammatory cytokines, including IL-2, IL-17 and IFN-γ6,7,10,13,14. Activated Treg cells, including CD4+ CD25+ FoxP3+ cells that can express the cannabinoid receptor CB1, CD152 (CTLA-4) or glycoprotein A repetitions predominant (GARP)7,10,13,14,15, promote CIRS by regulating the immune response and/or generating anti-inflammatory cytokines such as IL-10 and TGF-β16. By modulating the production and secretion of pro-inflammatory cytokines17, Treg cells play a crucial role in preserving immune balance and promoting immune tolerance. In animal studies, Treg depletion may amplify immune-inflammatory pathways and result in autoimmunity18.

Using flow cytometry results, early precision medicine studies demonstrated that major depression, particularly melancholia, is qualitatively distinct from controls and minor depressed patients19. The most significant CD markers were elevated CD25+ and HLA-DR+ expression on CD4+ T cells, indicating T cell activation in depression and particularly melancholia19. In a separate study, Maes et al.20 determined that 64% of patients with MDD had increased expression of CD7+ CD25+ and CD2+ HLADR+ cells with a specificity of 91%. According to these findings, numerous patients with MDD exhibit T cell activation. The soluble transferrin receptor (sTfR or sCD71) and sIL-2R (CD25+) are elevated in the serum of patients with MDD and BD21,22,23, corroborating these findings.

Through their involvement in modulating cytokine production, Treg cells have been found to potentially impact mood disorders. Individuals with MDD have lower Treg levels than healthy individuals17, and antidepressant-treated MDD patients have a greater number of CD4+ CD25+ and CD4+ CD25+ FOXP3+ Treg cells24. The remitted phase of BD is characterized by a suppression of Teff cells and an activation of Treg cells, while increasing severity of BD, as measured by duration or number of manic episodes, is associated with Teff cell and Treg cell aberrations10. Mice with depleted CD4+ CD25+ Treg cells exhibit elevated levels of despair behaviors and decreased 5-HT in the hippocampus25.

Nevertheless, no studies have examined whether activation of T cells and depletion of Treg cells are associated with increased severity of depression in an acute phase of MDD and whether MDD is characterized by very early (CD69+), early (CD71+), late (HLADR+) activation markers, and/or increases in CD40L+, a key player in T-cell dependent effector functions and humoral immunity26,27,28,29. Consequently, the purpose of this study is to quantify the levels of CD69, CD71, CD40L an HLADR- bearing CD3 (pan T), CD4 and CD8 cells, and CB1, CD152 and GARP-bearing CD25+ FoxP3+ Treg cells in the acute phase of MDD, and to determine if activation of T cells and depletion of Treg are associated with the severity of the acute phase. In addition, the relationships between these T cell subsets and immune cell profiles, based on cytokine production, are assessed.



The investigation included participants from Bangkok’s Chulalongkorn University. The Department of Psychiatry’s senior psychiatrists recruited outpatients with MDD. Posters and word-of-mouth were used to recruit healthy controls. Before taking part in the study, all participants were required to sign a written consent form. Once participants were recruited, they were required to complete standard questionnaires and provide blood samples. Standard questionnaires and a semi-structured interview were utilized to collect demographic data, including age, gender, and body mass index. Patients between the ages of 18 and 65 who comprehended Thai, had been diagnosed with MDD by a psychiatrist using the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria, and had a Hamilton Depression Rating Scale (HAMD) score greater than 17, indicating moderate to severe depression, met the inclusion criteria. Patients with other DSM-5 axis 1 disorders, including schizophrenia, schizoaffective disorders, obsessive compulsive disorder, post-traumatic stress disorder, psycho-organic disorders, or substance misuse disorders, were excluded from the study. Healthy volunteers were recruited from the same catchment area as family or friends of staff members or friends of patients. We excluded individuals with a diagnosis of any axis 1 DSM-5 disorder or a positive family history of MDD, bipolar disorder, psychosis or suicide. Patients and controls were ineligible if they had experienced any allergic or inflammatory responses in the previous three months, if they had neuroinflammatory, neurodegenerative, or neurological disorders such as epilepsy, Alzheimer’s disease, multiple sclerosis, or Parkinson’s disease, if they had (auto)immune diseases such as COVID-19 infection (lifetime), chronic obstructive pulmonary disease, inflammatory bowel disease, psoriasis, diabetes type 1, asthma, or rheumatoid arthritis, had a history of receiving immunomodulatory drugs, had taken therapeutic doses of antioxidants or omega-3 polyunsaturated fatty acid supplements within three months before the study, had used anti-inflammatory drugs (NSAID or steroids) within one month of the study. We did not include lactating or expectant women. Some patients were taking psychotropic medications, such as sertraline (18 patients), other antidepressants (8 patients, such as fluoxetine, venlafaxine, escitalopram, bupropion, and mirtazapine), benzodiazepines (22 patients), atypical antipsychotics (14 patients), and mood stabilizers (4 patients). In the statistical analysis, the potential effects of these drug variables were taken into consideration.

Before participating in the study, all controls and patients gave their written informed consent. The study was carried out following both international and Thai ethics and privacy regulations. The Institutional Review Board of the Faculty of Medicine at Chulalongkorn University in Bangkok, Thailand approved the study (#528/63), which was in accordance with the international guidelines for human research protection, such as the Declaration of Helsinki, the Belmont Report, CIOMS Guideline, and the International Conference on Harmonization in Good Clinical Practice (ICH-GCP).


A research assistant conducted a semi-structured interview and an experienced psychiatrist administered the 17-item version of the HAMD to assess the intensity of depressive symptoms30. The Mini-International Neuropsychiatric Interview (M.I.N.I.) was used to make the diagnosis of psychiatric disorders31.

After an overnight fast, between 8:00 and 9:00 a.m., participants’ blood (20 mL) was collected. The blood was collected in BD Vacutainer® EDTA (10 mL) and BD Vacutainer® SST™ (5 mL) tubes provided by BD Biosciences (Franklin Lakes, NJ, USA). The serum-separating tubes were left to clot at room temperature for 30 min to obtain serum. The tubes were spun at 1100×g for 10 min at 4 °C. Peripheral blood mononuclear cells (PBMCs) were separated from the blood sample through density gradient centrifugation (30 min at 900×g) using Ficoll® Paque Plus (GE Healthcare Life Sciences, Pittsburgh, PA, USA). Utilizing a hemocytometer and trypan blue, 0.4% solution, pH 7.2–7.3 (Sigma-Aldrich Corporation, St. Louis, Missouri, United States), the cell count and viability were determined. Counting total and blue-stained cells ensured that greater than 95% of cells were viable under all conditions. To activate PBMCs, the 96-well plates were coated with 5 µg/mL of the anti-human CD3 antibody (OKT3, from eBioscience), overnight. 3 × 105 PBMCs were added to each well along with 5 μg/mL of the anti-human CD28 antibody (CD28.2, eBioscience). The cells were cultured in RPMI 1640 medium with L-glutamine supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin (Gibco Life Technologies, Carlsbad, CA, USA) for 3 days at 37 °C in an incubator with 5% CO2. The unstimulated PBMCs cultured for the same period were used as a negative control. After 3 days, the lymphocyte immunophenotypes were determined through flow cytometry.

To study lymphocyte immunophenotypes, 3 × 105 PBMCs were labeled to surface markers with monoclonal antibodies for 30 min, including CD3-PEcy7, CD4-APCcy7, CD8-APC, CD40L-FITC, CD69-AF700, CD71-PerCPcy5.5, HLA-DR-PE594, CD25-APC, CD152 PE Dazzle594, GARP-PE, CB1-AF700 and 7-AAD (Biolegend, BD Biosciences and R&D Systems). For Treg cells, we firstly stained cells with surface markers (CD3, CD4, CD25) followed by intracellular FoxP3 staining, which was performed using the FoxP3/Transcription Factor Staining Buffer Set (eBioscience) for fixation and permeabilization before staining with antibody to Foxp3-FITC (Biolegend). Electronic Supplementary File (ESF1), Figs. 1 and 2 show our gating strategies. The MFI of the activation markers expression was defined as the MFI of the activation markers expression (e.g. CD71+, HLA-DR+) gated within the CD4+ /CD8+ subpopulations. The flow cytometry was performed using LSRII flow cytometer (BD Biosciences) to evaluate lymphocyte immunophenotypes. All data were analyzed using FlowJo X software (Tree Star Inc., Ashland, OR, USA). Our study focused on characterizing different types of T cells. We did not utilize a Fc receptor blocking agent in our study. Table 1 of the ESF2 shows the different CD and cell surface markers, as well as their key functions, measured in the present study. We assessed both the percentages of T cells and the median fluorescence intensity (MFI) of the markers.

As previously described15, we measured cytokines and chemokines in stimulated diluted whole blood culture supernatant using the same blood samples as those employed for the flow cytometry. We used RPMI-1640 medium supplemented with L-glutamine, phenol red, and 1% penicillin (Gibco Life Technologies, USA), with or without 5 µg/mL PHA (Merck, Germany)+ 25 µg/mL lipopolysaccharide (LPS; Merck, Germany). On sterile 24-well plates, 1.8 mL of each medium was combined with 0.2 mL of 1/10-diluted whole blood. Using the Bio-Plex Pro human cytokine 27-plex assay kit (BioRad, Carlsbad, California, United States of America) and the LUMINEX 200 instrument (BioRad, Carlsbad, California, United States of America), the unstimulated and stimulated levels of cytokines and chemokines were measured. The intra-assay CV values were less than 11%. ESF2, Table 2 shows the cytokines and chemokines that were measured to compute the immune profiles15. ESF2, Table 3 displays how the immune profiles were computed as z unit composite scores using the cytokine and chemokine assessments15. In the current study, we examined whether the immune cell percentages or MFI values could predict the stimulated M1, Th-1, IRS, T cell growth, and neurotoxicity profiles.


Among study groups, we compared nominal variables via analysis of contingency tables (χ2-test) and continuous variables via analysis of covariance (ANCOVA). We employed Pearson’s product moment correlation coefficients for the purpose of examining the relationships between scale variables. A pre-specified generalized estimating equation (GEE), repeated measures, included fixed categorical effects of diagnosis (differences between patients and controls), cell types (T cell populations, namely CD69, CD71, CD40L, and HLADR-bearing CD3+, CD4+ and CD8+ cells), and the responsivity of the cells to in vitro administration of anti-CD3/CD28 antibodies (unstimulated versus stimulated condition), while allowing for the effects of age, sex, the drug state, BMI and smoking. Using pairwise comparisons, predefined comparisons were examined. Prespecified GEE, repeated measures, were used to examine the fixed categorical effects of diagnosis, treatment (unstimulated versus stimulated), cell type, and their interactions on Treg cells (CD25+ FoxP3+ CB1+, CD25+ FoxP3+ CD152+, and CD25+ FoxP3+ GARP+), while allowing for the effects of the aforementioned background variables. The primary statistical analyses examined the effects of Teff and Treg subsets on the HAMD score using multiple regression analyses. These analyses are the most important because a) a quantitative outcome score such as the HAMD is much more informative or correct than the binary diagnosis MDD32,33, and b) this method allows to examine the combined effects of both activated T and Treg cells on severity of illness. Furthermore, the effects of immune cell numbers or MFI levels on the stimulated production of cytokine profiles (M1, Th-1, IRS, T cell growth and neurotoxicity) were examined using multiple regression analysis with the stimulated production as dependent variables, and GEE analysis, repeated measures, with both unstimulated and stimulated immune profiles as dependent variables. Using manual multiple regression analysis, the effects of cell populations on the HAMD score were evaluated. We utilized an automatic forward stepwise regression strategy with a p-to-enter of 0.05 and a p-to-remove of 0.06 to determine which variables would be included and which would be excluded from the final regression model. We estimated the standardized β coefficients with t-statistics and exact p-values for each explanatory variable, F-statistics (and p-values), and partial eta squared (effect size) of the model. Using the modified Breusch-Pagan test and the White test, heteroskedasticity was investigated. We examined the probability of multicollinearity and collinearity utilizing the tolerance (cut-off value: 0.25) and variance inflation factor (cut-off value: > 4), as well as the condition index and variance proportions from the collinearity diagnostics table. Residuals, residual plots, and data quality were always evaluated in the final model. We also computed partial regression analyses, including partial regression plots, based on the results of the linear modeling analyses. When required, we utilized transformations such as Log10, square-root, rank-based inversed normal, and Winsorization to normalize the data distribution of the indicators. At p = 0.05, all analyses were two-tailed. To conduct the aforementioned statistical analyses, we utilized IBM’s Windows version of SPSS 28. Given an effect size of 0.25 (equivalent to 20% explained variance), alpha = 0.05, power = 0.8, and 3 covariates, an a priori power analysis (G*Power for a linear multiple regression analysis (the primary analysis) suggests that the minimum sample size should be 48.

Ethical statement

All subjects gave their written informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki of 1975, revised in 2013, and the protocol was approved by the Ethics Committee of the Institutional Review Board of the Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand (#528/63).


Demographic data

Table 1 displays the sociodemographic and clinical characteristics of the study’s patients and controls. There were no significant differences between the study groups in terms of age, gender distribution, education, or smoking. Patients had a higher BMI than controls. In any case, we have controlled all results for possible effects of age, BMI, sex, smoking, and the drug state of the patients but could not find any significant effects. The average HAMD scores of patients were substantially higher than those of controls, indicating that the majority of patients suffered from moderate to severe clinical depression. The outcomes of GEE analyses with immune profiles as dependent variables and stimulation status (baseline vs LPS+ PHA treatment), diagnosis (MDD versus controls), and stimulation by diagnosis interaction as predefined fixed effects are presented in the same table. The Th-1, IRS, T cell growth, and neurotoxicity (but not M1) profiles were all higher in MDD individuals compared to controls.

Table 1 Demographic and clinical data of the depressed patients and healthy controls (HC) included in the present study.

Baseline and stimulated Teff frequencies and MFI in MDD

Table 1 shows that there were no significant differences in CD3+, CD4+, and CD8+ percentages between MDD patients and controls. GEE analysis performed on the frequency of activated T cells showed significant effects of diagnosis (W = 4.78, df = 1, p = 0.029), stimulation (Wald = 596.28, df = 1, p < 0.001)), cell type (W = 1121.92, df = 1, p < 0.001)), diagnosis X cell type (W = 30.59, df = 11, p = 0.001), stimulation X cell type (W = 542.25, df = 1, p < 0.001), and the three-way interaction cell type X diagnosis X stimulation (W = 25.72, df = 11. P = 0.007). There were no significant effects of age, sex, BMI, smoking or drug state of the patients. Table 2 shows the frequencies of the T cell subtypes according to diagnosis and stimulation status. Anti-CD3/CD28 stimulation significantly increased the frequencies of all subsets. Activation markers were significantly higher in MDD (29.34 ± 0.98%) than in controls (27.6 ± 0.93%). ESF2, Table 4 shows the measurements of the 4 activation markers in CD3+, CD4+ and CD8+ T cells in baseline and stimulated conditions in MDD patients and controls. Pairwise comparisons showed that the number of activated T cells in the baseline condition was greater in patients than controls (p = 0.041). Table 2 shows that in the baseline condition, MDD patients showed higher frequencies of baseline CD3+ CD71+, CD3+ CD40L+, CD4+ CD71+, CD4+ CD40L+, CD4+ HLA-DR+ and CD8+ CD40L+ than normal controls. After stimulation with anti CD3/CD28 antibodies, there were no significant differences in any of the cell types between MDD and controls, except CD8+ CD40L+ which was elevated in MDD.

Table 2 Differences in unstimulated (UNST) and stimulated (STIM) changes in the percentage of T lymphocyte populations in healthy controls (HC) and major depressed patients (MDD).

ESF2, Table 5 shows the outcome of a similar GEE analysis performed on the MFI values of the activated T cells. This GEE analysis showed significant effects of anti-CD3/CD28 stimulation (Wald = 825.31, df = 1, p < 0.001), cell type (W = 1875.74, df = 1, p < 0.001), diagnosis X cell type (W = 23.00, df = 11, p = 0.018), stimulation X cell type (W = 1133.54, df = 1, p < 0.001), and the three-way interaction cell type X diagnosis X stimulation (W = 24.12, df = 11. P = 0.012). ESF2, Table 5 displays that MDD patients exhibit increased baseline MFI values of CD3+ CD71+, CD4+ CD71+, and CD4+ HLADR+ and a trend toward increased CD3+ HLADR+ MFI values as compared with controls. In addition, pairwise comparisons showed significant increases in stimulated CD3+ CD71+ and CD4+ HLADR+ MFI values between MDD patients and controls (p < 0.05). There were no significant effects of age, sex, BMI, smoking or drug state of the patients on any of the above measurements (these putative confounders were entered as covariates in the analyses).

GEE analysis performed on the frequencies of the three CD25+ FoxP3+ subtypes showed a significant effect of diagnosis (W = 3.96, df = 1, p = 0.047), stimulation (W = 211.51, df = 1, p < 0.001), cell type (W = 529.81, df = 1, p < 0.001) and cell type X simulation (W = 243.14, df = 1, p < 0.001). The CD25+ FoxP3+ cells were significantly lower in MDD (z score ± SE: -0.073 ± 0.055) than in controls (0.111 ± 0.074). As can be seen in Table 3, in vitro administration of anti-CD3/CD28 antibodies increased the number of all three FoxP3 cell types significantly. There was a significant group (diagnosis) by time (unstimulated versus stimulated condition) interaction for CD25+ FoxP3+ CB1, indicating that patients with MDD showed significantly lower CD25+ FoxP3+ CB1+ levels than controls in the unstimulated condition. The results of GEE, repeated measures, performed on the MFI of the CD25+ FoxP+ cells showed a significant effect of stimulation (W = 75.58, df = 1, p < 0.001) and cell types (W = 39.52, df = 1, p < 0.001). ESF2, Table 6 shows the interaction between diagnosis and stimulation with anti-CD3/CD28. Overall, stimulation significantly increased the MFI of FoxP3 cells. Pairwise comparisons show that the unstimulated expression of FoxP3 was significantly lower in MDD as compared with controls (p = 0.027), whereas there were no significant intergroup differences in the stimulated condition (p = 0.772).

Table 3 Differences in unstimulated (UNST) and stimulated (STIM) changes in the percentage of T regulatory (Treg) lymphocyte populations in healthy controls (HC) and major depressed patients (MDD).

Regression analyses of severity of illness on T cell subtype percentages and MFI values

Consequently, we performed the primary outcome statistical analyses, namely multiple regression analyses with the HAMD as outcome variable and the different cell populations (either frequencies or MFI) as explanatory variables. Table 4, regressions #1 and #2 show the outcome of these regressions using the percentages of the cell types, whereas regressions #4, #5 and #6 show the outcome of regressions on the MFI values. We examined the effects of the unstimulated and stimulated cell types. Regression #1 shows that 37.9% of the variance in the HAMD could be explained by unstimulated CD25+ FoxP3+ CB1% (inversely), unstimulated CD4+ CD40L+ % and stimulated CD4+ CD69+ % (positively). Figure 1 shows the partial regression of the total HAMD score on the baseline number of CD4+ CD40L+ cells. Nevertheless, also CD4+ CD71+ % was significantly associated with the HAMD score as shown in Fig. 2. Regression #2 (we omitted the CD4+ populations) shows part of the variance in the HAMD score could be explained by the regression on CD8+ CD40+ % (positively) and age (inversely).

Table 4 Results of multiple regression analyses with the Hamilton Depression rating Scale (HAMD) score as a dependent variable.
Figure 1
figure 1

Partial regression of the total Hamilton Depression Rating Scale (HAMD) score on baseline number of CD4+ CD40L+ cells.

Figure 2
figure 2

Partial regression of the total Hamilton Depression Rating Scale (HAMD) score on baseline number of CD4+ CD71+ cells.

Regarding the MFI data, we found that (regression #3) 35.5% of the variance in the HAMD could be explained by 3 stimulated MFI values (CD8+ CD40L+, CD3+ HLADR+, CD4+ CD69+, all positively) and age (inversely). Regression #4 shows that this prediction could be slightly improved by adding the unstimulated CD25+ FoxP3+ CB1 MFI data. Regression #5 displays that unstimulated CD25+ FoxP3+ GARP+ MFI was inversely associated with the HAMD.

Regression analyses of immune profiles on T cell subtype percentages and MFI values

ESF2, Table 7 shows the measurements of the unstimulated and stimulated M1, Th-1, IRS, T cell growth and immune-associated neurotoxicity profiles in MDD patients and controls. The stimulated M1, IRS, T cell growth and neurotoxicity profiles were significantly increased in MDD than in controls, whereas there was a trend towards an increased Th-1 profile. Table 5 shows the results of multiple regression analyses with stimulated immune profiles (M1, Th-1, IRS, T cell growth and neurotoxicity) as dependent variables and T cells, either percentage (regressions 1–5) or MFI (regressions 6–10) as explanatory variables, while allowing for the effects of extraneous variables. Regression #1 shows that 22.4% of the variance in M1 could be explained by CD8+ CD69+ % (positively) and CD25+ FoxP3+ CD152+ % (inversely). Th-1 was best predicted by CD25+ FoxP3+ CD152+ % (inversely). Regression #3 shows that a larger part (28.8%) of the variance in IRS was explained by baseline CD25+ FoxP3+ CD152+ % (inversely) and stimulated CD4+ CD69+ % and CD4+ CD71+ % (positively). Regression #4 shows that the T cell growth profile was predicted by three different cell types, namely CD3+ CD71+ %, CD4+ CD69+ % (positively) and CD25+ FoxP3+ CD152+ (inversely). Neurotoxicity was best predicted by CD25+ FoxP3+ GARP+ % (inversely).

With regard to the MFI values, we detected that M1 (regression #6) was best predicted by CD25+ FoxP3+ CD152+ MFI (inversely) and CD8+ CD69+ MFI (positively). Regression #7 shows that 12.9% of the variance in Th-1 was predicted by CD4+ CD71+ MFI (positively). IRS was strongly associated with CD8+ CD71+ MFI (regression #8). T cell growth was predicted (33.8% of the variance) by CD3+ CD71+ MFI, CD4+ CD69+ MFI (positively), and CD25+ FoxP3+ CD152+ MFI (inversely). Regression #10 shows that 24.4% of the variance in neurotoxicity was explained by the unstimulated CD25+ FoxP3+ CD152+ MFI (inversely) and CD8+ CD40L+ MFI values.

Table 6 shows that the changes from baseline to the LPS+ PHA-stimulated IRS and neurotoxicity scores were predicted by lowered numbers or expression of CD25+ FoxP3+ CD152+ T reg cells combined with baseline or stimulated CD8+ and CD4+ T activation subsets (either prevalence or MFI values).

Table 5 Results of multiple regression analyses with immune profiles as dependent variables and T cell subsets or expression markers as explanatory variables.
Table 6 Results of general equating estimation, repeated measurement.


Associations between T cells and immune profiles

The first major finding of this study is that there is a strong correlation between the flow cytometric assessments of T cell subtypes and cytokine-based immune-inflammatory profiles, including M1, Th-1, IRS, T cell growth, and neuroimmunotoxic profiles. Thus, these immune profiles were a) positively associated with the number or expression of diverse activated T cells, primarily CD4+ CD71+ and CD8+ CD69+ T cells, but also CD4+ CD69+ T, CD8+ CD71+ and CD8+ CD40L+ cells, and b) negatively associated with the number or expression of diverse Treg cells, including CD152- and GARP-bearing CD25+ FoxP3+ cells.

During T cell activation, distinct phases can be distinguished by the presence of distinct markers. CD69 is a cell surface marker that is expressed immediately after T cell activation and rapidly disappears after the initial stimulation14. CD69 is a costimulatory molecule for T cell differentiation, activation, and proliferation34 as well as lymphocyte retention in tissues14. CD69 activation stimulates the production of IL2 and TGF-β1 and controls the production of pro-inflammatory cytokines such as IL-17, IFN-γ, and IL-22, as well as Treg cells14. CD40L is present in the early phases of T cell activation35 and binds to CD40 on macrophages, CD8+ T cells, B cells, and dendritic cells, thereby promoting the differentiation and proliferation of immune cells, immune-inflammatory responses, and activation and maturation of B cells, which ultimately results in the production of autoantibodies35. After antigen presentation, preformed CD40L molecules are mobilized from the lysosomes of Th-1 cells36, while colligation of CD3 and CD40L increases the production of IFN-γ, TNF-α, and IL-1026. During a later phase of T cell activation, CD71 (TfR) is upregulated, resulting in an increase in iron uptake into the cells37, which stimulates T cell activation and proliferation38. Within twenty-four to forty-eight hours of T cell activation, the late phase activation marker HLA-DR is expressed. This latter molecule is a class II MHC, and its upregulation is associated with an increase in IFN-γ production39.

Conversely, the inverse associations of the immune profiles with CD152 (CTLA-4) and GARP-bearing Treg cells indicate immunoregulatory Treg effects on different aspects of the immune response. Depending on the surface expression of the proteins, Tregs regulate the immune system through a variety of mechanisms. CTLA-4 also known as CD152 is expressed on Treg cells (and activated T cells) and functions as a negative immune checkpoint that regulates immune responses40. By competing with CD28 for binding to the same ligands on antigen-presenting cells, CD152 inhibits the activation and proliferation of T cells41. GARP is a Treg cell surface molecule that indicates Treg activation and the release of TGF-β1, an anti-inflammatory cytokine. As a result, Treg cells promote immune tolerance41,42. Treg cells, may express CB1, which inhibits the production of pro-inflammatory cytokines and exerts negative immunoregulatory effects43. Previously, we discovered that the number of CD25+ FoxP3+ CB1+ cells was inversely associated with IRS responses15. However, after introducing the effects of GARP- and CD152-bearing Tregs, the effects of CB1 were no longer significant. Intriguingly, the current study detected that decreased numbers of baseline Treg cells predict increased neurotoxic potential, indicating that diminished fitness of Treg cells is associated not only with exaggerated immune responses, but also with increased neuroimmunotoxicity.

T cell activation in MDD

The second major finding of this study is that the numbers of baseline CD3+ CD71+, CD3+ CD40L+, CD4+ CD71+, CD4+ CD40L+, CD4+ HLADR+, and CD8+ CD40L+ and the expression of CD3+ CD71+, CD4+ CD71+, and CD4+ HLADR+ are significantly greater in MDD than in controls. Furthermore, the study demonstrates a strong correlation between the severity of depression and elevated numbers of CD4+ CD40L+ and CD4+ CD69+ cells, as well as increased expression of CD40L on CD8+ cells and HLADR on CD3+ cells.

These results suggest that patients with MDD have an abnormal distribution of CD4+ and CD8+ T cell subsets with elevated T cell activation markers. Maes et al. demonstrated in previous investigations that depressed patients had a higher expression of T cell activation markers, specifically CD25+ and HLADR+ 20. A recent meta-analysis also reported an increase in the mean absolute number of activated T cells expressing CD25+ and CD3+ HLADR+ 44. Moreover, Maes et al.10 found that in the symptomatic remission phase of BD, the frequency of unstimulated CD3+ CD8+ CD71+ cells is lower than in healthy controls. This suggests that the acute phase of major mood disorders is characterized by T cell activation and that this T cell activation is suppressed when the symptoms ameliorate. Therefore, future research on T cell activation (and the immune system in general) should always differentiate between the acute, partially remitted, and remitted phases of mood disorders.

Our findings indicate that the pathogenesis of the acute phase of severe MDD is associated with CD4+ and CD8+ T cell activation. Interaction with the antigen-MHC complex activates CD4+ T cells, which consequently differentiate into distinct subtypes under the influence of cytokines39. Almulla et al.45 reported that increased IL-16, Th-1 activation and Th-1 polarization are hallmarks of the acute phase of severe MDD. IL-16 signaling occurs via the CD4 molecule on Th cells, activating CD4+ cells and upregulating activation markers and the IRS46,47,48. Ligation of CD40 via CD40L stimulates the secretion of the Th-1-polarizing cytokine IL-1249,50. Moreover, CD40L expression by CD4+ cells is a crucial step in the activation of CD8 T cells29, which may express CDL40, perform T helper functions, and promote their own growth27,29. In addition, CD40L-CD40 interactions mediate the effect of CD4+ cells on CD8+ cytotoxic T lymphocytes51.

Consequently, our findings suggest that CD4+ and CD8+ T cell activation, along with increased Th-1 polarization and IL-16 production, are key phenomena that drive increased T effector, T helper, and T cytotoxic functions, as well as enhanced M1, Th-17, IRS, and neurotoxic immune profiles in the acute phase of severe MDD. These activation markers may contribute to the immune, autoimmune, and neurotoxic processes seen in MDD. First, CD69 is involved in systemic lupus erythematosus (SLE), rheumatoid arthritis, atopic dermatitis, and systemic sclerosis, as well as animal models of arthritis, myocarditis, and inflammatory bowel disease14. Second, increased expression of CD71 on T cells is associated with SLE disease activity and an enhanced Th-17 profile, whereas CD71 may cause autoimmunity in an SLE model52. Targeting CD71 reduces autoimmunity and pathology in a mouse model of SLE and increases the secretion of IL-1052. Third, CD40L-CD40 interactions play a significant role in neurodegenerative and neurological disorders, such as Parkinson’s disease, Alzheimer’s disease, stroke, multiple sclerosis, and epilepsy53. Increased permeability of the blood–brain barrier (BBB), injury to neuronal and glial cells, neuroinflammation, and the formation of microthrombi are mediated by CD40L-CD40 signaling53. In an animal model of autoimmune disease, CD4+ immune cells, attracted and activated by IL-16, are associated with neuroinflammation54. CD40L+ T cells may cross the BBB and activate microglial CD40, resulting in the production of pro-inflammatory cytokines and reactive nitrogen species, which may cause demyelination of axons53. In addition, activated CD8+ cells may induce IRS responses and induce neurotoxicity55.

Breakdown of immune tolerance in MDD

The third major finding of this study is that lowered numbers/expression of CD3+ CD4+ CD25+ FoxP3+ GARP+ and CD3+ CD4+ CD25+ FOXP3+ CB1+ coupled with T cell activation, is inversely associated with the severity of illness. Jahangard et al. showed that the frequency of FOXP3+ Tregs was decreased in untreated MDD patients compared to healthy controls, while the proliferation of circulating CD4+ T cells was increased in the patients56. Ellul et al.57 found an association between increased risk of MDD and decreased Treg numbers in association with increased inflammation57. Importantly, during the remission phase of mood disorders, Treg functions may be upregulated10.

Treg cells and CIRS activities ensure immune homeostasis by regulating innate and adaptive immune effector cells and restraining the detrimental activities of these effector cells. The transcription factor FoxP3 is a key regulatory factor in Treg activities by producing CIRS cytokines, such as IL-10 and TGF-β1, depleting growth factors, and expressing the co-inhibitory CTLA-4 molecule58. Furthermore, CD25 signaling and the expression levels of FoxP3 are important for the function and stability of these cells59,60,61,62. Moreover, CD25+ (or IL-2R) Treg cells may absorb IL-2 precluding IL-2 for the growth of CD8+ T effector and Th-1 cells63. Depletion of Treg cells may be the consequence of reduced fitness, altered apoptosis, and increased inhibitory activities of CD7152.

Overall, the findings indicate that the acute phase of severe MDD is characterized by a breakdown of peripheral immune tolerance. Loss of Treg-associated immune tolerance can lead to immunopathologies and autoimmune diseases64. In addition, Treg cells have neuroprotective properties by suppressing neurotoxic T effector responses and production of neurotoxic cytokines, and favoring anti-inflammatory microglia polarization6,63. When the frequency of Tregs is decreased or their immunosuppressive effects are diminished, the Teff/Treg balance is shifted towards the production of pro-inflammatory Th-1 and Th-7 phenotypes63. As such, diminished Treg functions predispose to neuroinflammatory disorders including multiple sclerosis, stroke, and Alzheimer’s disease. These findings further support the IRS/CIRS theory of mood disorders and previous studies that have highlighted the imbalance of IRS and CIRS activity in patients with MDD57,65,66,67.


This study would have been more interesting if we had measured other immune-related biomarkers of depression including oxidative and nitrosative stress markers. We stimulated PBMCs with anti-CD3/anti-CD28 antibodies and detected Treg cells by gating CD3+ CD4+ CD25+ FoxP3+ cells. Nevertheless, activated conventional T cells can be induced to transiently express Foxp368. However, GARP is only expressed on Treg cells and upregulated in activated Treg cells. Moreover, GARP+ Treg cells clearly show suppressive functions directed at target cells42,69. In our study, we detected surface GARP on Treg cells (CD3+ CD4+ CD25+ FoxP3+ cells) that is likely separated from FoxP3-expression by conventional T cells upon stimulation. In addition, we gated CD25++ which marks Treg cells and excluded activated T cells as shown in our gating strategy figures (in ESF1). Our gating strategy as shown in the ESF1 indicates that no spillover effects were present. Moreover, our unstimulated samples were used to separate negative from positive effects and to determine the best gate placement for the stimulated activation markers70. In addition, the current study also used the fluorescence minus one (FMO) method to identify positive populations.


The acute phase of severe MDD is characterized by a breakdown of immune tolerance, and CD40L activation, which are associated with increased neuroimmunotoxic potential. Diminished Treg neuroprotection, T cell activation, Th-1 polarization, increased CD71+ and CD40L expression and IL-16 production are new drug targets to treat the T effector-associated neurotoxicity of the acute phase of severe MDD.

There is now some evidence that anti-inflammatory agents may be useful to treat MDD patients with an elevated immune profile71. A recent meta-analysis shows that different types of anti-inflammatory agents have some clinical efficacy in treating MDD72. Nevertheless, our results show that MDD should not be treated with repurposing anti-inflammatory drugs (e.g., Cox-2 inhibitors) and antioxidant supplements (e.g., curcumin), but rather should boost Treg functions and target CD40L or CD71 and IL-16 to attenuate T cell activation.