Gene Correlation Network Analysis to Identify Biomarkers of Peri-Implantitis

Background and Objectives: The histopathological and clinical conditions for transforming peri-implant mucositis into peri-implantitis (PI) are not fully clarified. We aim to uncover molecular mechanisms and new potential biomarkers of PI. Materials and Methods: Raw GSE33774 and GSE57631 datasets were obtained from the Gene Expression Omnibus (GEO) database. The linear models for microarray data (LIMMA) package in R software completes differentially expressed genes (DEGs). We conducted a weighted gene co-expression network analysis (WGCNA) on the top 25% of altered genes and identified the key modules associated with the clinical features of PI. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the R software. We constructed a protein-protein interaction (PPI) network through the STRING database. After that we used Cytohubba plug-ins of Cytoscape to screen out the potential hub genes, which were subsequently verified via receiver operating characteristic (ROC) curves in another dataset, GSE178351, and revalidation of genes through the DisGeNET database. Results: We discovered 632 DEGs (570 upregulated genes and 62 downregulated genes). A total of eight modules were screened by WGCNA, among which the turquoise module was most correlated with PI. The Cytohubba plug-ins were used for filtering hub genes, which are highly linked with PI development, from the candidate genes in the protein-protein interaction (PPI) network. Conclusions: We found five key genes from PI using WGCNA. Among them, ICAM1, CXCL1, and JUN are worthy of further study of new target genes, providing the theoretical basis for further exploration of the occurrence and development mechanism of PI.