Clinical and Biological Factors Are Associated with Treatment-Resistant Depression

Abstract

Background: Treatment-resistant depression (TRD) is a debilitating condition associated with unmet clinical needs. Few studies have explored clinical characteristics and serum biomarkers associated with TRD. Aims: We investigated whether there were differences in clinical and biochemical variables between patients affected by TRD than those without. Methods: We recruited 343 patients (165 males and 178 females) consecutively hospitalized for MDD to the inpatient clinics affiliated to the Fondazione IRCCS Policlinico, Milan, Italy (n = 234), and ASST Monza, Italy (n = 109). Data were obtained through a screening of the clinical charts and blood analyses conducted during the hospitalization. Results: TRD versus non-TRD patients resulted to be older (p = 0.001), to have a longer duration of illness (p < 0.001), to be more currently treated with a psychiatric poly-therapy (p < 0.001), to have currently more severe depressive symptoms as showed by the Hamilton Depression Rating Scale (HAM-D) scores (p = 0.016), to have lower bilirubin plasma levels (p < 0.001). In addition, more lifetime suicide attempts (p = 0.035), more antidepressant treatments before the current episode (p < 0.001), and a lower neutrophil to lymphocyte ratio at borderline statistically significant level (p = 0.060) were all associated with the TRD group. Conclusion: We identified candidate biomarkers associated with TRD such as bilirubin plasma levels and NLR, to be confirmed by further studies. Moreover, TRD seems to be associated with unfavorable clinical factors such as a predisposition to suicidal behaviors. Future research should replicate these results to provide robust data in support of the identification of new targets of treatment and implementation of prevention strategies for TRD.

1. Introduction

Major Depressive Disorder (MDD) is a highly prevalent mental disorder worldwide. It is a complex and heterogeneous illness, causing individual suffering, important health care costs, loss of productivity and increased suicide risk [1]. Furthermore, up to 60% of individuals with MDD may show an incomplete response to their first antidepressant treatment, whilst approximately 30% of patients respond poorly to various antidepressant treatments, eventually falling into the condition of treatment-resistant depression (TRD) [2]. Even though variously defined, the most common TRD definition requires a non-response to two or more consecutive antidepressant medications in a single major depressive episode [3]. Available data indicate that TRD differs from MDD without treatment resistance in several aspects. Patients with TRD have higher rates of relapse and account therefore for a greater economic burden than those without TRD [4]. In addition, TRD represents an important risk factor for the occurrence of several medical comorbidities, and it is positively associated with decreased life expectancy as a result of both natural and unnatural causes [5]. A number of clinical variables including older age, comorbid personality and anxiety spectrum disorders, longer duration of illness, more frequent and recurrent depressive episodes, greater severity of depression, suicidal behavior, higher rate of hospitalizations and medical comorbidity (e.g., cardiovascular diseases) were all detected as possible predictors of TRD [6,7]. However, in light of specific sociodemographic and clinical correlates associated with TRD, some authors hypothesized that TRD patients may differ in biological characteristics from non-TRD ones [8]. In support of this speculation in the last years, different research groups have tried to better understand the pathophysiology of TRD, focusing on the identification of candidate biomarkers associated with this condition [9]. Particularly, many studies and meta-analyses reported that, on average, TRD patients have increased serum levels of pro-inflammatory molecules including interleukin-6 (IL-6), tumor necrosis factor (TNF) and C-reactive protein (CRP) [10,11]. Of note, patients affected by MDD seem to benefit less from antidepressant treatments in the presence of hyperactivity of the inflammatory system [12]. In this regard, the neutrophil–lymphocyte ratio (NLR) was suggested as a reproducible biomarker of the degree of systemic inflammation in patients with different medical conditions including psychiatric ones [13]. Higher values of NLR may be associated with a better response to some antidepressants [14]. Moreover, some authors reported that over-inflammation may hamper the effectiveness of antidepressants as a result of lipid peroxidation, impaired mitochondrial function, insulin resistance, oxidative damage to neuron membranes and DNA [15,16]. Indeed, some antioxidant agents, including serum levels of uric acid, albumin and high-density lipoprotein cholesterol (HDL-C) were found to be lower among depressed patients than controls, whilst some oxidative damage products, including malondialdehyde, were higher than in healthy subjects [17]. In addition, bilirubin, an endogenous antioxidant, may be impaired by excessive oxidative stress [18]. It is not surprising that lower total serum bilirubin levels were proposed as a risk factor of vulnerability for different neurological and psychiatric disorders, including MDD [19]. Similarly, in light of the association between obesity and inflammation, metabolic syndrome and related biomarkers (e.g., dysregulated cortisol release) were hypothesized to be related to TRD [20,21].

Despite the available literature, the overall picture of TRD still seems conflicting and fragmentary. Indeed, few studies evaluated concomitantly different clinical variables and laboratory markers, especially among patients with TRD compared to non-TRD ones [22]. Nevertheless, it is noteworthy that most findings, in term of biomarkers in pathophysiology of MDD/TRD, are inconsistent in light of different confounding factors such as sample characteristics, ongoing medications, medical comorbidities and definition of TRD [9]. Thus, in order to gain more information useful for clinical practice and identify any specific feature of TRD patients, this study had the aim to explore the possible roles of several clinical variables and biochemical markers on the occurrence of TRD in a sample of real-world inpatients affected by unipolar MDD.

Therefore, our main research questions were:

• What is the prevalence of TRD in hospitalized patients with MDD?

• Are there some clinical and biological variables for discriminating TRD individuals from depressed patients without this characteristic?

• Is it possible to hypothesize some interactions between the different biological systems implicated in TRD?

Given the exploratory nature of the selected clinical and biomarker panel based on previous available literature [13,17,23], it was expected that TRD patients may differ from non-TRD patients for some of the investigated variables.

2. Materials and Methods

An overall sample of 343 patients with a diagnosis of MDD according to DSM (Diagnostic and Statistical Manual of Mental Disorders) criteria [24] was enrolled among patients admitted to the inpatient clinics affiliated to (1) the Fondazione IRCCS Policlinico, Milan, Italy (n = 234), and (2) ASST Monza, Italy (n = 109). The protocol of this study was approved by the local Ethical Committees (approval number 1789). Clinical and biological data were obtained through a screening of the clinical charts, interviews with patients and their relatives, and blood analyses conducted during the first day of hospitalization. Exclusion criteria were the following: (1) an age < 18 years, (2) peripartum as perinatal depression presents a different clinical presentation associated with specific biological abnormalities [25,26]; (3) pharmacological treatments that could significantly change mood (e.g., steroids); (4) medical comorbidities that could significantly affect mood or biochemical parameters (e.g., rheumatoid arthritis) [27]. In case of psychiatric comorbidity, MDD represented the disorder that affected the patients for a longer time or/and was responsible for more disability. TRD was defined as the presence of at least two prior treatment failures at adequate dose and duration [28]. This criterion was adopted to identify TRD patients in our sample.

The following variables were collected on the first day of hospitalization:

• -. clinical variables: age, age at onset, duration of untreated illness (months), duration of illness (months), duration of hospitalization (days), duration of last antidepressant therapy (months), number of antidepressants assumed before the current episode, the presence of a current depressive episode with mixed features, presence of lifetime psychotic symptoms, presence of lifetime history of substance misuse, presence of lifetime history of misuse of more than one substance, presence of family history of psychiatric disorders, presence of family history of substance misuse, current treatment with more than one psychotropic drug (including benzodiazepines), comorbidity with at least one psychiatric diagnosis, comorbidity with more than one psychiatric diagnosis, presence of lifetime suicidal attempts, number of lifetime suicidal attempts, comorbidity with hypothyroidism, comorbidity with diabetes, comorbidity with obesity, the body mass index (BMI), and scores of Global Assessment of Functioning (GAF), Young Mania Rating Scale (YMRS), Hamilton Depression Rating Scale 21 items (HAM-D), Montgomery–Asberg Depression Rating Scale (MADRS), Brief Psychiatric Rating Scale (BPRS), Hamilton Anxiety Rating Scale (HAM-A);

