Alteration of transthyretin and thyroxine-binding globulin in major depressive disorder: multiple reaction monitoring-based proteomic analysis

Background Major depressive disorder (MDD), common mental disorder, lacks objective diagnostic and prognosis biomarkers. The objective of this study was to perform proteomic analysis to identify proteins with changed expression levels after antidepressant treatment and investigate differences in protein expression between MDD patients and healthy individuals. Methods A total of 111 proteins obtained from literature review were subjected to multiple reaction monitoring (MRM)-based protein quantitation. Finally, seven proteins were quantified for plasma specimens of 10 healthy controls and 78 MDD patients (those at baseline and at 6 weeks after antidepressant treatment of either selective serotonin reuptake inhibitors (SSRIs) or mirtazapine). Results Among 78 MDD patients, 35 patients were treated with SSRIs and 43 patients were treated with mirtazapine. Nineteen (54.3%) and 16 (37.2%) patients responded to SSRIs and mirtazapine, respectively. Comparing MDD patients with healthy individuals, alteration of transthyretin was observed in MDD (P = 0.026). A few differences were observed in protein levels related to SSRIs treatment, although they were not statistically significant. Plasma thyroxine-binding globulin (TBG) was different between before and after mirtazapine treatment only in responders (P = 0.007). Conclusions In proteomic analysis of plasma specimens from MDD patients, transthyretin and TBG levels were altered in MDD and changed after antidepressant treatment.

endocrine system that include a dysregulated hypothalamic-pituitary-adrenal axis [6,7] and reproductive endocrine changes [8]. However, etiology and consequences of MDD remain unclear due to innate complexity of MDD and multifactored etiologies such as genetic and environmental factors including trauma during childhood or adulthood and stressful life events [9]. Currently there are no firm biomarkers for diagnosis of MDD [10], leading to the lack of objective diagnostic tool and misdiagnosis of MDD, especially in primary care setting [10]. Although at least one-third of patients treated with second-generation antidepressants do not achieve treatment responses [11,12], there is no relevant biomarker to predict antidepressant treatment response either.
A few proteomic studies related to prediction of treatment response and changes by antidepressant treatment have been reported, including increase of pigment epithelium-derived factor level after electroconvulsive therapy (ECT) [27], decrease of proteasome subunit α type-2 in antidepressant responders [28], and association of apolipoprotein A-IV with antidepressant response [29]. However, these studies are focused on changes of proteins according to mixed antidepressants [28,29] and ECT [27], or on proteins involved in immune, endocrine, and metabolic processes [29]. More proteomic studies using precise analytic method to identify predictive protein markers, especially peripheral blood markers, are required because they are non-invasive and practical for clinical use [10].
Therefore, the objective of this study was to perform multiple reaction monitoring (MRM)-based proteomic analysis to identify proteins whose expression levels were changed after antidepressant treatment, that might be able to predict treatment response, and assess differences in protein expression between MDD patients and healthy individuals. We screened 111 proteins using MRM method and performed quantitative proteomic analysis usingstable isotope labeled peptide in plasma from MDD patients. To the best of our knowledge, this is the first report describing changes in protein levels after antidepressant treatment in MDD patients using MRM-based proteomic analysis.

