Identification of potential biomarkers and pathways related to major depressive disorder by integrated bioinformatic analysis and experimental validation

Abstract

Objective: To identify promising biomarkers for the pathogenesis of major depressive disorder (MDD).

Methods: Microarray chips of MDD patients, including the GSE98793, GSE52790, and GSE39653 datasets, were obtained from the Gene Expression Omnibus database. The biological processes and pathways related to MDD were investigated using the GO and KEGG pathway tools. Weighted gene coexpression network analysis was conducted to identify modules related to MDD. The hub genes associated with MDD were obtained via protein-protein interaction analysis. Finally, the expression of hub genes in the hippocampal tissues of depression-like rats was detected by reverse transcription-polymerase chain reaction and Western blotting.

Results: A total of 658 differentially expressed genes were identified from the Gene Expression Omnibus datasets; thus, these genes and the GSE98793 dataset were used to conduct weighted gene coexpression network analysis. A total of 244 module-related genes were identified and these genes were highly correlated with MDD. These genes were involved in the Ras signaling pathway, regulation of the actin cytoskeleton, and axon guidance according to the KEGG analysis. Hub genes, including MAPK14, SOCS1, TLR2, PTK2B, and GRB2, were obtained via protein-protein interaction analysis. All these hub genes showed better diagnostic efficiency in the GSE52790, GSE39653, and GSE98793 datasets. In vivo experiments revealed that compared with those in control rats, SOCS1 and MAPK14 expression was significantly decreased; while GRB2, TLR2, and PTK2B expression was increased in the hippocampi of depression-like rats.

Conclusions: Our study demonstrates that GRB2, TLR2, SOCS1, PTK2B, and MAPK14 are promising hub genes, and targeting these five genes may be an effective treatment strategy for MDD.

Keywords

Major depressive disorder / Bioinformatic / Biomarkers / Microarray / Hub genes / Weighted gene coexpression network analysis

Conflict of interest statement

The authors declare that there is no conflict of interest.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81774281 and No.82474303), the Natural Science Foundation of Hunan Province (2023JJ30888), and the leading national joint discipline of Chinese and Western medicines to the Chinese Medicine Department, Xiangya Hospital, CSU.

Data availability statement

The data supporting the findings of this study are available from the corresponding author upon request.

Authors’ contributions

YZ performed the experiments and wrote the paper. MZ, LQP, RZ, and KCN contributed to data analysis and interpretation, and provided experimental reagents and materials. WH was responsible for designing the experimental framework and supervising the entire project. All authors reviewed and approved the final manuscript.

Publisher’s note

The Publisher of the Journal remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

[1] Kennis M, Gerritsen L, van Dalen M, Williams A, Cuijpers P, Bockting C. Prospective biomarkers of major depressive disorder: A systematic review and meta-analysis. Mol Psychiatry 2020;25:321-338.

[2] Rehm J, Shield KD. Global burden of disease and the impact of mental and addictive disorders. Curr Psychiatry Rep 2019;21:10.

[3] Almeida OP, Ford AH, Hankey GJ, Golledge J, Yeap BB, Flicker L. Depression, antidepressants and the risk of cardiovascular events and death in older men. Maturitas 2019;128:4-9.

[4] Liu L, Wang H, Chen X, Zhang Y, Zhang H, Xie P. Gut microbiota and its metabolites in depression: From pathogenesis to treatment. Ebiomedicine 2023;90:104527.

[5] Zhang G, Xu S, Yuan Z, Shen L. Weighted gene coexpression network analysis identifies specific modules and hub genes related to major depression. Neuropsychiatr Dis Treat 2020;16:703-713.

[6] Geschwind DH, Flint J. Genetics and genomics of psychiatric disease. Science 2015;349:1489-1494.

[7] Okazawa H, Ikawa M, Tsujikawa T, Kiyono Y, Yoneda M. Brain imaging for oxidative stress and mitochondrial dysfunction in neurodegenerative diseases. Q J Nucl Med Mol Imag 2014;58:387-397.

[8] Tobore TO. On the neurobiological role of oxidative stress in alcohol-induced impulsive, aggressive and suicidal behavior. Subst Use Misuse 2019;54:2290-2303.

