Multi-omics integration identifies NK cell dysregulation and a five-gene diagnostic signature in major depressive disorder

Jia Wang¹Ye Kuang¹Chuanmei Peng¹Yong Yuan¹Yong Ji¹Sulian Chen¹Jinrong Tian¹Yuanyuan Zhou¹Xingying Chen¹Jing Li¹Lei Feng^1*Shengjie Nie^2*

• ¹Yan'an Hospital Affiliated To Kunming Medical University, Kunming, China
• ²Kunming Medical University, Kunming, China

Background: Major Depressive Disorder (MDD), a leading global disability affecting 280 million people, has poor treatment efficacy due to persistent biological variability involving cell-type-specific transcriptomic dysregulation and immune dysfunction, and integrated multi-omics approaches are vital to uncover pathways and therapeutic targets. Methods: This research utilized a comprehensive multi-omics approach, merging bulk RNA sequencing (RNA-seq) data from the GSE39653 dataset with single-cell RNA sequencing (scRNA-seq) data derived from peripheral blood mononuclear cells (PBMCs) of three MDD patients and three healthy controls. Analysis of differential gene expression (DEGs1) and identification of genes inside Weighted Gene Co-expression Network Analysis (WGCNA) modules were conducted using bulk RNA-seq data. Analysis of differential cell population abundance and DEGs2 was performed on the scRNA-seq data. Detection of CD3⁻CD56⁺ or CD3⁻CD16⁺ NK cells in human peripheral blood samples by flow cytometry. Candidate genes were subsequently identified from the intersection of DEGs1, WGCNA module genes, and DEGs2. Subsequently, machine learning methods were employed to discern key genes from these candidates. The functional characterisation of essential cell populations was accomplished via pseudotime trajectory analysis, Gene Set Variation Analysis (GSVA), metabolic pathway analysis (scMetabolism), and transcription factor inference (SCENIC). Ultimately, diagnostic models, regulatory networks, and compound screenings were developed based on the key genes. Results: In the RNA-seq analysis, 803 DEGs1 and 2080 WGCNA module genes were identified. scRNA-seq analysis revealed 1,539 DEGs2 and identified natural killer (NK) cells as a major dysfunctional immune cell subpopulation in MDD, exhibiting a significantly increased proportion (CD3⁻CD56⁺ or CD3⁻CD16⁺, p < 0.05) in the MDD patients. The intersection of DEGs1, WGCNA module genes, and DEGs2 yielded 26 candidate genes. Subsequent machine learning analysis identified five key genes: CSPP1, ZNF84, HLA-DPA1, CCZ1, and LRRC8D. A diagnostic nomogram constructed using these key genes demonstrated robust discriminatory performance in distinguishing MDD patients. Mechanistic investigations implicated these five key genes in MDD pathogenesis through neurodegenerative signaling pathways. Conclusion: Our study establishes NK cell dysfunction as a core pathophysiological mechanism in MDD. The identified key genes serve as robust diagnostic biomarkers and therapeutic targets. Elucidation of their regulatory networks provides critical insights for precision psychiatry interventions.