MDD patients fulfilling the Diagnostic and Statistical
Manual of Mental Disorders, Fourth Edition, criteria for major depressive episode who were treated with selective serotonin reuptake inhibitors (SSRIs: escitalopram, fluoxetine, paroxetine, sertraline), or non-SSRI (mirtazapine) were recruited. Diagnosis was confirmed by board certified psychiatrists based on Samsung Psychiatric Evaluation Schedule, case review notes, and SCID (structured clinical interview for DSM-IV) to diagnose depression. A minimum baseline score of 15 for 17-item Hamilton Rating Scale for Depression (HAM-D) was required [30]. Study participants were excluded if they had pregnancy, significant medical conditions, abnormal laboratory baseline values, unstable psychiatric features (e.g., suicidal), history of alcohol or drug dependence, seizures, head trauma with loss of consciousness, neurological illness, or concomitant Axis I psychiatric disorder.
Patients received monotherapy for 6 weeks with one of four commonly used SSRI or mirtazapine antidepressants by clinician's choice. In this study, choice of drug was driven by the preference of the physician, with consideration of anticipated side effects in at-risk individuals. Dose titration was completed within two weeks. Trough plasma samples were drawn at the end of week 6 for plasma drug concentrations. Lorazepam 0.5-1 mg was allowed at bedtime for insomnia. Patients were seen by a psychiatrist, who monitored their adverse events and severity of depression. The HAM-D was administered by a single trained rater every two weeks [30]. The rater and personnel of protein quantitation were blinded to the hypotheses, drug assignment, and HAM-D data. Therapeutic response was defined as 50% or more reduction of HAM-D score by 6 weeks after initiation of antidepressant treatment.
Clinical data of depressant patients including age, gender, HAM-D score, antidepressant, family history, onset age, number of episode, and duration of current episode were collected for each individual. Healthy individuals without known past medical or psychiatric history or family history of MDD were included as controls. This study was approved by Samsung Medical Center Institutional Review Board. Informed consent was obtained from all participants.
Peripheral blood specimens were obtained from patients with MDD at baseline and at 6 weeks after initiation of antidepressant treatment. Patients were categorized into four groups according to antidepressant type used and treatment responsiveness: SSRIs responder, SSRIs nonresponder, mirtazapine responder, and mirtazapine nonresponder. For marker prioritization, we prepared pooled plasma by pooling equal amounts of plasma specimen from five individuals for each group. To investigate whether proteins detected in pooled plasma were detectable in each plasma, we used these 45 specimens individually. For protein quantitation, we used additional 166 plasma specimens at baseline and at 6 weeks after initiation of either SSRIs or mirtazapine treatment from 78 MDD patients and 10 healthy controls (Fig. 1).

Marker surrogate selection
We searched proteins that were either related to depressive disorders or located on antidepressant action sites or metabolic pathways in previous individual studies and meta-analysis. In addition to protein markers, genetic markers, DNA, and mRNA were also selected. We manually reviewed them and obtained proteins encoded by genes as well as proteins described in previous literatures. Among selected proteins, we gave priorities to proteins present in serum or plasma based on public databases including Healthy Human Individuals Integrated Plasma Proteome Database and Sys-Body Fluid database as they were studied much more in literatures. Finally, a total of 111 proteins were selected for MRM-based marker prioritization (Additional file 1: Table S1).

Sample preparation
For marker prioritization, 50 μL of plasma was resolved on 4 ~ 12% NuPAGE gel which was cut into 10 bands and subjected to in-gel tryptic digestion prior to MRM analysis. In-gel digestion was achieved at a 50:1 ratio for 16 h at 37 °C. Tryptic digests were recovered by extraction with 50% ACN/0.1% FA and purified using an OMIX C-18 tip (Agilent, Technologies, Santa Clara, CA, USA).
For protein quantitation, 1 μL of plasma was mixed with 40 μL of 10 mM of dithiothreitol and incubated at 60 °C for 45 min for protein reduction. Next, 5 mM iodoacetamide was added. Samples were incubated at room temperature for 30 min in the dark to induce alkylation followed by digestion with trypsin (1:50) overnight at 37 °C, after which digested samples were dried using a centrifugal evaporator. Each tryptic digest sample was spiked with stable isotope labeled peptide. Samples were dried and reconstituted with 20 μL of 0.1% FA in water.

MRM-based marker prioritization
We generated MRM transitions using MRMPilot ™ v2.0 (AB Sciex, Framingham, MA, USA) against 111 proteins selected by literature review. We then monitored these MRM methods using pooled plasma. We performed MRM analysis using a QTRAP 5500 hybrid triple quadrupole/linear ion trap mass spectrometer (AB Sciex) equipped with a nano-electrospray ion source. MRM mode setting was as follows: curtain gas and spray gas at 20 and 25 psi, respectively; collision gas set on high level; and declustering potential set at 100 V. Among 111 proteins, 33 proteins showed signal-to-noise ratio (S/N) above 8 in pooled plasma specimen. For these 33 proteins, we performed MRM analysis using individual plasma specimens to select marker candidates for protein quantitation using stable isotope labeled peptide. The MRM method was identical to that was used for protein measurement of pooled plasma.

