Identification of molecular mechanism and key biomarkers in membranous nephropathy by bioinformatics analysis

Objectives: Membranous nephropathy (MN) is an autoimmune nephropathy. The incidence of MN is increasing gradually in recent years. Previous studies focused on antibody production, complement activation and podocyte injury in MN. However, the etiology and underlying mechanism of MN remain to be further studied. Methods: GSE104948 and GSE108109 of glomerular expression profile were downloaded from Gene Expression Omnibus (GEO) database, GSE47184, GSE99325, GSE104954, GSE108112, GSE133288 of renal tubule expression profile, and GSE73953 of peripheral blood mononuclear cells (PBMCs) expression profile. After data integration by Network analyst, differentially expressed genes (DEGs) between MN and healthy samples were obtained. DEGs were enriched in gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) networks of these genes were constructed through Metascape, etc. We further understood the function of hub genes through gene set enrichment analysis (GSEA). The diagnostic value of DEGs in MN was evaluated by receiver operating characteristic (ROC) analysis. Results: A total of 3 genes (TP53, HDAC5, and SLC2A3) were screened out. Among them, the up-regulated TP53 expression may be closely related to MN renal pathological changes. However, the expression of MN podocyte target antigen was not significantly different from that of healthy controls. In addition, the changes of Wnt signaling pathway in PBMCs and the effects of SLC2A3 on the differentiation of M2 monocyte need further study. Conclusion: It is difficult to unify a specific mechanism for the changes of glomerulus, renal tubules and PBMCs in MN patients. This may be related to the pathogenesis, pathology and immune characteristics of MN. MN podocyte target antigen may not be the root cause of the disease, but a stage result in the pathogenesis process.