Abstract
This study aimed at identifying molecular biomarkers of inflammation-related depression in order to improve diagnosis and treatment. For this, we performed whole-genome expression profiling from peripheral blood in a naturalistic model of inflammation-associated major depressive disorder (MDD) represented by comorbid depression in obese patients. We took advantage of the marked reduction of depressive symptoms and inflammation following bariatric surgery to test the robustness of the identified biomarkers. Depression was assessed during a clinical interview using Mini-International Neuropsychiatric Interview and the 10-item, clinician-administered, Montgomery–Asberg Depression Rating Scale. From a cohort of 100 massively obese patients, we selected 33 of them for transcriptomic analysis. Twenty-four of them were again analyzed 4–12 months after bariatric surgery. We conducted differential gene expression analyses before and after surgery in unmedicated MDD and non-depressed obese subjects. We found that TP53 (Tumor Protein 53), GR (Glucocorticoid Receptor), and NFκB (Nuclear Factor kappa B) pathways were the most discriminating pathways associated with inflammation-related MDD. These signaling pathways were processed in composite z-scores of gene expression that were used as biomarkers in regression analyses. Results showed that these transcriptomic biomarkers highly predicted depressive symptom intensity at baseline and their remission after bariatric surgery. While inflammation was present in all patients, GR signaling over-activation was found only in depressed ones where it may further increase inflammatory and apoptosis pathways. In conclusion, using an original model of inflammation-related depression and its remission without antidepressants, we provide molecular predictors of inflammation-related MDD and new insights in the molecular pathways involved.
Similar content being viewed by others
Introduction
Over the past two decades, compelling evidence has emerged to highlight the role of inflammation in the onset and perpetuation of major depressive disorders (MDD). Epidemiological and clinical studies have repeatedly shown the high prevalence of depressive symptoms in chronic inflammatory disorders. Furthermore, chronic treatment with the pro-inflammatory cytokine, interferon-alpha, in medically ill patients was found to be responsible for the development of major depression in over 30% of patients1,2. Consistent with these data, a recent meta-analysis confirmed that MDD is associated with elevated levels of pro-inflammatory cytokines, such as interleukin (IL)-6 and tumor necrosis factor (TNF)-α3. Additionally, genetic studies have revealed the role of immune genetic variants in MDD4. Interestingly, inflammation has been recently hypothesized to also have a role in antidepressant treatment resistance, as patients presenting signs of systemic low-grade inflammation display a poor response to classical antidepressants, such as selective serotonin reuptake inhibitors5. Thus, a better understanding of the molecular mechanisms underlying inflammation-related depression is needed. In particular, the identification of biomarkers specific to this depression subtype is required for improving its diagnosis and treatment.
Transcriptomic analysis from peripheral blood has successfully been used to reveal biological pathways linked to MDD6,7. A recent study examined the implication of specific molecular pathways in inflammation-related depression by measuring self-reported depressive symptoms 2 h after endotoxin injection in healthy subjects8. However, hypothesis-free transcriptomic analyses have not been performed specifically in diagnosed inflammation-associated MDD patients. To address this issue, we used a naturalistic model of inflammation-related depression represented by comorbid depression in obese patients. Indeed, obesity is a chronic medical condition characterized by a low-grade inflammatory state associated with an increased prevalence of depression9. Remarkably, systemic inflammation, reflected by increased levels of C-reactive protein levels, predicts obesity-related depressive symptoms better than metabolic health per se10. Interestingly, bariatric surgery-induced weight loss was found to correlate with a reduction of both inflammation and depressive symptoms11,12.
In this study, we sought to identify transcriptional control pathways associated with inflammation-related depression in unmedicated obese patients. For this purpose, gene expression profiling of peripheral blood was performed and compared between depressed and non-depressed severely obese patients before and after bariatric surgery, with the hypothesis that these biological pathways would be changed after surgery-induced weight loss and subsequent decrease in inflammation. Upstream regulator analyses were employed to uncover the molecular signaling pathways that differentiate depressed from non-depressed subgroups of patients. Finally, the most discriminating signaling pathways identified were processed in composite z-scores of gene expression in order to be used as biomarkers predicting depressive symptoms.
Methods
Patients
One hundred severely or morbidly obese patients [body mass index (BMI) > 35 kg/m2] awaiting bariatric surgery were recruited from the services of digestive and bariatric surgery at two private clinics (Tivoli and Jean Villar clinics, Bordeaux, France). Patients were scheduled to receive either a sleeve gastrectomy or a gastric bypass. Patients were chosen on the basis of their levels of depressive symptoms in order to cover a range from non-depressed to clinically depressed patients. Seventy-three percent of patients (n = 24) were followed up at 4–12 months after surgery (mean time at follow-up 6.5 months). The variability in the follow-up interval times is due to the patient not coming at the first appointments. Exclusion criteria were: age >65 years old; acute or chronic inflammatory conditions (other than obesity); current treatment with antidepressants or any other psychotropic drug; current diagnosis of psychiatric disease (except for major depression); and/or severe medical illness including neurological disorders and cancer. Patients using anti-inflammatory medications (checked at recruitment and at medical consultation) were excluded. The study was approved by the local Committee for the Protection of Persons (Bordeaux, France). All patients provided written informed consent after reading a complete description of the study.
Clinical evaluation
DSM criteria for current major depression were determined using the Mini-International Neuropsychiatric Interview (MINI) administered during a semi-structured interview by trained raters at baseline and after bariatric surgery13. Concomitantly, the intensity of depressive symptoms was assessed using the 10-item, clinician-administered, Montgomery–Asberg Depression Rating Scale (MADRS)14.
Genome-wide transcriptional profiling
Blood samples were collected twice, before (n = 33) and after surgery (n = 24) in study participants. Five ml of venous blood was collected per patient on vacuum tubes (PAXgene Blood RNA system; PreAnalytiX GmbH, Hombrechtikon, Switzerland). The samples were maintained at room temperature for >2 h as required for stabilization of RNA and then kept at −20 °C until RNA extraction. Total RNA was extracted using the PAXgene Blood RNA Kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s protocol. The quality of the total RNA was assessed using RNA Nano chips on a Bioanalyser 2100 (Agilent, Boeblingen, Germany). All samples had an RNA Integrity Number (RIN) score >8.0. RNA concentrations were measured with a Nanodrop spectrophotometer.
Gene expression profiles of the 57 samples (33 before and 24 after surgery) were analyzed in a single batch at the GeT‐TRiX facility (GenoToul, Génopole Toulouse, France) using Agilent SurePrint G3 Human GE V2 8x60K (Design 039494) microarrays following the manufacturer’s instructions. Microarray data and experimental details are available in NCBI’s Gene Expression Omnibus15 and are accessible through GEO Series accession number GSE99725.
sPLS regression
Sparse Partial Least Square (sPLS) regression16 using the R package, mixOmics (http://www.mixOmics.org) was performed to select the expressed genes that correlated the best to MDD diagnosis or MADRS items. From the 2000 probes found highly correlated with a correlation coefficient r > 0.5, 1474 unique genes were identified.
Differential gene expression analyses
Raw data (median of pixel intensity) were analyzed using R (www.r-project.org, R v. 3.1.2) Bioconductor packages (www.bioconductor.org, v 3.017) as described in GEO entry GSE99725. Briefly, raw data were filtered, log2 transformed, corrected for bath effects (washing and labeling serials), and normalized using the quantile method18. A first exploratory and statistical analysis showed a possible correlational structure among gene expression, which could negatively impact the multiple testing procedures. We applied the FAMT method19 to reduce the dependence structure. A model was fitted using the limma lmFit function20 considering “patients” as a blocking factor for the time pairs using a duplicate correlation function. A correction for multiple testing was then applied using the Benjamini–Hochberg procedure21 for false discovery rate (FDR). Probes with FDR ≤ 0.05 were considered differentially expressed between conditions.
Impact of blood cell proportions in gene expression profiles
We assessed whether the observed gene expression changes in the different analyses were related to changes in cell proportions in the blood samples using the Cell-type Computational Differential Estimation CellCODE R package22.
Transcript origin analysis
To identify specific cell types that contributed to the empirically observed differences in gene expression, we applied transcript origin analysis23 to separate lists of upregulated and downregulated genes. TOA tested for significant overrepresentation of genes preponderantly expressed by a specific subset of blood cells (i.e., monocytes, dendritic cells, CD4-positive T cells, CD8-positive T cells, B cells, or NK cells) using cell-specific reference transcriptomes derived by flow cytometric isolation23,24.
Transcriptional pathways identification, comparison analysis, and composite z-scores analyses
The biological pathways involved in MDD were predicted from lists of differentially expressed genes using the “Upstream Regulator Analysis” from Ingenuity Pathway Analysis (IPA, Qiagen, USA).
To calculate the composite z-scores representing a given signaling pathway, we first extracted the list of target genes for a given transcription factor from ChIP-seq experiments available in ChEA 2016, an integrative database of public ChIP-seq data (available in http://amp.pharm.mssm.edu/Enrichr/). Then, from the list of target genes for a given transcription factor, we selected the ones that were expressed in our data set. For each of the patients, the expression data for each gene was first standardized by calculating a z-score, z-score = \(\frac{{x - \mu }}{\sigma }\), where x is the gene expression data, μ the mean of the gene expression for the population (all patients of the study) and σ the standard deviation for the population. For each patient, the composite z-score for a given transcription factor was then calculated as the mean of the target genes’ z-scores.
Cytokines measures in serum
High-sensitivity C-reactive protein (CRP) was measured from fasting blood samples as described in ref. 10 using an ELISA kit (Millipore, Billerica, MA). Similarly, high-sensitivity Interleukin-6 (IL-6) was measured from an ELISA kit (R&D Systems, Minneapolis, MN) according to manufacturer instructions.
Statistics
Two-tailed t-tests, ANOVA, and Multivariate regression analyses were done with the software Statistica and the graphs designed with GraphPadPrism. A p-value of p < 0.05 was considered significant.
Results
Characteristics of participants
The general characteristics of the participants are described in Table 1. Since study participants were chosen at baseline to cover a range from non-depressed to clinically depressed subjects, the percentage of patients with MDD in the present study (45.45%) is higher than what is generally found in the obese population. There were no statistical differences in terms of age, gender, BMI (before or after surgery), CRP, IL-6 or interval time after surgery between non-depressed versus depressed patients (unpaired t-tests except Fisher’s exact test for % women). As expected, MADRS scores were significantly higher in depressed patients at baseline (p < 0.0001) and remained significantly higher after surgery (p = 0.0068) although all patients were no more depressed. The number of patients that were lost after follow-up was balanced between the 2 groups, 4 depressed and 5 non-depressed patients. These patients did not come to the follow-up medical consultation probably by lack of interest but not because of health complications. BMI was strongly reduced after surgery in the 24 patients with available data at follow-up either as a whole or within depressed or non-depressed subgroups (paired t-test p < 0.0001 for each analysis). Similarly, the percentage of patients with MDD was significantly decreased after surgery compared to baseline (from 45.45 to 0%, Fisher exact test p < 0.0006) as well as the MADRS scores in the depressed subgroup (paired t-test, p < 0.0002).
Gene expression analyses
As a first exploration of the peripheral blood gene expression data from the 33 obese patients at baseline, we performed an sPLS regression analysis with the variables related to depressive symptoms. The heatmap (Fig. 1A) illustrates the strength of correlation between the 2000 most correlated probes and the variables related to depressive symptomatology. Strong correlations with gene expression were found with total MADRS score and MDD diagnosis. Not surprisingly, the correlation was the weakest with the MADRS item measuring reduced appetite. When the selected variable was the total MADRS score, the 2000 most correlated probes corresponded to 1474 unique genes, 428 downregulated and 1046 upregulated. When this gene data set was submitted to IPA, the top canonical pathway identified was the Glucocorticoid Receptor (GR) signaling pathway (p-value 2.13 × 10−5 with an overlap of 12.5%), while the top upstream regulators were cystatin D (CST5), tumor protein 53 (TP53) and hepatocyte nuclear factor4 (HNF4A) (p = 1.09 × 10−8, p = 9.8 × 10−8 and p = 6.10 × 10−7, respectively). The IPA analysis could not predict whether these transcription factors were activated or inhibited in patients.
Then, we performed 4 differential gene expression analyses where we compared the gene expression of MDD to non-depressed patients at baseline, MDD patients after versus before bariatric surgery, non-depressed patients after versus before surgery, and lastly, MDD versus non-depressed (as they were at baseline) after surgery.
At a false discovery rate of 5%, the three first analyses provided a number of differentially expressed genes (DEG) as illustrated in Fig. 1B, with a Venn diagram proportional to the number of DEG. For the later analysis (MDD versus non-MDD after surgery), no significant DEGs were detected as expected since patients were no longer depressed in either group. We found no significant changes in immune cell subtypes’ proportion in the different analyses with significant DEG using CellCODE (Supplementary Table 1). For the MDD versus non-MDD comparison at baseline, the transcript origin analysis revealed that upregulated transcripts derived preferentially from NK cells, CD4+ T cells, and CD8+ T cells. For the MDD after versus before surgery comparison, the upregulated transcript derived from monocytes and the downregulated genes from monocytes, dendritic cells, and NK cells. For the non-MDD after versus before surgery comparison, the upregulated transcripts originated preferentially from CD4+ T cells and CD8+ T cells and the downregulated genes from NK cells. These results are presented in Supplementary Fig. 1. The list of DEG for each analysis, and the corresponding fold change, are presented in Supplementary Table 2.
Upstream regulator analyses
To identify the transcriptional pathways associated with MDD in these patients, we perform upstream regulator analyses from each list of DEG using IPA. The lists of upstream regulators for each analysis are provided in Supplementary Table 3. At a cutoff p-value of overlap p < 0.001, only one upstream regulator, TP53, was significant for the MDD versus non-MDD at baseline analysis, with 17 downstream target genes. For the MDD after versus before surgery analysis, 67 upstream regulators were detected, the three most significant were, respectively, dehydrotestosterone (p = 4.25 × 10−7, an androgen receptor agonist), TP53 (p = 1.6 × 10−6), and dexamethasone, a specific glucocorticoid receptor agonist (p = 2.84 × 10−6). TP53 and dexamethasone showed the highest number of downstream target genes with 67 and 71 target genes, respectively. Interestingly, a mechanistic network was revealed that linked dexamethasone and TP53 along with other upstream regulators of this data set (Fig. 2A). The NFκB (Nuclear Factor kappa B) complex appears as an important hub between dexamethasone and TP53 in this network. For the non-MDD after versus before surgery analysis, 7 upstream regulators were found significant but all of them with a low number of target genes, GATA1 (GATA binding protein 1) having the highest number with 9 target genes. To identify the upstream regulators that were in common or different between MDD and non-depressed patients and explaining the changes after surgery, a comparison analysis was performed (Fig. 2B). This analysis highlighted that among the 10 regulators sorted by activation z-score, several cytokines pathways (IL-4, IL15, IL2, IL-1β) were downregulated in both MDD and non-depressed patients after surgery. Remarkably, dexamethasone was found downregulated in the MDD after versus before surgery analysis only.
Transcriptomic pathway biomarkers of inflammation-related depression
The upstream regulator analyses above, in accordance with the sPLS regression, pinpointed dexamethasone/GR, NFκB, and TP53 as the most important and convergent transcriptional regulators associated with MDD. To assess whether these pathways could serve as biomarkers predicting inflammation-related depression in our sample, we calculated for each of these regulators a composite z-score of gene expression using all their target genes, identified from ChIP-seq experiments and expressed on the microarrays (see “Methods” section for details). The lists of target genes representing NR3C1 (for dexamethasone/GR), RELA (for NFκB), and TP53 pathways and used to calculate the composite z-scores are given in Supplementary Table 4. Results from a two-way ANOVA indicated that MDD patients could be differentiated from non-depressed patients at baseline by the composite z-scores of TP53, NR3C1, and RELA, which were normalized after surgery when patients were no longer depressed (Fig. 3) (for each graph two-way ANOVA interaction (MDD x Surgery) p < 0.01, Bonferroni post hoc tests: MDD versus non-depressed patients before surgery p < 0.001, MDD versus non-depressed patients after surgery p = ns, MDD patients before surgery versus after surgery p < 0.001).
Furthermore, we tested by linear regression analysis whether these composite z-scores could predict the intensity of depressive symptoms at baseline using total MADRS scores (Fig. 4A). Indeed, highly significant linear relationships were detected for each of the three pathways. TP53 showed the strongest relationship (β = 0.80, F(1,31) = 56.8, p < 10−7), followed by NR3C1 (β = 0.79, F(1,31) = 50.3, p < 10−7) and RELA (β = 0.78.5, F(1,31) = 49.8, p < 10−7). These associations remained significant after adjusting for gender, age, and BMI (Supplementary Table 5). Using the same composite z-scores in regression analyses with BMI, we found no relationship (Supplementary Fig. 2) suggesting the specificity of the transcriptomic biomarkers for depressive symptoms. Since depressive symptoms were considerably reduced in most obese patients after surgery, we hypothesized that the difference in composite z-scores before and after surgery should predict the decrease in MADRS scores from baseline to post-surgery. Consistent with this hypothesis, multivariate linear regression analyses showed a significant linear relationship for the three transcription factor composite z-scores (β = 0.65, F(1,22) = 16.4, p < 0.001 for TP53, β = 0.47, F(1,22) = 6.34, p < 0.02 for NR3C1 and β = 0.49, F(1,22) = 7.19, p < 0.015 for RELA) (Fig. 4B).
Discussion
Low-grade-inflammation in obesity originates primarily from the adipose tissue, in which immune cells accumulate and secrete inflammatory factors9,25. The resulting systemic inflammation is thought to contribute to increased neuroinflammatory processes that subsequently drive behavioral alterations, including depressive symptoms11. In the present study, we used this model of inflammation-related depression to reveal the transcriptomic pathways involved, taking advantage of the marked reduction of depressive symptoms and inflammation following bariatric surgery. Our transcriptomic and bioinformatic analyses identify TP53, NR3C1, and NFκB as the most important upstream regulators predicting depression in obese patients. Indeed, NR3C1/Glucocorticoid Receptor and TP53 signaling pathways were revealed as linked to depressive symptoms by an sPLS analysis using the 33 patients at baseline, and the involvement of these two upregulators were confirmed in the differential gene expression analysis comparing post-surgery versus baseline in MDD patients. Regarding NFΚB, it was first identified as an important intermediary in a mechanistic network including NR3C1 and TP53.
TP53 is known as a tumor suppressor, pro-apoptotic protein, which levels are altered in response to cellular stress. Increased apoptotic stress was reported previously in animal models of stress-induced depression and in MDD patients26,27. For example, a transcriptomic analysis of postmortem prefrontal cortex tissues from patients with a history of MDD revealed altered expression of apoptosis factors together with increased cytokines expression28. The importance of the TP53 pathway itself was described in an animal model of trauma by transcriptional analysis from blood, amygdala, and hippocampus29, and in a genetic study where DNA polymorphisms within the TP53 gene were found associated with MDD in a Slovak population30. Interestingly, TP53 was also shown to be increased in the adipose tissue of obese mice and responsible for cytokine production and senescence-like changes thereby contributing to insulin resistance31. Our finding that the TP53 pathway is involved in inflammation-related depression in obese individuals is consistent with these data. It is tempting to postulate that TP53 induces the secretion of pro-inflammatory factors in adipose tissue, which in turn triggers neuroinflammation and subsequent behavioral alterations.
Since MDD patients were not different from non-depressed patients in terms of BMI, it is very likely that other factors besides adiposity explain increased TP53 and NFκB signaling in MDD patients. Our analysis put forward NR3C1/GR as one possible factor since it came out in both sPLS and differential gene expression analyses as an important upstream regulator associated with MDD. Furthermore, in the comparison analysis, dexamethasone appeared as a significant upstream regulator in obese MDD patients but not in obese non-depressed patients. Evidence from the literature strongly suggests a role for GR in the development of MDD and its related neurobiological disturbances. Indeed, MDD is often preceded by a history of stress and chronically elevated glucocorticoid levels, which contribute to reduced hippocampal volume and impaired hippocampal neurogenesis commonly found in MDD32. Recently, GR signaling was found as the most significant pathway associated with trauma-related individual differences in a preclinical model as revealed by gene expression analysis. It was also the most significant transcription factor convergent across blood and brain tissues. Furthermore, a translational study discovered that GR sensitivity had a crucial role in antidepressant treatment response in MDD patients and in a mouse model of depression33. As for inflammation-related depression, a recent study examined candidate transcriptomic pathways contributing to depressed mood in healthy subjects treated acutely with endotoxin. NFκB and CREB (cAMP Responsive Element Binding Protein) signaling were found increased in subjects who developed depression but GR signaling was decreased8. This study confirmed earlier reports showing for example that activation of inflammatory processes in MDD is associated with a reduced GRα/β expression ratio in monocytes34,35 or that dexamethasone-stimulated gene expression in peripheral blood is lower in depressed patients than in healthy subjects36. Recently, decreased GR mRNA together with increased FKBP5 mRNA levels were found in drug-free depressed and treatment-resistant patients compared to healthy controls or antidepressant responsive patients, suggesting GR resistance37. Although these data may seem contradictory to ours for GR, it could be explained by the anti and pro-inflammatory actions of glucocorticoids, which depend on the timing of exposure38. Indeed, exposure to glucocorticoids prior to an immune challenge primes the neuroinflammatory responses in rodents39 as well as in human hippocampal progenitor cells40. Along this line, we can hypothesize that obese patients with depressive symptoms may have been exposed to more stress and elevated glucocorticoid levels before developing obesity, which would result in increased inflammation during fat accumulation. In support of this hypothesis, we showed in a previous study that patients who developed MDD during chronic interferon-α treatment displayed, as soon as the first administration of the cytokine (at the time where they were not depressed yet), an exaggerated ACTH and cortisol responses, suggestive of a sensitized stress response system41.
Interaction between GR and TP53 has been reported in neuronal cell cultures where dexamethasone through GR enhances TP53 activity by increasing its nuclear translocation and transcriptional activity42. Induction of TP53 by dexamethasone was also found in cultured tenocytes and in tendon biopsies from patients treated with dexamethasone43 supporting the role of GR in TP53 induction in vivo.
Besides these three main pathways, we found that several cytokines pathways were also reduced after surgery in both MDD and non-MDD patients (e.g., IL-4, IL-6). Of note, the TNF pathway was found significantly reduced in MDD patients after surgery, which is interesting in view of the efficiency of the anti-TNF agent infliximab to reduce depressive symptoms in antidepressant treatment-resistant patients44. HNF4A was also an interesting candidate from this study, being associated with glucose and lipid metabolism as well as with inflammation-related depression and antidepressant treatment44. Although the HNF4A pathway was detected in the sPLS analysis, we found no change in this pathway from the differentially expressed gene analyses may be due to lack of statistical power.
A number of studies have utilized gene expression profiles to investigate the determinants of MDD. Although few associated genes are found in common across studies, all of them implicate immune pathways, such as NFκB45,46,47, IL-648, IFN-α48,49,50, oxidative stress, or apoptosis28,47,51 signaling pathways. These convergent findings confirm that transcriptomic pathways represent stronger biomarkers than individual genes, probably due to statistical issues, effects of cellular heterogeneity and temporal dynamics, as well as multifactorial gene regulation, as discussed before for stress genomics52, and heterogeneity of depression types. Indeed, the present transcript origin analyses implicated lymphoid lineage cells (particularly CD4+ T cells, CD8+ T cells, and NK cells) as primary cellular mediators of differential gene expression in MDD versus non-depressed controls under baseline conditions, as well as pre- to post-surgery transcriptome changes in non-depressed individuals. In MDD patients, by contrast, pre- to post-surgery changes were dominated by transcripts typically expressed in myeloid lineage monocytes and dendritic cells (as well as NK cells). These data suggest the potential for distinct immunoregulatory dynamics in MDD (e.g., due to differential activity of the glucocorticoid signaling system as noted above).
In most studies, the transcriptomic pathways are identified from the differentially expressed genes by bioinformatic prediction based on transcription binding motifs using various softwares (e.g., DAVID53, TELiS54) or literature curation (IPA). Here, we used a new approach relying on ChIP-seq experiments. By this method, we believe that we captured a closer reflection of the transcriptional pathways involved. Indeed, the composite z-score of gene expression representing the pathways of TP53, GR, and RELA was highly predictive of MADRS scores in our patient sample. Furthermore, the same composite z-scores were also predictive of the remission of depressive symptoms after bariatric surgery, supporting the validity of these biomarkers.
The following limitations of our study must be considered. The main one resides in the low number of patients, especially after bariatric surgery. Another limitation is that the findings may be restricted to obesity-related depression and not relevant to all inflammation-related depressions. Clearly, replication of these findings in larger cohorts and in non-obese patients with inflammation-associated depression is needed. The variability in the interval times after surgery between patients is another difficulty but as these data did not differ on average between depressed and non-depressed patients, it had no direct implication on the results presented here. The mechanism of GR over-activation versus GR resistance in inflammation-related depression remains an open question that requires further investigation. Finally, we focused on the main transcriptional pathways detected in the study while many other pathways were found significant and others may have been missed by lack of statistical power.
In conclusion, this study highlights the importance of TP53, NFκB, and GR signaling pathways in a clinical model of inflammation-related MDD and their potential as predictive biomarkers to identify individuals vulnerable to this type of depression. The present data provide insights of the molecular mechanisms involved in inflammation-associated depression that can be used to improve diagnosis and treatment.
References
Capuron, L. & Miller, A. H. Immune system to brain signaling: neuropsychopharmacological implications. Pharm. Ther. 130, 226–238 (2011).
Pape, K., Tamouza, R., Leboyer, M. & Zipp, F. Immunoneuropsychiatry—novel perspectives on brain disorders. Nat. Rev. Neurol. 15, 317–328 (2019).
Goldsmith, D. R., Rapaport, M. H. & Miller, B. J. A meta-analysis of blood cytokine network alterations in psychiatric patients: comparisons between schizophrenia, bipolar disorder and depression. Mol. Psychiatry 21, 1696–1709 (2016).
Barnes, J., Mondelli, V. & Pariante, C. M. Genetic contributions of inflammation to depression. Neuropsychopharmacology 42, 81–98 (2017).
Strawbridge, R. et al. Inflammation and clinical response to treatment in depression: a meta-analysis. Eur. Neuropsychopharmacol. 25, 1532–1543 (2015).
Hepgul, N., Cattaneo, A., Zunszain, P. A. & Pariante, C. M. Depression pathogenesis and treatment: what can we learn from blood mRNA expression? BMC Med. 11, 28 (2013).
Lin, E. & Tsai, S.-J. Genome-wide microarray analysis of gene expression profiling in major depression and antidepressant therapy. Prog. Neuropsychopharmacol. Biol. Psychiatry 64, 334–340 (2016).
Cho, J. H.-J., Irwin, M. R., Eisenberger, N. I., Lamkin, D. M. & Cole, S. W. Transcriptomic predictors of inflammation-induced depressed mood. Neuropsychopharmacology 44, 923–929 (2019).
Castanon, N., Lasselin, J. & Capuron, L. Neuropsychiatric comorbidity in obesity: role of inflammatory processes. Front. Endocrinol. 5, 74 (2014).
Delgado, I. et al. Depressive symptoms in obesity: relative contribution of low-grade inflammation and metabolic health. Psychoneuroendocrinology 91, 55–61 (2018).
Capuron, L., Lasselin, J. & Castanon, N. Role of adiposity-driven inflammation in depressive morbidity. Neuropsychopharmacology 42, 115–128 (2017).
Emery, C. F. et al. Gastric bypass surgery is associated with reduced inflammation and less depression: a preliminary investigation. Obes. Surg. 17, 759–763 (2007).
Sheehan, D. V. et al. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J. Clin. Psychiatry 59(Suppl 20), 22–33 (1998).
Montgomery, S. A. & Asberg, M. A new depression scale designed to be sensitive to change. Br. J. Psychiatry 134, 382–389 (1979).
Barrett, T. & Edgar, R. Gene expression omnibus: microarray data storage, submission, retrieval, and analysis. Methods Enzymol. 411, 352–369 (2006).
Lê Cao, K.-A., Rossouw, D., Robert-Granié, C., & Besse, P. A sparse PLS for variable selection when integrating omics data. Stat. Appl. Genet. Mol. Biol. 7, Article 35 (2008).
Gentleman, R. C. et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5, R80 (2004).
Bolstad, B. M., Irizarry, R. A., Astrand, M. & Speed, T. P. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19, 185–193 (2003).
Friguet, C., Kloareg, M. & Causeur, D. A factor model approach to multiple testing under dependence. J. Am. Stat. Assoc. 104, 1406–1415 (2009).
Smyth, G. K. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3, 1–25 (2004).
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful. approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300 (1995).
Chikina, M., Zaslavsky, E. & Sealfon, S. C. CellCODE: a robust latent variable approach to differential expression analysis for heterogeneous cell populations. Bioinformatics 31, 1584–1591 (2015).
Cole, S. W., Hawkley, L. C., Arevalo, J. M. G. & Cacioppo, J. T. Transcript origin analysis identifies antigen-presenting cells as primary targets of socially regulated gene expression in leukocytes. Proc. Natl Acad. Sci. USA 108, 3080–3085 (2011).
Black, D. S., Cole, S. W., Christodoulou, G. & Figueiredo, J. C. Genomic mechanisms of fatigue in survivors of colorectal cancer. Cancer 124, 2637–2644 (2018).
Wellen, K. E. & Hotamisligil, G. S. Inflammation, stress, and diabetes. J. Clin. Invest. 115, 1111–1119 (2005).
Liu, W. et al. The role of neural plasticity in depression: from hippocampus to prefrontal cortex. Neural Plast. 2017, 6871089 (2017).
Lucassen, P. J. et al. Stress, depression and hippocampal apoptosis. CNS Neurol. Disord. Drug Targets 5, 531–546 (2006).
Shelton, R. C. et al. Altered expression of genes involved in inflammation and apoptosis in frontal cortex in major depression. Mol. Psychiatry 16, 751–762 (2011).
Daskalakis, N. P., Cohen, H., Cai, G., Buxbaum, J. D. & Yehuda, R. Expression profiling associates blood and brain glucocorticoid receptor signaling with trauma-related individual differences in both sexes. Proc. Natl Acad. Sci. USA 111, 13529–13534 (2014).
Mahmood, S., Evinová, A., Škereňová, M., Ondrejka, I. & Lehotský, J. Association of EGF, IGFBP-3 and TP53 gene polymorphisms with major depressive disorder in Slovak population. Cent. Eur. J. Public Health 24, 223–230 (2016).
Minamino, T. et al. A crucial role for adipose tissue p53 in the regulation of insulin resistance. Nat. Med. 15, 1082–1087 (2009).
Anacker, C., Zunszain, P. A., Carvalho, L. A. & Pariante, C. M. The glucocorticoid receptor: pivot of depression and of antidepressant treatment? Psychoneuroendocrinology 36, 415–425 (2011).
Carrillo-Roa, T. et al. Common genes associated with antidepressant response in mouse and man identify key role of glucocorticoid receptor sensitivity. PLoS Biol. 15, e2002690 (2017).
Carvalho, L. A. et al. Inflammatory activation is associated with a reduced glucocorticoid receptor alpha/beta expression ratio in monocytes of inpatients with melancholic major depressive disorder. Transl. Psychiatry 4, e344–e344 (2014).
Grosse, L. et al. Clinical characteristics of inflammation-associated depression: Monocyte gene expression is age-related in major depressive disorder. Brain Behav. Immun. 44, 48–56 (2015).
Menke, A. et al. Dexamethasone stimulated gene expression in peripheral blood is a sensitive marker for glucocorticoid receptor resistance in depressed patients. Neuropsychopharmacology 37, 1455–1464 (2012).
Cattaneo, A. et al. Whole-blood expression of inflammasome- and glucocorticoid-related mRNAs correctly separates treatment-resistant depressed patients from drug-free and responsive patients in the BIODEP study. Transl. Psychiatry 10, 232 (2020).
Horowitz, M. A. & Zunszain, P. A. Neuroimmune and neuroendocrine abnormalities in depression: two sides of the same coin. Ann. N. Y Acad. Sci. 1351, 68–79 (2015).
Frank, M. G., Miguel, Z. D., Watkins, L. R. & Maier, S. F. Prior exposure to glucocorticoids sensitizes the neuroinflammatory and peripheral inflammatory responses to E. coli lipopolysaccharide. Brain Behav. Immun. 24, 19–30 (2010).
Horowitz, M. A. et al. Glucocorticoids prime the inflammatory response of human hippocampal cells through up-regulation of inflammatory pathways. Brain Behav. Immun. 87, 777–794 (2020).
Capuron, L. et al. Association of exaggerated HPA axis response to the initial injection of interferon-alpha with development of depression during interferon-alpha therapy. Am. J. Psychiatry 160, 1342–1345 (2003).
Crochemore, C., Michaelidis, T. M., Fischer, D., Loeffler, J.-P. & Almeida, O. F. X. Enhancement of p53 activity and inhibition of neural cell proliferation by glucocorticoid receptor activation. FASEB J. 16, 761–770 (2002).
Poulsen, R. C. et al. Glucocorticoids induce senescence in primary human tenocytes by inhibition of sirtuin 1 and activation of the p53/p21 pathway: in vivo and in vitro evidence. Ann. Rheum. Dis. 73, 1405–1413 (2014).
Mehta, D. et al. Transcriptional signatures related to glucose and lipid metabolism predict treatment response to the tumor necrosis factor antagonist infliximab in patients with treatment-resistant depression. Brain Behav. Immun. 31, 205–215 (2013).
Mellon, S. H. et al. Alterations in leukocyte transcriptional control pathway activity associated with major depressive disorder and antidepressant treatment. Transl. Psychiatry 6, e821 (2016).
Chiang, J. J. et al. Depressive symptoms and immune transcriptional profiles in late adolescents. Brain Behav. Immun. 80, 163–169 (2019).
Malki, K. et al. Identification of genes and gene pathways associated with major depressive disorder by integrative brain analysis of rat and human prefrontal cortex transcriptomes. Transl. Psychiatry 5, e519 (2015).
Jansen, R. et al. Gene expression in major depressive disorder. Mol. Psychiatry 21, 339–347 (2016).
Mostafavi, S. et al. Type I interferon signaling genes in recurrent major depression: increased expression detected by whole-blood RNA sequencing. Mol. Psychiatry 19, 1267–1274 (2014).
Felger, J. C. et al. Molecular signatures of peripheral blood mononuclear cells during chronic interferon-α treatment: relationship with depression and fatigue. Psychol. Med. 42, 1591–1603 (2012).
Savitz, J. et al. Inflammation and neurological disease-related genes are differentially expressed in depressed patients with mood disorders and correlate with morphometric and functional imaging abnormalities. Brain Behav. Immun. 31, 161–171 (2013).
Cole, S. W. Elevating the perspective on human stress genomics. Psychoneuroendocrinology 35, 955–962 (2010).
Huang, D. W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57 (2009).
Cole, S. W., Yan, W., Galic, Z., Arevalo, J. & Zack, J. A. Expression-based monitoring of transcription factor activity: the ℡iS database. Bioinformatics 21, 803–810 (2005).
Acknowledgements
We thank Claire Naylies and Yannick Lippi for their contribution to microarray fingerprints acquisition and microarray data analysis carried out at GeT-TRiX Genopole Toulouse Midi-Pyrénées facility, Dr. J.J. Peroua, Dr. M. Salzmann-Lovato, and the medical assistants for their help with patient recruitment, and the patients for their participation in the study. This work was supported by funds from the French National Research Agency, ANR (ANR-11-JSV1-0006, LC), and the Fondation pour la Recherche Médicale (DPP20151033943, LC).
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Moisan, MP., Foury, A., Dexpert, S. et al. Transcriptomic signaling pathways involved in a naturalistic model of inflammation-related depression and its remission. Transl Psychiatry 11, 203 (2021). https://doi.org/10.1038/s41398-021-01323-9
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41398-021-01323-9