Identification of hub genes and transcription factors involved in periodontitis on the basis of multiple microarray analysis.

Objectives: To identify the differentially expressed genes (DEGs) during the pathogenesis of periodontitis by bioinformatics analysis.

Methods: GEO2R was used to screen DEGs in GSE10334 and GSE16134. Then, the overlapped DEGs were used for further analysis. g:Profiler was used to perform Gene Ontology analysis and pathway analysis for upregulated and downregulated DEGs. The STRING database was used to construct the protein-protein interaction (PPI) network, which was further visua-lized and analyzed by Cytoscape software. Hub genes and key modules were identified by cytoHubba and MCODE plug-ins, respectively. Finally, transcription factors were predicted via iRegulon plug-in.

Results: A total of 196 DEGs were identified, including 139 upregulated and 57 downregulated DEGs. Functional enrichment analysis showed that the upregulated DEGs were mainly enriched in immune-related pathways including immune system, viral protein interaction with cytokine and cytokine receptor, cytokine-cytokine receptor interaction, leukocyte transendothelial migration, and chemokine receptors bind chemokines. On the contrary, the downregulated DEGs were mainly related to the formation of the cornified envelope and keratinization. The identified hub genes in the PPI network were CXCL8, CXCL1, CXCR4, SEL, CD19, and IKZF1. The top three modules were involved in chemokine response, B cell receptor signaling pathway, and interleukin response, respectively. iRegulon analysis revealed that IRF4 scored the highest.

Conclusions: The pathogenesis of periodontitis was closely associated with the expression levels of the identified hub genes including CXCL8, CXCL1, CXCR4, SELL, CD19, and IKZF1. IRF4, the predicted transcription factor, might serve as a dominant upstream regulator.