A correction has been applied to this article in:
Corrigendum: Isobaric Tags for Relative and Absolute Quantitation Identification of Blood Proteins Relevant to Paroxetine Response in Patients With Major Depressive Disorder
Chin-Chuen Lin
Hung Su2
Jentaie Shiea
Tiao-Lai Huang- 1Department of Psychiatry, Kaohsiung Chang Gung Memorial Hospital, Chang Gung University College of Medicine, Kaohsiung, Taiwan
- 2Department of Chemistry, National Sun Yat-sen University, Kaohsiung, Taiwan
- 3Genomic and Proteomic Core Laboratory, Department of Medical Research, Kaohsiung Chang Gung Memorial Hospital, Kaohsiung, Taiwan
Objectives: Isobaric tags for relative and absolute quantitation (iTRAQ) is a proteomic investigation that could be utilized for rapid identification and quantification of proteins, which we would use to identify differentially expressed proteins in treatment responsive patients with major depressive disorder (MDD).
Methods: Six treatment responsive patients of MDD were recruited, and their peripheral blood mononuclear cell (PBMC) were collected before and after 4 weeks of paroxetine treatment. iTRAQ and Mascot search engine were used to detect differentially expressed proteins, which were then validated by Western blot.
Results: Two thousand one hundred and fifty three proteins were screened, and seven proteins showed differences of more than two-fold and 62 proteins with a differences of less than two-fold. Six proteins with commercially available antibodies were identified, and were validated by Western blot in 10 paroxetine responsive MDD patients. Putative hydroxypyruvate isomerase (HYI), eukaryotic translation initiation factor 4H (eIF4H), and RNA binding motif 8A (RBM8A) had statistically significant differences before and after treatment in the validation. Data are available via ProteomeXchange with identifier PXD028947.
Conclusions: By using iTRAQ and Western blot, we were able to identify HYI, eIF4H, and RAM8a to be the potential predictors of paroxetine treatment response in patients with MDD. This finding could help establish future individualized medicine.
Introduction
Major depression disorder (MDD) is a mental illness characterized by pathological low mood, loss of interest or pleasure, and pessimistic thoughts that negatively impact the quality of life of a significant proportion of the world's population, and has a lifetime incidence ranging from 4.4 to 20% (1). In several decades, many studies have tried to identify the underlying pathogenic mechanisms for MDD (2). However, the pathogenesis of the MDD and related biological regulatory mechanisms are still largely unknown. Therefore, utilization of new technologies to search for potential biomarkers in MDD is mandated.
After the successful deconstruction of the human genome, it was first thought that decoding the gene would be sufficient to unravel the mystery of human brain. However, it was found that genes alone do not necessarily represent the phenotypes, thus entered the post-genome exploration, where the human proteomics program aims at the systematic study of protein expressions.
Earlier proteomic studies frequently utilized cerebrospinal fluid (CSF) or postmortem brain tissues to identify the disease relevant molecules in MDD. The known pathogenesis of MDD involves many brain regions, but most studies focus on one regulatory pathway at a time (3, 4). In clinical population, CSF is challenging to collect, let alone brain tissues. Therefore, investigations of peripheral blood of patients of MDD allowed easier access and flexibility, such as comparing the proteomic changes before and after antidepressant treatment (5). Peripheral blood mononuclear cell (PBMC) has been suggested to be an ideal peripheral meidrum in searching of biomarker (5). Functional lymphatic vessels found in dural sinuses (6) suggested a possible crosstalk between central nervous system (CNS) and the peripheral (7). PBMC has been utilized in both animal (5) and clinical studies of MDD (8, 9).
In the past, we analyzed the peripheral blood of patients with MDD, using acid hydrolysis and modern mass spectrometry, and found that the ratio of transferrin to fibrinogen was significantly different between patients and healthy controls (10, 11). The same technique has been successfully applied to serum samples of patients with schizophrenia, before and after antipsychotic treatment. In that study we found that the serum IgG content changed significantly after risperidone, a frequently prescribed antipsychotic, treatment (12). Furthermore, we also observed that disulfide-isomerase A3 and F-actin-capping protein subunit beta in platelet samples of patients with MDD have higher performance than those of healthy controls using two-dimensional differential gel electrophoresis (2D-DIGE) (13).
2D-DIGE, however, is limited by the range of isoelectric point and molecular weight, and the development of isobaric tags for relative and absolute quantitation (iTRAQ) could be utilized for rapid identification and quantification of proteins. Wang et al. used the iTRAQ method to study relevant proteins in human serum and quickly screened 472 proteins, of which 154 were found to be differentially expressed between patients with depression and healthy conrtols (14). In 2015, Han et al. used iTRAQ technology to study hippocampal postsynaptic density-associated proteins in animal models, screened 1,500 proteins, and found that 74 membrane proteins were differentially expressed (15). Palmfeldt et al. investigated the prefrontal cortices of rat models of depression with iTRAQ and found several pathways that may contribute to the development of new therapeutic directions (16).
Another clinical challenge for managing MDD is that at least a third of patients experienced inadequate response from antidepressant treatment (17). Pharmacogenetics could help identify the most suitable antidepressant based on an individual's genetic profile on pharmacokinetics and pharmacodynamics (18), though individual genotyping remains expensive. Commercial iTRAQ for discovery of target proteins is also expensive, but if the discovered peripheral target proteins are validated, the cost of measuring only identified peripheral markers becomes more affordable.
In this study, iTRAQ was used to analyze PBMC in patients with MDD, before and after 4 weeks of antidepressant paroxetine, a selective serotonin reuptake inhibitor (SSRI), treatment, in order to identify proteins relevant to paroxetine treatment response.
Materials and Methods
Study Samples
The patients were recruited from a tertiary medical center over an 1-year period. Their age was between 20 and 65 years old, and they had no systemic diseases such as diabetes, hypertension, and hyperthyroidism. They also had no smoking or alcohol addiction. All participants had the ability to sign the informed consent, after the investigators thoroughly explained the details of the study. The diagnosis of MDD was established by the same senior psychiatrist, according to the Diagnostic and Statistical Manual of Mental Disorders (DSM)-5 criteria. Patients were free of antidepressants or other medications for at least 1 month prior to entering the study. The severity of MDD was assessed with 17-item Hamilton Depression Rating Scale (HAM-D) (19). All patients had scores greater or equal to 18 prior to treatment. All patients received 4 weeks of paroxetine treatment. Treatment response was defined by HAM-D score reduction of >50%. PBMC of patients were collected in the morning, after fasting for 6 h, before and after 4 weeks of antidepressant treatment. Only those showing treatment response were analyzed (20–22). Written informed consent was provided by all participants after the content and context of the study was fully explained. The institutional review board (IRB) of Chang Gung Memorial Hospital approved the study design (IRB number 201601596B0).
Isobaric Tags for Relative and Absolute Quantification Labeling
Equal amounts of samples at baseline were pooled, and equal amounts of the samples after treatment were pooled, to minimize the effect of individual variation. One hundred micrograms of protein was reduced, blocked on cysteine, alkylated, and subsequently digested with trypsin overnight at 37°C. The resulting products were than labeled with isobaric tags using a 4-plex iTRAQ kit. Isobaric tags 114 and 115 were used to label lithium responsive groups before and after treatment, respectively, and isobaric tags 116 and 117 were used to label paroxetine responsive groups before and after treatment, respectively. A similar design decision in which only reagents 114 and 116 were used had been seen before to cut down experimental cost (23).
Peptide Fractionation
The peptides labeled with different channels of isobaric tags were pooled and desalted by using Sep-Pak C18 cartridges (Waters). The desalted peptides were dried under vacuum and then re-suspended in 0.5% trifluoroacetic acid for further fractionized by reversed phase C18 under high pH condition (Pierce High pH Reversed-Phase Peptide Fractionation Kit, Thermo Fisher Scientific, USA) as manufacturer's suggested protocol. The resulting bound peptides were eluted in 1% ammonia solution with various amounts of acetonitrile (5–50 %) resulting in about 10 fractions. Each eluted fraction was then dried under vacuum and re-suspended in 0.1% formic acid solution for the liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) analysis.
Tandem Mass Spectrometry-Based Protein Identification
Tandem mass spectrometry analysis was performed on Q ExactiveTM HF mass spectrometer (Thermo Fisher, San Jose, California, USA) in combination with a Thermo ScientificTM UltiMateTM 3000 RSLCnano HPLC system. The peptide mixtures were directly loaded onto a 50-cm analytic column (EASY-Spray™ C18 Column), and separated by a gradient with gradually increased of buffer B (80% acetonitrile in 0.1% formic acid) at a flow rate of 250 nl/min over about 165 min. The peptide spectra were acquired in the positive ion mode with a data-dependent acquisition. The top abundant 15 precursor ions within 375-1,400 m/z scan range were then dynamically selected for further fragmented in high collision dissociation (HCD) mode with normalized collision energy set to 33 ± 1. In Full MS scan, the resolution was set to 60,000 at m/z 200, AGC target to 3e6, maximum inject time to 50 ms. In MS/MS scan, the resolution was set to 15,000, with AGC target to 5e4, and maximum inject time set to 100 ms. The release of dynamic exclusion of the selected precursor ions was set to 20 s.
Mass Spectrometry-Based Proteomics Data Analysis by Proteome Discoverer Software and Mascot Search Engine
The MS files were uploaded into Proteome Discoverer (version 2.1, Thermo Scientific) to generate peak list for protein identification using Mascot search algorithm (version 2.5, Matrix Science) against Swiss-Prot human protein database. For protein identification, the following parameteres were used for database search: carbamidomethylation at Cys as the fixed modification, oxidation at Met, acetylation at protein N-terminus, iTRAQ-labeled at peptide N-terminus and K residue as dynamic modifications, maximum missing cleavage sites with 2, 10 ppm for MS tolerance, and 0.02 Da for MS/MS tolerance. The acceptable false discovery rate of peptide and protein identifications was set to <1%. A minimum of two unique peptides and quantification with two spectrum ratio counts from the unique peptides were required to improve the confidence of the protein identification. The median value of the spectrum ratios was calculated as protein abundance. The global median normalization was applied to recalculate the protein abundance to generate the normalized protein ratios to reduce the system error from sample preparation in each experiment. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (24) partner repository with the dataset identifier PXD028947 and 10.6019/PXD028947.
Western Blotting Validation
The volume of the loading buffer was one-fourth of the total volume. The proteins were denatured at 95°C for 10 min. Following SDS-PAGE, the proteins were transferred to PVDF membranes. After blocking in 2% BSA in PBST for 1 h at room temperature, the PVDF membrane was incubated for 10 h at 4°C with the following primary antibodies: anti-hydroxypyruvate isomerase (HYI) antibody (Abcam, Cambridge, UK), anti-alpha cardiac muscle 1 actin (ACTC1) antibody (R&D Systems, Minneapolis, MN, USA), anti-eukaryotic translation initiation factor 4H (eIF4H) antibody (Thermo Fisher Scientific, Waltham, MA, USA), anti-adenylate kinase isoenzyme 1 (AK1) antibody (Thermo Fisher Scientific, Waltham, MA, USA), anti-mitochondrial dynamin-like 120 kDa protein (OPA1) antibody (R&D Systems, Minneapolis, MN, USA), and anti-RNA binding motif 8A (RBM8A) antibody (Thermo Fisher Scientific, Waltham, MA, USA). Anti-cytochrome c oxidase (COX IV) antibody was used as controls (25). After three washes with PBST, the membrane was incubated with anti-goat or anti-rabbit secondary antibody (1:10,000) at room temperature for 1 h, and then the membrane was washed another five times with PBST (5 min each). The Western blot bands were photographed with UVP BioSpectrum 810 Imaging System (Thermo Fisher Scientific, Waltham, MA, USA). The relative intensities of the protein bands were calculated with ImageJ software (National Institutes of Health), using background subtraction method (26).
Protein-Protein Interactions
Protein-protein interactions were analyzed using BioGRID database (27). The interested protein names were inputed in the main page (https://thebiogrid.org/), searching with “By Protein/Gene” and “Homo sapiens.”
Statistical Analysis of iTRAQ Data
The statistical method utilizing peptide-level intensities of iTRAQ reporter ions was used to determine which proteins are significantly changed (28).
Results
Study Participants
We recruited six patients with MDD showing response after 4 weeks of paroxetine treatment. Their demographic data were summarized in Table 1. Samples from six patients were analyzed using iTRAQ-based quantitative proteomics.
TABLE 1
Table 1. Demographic data of paroxetine responsive patients with MDD.
iTRAQ-Based Quantitative Proteomics
Protein and peptide identifications with false determinate rate (FDR) < 0.01 and with at least one unique peptide screened 2,153 proteins. Seven proteins showed more than two-fold difference before and after treatment, and 62 proteins showed differences less than two-fold (Table 2). Six proteins with significant fold changes and commercially available antibodies were identified for further Western blot validation, and they are: putative hydroxypyruvate isomerase (HYI, fold change −311.80), alpha cardiac muscle 1 actin (ACTC1, fold change −14.27), eukaryotic translation initiation factor 4H (eIF4H, fold change −6.42), adenylate kinase isoenzyme 1 (AK1, fold change −3.64), mitochondrial dynamin-like 120 kDa protein (OPA1, fold change −3.55), and RNA binding motif 8A (RBM8A, fold change 2.000). Their representative MS/MS spectra annotated with b and y ions identified are shown in Supplementary Figures 1–6.
TABLE 2
Table 2. Western blot validations of potential biomarkers for paroxetine treatment response in MDD patients.
Western Blot Validation
Ten pairs of PBMC samples from 10 MDD patients with paroxetine response were analyzed with Western blot analysis (Figure 1). The blot intensities were quantified using ImageJ (Supplementary Table 2). The relative intensities were then calculated with background subtraction method (26) by dividing the intensity of target band by the control (COX-IV) intensity. The relative intensities of HYI, ACTC1, eIF4H, AK1, OPA1, and RBM8A before and after 4 weeks of paroxetine treatment are summarized in Supplementary Table 3. Using paired t-test, HYI, eIF4H, and RBM8A showed statistically significant difference before and after paroxetine treatment (p = 0.005, 0.019, and 0.013, respectively). Using Wilcoxon signed-rank test, HYI, eIF4H, and RBM8A showed statistically significant difference before and after paroxetine treatment (p = 0.013, 0.022, and 0.022, respectively) (Table 2).
FIGURE 1
Figure 1. Western blots of 10 MDD patients with paroxetine reponse, before and after antidepressant treatment. A, after treatment; ACTC1, alpha cardiac muscle 1 actin; AK1, adenylate kinase isoenzyme 1; B, before treatment; COX-IV, cytochrome c oxidase subunit IV, as control; EIF4H, eukaryotic translation initiation factor 4H; HYI, putative hydroxypyruvate isomerase; MDD, major depressive disorder; OPA1, mitochondrial dynamin-like 120 kDa protein; RBM8A, RNA binding motif 8A.
Protein-Protein Interactions
eIF4H (https://thebiogrid.org/113297/summary/homo-sapiens/eif4h.html) showed 88 unique interactors, among them eIF4A1. eIF4A1 (https://thebiogrid.org/108289/summary/homo-sapiens/eif4a1.html) had 241 unique interactors, eIF4A3 included. RBM8a (https://thebiogrid.org/115265/summary/homo-sapiens/rbm8a.html) interacted with eIF4A3, among 171 unique interactors, showing that eIF4H and RBM8a can be linked via eukaryotic translation initiation factor family. Those html pages were converted to PDF using a free online converter (https://www.sejda.com/html-to-pdf), and the interested interactors were marked using PDF-XChange Viewer (Version 2.5: Build 322.10). Those html and PDF files were uploaded to a public data depository.
Peptide-Level Intensities of iTRAQ Reporters
A two-fold cut off was initially used, but after the western blots data we were suggested to evaluate the candidates using a statistical method utilizing peptide-level intensities of iTRAQ reporters to determine which proteins are significantly changed (28). This method could not be applied to HYI, eIF4H, and RBM8A due to their low numbers of peptide spectrum matches (and they all have only one unique peptide), but ACTC1, AK1, and OPA1 all had p values lower than 0.05 (p = 1.33 x 10−55, 0.00162, and 0.0125, respectively). p-values for proteins with two-fold increases or decreases were also added to Supplementary Table 1.
Discussion
The major finding of this study was that by using iTRAQ, we were able to identify three proteins: putative hydroxypyruvate isomerase (HYI), eukaryotic translation initiation factor 4H (eIF4H), and RNA binding motif 8A (RBM8A), which were validated by Western blot, as the potential biomarkers for paroxetine treatment response in MDD patients.
It is the first time HYI and eIF4H were found to associate with psychiatric disorders. HYI catalyzes the reversible isomerization between hydroxypyruvate and 2-hydroxy-3-oxopropanoate (tartronate semialdehyde). HYI had not previously been associated with human diseases. eIF4H, in contrast, had been associated with various diseases, most prominently cancer. In a proteomic investigation of attached (surviving) and detached (dying) cancer cells, eIF4H was identified among six consistently differentially expressed proteins, and silencing it had an anti-proliferative effect (29). In lung adenocarcinoma, eIF4H was the target of miRNA519d to inhibit cell proliferation and invasion (30). eIF4H is overexpressed in lung cancer cells and can be used to predict the response of chemotherapy (31). eIF4h has been found to be associated with neurodegenerative diseases as well. In the post-mortem tissues of patients with amyotrophic lateral sclerosis and frontotemporal degeneration, eIF4H was downregulated in patients with G4C2 mutation compared to patients without the particular mutation and healthy controls, suggesting eIF4H as a disease modifier (32). In a study investigating spatial memory performance and associated transcriptomic alternations in hippocampus of mice underwent chronically unpredictable stress, eIF4H was among identified genes (33). On one hand, while there had been no direct association of HYI and eIF4H with MDD, our findings suggest further investigations would be warranted.
On the other hand, RBM8a had been associated with brain development and animal model of depression. RBM8a is one of the four core proteins of exon junction complex (EJC), along with Magoh, eIF4A3 (DDX48), and CASC3 (MLN51) (34). Microdeletions of 1q21.1 chromosome, which include Rbm8a gene, were associated with abnormal brain developments such as microcephaly or macrocephaly (35) as well as moderate mental retardation (36). In an animal study, over-expression of RBM8a in the dentate gyrus of the hippocampus of mice induced abnormal behaviors, in anxiety and depression models such as open field test, social interaction test, and forced swimming test (37). In terms of protein-protein interaction analysis using BioGRID database (27), RBM8a interacts with eIF4A3, which is linked with eIF4A1, which is further connected with eIF4H, suggesting the involvement of eukaryotic translation initiation factor family. Our data suggest RBM8a and eIF4H could be involved in the paroxetine treatment response in MDD patients.
Three other proteins identified through iTRAQ [alpha cardiac muscle 1 actin (ACTC1), adenylate kinase isoenzyme 1 (AK1), and mitochondrial dynamin-like 120 kDa protein (OPA1)] were not verified by the Western blot validations, but they could be linked to MDD in interesting ways.
Mutations of alpha cardiac muscle 1 actin (ACTC1) had been associated with hypertrophic cardiomyopathy (HCM) (38). Studies had shown that patients with HCM were more depressed (39, 40). Furthermore, ACTC1 upregulation could be induced by over-expression of dysbindin (41), a protein found to be reduced in hippocampus in patients with schizophrenia (42, 43). Those data suggested ACTC1 could be affected in psychiatric disorders.
Adenylate kinase (AK) plays a crucial role in energetic metabolism (44). Its isoenzyme AK1 had been found to associate with neurodegenerative diseases. Upregulation of AK1 in the frontal cortex was associated with stages 5-6 of Parkinson's disease (45). Immunization of recombinant AK1 in mice could induce Th1-response as an immune protection (46). MDD had been associated with an activation of the immune/inflammatory system, including changes in serum acute phase protein (47–50) and cytokine levels (51–55). Our previous data also showed that cytokine levels differed between antidepressant responders and non-responders (56). The role of AK1 in MDD treatment response should be followed up in the future.
Mitochondrial dynamin-like 120 kDa protein (OPA1) is a GTPase that regulates mitochondrial fusion and cristae structure in the inner mitochondrial membrane (57). Mutation of OPA1 gene has long been associated with dominant optic atrophy (58), but neurodegenerative symptoms had also been reported, as such parkinsonism and dementia (59), multiple sclerosis-like disorder (60), and fatal infantile mitochondrial encephalomyopathy (61). Disturbances in mitochondrial structure and function could impair neuroplasticity (62), while we did not find any report directly linking OPA1 to MDD at the moment, but a follow-up with more clinical samples could be interesting.
Proteomics had been used in identifying novel biomarkers for MDD, mostly comparing the proteomic profiles of patients with MDD with healthy controls (63). Proteomic studies focusing on antidepressant response were far more scarce, let alone ones focusing on one particular drug (64). Utilization of iTRAQ was also relatively novel. In animal model of depression, iTRAQ analysis identified 241 differentially expressed proteins in hippocampi after Chinese herbal medicine xiaochaihutang (XCHT) treatment, and protein-protein interaction network analysis predicted cell division cycle and apoptosis regulator protein 1 (Ccar1) and Calretinin (Calb2) might play the central roles in XCHT's antidepressant network (65). Our investigation could be the first to use iTRAQ to identify differentially expressed proteins in paroxetine response in patients with MDD.
There were several limitations in our study. The study sample size was small. Samples from six paroxetine responsive patients with MDD were analyzed in the iTRAQ analysis, and only 10 were analyzed in the Western blot validation. Ideally, an independent large cohort should be used to validate the results obtained from iTRAQ analysis, and further analyses such as receiver operating characteristic (ROC) analysis and principal component analysis (PCA) should be performed. Secondly, while significant changes were shown in both iTRAQ and Western blot validation, the directions of changes were not always consistent, which may indicate the possibility of false positive results. Thirdly, only one iTRAQ experiment was performed to compare data before and after paroxetine treatment without experimental, biological, or technical replicate. At least one biological replcate was highly recommended in iTRAQ experiments to minimize variations (66), though by combining our samples we hoped to reduce random biological variations. Fourthly, p-values calculation should be performed prior to candidate selection, and proteins identified and quantified by a single peptide are often discarded from further analysis. The readers are warned against over-interpret our study results because of those limitations, and a larger sample size with more controlled variables will be needed before a firm conclusion could be made.
iTRAQ could rapidly identify and quantify differentially expressed proteins, despite being costly. In patients with MDD showing response to paroxetine treatment, putative hydroxypyruvate isomerase (HYI), eukaryotic translation initiation factor 4H (eIF4H), and RNA binding motif 8A (RBM8A) were differentially expressed before and after treatment, as identified by iTRAQ and validated by Western blot. Those identified proteins could be potential predictors of paroxetine treatment response in patients with MDD. This finding could help establish future individualized medicine.
Data Availability Statement
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (24) partner repository with the dataset identifier PXD028947 and 10.6019/PXD028947.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Review Board (IRB) of Chang Gung Memorial Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
T-LH and JS contributed substantially to conception and design. T-LH and HS contributed to acqusition of data. HS, C-CL, and T-LH helped with analysis and interpretation of data. C-CL drafted the article and revised it critically for important intellectual content. T-LH approved the final draft. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by a clinical research grants CMRPG8G1161 and CMRPG8J1511 from Kaohsiung Chang Gung Memorial Hospital, Taiwan and Ministry of Science and Technology, Taiwan (NMRPG8H0101).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The mass spectrometry-based proteomics analysis was done by the Clinical Proteomics Core Laboratory of Chang Gung Memorial Hospital at Linkou, Taiwan.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.577857/full#supplementary-material
References
1. Bakish D. New standard of depression treatment: remission and full recovery. J Clin Psychiatry. (2001) 62(Suppl 26):5–9.
2. Belmaker RH, Agam G. Major depressive disorder. N Engl J Med. (2008) 358:55–68. doi: 10.1056/NEJMra073096
3. Rao C, Shi H, Zhou C, Zhu D, Zhao M, Wang Z, et al. Hypothalamic proteomic analysis reveals dysregulation of glutamate balance and energy metabolism in a mouse model of chronic mild stress-induced depression. Neurochem Res. (2016) 41:2443–56. doi: 10.1007/s11064-016-1957-2
4. Cheng K, Li J, Yang D, Yang Y, Rao C, Zhang S, et al. 2D-gel based proteomics unravels neurogenesis and energetic metabolism dysfunction of the olfactory bulb in CUMS rat model. Behav Brain Res. (2016) 313:302–9. doi: 10.1016/j.bbr.2016.05.026
5. Li J, Zhang SX, Wang W, Cheng K, Guo H, Rao CL, et al. Potential antidepressant and resilience mechanism revealed by metabolomic study on peripheral blood mononuclear cells of stress resilient rats. Behav Brain Res. (2017) 320:12–20. doi: 10.1016/j.bbr.2016.11.035
6. Louveau A, Smirnov I, Keyes TJ, Eccles JD, Rouhani SJ, Peske JD, et al. Structural and functional features of central nervous system lymphatic vessels. Nature. (2015) 523:337–41. doi: 10.1038/nature14432
7. Miller AH, Raison CL. The role of inflammation in depression: from evolutionary imperative to modern treatment target. Nat Rev Immunol. (2016) 16:22–34. doi: 10.1038/nri.2015.5
8. Wolkowitz OM, Mellon SH, Epel ES, Lin J, Reus VI, Rosser R, et al. Resting leukocyte telomerase activity is elevated in major depression and predicts treatment response. Mol Psychiatry. (2012) 17:164–72. doi: 10.1038/mp.2010.133
9. Belzeaux R, Formisano-Treziny C, Loundou A, Boyer L, Gabert J, Samuelian JC, et al. Clinical variations modulate patterns of gene expression and define blood biomarkers in major depression. J Psychiatr Res. (2010) 44:1205–13. doi: 10.1016/j.jpsychires.2010.04.011
10. Huang TL, Cho YT, Su H, Shiea J. Principle component analysis combined with matrix-assisted laser desorption ionization mass spectrometry for rapid diagnosing the sera of patients with major depression. Clin Chim Acta. (2013) 424:175–81. doi: 10.1016/j.cca.2013.06.003
11. Lo LH, Huang TL, Shiea J. Acid hydrolysis followed by matrix-assisted laser desorption/ionization mass spectrometry for the rapid diagnosis of serum protein biomarkers in patients with major depression. Rapid Commun Mass Spectrom. (2009) 23:589–98. doi: 10.1002/rcm.3908
12. Lo LH, Shiea J, Huang TL. Rapid detection of alteration of serum IgG in patients with schizophrenia after risperidone treatment by matrix-assisted laser desorption ionization/time-of-flight mass spectrometry. Rapid Commun Mass Spectrom. (2016) 30:2645–9. doi: 10.1002/rcm.7753
13. Huang TL, Sung ML, Chen TY. 2D-DIGE proteome analysis on the platelet proteins of patients with major depression. Proteome Sci. (2014) 12:1. doi: 10.1186/1477-5956-12-1
14. Wang Q, Su X, Jiang X, Dong X, Fan Y, Zhang J, et al. iTRAQ technology-based identification of human peripheral serum proteins associated with depression. Neuroscience. (2016) 330:291–325. doi: 10.1016/j.neuroscience.2016.05.055
15. Han X, Shao W, Liu Z, Fan S, Yu J, Chen J, et al. iTRAQ-based quantitative analysis of hippocampal postsynaptic density-associated proteins in a rat chronic mild stress model of depression. Neuroscience. (2015) 298:220–92. doi: 10.1016/j.neuroscience.2015.04.006
16. Palmfeldt J, Henningsen K, Eriksen SA, Muller HK, Wiborg O. Protein biomarkers of susceptibility and resilience to stress in a rat model of depression. Mol Cell Neurosci. (2016) 74:87–95. doi: 10.1016/j.mcn.2016.04.001
17. Fleurence R, Williamson R, Jing Y, Kim E, Tran QV, Pikalov AS, et al. A systematic review of augmentation strategies for patients with major depressive disorder. Psychopharmacol Bull. (2009) 42:57–90.
18. Bousman CA, Forbes M, Jayaram M, Eyre H, Reynolds CF, Berk M, et al. Antidepressant prescribing in the precision medicine era: a prescriber's primer on pharmacogenetic tools. BMC Psychiatry. (2017) 17:60. doi: 10.1186/s12888-017-1230-5
19. Hamilton M. Development of a rating scale for primary depressive illness. Br J Soc Clin Psychol. (1967) 6:278–96. doi: 10.1111/j.2044-8260.1967.tb00530.x
20. Beck AT, Kovacs M, Weissman A. Assessment of suicidal intention: the Scale for Suicide Ideation. J Consult Clin Psychol. (1979) 47:343–52. doi: 10.1037/0022-006X.47.2.343
21. Butler AM, Malone K. Attempted suicide v. non-suicidal self-injury: behaviour, syndrome or diagnosis? Br J Psychiatry. (2013) 202:324–5. doi: 10.1192/bjp.bp.112.113506
22. Kageyama Y, Kasahara T, Morishita H, Mataga N, Deguchi Y, Tani M, et al. Search for plasma biomarkers in drug-free patients with bipolar disorder and schizophrenia using metabolome analysis. Psychiatry Clin Neurosci. (2017) 71:115–23. doi: 10.1111/pcn.12461
23. Xu HB, Zhang RF, Luo D, Zhou Y, Wang Y, Fang L, et al. Comparative proteomic analysis of plasma from major depressive patients: identification of proteins associated with lipid metabolism and immunoregulation. Int J Neuropsychopharmacol. (2012) 15:1413–25. doi: 10.1017/S1461145712000302
24. Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. (2019) 47:D442–50. doi: 10.1093/nar/gky1106
25. Haugse R, Langer A, Gullaksen SE, Sundoy SM, Gjertsen BT, Kotopoulis S, et al. Intracellular signaling in key pathways is induced by treatment with ultrasound and microbubbles in a leukemia cell line, but not in healthy peripheral blood mononuclear cells. Pharmaceutics. (2019) 11:319. doi: 10.3390/pharmaceutics11070319
26. Gallo-Oller G, Ordonez R, Dotor J. A new background subtraction method for Western blot densitometry band quantification through image analysis software. J Immunol Methods. (2018) 457:1–5. doi: 10.1016/j.jim.2018.03.004
27. Oughtred R, Rust J, Chang C, Breitkreutz BJ, Stark C, Willems A, et al. The BioGRID database: a comprehensive biomedical resource of curated protein, genetic, and chemical interactions. Protein Sci. (2021) 30:187–200. doi: 10.1002/pro.3978
28. Pham TK, Roy S, Noirel J, Douglas I, Wright PC, Stafford GP. A quantitative proteomic analysis of biofilm adaptation by the periodontal pathogen Tannerella forsythia. Proteomics. (2010) 10:3130–41. doi: 10.1002/pmic.200900448
29. Saei AA, Sabatier P, Tokat UG, Chernobrovkin A, Pirmoradian M, Zubarev RA. Comparative proteomics of dying and surviving cancer cells improves the identification of drug targets and sheds light on cell life/death decisions. Mol Cell Proteomics. (2018) 17:1144–55. doi: 10.1074/mcp.RA118.000610
30. Bai Y, Lu C, Zhang G, Hou Y, Guo Y, Zhou H, et al. Overexpression of miR-519d in lung adenocarcinoma inhibits cell proliferation and invasion via the association of eIF4H. Tumour Biol. (2017) 39:1010428317694566. doi: 10.1177/1010428317694566
31. Vaysse C, Philippe C, Martineau Y, Quelen C, Hieblot C, Renaud C, et al. Key contribution of eIF4H-mediated translational control in tumor promotion. Oncotarget. (2015) 6:39924–40. doi: 10.18632/oncotarget.5442
32. Goodman LD, Prudencio M, Srinivasan AR, Rifai OM, Lee VM, Petrucelli L, et al. eIF4B and eIF4H mediate GR production from expanded G4C2 in a Drosophila model for C9orf72-associated ALS. Acta Neuropathol Commun. (2019) 7:62. doi: 10.1186/s40478-019-0711-9
33. Jung SH, Brownlow ML, Pellegrini M, Jankord R. Divergence in morris water maze-based cognitive performance under chronic stress is associated with the hippocampal whole transcriptomic modification in mice. Front Mol Neurosci. (2017) 10:275. doi: 10.3389/fnmol.2017.00275
34. Bartkowska K, Tepper B, Turlejski K, Djavadian RL. Roles of the exon junction complex components in the central nervous system: a mini review. Rev Neurosci. (2018) 29:817–24. doi: 10.1515/revneuro-2017-0113
35. Brunetti-Pierri N, Berg JS, Scaglia F, Belmont J, Bacino CA, Sahoo T, et al. Recurrent reciprocal 1q21.1 deletions and duplications associated with microcephaly or macrocephaly and developmental and behavioral abnormalities. Nat Genet. (2008) 40:1466–71. doi: 10.1038/ng.279
36. Gamba BF, Zechi-Ceide RM, Kokitsu-Nakata NM, Vendramini-Pittoli S, Rosenberg C, Krepischi Santos AC, et al. Interstitial 1q21.1 microdeletion is associated with severe skeletal anomalies, dysmorphic face and moderate intellectual disability. Mol Syndromol. (2016) 7:344–8. doi: 10.1159/000450971
37. Alachkar A, Jiang D, Harrison M, Zhou Y, Chen G, Mao Y. An EJC factor RBM8a regulates anxiety behaviors. Curr Mol Med. (2013) 13:887–99. doi: 10.2174/15665240113139990019
38. Despond EA, Dawson JF. Classifying cardiac actin mutations associated with hypertrophic cardiomyopathy. Front Physiol. (2018) 9:405. doi: 10.3389/fphys.2018.00405
39. Igoumenou A, Alevizopoulos G, Anastasakis A, Stavrakaki E, Toutouzas P, Stefanadis C. Depression in patients with hypertrophic cardiomyopathy: is there any relation with the risk factors for sudden death? Heart Asia. (2012) 4:44–8. doi: 10.1136/heartasia-2012-010099
40. Serber ER, Sears SF, Nielsen CD, Spencer WH 3rd, Smith KM. Depression, anxiety, and quality of life in patients with obstructive hypertrophic cardiomyopathy three months after alcohol septal ablation. Am J Cardiol. (2007) 100:1592–7. doi: 10.1016/j.amjcard.2007.06.063
41. Rangrez AY, Bernt A, Poyanmehr R, Harazin V, Boomgaarden I, Kuhn C, et al. Dysbindin is a potent inducer of RhoA-SRF-mediated cardiomyocyte hypertrophy. J Cell Biol. (2013) 203:643–56. doi: 10.1083/jcb.201303052
42. Weickert CS, Rothmond DA, Hyde TM, Kleinman JE, Straub RE. Reduced DTNBP1 (dysbindin-1) mRNA in the hippocampal formation of schizophrenia patients. Schizophr Res. (2008) 98:105–10. doi: 10.1016/j.schres.2007.05.041
43. Talbot K, Eidem WL, Tinsley CL, Benson MA, Thompson EW, Smith RJ, et al. Dysbindin-1 is reduced in intrinsic, glutamatergic terminals of the hippocampal formation in schizophrenia. J Clin Invest. (2004) 113:1353–63. doi: 10.1172/JCI200420425
44. Ionescu MI. Adenylate Kinase: A ubiquitous enzyme correlated with medical conditions. Protein J. (2019) 38:120–33. doi: 10.1007/s10930-019-09811-0
45. Garcia-Esparcia P, Hernandez-Ortega K, Ansoleaga B, Carmona M, Ferrer I. Purine metabolism gene deregulation in Parkinson's disease. Neuropathol Appl Neurobiol. (2015) 41:926–40. doi: 10.1111/nan.12221
46. Gao Y, Zhou X, Wang H, Liu R, Ye Q, Zhao Q, et al. Immunization with recombinant schistosome adenylate kinase 1 partially protects mice against Schistosoma japonicum infection. Parasitol Res. (2017) 116:1665–74. doi: 10.1007/s00436-017-5441-y
47. Huang T-L. Lower serum albumin levels in patients with mood disorders. Chang Gung Med J. (2002) 25:509–13.
48. Maes M, Delange J, Ranjan R, Meltzer HY, Desnyder R, Cooremans W, et al. Acute phase proteins in schizophrenia, mania and major depression: modulation by psychotropic drugs. Psychiatry Res. (1997) 66:1–11. doi: 10.1016/S0165-1781(96)02915-0
49. Maes M, Wauters A, Neels H, Scharpe S, Van Gastel A, D'Hondt P, et al. Total serum protein and serum protein fractions in depression: relationships to depressive symptoms and glucocorticoid activity. J Affect Disord. (1995) 34:61–9. doi: 10.1016/0165-0327(94)00106-J
50. Song C, Dinan T, Leonard BE. Changes in immunoglobulin, complement and acute phase protein levels in the depressed patients and normal controls. J Affect Disord. (1994) 30:283–8. doi: 10.1016/0165-0327(94)90135-X
51. Huang T-L, Chen J-F. Cholesterol and lipids in depression: stress, hypothalamo-pituitary-adrenocortical axis, and inflammation/immunity. Adv Clin Chem. (2005) 39:81–8105. doi: 10.1016/S0065-2423(04)39003-7
52. Maes M, Bosmans E, Suy E, Vandervorst C, DeJonckheere C, Raus J. Depression-related disturbances in mitogen-induced lymphocyte responses and interleukin-1 beta and soluble interleukin-2 receptor production. Acta Psychiatr Scand. (1991) 84:379–86. doi: 10.1111/j.1600-0447.1991.tb03163.x
53. Seidel A, Arolt V, Hunstiger M, Rink L, Behnisch A, Kirchner H. Cytokine production and serum proteins in depression. Scand J Immunol. (1995) 41:534–8. doi: 10.1111/j.1365-3083.1995.tb03604.x
54. Sluzewska A, Rybakowski J, Bosmans E, Sobieska M, Berghmans R, Maes M, et al. Indicators of immune activation in major depression. Psychiatry Res. (1996) 64:161–7. doi: 10.1016/S0165-1781(96)02783-7
55. Hirano T, Akira S, Taga T, Kishimoto T. Biological and clinical aspects of interleukin 6. Immunol Today. (1990) 11:443–9. doi: 10.1016/0167-5699(90)90173-7
56. Lin CC, Lee CT, Huang TL. Cytokines, treatment responses, and body weight changes in major depressive patients treated with antidepressants. Taiwan J Psychiatry. (2014) 28:95–102.
57. Del Dotto V, Fogazza M, Lenaers G, Rugolo M, Carelli V, Zanna C. OPA1: how much do we know to approach therapy? Pharmacol Res. (2018) 131:199–210. doi: 10.1016/j.phrs.2018.02.018
58. Hudson G, Amati-Bonneau P, Blakely EL, Stewart JD, He L, Schaefer AM, et al. Mutation of OPA1 causes dominant optic atrophy with external ophthalmoplegia, ataxia, deafness and multiple mitochondrial DNA deletions: a novel disorder of mtDNA maintenance. Brain. (2008) 131:329–37. doi: 10.1093/brain/awm272
59. Carelli V, Musumeci O, Caporali L, Zanna C, La Morgia C, Del Dotto V, et al. Syndromic parkinsonism and dementia associated with OPA1 missense mutations. Ann Neurol. (2015) 78:21–38. doi: 10.1002/ana.24410
60. Verny C, Loiseau D, Scherer C, Lejeune P, Chevrollier A, Gueguen N, et al. Multiple sclerosis-like disorder in OPA1-related autosomal dominant optic atrophy. Neurology. (2008) 70:1152–3. doi: 10.1212/01.wnl.0000289194.89359.a1
61. Spiegel R, Saada A, Flannery PJ, Burte F, Soiferman D, Khayat M, et al. Fatal infantile mitochondrial encephalomyopathy, hypertrophic cardiomyopathy and optic atrophy associated with a homozygous OPA1 mutation. J Med Genet. (2016) 53:127–31. doi: 10.1136/jmedgenet-2015-103361
62. Cheng A, Hou Y, Mattson MP. Mitochondria and neuroplasticity. ASN Neuro. (2010) 2:e00045. doi: 10.1042/AN20100019
63. Woods AG, Wormwood KL, Iosifescu DV, Murrough J, Darie CC. Protein biomarkers in major depressive disorder: an update. Adv Exp Med Biol. (2019) 1140:585–600. doi: 10.1007/978-3-030-15950-4_35
64. Gadad BS, Jha MK, Czysz A, Furman JL, Mayes TL, Emslie MP, et al. Peripheral biomarkers of major depression and antidepressant treatment response: current knowledge and future outlooks. J Affect Disord. (2018) 233:3–14. doi: 10.1016/j.jad.2017.07.001
65. Zhang K, He M, Su D, Pan X, Li Y, Zhang H, et al. Quantitative proteomics reveal antidepressant potential protein targets of xiaochaihutang in corticosterone induced model of depression. J Ethnopharmacol. (2019) 231:438–45. doi: 10.1016/j.jep.2018.11.020
Keywords: eukaryotic translation initiation factor 4H (eIF4H), iTRAQ, major depressive disorder, putative hydroxypyruvate isomerase (HYI), RNA binding motif 8A (RBM8A)
Citation: Lin C-C, Su H, Shiea J and Huang T-L (2022) Isobaric Tags for Relative and Absolute Quantitation Identification of Blood Proteins Relevant to Paroxetine Response in Patients With Major Depressive Disorder. Front. Psychiatry 13:577857. doi: 10.3389/fpsyt.2022.577857
Received: 30 June 2020; Accepted: 24 March 2022;
Published: 18 April 2022.
Edited by:
Gianluca Serafini, San Martino Hospital (IRCCS), ItalyReviewed by:
Caroline Evans, The University of Sheffield, United KingdomKuanjun He, Inner Mongolia University for Nationalities, China
Copyright © 2022 Lin, Su, Shiea and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tiao-Lai Huang, a540520@adm.cgmh.org.tw