Abstract
Major depressive disorder and a major depressive episode (MDD/MDE) are characterized by activation of the immune-inflammatory response system (IRS) and the compensatory immune-regulatory system (CIRS). In MDD/MDE, recent precision nomothetic psychiatry studies discovered a new endophenotype class namely Major Dysmood Disorder, a new pathway phenotype namely reoccurrence of illness (ROI), and a new model of the phenome of depression. The aim of the present study is to examine the association between ROI, the phenome of depression, and MDMD’s features and IRS, CIRS, M1, T helper (Th)1, Th2, Th17, Tregulatory, and growth factor (GF) profiles. Culture supernatants of unstimulated and stimulated (5 μg/mL of PHA and 25 μg/mL of LPS) diluted whole blood of 30 MDD/MDE patients and 20 controls were assayed for cytokines/GF using the LUMINEX assay. MDMD was characterized by increased M1, Th1, Th2, Th17, Treg, IRS, CIRS, neurotoxicity and GF profiles. Factor analysis shows that ROI features and immune-GF profiles may be combined into one a new pathway phenotype (an extracted latent vector). ROI, lifetime and recent suicidal behaviors, and severity of depression are significantly associated with immunotoxicity and GF profiles. Around 80.0% of the variance in the phenome is predicted by ROI and neurotoxicity or the IRS / CIRS ratio. The molecular pathways underpinning ROI-associated sensitization of immune/growth networks are transmembrane receptor protein kinase triggered STAT protein phosphorylation, TLR/NF-κB, JAK-STAT, and the main proliferation/survival PI3K/Akt/RAS/MAPK pathway. In conclusion, MDMD’s heightened immune responses are the consequence of ROI-associated sensitization combined with immunostimulatory triggers.
Introduction
Cytokines have been implicated in major depressive disorder (MDD) and a major depressive episode (MDD) in bipolar disorder (BD) in recent reviews and meta-analyses [1,2]. Maes et al. discovered increased levels of cytokine or cytokine receptor production in the culture supernatant of stimulated peripheral blood mononuclear cells (PBMCs) or in the serum of MDD and BD in the 1990s, including interleukin(IL)-2, soluble IL-2 receptor (sIL-2R), IL-1, the sIL-1R antagonist (sIL-1RA), IL-6, sIL-6R, and tumor necrosis factor (TNF)-α [review: 3]. These findings were corroborated by the same laboratory’s findings that MDD and BD are associated with an inflammatory or acute phase (APR) response characterized by elevated levels of positive AP reactants (including haptoglobin) and complement factors, and decreased levels of negative AP reactants (including albumin) [3]. Additionally, early machine learning findings indicated that clinical depression is associated with an increase in the expression of T cell activation markers including CD25+ (IL-2R) [3]. The cytokine, monocyte-T lymphocyte, and immune-inflammatory responses system (IRS) theories of depression were developed in 1995 because of these observations [3].
According to a recent review [1], MDD and BD are associated with activation of the IRS and the compensatory immune regulatory system (CIRS), which includes increased levels of immunoregulatory products that downregulate the IRS and prevent hyperinflammation, such as sIL-1RA, sIL-2R, IL-4, IL-10, and some AP proteins [1]. MDD and MDE are both characterized by activation of macrophage M1 (IL-1, IL-6, TNF-α), T helper (Th)1 (IL-2, IFN-γ), Th2 (Il-4, IL-5), Th17 (IL-17), and T regulatory (Treg) cells (IL-10) [1]. Nonetheless, the IRS is more active than the CIRS during the acute phases of MDD/MDE, resulting in a net activation of the immune system. Moreover, also growth factors including platelet derived growth factor (PDGF), vascular endothelial growth factor (VEGF), and fibroblast growth factor (FGF) are higher in MDD/MDE than in controls [4;5]. It is critical to note that an immune profile characterized by neurotoxic cytokines (M1, Th1, and Th17) is another hallmark of MDD/MDE leading to the conclusion that this profile may result in neuro-affective toxicity, resulting in the affective and cognitive symptoms of depression [1]. The immunological response is intimately linked to redox mechanisms including increased formation of reactive oxygen and nitrogen species (RONS), diminished antioxidant defenses, and indicators of enhanced nitro-oxidative stress toxicity (OSTOX), all of which are found in MDD and BD [6–8]. Suicidal behaviors (SB), both recent and lifetime suicidal ideation and attempts, are associated with activated IRS, CIRS, RONS, and OSTOX pathways [9].
The staging of illness (conceptualized as the reoccurrence of episodes and SBs) is significantly associated with IRS, CIRS, and redox pathways. Thus, in MDD/BD, the number of prior depressive and manic episodes, as well as SBs, is associated with an increase in plasma IL-1β, sIL-1RA, IL-6, TNF-α, and neopterin, lowered antioxidant enzyme activities and damage to lipids and proteins [10–13], while other findings indicate that immune-inflammatory responses are more pronounced in later stages of illness [14,15]. In BD, we detected that the frequency of episodes is inversely related to the proportions of stimulated antigen specific activated CD3+CD4+ T cells, and the expression of early activation markers and the transferrin receptor on CD4+ and CD8+ cells, whereas later stages of illness are characterized by decreased frequencies of activated Treg cells [16]. The association between staging and immunological pathways was explained by the fact that pro-inflammatory signals may sensitize the immune system, as well as by abnormalities in CIRS and proliferative responses [12,16].
In this respect, using a new precision nomothetic psychiatry approach, Maes et al. [17–20] established a) a new reoccurrence of illness index (ROI) of MDD/BD, namely a latent vector (LV) extracted from number of depressive and manic episodes and SB as well; b) a novel model of affective disorders based on antioxidant gene variants, adverse outcome pathways (immune-redox pathways), ROI, and the phenome, which was conceptualized as a LV extracted from depression, anxiety, clinical, SB, disability, and quality of life ratings; c) a new ROI-redox pathway phenotype, which was conceptualized as a LV extracted from ROI and redox biomarkers; and d) a new diagnostic class namely Major DysMood Disorder (MDMD), characterized by aberrations in immune-redox biomarkers, increased ROI, and phenome scores. Importantly, we found that MDMD cuts across unipolar depression and bipolar disorder and is more prominent than the latter diagnoses indicating that both are phenotypes of the same disorder [16,18,20]. Nevertheless, there are no data whether MDMD and increased phenome scores are accompanied by aberrations in IRS/CIRS and an increased immunoneurotoxicity and whether a ROI-immune pathway phenotype may be constructed by linking ROI and cytokines/growth factors.
In previous studies, we used not only serum/plasma but also 72 hours culture supernatant of unstimulated and LPS + PHA stimulated diluted whole blood to measure cytokine/growth factor production [4,21–23]. While the serum levels of some cytokines (IFN-γ, IL-2, IL-4, and IL-5) are difficult to measure, their ex vivo production in stimulated whole blood is well measurable [21–24]. Moreover, the latter method adequately reflects the in vivo cytokine production especially that of cytokines/growth factors, which are produced by different immune profiles [21–24]. Nevertheless, there are no data on the cytokine and growth factor profiles of MDMD, SBs and ROI either in unstimulated or stimulated diluted whole blood cultures.
Hence, the aim of the study is to examine a) whether a MDMD class may be constructed that is associated with greater (un)stimulated immune (M1, Th1, Th2, Th17, Treg, IRS, CIRS, T cell growth, growth factor, and neurotoxicity) profiles and a IRS/CIRS ratio; b) whether ROI is associated with the immune profiles and whether a ROI-immune pathway phenotype may be constructed; and c) whether ROI and associated IRS and neurotoxicity profiles predict the phenome of MDD/MDE.
Methods and Participants
Participants
In this study, we included 20 normal controls and 30 depressed patients recruited from the outpatient clinic of the King Chulalongkorn Memorial Hospital’s Department of Psychiatry in Bangkok, Thailand. According to DSM-5 criteria, the patients were diagnosed with MDD or MDE and had moderate to severe depression as measured by the Hamilton Depression Rating Scale (HDRS). We recruited normal volunteers of both sexes, ranging in age from 18 to 65 years from the same catchment area, namely Bangkok, Thailand. The controls were recruited by word of mouth. Both patients and controls were excluded if they a) had neuroinflammatory, neurodegenerative, or neurological disorders such as multiple sclerosis, epilepsy, Alzheimer’s disease, stroke, or Parkinson’s disease; b) had (auto)immune diseases such as chronic obstructive pulmonary disease, cancer, psoriasis, type 1 diabetes, asthma, and inflammatory bowel disease, c) had inflammatory or allergic reactions three months prior to the study; d) were treated with immunomodulatory drugs (lifetime history) including glucocorticoids); e) were treated with therapeutic doses of omega-3 or antioxidant supplements or anti-inflammatory medication the month prior to the study; f) were pregnant or lactating women. Other DSM-5 axis 1 illnesses such as psycho-organic disorders, schizoaffective disorders, schizophrenia, obsessive compulsive disorder, post-traumatic stress disorder, and drug abuse disorders were excluded as exclusion criteria for MDD/MDE patients. Healthy participants were excluded if they had a diagnosis of any DSM-5 axis 1 condition or a positive family history of MDD or BD. We statistically adjusted for the possible impacts of the drug state of the patients, namely sertraline (n=18), other antidepressants (n=8, including fluoxetine, venlafaxine, escitalopram, bupropion, and mirtazapine), benzodiazepines (n=22), atypical antipsychotics (n=14), and mood stabilizers (n=4). Prior to participating in this study, all controls and patients submitted written informed consent. The research adhered to international and Thai ethical standards and privacy legislation. The study was approved by the Institutional Review Board of Chulalongkorn University’s Faculty of Medicine in Bangkok, Thailand (#528/63), in accordance with the International Guidelines for the Protection of Human Subjects as required by the Declaration of Helsinki, The Belmont Report, the CIOMS Guideline, and the International Conference on Harmonization in Good Clinical Practice (ICH-GCP).
Clinical measurements
A research assistant with expertise in mood disorders performed semi-structured interviews. To assess the severity of depression symptoms, we utilized the HDRS, 17-item version, given by an expert psychiatrist [25]. The Thai state version of the State-Trait Anxiety Assessment (STAI) is a psychological inventory designed to assess the intensity of state anxiety [26]. To evaluate psychiatric axis-1 diagnoses, the Mini-International Neuropsychiatric Interview (M.I.N.I.) was utilized [27]. To compute the ROI, we registered the number of depressive and (hypo)manic episodes and assessed recent and lifetime suicidal behaviors (SB) using the Columbia-Suicide Severity Rating Scale (C-SSRS) lifeline version [28]. Lifetime suicidal behaviors were assessed as C-SSRS item 1 (lifetime suicidal ideation, namely wish to be dead), and item C-SSRS suicidal behavior, number of actual attempts. Recent SB was conceptualized as a PC (labeled “PC recent SB”) extracted from nine C-SSRS items, namely wish to be dead, non-specific active suicidal thoughts, active suicidal ideation with any methods, active suicidal ideation with some intent to act, active suicidal ideation with specific plan/intent, frequency and duration of suicidal ideation, actual attempts, and total number of actual attempts (all past month). Lifetime SB was conceptualized as a PC (labeled “PC lifetime SB”) extracted from 11 C-SSRS items, namely wish to be death lifetime, non-specific active suicidal thoughts, active suicidal ideation with any methods, active suicidal ideation with some intent to act, active suicidal ideation with specific plan/intent, frequency and duration of ideation, number of actual attempts, preparatory acts or behavior and total number of preparatory acts (all lifetime).
Assays
At 8:00 a.m., after an overnight fast (at least 10 hours), blood was collected in BD Vacutainer® EDTA (10 mL) tubes (BD Biosciences, Franklin Lakes, NJ, USA). In the present study, we measured cytokines/growth factors in unstimulated and stimulated (PHA and LPS) whole blood culture supernatant [21–23]. We used RPMI-1640 medium (Gibco Life Technologies, USA) supplemented with L-glutamine and phenol red and containing 1% penicillin (Gibco Life Technologies, USA) with or without 5 g/mL PHA (Merck, Germany) + 25 g/mL lipopolysaccharide (unstimulated) (LPS; Merck, Germany). 1.8 mL of each of these two mediums was added to 0.2 mL of whole blood, 1/10 diluted, on 24-well sterile plates. Whole blood was seeded on 24-well culture plates. Each subject’s specimens were separated into unstimulated and stimulated conditions and were incubated for 72 hours at 37°C, 5% CO2 in a humidified environment. After incubation, the plates were centrifuged for 8 minutes at 1500 rpm. Supernatants were carefully removed under sterile circumstances, split into Eppendorf tubes, and promptly frozen at −70°C until thawed for cytokine/growth factor assays. Electronic Supplementary File (ESF), Table 1 shows the names, acronyms, and official gene symbols of all cytokines/growth factors measured in the current study. ESF, Table 2 offers a list of the various immune profiles studied here.
Demographic and clinical data of the healthy controls (HC) and depressed patients divided into those with simple depression and Major Dysmood Disorder (MDMD)
Differences in unstimulated (UNST) and lipopolysaccharide + phytohaemagglutinin-stimulated (STIM) changes in various immune profiles in healthy controls (HC) and depressed patients divided into those with simple depression and Major Dysmood Disorder (MDMD)
The cytokines/growth factor were quantified using the LUMINEX 200 equipment (BioRad, Carlsbad, California, United States of America), a multiplex approach. In summary, supernatants were diluted fourfold with medium and incubated for 30 minutes with linked magnetic beads. After adding detection antibodies and streptavidin-PE for 30 minutes and 10 minutes, respectively, the fluorescence intensities (FI) were measured. We chose the (blank analyte subtracted) FI values in the current research for statistical analyses since FI are often a better option than absolute concentrations, particularly when numerous plates are employed [29]. The range of FI values that fall inside the concentration curve is shown in Table 1. All samples of all cytokines were quantifiable, except for IL-7, which displayed an abnormally high number of results below the assay’s sensitivity and was thus removed [26]. IL-13 demonstrated an acceptable lower limit of 30% of the sensitivity and so could be included in the statistical analysis. The CV values between analyses are less than 11% for all studies.
Statistical analysis
Analysis of variance (ANOVA) was used to compare scale variables, whereas chi-square tests or Fisher’s Freeman-Halton test were used to compare nominal variables across categories. We used generalized estimating equations (GEE) analysis to investigate the association between MDMD and the immune profiles and cytokines/growth factors. The pre-specified GEE analysis, which used repeated measures (unstructured working correlation matrix, linear scale response, and maximum likelihood estimation as scale parameter method), included fixed categorical effects of time (unstimulated versus stimulated), groups (MDMD versus simple depression and controls), and time x groups or time-by-continuous variable interactions (e.g. HDRS, STAI, staging, PC lifetime and recent SV) as well as sex, smoking, age, and BMI as covariates. The primary outcome variables in the GEE analyses were the immune profiles, and if these revealed substantial results, we additionally investigated the individual cytokines/growth factors. Multiple effects of time or group on the immune profiles were corrected using the false discovery rate (FDR) p-value [30]. Additionally, we incorporated the drug state of the patients as additional predictors in the GEE analysis to rule out any influence of these potential confounders. The GEE approach allows us to account for significant interactions and confounders while avoiding biased imputations induced by incomplete assessments. Nonetheless, there were no missing values in any of the demographic, clinical, or cytokine/growth factor data analyzed in this investigation (except for IL-7, which was excluded from the analysis). We calculated the estimated marginal means for the groups as well as the time X group interactions and used (protected) pairwise contrasts (least significant difference at p=0.05) to examine differences between groups and time x group interactions. We employed multiple regression analysis (automatic technique with a p-to-entry of 0.05 and a p-to-remove of 0.06 while evaluating the change in R2) to identify the biomarkers that predict the phenome or ROI scores. Multicollinearity was determined using tolerance and VIF and multivariate normality was determined using Cook’s distance and leverage, and homoscedasticity was determined using the White and modified Breusch-Pagan tests. The results of these regression analyses were always bootstrapped using 5.000 bootstrap samples, and the latter are presented if the findings were not concordant. Principal component analysis (PCA) was used to reduce the number of suicidal behaviors features and to summarize the information in summary indices or patterns. The first PC was considered to indicate a valid pattern when it explained at least 50.0% of the total variance and all loadings were greater than 0.6. To construct latent vectors (factors) underpinning several indicators we used exploratory factor analysis (unweighted least squares) with the same quality criteria as described above. In addition, we always estimated the factorability using the Kaiser-Meyer-Olkin test for sampling adequacy (should be >0.6) and Bartlett’s test of sphericity. In addition, we performed Partial least squares (PLS)-SEM analysis to derive latent variable (LV) scores only when the factors complied with prespecified quality criteria, namely all factor loadings should be >0.6 at p<0.0001, average variance extracted (AVE) >0.5, Cronbach alpha >0.7, composite reliability >0.8, and rho_A>0.8 [20]. All statistical analyses (except PLS) were conducted using IBM SPSS windows version 28. Two-tailed tests were used, and statistical significance was set at p<0.05. Using a two-tailed test with a significance threshold of 0.05 and assuming an effect size of 0.23 and a power of 0.80 while considering three groups and intercorrelations of around 0.6, the estimated sample size for a repeated measurement design ANOVA is approximately forty-two.
We developed seed-gene based protein-protein interaction (PPI) networks based on the differentially expressed proteins (DEPs) in MDMD versus controls. The networks were constructed using STRING version 11.0 (https://string-db.org), a predictive database, and IntAct Molecular Interaction Database (https://www.ebi.ac.uk/intact/), a database based on peer-reviewed articles. We developed zero-order PPIs (consisting entirely of seed proteins), first-order PPI network (consisting of 50 interactions in the first shell and none in the second shell; set organism: homo sapiens, and a minimum required interaction score of 0.400), and expanded networks including using OmicsNet (OmicsNet). Markov Clustering (MCL) analysis was performed using STRING to discover DEP communalities. STRING and the Cytoscape (https://cytoscape.org) plugin Network Analyzer were used to examine the network topology. The network’s backbone was defined as a collection of top hubs (nodes with the greatest degree) and non-hub bottlenecks (nodes with the highest betweenness centrality). STRING was utilized to show the physical interactions between the DEPs, and Gene Ontology (GO)net (GOnet (dice-database.org) to make graphs which contain GO terms and genes. The PPI networks were analyzed for their enrichment scores and annotated terms using the following tools: a) STRING to establish GO biological processes and molecular functions, diseases, as well as KEGG (https://genome.jp/kegg/) and WIKI (WikiPathways - WikiPathways) pathways; b) OmicsNet (using InAct,) for establishing REACTOME (European Bioinformatics Institute Pathway Database; https://reactome.org) and PANTHER (www.pantherdb.org/pathway/) pathways; Enrichr (Enrichr (maayanlab.cloud)) to establish Elsevier Pathways visualized using Appyter; and MetaScape (Metascape) to establishmolecular complex detection (MCODE) components based on GO terms. The enrichment analysis results are always shown using FDR-corrected p-values or q-values.
Results
Construction of ROI, phenome scores and the MDMD phenotype
The ROI was computed as the first LV extracted from number of depressive episodes, total number of episodes, C-SSRS lifetime suicidal ideation, C-SSRS number of lifetime suicidal attempts, and PC lifetime SB (the latter was computed as the first PC extracted from 11 lifetime SB C-SSRS items and this PC explained 62.21% of the variance while all the items showed loadings > 0.740). The first ROI LV explained 75.6% of the total variance and all loadings were >0.6, namely number of depressive episodes: 0.891, total number of episodes: 0.909, lifetime suicidal ideation: 0.931, number of lifetime suicidal attempts: 0.664, and PC lifetime SB: 923, and showed excellent quality criteria including Cronbach alpha: 0.916, composite reliability: 0.939, rho_A: 0.939, and AVE: 0.756. In order to compute the phenome score we a) computed the PC recent SB by extracting the first PC from the 9 recent SB C-SSRS items, which explains 60.54% of the variance while all those items showed loadings > 0.6; and b) computed the first LV extracted from the PC recent SB (0.862), diagnosis (controls: 0, MDD/MDE: 1, with psychotic/melancholia features: 2; loading: 0.906), HDRS (0.942), and STAI (0.805). This LV showed excellent quality criteria including Cronbach alpha: 0.902, composite reliability: 0.932, rho_A: 0.920, and AVE: 0.775. The combined ROI-phenome score was computed as the first LV (explaining 72.5% of the variance) extracted from all 5 ROI indicators and all four phenome indicators. This LV showed excellent quality criteria with factor loadings that were all higher than 0.646, and Cronbach alpha: 0.951, composite reliability: 0.959, rho_A: 0.957, and AVE: 0.725. We divided the study sample into three non-overlapping groups using the ROI-phenome score and a visual binning procedure (analysis of apparent modes and local minima in the frequency histogram with two cutoff points, namely −0.6 and 0.6). As such, we obtained three study groups, namely controls and depressed patients divided into these with lower (labeled: simple depression) versus high (labeled: MDMD) ROI-phenome scores.
Demographic and clinical data of the study groups
Table 1 shows the features of the normal controls, simple depression and MDMD. There were no significant differences in age, sex, education between the study groups, while depressed patients showed a somewhat higher BMI than controls. The MDMD group showed a higher prevalence of depression with melancholia/psychotic features, number of depression and all episodes, PC recent and lifetime SB scores.
Differences in immune biomarkers between MDMD and SD
Table 2 shows the measurement of unstimulated and stimulated immune profiles in controls, simple depression and MDMD. The time x group interactions of all immune profiles were significant, and these differences remained significant after FDR p-correction (at p<0.04). The interaction patterns showed that the stimulated production of all those profiles was significantly higher than the unstimulated production (all p<0.001) and that the production was always higher in MDMD than in controls while there were no differences between simple depression and controls. Moreover, the M1, Th1, T cell growth, growth factor, and neurotoxicity profiles were significantly higher in MDMD than in simple depression. In the GEE analyses there were no significant effects of age, sex, BMI and TUD. In addition, no significant effects of any of the drugs could be found on any of the immune profiles or single cytokines/growth factors even without FDR p-correction. For example, there were no significant effects of sertraline (W=1.01, df=1, p=0.314), other antidepressants (W=0.895, p=0.344), benzodiazepines (W=0.91, p=0.340), mood stabilizers (W=0.10, p=0.745) and atypical antipsychotics (W=0.59, p=0.443) on the neurotoxicity index.
Table 3 shows the secondary analyses performed on all separate cytokines/growth factors. We found significant interaction patterns between time x group with significantly increased stimulated production of sIL-1RA, IL-5, CXCL-8, IL-9, IL-15, IL-17, FGF, IFN-γ, CXCL10, PDGF, CCL5, TNF-α, VGEF, and G-CSF in MDMD versus controls, while there were no significant differences between controls and simple depression. Moreover, the interaction terms showed that the stimulated production of sIL-1RA, FGF, PDGF, and IFN-γ was significantly greater in MDMD than in simple depression.
Differences in lipopolysaccharide + phytohaemagglutinin-stimulated changes in cytokines/growth factors in healthy controls (HC) and depressed patients divided into those with simple depression and Major Dysmood Disorder (MDMD)
Associations with the features of MDMD
To examine the associations between the immune profiles and the key features of MDMD we performed GEE analyses and analyzed the interactions between time x HDRS, time x ROI, and time x PCs lifetime/recent SB. Table 4 shows significant time x ROI interactions for all nine immune profiles and these effects remained significant after FDR p-correction (at p<0.043). GEE analyses performed on the solitary cytokines/growth factors showed significant time x ROI interactions for sIL-1RA (W=8.18, p=0.004), IL-5 (W=7.08, p=0.008), CXCL8 (W=6.22, p=0.013), IL-9 (w=5.95, p=0.017), IL-15 (W=9.86, p=0.002), IL-17 (W=4.03, p=0.045), FGF (W=6.06, p=0.014), G-CSF (W=7.81, p=0.005), GM-CSF (W=5.75, p=0.017), IFN-γ (W=5.60, p=0.018), PDGF (W=6.58, p=0.010), CCL5 (W=5.22, p=0.022), TNF-α (W=5.57, p=0.018), and VGEF (W=5.52, p=0.019).
Associations between immune/growth factor profiles and the features of Major Dysmood Disorder, namely reoccurrence of illness (ROI), Hamilton Depression rating Scale (HDRS) score, and lifetime and recent suicidal behaviors (SB)
Table 4 shows significant time x HDRS interactions for all immune profiles (except Th17 and CIRS) and these differences remained significant after FDR p-correction (at p<0.036). GEE analyses performed on the solitary cytokines/growth factor showed significant time x HDRS interaction effects for sIL-1RA (W=5.40, p=0.020), IL-5 (W=7.32, p=0.007), CXCL8 (W=6.78, p=0.009), IL-9 (w=7.32, p=0.007), IL-12 (W=4.37, p=0.037), IL-15 (W=17.73, p<0.001), IL-17 (W=7.18, p=0.007), G-CSF (W=8.05, p=0.005), IFN-γ (W=5.71, p=0.017), PDGF (W=7.39, p=0.007), CCL5 (W=5.61, p=0.018), TNF-α (W=4.08, p=0.043), and VGEF (W=8.17, p=0.004).
Table 4 shows that all time x PC lifetime SB interactions were significant, except Th17 and Th2, and that the effects remained significant after FDR p-correction (at p<0.042). GEE analyses performed on the cytokines/growth factors showed significant time x PC lifetime SB interactions for sIL-1RA (W=12.34, p<0.001), IL-15 (W=6.75, p=0.009), FGF (W=4.55, p=0.039), IFN-γ (w=4.09, p=0.043), CCL5 (W=4.50, p=0.039), and VGEF (W=7.15, p=0.006). In addition, there were significant interactions between time x PC recent SB and M1, Th1, Th2, IRS, Tcell growth, and NT profiles. GEE analyses performed on the cytokines/growth factors showed significant time x recent SB interactions for sIL-1RA (W=8.0, p=0.005), IL-2 (W=4.30, p=0.038), GM-CSF (W=5.44, p=0.020), CXCL10 (w=4.29, p=0.038), CCL3 (W=4.34), p=0.037), and VGEF (W=3.97, p=0.046).
Construction of a ROI-immune score and prediction of the phenome
To delineate the prediction of the phenome using the immune profiles we performed multiple regression analysis with the phenome score as dependent variable and the immune profiles as explanatory variables while allowing for the effects of age, sex, BMI, education, and smoking. The phenome score was associated with the residualized M1 (partial correlation coefficient=0.438, p=0.002), Th1 (r=0.432, p=0.002), Th2 (r=0.404, p=0.005), Th17 (r=0.382, p=0.008), IRS (r=0.473, p<0.001), CIRS (r=0.302, p=0.039), neurotoxicity (r=0.459, p=0.001), Tcell growth (r=0.477, p=0.001), and growth factor (r=0.450, p=0.001) values after partialling out the effects of the unstimulated levels. These effects remained significant after FDR p-correction. We found that 82.3% of the variance in the phenome score was predicted (F=23.88, df=8/41, p<0.001) by ROI (β=0.747, t=10.09, p<0.001), neurotoxicity (β=0.434, t=3.67, p<0.001), CIRS (β=-0.275, t=-2.31, p=026), and age (β=-0.226, t=-3.16, p=0.003), while education (p=0.234), sex (p=0.095), BMI (p=0.393) and TUD (p=0.540) were all not significant. Figures 1–3 show the partial regressions of the phenome on staging, residualized neurotoxicity, and CIRS values, respectively. Moreover, substituting the neurotoxicity and CIRS in this regression with the IRS/CIRS ratio (computed as z IRS – z CIRS) showed that the IRS/CIRS ratio was significantly associated with the phenome score (β=0.418, t=2.98, p=0.005). Figure 4 shows the partial regression plot of the phenome score on the z IRS/CIRS ratio.
Partial regression of the phenome score on the recurrence of illness (ROI)
Partial regression of the phenome score on the immune-neurotoxicity index
Partial regression of the phenome score on an index of the compensatory immune-regulatory system
Partial regression plot of the phenome score on the immune-inflammatory response system (IRS) / compensatory immune-regulatory system (CIRS) ratio
We found that one general factor could be extracted from the residualized neurotoxicity (loading: 0.645), T cell growth (0.719), growth factor (0.702) scores, number of depressive episodes (0.884), number of all episodes (0.976), and PC lifetime SB (0.777). This factor explained 57.6% of the variance and showed adequate Cronbach alpha (0.894), composite reliability (0.904), and rho_A (0.932) values. This ROI-IMMUNE construct was strongly associated (t=18.76, p<0.001) with the phenome score with a zero-order correlation coefficient of 0.825 and a partial correlation coefficient of 0.799, i.e., after adjusting for the effects of age (t=-2.70, p=0.010), sex, BMI, education (all non-significant). There were no significant effects of the five drug state variables on the phenome scores.
Results of network, annotation, and enrichment analysis
The protein network of the 14 DEPs that are elevated in MDMD is shown in Figure 5A. It consists of 14 nodes with 79 edges, exceeding the predicted number (n=11) with a p-enrichment value of <1.0E-16. This network has an average node degree of 11.3 and an average local clustering coefficient of 0.967. Additionally, we built a first-order PPI network with 50 interactions in the first shell and none in the second and this network has 64 nodes and 760 edges, exceeding the predicted amount (n=177; p<1.0E-16). This network has an average node degree of 23.8, a mean local clustering coefficient of 0.707, a network diameter of 3 and a radius of 2, a typical path length of 1.651, a network density of 0.377, and a heterogeneity of 0.511. In decreasing order of significance, the top five seed hubs were TNF (degree=50), VEGFA (45), CXCL8 (42), IFNG (39), and CSF3 (37). FGF2 (betweenness centrality=0.0275) and PDGFA were the top two non-hub bottlenecks (0.0147). As a result, the first-order PPI network’s backbone is composed of four immune DEPs and three growth factor DEPs. MCL cluster analysis (inflation parameter: 3) indicated two communalities, one focusing on immunological DEPs and the other on growth factors (see Figure 5B for a first-order network).
A. The protein-protein interaction network (PPIN) of the 14 differentially expressed immune/growth factors of Major DysMood Disorder (MDMD); 5B: the first order PPIN of the growth factor cluster of MDMD
Table 5 summarizes the results of an enrichment study performed on all DEPs in MDMD’s first-order PPI network and the KEGG pathway classifications determined using STRING. The most over-represented KEGG pathways were immune-inflammatory pathways and viral infections, as well as JAK-STAT and MAPK pathways. Additionally, this table summarizes the KEGG pathways that were significantly enriched in the second growth factor cluster, including the RAS, MAPK, Rap1, and PI3K-Akt signaling pathways. Additionally, the same table includes the diseases that were enriched in all immunological DEPs associated with MDMD. Figure 6 shows the top ten Elsevier Pathways that are overrepresented in the first order protein-protein interaction network of MDMD.
Bar chart with the top ten Elsevier Pathways that are overrepresented in the first order protein-protein interaction network of Major DysMood Disorder
KEGG pathway classifications of the differently expressed proteins (DEPs) in Major Dysmood Disorder
The results of MCODE analysis for all first order DEPs utilizing GO biological and molecular terms are shown in Table 6. We identified three molecular complexes that represent cytokine responses, nuclear factor (NF)-B signaling regulation, and necrotic cell death. MCODE analysis of growth factor DEPs revealed the presence of a single cluster of transmembrane receptor protein kinase activity. Table 7 summarizes the top ten REACTOME pathways that were enriched in the first-order PPI network (as evaluated using OmicsNet), mostly involving Toll-Like Receptor signaling. Additionally, the same table illustrates the PANTHER biological processes that are enriched in the PPI network, with a particular emphasis on transcription through RNA polymerase, viral activities, and angiogenesis.
Results of Molecular Complex Detection (MCODE) analysis performed on differentially expressed proteins (DEPs) of major DysMood Disorder
REACTOME pathways and PANTHER biological processes statistically over-represented in the differentially expressed proteins of Major Dysmood Disorder
The GOnet enrichment analysis (biological process, q value threshold of <0.0001; p value threshold of <2.5E-9) is shown in Figure 7, highlighting the 14 seed MDMD DEPs and the significant GO annotations. The GO terms that are substantially overrepresented in the DEP list include inflammation, proliferation, STAT protein phosphorylation, JAK-STAT pathway regulation, and cellular responses to stress, lipid and LPS.
Results of GOnet enrichment analysis showing the 14 seed differentially expressed proteins of Major DysMood Disorder and their significant GO annotations
Discussion
Activated immune profiles in MDMD
The first key conclusion of this research is that an acute episode of MDMD is characterized by enhanced activation of all immune profiles and 14 of the 23 cytokine/growth factors in stimulated (but not unstimulated) diluted whole blood as compared with controls. Notably, immunological profiles and single cytokines/growth factors did not differ between simple depression and controls but were consistently greater in MDMD than in simple depression, with M1, Th1, T cell growth, and neurotoxicity profiles, sIL-1RA, FGF, and IFN-γ being significantly higher. These results support the IRS/CIRS hypothesis of depression, which states that in depression, IRS profiles such as M1, Th1, Th17, as well as CIRS profiles, particularly Treg and Th2, are considerably active [1]. We reviewed the many previous data on the serum cytokines when MDD/MDE was compared to controls [1,3,4]. Nonetheless, most of these findings are not entirely comparable to ours since we focused on MDMD, a subgroup of MDD/MDE. One significant contrast with past findings concerns the Th1 profile, as a meta-analysis [2] revealed decreased IFN-γ production in depression, but the present and earlier investigations found enhanced stimulated IFN-γ production [31]. Nonetheless, serum IFN-γ levels are often undetectable, and hence it is contradictory to assume that depressed individuals have lower blood levels of this cytokine than the barely detectable control values. Additionally, increased IFN-γ production was identified as a characteristic of MDMD in the present investigation, highlighting the critical involvement of the Th1 phenotype and cell-mediated immunity in that illness [32].
Notably, the current research revealed that MDMD was associated with a significantly elevated neurotoxicity profile, supporting the neurotoxicity theory of depression [1]. MDMD exhibits enhanced neurotoxicity because of increased production of IL-15, IL-17, TNF-α, IFN-γ, CXCL8, CXCL10, and CCL5, all of which have neurotoxic properties [1,33]. Notably, the adverse effects of those cytokines may exacerbate the neurotoxic consequences of increased RONS, lipid and protein oxidation, aldehyde formation, hypernitrosylation, and decreased antioxidant and neurotrophic defenses in MDMD [18,19]. Additionally, this research showed that MDMD is associated with elevated levels of the growth factors VEGF, PGDF, and FGF. Previously, a meta-analysis revealed elevated FDF levels in depression, however, results on VEGF and PDGF levels were more controversial [5,34–38]. Nonetheless, the present study’s findings are more appropriate since stimulated cultures of diluted whole blood more accurately mirror the in vivo scenario, and MDMD is a more accurate model than MDD/MDE. Interestingly, these three growth factors regulate cell division, endothelial cell chemotaxis, mitogen-activated protein kinase (MAPK) signaling pathways and angiogenesis (3 items (human) - STRING interaction network (string-db.org) and as such these growth factors may contribute to the increased immune responses.
Our results confirm that the clinical study population for major depression is divided into two clinically and biologically different subtypes, i.e., MDMD and simple depression. We generated the diagnosis MDMD in this investigation using previously described clinical characteristics [18,19]. As such, immune profile tests were used as a proof of concept of the MDMD construct. The primary goal of precision psychiatry should be to define the correct (reliable and cross-validated) model of depression as a major psychosis or a severe medical condition and to distinguish this class from depressive-like emotional reactions [20]. Thus, MDMD is clinically defined by an increase in ROI and phenome scores, as well as an increase in melancholia and psychotic symptoms, suicidal behaviors, and disabilities, lowered quality of life, as well as activated IRS/CIRS/neurotoxin/growth factor and redox pathways [this study, 19,18,20]. As a result, future research on depression should always take MDMD into account rather than MDD or MDE [20,39].
Immune profiles and the features of MDMD
The second main conclusion of this research is that activated immune profiles are strongly linked with MDMD characteristics, including increased suicidal behaviors, staging, and severity of the phenome. To begin, both lifetime and current suicide behaviors are related with increases in most immunological profiles, most notably the IRS and neurotoxicity profiles. Again, in light of the IRS/neurotoxicity hypotheses of suicidal ideation and attempts [9], the neurotoxic profile and neurotoxic cytokines (e.g., IL-2, IL-15, CCL5, CXCL10) may exacerbate the adverse effects of RONS/OSTOX on suicidal behaviors [9].
Second, by integrating the number of depressive and total episodes and lifetime suicidal behaviors into a ROI or staging index, we were able to corroborate the ROI model [10,18,20]. Furthermore, our findings reveal that ROI and activated immune-redox pathways are not only intrinsically linked but also generate a novel ROI-immune pathway phenotype, indicating that ROI and IRS, neurotoxicity and T cell growth profiles are manifestations of a common core in MDMD. These results corroborate our prior findings that ROI and immune-redox indicators belong to a common core [18,20]. As discussed in the introduction, there is some evidence that the frequency of depressive and manic episodes is related to serum IRS (sIL-1RA, IL-6) and CIRS (IL-10) biomarkers and to the stimulated expression of activation markers [16]. It is worth comparing this to “epileptogenic basolateral amygdala kindling,” a model for the recurrence of epileptic seizures in which a seizure raises the probability of subsequent seizures [40]. Notably, a shared core underpins the frequency and uncontrollability of seizures, comorbid depression and anxiety, and redox pathways, showing that this pathway phenotype identifies a particularly severe form of TLE [20]. Similarly, the kindling theory of affective disorders states that as episodes recur, the latter become increasingly sensitized [41], and our findings indicate that this “affective kindling” is a central feature of MDMD and is caused by aberrations in immune-redox pathways [18,20].
Thirdly, the MDMD phenome might be understood as a common core comprised of the intensity of depression and anxiety, the presence of melancholia or psychotic features, and recent suicidal behaviors which incorporate various aspects of suicidal ideation and attempts. Most notably, a considerable proportion (around 80%) of the variance in the phenome score may be attributed to ROI and the IRS or neurotoxicity profiles (positively related) and CIRS and age (inversely associated). These results corroborate the IRS/CIRS hypothesis of affective disorders [1], demonstrating that the neurotoxicity profile may result in affective disorders, particularly in people with weakened or inadequate CIRS defenses. Additionally, the results corroborate our prior findings that the ROI-redox pathway phenotype substantially predicts the phenome of mood disorders [18–20].
Increased immune responsivity in MDMD
The third important conclusion of this research is that MDMD and its features relate to the stimulated immune responses and not the unstimulated activity levels of the immune profiles. As such, the stimulated immune/growth factor production is part of the ROI and may, therefore, be ascribed to sensitization (kindling) processes [10–13,16]. Since ROI had no effect on unstimulated profiles, we may infer that not only ROI-induced sensitization but also continuous immunostimulatory triggers are required to account for the enhanced immune and growth factor responses. By inference, because our ex vivo stimulated whole blood assay reflects the measurement of serum cytokines in vivo, we can deduce that the elevated serum cytokine levels frequently observed in depression [1,2] are a result of both sensitization and continuous presence of immunostimulatory triggers during the acute phase of illness. Recently, a variety of trigger factors have been identified that may account for these effects, including increased LPS load from increased bacterial translocation and leaky gut [42], LPS particles on the membrane of outer membrane vesicles that circulate in the serum [43], increased LPS contents due to apical periodontitis [44], continuous stimulation of the Toll-Like Receptor (TLR)-Redox cycle by damage-associated molecular patterns including oxidatively modified neoepitopes which are abundant in affective disorders [45], other pathogen-associated molecular patterns including due to latent viral infection such as with cytomegalovirus [16], and immunostimulatory effects of psychological stressors and medical disorders that activate the cytokine network [46]. The impact of the latter is further corroborated by our disease enrichment analysis which shows that a multitude of (auto)immune diseases are enriched in the PPI network built using the DEPs of MDMD (see listings in Figure 6 and Table 5).
Notably, we observed no impact of antidepressant intake on the ex vivo (un)stimulated immunological profiles, even though these medications have a strong negative immunoregulatory effect in vitro, attenuating cell-mediated immune responses by lowering IFN-γ production while raising IL-10 production [21,22]. By extrapolation, despite antidepressants’ inherent immunoregulatory properties, IRS and neurotoxicity profiles in serum and stimulated whole blood cultures remain elevated, because of the continual stimulation throughout the acute phase of MDMD. By inference, these trigger factors as well as the molecular networks that underpin the sensitized immune-growth responses are new drug targets to treat affective disorders.
Using annotation and enrichment analysis we were able to pinpoint the most important molecular pathways and functions that play a role in MDMD, namely crosstalks between transmembrane receptor protein kinase activated STAT protein phosphorylation, hyperresponsive TLR/NF-κB and JAK-STAT pathways [46], and hyperresponsivity of the main proliferation/survival pathway, namely PI3K/Akt/RAS/MAPK signaling [47].
Conclusions
Increased M1, Th1, Th2, Th17, Treg, IRS, CIRS, neurotoxicity, and growth factor profiles are seen in MDMD. ROI, suicidal behaviors, and the phenome of depression are all linked to immunotoxicity and growth factor profiles. Factor analysis reveals that ROI characteristics and immune-growth factor profiles may be used to construct a novel pathway phenotype. ROI and immune profiles account for around 80.0 percent of the variance in the phenome. The exaggerated immune responses in MDMD may be explained by ROI-induced sensitization coupled with continuous immunostimulatory triggers. This study delineated new drug targets to treat affective disorders, namely ROI-associated sensitization of immune and growth factors and the underlying molecular pathways, namely the transmembrane receptor protein kinase triggered STAT protein phosphorylation, TLR/NF-B, JAK-STAT, and the proliferation/survival PI3K/Akt/RAS/MAPK pathway.
Author’s contributions
All the contributing authors have participated in the manuscript. M.M., M.R. and P.S. designed the study and M.M. performed the statistical, network, enrichment, and annotation analyses. M.R. and K.J. recruited patients and controls. Assays were performed by P.S. and S.K. A.S. secured financial support. M.M., M.R., K.J., P.S., S.K., P.G., and A.S. contributed to the interpretation of the data and writing of the manuscript and agreed to publish the final version of the manuscript.
Funding Statement
AMERI-ASIA MED CO, Ltd, supported this work
Institutional Review Board Statements
The study was conducted according to International and Thai ethics and privacy laws. Approval for the study was obtained from the Institutional Review Board of the Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand (#528/63).
Informed Consent Statement
All controls and patients provided written informed consent prior to participation in this study.
Data Availability Statement
The dataset generated during and/or analyzed during the current study will be available from the corresponding author (M.M.) upon reasonable request and once the dataset has been fully exploited by the authors.
Conflict of interest
The authors have no conflict of interest with any commercial or other association in connection with the submitted article.
Data Availability
The dataset generated during and/or analyzed during the current study will be available from the corresponding author (M.M.) upon reasonable request and once the dataset has been fully exploited by the authors.
Footnotes
e-mail addresses: dr.michaelmaes{at}hotmail.com, muanpetch.mp{at}gmail.com, ket.kett{at}hotmail.com, yokpim{at}gmail.com, siriwanklinchanhom{at}gmail.com, piotr.galecki{at}umed.lodz.pl, atapol.s{at}gmail.com
References