• -. biological parameters: number of red blood cells (10^12/L), mean corpuscular volume (MCV) (fL), hemoglobin (HB) (g/dL), number of white blood cells (WBC) (10^9/L), number of lymphocytes (10^9/L), neutrophil-to-lymphocyte ratio (NLR), number of platelets (10^9/L), mean platelet volume (MPV) (fL), C-reactive protein (CRP) (mg/dL), total plasmatic proteins (g/dL), albumin (g/dL), bilirubin (mg/dL), uric acid (mg/dL), cholesterol (mg/dL), glycaemia (fasting) (mg/dL), creatine phosphokinase (CPK) (U/L), high-density lipoproteins (HDL) (mg/dL), low-density lipoproteins (LDL) (mg/dL), triglicerydes (mg/dL) and thyroid-stimulating hormone (TSH) (mcU/mL).

Duration of untreated illness was considered as the time elapsing between first episode of MDD and the prescription of a proper pharmacological treatment (antidepressant or quetiapine) [29]. Obesity was defined as a BMI ≥ 30 [30].

With regard to rating scales, GAF measures patients’ global functioning with total scores ranging from 1 to 100 and higher scores indicating less disability [31]. YMRS is a tool used to assess the severity of manic symptoms, and in the case of patients affected by MDD it is useful to identify mixed features (YMRS scores ≥ 10) [32]. Both HAM-D and MADRS measure the severity of depressive symptoms with the first focusing much more on anxiety-somatization aspects (HAM-D ≥ 8 indicating clinically significant depressive symptoms) [33] and the latter more specific to identify core symptoms of MDD (e.g., guilt feelings and anhedonia). MADRS scores ≥ 10 reveal clinically significant depressive symptoms [34]. BPRS assesses the severity of general psychopathology with a score ≥ 31 indicating the presence of clinically significant psychiatric symptoms [35] and, finally, HAM-A measures the severity of anxiety symptoms with a score ≥ 8 requiring clinical attention [36].

Given that a difference of at least 5 points in age is expected between TRD and non-TRD patients according to previous literature [37] and that a p value = 0.05 is considered statistically significant, for a power of 80% a sample of 212 subjects with at least 106 individuals for each subgroup is calculated to be reliable.

Statistical analyses were performed through The Statistical Package for Social Sciences (SPSS) for Windows (version 26.0). Descriptive analyses of the total sample were initially performed, then groups identified according to the presence of TRD were compared using Student’s t-tests for means of quantitative variables and Chi-square tests for qualitative variables (Table 1—clinical variables and Table 2—biological variables). Owing to the large number of factors statistically related to the dependent variable (the presence of TRD) at the univariate analyses, two multiple logistic regression analyses (one for clinical variables—Table 3 and one for biological ones—Table 4) were performed including only the statistically significant variables. This approach was followed as recommended by previous publications on purposeful selection of variables in logistic regression [38] and similarly to the method used by our group in other articles for the analyses of data from samples of patients affected by mood disorders [39]. The choice to initially perform two separate models for clinical and biological variables was driven by the hypothesis that no inter-relation was present between these two sets of factors [40]. A final model entering both clinical and biological variables was then performed to verify the stability of the results of the previous models [40] (Table 5). Finally, correlation analyses (Pearson’s r) were performed to detect eventual interactions among all the independent variables inserted in the final logistic regression model and explaining eventual differences in the findings of separate versus final regression models (Table 3 and Table 4 versus Table 5).

The goodness of fit of the models was assessed by the Omnibus and Hosmer–Lemeshow tests. The level of statistical significance was set at p ≤ 0.05, and confidence intervals at 95% for odds ratios were calculated.

3. Results

The total sample included 343 patients, of whom 165 were males (48.1%) and 178 females (52.1%). The mean age was 49.27 ± 14.98 years, with a minimum of 18 and a maximum of 86. Descriptive analyses and p values of univariate analyses are reported in Table 1 (demographic and clinical variables) and Table 2 (biological variables). With regard to antidepressant therapy assumed during the last depressive episode, 17 patients were treated with tricyclic antidepressants (TCAs), 69 with selective serotonin reuptake inhibitors (SSRIs), 42 with serotonin-norepinephrine reuptake inhibitors (SNRIs), 12 with mirtazapine, 12 with trazodone, 10 with vortioxetine, 7 with bupropion and 13 with other antidepressants. Patients without TRD were more frequently treated with SSRIs than TRD ones (36.8% versus 19.3%, p < 0.05).

First of all, considering demographical and clinical variables, there were not significant differences between TRD and non-TRD patients (p > 0.05) in: age at onset, duration of untreated illness, duration of hospitalization, duration of last antidepressant therapy, presence of mixed features during the current episode, presence of lifetime psychotic symptoms, presence of lifetime history of substance misuse, presence of lifetime history of more than one substance, family history of psychiatric disorders, family history of substance abuse, frequency of current treatment with statins, comorbidity with at least one psychiatric diagnosis, comorbidity with more than one psychiatric diagnosis, history of lifetime suicidal attempts, comorbidity with hypothyroidism, comorbidity with obesity, BMI, comorbidity with diabetes, GAF score, YMRS score, MADRS score, BPRS score, HAM-A score. In contrast, TRD versus non-TRD patients were older (p = 0.001), and they showed: a longer duration of illness (p < 0.001), a higher number of antidepressants assumed before the current episode (p < 0.001), a higher frequency of a current therapy with more than one psychotropic drug (p = 0.001), more lifetime suicidal attempts (p = 0.019) and higher HAM-D scores at the moment of hospitalization (p = 0.016).

With regards to biological parameters, there were no differences between the two groups (p > 0.05) in: number of red blood cells (RBC), hemoglobin, number of white blood cells (WBC), number of lymphocytes, number of platelets, mean platelet volume (MPV), C-reactive protein (CRP), total plasmatic proteins, albumin, uric acid, total cholesterol, HDL LDL, triglycerides, glycaemia (fasting), CPK, TSH. In contrast, TRD versus non-TRD patients showed lower NLR (p = 0.012) and bilirubin (p < 0.001) levels and a smaller MCV at a borderline statistically significant level (p = 0.069).

With regard to the binary logistic regression with clinical variables (Table 3), the goodness-of-fit test (Hosmer and Lemeshow Test: χ² = 6.773, p = 0.561) showed that the model was reliable, allowing for a correct classification of 83.1% of the cases. In addition, the model was overall significant (Omnibus test: χ² = 68.612, df = 5, p < 0.001). TRD patients (versus non-TRD ones) were confirmed to have received more antidepressants before the current episode (p < 0.001) and to show more lifetime suicidal attempts (p = 0.035).

With regard to the binary logistic regression with biological variables including MCV in light of the borderline statistically significant difference between groups (Table 4), the goodness-of-fit test (Hosmer and Lemeshow Test: χ² = 7.440, p = 0.490) showed that the model was reliable, allowing for a correct classification of 66.7% of the cases. In addition, the model was overall significant (Omnibus test: χ² = 14.585, df = 3, p = 0.002). A lower NLR was confirmed to be associated with TRD patients at the borderline statistically significant level (p = 0.060).

The final model (Table 5) confirmed substantially the results of the previous model with the exception of a loss of statistical strength about the association between NLR and TRD and a more robust association between more lifetime suicide attempts and TRD. Duration of illness reached statistical significance, but with an opposite trend with respect to univariate analyses. Of note, correlation analyses revealed a significant inverse interaction between bilirubin levels and number of suicide attempts (r = −0.189, p < 0.01) and between duration of illness and NLR (r = −0.195, p < 0.01). In addition, NLR and bilirubin resulted to have a direct correlation (r = 0.233, p < 0.01).

