Abstract:Objective: To identify the gene coexpression network and hub genes of 3 brain regions including NAc, PAG and mPFC in chronic pain induced mice model of depression, weighted gene coexpression network analysis (WGCNA) was processed. Methods: The RNA-seq dataset of GSE91396 was downloaded from the gene expression omnibus (GEO) database. WGCNA was used to construct the gene coexpression network, and identify the key modules associated with NAc, PAG and mPFC respectively. Gene ontology and pathway enrichment analyses were performed on the key module genes. Then, Cytoscape software was used to screen the hub genes in key modules. Results: There are 9 coexpression modules identified by WGCNA of which the brown, turquoise and yellow modules were the key module related to NAc, PAG and mPFC respectively. Module genes were mainly involved in synapse structure, membrane potential, neurotransmitter and so on. The hub genes with high connectivity were identified in 3 key modules, such as Drd1a, Pdyn and Ppp1r1b in brown module, Cd99l2 and Pcyt1a in turquoise module, and Neurod6, Dlgap1 and Adcy1 in yellow module. Conclusion: WGCNA can identify key modules and hub genes which are biologically relevant to the 3 depression associated brain regions. This study will provide a novel insight into the understanding of the depression induced by chronic pain.