MRM-based protein quantitation
Through marker prioritization, seven proteins were detectable (S/N > 8) in individual plasma. For these seven proteins, we performed protein quantitation using stable isotope labeled peptide. Each tryptic digest sample was spiked with stable isotope labeled peptide. Samples were dried and reconstituted with 20 μL of 0.1% FA in water. MRM analysis was performed using Agilent 6490 Triple Quadrupole mass spectrometer equipped with Agilent 1260 Infinity LC system (Agilent Technologies Inc.,

Patient characteristics and protein quantitation
Clinical and demographic characteristics of MDD patients at baseline and healthy individuals are summarized in Table 2. Age, gender, HAM-D score, onset age, number of episode, duration of current episode, and antidepressant used in study population were investigated. Overall, SSRIs-or mirtazapine-treated patients and healthy individuals and MDD patients showed no statistically significant differences. Responders and nonresponders did not differ in any variables at baseline. Among 111 proteins selected by literature review, 33 proteins showed S/N ratio above 8 in pooled plasma specimen. Among these 33 proteins, seven proteins (AGP1, HGFA, LBP, prothrombin, SeP, TBG, and transthyretin) showed S/N ratio above 8 in individual plasma specimens. For these seven proteins, we performed quantitative analysis using stable isotope labeled peptide in 166 plasma specimens from 78 MDD patients and 10 healthy individuals and compared protein expressions for the following groups: patients vs. controls, responders vs. nonresponders, and at baseline vs. at 6 weeks after treatment.

Differentially expressed proteins between MDD patients and controls
After comparison between patients at baseline and controls, transthyretin showed lower concentration in MDD patients than that in controls (mean, 6338 vs. 7718 fmol/ μL, uni-P = 0.026, Table 3 and Fig. 2a). This statistical significance was maintained after adjusting for age and gender (multi-P = 0.038). Although AGP1 concentration was also lower in MDD patients than that in controls, such difference did not reach statistical significance (7874 vs. 9448 fmol/μL, uni-P = 0.116).

Plasma protein concentrations in MDD patients treated with SSRIs
In patients treated with SSRIs, the concentrations of HGFA and LBP were numerically different between responders and nonresponders at baseline (5.44 vs. 4.00 fmol/μL, uni-P = 0.141, multi-P = 0.193 for HGFA; 28.8 vs. 33.6 fmol/μL, uni-P = 0.362, multi-P = 0.053 for LBP; Table 4 and Fig. 2b). TBG levels were increased in responders after SSRI treatment (118 vs. 129 fmol/μL, uni-P = 0.090). However, these differences according to responsiveness to SSRIs and changes according to SSRIs treatment of plasma protein did not reach statistical significance.

Plasma protein concentrations in MDD patients treated with mirtazapine
In patients treated with mirtazapine, although plasma protein concentration was higher in responders than that in nonresponders for HGFA (4.48 vs. 3.79 fmol/μL, uni-P = 0.191) and lower in responders for prothrombin (1450 vs. 1528 fmol/μL, uni-P = 0.057) at baseline, these differences did not show statistical significance (Table 5 and Fig. 2c).  In analysis for changes according to mirtazapine treatment, TBG was increased after mirtazapine treatment in responders (120 vs. 139 fmol/μL, uni-P = 0.007). On the other hand, TBG did not show statistically significant changes according to mirtazapine treatment in nonresponders.

Discussion
In this study, we performed MRM-based proteomic analysis using peripheral blood specimens at baseline and at 6 weeks after antidepressant treatment in MDD patients and from healthy individuals. We also quantified proteins using stable isotope labeled peptide. We identified a difference in plasma transthyretin concentration between MDD patients and healthy individuals, and a change in plasma TBG concentration related to mirtazapine treatment.
The diagnosis of MDD is primarily based on clinical features, not objective biomarkers. Biomarkers of diagnosis are potential candidates for disease monitoring. Our study showed lower transthyretin level in MDD patients than that in healthy individuals. Transthyretin is a carrier protein which transports thyroxin (T4) and retinol in plasma and CSF [31]. It has been reported that transthyretin is associated with neurodegenerative disorders and psychiatric condition including major depressive disorder [31]. Transthyretin level has been found to be low in CSF from patients with Alzheimer's Disease [31][32][33]. Correlation between plasma transthyretin level and severity and progression of Alzheimer's disease has also been reported [34]. Previous study has also reported that CSF transthyretin level is lower in depressive patients compared to that in healthy controls [35,36]. Serum transthyretin level is also lower in patients with poststroke depression than that in patients without poststroke depression [37]. In preclinical studies using animal model of depression, administration of sodium butyrate, a histone deacetylase inhibitor with antidepressant-like effect, has resulted upregulation of transthyretin RNA level [38,39]. Low plasma transthyretin level in MDD patients in our study could be meaningful in the aspect of that this is the first demonstration using plasma from MDD patients. Although conflicting findings have been reported in preclinical studies [40] and whether transthyretin alteration is a cause or consequence of MDD remains unclear, it has been suggested that low transthyretin level can lower thyroid hormone availability which may lead to depression [35,39].
In analysis for identifying proteins whose levels are changed after antidepressant treatment, TBG level was found to be elevated after mirtazapine treatment only in responders. TBG has not been measured in MDD yet. Its changes in response to antidepressant treatment has not been investigated in previous studies either. Only one study has evaluated TBG as a predictor of perinatal syndromal depression based on relationship between hypothalamic-pituitary-thyroid axis abnormality and pregnancy-related depression [41]. TBG is one of the major thyroid-hormone binding proteins, along with transthyretin and albumin. In the present study, TBG level was increased in responders after treatment with SSRIs (uni-P = 0.090) and transthyretin level was increased in responders after treatment with mirtazapine (uni-P = 0.095). Although these increases did not show statistical significance, they were in line with the difference of transthyretin levels between MDD patients and controls and changes of TBG levels in responders to mirtazapine. These findings suggest that alterations and changes of transportation of thyroid hormones can mediate the effect of antidepressants and clinical manifestation of MDD. These findings have not been reported in previous studies. They need to be confirmed by further studies.
When compared between responders and nonresponders to both SSRIs and mirtazapine, HGFA showed higher concentrations in responders. HGFA is involved in the HGF-Met pathway correlated in tissue protection, regeneration, and anti-fibrosis/inflammation [42]. With regard to MDD, it was suggested that the HGFA concentration correlated with severity and symptoms of depression in previous studies [43,44]. HGFA as a marker that predicts treatment response should be demonstrated further analysis because the correlation between HGFA concentration and treatment responsiveness was not statistically significant in our study.
We performed proteomic analysis using pre-and post-plasma specimens after SSRIs and mirtazapine treatment from MDD patients and healthy individuals. Our results showed alterations and changes of plasma proteins including transthyretin and TBG. These (See figure on next page.) Fig. 2 Plasma levels of thyroxine-binding globulin (TBG) and transthyretin in healthy controls and major depressive disorder (MDD) patients treated with either selective serotonin reuptake inhibitors (SSRIs) or mirtazapine. a Comparison between healthy individuals and patients with MDD. Transthyretin was higher in controls than that in MDD patients. b Plasma protein levels in MDD patients treated with SSRIs. TBG levels were increased in responders after SSRIs treatment. However, the difference did not reach statistical significances. c Plasma protein levels in MDD patients treated with mirtazapine. TBG level was increased after mirtazapine treatment. Transthyretin levels tended to be lower in responders. However, this difference did not reach the statistical significance. NR, nonresponder; R, responder findings are meaningful because there are few studies about changes of protein expression after antidepressant treatment or differences at baseline associated with responsiveness. Our findings have some limitations. First, performed evaluation of our method to detect plasma proteins through several experimental stages. However, we did not validate our results using independent study population. Therefore, findings of this study should be considered as exploratory findings. Further validation study is required.

Conclusions
We identified proteins including transthyretin and TBG, which were altered in MDD after antidepressant treatment using proteomic analysis of plasma specimens. Our findings suggest that alteration and change in transportation of thyroid hormones can be related to the effect of antidepressants and clinical manifestation of MDD.

Supplementary Information
The online version contains supplementary material available at https ://doi. org/10.1186/s1296 7-021-02702 -y.  Table 5 Plasma protein concentrations (fmol/μL) at baseline and at 6 weeks after mirtazapine treatment HAM-D Hamilton Depression Rating Scale, SD standard deviation a P-value from Wilcoxon signed rank test or paired T-test between at baseline and at 6 weeks in each response group b P-value from Mann-Whitney test or T-test between responders and nonresponders at baseline c Multivariable P-value from partial Spearman correlation analysis including clinical variables with P < 0.020 (HAM-D score at baseline, and duration of current episode)

Protein Responder Nonresponder Responsiveness
Mean ± SD Uni-P a Mean ± SD Uni-P a Uni-P b Multi-P c