4. Discussion

Our study focused on the entanglement of clinical variables and biochemical parameters in the occurrence of TRD in quite a large sample of patients hospitalized for MDD. In line with literature [37,41], TRD was present in slightly more than one-third (35.9%) of the recruited sample. In addition, with regard to rating scales, our sample showed scores that are reliable for inpatients with a diagnosis of MDD (HAM-D and MADRS mean total scores with standard deviations corresponding to a moderate-to-severe major depressive episode) [34,42].

Either clinical or biological differences emerged between patients suffering from TRD and non-TRD ones. First of all, it is not surprising that patients without TRD were more frequently treated with SSRIs than TRD ones (36.8% versus 19.3%). Indeed, the majority of the clinical practice guidelines recommend SSRI as first-line treatment thanks to their favourable risk–benefit ratio [43]. On the other hand, SNRIs and TCAs are commonly considered as a second-line option or for patients with more severe depressive symptoms [44]. Specifically, TCAs have a less favourable safety profile with respect to SSRIs and SNRIs, and they are more associated with troublesome side effects (e.g., constipation) that prevent an optimal adherence to pharmacological treatment [45].

With regard to clinical variables, TRD patients showed features typically associated with less clinical stabilization and overall more severity of MDD. Indeed, albeit revealed only by the univariate analyses, TRD versus non-TRD patients had a significantly longer duration of illness and a higher frequency of a current therapy with more than one psychotropic drug. The prescription of a polytherapy and a chronic course of illness are recognized factors associated with poor response to pharmacological treatment in MDD [37,46,47]. Specifically, patients with TRD are often treated with combined treatments (e.g., antidepressants plus atypical antipsychotics or lithium) associated with an increased risk of side effects and less adherence to medical prescription, resulting in a further worsening of symptoms [48]. This is consistent with our result, in line with previous literature, that patients with TRD experience more different antidepressant treatments with respect to non-TRD ones [49,50]. Furthermore, in our sample, a more severe suicidal behaviour resulted to be robustly associated with TRD, as showed by the binary logistic regression models. Several authors underlined the association between suicidal behaviour and poor response to treatment in MDD, suggesting that history of self-harm [50,51,52] or number of suicidal attempts are strong predictors of TRD [53,54,55].

In terms of biological parameters, even though not confirmed by logistic regression, TRD versus non-TRD patients seem to have lower bilirubin levels. Lower bilirubin levels were also found to be correlated with more lifetime suicide attempts, thus confirming that antioxidant defenses can be compromised in severe patients. Based on the authors’ knowledge, this is the first study reporting a potential role between bilirubin levels and MDD severity. Bilirubin has antioxidant and anti-inflammatory properties, and it is assumed that increased oxidative stress might impair its functioning [19]. Although in recent times a growing number of studies have explored the involvement of anti-oxidative defences in the pathophysiology of MDD [56,57,58], there is little evidence concerning the role played by bilirubin, suggesting that lower levels might represent a risk to more severe forms of illness [19,59]. Of note, some authors proposed bilirubin as a promising therapy for the management of neurodegenerative diseases like Parkinson’s Disease [60]. Similarly, a lower NLR resulted to be associated with TRD. NLR is commonly considered as a reliable measure of subclinical inflammation [61], as well as a potential biomarker associated to the severity of MDD and propensity to suicidal behaviour [62,63]. In contrast with our findings, most of the literature data highlighted higher NLR levels in depressed patients compared to healthy controls [64,65]. Similarly, other studies reported higher NLR in suicide attempters than individuals without suicidal behaviour [63,66]. Conversely, there is evidence that antidepressants treatment might influence inflammatory processes by affecting proliferation of lymphocyte, hence decreasing NRL [67,68]. In addition, the finding of a lower NLR in TRD versus non-TRD patients is consistent with data reported by our group on a different sample showing that lower NLR levels are associated with violent suicide attempts [13]. Even though we did not focus on the modalities of self-harm, TRD patients of this sample resulted to have more lifetime suicidal attempts, consistently with the result of a lower NLR with respect to patients without TRD.

Some limitations should be considered in the interpretation of the aforementioned results. First, even though the sample was quite large, the recruitment in only two centers might partly affect the generalization of our results. Second, last pharmacological treatment with different compounds could have influenced either clinical or biological parameters, as non-TRD patients were more often treated with SSRI. Third, the information partly obtained interviewing patients and their relatives (e.g., family history of psychiatric disorders) may not be always accurate. Fourth, with regard to some variables, a lot of data are missing either because this information is not routinely collected in one of the two inpatient clinics or because it cannot be derived from medical records. Multiple imputation was not performed in light of the retrospective naturalistic design of our study.

From a clinical point of view, our data indicate that patients with TRD should be managed considering the increased risk of suicide and a probable weakness of antioxidant systems compared to non-TRD patients. With regard to suicidal risk, clinicians may consider pharmacological approaches, prescribing compounds with a clear anti-suicidal effect (e.g., lithium) [69], psychosocial interventions (e.g., psychoeducation) [70] or neuromodulation techniques (e.g., vagus nerve stimulation) [71]. However, it is necessary to underline that a recent review highlighted that these strategies showed poor efficacy in samples of patients with TRD [53]. Studies assessing the effectiveness of different treatments to mitigate self-harm risk in TRD patients should be therefore absolutely promoted for the best management of subjects affected by this condition. In this context, as regards the enhancement of antioxidative systems, it would be important to understand whether the administration of antioxidant agents such as vitamin E [72] could improve some clinical aspects of patients with TRD.

In light of the previous considerations, the results of the present study corroborate the hypothesis that TRD patients differ from non-TRD ones for some clinical and biological variables, thus potentially providing insight into the pathophysiology and treatment of this severe condition. Future research should replicate these results to provide robust data in support of the identification of new targets of treatment and implementation of prevention strategies for TRD. Furthermore, these data stimulate a greater deepening of the relationships between the different clinical and biological markers in explaining the onset and course of mood disorders as well as the search for treatments that can promote antioxidative defenses in the most serious patients. These research directions will favor an approach to the patients with TRD based on precision medicine.

Author Contributions

Conceptualization, M.B., E.C. and A.D.; methodology, G.N.; software, M.B.; validation, M.B.; formal analysis, M.B., E.C. and A.C (Alice Caldiroli).; investigation, M.C. (Martina Capellazzi), T.S. and C.P.; resources, L.C.; data curation, A.C. (Alessandro Ceresa), C.M.E., A.M.A., I.T., F.L. and M.D.P.; writing—original draft preparation, M.B., E.C. and A.C. (Alice Caldiroli); writing—review and editing, A.C. (Alessandro Ceresa) and C.M.E.; visualization, F.Z.Q. and F.C.; supervision, M.B., A.D. and M.C. (Massimo Clerici); project administration, A.D.; funding acquisition, A.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by a grant (2019-3396–title: “shedding light on treatment-resistant depression and unraveling ketamine’s fast-acting antidepressant mechanisms of action through patient-derived induced pluripotent stem cells”) from the Italian Cariplo Foundation, which has no involvement in any stage of research and dissemination of results.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of the Fondazione IRCCS Ca Granda Ospedale Maggiore Policlinico (protocol code 1789).

Informed Consent Statement

Patient consent was waived due to the retrospective design of this study.

Data Availability Statement

The data presented in this study are available upon request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Tables

Table 1

Clinical variables of the total sample and of the two groups divided according to the presence of TRD.

