Investigation of molecular biomarker candidates for diagnosis and prognosis of chronic periodontitis by bioinformatics analysis of pooled microarray gene expression datasets in Gene Expression Omnibus (GEO).

Background: Chronic periodontitis (CP) is a multifactorial inflammatory disease. For the diagnosis of CP, it is necessary to investigate molecular biomarkers and the biological pathway of CP. Although analysis of mRNA expression profiling with microarray is useful to elucidate pathological mechanisms of multifactorial diseases, it is expensive. Therefore, we utilized pooled microarray gene expression data on the basis of data sharing to reduce hybridization costs and compensate for insufficient mRNA sampling. The aim of the present study was to identify molecular biomarker candidates and biological pathways of CP using pooled datasets in the Gene Expression Omnibus (GEO) database.

Methods: Three pooled transcriptomic datasets (GSE10334, GSE16134, and GSE23586) of gingival tissue with CP in the GEO database were analyzed for differentially expressed genes (DEGs) using GEO2R, functional analysis and biological pathways with the Database of Annotation Visualization and Integrated Discovery database, Protein-Protein Interaction (PPI) network and hub gene with the Search Tool for the Retrieval of Interaction Genes database, and biomarker candidates for diagnosis and prognosis and upstream regulators of dominant biomarker candidates with the Ingenuity Pathway Analysis database.

Results: We shared pooled microarray datasets in the GEO database. One hundred and twenty-three common DEGs were found in gingival tissue with CP, including 81 upregulated genes and 42 downregulated genes. Upregulated genes in Gene Ontology were significantly enriched in immune responses, and those in the Kyoto Encyclopedia of Genes and Genomes pathway were significantly enriched in the cytokine-cytokine receptor interaction pathway, cell adhesion molecules, and hematopoietic cell lineage. From the PPI network, the 12 nodes with the highest degree were screened as hub genes. Additionally, six biomarker candidates for CP diagnosis and prognosis were screened.

Conclusions: We identified several potential biomarkers for CP diagnosis and prognosis (e.g., CSF3, CXCL12, IL1B, MS4A1, PECAM1, and TAGLN) and upstream regulators of biomarker candidates for CP diagnosis (TNF and TGF2). We also confirmed key genes of CP pathogenesis such as CD19, IL8, CD79A, FCGR3B, SELL, CSF3, IL1B, FCGR2B, CXCL12, C3, CD53, and IL10RA. To our knowledge, this is the first report to reveal associations of CD53, CD79A, MS4A1, PECAM1, and TAGLN with CP.