4. Results
The study found further associations between inflammation and depression. Inflammatory markers (IL-6, IL-1β, and TNF-α) were measured in a group of 31 patients and compared with a control group without psychiatric disorders or associated inflammatory pathologies of 19 individuals. The analysis of the distribution of patients in the study group provides an insight into the diversity of the sample in terms of gender and background. The gender distribution in the study group reflects potential differences in access to care or the prevalence of the health conditions studied.
Table 1.
Distribution of patients according to gender, BMI and level of inflammation.
Table 1.
Distribution of patients according to gender, BMI and level of inflammation.
| |
BMI |
| GENDER |
< 18.5 kg/m2 (Underweight) |
18.5 – 25 kg/m2 (Normal) |
25 – 30 kg/m2 (Overweight) |
BMI ≥ 30 kg/m2 (Obesity) |
| Men |
2 (6.45%) |
4 (12.90%) |
3 (9.67%) |
1 (3.22%) |
| Women |
5 (16.12%) |
10 (32.25%) |
3 (9.67%) |
3 (9.67%) |
| The presence of inflammation |
2 (6.45%) |
4 (12.90%) |
2 (6.45%) |
2 (6.45%) |
These data offer an intriguing perspective on the distribution of patients by gender, Body Mass Index (BMI), and the presence of elevated inflammatory markers. Interpreting these data brings to the forefront complex aspects that can provide a solid basis for further research and a deeper understanding of these relationships. We will explore these elements without resorting to enumerations, using appropriate academic language.
Firstly, we observe significant variations in the body composition of subjects based on BMI and gender. Regarding men, a significant majority fall into the normal BMI category. This observation suggests a relatively uniform distribution of BMIs among men in this sample. In contrast, women exhibit a greater diversity in BMIs, which may reflect not only biological variability but also behavioural influences on body weight.
Table 2.
Distribution of patients according to gender, family medical history and inflammatory markers.
Table 2.
Distribution of patients according to gender, family medical history and inflammatory markers.
| Gender |
Family medical history |
No. of Patients |
Increased inflammatory markers |
| Men |
Present |
3 (9.67%) |
1 (3.22%) |
| Absent |
7 (22.58%) |
4 (12.90%) |
| Women |
Present |
9 (29.03%) |
5 (16.12%) |
| Absent |
12 (38.70%) |
0 (0.00%) |
From a broader perspective, when analysing the presence of elevated inflammatory markers, we can easily observe that these markers are significantly present across all BMI categories. This indicates an important aspect of this analysis: inflammation appears to be a relevant issue among these patients, regardless of BMI. However, it is important to mention that inflammation can have multiple origins, such as underlying medical conditions, genetic factors, or behavioural factors, and its intensity can vary from one individual to another.
The relevant data presented in
Table 2 highlight a series of noteworthy aspects. Firstly, it is evident that patients were categorized based on gender and the presence or absence of hereditary antecedents. This may reflect a patient classification approach based on biological factors and family history. In this context, there is significant diversity in the distribution of patients across the four categories.
Furthermore, the presence or absence of elevated inflammatory markers in these groups is notable. This may suggest complex relationships between genetic predisposition, family history, and the inflammatory status of patients. However, we must be aware that inflammation is a multifactorial process and can have multiple causes, including environmental factors and lifestyle. This could be an important direction for future research and, overall, these data represent just a "window" into a much more complex medical reality.
In the graph below, we have the correlation between the inflammatory marker IL-1 and the HAM-D test score (
Figure 1). As we can see, we have a Pearson correlation coefficient r = -0.3720, which indicates a negative correlation but with weak statistical significance (p > 0.005), from which we cannot conclude that these variables are correlated and statistically significant.
Analysing the next graph (
Figure 2), which represents the correlation between the inflammatory marker IL-6 and the HAM-D test score, we have a Pearson correlation coefficient r = 0.580, which represents a reasonable positive correlation and statistical significance where p = 0.0390 < 0.05, indicating that the two variables are correlated. The results indicate a statistically significant relationship between the inflammatory marker IL-6 and the severity of depression, as measured by HAM-D scores. This suggests that inflammation, as indicated by IL-6 levels, could play a role in the pathogenesis or severity of depression.
Next,
Figure 3 shows the correlation between the inflammatory marker IL-6 and the HAM-D test in the control group. We can observe that the Pearson correlation coefficient is r = 0.008 and p = 0.0021. An r value of 0.008 indicates a very weak correlation between IL-6 and HAM-D scores. Practically, this coefficient close to zero suggests that there is no significant linear relationship between IL-6 levels and the severity of depression as measured by the HAM-D test in the control group. The p-value of 0.0021 is well below the standard threshold for statistical significance. The mentioned p-value shows that, although the correlation coefficient suggests a lack of linear correlation between the two variables, the statistical test results are significant, suggesting that any observed association between IL-6 and HAM-D scores is unlikely to be due to chance.
It is important to understand that a small p-value in the context of a correlation coefficient very close to zero can be somewhat confusing. The statistical significance indicated by the p-value refers to testing the hypothesis that there is no relationship between the two variables. In this case, the small p-value might indicate the detection of a statistical difference, but the near-zero correlation coefficient suggests that the relationship between the variables is practically negligible.
We can state that the inflammatory marker IL-6 can be considered an indicator in assessing depression, supported by the fact that the two are in a reasonable positive correlation, with p = 0.0021, reaching the statistical significance threshold. In this situation, the inflammatory hypothesis in depression is confirmed, with a very statistically significant correlation, so that IL-6 values increase as the intensity of depression increases.
Figure 3.
Correlation between the HAM-D and IL-6 - control lot.
Figure 3.
Correlation between the HAM-D and IL-6 - control lot.
In
Figure 3, we have a Pearson correlation coefficient of r = 0.06310, which indicates that there is no correlation between the two variables, and a p-value > 0.005, which does not attest to statistical significance.
Figure 4 represents the correlation between the inflammatory marker TNF-α and the HAM-D test score.
Regarding study hypothesis 1, we can state that inflammation is associated to some extent with the highest HAM-D test scores. This hypothesis is confirmed only in the case of IL-6, while for IL-1β and TNF-α, the correlations were not statistically significant.
We further examined the extent to which study hypothesis 4 is confirmed, a hypothesis proposing that marital status influences the depression-inflammation correlation.
Figure 5.
Pearson correlations for IL-1 and HAM-D variables in married patients.
Figure 5.
Pearson correlations for IL-1 and HAM-D variables in married patients.
In married patients, where we correlated the HAM-D score with the IL-1 value, we have a Pearson correlation coefficient of r = -0.535, indicating a moderate negative correlation and statistical significance with p = 0.018262 < 0.05, attesting that the two variables are interconnected. The results indicate a statistically significant relationship between IL-1 inflammatory marker levels and the severity of depression as measured by HAM-D scores in married patients. The moderate negative correlation between IL-1 and HAM-D suggests that higher levels of IL-1 might be associated with lower severity of depressive symptoms in this subgroup of patients, or that more severe depression is associated with lower levels of IL-1.
Similarly, we have a strong negative correlation in divorced and cohabiting patients.
Figure 6.
Pearson correlations for IL-1 and HAM-D variables in divorced patients.
Figure 6.
Pearson correlations for IL-1 and HAM-D variables in divorced patients.
We identified a Pearson correlation coefficient of r = -0.944 for divorced patients (
Figure 6) and r = -0.809 for those cohabiting with partners, indicating a strong association but without statistical significance, with p = 0.214061 > 0.05 for divorced patients and p = 0.400018 for cohabiting patients where IL-1 was measured in conjunction with the HAM-D score. Analyzing the results presented for the correlations between IL-1 inflammatory marker levels and HAM-D scores by marital status, we have two different sets of Pearson coefficients and p values for two subgroups: divorced patients and those cohabiting with partners.
For divorced patients: the Pearson coefficient r = -0.944 indicates a very strong negative correlation between IL-1 and the HAM-D score. This suggests that in the group of divorced patients, higher IL-1 levels are associated with lower HAM-D scores and vice versa, signifying a strong inverse relationship between inflammation and depression severity. The p-value of 0.214061 is greater than the standard threshold for statistical significance (0.05), indicating that although the correlation is strong, it is not statistically significant. This means we cannot exclude the possibility that the observed relationship is due to chance.
For cohabiting patients: the Pearson coefficient r = -0.809 also shows a strong negative correlation between IL-1 and the HAM-D score, suggestive of an inverse association between inflammation and depression severity in this group as well. The p-value of 0.400018, like the divorced group, indicates a lack of statistical significance for the observed correlation. The p value > 0.05 suggests that the observed strong relationship cannot be considered statistically significant, raising the possibility that it is the result of random variability.
The results indicate strong negative correlations between IL-1 and HAM-D scores for both divorced and cohabiting patients, suggesting an inverse association between inflammation levels and depression severity in these groups. However, the lack of statistical significance in both cases indicates that these strong correlations may not be robust or reflect a causal relationship.
Figure 7.
Pearson correlations for IL-1 and HAM-D variables in cohabiting patients.
Figure 7.
Pearson correlations for IL-1 and HAM-D variables in cohabiting patients.
Regarding the association of IL-1 with HAM-D scores and marital status, we obtained negative correlations in all three categories. In this case, we can state that the inflammation-depression relationship, although present, is inversely proportional and confirms the hypothesis that marital status influences the inflammation-depression correlation, with inflammation viewed from the perspective of IL-1 values.
While the correlation is negative and reached statistical significance in married patients, it was not statistically significant in divorced and cohabiting patients.
We further analysed the correlation between IL-6 and HAM-D scores in the three categories of patients: married, divorced, and cohabiting. In the case of married patients (
Figure 8), we identified a Pearson correlation coefficient of r = -0.534, indicating a moderate negative correlation and statistical significance, with p = 0.018443 < 0.05, showing that the two variables are correlated.
Although in the case of the patients from
Figure 8, we identified a good but negative correlation between IL-6 and HAM-D scores, in the case of divorced and cohabiting patients, the correlations are strong and positive, with a good degree of association, with a Pearson correlation coefficient of r = 0.972 for divorced patients and r = 0.642 for cohabiting patients.
The correlation is stronger in divorced patients than in cohabiting patients, being an almost perfectly positive correlation.
Figure 9.
Pearson correlations for IL-6 and HAM-D variables in divorced patients.
Figure 9.
Pearson correlations for IL-6 and HAM-D variables in divorced patients.
In both situations (
Figure 8 and
Figure 9), the threshold of statistical significance was not reached, in divorced patients, the p-value was 0.149645, and in cohabiting patients p = 0.55621.
Figure 10.
Pearson correlations for IL-1 and HAM-D variables in cohabiting patients.
Figure 10.
Pearson correlations for IL-1 and HAM-D variables in cohabiting patients.
To identify and study the relationship between TNF-α and HAM-D scores in the three groups, the Pearson correlation coefficient was calculated, which in this case shows extreme values with different correlations. For married patients, the Pearson coefficient is r = 0.009, suggesting that there is no correlation between the two variables, with a p-value of 0.400018 > 0.05, which is not statistically significant.
The results suggest that, in married patients, there is no significant or practically relevant relationship between TNF-α levels and the severity of depression as evaluated by HAM-D scores. The absence of a significant correlation between these two variables indicates that, at least in this subgroup of patients, TNF-α does not appear to be a predictive marker or directly associated with the severity of depression.
Figure 11.
Pearson correlations for TNF-α and HAM-D variables in married patients.
Figure 11.
Pearson correlations for TNF-α and HAM-D variables in married patients.
For divorced patients (
Figure 12), the Pearson correlation coefficient is r = 0.901, indicating a strong, almost perfect positive correlation and a good degree of association between TNF-α levels and HAM-D scores, with a p-value of 0.004158 < 0.005, confirming that the two variables are correlated.
On the other side, in
Figure 12, is the group of divorced patients, where we identified a Pearson correlation coefficient value of r = -0.109, indicating a weak negative correlation, and a p-value of 0.93047, which is greater than the statistical significance threshold of 0.05.
Figure 13.
Pearson correlations for TNF-α and HAM-D variables in cohabiting patients.
Figure 13.
Pearson correlations for TNF-α and HAM-D variables in cohabiting patients.
In the study hypothesis where we assumed that patients from urban areas would specifically have elevated TNF-α levels, we calculated the mean with the corresponding standard deviation and analysed the minimum and maximum components. Here, the mean in the urban area is lower than in the rural area. The hypothesis is refuted, but still, the maximum value is higher for the urban environment than for the rural one (
Table 3).
Similarly, we tested and confirmed the hypothesis that the intensity of depression in women is higher than in men. We calculated the mean HAM-D score for women and men, along with the standard deviation and minimum and maximum components. As shown in the table below (
Table 4), the hypothesis is confirmed, with women having an average HAM-D score of 30.83 with a standard deviation of 3.1, compared to an average score of 30.05 for men with a standard deviation of 2.99.
Next, we evaluated the relationship between the intensity of depression measured by the HAM-D scale and the value of inflammatory markers IL-1, IL-1β, and TNF-α, based on the gender of the patients, to see the extent to which study hypotheses 1 and 2, established earlier, are confirmed.
The Pearson correlations applied to the group of women (
Figure 14) with major recurrent depressive disorder, in terms of the variables IL-1 and HAM-D, indicate a moderate, satisfactory correlation, with r = 0.372536716, which reaches the threshold of statistical significance p = 0.03902557 < 0.05.
We evaluated the relationship between the intensity of depression in women (
Figure 14) and the intensity of inflammation in terms of IL-6 values, applying Pearson correlations, and we obtained a strong positive correlation (r = 0.580456) and a very statistically significant threshold (p = 0.000619 < 0.05). The results indicate a clear and statistically significant relationship between the intensity of inflammation and the severity of depression among the women studied. This finding is consistent with the theory that inflammatory processes may contribute to the development or exacerbation of depressive symptoms.
The strong positive correlation and high statistical significance suggest that measuring IL-6 levels could provide valuable insights into the biological mechanisms underlying depression and could help in identifying women at higher risk of severe depression due to inflammation.
Figure 15.
Pearson correlations in the case of women from the perspective of IL-6 and HAM-D variables.
Figure 15.
Pearson correlations in the case of women from the perspective of IL-6 and HAM-D variables.
In the case of the association between the intensity of depression and the degree of inflammation associated with TNF-α (
Figure 16), the correlation is insignificant, minor (r = 0.063099591), without reaching the established threshold of statistical significance p < 0.01 (p = 0.735949118). Regarding the inflammation-depression correlation, at least in the group of women, for the TNF-α component, hypothesis 1 is not confirmed.
From the association between IL-1 and the HAM-D score, a moderate correlation is observed in
Figure 17 (r = 0.392211), but it is statistically significant with p = 0.047507. In other words, as the HAM-D score increases in men with recurrent major depressive disorder, higher values of IL-1 are encountered, which supports our hypothesis.
Similar to the case of IL-1, we applied Pearson correlations in evaluating the degree of correlation between IL-6 and the intensity of associated depression (HAM-D score). This time, the threshold for statistical significance is very high (p = 0.00897), and the correlation is also strong and satisfactory (r = 0.546145), indicating a strong link between depression and inflammation, with inflammation being expressed through high IL-6 values.
Figure 18.
Pearson correlations in the case of men from the perspective of IL-6 and HAM-D variables.
Figure 18.
Pearson correlations in the case of men from the perspective of IL-6 and HAM-D variables.
Finally, the relationship between inflammation and depression was evaluated in the male group, considering the HAM-D score and TNF-α value as variables. The result was not significant, with a small correlation (r = 0.072294), lacking statistical significance (p = 0.725618 > 0.05).
As observed in the figure below (
Figure 19), regarding TNF-α, study hypothesis 1 is not confirmed for either men or women, as there is no correlation between the intensity of depression and inflammation. The strongest correlations were found between HAM-D and IL-6 in both men and women, and for IL-1, there are satisfactory correlations with the intensity of depression in both men and women.
Analyzing the above graphs (
Figure 18 and
Figure 19) from the perspective of the statistical analysis performed, we can conclude that study hypothesis 1 is partially confirmed and that depression is specifically associated with an inflammatory status, as evidenced by the statistically significant correlations between high HAM-D scores and high IL-1 and IL-6 values.
Regarding study hypothesis 2, where we assumed that there are stronger correlations between inflammation and depression in women than in men, as observed in the table below (Table 5), this aspect is confirmed only in the case of IL-6, where the value of r in women (r = 0.580456) is higher than in men (r = 0.546145) (
Table 5).
Although there are numerous scales for measuring depression and the associated risk of suicide, none of them provide a score that is sufficiently reliable or clinically useful for accurately estimating a patient's real risk of self-harm.
Healthcare professionals should consider the social constructs and norms that influence how suicide is addressed and indirectly minimized using these terms to describe suicidal ideation.
In this study, we paid particular attention to the clinical interview, anamnesis, and psychological evaluation using the Columbia scale or the HAM-A anxiety scale to identify elements that could be indicative of suicide risk. However, in the patients selected to participate in this scientific research, we did not identify such elements.
As the study aims at evaluating the correlation between the inflammatory markers present in the study group with the severity of depressive symptoms assessed, a regresion analysis was performed to determine if there are statistically significant differences between the mean scores of both groups. Also, the analysis of variance ANOVA was used to analyse the impact of multiple independent variables (inflammatory markers) on dependant variable (HAM-D scores). The p-values obtained through ANOVA indicate which markers are statistically significant in predicting depression severity and explains the relationship between the inflammatory markers and depression in multiple groups.
The regression model analysis performed with control group data suggests that IL-1 and TNF significantly predict HAM-D scores, suggesting a corelation between these two inflammatory markers and depression severity, whereas IL-6 does not appear to be a significant predictor in this group.
Table 6.
Control group ANOVA.
Table 6.
Control group ANOVA.
| ANOVA a |
| Model |
Sum of Squares |
df |
Mean Square |
F |
Sig. |
| 1 |
Regression |
23,538 |
3 |
7,846 |
2,901 |
,069 b
|
| Residual |
40,567 |
15 |
2,704 |
|
|
| Total |
64,105 |
18 |
|
|
|
| a. Dependent Variable: HAMD |
| b. Predictors: (Constant), TNF, IL6, IL1 |
The ANOVA analysis for the study group suggests that the regression model the predictors, respectively the inflammatory markers IL-1, IL-6 and TNF, significantly explain the variance in the HAM-D scores.
Table 7.
Study group ANOVA.
Table 7.
Study group ANOVA.
| ANOVA a |
| Model |
Sum of Squares |
df |
Mean Square |
F |
Sig. |
| 1 |
Regression |
132,352 |
3 |
44,117 |
7,183 |
,001 b
|
| Residual |
165,841 |
27 |
6,142 |
|
|
| Total |
298,194 |
30 |
|
|
|
| a. Dependent Variable: HAMD |
| b. Predictors: (Constant), TNF, IL6, IL1 |
To better observe the relationship between the presence of inflammatory markers and demographic data with the severity of depression assessed, we analyse the correlation matrix of all variables for the study group (
Table 7). While there are no statistically significant correlations observed in regard to age, gender and area of residence to HAM-D scores, we also do not find a significant correlation of TNF marker to the variables provided. However, we have strong positive correlation between HAM-D scores and IL-6 (0.581, p = 0.001), suggesting that higher IL6 levels are associated with higher HAM-D scores. Also, the analysis provides that negative correlation (-0.409, p = 0.022) between IL-1 and HAM-D scores, suggest that higher IL-1 levels are associated with lower depression scores.
Table 8.
Correlation matrix all variables.
Table 8.
Correlation matrix all variables.
| Correlations |
| |
Area |
Status |
Gender |
Age |
HAMD |
TNF |
| Area |
Pearson Correlation |
1 |
,555 ** |
-,065 |
-,708 ** |
-,048 |
,027 |
| Sig. (2-tailed) |
|
,001 |
,728 |
,000 |
,796 |
,887 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| Status |
Pearson Correlation |
,555 ** |
1 |
,226 |
-,495 ** |
,211 |
,096 |
| Sig. (2-tailed) |
,001 |
|
,222 |
,005 |
,256 |
,608 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| Gender |
Pearson Correlation |
-,065 |
,226 |
1 |
-,095 |
,075 |
,232 |
| Sig. (2-tailed) |
,728 |
,222 |
|
,611 |
,687 |
,208 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| Age |
Pearson Correlation |
-,708 ** |
-,495 ** |
-,095 |
1 |
-,154 |
-,262 |
| Sig. (2-tailed) |
,000 |
,005 |
,611 |
|
,408 |
,154 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| HAMD |
Pearson Correlation |
-,048 |
,211 |
,075 |
-,154 |
1 |
,063 |
| Sig. (2-tailed) |
,796 |
,256 |
,687 |
,408 |
|
,736 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| TNF |
Pearson Correlation |
,027 |
,096 |
,232 |
-,262 |
,063 |
1 |
| Sig. (2-tailed) |
,887 |
,608 |
,208 |
,154 |
,736 |
|
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| IL6 |
Pearson Correlation |
-,011 |
,106 |
-,013 |
-,134 |
,581 ** |
,012 |
| Sig. (2-tailed) |
,954 |
,569 |
,947 |
,472 |
,001 |
,947 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
| IL1 |
Pearson Correlation |
-,153 |
-,047 |
,195 |
,082 |
-,409 * |
-,094 |
| Sig. (2-tailed) |
,411 |
,803 |
,292 |
,661 |
,022 |
,616 |
| N |
31 |
31 |
31 |
31 |
31 |
31 |
For the purpose of our study, we look at the confounding variables to accurately evaluate the correlation of the depression severity and the levels of inflammation present: age could be a confounding variable, as older people tend to have higher levels of inflammatory markers compared to younger patients, which could be impacting the correlation. Also, gender differences are relevant, where women show stronger correlations between inflammatory levels and depression than men, whereas marital status is impactful, since divorced patients show a significant positive correlation between IL-6 and depressive scores, while married people did now emphasize such relationship. Also, the presence or absence of comorbidities or use of medication can influence the levels of inflammatory markers and thus confound the relationship to depression scores. The interpretation of regression coefficients (
Table 8) suggests that IL-6 has a strong positive impact on HAM-D scores, the higher IL-6 levels, the higher the depression severity. Another significant coefficient is observed at IL-1, with moderate negative impact on HAM-D scores, the higher IL-1 levels, the lower depression severity. Age, status, living area and TNF do not suggest statistically significant relationship to HAM-D scores.
Table 9.
Regression coefficients all variables.
Table 9.
Regression coefficients all variables.
| Coefficients a |
| Model |
Unstandardized Coefficients |
Standardized Coefficients |
t |
Sig. |
| B |
Std. Error |
Beta |
| 1 |
(Constant) |
33,047 |
4,864 |
|
6,794 |
,000 |
| Age |
-,067 |
,064 |
-,233 |
-1,040 |
,309 |
| Gender |
,478 |
1,050 |
,072 |
,455 |
,653 |
| Status |
1,240 |
,962 |
,236 |
1,289 |
,210 |
| Area |
-2,545 |
1,518 |
-,393 |
-1,677 |
,107 |
| TNF |
,000 |
,000 |
-,069 |
-,441 |
,663 |
| IL6 |
,005 |
,002 |
,462 |
3,095 |
,005 |
| IL1 |
-,003 |
,001 |
-,389 |
-2,578 |
,017 |