Variables	Total Samplen = 343	Non-TRDn = 220 (64.1%)	TRDn = 123 (35.9%)	p Value
Age (years)	49.27 (±14.98)	47.20 (±15.34)	52.99 (±13.60)	0.001
Age at onset of MDD (years) (Missing = 35)	37.31 (±16.41)	37.61 (±16.78)	36.79 (±15.79)	0.677
Duration of untreated illness (months) (Missing = 125)	52.69 (±104.44)	48.05 (±102.84)	60.69 (±107.33)	0.390
Duration of illness (months) (Missing = 38)	138.01 (±157.88)	107.90 (±140.99)	190.71 (±172.08)	<0.001
Duration of hospitalization (days) (Missing = 113)	12.90 (±9.08)	12.98 (±8.63)	12.75 (±9.85)	0.856
Duration of last antidepressant therapy (months) (Missing = 168)	9.61 (±19.58)	11.23 (±21.56)	7.12 (±15.88)	0.149
Number of antidepressants assumed before the current episode (Missing = 154)	1.32 (±1.76)	0.68 (±0.85)	2.71 (±2.37)	<0.001
Current depressive episode with mixed features (Missing = 114)	36 (15.7%)	23 (15.9%)	13 (15.5%)	1.000
Presence of lifetime psychotic symptoms (Missing = 116)	51 (22.5%)	28 (19.4%)	23 (27.7%)	0.186
Presence of lifetime history of substance misuse (Missing = 7)	89 (26.5%)	57 (26.3%)	32 (26.9%)	1.000
Presence of lifetime history of misuse of more than one substance (Missing = 7)	23 (6.8%)	17 (7.8%)	6 (5.0%)	0.376
Presence of family history of psychiatric disorders (Missing = 140)	92 (45.3%)	54 (41.9%)	38 (51.4%)	0.241
Presence of family history of substance misuse (Missing = 149)	16 (8.2%)	10 (7.9%)	6 (8.8%)	1.000
Current therapy with more than one psychotropic drug (Missing = 116)	177 (78.0%)	102 (70.8%)	75 (90.4%)	0.001
Current treatment with statins(Missing = 222)	9 (7.4%)	3 (4.2%)	6 (12.0%)	0.160
Comorbidity with at least one psychiatric diagnosis (Missing = 115)	78 (34.2%)	46 (32.4%)	32 (37.2%)	0.474
Comorbidity with more than one psychiatric diagnosis (Missing = 113)	9 (3.9%)	6 (4.2%)	3 (3.5%)	1.000
Presence of lifetime suicidal attempts (Missing = 66)	128 (46.2%)	72 (42.6%)	56 (51.9%)	0.140
Number of lifetime suicidal attempts (Missing = 69)	0.73 (±1.02)	0.61 (±0.88)	0.92 (±1.19)	0.019
Comorbidity with hypothyroidism(Missing = 113)	14 (6.1%)	6 (4.2%)	8 (9.3%)	0.154
Comorbidity with obesity(Missing = 246)	26 (26.8%)	10 (19.6%)	16 (34.8%)	0.111
BMI (Missing = 261)	23.44 (±3.07)	23.21 (±3.11)	23.80 (±3.03)	0.399
Comorbidity with diabetes(Missing = 115)	15 (6.6%)	6 (4.2%)	9 (10.5%)	0.096
GAF score (Missing = 116)	46.48 (±15.30)	47.33 (±15.50)	45.06 (±14.93)	0.280
YMRS score (Missing = 122)	3.59 (±4.75)	3.61 (±4.75)	3.55 (±4.79)	0.925
HAM-D score (Missing = 124)	15.95 (±5.57)	15.25 (±5.28)	17.12 (±5.86)	0.016
MADRS score (Missing = 128)	24.87 (±6.77)	24.50 (±6.87)	25.53 (±6.59)	0.285
BPRS score (Missing = 104)	36.75 (±7.82)	36.47 (±7.83)	37.21 (±7.84)	0.481
HAM-A score (Missing = 133)	8.98 (±4.75)	8.61 (±4.38)	9.62 (±5.31)	0.137

BMI = Body Mass Index; BPRS = Brief Psychiatric Rating Scale; GAF = Global Assessment of Functioning; HAM-A = Hamilton Rating Scale for Anxiety; HAM-D = Hamilton Rating Scale for Depression; MADRS = Montgomery–Asberg Depression Rating Scale; MDD = Major Depressive Disorder; TRD = Treatment-Resistant Depression; YMRS = Young Mania Rating Scale. Note: For some variables, more than 100 values are missing as they are not routinely collected or are not available in one of the involved inpatient clinics. Means for quantitative variables and frequencies for qualitative ones are reported. Standard deviations for quantitative variables and percentages for qualitative variables are reported in brackets. In bold, statistically significant p resulting from χ² or unpaired Student’s t-tests.

Table 2

Biological markers in study sample stratified by TRD.

Variables	Total Samplen = 343	Non-TRDn = 220 (64.1%)	TRDn = 123 (35.9%)	p Value
Number of RBC (10^12/L) (Missing = 44)	4.65 (±0.52)	4.61 (±0.53)	4.65 (±0.51)	0.494
MCV (fL) (Missing = 153)	86.63 (±6.47)	87.29 (±5.94)	85.53 (±7.18)	0.069
HB (g/dL) (Missing = 44)	13.73 (±1.72)	13.81 (±1.65)	13.58 (±1.82)	0.281
Number of WBC (10^9/L) (Missing = 45)	7.55 (±2.72)	7.67 (±2.94)	7.34 (±2.30)	0.321
Number of RBC (10^12/L) (Missing = 44)	4.65 (±0.52)	4.61 (±0.53)	4.65 (±0.51)	0.494
Number of lymphocytes (10^9/L) (Missing = 115)	2.30 (±0.83)	2.24 (±0.85)	2.39 (±0.79)	0.187
NLR (Missing = 114)	2.23 (±1.86)	2.46 (±2.14)	1.89 (±1.26)	0.012
Number of PLT (10^9/L) (Missing = 150)	243.90 (±69.00)	245.13 (±73.65)	241.83 (±60.82)	0.749
MPV (fL) (Missing = 157)	12.11 (±17.88)	12.89 (±22.52)	10.78 (±1.08)	0.438
CRP (mg/dL) (Missing = 250)	1.37 (±2.77)	1.51 (±3.33)	1.16 (±1.69)	0.551
Total plasmatic proteins (g/dL) (Missing = 163)	6.59 (±0.56)	6.59 (±0.61)	6.60 (±0.48)	0.925
Albumin (g/dL) (Missing = 150)	4.19 (±0.40)	4.18 (±0.43)	4.19 (±0.35)	0.872
Bilirubin (mg/dL) (Missing = 63)	0.62 (±0.53)	0.70 (±0.60)	0.49 (±0.36)	<0.001
Uric acid (mg/dL) (Missing = 159)	4.61 (±1.36)	4.62 (±1.44)	4.60 (±1.26)	0.934
Cholesterol (mg/dL) (Missing = 84)	185.06 (±42.61)	181.91 (±42.27)	190.14 (±42.87)	0.131
HDL (mg/dL) (Missing = 220)	52.41 (±18.59)	52.58 (±21.00)	52.11 (±13.86)	0.891
LDL (mg/dL) (Missing = 230)	112.97 (±36.04)	115.75 (±36.07)	108.10 (±35.91)	0.280
Triglicerydes (mg/dL) (Missing = 216)	115.77 (±63.10)	113.38 (±63.40)	119.57 (±63.09)	0.593
Glycaemia (fasting) (mg/dL) (Missing = 53)	92.03 (±19.82)	91.78 (±17.97)	92.48 (±22.77)	0.771
CPK (U/L) (Missing = 169)	116.40 (±165.57)	125.71 (±174.60)	102.90 (±151.70)	0.373
TSH (mcU/mL) (Missing = 126)	1.93 (±1.73)	1.97 (±1.99)	1.87 (±1.20)	0.670