[9] Hodes GE, Kana V, Menard C, Merad M, Russo SJ. Neuroimmune mechanisms of depression. Nat Neurosci 2015;18:1386-1393.

[10] Myint AM, Schwarz MJ, Steinbusch HW, Leonard BE. Neuropsychiatric disorders related to interferon and interleukins treatment. Metab Brain Dis 2009;24:55-68.

[11] Deng Z, Cai W, Liu J, Deng A, Yang Y, Tu J, et al. Co-expression modules construction by WGCNA and identify potential hub genes and regulation pathways of postpartum depression. Front Biosci (Landmark Ed) 2021;26:1019-1030.

[12] Zhou Z, Wang Y, Sun S, Zhang K, Wang L, Zhao H, et al. Paeonia lactiflora pall. polysaccharide alleviates depression in CUMS mice by inhibiting the NLRP3/ASC/Caspase-1 signaling pathway and affecting the composition of their intestinal flora. J Ethnopharmacol 2023;316:116716.

[13] Zuo X, Chen Z, Cai J, Gao W, Zhang Y, Han G, et al. 5-Hydroxytryptamine receptor 1d aggravates hepatocellular carcinoma progression through FoxO6 in AKT-dependent and independent manners. Hepatology 2019;69:2031-2047.

[14] Rotenstein LS, Ramos MA, Torre M, Segal JB, Peluso MJ, Guille C, et al. Prevalence of depression, depressive symptoms, and suicidal ideation among medical students: A systematic review and meta-analysis. JAMA 2016;316:2214-2236.

[15] Silverman-Gavrila LB, Charlton MP. Calcineurin and cytoskeleton in low-frequency depression. J Neurochem 2009;109:716-732.

[16] Myrum C, Soule J, Bittins M, Cavagnini K, Goff K, Ziemek SK, et al. Arc interacts with the integral endoplasmic reticulum protein, Calnexin. Front Cell Neurosci 2017;11:294.

[17] Li S, Zhu Z, Lan T, Wu Y, Li Y, Wang C, et al. Levomilnacipran ameliorates lipopolysaccharide-induced depression-like behaviors and suppressed the TLR4/Ras signaling pathway. Int Immunopharmacol 2023;122:110595.

[18] Li H, Linjuan-Li, Wang Y. G-CSF improves CUMS-induced depressive behaviors through downregulating Ras/ERK/MAPK signaling pathway. Biochem Biophys Res Commun 2016;479:827-832.

[19] Stoeckli ET. Understanding axon guidance: Are we nearly there yet? Development 2018;145.

[20] He J, Xie J, Zhou G, Jia C, Han D, Li D, et al. Active fraction of Polyrhachis vicina roger (AFPR) ameliorate depression induced inflammation response by FTO/miR-221-3p/SOCS1 axis. J Inflamm Res 2023;16:6329-6348.

[21] Sun Z, Zhan H, Wang C, Guo P. Shanzhiside methylester protects against depression by inhibiting inflammation via the miRNA-155-5p/SOCS1 axis. Psychopharmacology (Berl) 2022;239:2201-2213.

[22] Liu Y, Yu J, Wang X, Dong J. MicroRNA-345-5p regulates depression by targeting suppressor of cytokine signaling 1. Brain Behav 2020;10:e1653.

[23] Liu Q, Sun NN, Wu ZZ, Fan DH, Cao MQ. Chaihu-Shugan-San exerts an antidepressive effect by downregulating miR-124 and releasing inhibition of the MAPK14 and Gria3 signaling pathways. Neural Regen Res 2018;13:837-845.

[24] Su W, Ju J, Gu M, Wang X, Liu S, Yu J, et al. SARS-CoV-2 envelope protein triggers depression-like behaviors and dysosmia via TLR2-mediated neuroinflammation in mice. J Neuroinflammation 2023;20:110.

[25] Medina-Rodriguez EM, Cheng Y, Michalek SM, Beurel E, Jope RS. Toll-like receptor 2 (TLR2)-deficiency impairs male mouse recovery from a depression-like state. Brain Behav Immun 2020;89:51-58.

[26] Amiresmaeili A, Roohollahi S, Mostafavi A, Askari N. Effects of oregano essential oil on brain TLR4 and TLR2 gene expression and depressive-like behavior in a rat model. Res Pharm Sci 2018;13:130-141.