OBJECTIVE: We aim to explore the detailed molecular mechanisms of membrane nephropathy (MN) related genes by bioinformatics analysis.METHODS: Two microarray datasets (GSE108109 and GSE104948) with glomerular gene expression data from 65 MN patients and 9 healthy donors were obtained from the Gene Expression Omnibus (GEO) database. After processing the raw data, DEGs screening was conducted using the LIMMA (linear model for microarray data) package and Gene set enrichment analysis (GSEA) was performed with GSEA software (v. 3.0), followed by gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. The protein-protein interaction (PPI) network analysis was carried out to determine the hub genes, by applying the maximal clique centrality (MCC) method, which was visualized by Cytoscape. Finally, utilizing the Nephroseq v5 online platform, we analyzed subgroups associated with hub genes. The findings were further validated by immunohistochemistry (IHC) staining in renal tissues from MN or control patients.RESULTS: A sum of 370 DEGs (188 up-regulated genes, 182 down-regulated genes) and 20 hub genes were ascertained. GO and KEGG enrichment analysis demonstrated that DEGs of MN were preponderantly associated with cell damage and complement cascade-related immune responses. Combined with literature data and hub gene-related MN subset analysis, CTSS, ITGB2, and HCK may play important roles in the pathological process of MN.CONCLUSION: This study identified novel hub genes in MN using bioinformatics. We found that some hub genes such as CTSS, ITGB2, and HCK might contribute to MN immunopathological process, providing new insights for further study of the molecular mechanisms underlying glomerular injury of MN.