CPK = creatine phosphokinase; CRP = C-reactive protein; HB = haemoglobin; RBC = red blood cells; MCV = mean corpuscular volume; MPV = mean platelet volume; NLR = neutrophil-to-lymphocyte ratio; PLT = platelets; TRD = Treatment-Resistant Depression; TSH = thyroid-stimulating hormone; WBC = white blood cells. Note: For some variables, more than 100 values are missing as they are not routinely collected or are not available in one of the involved inpatient clinics. Means and standard deviations (in brackets) are reported. In bold, statistically significant p resulting from unpaired Student’s t-tests.

Table 3

Summary of the statistics of best-fit binary regression model for quantitative clinical variables.

Variables	B	S.E.	Wald	p	OR	95% CI for OR
Age (years)	0.025	0.013	3.509	0.061	1.025	0.999–1.052
Duration of illness (months)(Missing = 38)	−0.001	0.002	0.772	0.380	0.999	0.996–1.002
Number of antidepressants assumed before the current episode (Missing = 154)	1.110	0.211	27.624	<0.001	3.034	2.006–4.589
Number of lifetime suicidal attempts (Missing = 69)	0.411	0.195	4.431	0.035	1.509	1.029–2.213
HAM-D score (Missing = 124)	0.069	0.039	3.163	0.075	1.071	0.993–1.156

In this analysis, the dependent variable was presence of Treatment-Resistant Depression (TRD) versus absence of TRD. B = regression coefficient; CI = confidence interval; HAM-D = Hamilton Rating Scale for Depression; OR = odds ratio; S.E. = standard error of B; Wald = Wald statistics. In bold, statistically significant p resulting from the model.

Table 4

Summary of the statistics of best-fit binary regression model for quantitative biological variables.

Variables	B	S.E.	Wald	p	OR	95% CI for OR
NLR (Missing = 114)	−0.235	0.125	3.529	0.060	0.790	0.618–1.010
Bilirubin (mg/dL)(Missing = 63)	−0.695	0.411	2.852	0.091	0.499	0.223–1.118
MCV (fL) (Missing = 153)	−0.039	0.025	2.496	0.114	0.962	0.917–1.009

In this analysis, the dependent variable was presence of Treatment Resistant Depression (TRD) versus absence of TRD. B = regression coefficient; CI = confidence interval; MCV = mean corpuscular volume; NLR = neutrophil to lymphocyte ratio; OR = odds ratio; S.E. = standard error of B; Wald = Wald statistics. In bold, statistically significant p resulting from the model.

Table 5

Regression analysis predictors of TRD.

Variables	B	S.E.	Wald	p	OR	95% CI for OR
Age (years)	0.035	0.020	3.233	0.072	1.036	0.997–1.076
Duration of illness (months)(Missing = 38)	−0.005	0.002	4.213	0.040	0.995	0.991–1.000
Number of antidepressants assumed before the current episode (Missing = 154)	1.485	0.318	27.791	<0.001	4.416	2.367–8.238
Number of lifetime suicidal attempts (Missing = 69)	0.662	0.248	7.112	0.008	1.939	1.192–3.154
HAM-D score (Missing = 124)	0.110	0.057	3.698	0.054	1.16	0.998–1.248
NLR (Missing = 114)	−0.294	0.192	2.348	0.125	0.512	0.618–1.085
Bilirubin (mg/dL)(Missing = 63)	−0.337	0.464	0.529	0.467	0.714	0.287–1.772
MCV (fL) (Missing = 153)	−0.076	0.057	1.777	0.183	0.927	0.828–1.037

In this analysis, the dependent variable was presence of Treatment Resistant Depression (TRD) versus absence of TRD. B = regression coefficient; CI = confidence interval; HAM-D = Hamilton Rating Scale for Depression; MCV = mean corpuscular volume; NLR = neutrophil to lymphocyte ratio; OR = odds ratio; S.E. = standard error of B; Wald = Wald statistics. In bold, statistically significant p resulting from the model.

References

