Identification of key genes in late-onset major depressive disorder through a co-expression network module

Late-onset major depressive disorder (LOD) increases the risk of disability and suicide in elderly patients. However, the complex pathological mechanism of LOD still remains unclear. We selected 10 LOD patients and 12 healthy control samples from the GSE76826 dataset for statistical analysis. Under the screening criteria, 811 differentially expressed genes (DEGs) were screened. We obtained a total of two most clinically significant modules through the weighted gene co-expression network analysis (WGCNA). Functional analysis of the genes in the most clinically significant modules was performed to explore the potential mechanism of LOD, followed by protein-protein interaction (PPI) analysis and hub gene identification in the core area of the PPI network. Furthermore, we identified immune infiltrating cells using the cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm between healthy subjects and LOD patients with the GSE98793 dataset. Next, six hub genes (CD27, IL7R, CXCL1, CCR7, IGLL5, and CD79A) were obtained by intersecting hub genes with DEGs, followed by verifying the diagnostic accuracy with the receiver operating characteristic curve (ROC). In addition, we constructed the least absolute shrinkage and selection operator (LASSO) regression model for hub gene cross-validation. Finally, we found that CD27 and IGLL5 were good diagnostic indicators of LOD, and CD27 may be the key gene of immune function change in LOD. In conclusion, our research shows that the changes in the immune function may be an important mechanism in the development of LOD, which can provide some guidance for the related research of LOD in the future.