1. Patel, V.; Chisholm, D.; Parikh, R.; Charlson, F.; Degenhardt, L.; Dua, T.; Ferrari, A.; Hyman, S.; Laxminarayan, R.; Levin, C. et al. Addressing the burden of mental, neurological, and substance use disorders: Key messages from Disease Control Priorities, 3rd edition. Lancet; 2016; 387, pp. 1672-1685. [DOI: https://dx.doi.org/10.1016/S0140-6736(15)00390-6]

2. Voineskos, D.; Daskalakis, Z.J.; Blumberger, D.M. Management of Treatment-Resistant Depression: Challenges and Strategies. Neuropsychiatr. Dis. Treat.; 2020; 16, pp. 221-234. [DOI: https://dx.doi.org/10.2147/NDT.S198774] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32021216]

3. McIntyre, R.S.; Filteau, M.J.; Martin, L.; Patry, S.; Carvalho, A.; Cha, D.S.; Barakat, M.; Miguelez, M. Treatment-resistant depression: Definitions, review of the evidence, and algorithmic approach. J. Affect. Disord.; 2014; 156, pp. 1-7. [DOI: https://dx.doi.org/10.1016/j.jad.2013.10.043] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24314926]

4. Johnston, K.M.; Powell, L.C.; Anderson, I.M.; Szabo, S.; Cline, S. The burden of treatment-resistant depression: A systematic review of the economic and quality of life literature. J. Affect. Disord.; 2019; 242, pp. 195-210. [DOI: https://dx.doi.org/10.1016/j.jad.2018.06.045]

5. Li, G.; Fife, D.; Wang, G.; Sheehan, J.J.; Bodén, R.; Brandt, L.; Brenner, P.; Reutfors, J.; DiBernardo, A. All-cause mortality in patients with treatment-resistant depression: A cohort study in the US population. Ann. Gen. Psychiatry; 2019; 18, 23. [DOI: https://dx.doi.org/10.1186/s12991-019-0248-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31583010]

6. De Carlo, V.; Calati, R.; Serretti, A. Socio-demographic and clinical predictors of non-response/non-remission in treatment resistant depressed patients: A systematic review. Psychiatry Res.; 2016; 240, pp. 421-430. [DOI: https://dx.doi.org/10.1016/j.psychres.2016.04.034] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27155594]

7. Halaris, A.; Sohl, E.; Whitham, E.A. Treatment-Resistant Depression Revisited: A Glimmer of Hope. J. Pers. Med.; 2021; 11, 155. [DOI: https://dx.doi.org/10.3390/jpm11020155]

8. Klok, M.; van Eijndhoven, P.F.; Argyelan, M.; Schene, A.H.; Tendolkar, I. Structural brain characteristics in treatment-resistant depression: Review of magnetic resonance imaging studies. BJPSych Open; 2019; 5, e76. [DOI: https://dx.doi.org/10.1192/bjo.2019.58]

9. Yang, C.; Wardenaar, K.J.; Bosker, F.J.; Li, J.; Schoevers, R.A. Inflammatory markers and treatment outcome in treatment resistant depression: A systematic review. J. Affect. Disord.; 2019; 257, pp. 640-649. [DOI: https://dx.doi.org/10.1016/j.jad.2019.07.045]

10. Strawbridge, R.; Arnone, D.; Danese, A.; Papadopoulos, A.; Herane Vives, A.; Cleare, A.J. Inflammation and clinical response to treatment in depression: A meta-analysis. Eur. Neuropsychopharmacol.; 2015; 25, pp. 1532-1543. [DOI: https://dx.doi.org/10.1016/j.euroneuro.2015.06.007]

11. Chamberlain, S.R.; Cavanagh, J.; de Boer, P.; Mondelli, V.; Jones, D.N.; Drevets, W.C.; Cowen, P.J.; Harrison, N.A.; Pointon, L.; Pariante, C.M. et al. Treatment-resistant depression and peripheral C-reactive protein. Br. J. Psychiatry; 2019; 214, pp. 11-19. [DOI: https://dx.doi.org/10.1192/bjp.2018.66] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29764522]

12. Rizzo, S.J.S.; Neal, S.J.; Hughes, Z.; Beyna, M.; Rosenzweig-Lipson, S.; Moss, S.J.; Brandon, N. Evidence for sustained elevation of IL-6 in the CNS as a key contributor of depressive-like phenotypes. Transl. Psychiatry; 2012; 2, e199. [DOI: https://dx.doi.org/10.1038/tp.2012.120] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23212583]

13. Capuzzi, E.; Caldiroli, A.; Capellazzi, M.; Tagliabue, I.; Auxilia, A.; Ghilardi, G.; Buoli, M.; Clerici, M. Exploring the role of serum lipid profile and neutrophil-to-lymphocyte ratio in violent suicide attempters: A cross sectional study. CNS Spectr.; 2020; 10, pp. 1-7. [DOI: https://dx.doi.org/10.1017/S1092852920002199] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33298214]

14. Llorca-Bofí, V.; Palacios-Garrán, R.; Routo, D.R.; Buil-Reiné, E.; Adrados-Pérez, M.; Gich, I.; Bioque, M.; Mur, M. High neutrophil-lymphocyte ratio upon admission is associated with better response in psychotic depression. J. Psychiatr. Res.; 2021; 143, pp. 38-42. [DOI: https://dx.doi.org/10.1016/j.jpsychires.2021.08.021]

15. Maes, M.; Galecki, P.; Chang, Y.S.; Berk, M. A review on the oxidative and nitrosative stress (O&NS) pathways in major depression and their possible contribution to the (neuro)degenerative processes in that illness. Prog. Neuropsychopharmacol. Biol. Psychiatry; 2011; 35, pp. 676-692. [DOI: https://dx.doi.org/10.1016/j.pnpbp.2010.05.004]

16. Lindqvist, D.; Dhabhar, F.S.; James, S.J.; Hough, C.M.; Jain, F.A.; Bersani, F.S.; Reus, V.; Verhoeven, J.E.; Epel, E.S.; Mahan, L. et al. Oxidative stress, inflammation and treatment response in major depression. Psychoneuroendocrinology; 2017; 76, pp. 197-205. [DOI: https://dx.doi.org/10.1016/j.psyneuen.2016.11.031]

17. Liu, T.; Zhong, S.; Liao, X.; Chen, J.; He, T.; Lai, S.; Jia, Y. A Meta-Analysis of Oxidative Stress Markers in Depression. PLoS ONE; 2015; 10, e0138904. [DOI: https://dx.doi.org/10.1371/journal.pone.0138904]

18. Dipnall, J.F.; Pasco, J.A.; Berk, M.; Williams, L.J.; Dodd, S.; Jacka, F.N.; Meyer, D. Fusing Data Mining, Machine Learning and Traditional Statistics to Detect Biomarkers Associated with Depression. PLoS ONE; 2016; 11, e0148195. [DOI: https://dx.doi.org/10.1371/journal.pone.0148195]

19. Jayanti, S.; Moretti, R.; Tiribelli, C.; Gazzin, S. Bilirubin and inflammation in neurodegenerative and other neurological diseases. Neuroimmunol. Neuroinflamm.; 2020; 7, pp. 92-108. [DOI: https://dx.doi.org/10.20517/2347-8659.2019.14]

20. Veltman, E.M.; Lamers, F.; Comijs, H.C.; Stek, M.L.; van der Mast, R.C.; Rhebergen, D. Inflammatory markers and cortisol parameters across depressive subtypes in an older cohort. J. Affect. Disord.; 2018; 234, pp. 54-58. [DOI: https://dx.doi.org/10.1016/j.jad.2018.02.080]

21. Moradi, Y.; Albatineh, A.N.; Mahmoodi, H.; Gheshlagh, R.G. The relationship between depression and risk of metabolic syndrome: A meta-analysis of observational studies. Clin. Diabetes Endocrinol.; 2021; 7, 4. [DOI: https://dx.doi.org/10.1186/s40842-021-00117-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33648597]

22. Sowa-Kućma, M.; Styczeń, K.; Siwek, M.; Misztak, P.; Nowak, R.J.; Dudek, D.; Rybakowski, J.K.; Nowak, G.; Maes, M. Lipid Peroxidation and Immune Biomarkers Are Associated with Major Depression and Its Phenotypes, Including Treatment-Resistant Depression and Melancholia. Neurotox. Res.; 2018; 33, pp. 448-460. [DOI: https://dx.doi.org/10.1007/s12640-017-9835-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29103192]

23. Buoli, M.; Caldiroli, L.; Guenzani, D.; Carnevali, G.S.; Cesari, M.; Turolo, S.; Barkin, J.L.; Messa, P.; Agostoni, C.; Vettoretti, S. et al. Associations Between Cholesterol and Fatty Acid Profile on the Severity of Depression in Older Persons With Nondialysis Chronic Kidney Disease. J. Ren. Nutr.; 2021; 31, pp. 537-540. [DOI: https://dx.doi.org/10.1053/j.jrn.2020.08.017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33386204]

24. American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders; 5th ed. American Psychiatric Association: Arlington, VA, USA, 2013.

25. Serati, M.; Redaelli, M.; Buoli, M.; Altamura, A.C. Perinatal Major Depression Biomarkers: A systematic review. J. Affect. Disord.; 2016; 193, pp. 391-404. [DOI: https://dx.doi.org/10.1016/j.jad.2016.01.027]

26. Serati, M.; Esposito, C.M.; Grassi, S.; Bollati, V.; Barkin, J.L.; Buoli, M. The Association between Plasma ERVWE1 Concentrations and Affective Symptoms during Pregnancy: Is This a Friendly Alien?. Int. J. Environ. Res. Public Health; 2020; 17, 9217. [DOI: https://dx.doi.org/10.3390/ijerph17249217]

27. Ingegnoli, F.; Schioppo, T.; Ubiali, T.; Ostuzzi, S.; Bollati, V.; Buoli, M.; Caporali, R. Patient Perception of Depressive Symptoms in Rheumatic Diseases: A Cross-sectional Survey. J. Clin. Rheumatol.; 2020; 28, pp. e18-e22. [DOI: https://dx.doi.org/10.1097/RHU.0000000000001564]

28. Gaynes, B.N.; Lux, L.; Gartlehner, G.; Asher, G.; Forman-Hoffman, V.; Green, J.; Boland, E.; Weber, R.; Randolph, C.; Bann, C. et al. Defining treatment-resistant depression. Depress Anxiety; 2020; 37, pp. 134-145. [DOI: https://dx.doi.org/10.1002/da.22968]

29. Altamura, A.C.; Buoli, M.; Albano, A.; Dell’Osso, B. Age at onset and latency to treatment (duration of untreated illness) in patients with mood and anxiety disorders: A naturalistic study. Int. Clin. Psychopharmacol.; 2010; 25, pp. 172-179. [DOI: https://dx.doi.org/10.1097/YIC.0b013e3283384c74]

30. Iodice, S.; Ceresa, A.; Esposito, C.; Mucci, F.; Conti, D.; Pergoli, L.; Tarantini, L.; Vigna, L.; Bollati, V.; Buoli, M. et al. The Independent Role of Body Mass Index (BMI) and Severity of Depressive Symptoms on Biological Changes of Women Affected by Overweight/Obesity. Int. J. Environ. Res. Public Health; 2021; 18, 2923. [DOI: https://dx.doi.org/10.3390/ijerph18062923]

31. Moos, R.H.; Nichol, A.C.; Moos, B.S. Global Assessment of Functioning ratings and the allocation and outcomes of mental health services. Psychiatr. Serv.; 2002; 53, pp. 730-737. [DOI: https://dx.doi.org/10.1176/appi.ps.53.6.730]

32. Buoli, M.; Caldiroli, A.; Caletti, E.; Zugno, E.; Altamura, A.C. The impact of mood episodes and duration of illness on cognition in bipolar disorder. Compr. Psychiatry; 2014; 55, pp. 1561-1566. [DOI: https://dx.doi.org/10.1016/j.comppsych.2014.06.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25011691]

33. Altamura, A.C.; Dell’Osso, B.; Buoli, M.; Zanoni, S.; Mundo, E. Intravenous augmentative citalopram versus clomipramine in partial/nonresponder depressed patients: A short-term, low dose, randomized, placebo-controlled study. J. Clin. Psychopharmacol.; 2008; 28, pp. 406-410. [DOI: https://dx.doi.org/10.1097/JCP.0b013e31817d5931] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18626267]

34. Hawley, C.J.; Gale, T.M.; Sivakumaran, T. Hertfordshire Neuroscience Research group. Defining remission by cut off score on the MADRS: Selecting the optimal value. J. Affect. Disord.; 2002; 72, pp. 177-184. [DOI: https://dx.doi.org/10.1016/S0165-0327(01)00451-7]

35. Leucht, S.; Kane, J.M.; Kissling, W.; Hamann, J.; Etschel, E.; Engel, R. Clinical implications of Brief Psychiatric Rating Scale scores. Br. J. Psychiatry; 2005; 187, pp. 366-371. [DOI: https://dx.doi.org/10.1192/bjp.187.4.366]

36. Matza, L.S.; Morlock, R.; Sexton, C.; Malley, K.; Feltner, D. Identifying HAM-A cutoffs for mild, moderate, and severe generalized anxiety disorder. Int. J. Methods Psychiatr. Res.; 2010; 19, pp. 223-232. [DOI: https://dx.doi.org/10.1002/mpr.323]

37. Di Salvo, G.; Cuniberti, F.; Bianco, M.; Ramello, D.; Maina, G.; Rosso, G. Features and pharmacotherapy of treatment-resistant depression: An observational study on hospitalized patients. J. Psychopathol.; 2020; 26, pp. 257-263. [DOI: https://dx.doi.org/10.36148/2284-0249-396]

38. Bursac, Z.; Gauss, C.H.; Williams, D.K.; Hosmer, D.W. Purposeful selection of variables in logistic regression. Source Code Biol. Med.; 2008; 3, 17. [DOI: https://dx.doi.org/10.1186/1751-0473-3-17]

39. Buoli, M.; Cesana, B.M.; Bolognesi, S.; Fagiolini, A.; Albert, U.; Di Salvo, G.; Maina, G.; de Bartolomeis, A.; Pompili, M.; Palumbo, C. et al. Factors associated with lifetime suicide attempts in bipolar disorder: Results from an Italian nationwide study. Eur. Arch. Psychiatry Clin. Neurosci.; 2021; [DOI: https://dx.doi.org/10.1007/s00406-021-01343-y]

40. Heinze, G.; Wallisch, C.; Dunkler, D. Variable selection-A review and recommendations for the practicing statistician. Biom. J.; 2018; 60, pp. 431-449. [DOI: https://dx.doi.org/10.1002/bimj.201700067]

41. Zhdanava, M.; Pilon, D.; Ghelerter, I.; Chow, W.; Joshi, K.; Lefebvre, P.; Sheehan, J.J. The Prevalence and National Burden of Treatment-Resistant Depression and Major Depressive Disorder in the United States. J. Clin. Psychiatry; 2021; 82, 20m13699. [DOI: https://dx.doi.org/10.4088/JCP.20m13699]

42. Hamilton, M. Development of a rating scale for primary depressive illness. Br. J. Soc. Clin. Psychol.; 1967; 6, pp. 278-296. [DOI: https://dx.doi.org/10.1111/j.2044-8260.1967.tb00530.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/6080235]

43. Gabriel, F.C.; de Melo, D.O.; Fráguas, R.; Leite-Santos, N.C.; Mantovani da Silva, R.A.; Ribeiro, E. Pharmacological treatment of depression: A systematic review comparing clinical practice guideline recommendations. PLoS ONE; 2020; 15, e0231700. [DOI: https://dx.doi.org/10.1371/journal.pone.0231700] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32315333]

44. Bauer, M.; Pfennig, A.; Severus, E.; Whybrow, P.C.; Angst, J.; Möller, H.J. on behalf of the Task Force on Unipolar Depressive Disorders. World Federation of Societies of Biological Psychiatry (WFSBP) guidelines for biological treatment of unipolar depressive disorders, part 1: Update 2013 on the acute and continuation treatment of unipolar depressive disorders. World J. Biol. Psychiatry; 2013; 14, pp. 334-385. [DOI: https://dx.doi.org/10.3109/15622975.2013.804195] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23879318]

45. Von Wolff, A.; Hölzel, L.P.; Westphal, A.; Härter, M.; Kriston, L. Selective serotonin reuptake inhibitors and tricyclic antidepressants in the acute treatment of chronic depression and dysthymia: A systematic review and meta-analysis. J. Affect. Disord.; 2013; 144, pp. 7-15. [DOI: https://dx.doi.org/10.1016/j.jad.2012.06.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22963896]

46. Fabbri, C.; Hagenaars, S.P.; John, C.; Williams, A.T.; Shrine, N.; Moles, L.; Hanscombe, K.B.; Serretti, A.; Shepherd, D.J.; Free, R.C. et al. Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts. Mol. Psychiatry; 2021; 26, pp. 3363-3373. [DOI: https://dx.doi.org/10.1038/s41380-021-01062-9]

47. Murphy, J.A.; Sarris, J.; Byrne, G.J. A Review of the Conceptualisation and Risk Factors Associated with Treatment-Resistant Depression. Depress. Res. Treat.; 2017; 2017, 4176825. [DOI: https://dx.doi.org/10.1155/2017/4176825]

48. Taylor, R.W.; Marwood, L.; Oprea, E.; Deangel, V.; Mather, S.; Valentini, B.; Zahn, R.; Young, A.H.; Cleare, A.J. Pharmacological Augmentation in Unipolar Depression: A Guide to the Guidelines. Int. J. Neuropsychopharmacol.; 2020; 23, pp. 587-625. [DOI: https://dx.doi.org/10.1093/ijnp/pyaa033] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32402075]

49. Döme, P.; Kunovszki, P.; Takács, P.; Fehér, L.; Balázs, T.; Dede, K.; Mulhern-Haughey, S.; Barbreau, S.; Rihmer, Z. Clinical characteristics of treatment-resistant depression in adults in Hungary: Real-world evidence from a 7-year-long retrospective data analysis. PLoS ONE; 2021; 16, e0245510. [DOI: https://dx.doi.org/10.1371/journal.pone.0245510]

50. Kautzky, A.; Dold, M.; Bartova, L.; Spies, M.; Kranz, G.S.; Souery, D.; Montgomery, S.; Mendlewicz, J.; Zohar, J.; Fabbri, C. et al. Clinical factors predicting treatment resistant depression: Affirmative results from the European multicenter study. Acta Psychiatr. Scand; 2019; 139, pp. 78-88. [DOI: https://dx.doi.org/10.1111/acps.12959]

51. Balestri, M.; Calati, R.; Souery, D.; Kautzky, A.; Kasper, S.; Montgomery, S.; Zohar, J.; Mendlewicz, J.; Serretti, A. Socio-demographic and clinical predictors of treatment resistant depression: A prospective European multicenter study. J. Affect. Disord.; 2016; 189, pp. 224-232. [DOI: https://dx.doi.org/10.1016/j.jad.2015.09.033]

52. Fornaro, M.; Fusco, A.; Novello, S.; Mosca, P.; Anastasia, A.; De Blasio, A.; Iasevoli, F.; de Bartolomeis, A. Predictors of Treatment Resistance Across Different Clinical Subtypes of Depression: Comparison of Unipolar vs. Bipolar Cases. Front. Psychiatry; 2020; 11, 438. [DOI: https://dx.doi.org/10.3389/fpsyt.2020.00438] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32670098]

53. Bergfeld, I.O.; Mantione, M.; Figee, M.; Schuurman, P.R.; Lok, A.; Denys, D. Treatment-resistant depression and suicidality. J. Affect. Disord.; 2018; 235, pp. 362-367. [DOI: https://dx.doi.org/10.1016/j.jad.2018.04.016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29665520]

54. Congio, A.C.; Norcia, M.; Urbano, M.R.; Verri, W.A.; Vargas Nunes, S.O. Association of clinical features and biomarkers with treatment-resistant depression. Neurol. Psychiatry Brain Res.; 2020; 36, pp. 32-38. [DOI: https://dx.doi.org/10.1016/j.npbr.2020.02.004]

55. Sagud, M.; Mihaljevic-Peles, A.; Uzun, S.; Cusa, B.V.; Kozumplik, O.; Kudlek-Mikulic, S.; Mustapic, M.; Barisic, I.; Muck-Seler, D.; Pivac, N. The lack of association between components of metabolic syndrome and treatment resistance in depression. Psychopharmacology; 2013; 230, pp. 15-21. [DOI: https://dx.doi.org/10.1007/s00213-013-3085-x]

56. Black, C.N.; Bot, M.; Scheffer, P.G.; Cuijpers, P.; Penninx, B.W.J.H. Is depression associated with increased oxidative stress? A systematic review and meta-analysis. Psychoneuroendocrinology; 2015; 51, pp. 164-175. [DOI: https://dx.doi.org/10.1016/j.psyneuen.2014.09.025]

57. Camkurt, M.A.; Findikli, E.; Izci, F.; Kurutaş, E.B.; Tuman, T.C. Evaluation of malondialdehyde, superoxide dismutase and catalase activity and their diagnostic value in drug naïve, first episode, non-smoker major depression patients and healthy controls. Psychiatry Res.; 2016; 238, pp. 81-85. [DOI: https://dx.doi.org/10.1016/j.psychres.2016.01.075]

58. Hacimusalar, Y.; Eşel, E. Suggested Biomarkers for Major Depressive Disorder. Arch. Neuropsychiatry; 2018; 55, pp. 280-290. [DOI: https://dx.doi.org/10.5152/npa.2017.19482]

59. Peng, Y.F.; Xiang, Y.; Wei, Y.S. The significance of routine biochemical markers in patients with major depressive disorder. Sci. Rep.; 2016; 6, 34402. [DOI: https://dx.doi.org/10.1038/srep34402]

60. Jayanti, S.; Moretti, R.; Tiribelli, C.; Gazzin, S. Bilirubin: A Promising Therapy for Parkinson’s Disease. Int. J. Mol. Sci.; 2021; 22, 6223. [DOI: https://dx.doi.org/10.3390/ijms22126223]

61. Faria, S.S.; Fernandes, P.C.; Silva, M.J.B.; Lima, V.C.; Fontes, W.; Freitas, R.; Eterovic, A.K.; Forget, P. The neutrophil-to-lymphocyte ratio: A narrative review. Ecancermedicalscience; 2016; 10, 702. [DOI: https://dx.doi.org/10.3332/ECANCER.2016.702]

62. Sunbul, E.A.; Sunbul, M.; Yanartas, O.; Cengiz, F.; Bozbay, M.; Sari, I.; Gulec, H. Increased Neutrophil/Lymphocyte Ratio in Patients with Depression is Correlated with the Severity of Depression and Cardiovascular Risk Factors. Psychiatry Investig.; 2016; 13, pp. 121-126. [DOI: https://dx.doi.org/10.4306/pi.2016.13.1.121] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26766954]

63. Velasco, Á.; Rodríguez-Revuelta, J.; Olié, E.; Abad, I.; Fernández-Peláez, A.; Cazals, A.; Guillaume, S.; De La Fuente-Tomás, L.; Jiménez-Treviño, L.; Gutiérrez, L. et al. Neutrophil-to-lymphocyte ratio: A potential new peripheral biomarker of suicidal behavior. Eur. Psychiatry; 2020; 63, E14. [DOI: https://dx.doi.org/10.1192/j.eurpsy.2019.20] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32093807]

64. Brinn, A.; Stone, J. Original research: Neutrophil–lymphocyte ratio across psychiatric diagnoses: A cross-sectional study using electronic health records. BMJ Open; 2020; 10, e036859. [DOI: https://dx.doi.org/10.1136/bmjopen-2020-036859] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32690528]

65. Mazza, M.G.; Lucchi, S.; Tringali, A.G.M.; Rossetti, A.; Botti, E.R.; Clerici, M. Neutrophil/lymphocyte ratio and platelet/lymphocyte ratio in mood disorders: A meta-analysis. Prog. Neuro-Psychopharmacol. Biol. Psychiatry; 2018; 84, pp. 229-236. [DOI: https://dx.doi.org/10.1016/j.pnpbp.2018.03.012]

66. Ekinci, O.; Ekinci, A. The connections among suicidal behavior, lipid profile and low-grade inflammation in patients with major depressive disorder: A specific relationship with the neutrophil-to-lymphocyte ratio. Nord. J. Psychiatry; 2017; 71, pp. 574-580. [DOI: https://dx.doi.org/10.1080/08039488.2017.1363285]

67. Demircan, F.; Gözel, N.; Kılınç, F.; Ulu, R.; Atmaca, M. The Impact of Red Blood Cell Distribution Width and Neutrophil/Lymphocyte Ratio on the Diagnosis of Major Depressive Disorder. Neurol. Ther.; 2016; 5, pp. 27-33. [DOI: https://dx.doi.org/10.1007/s40120-015-0039-8]

68. Öztürk, A.; Sahan, E.; Mirçik, A.B.; Deveci, E.; Yilmaz, O.; Kirpinar, I. Mean platelet volume and neutrophil to lymphocyte ratio decrease in patients with depression with antidepressant treatment. Arch. Clin. Psychiatry; 2019; 46, pp. 9-13. [DOI: https://dx.doi.org/10.1590/0101-60830000000184]

69. Guzzetta, F.; Tondo, L.; Centorrino, F.; Baldessarini, R.J. Lithium treatment reduces suicide risk in recurrent major depressive disorder. J. Clin. Psychiatry; 2007; 68, pp. 380-383. [DOI: https://dx.doi.org/10.4088/JCP.v68n0304]

70. Wasserman, D.; Carli, V.; Iosue, M.; Javed, A.; Herrman, H. Suicide prevention in psychiatric patients. Asia Pac Psychiatry; 2021; 13, e12450. [DOI: https://dx.doi.org/10.1111/appy.12450]

71. Bottomley, J.M.; LeReun, C.; Diamantopoulos, A.; Mitchell, S.; Gaynes, B.N. Vagus nerve stimulation (VNS) therapy in patients with treatment resistant depression: A systematic review and meta-analysis. Compr. Psychiatry; 2019; 98, 152156. [DOI: https://dx.doi.org/10.1016/j.comppsych.2019.152156]

72. Manosso, L.M.; Camargo, A.; Dafre, A.L.; Rodrigues, A.L.S. Vitamin E for the management of major depressive disorder: Possible role of the anti-inflammatory and antioxidant systems. Nutr. Neurosci.; 2020; [DOI: https://dx.doi.org/10.1080/1028415X.2020.1853417] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33314993]