# Shared Molecular Mechanisms between Alzheimer’s Disease and Periodontitis Revealed by Transcriptomic Analysis

Authors: Jieqi Jin; Mengkai Guang; Anthony Chukwunonso Ogbuehi; Simin Li; Kai Zhang; Yihong Ma; Aneesha Acharya; Bihan Guo; Zongwu Peng; Xiangqiong Liu; Yupei Deng; Zhaobi Fang; Xiongjie Zhu; Shiting Hua; Cong Li; Rainer Haak; Dirk Ziebolz; Gerhard Schmalz; Lei Liu; Baohua Xu; Xiaofeng Huang
Journal: BioMed Research International
 (2021)

---

## Abstract

Objective. To investigate the genetic crosstalk mechanisms that link periodontitis and Alzheimer’s disease (AD). Background. Periodontitis, a common oral infectious disease, is associated with Alzheimer’s disease (AD) and considered a putative contributory factor to its progression. However, a comprehensive investigation of potential shared genetic mechanisms between these diseases has not yet been reported. Methods. Gene expression datasets related to periodontitis were downloaded from the Gene Expression Omnibus (GEO) database, and differential expression analysis was performed to identify differentially expressed genes (DEGs). Genes associated with AD were downloaded from the DisGeNET database. Overlapping genes among the DEGs in periodontitis and the AD-related genes were defined as crosstalk genes between periodontitis and AD. The Boruta algorithm was applied to perform feature selection from these crosstalk genes, and representative crosstalk genes were thus obtained. In addition, a support vector machine (SVM) model was constructed by using the scikit-learn algorithm in Python. Next, the crosstalk gene-TF network and crosstalk gene-DEP (differentially expressed pathway) network were each constructed. As a final step, shared genes among the crosstalk genes and periodontitis-related genes in DisGeNET were identified and denoted as the core crosstalk genes. Results. Four datasets (GSE23586, GSE16134, GSE10334, and GSE79705) pertaining to periodontitis were included in the analysis. A total of 48 representative crosstalk genes were identified by using the Boruta algorithm. Three TFs (FOS, MEF2C, and USF2) and several pathways (i.e., JAK-STAT, MAPK, NF-kappa B, and natural killer cell-mediated cytotoxicity) were identified as regulators of these crosstalk genes. Among these 48 crosstalk genes and the chronic periodontitis-related genes in DisGeNET, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3 were shared and identified as the most pivotal candidate links between periodontitis and AD. Conclusions. Exploration of available transcriptomic datasets revealed C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3 as the top candidate molecular linkage genes between periodontitis and AD.


---

## Body

## 1. Introduction

An association between periodontitis and Alzheimer’s disease (AD) has been demonstrated, and periodontitis reportedly confers risk for the incidence and progression of Alzheimer’s disease (AD) [1, 2]. Alzheimer’s disease is a neurodegenerative disease characterized by the formation of amyloid-β peptide (AβP) plaques and intraneuronal neurofibrillary tangles (NFTs), which drives neuroinflammation in the brain [3]. Periodontitis, a chronic, immunoinflammatory disease affecting supporting structures of teeth, is multifactorial in nature, driven by polymicrobial dysbiosis and unfavorable shifts in the plaque biofilm composition, which disrupts the host-microbial homeostasis [4, 5]. Inflammation is considered the key connecting link between both these diseases [6]. Two purported mechanistic links have been highlighted. First, proinflammatory mediators including specific cytokines or chemokines in the periodontal milieu that enter systemic circulation impose a systemic inflammatory burden, propagating inflammatory responses by microglial cells in the brain [6]. Second, periodontal pathogens may directly enter the brain via blood circulation or peripheral nerves, as evidenced by the discovery of the keystone periodontal pathogen Porphyromonas gingivalis (Pg) in AD patients’ brains [7, 8].At the same time, genetic susceptibility and gene dysregulation have been identified in the context of both AD and periodontitis. In AD, frequently reported dysregulated genes include amyloid-beta precursor protein (APP), presenilin 1 (PSEN1), and presenilin 2 (PSEN2), apolipoproteins, and lipid homeostasis, genes involved in endocytosis, and membrane-spanning 4 (MS4) family [9]. In periodontitis, aberrant genes highlighted include interleukin-1, interleukin-6, interleukin-10, transforming growth factor-beta (TGF-β), tumor necrosis factor-α (TNF-α), interferon-gamma (IFN-γ), and matrix metalloproteinases (MMPs) among others [10]. It is plausible that gene dysregulation in periodontitis could contribute to its association with AD, and crosstalk genes may biologically link AD and periodontitis serving as either shared susceptibility factors or molecular links.Here, we designed a bioinformatic study of existing experimental datasets to understand putative molecular links between periodontitis and Alzheimer’s disease by identifying crosstalk genes, transcription factors, and signaling pathways involved in both disorders. The molecular mechanisms identified through this approach could suggest potential therapeutic targets particularly relevant to drug development and personalized medicine approaches.


## 2. Materials and Methods

### 2.1. Study Design

Figure1 depicts a flowchart outlining the study workflow. DEGs dysregulated in periodontitis- and AD-related genes were obtained from the GEO database and DisGeNET database, respectively. Crosstalk genes linking periodontitis and AD were identified as the AD-related genes that overlapped with significantly up- or downregulated DEGs in periodontitis. Thereafter, feature selection from the crosstalk genes was performed using a conventional recursive feature elimination (RFE) algorithm and the Boruta algorithm. The crosstalk genes obtained by feature selection were used to construct two networks to identify the transcription factors and the differentially expressed pathways that target these crosstalk genes. In the next step, “core” crosstalk genes were identified as the crosstalk genes obtained by feature selection that were overlapping with chronic periodontitis-related genes in the DisGeNET database.Figure 1
Flowchart depicting study workflow.


### 2.2. Procurement of Periodontitis-Related Datasets

Sample-matched whole-genome gene expression datasets from periodontitis were sourced and downloaded from the NCBI Gene Expression Omnibus (GEO). The eligibility criteria for these datasets were as follows: datasets that included established periodontitis samples as the experimental group and healthy gingival samples as the control group, where periodontitis was defined based in accordance with the case definition presented in the 2017 World Workshop: (1) interdental CAL detectable at ≥2 nonadjacent teeth or (2) buccal or oralCAL≥3mm with pocketing >3 mm detectable at ≥2 teeth [11].


### 2.3. Differential Gene Expression Analysis

Differential gene expression analysis of periodontitis-related datasets was carried out using the Linear Models for Microarray (limma) package [12] in the R project (version 3.0.1, http://www.r-project.org/) [13]. Three such datasets, GSE23586, GSE16134, and GSE10334, were sourced and analyzed. Genes with p value < 0.05 and ∣logFCfoldchange∣≥1 were regarded as significant differentially expressed genes (DEGs). For another dataset, GSE79705, the screening range of DEGs was broadened by extending the thresholds and settings to p value < 0.05 and ∣logFC∣>0 as DEGs.Next, a Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/) was drawn to identify shared genes within the DEGs identified from the four datasets. The common up/downregulated DEGs in four datasets were used for the following analyses, and DEGs that were not common to all the datasets were excluded.


### 2.4. Functional Enrichment Analysis

Functional enrichment analysis of up/downregulated periodontitis-related DEGs was performed with DAVID (Database for Annotation, Visualization and Integrated Discovery, v6.8) [14, 15]. p<0.05 was set as the threshold. GO (Gene Ontology) and pathway enrichment analysis of the identified DEGs was performed [16, 17].


### 2.5. Construction of the Protein-Protein Interaction Network of Periodontitis-Related DEGs

280,826 PPI (protein-protein interaction) pairs including 19,610 genes were downloaded from HPRD (Human Protein Reference Database) [18], BioGRID (Biological General Repository for Interaction Datasets) [19], DIP (Database of Interacting Proteins) [20], MINT (Molecular INTeraction database) [21], PINA (Protein Interaction Network Analysis) [22], InnateDB (a knowledge resource for innate immunity interactions and pathways) [23], and INstruct (3D protein interactome networks with structural resolution) [24]. The Cytoscape platform was used to visualize the network derived from PPI pairs of DEGs and conduct a topological analysis [25].


### 2.6. Identification of Crosstalk Genes and Construction of the Crosstalk Gene-Related PPI Network

AD-related genes were downloaded from the DisGeNET database [26]. The AD-related genes that overlapped with the up- and downregulated periodontitis-related DEGs were identified. These overlapping genes were regarded as “crosstalk” genes linking AD and periodontitis and further used for constructing a crosstalk gene-related PPI network. The crosstalk gene-related PPI network consisted of four types of nodes, namely, (1) DEGs dysregulated in periodontitis (not related to AD), (2) crosstalk genes or AD-related genes which were also DEGs dysregulated in periodontitis, (3) AD-related genes (not dysregulated in periodontitis), and (4) other genes (neither related to AD nor dysregulated in periodontitis).


### 2.7. Feature Selection from Crosstalk Genes

Since GSE16134 had the largest sample size among the periodontitis-related datasets, it was used as the test set. The other three datasets (GSE23586, GSE10334, and GSE79705) were used as validation sets. Firstly, expression values of the crosstalk genes (identified in the previous step) from GSE16134 were used as input for the Boruta algorithm in the R project [27] and the conventional recursive feature elimination (RFE) algorithm [28] and feature selection was performed. Each gene was regarded as a feature.


### 2.8. Support Vector Machine (SVM) Modeling Using Feature-Selected Crosstalk Genes

The expression values in GSE16134 and GSE10334 were scale-standardized. Next, it was examined if the crosstalk genes obtained by feature selection were found in the four periodontitis-related datasets (i.e., GSE23586, GSE16134, GSE10334, and GSE79705). The gene expression values of these feature selection-obtained crosstalk genes were extracted from these datasets. If the number of expression profile genes in a certain dataset after extraction was lower than the number of feature selection-obtained crosstalk genes, the expression values of the missing genes were considered missing values and represented by the NA symbol. The missing values were processed by using the DMwR package [29] in R, and the K-Nearest Neighbors (KNN) algorithm was used to impute these missing values (NA) in that dataset. By imputing these missing values, all the four periodontitis-related datasets presented all the feature selection-obtained crosstalk genes. Thereafter, the scikit-learn package [30] was used to perform a grid search, and the best hyperparameters of a support vector machine (SVM) model were found by using 5-fold cross-validation (CV) [31]. A SVM classifier model was established by using data from GSE16134 as the training set and test set, where the samples of the GSE16134 dataset were split into 60% : 40% for the training set and test set, respectively. Data from the other three datasets (GSE23586, GSE10334, and GSE79705) were used as the validation set. The decision function method was used to obtain the score for each sample. Next, receiver operating characteristic (ROC) curves for the four datasets were generated by using the pROC package and displayed using the ggplot2 package in R.


### 2.9. Targeting Relationships between Transcription Factors (TFs) and Crosstalk Genes

Transcription factor- (TF-) target gene regulation pairs were obtained and downloaded from TRRUST [32], cGRNB [33], HTRIdb [34], ORTI [35], and TRANSFAC [36] databases. The TF-target gene interaction pairs corresponding to the feature selection-obtained crosstalk genes were extracted and used for constructing a TF-target gene interaction network with visualization using Cytoscape software [37].


### 2.10. Pathway Analysis of the Crosstalk Genes

Human data describing relationships between signaling pathways and genes were downloaded from the KEGG database, and all pathways related to the feature selection-obtained crosstalk genes were extracted. The expression levels of these crosstalk gene-related pathways were plotted as a heatmap, and differential expression analysis was performed and applied to the four periodontitis-related datasets to identify the DEPs (differentially expressed pathways) using the R package limma. For the three datasets (GSE23586, GSE16134, and GSE79705), the pathways withp value < 0.05 and ∣logFC∣≥1 were regarded as differentially expressed pathways (DEPs), while for the dataset GSE79705, the pathways with p value < 0.05 and ∣logFC∣>0 were regarded as DEPs.


### 2.11. Classification Performance of the Core Crosstalk Genes

The genes related to periodontitis were downloaded from the DisGeNET database [26]. The overlap between the periodontitis-related genes obtained from the DisGeNET database and the feature selection-obtained crosstalk genes was analyzed, and the overlapping genes were termed core crosstalk genes. The corresponding expression values of these overlapping genes in the four periodontitis-related datasets were obtained, and ROC curves were drawn.


## 2.1. Study Design

Figure1 depicts a flowchart outlining the study workflow. DEGs dysregulated in periodontitis- and AD-related genes were obtained from the GEO database and DisGeNET database, respectively. Crosstalk genes linking periodontitis and AD were identified as the AD-related genes that overlapped with significantly up- or downregulated DEGs in periodontitis. Thereafter, feature selection from the crosstalk genes was performed using a conventional recursive feature elimination (RFE) algorithm and the Boruta algorithm. The crosstalk genes obtained by feature selection were used to construct two networks to identify the transcription factors and the differentially expressed pathways that target these crosstalk genes. In the next step, “core” crosstalk genes were identified as the crosstalk genes obtained by feature selection that were overlapping with chronic periodontitis-related genes in the DisGeNET database.Figure 1
Flowchart depicting study workflow.


## 2.2. Procurement of Periodontitis-Related Datasets

Sample-matched whole-genome gene expression datasets from periodontitis were sourced and downloaded from the NCBI Gene Expression Omnibus (GEO). The eligibility criteria for these datasets were as follows: datasets that included established periodontitis samples as the experimental group and healthy gingival samples as the control group, where periodontitis was defined based in accordance with the case definition presented in the 2017 World Workshop: (1) interdental CAL detectable at ≥2 nonadjacent teeth or (2) buccal or oralCAL≥3mm with pocketing >3 mm detectable at ≥2 teeth [11].


## 2.3. Differential Gene Expression Analysis

Differential gene expression analysis of periodontitis-related datasets was carried out using the Linear Models for Microarray (limma) package [12] in the R project (version 3.0.1, http://www.r-project.org/) [13]. Three such datasets, GSE23586, GSE16134, and GSE10334, were sourced and analyzed. Genes with p value < 0.05 and ∣logFCfoldchange∣≥1 were regarded as significant differentially expressed genes (DEGs). For another dataset, GSE79705, the screening range of DEGs was broadened by extending the thresholds and settings to p value < 0.05 and ∣logFC∣>0 as DEGs.Next, a Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/) was drawn to identify shared genes within the DEGs identified from the four datasets. The common up/downregulated DEGs in four datasets were used for the following analyses, and DEGs that were not common to all the datasets were excluded.


## 2.4. Functional Enrichment Analysis

Functional enrichment analysis of up/downregulated periodontitis-related DEGs was performed with DAVID (Database for Annotation, Visualization and Integrated Discovery, v6.8) [14, 15]. p<0.05 was set as the threshold. GO (Gene Ontology) and pathway enrichment analysis of the identified DEGs was performed [16, 17].


## 2.5. Construction of the Protein-Protein Interaction Network of Periodontitis-Related DEGs

280,826 PPI (protein-protein interaction) pairs including 19,610 genes were downloaded from HPRD (Human Protein Reference Database) [18], BioGRID (Biological General Repository for Interaction Datasets) [19], DIP (Database of Interacting Proteins) [20], MINT (Molecular INTeraction database) [21], PINA (Protein Interaction Network Analysis) [22], InnateDB (a knowledge resource for innate immunity interactions and pathways) [23], and INstruct (3D protein interactome networks with structural resolution) [24]. The Cytoscape platform was used to visualize the network derived from PPI pairs of DEGs and conduct a topological analysis [25].


## 2.6. Identification of Crosstalk Genes and Construction of the Crosstalk Gene-Related PPI Network

AD-related genes were downloaded from the DisGeNET database [26]. The AD-related genes that overlapped with the up- and downregulated periodontitis-related DEGs were identified. These overlapping genes were regarded as “crosstalk” genes linking AD and periodontitis and further used for constructing a crosstalk gene-related PPI network. The crosstalk gene-related PPI network consisted of four types of nodes, namely, (1) DEGs dysregulated in periodontitis (not related to AD), (2) crosstalk genes or AD-related genes which were also DEGs dysregulated in periodontitis, (3) AD-related genes (not dysregulated in periodontitis), and (4) other genes (neither related to AD nor dysregulated in periodontitis).


## 2.7. Feature Selection from Crosstalk Genes

Since GSE16134 had the largest sample size among the periodontitis-related datasets, it was used as the test set. The other three datasets (GSE23586, GSE10334, and GSE79705) were used as validation sets. Firstly, expression values of the crosstalk genes (identified in the previous step) from GSE16134 were used as input for the Boruta algorithm in the R project [27] and the conventional recursive feature elimination (RFE) algorithm [28] and feature selection was performed. Each gene was regarded as a feature.


## 2.8. Support Vector Machine (SVM) Modeling Using Feature-Selected Crosstalk Genes

The expression values in GSE16134 and GSE10334 were scale-standardized. Next, it was examined if the crosstalk genes obtained by feature selection were found in the four periodontitis-related datasets (i.e., GSE23586, GSE16134, GSE10334, and GSE79705). The gene expression values of these feature selection-obtained crosstalk genes were extracted from these datasets. If the number of expression profile genes in a certain dataset after extraction was lower than the number of feature selection-obtained crosstalk genes, the expression values of the missing genes were considered missing values and represented by the NA symbol. The missing values were processed by using the DMwR package [29] in R, and the K-Nearest Neighbors (KNN) algorithm was used to impute these missing values (NA) in that dataset. By imputing these missing values, all the four periodontitis-related datasets presented all the feature selection-obtained crosstalk genes. Thereafter, the scikit-learn package [30] was used to perform a grid search, and the best hyperparameters of a support vector machine (SVM) model were found by using 5-fold cross-validation (CV) [31]. A SVM classifier model was established by using data from GSE16134 as the training set and test set, where the samples of the GSE16134 dataset were split into 60% : 40% for the training set and test set, respectively. Data from the other three datasets (GSE23586, GSE10334, and GSE79705) were used as the validation set. The decision function method was used to obtain the score for each sample. Next, receiver operating characteristic (ROC) curves for the four datasets were generated by using the pROC package and displayed using the ggplot2 package in R.


## 2.9. Targeting Relationships between Transcription Factors (TFs) and Crosstalk Genes

Transcription factor- (TF-) target gene regulation pairs were obtained and downloaded from TRRUST [32], cGRNB [33], HTRIdb [34], ORTI [35], and TRANSFAC [36] databases. The TF-target gene interaction pairs corresponding to the feature selection-obtained crosstalk genes were extracted and used for constructing a TF-target gene interaction network with visualization using Cytoscape software [37].


## 2.10. Pathway Analysis of the Crosstalk Genes

Human data describing relationships between signaling pathways and genes were downloaded from the KEGG database, and all pathways related to the feature selection-obtained crosstalk genes were extracted. The expression levels of these crosstalk gene-related pathways were plotted as a heatmap, and differential expression analysis was performed and applied to the four periodontitis-related datasets to identify the DEPs (differentially expressed pathways) using the R package limma. For the three datasets (GSE23586, GSE16134, and GSE79705), the pathways withp value < 0.05 and ∣logFC∣≥1 were regarded as differentially expressed pathways (DEPs), while for the dataset GSE79705, the pathways with p value < 0.05 and ∣logFC∣>0 were regarded as DEPs.


## 2.11. Classification Performance of the Core Crosstalk Genes

The genes related to periodontitis were downloaded from the DisGeNET database [26]. The overlap between the periodontitis-related genes obtained from the DisGeNET database and the feature selection-obtained crosstalk genes was analyzed, and the overlapping genes were termed core crosstalk genes. The corresponding expression values of these overlapping genes in the four periodontitis-related datasets were obtained, and ROC curves were drawn.


## 3. Results

### 3.1. Included Periodontitis-Related Datasets

Four datasets pertaining to periodontitis (i.e., GSE23586, GSE16134, GSE10334, and GSE79705) were included and analyzed. Table1 provides key details regarding the included datasets.Table 1
Details of the included periodontitis-related GEO datasets.
Included four datasetsType of periodontitisExperimental platformNumber of examined genesNumber of inflamed gingival tissue samplesNumber of healthy control samplesNumber of total samplesGSE235863 patients with severe chronic periodontitisGPL57023,518336GSE16134120 patients (65 with chronic periodontitis and 55 with aggressive periodontitis)GPL57024,44124169310GSE1033463 with chronic periodontitis and 27 with aggressive periodontitisGPL57024,44118364247GSE79705Generalized aggressive periodontitis (GAgP):n=4; chronic periodontitis (CP): n=4GPL1873419,3058412


### 3.2. Identification of Periodontitis-Related DEGs and Their Functions

Table2 shows the number of up- and downregulated DEGs that were identified in each of the four datasets. Figure 2 depicts a Venn diagram showing the overlap of DEGs identified from the four datasets. A total of 816 DEGs including 405 upregulated genes and 411 downregulated genes were used for the following analyses. Figure 3 shows the biological processes and signaling pathways in which the up- and downregulated DEGs were significantly enriched.Table 2
The number of up/downregulated DEGs identified in the four periodontitis-related datasets GSE23586, GSE16134, GSE10334, and GSE79705 at thresholds set for defining DEGs.
DataNumber of upregulated DEGsNumber of downregulated DEGsNumber of total DEGs∣log2FC∣p valueExpression scale (if the data normalization was performed)GSE2358676976>1<0.05TrueGSE1613418848236>1<0.05FalseGSE1033415242194>1<0.05FalseGSE79705219291510>0<0.05TrueFigure 2
Venn diagram showing the overlap between DEGs identified in the four periodontitis-related datasets.Figure 3
Functional enrichment analysis of up- and downregulated periodontitis-related DEGs. (a) The functional terms (the biological processes (A) and KEGG pathways (B)) enriched by upregulated periodontitis-related DEGs. (b) The functional terms (the biological processes (A) and KEGG pathways (B)) enriched by downregulated periodontitis-related DEGs.
(a)(b)


### 3.3. The Hub Genes Identified by the Periodontitis-Related PPI Network

Figure4 depicts a PPI network based on the periodontitis-related DEGs, and Table 3 shows the topological characteristics of the top 30 nodes in the PPI network. As seen in Table 3, several genes with the highest degree were identified as hub genes. These upregulated DEGs included SMAD3, TRIM27, VIM, YWHAH, and FOS, and downregulated DEGs included MYC, HSPB1, DDB1, RPS3, KAT5, SMARCA4, RPL13, and PLCG1.Figure 4
Protein-protein interaction (PPI) network associated with periodontitis. Red nodes represent the DEGs which were upregulated in periodontitis, green nodes represent the DEGs which were downregulated in periodontitis, and blue nodes represent other genes which were not DEGs but interacted with DEGs.Table 3
Topological characteristics of the top 30 nodes in the periodontitis-related PPI network, ranked in descending order of degree.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientMYCDEG_down10272.4327270.1772170.4110610.0018960.002805HSPB1DEG_down4332.6380610.0599080.3790660.002780.005213SMAD3DEG_up3962.6676360.0521670.3748640.0042830.005425UBC3952.1120.2416730.4734850.0045240.007004TRIM27DEG_up3342.8647270.0467920.3490730.0007910.006306VIMDEG_up3202.7052120.0376070.3696570.0052510.007619YWHAHDEG_up2822.8387880.0341150.3522630.0018680.00865DDB1DEG_down2762.764970.0296770.3616680.0035310.008012RPS3DEG_down2752.8448480.0150360.3515130.0193230.014779KAT5DEG_down2452.7176970.0282240.3679590.007160.007732SMARCA4DEG_down2422.7157580.0216460.3682210.0101160.008636RPL13DEG_down2342.7551520.0157730.3629560.0245410.011238PLCG1DEG_down2272.8850910.0300110.346610.0021440.009614FOSDEG_up2262.8807270.0245380.3471350.0032650.009974AP2M1DEG_down2142.8501820.0255820.3508550.0076350.010848SF3B3DEG_down2032.7655760.0151880.3615880.0145830.010661TP53RKDEG_up2002.8713940.0139720.3482630.016080.012798MMS19DEG_down1952.8069090.0223350.3562640.0047050.008763RPS7DEG_down1932.883030.0081260.3468570.0283350.016367MEPCEDEG_down1892.9052120.01450.3442090.0050660.012337TUBB4BDEG_down1732.8803640.0081670.3471780.0075950.019412CDH1DEG_down1712.8884850.0161860.3462020.0058480.011775EZRDEG_down1702.9066670.0162950.3440370.0041070.01554DDX41DEG_down1682.9500610.0095450.3389760.0070570.017509NUDCDEG_up1602.9437580.0168860.3397020.002830.013954SF3A2DEG_down1562.9498180.0090750.3390040.0099260.01593RPL27DEG_down1472.9127270.0039490.3433210.0439850.020057NIFKDEG_up1453.1078790.0111110.3217630.0103450.016053SRPRBDEG_up1432.9453330.0175450.339520.0016740.015983ATP1A1DEG_up1402.8166060.0112690.3550370.0091470.013032


### 3.4. Crosstalk Genes Bridging Alzheimer’s Disease and Periodontitis

In total, 51 upregulated crosstalk genes and 41 downregulated crosstalk genes were identified and are listed in Table4. The crosstalk gene-related PPI network as shown in Figure 5 consisted of 3496 nodes and 5141 edges. The topological characteristics of the top 30 nodes in this network are presented in Table 5. Hub crosstalk genes with the highest degree included MYC, HSPB1, VIM, KAT5, RPL13, FOS, and CDH1.Table 4
The 51 crosstalk genes that were upregulated and 41 crosstalk genes which were downregulated in periodontitis.
Regulation pattern in periodontitisCrosstalk genes linking periodontitis and Alzheimer’s diseaseUpregulatedAMIGO1 (adhesion molecule with Ig-like domain 1; gene ID: 57463)ARG2 (arginase 2; gene ID: 384)BDKRB2 (bradykinin receptor B2; gene ID: 624)BRI3 (brain protein I3; gene ID: 25798)C1D (C1D nuclear receptor corepressor; gene ID: 10438)C3 (complement C3; gene ID: 718)C4A (complement C4A (Rodgers blood group); gene ID: 720)C4B (complement C4B; gene ID: 721)C4B_2 (complement component 4B (Chido blood group), copy 2; gene ID: 100293534)CASP7 (caspase 7; gene ID: 840)CD177 (CD177 molecule; gene ID: 57126)CD38 (CD38 molecule; gene ID: 952)CHCHD10 (coiled-coil-helix-coiled-coil-helix domain containing 10; gene ID: 400916)COX8A (cytochrome c oxidase subunit 8A; gene ID: 1351)CSF3 (colony-stimulating factor 3; gene ID: 1440)CTGF (cellular communication network factor 2; gene ID: 1490)CTSS (cathepsin S; gene ID: 1520)CXCL1 (C-X-C motif chemokine ligand 1; gene ID: 2919)CXCL12 (C-X-C motif chemokine ligand 12; gene ID: 6387)CXCL8 (C-X-C motif chemokine ligand 8; gene ID: 3576)CXCR4 (C-X-C motif chemokine receptor 4; gene ID: 7842)ENPP2 (ectonucleotide pyrophosphatase/phosphodiesterase 2; gene ID: 5168)FCGR3A (Fc fragment of IgG receptor IIIa; gene ID: 2214)FCGR3B (Fc fragment of IgG receptor IIIb; gene ID: 2215)FHL2 (four and a half LIM domains 2; gene ID: 2274)FOS (Fos proto-oncogene, AP-1 transcription factor subunit; gene ID: 2353)FOXO1 (forkhead box O1; gene ID: 2308)GSTK1 (glutathione S-transferase kappa 1; gene ID: 373156)GSTO1 (glutathione S-transferase omega 1; gene ID: 9446)HCLS1 (hematopoietic cell-specific Lyn substrate 1; gene ID: 3059)IGFBP7 (insulin-like growth factor binding protein 7; gene ID: 3490)IL1B (interleukin-1 beta; gene ID: 3553)LYZ (lysozyme; gene ID: 4069)MAK16 (MAK16 homolog; gene ID: 84549)MEF2C (myocyte-specific enhancer factor 2C; gene ID: 4208)MME (membrane metalloendopeptidase; gene ID: 4311)MMP1 (matrix metallopeptidase 1; gene ID: 4312)MMP3 (matrix metallopeptidase 3; gene ID: 4314)MS4A1 (membrane-spanning 4-domains A1; gene ID: 931)MSRA (methionine sulfoxide reductase A; gene ID: 4482)MZB1 (marginal zone B and B1 cell-specific protein; gene ID: 51237)NDUFB8 (NADH:ubiquinone oxidoreductase subunit B8; gene ID: 4714)PECAM1 (platelet endothelial cell adhesion molecule 1; gene ID: 5175)PLAT (plasminogen activator, tissue type; gene ID: 5327)PSENEN (presenilin enhancer, gamma-secretase subunit; gene ID: 55851)SEL1L (SEL1L adaptor subunit of ERAD E3 ubiquitin ligase; gene ID: 6400)SERPINI1 (serpin family I member 1; gene ID: 5274)ST6GAL1 (ST6 beta-galactoside alpha-2,6-sialyltransferase 1; gene ID: 6480)VCAN (versican; gene ID: 1462)VIM (vimentin; gene ID: 7431)WARS2 (tryptophanyl tRNA synthetase 2, mitochondrial; gene ID: 10352)DownregulatedABCA12 (ATP binding cassette subfamily A member 12; gene ID: 26154)ACE (angiotensin I converting enzyme; gene ID: 1636)ADRB1 (adrenoceptor beta 1; gene ID: 153)BCL2L2 (BCL2-like 2; gene ID: 599)BRSK1 (BR serine/threonine kinase 1; gene ID: 84446)CALML5 (calmodulin-like 5; gene ID: 51806)CALR (calreticulin; gene ID: 811)CD36 (CD36 molecule; gene ID: 948)CDH1 (cadherin 1; gene ID: 999)CDK5R1 (cyclin-dependent kinase 5 regulatory subunit 1; gene ID: 8851)CTDSP2 (CTD small phosphatase 2; gene ID: 10106)DNM1 (dynamin 1; gene ID: 1759)DSG1 (desmoglein 1; gene ID: 1828)FLG (filaggrin; gene ID: 2312)GJA1 (gap junction protein alpha 1; gene ID: 2697)GRHL3 (grainyhead-like transcription factor 3; gene ID: 57822)GRIN3B (glutamate ionotropic receptor NMDA type subunit 3B; gene ID: 116444)HM13 (histocompatibility minor 13; gene ID: 81502)HMGCR (3-hydroxy-3-methylglutaryl-CoA reductase; gene ID: 3156)HSPB1 (heat shock protein family B (small) member 1; gene ID: 3315)KAT5 (lysine acetyltransferase 5; gene ID: 10524)MAP3K12 (mitogen-activated protein kinase kinase kinase 12; gene ID: 7786)MED12 (mediator complex subunit 12; gene ID: 9968)MFN2 (mitofusin 2; gene ID: 9927)MYC (MYC proto-oncogene, bHLH transcription factor; gene ID: 4609)NEFL (neurofilament light; gene ID: 4747)NES (exportin 1; gene ID: 7514)NGF (nerve growth factor; gene ID: 4803)NOS3 (nitric oxide synthase 3; gene ID: 4846)NPTXR (neuronal pentraxin receptor; gene ID: 23467)PLXNA3 (plexin A3; gene ID: 55558)PPIL2 (peptidylprolyl isomerase-like 2; gene ID: 23759)PTGS1 (prostaglandin-endoperoxide synthase 1; gene ID: 5742)RPL13 (ribosomal protein L13; gene ID: 6137)SCIMP (SLP adaptor and CSK interacting membrane protein; gene ID: 388325)SPPL2B (signal peptide peptidase-like 2B; gene ID: 56928)TIAF1 (TGFB1-induced antiapoptotic factor 1; gene ID: 9220)TP73 (tumor protein P73; gene ID: 7161)TYRP1 (tyrosinase-related protein 1; gene ID: 7306)USF2 (upstream transcription factor 2, C-Fos interacting; gene ID: 7392)VEGFA (vascular endothelial growth factor A; gene ID: 7422)Figure 5
The crosstalk gene-related PPI network.Table 5
Topological characteristics of the top 30 nodes in the crosstalk gene-related PPI network, ranked in descending order of degree.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientMYCDEG_AD_down10272.1275250.4666980.470035.03E−040.001842HSPB1DEG_AD_down4332.369590.1845980.4220140.0016360.003351VIMDEG_AD_up3202.6633010.1166110.3754740.0017630.005783KAT5DEG_AD_down2452.5663590.0894760.3896570.0030110.005621RPL13DEG_AD_down2342.5897290.0897080.3861410.0013940.005672FOSDEG_AD_up2262.8831510.0836420.3468431.97E−040.017606CDH1DEG_AD_down1712.931910.0584280.34107500.025911TP73DEG_AD_down1372.614830.0467810.3824340.0063330.009173FHL2DEG_AD_up1362.6410850.0571920.3786320.0040310.00889CALRDEG_AD_down1352.8124640.0461420.355560.0027640.010949GSTK1DEG_AD_up1112.9884590.0351520.33462100.026289MMEDEG_AD_up1042.8369880.0403860.3524870.0013070.011694MED12DEG_AD_down912.6956150.0306820.3709730.0061050.012517CASP7DEG_AD_up902.9010390.0271310.3447040.0052430.015603DNM1DEG_AD_down842.9601850.0290860.3378170.0011470.023716FOXO1DEG_AD_up672.9777840.0174030.335820.0040710.026586CXCR4DEG_AD_up652.9985570.0208510.3334940.0038460.028808NOS3DEG_AD_down622.9636470.0187860.3374220.0026440.027957NDUFB8DEG_AD_up563.0424120.021360.32868700.027447UBCOther gene532.0848240.1492650.4796570.0159650.027336NESDEG_AD_down532.9209460.0170120.3423550.0079830.022954CDK5R1DEG_AD_down503.01760.01330.3313890.0024490.037222NEFLDEG_AD_down492.9134450.0148860.3432360.0085030.024673C3DEG_AD_up433.0432780.017050.32859300.042101SEL1LDEG_AD_up423.0360650.0130660.32937400.045635DSG1DEG_AD_down413.0461630.008570.32828200.068464GJA1DEG_AD_down413.0351990.0157030.32946800.04037MEF2CDEG_AD_up413.0392380.0086760.3290300.068892LYZDEG_AD_up403.0421230.0085780.32871800.065086GSTO1DEG_AD_up382.8202540.0109050.3545780.0199150.030356


### 3.5. Crosstalk Genes Obtained by Feature Selection

A total of 48 crosstalk genes were selected by using the Boruta algorithm (Figure6(a)). In addition, 62 crosstalk genes were selected by using the RFE algorithm (Figure 6(b)). All 48 genes obtained by using the Boruta algorithm were included in the 62 genes obtained by the RFE algorithm, indicating that this 48-gene set was representative of the characteristics of all 92 crosstalk genes.Figure 6
48 crosstalk genes selected by using the Boruta algorithm (a) and 62 crosstalk genes selected by using the RFE algorithm (b).
(a)
Boruta algorithm to obtain genes(b)
RFE algorithm to obtain genes


### 3.6. Classification Accuracy Using the Feature Selection-Obtained Crosstalk Genes

The 48 genes obtained by SVM feature selection were not shown in all of the four periodontitis-related datasets. The gene expression profiles of GSE16134 and GSE10334 included all of the 48 crosstalk genes, whereas GSE23586 included 44 crosstalk genes and GSE79705 included 45 crosstalk genes. The classification performance of these 48 crosstalk genes for the four datasets is shown in Table6. For the test set GSE16134 and the validation set GSE10334, the accuracy performance was high at 91.94% and 88.26%, respectively. By comparison, the performance for the other two datasets GSE23586 and GSE79705 was low at 50% and 66.66%, respectively.Table 6
Classification performance of the 48 feature selection-obtained crosstalk genes in the four periodontitis-related datasets.
GSE16134GSE10334GSE23586GSE79705Samples: test setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracyCase939291.94%18317788.26%3350%8866.66%Control312264413040Total12411424721863128


### 3.7. ROC Curves for the Four Periodontitis-Related Datasets

As shown in Figure7, the AUC (area under the curve) values for the GSE16134 test set and the GSE10334 validation set were high at 95.77% and 90.53%, respectively, congruent with the results in Table 6. It was thus inferred that the classifier performance was adequate only when the sample sizes of the validation sets were similar to those of the training and test datasets, and therefore, poor performance was noted for GSE23586 and GSE79705 having much lower sample numbers.Figure 7
The ROC curves of four periodontitis-related datasets.


### 3.8. The Identification of Transcription Factors Regulating the Crosstalk Genes

As shown in Figure8, the TF-crosstalk gene target network consisted of 388 nodes and 1178 edges. Several transcription factors which were also DEGs played critical roles by regulating the most number of crosstalk genes, for example, FOS, MEF2C, and USF2 (Table 7).Figure 8
The transcription factor-crosstalk gene target network.Table 7
Topological characteristics of the top 30 nodes in the TF-crosstalk gene target network, ranking in descending order.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientFOSDEG_up&TF892.1782950.2061010.4590750.0204290.03825VEGFADEG_down772.2118860.1972190.4521030.0109360.038456MEF2CDEG_up&TF762.3281650.1599450.4295230.0042110.044258CDH1DEG_down692.5400520.1593920.39369300.086634CXCR4DEG_up472.5968990.0751050.3850759.25E−040.083836MMEDEG_up442.6640830.043650.37536400.131313C4ADEG_up382.6614990.0297420.3757280.011380.118617VIMDEG_up382.5219640.0602930.3965160.0128020.065789CSF3DEG_up362.5503880.0317340.3920970.0190480.075137C4BDEG_up342.6795870.0188870.3731920.0124780.122037VCANDEG_up342.7260980.0305670.36682500.143717IL1BDEG_up332.7416020.057460.3647500.134602BRSK1DEG_down302.7416020.0304040.3647500.153788SERPINI1DEG_up292.7932820.0190140.35800200.127586USF2DEG_down&TF272.5917310.0275750.3858420.0455840.081607GRHL3DEG_down272.7751940.0270670.36033500.150732HSPB1DEG_down262.7803620.0179810.35966500.155678MMP3DEG_up262.6072350.0243770.3835480.0369230.074984CTGFDEG_up252.7881140.0206690.35866500.161905HMGCRDEG_down252.7674420.0411190.36134500.16GATA2TF242.0956070.0526870.4771890.0181160.104678PLATDEG_up242.7622740.0280750.36202100.173148ETS1TF222.2739020.037670.4397730.0129870.092022NPTXRDEG_down222.8268730.0204780.35374800.154545ARTF212.240310.0406980.4463670.023810.105465YBX1TF212.3152450.0378010.431920.0047620.098928BCL2L2DEG_down212.8139530.0181790.35537200.144048FCGR3ADEG_up212.811370.0205210.35569900.192799ABCA12DEG_down202.798450.0109610.35734100.2SP1TF202.1912140.0373620.4563680.0263160.109167CDK5R1DEG_down202.8242890.0178780.35407100.156098FLGDEG_down202.6537470.0122810.3768260.0421050.093043


### 3.9. Signaling Pathways Enriched in the Crosstalk Genes

From the 48 feature selection-obtained crosstalk genes, 37 crosstalk genes were found among gene-pathway interaction pair data in the KEGG database. 137 KEGG pathways corresponded to these 37 crosstalk genes. Figure9 shows the expression values of these 137 pathways in the four periodontitis-related datasets.Figure 9
Heatmaps depicting 137 signaling pathways enriched in the crosstalk genes in the four periodontitis-related datasets.
(a)
GSE10334(b)
GSE16134(c)
GSE23586(d)
GSE79705The numbers of crosstalk gene-related DEPs obtained from each of the four periodontitis-related datasets are listed in Table8. The interaction relationships between crosstalk genes and DEPs are depicted in Figure 10, showing that several DEPs were dysregulated in at least two datasets, including cytokine-related pathways (cytokine-cytokine receptor interaction, chemokine, and IL-17), immune cell-related pathways (T cell receptor, B cell receptor, Th1 and Th2 cell differentiation, Th17 cell differentiation, natural killer cell-mediated cytotoxicity, and osteoclast differentiation), JAK-STAT signaling, NOD-like receptor signaling, MAPK signaling, Toll-like receptor signaling, NF-kappa B signaling, and C-type lectin receptor signaling.Table 8
The number of crosstalk gene-related DEPs within the four periodontitis-related datasets.
DataSampleDEPDEPDEPLog FCp valueUpDownTotalAbsGSE23586611213>1<0.05GSE1613431084286>1<0.05GSE1033424777279>1<0.05GSE7970512448>0<0.05Figure 10
The crosstalk gene-differentially expressed pathway interaction network.


### 3.10. The Identification of Core Crosstalk Genes

Among the 48 feature selection-obtained crosstalk genes and periodontitis-related genes in the DisGeNET database, 12 common genes were identified as core crosstalk genes including C3, C4A, C4B, CXCL12, FCGR3A, FCGR3B, HSPB1, IL1B, MME, MMP3, PLAT, and VEGFA (Table9). Among these 12 genes, 6 genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, were found associated with chronic periodontitis.Table 9
The 12 overlapping genes between 48 feature selection-obtained crosstalk genes and periodontitis-related genes in the DisGeNET database. Among the 12 genes, 6 genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, were associated with chronic periodontitis.
Gene symbolGene IDDisease IDDisease nameC3718C0031099PeriodontitisC4A720C0266929Chronic periodontitisC4B721C0266929Chronic periodontitisCXCL126387C0266929Chronic periodontitisFCGR3A2214C0266929Chronic periodontitisFCGR3A2214C0031099PeriodontitisFCGR3A2214C0031106Periodontitis, juvenileFCGR3A2214C0399447Early-onset periodontitisFCGR3A2214C0031030Periapical periodontitisFCGR3B2215C0399447Early-onset periodontitisFCGR3B2215C0031106Periodontitis, juvenileHSPB13315C0031030Periapical periodontitisIL1B3553C1719494Periodontitis, localized aggressiveIL1B3553C1719495Aggressive periodontitis, generalizedIL1B3553C0399447Early-onset periodontitisIL1B3553C0001342Acute periodontitisIL1B3553C0031106Periodontitis, juvenileIL1B3553C0031030Periapical periodontitisIL1B3553C4025886Severe periodontitisIL1B3553C0031099PeriodontitisIL1B3553C0266929Chronic periodontitisMME4311C0031099PeriodontitisMMP34314C0031106Periodontitis, juvenileMMP34314C0001342Acute periodontitisMMP34314C0266929Chronic periodontitisMMP34314C0031099PeriodontitisPLAT5327C1719495Aggressive periodontitis, generalizedVEGFA7422C0031099PeriodontitisFigure11 shows ROC curves for the 6 chronic periodontitis-related genes for each of the four datasets. An AUC value of more than 80% was presented by three genes, C4A, C4B, and CXCL12, for the dataset GSE10334, 5 genes, C4A, C4B, CXCL12, FCGR3A, and IL1B, for GSE16134, 3 genes, CXCL12, IL1B, and MMP3, for GSE23586, and all 6 genes for GSE79705. The 3 genes, C4A, C4B, and CXCL1, had the highest classification accuracy when the datasets with small samples (GSE23586 and GSE79705) were not considered.Figure 11
The ROC curves for 6 chronic periodontitis-related genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, in the four periodontitis-related datasets.
(a)
GSE10334(b)
GSE16134(c)
GSE23586(d)
GSE79705


## 3.1. Included Periodontitis-Related Datasets

Four datasets pertaining to periodontitis (i.e., GSE23586, GSE16134, GSE10334, and GSE79705) were included and analyzed. Table1 provides key details regarding the included datasets.Table 1
Details of the included periodontitis-related GEO datasets.
Included four datasetsType of periodontitisExperimental platformNumber of examined genesNumber of inflamed gingival tissue samplesNumber of healthy control samplesNumber of total samplesGSE235863 patients with severe chronic periodontitisGPL57023,518336GSE16134120 patients (65 with chronic periodontitis and 55 with aggressive periodontitis)GPL57024,44124169310GSE1033463 with chronic periodontitis and 27 with aggressive periodontitisGPL57024,44118364247GSE79705Generalized aggressive periodontitis (GAgP):n=4; chronic periodontitis (CP): n=4GPL1873419,3058412


## 3.2. Identification of Periodontitis-Related DEGs and Their Functions

Table2 shows the number of up- and downregulated DEGs that were identified in each of the four datasets. Figure 2 depicts a Venn diagram showing the overlap of DEGs identified from the four datasets. A total of 816 DEGs including 405 upregulated genes and 411 downregulated genes were used for the following analyses. Figure 3 shows the biological processes and signaling pathways in which the up- and downregulated DEGs were significantly enriched.Table 2
The number of up/downregulated DEGs identified in the four periodontitis-related datasets GSE23586, GSE16134, GSE10334, and GSE79705 at thresholds set for defining DEGs.
DataNumber of upregulated DEGsNumber of downregulated DEGsNumber of total DEGs∣log2FC∣p valueExpression scale (if the data normalization was performed)GSE2358676976>1<0.05TrueGSE1613418848236>1<0.05FalseGSE1033415242194>1<0.05FalseGSE79705219291510>0<0.05TrueFigure 2
Venn diagram showing the overlap between DEGs identified in the four periodontitis-related datasets.Figure 3
Functional enrichment analysis of up- and downregulated periodontitis-related DEGs. (a) The functional terms (the biological processes (A) and KEGG pathways (B)) enriched by upregulated periodontitis-related DEGs. (b) The functional terms (the biological processes (A) and KEGG pathways (B)) enriched by downregulated periodontitis-related DEGs.
(a)(b)


## 3.3. The Hub Genes Identified by the Periodontitis-Related PPI Network

Figure4 depicts a PPI network based on the periodontitis-related DEGs, and Table 3 shows the topological characteristics of the top 30 nodes in the PPI network. As seen in Table 3, several genes with the highest degree were identified as hub genes. These upregulated DEGs included SMAD3, TRIM27, VIM, YWHAH, and FOS, and downregulated DEGs included MYC, HSPB1, DDB1, RPS3, KAT5, SMARCA4, RPL13, and PLCG1.Figure 4
Protein-protein interaction (PPI) network associated with periodontitis. Red nodes represent the DEGs which were upregulated in periodontitis, green nodes represent the DEGs which were downregulated in periodontitis, and blue nodes represent other genes which were not DEGs but interacted with DEGs.Table 3
Topological characteristics of the top 30 nodes in the periodontitis-related PPI network, ranked in descending order of degree.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientMYCDEG_down10272.4327270.1772170.4110610.0018960.002805HSPB1DEG_down4332.6380610.0599080.3790660.002780.005213SMAD3DEG_up3962.6676360.0521670.3748640.0042830.005425UBC3952.1120.2416730.4734850.0045240.007004TRIM27DEG_up3342.8647270.0467920.3490730.0007910.006306VIMDEG_up3202.7052120.0376070.3696570.0052510.007619YWHAHDEG_up2822.8387880.0341150.3522630.0018680.00865DDB1DEG_down2762.764970.0296770.3616680.0035310.008012RPS3DEG_down2752.8448480.0150360.3515130.0193230.014779KAT5DEG_down2452.7176970.0282240.3679590.007160.007732SMARCA4DEG_down2422.7157580.0216460.3682210.0101160.008636RPL13DEG_down2342.7551520.0157730.3629560.0245410.011238PLCG1DEG_down2272.8850910.0300110.346610.0021440.009614FOSDEG_up2262.8807270.0245380.3471350.0032650.009974AP2M1DEG_down2142.8501820.0255820.3508550.0076350.010848SF3B3DEG_down2032.7655760.0151880.3615880.0145830.010661TP53RKDEG_up2002.8713940.0139720.3482630.016080.012798MMS19DEG_down1952.8069090.0223350.3562640.0047050.008763RPS7DEG_down1932.883030.0081260.3468570.0283350.016367MEPCEDEG_down1892.9052120.01450.3442090.0050660.012337TUBB4BDEG_down1732.8803640.0081670.3471780.0075950.019412CDH1DEG_down1712.8884850.0161860.3462020.0058480.011775EZRDEG_down1702.9066670.0162950.3440370.0041070.01554DDX41DEG_down1682.9500610.0095450.3389760.0070570.017509NUDCDEG_up1602.9437580.0168860.3397020.002830.013954SF3A2DEG_down1562.9498180.0090750.3390040.0099260.01593RPL27DEG_down1472.9127270.0039490.3433210.0439850.020057NIFKDEG_up1453.1078790.0111110.3217630.0103450.016053SRPRBDEG_up1432.9453330.0175450.339520.0016740.015983ATP1A1DEG_up1402.8166060.0112690.3550370.0091470.013032


## 3.4. Crosstalk Genes Bridging Alzheimer’s Disease and Periodontitis

In total, 51 upregulated crosstalk genes and 41 downregulated crosstalk genes were identified and are listed in Table4. The crosstalk gene-related PPI network as shown in Figure 5 consisted of 3496 nodes and 5141 edges. The topological characteristics of the top 30 nodes in this network are presented in Table 5. Hub crosstalk genes with the highest degree included MYC, HSPB1, VIM, KAT5, RPL13, FOS, and CDH1.Table 4
The 51 crosstalk genes that were upregulated and 41 crosstalk genes which were downregulated in periodontitis.
Regulation pattern in periodontitisCrosstalk genes linking periodontitis and Alzheimer’s diseaseUpregulatedAMIGO1 (adhesion molecule with Ig-like domain 1; gene ID: 57463)ARG2 (arginase 2; gene ID: 384)BDKRB2 (bradykinin receptor B2; gene ID: 624)BRI3 (brain protein I3; gene ID: 25798)C1D (C1D nuclear receptor corepressor; gene ID: 10438)C3 (complement C3; gene ID: 718)C4A (complement C4A (Rodgers blood group); gene ID: 720)C4B (complement C4B; gene ID: 721)C4B_2 (complement component 4B (Chido blood group), copy 2; gene ID: 100293534)CASP7 (caspase 7; gene ID: 840)CD177 (CD177 molecule; gene ID: 57126)CD38 (CD38 molecule; gene ID: 952)CHCHD10 (coiled-coil-helix-coiled-coil-helix domain containing 10; gene ID: 400916)COX8A (cytochrome c oxidase subunit 8A; gene ID: 1351)CSF3 (colony-stimulating factor 3; gene ID: 1440)CTGF (cellular communication network factor 2; gene ID: 1490)CTSS (cathepsin S; gene ID: 1520)CXCL1 (C-X-C motif chemokine ligand 1; gene ID: 2919)CXCL12 (C-X-C motif chemokine ligand 12; gene ID: 6387)CXCL8 (C-X-C motif chemokine ligand 8; gene ID: 3576)CXCR4 (C-X-C motif chemokine receptor 4; gene ID: 7842)ENPP2 (ectonucleotide pyrophosphatase/phosphodiesterase 2; gene ID: 5168)FCGR3A (Fc fragment of IgG receptor IIIa; gene ID: 2214)FCGR3B (Fc fragment of IgG receptor IIIb; gene ID: 2215)FHL2 (four and a half LIM domains 2; gene ID: 2274)FOS (Fos proto-oncogene, AP-1 transcription factor subunit; gene ID: 2353)FOXO1 (forkhead box O1; gene ID: 2308)GSTK1 (glutathione S-transferase kappa 1; gene ID: 373156)GSTO1 (glutathione S-transferase omega 1; gene ID: 9446)HCLS1 (hematopoietic cell-specific Lyn substrate 1; gene ID: 3059)IGFBP7 (insulin-like growth factor binding protein 7; gene ID: 3490)IL1B (interleukin-1 beta; gene ID: 3553)LYZ (lysozyme; gene ID: 4069)MAK16 (MAK16 homolog; gene ID: 84549)MEF2C (myocyte-specific enhancer factor 2C; gene ID: 4208)MME (membrane metalloendopeptidase; gene ID: 4311)MMP1 (matrix metallopeptidase 1; gene ID: 4312)MMP3 (matrix metallopeptidase 3; gene ID: 4314)MS4A1 (membrane-spanning 4-domains A1; gene ID: 931)MSRA (methionine sulfoxide reductase A; gene ID: 4482)MZB1 (marginal zone B and B1 cell-specific protein; gene ID: 51237)NDUFB8 (NADH:ubiquinone oxidoreductase subunit B8; gene ID: 4714)PECAM1 (platelet endothelial cell adhesion molecule 1; gene ID: 5175)PLAT (plasminogen activator, tissue type; gene ID: 5327)PSENEN (presenilin enhancer, gamma-secretase subunit; gene ID: 55851)SEL1L (SEL1L adaptor subunit of ERAD E3 ubiquitin ligase; gene ID: 6400)SERPINI1 (serpin family I member 1; gene ID: 5274)ST6GAL1 (ST6 beta-galactoside alpha-2,6-sialyltransferase 1; gene ID: 6480)VCAN (versican; gene ID: 1462)VIM (vimentin; gene ID: 7431)WARS2 (tryptophanyl tRNA synthetase 2, mitochondrial; gene ID: 10352)DownregulatedABCA12 (ATP binding cassette subfamily A member 12; gene ID: 26154)ACE (angiotensin I converting enzyme; gene ID: 1636)ADRB1 (adrenoceptor beta 1; gene ID: 153)BCL2L2 (BCL2-like 2; gene ID: 599)BRSK1 (BR serine/threonine kinase 1; gene ID: 84446)CALML5 (calmodulin-like 5; gene ID: 51806)CALR (calreticulin; gene ID: 811)CD36 (CD36 molecule; gene ID: 948)CDH1 (cadherin 1; gene ID: 999)CDK5R1 (cyclin-dependent kinase 5 regulatory subunit 1; gene ID: 8851)CTDSP2 (CTD small phosphatase 2; gene ID: 10106)DNM1 (dynamin 1; gene ID: 1759)DSG1 (desmoglein 1; gene ID: 1828)FLG (filaggrin; gene ID: 2312)GJA1 (gap junction protein alpha 1; gene ID: 2697)GRHL3 (grainyhead-like transcription factor 3; gene ID: 57822)GRIN3B (glutamate ionotropic receptor NMDA type subunit 3B; gene ID: 116444)HM13 (histocompatibility minor 13; gene ID: 81502)HMGCR (3-hydroxy-3-methylglutaryl-CoA reductase; gene ID: 3156)HSPB1 (heat shock protein family B (small) member 1; gene ID: 3315)KAT5 (lysine acetyltransferase 5; gene ID: 10524)MAP3K12 (mitogen-activated protein kinase kinase kinase 12; gene ID: 7786)MED12 (mediator complex subunit 12; gene ID: 9968)MFN2 (mitofusin 2; gene ID: 9927)MYC (MYC proto-oncogene, bHLH transcription factor; gene ID: 4609)NEFL (neurofilament light; gene ID: 4747)NES (exportin 1; gene ID: 7514)NGF (nerve growth factor; gene ID: 4803)NOS3 (nitric oxide synthase 3; gene ID: 4846)NPTXR (neuronal pentraxin receptor; gene ID: 23467)PLXNA3 (plexin A3; gene ID: 55558)PPIL2 (peptidylprolyl isomerase-like 2; gene ID: 23759)PTGS1 (prostaglandin-endoperoxide synthase 1; gene ID: 5742)RPL13 (ribosomal protein L13; gene ID: 6137)SCIMP (SLP adaptor and CSK interacting membrane protein; gene ID: 388325)SPPL2B (signal peptide peptidase-like 2B; gene ID: 56928)TIAF1 (TGFB1-induced antiapoptotic factor 1; gene ID: 9220)TP73 (tumor protein P73; gene ID: 7161)TYRP1 (tyrosinase-related protein 1; gene ID: 7306)USF2 (upstream transcription factor 2, C-Fos interacting; gene ID: 7392)VEGFA (vascular endothelial growth factor A; gene ID: 7422)Figure 5
The crosstalk gene-related PPI network.Table 5
Topological characteristics of the top 30 nodes in the crosstalk gene-related PPI network, ranked in descending order of degree.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientMYCDEG_AD_down10272.1275250.4666980.470035.03E−040.001842HSPB1DEG_AD_down4332.369590.1845980.4220140.0016360.003351VIMDEG_AD_up3202.6633010.1166110.3754740.0017630.005783KAT5DEG_AD_down2452.5663590.0894760.3896570.0030110.005621RPL13DEG_AD_down2342.5897290.0897080.3861410.0013940.005672FOSDEG_AD_up2262.8831510.0836420.3468431.97E−040.017606CDH1DEG_AD_down1712.931910.0584280.34107500.025911TP73DEG_AD_down1372.614830.0467810.3824340.0063330.009173FHL2DEG_AD_up1362.6410850.0571920.3786320.0040310.00889CALRDEG_AD_down1352.8124640.0461420.355560.0027640.010949GSTK1DEG_AD_up1112.9884590.0351520.33462100.026289MMEDEG_AD_up1042.8369880.0403860.3524870.0013070.011694MED12DEG_AD_down912.6956150.0306820.3709730.0061050.012517CASP7DEG_AD_up902.9010390.0271310.3447040.0052430.015603DNM1DEG_AD_down842.9601850.0290860.3378170.0011470.023716FOXO1DEG_AD_up672.9777840.0174030.335820.0040710.026586CXCR4DEG_AD_up652.9985570.0208510.3334940.0038460.028808NOS3DEG_AD_down622.9636470.0187860.3374220.0026440.027957NDUFB8DEG_AD_up563.0424120.021360.32868700.027447UBCOther gene532.0848240.1492650.4796570.0159650.027336NESDEG_AD_down532.9209460.0170120.3423550.0079830.022954CDK5R1DEG_AD_down503.01760.01330.3313890.0024490.037222NEFLDEG_AD_down492.9134450.0148860.3432360.0085030.024673C3DEG_AD_up433.0432780.017050.32859300.042101SEL1LDEG_AD_up423.0360650.0130660.32937400.045635DSG1DEG_AD_down413.0461630.008570.32828200.068464GJA1DEG_AD_down413.0351990.0157030.32946800.04037MEF2CDEG_AD_up413.0392380.0086760.3290300.068892LYZDEG_AD_up403.0421230.0085780.32871800.065086GSTO1DEG_AD_up382.8202540.0109050.3545780.0199150.030356


## 3.5. Crosstalk Genes Obtained by Feature Selection

A total of 48 crosstalk genes were selected by using the Boruta algorithm (Figure6(a)). In addition, 62 crosstalk genes were selected by using the RFE algorithm (Figure 6(b)). All 48 genes obtained by using the Boruta algorithm were included in the 62 genes obtained by the RFE algorithm, indicating that this 48-gene set was representative of the characteristics of all 92 crosstalk genes.Figure 6
48 crosstalk genes selected by using the Boruta algorithm (a) and 62 crosstalk genes selected by using the RFE algorithm (b).
(a)
Boruta algorithm to obtain genes(b)
RFE algorithm to obtain genes


## 3.6. Classification Accuracy Using the Feature Selection-Obtained Crosstalk Genes

The 48 genes obtained by SVM feature selection were not shown in all of the four periodontitis-related datasets. The gene expression profiles of GSE16134 and GSE10334 included all of the 48 crosstalk genes, whereas GSE23586 included 44 crosstalk genes and GSE79705 included 45 crosstalk genes. The classification performance of these 48 crosstalk genes for the four datasets is shown in Table6. For the test set GSE16134 and the validation set GSE10334, the accuracy performance was high at 91.94% and 88.26%, respectively. By comparison, the performance for the other two datasets GSE23586 and GSE79705 was low at 50% and 66.66%, respectively.Table 6
Classification performance of the 48 feature selection-obtained crosstalk genes in the four periodontitis-related datasets.
GSE16134GSE10334GSE23586GSE79705Samples: test setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracySamples: validation setSamples: SVM model predicts accuratelyAccuracyCase939291.94%18317788.26%3350%8866.66%Control312264413040Total12411424721863128


## 3.7. ROC Curves for the Four Periodontitis-Related Datasets

As shown in Figure7, the AUC (area under the curve) values for the GSE16134 test set and the GSE10334 validation set were high at 95.77% and 90.53%, respectively, congruent with the results in Table 6. It was thus inferred that the classifier performance was adequate only when the sample sizes of the validation sets were similar to those of the training and test datasets, and therefore, poor performance was noted for GSE23586 and GSE79705 having much lower sample numbers.Figure 7
The ROC curves of four periodontitis-related datasets.


## 3.8. The Identification of Transcription Factors Regulating the Crosstalk Genes

As shown in Figure8, the TF-crosstalk gene target network consisted of 388 nodes and 1178 edges. Several transcription factors which were also DEGs played critical roles by regulating the most number of crosstalk genes, for example, FOS, MEF2C, and USF2 (Table 7).Figure 8
The transcription factor-crosstalk gene target network.Table 7
Topological characteristics of the top 30 nodes in the TF-crosstalk gene target network, ranking in descending order.
NodeLabelDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficientTopological coefficientFOSDEG_up&TF892.1782950.2061010.4590750.0204290.03825VEGFADEG_down772.2118860.1972190.4521030.0109360.038456MEF2CDEG_up&TF762.3281650.1599450.4295230.0042110.044258CDH1DEG_down692.5400520.1593920.39369300.086634CXCR4DEG_up472.5968990.0751050.3850759.25E−040.083836MMEDEG_up442.6640830.043650.37536400.131313C4ADEG_up382.6614990.0297420.3757280.011380.118617VIMDEG_up382.5219640.0602930.3965160.0128020.065789CSF3DEG_up362.5503880.0317340.3920970.0190480.075137C4BDEG_up342.6795870.0188870.3731920.0124780.122037VCANDEG_up342.7260980.0305670.36682500.143717IL1BDEG_up332.7416020.057460.3647500.134602BRSK1DEG_down302.7416020.0304040.3647500.153788SERPINI1DEG_up292.7932820.0190140.35800200.127586USF2DEG_down&TF272.5917310.0275750.3858420.0455840.081607GRHL3DEG_down272.7751940.0270670.36033500.150732HSPB1DEG_down262.7803620.0179810.35966500.155678MMP3DEG_up262.6072350.0243770.3835480.0369230.074984CTGFDEG_up252.7881140.0206690.35866500.161905HMGCRDEG_down252.7674420.0411190.36134500.16GATA2TF242.0956070.0526870.4771890.0181160.104678PLATDEG_up242.7622740.0280750.36202100.173148ETS1TF222.2739020.037670.4397730.0129870.092022NPTXRDEG_down222.8268730.0204780.35374800.154545ARTF212.240310.0406980.4463670.023810.105465YBX1TF212.3152450.0378010.431920.0047620.098928BCL2L2DEG_down212.8139530.0181790.35537200.144048FCGR3ADEG_up212.811370.0205210.35569900.192799ABCA12DEG_down202.798450.0109610.35734100.2SP1TF202.1912140.0373620.4563680.0263160.109167CDK5R1DEG_down202.8242890.0178780.35407100.156098FLGDEG_down202.6537470.0122810.3768260.0421050.093043


## 3.9. Signaling Pathways Enriched in the Crosstalk Genes

From the 48 feature selection-obtained crosstalk genes, 37 crosstalk genes were found among gene-pathway interaction pair data in the KEGG database. 137 KEGG pathways corresponded to these 37 crosstalk genes. Figure9 shows the expression values of these 137 pathways in the four periodontitis-related datasets.Figure 9
Heatmaps depicting 137 signaling pathways enriched in the crosstalk genes in the four periodontitis-related datasets.
(a)
GSE10334(b)
GSE16134(c)
GSE23586(d)
GSE79705The numbers of crosstalk gene-related DEPs obtained from each of the four periodontitis-related datasets are listed in Table8. The interaction relationships between crosstalk genes and DEPs are depicted in Figure 10, showing that several DEPs were dysregulated in at least two datasets, including cytokine-related pathways (cytokine-cytokine receptor interaction, chemokine, and IL-17), immune cell-related pathways (T cell receptor, B cell receptor, Th1 and Th2 cell differentiation, Th17 cell differentiation, natural killer cell-mediated cytotoxicity, and osteoclast differentiation), JAK-STAT signaling, NOD-like receptor signaling, MAPK signaling, Toll-like receptor signaling, NF-kappa B signaling, and C-type lectin receptor signaling.Table 8
The number of crosstalk gene-related DEPs within the four periodontitis-related datasets.
DataSampleDEPDEPDEPLog FCp valueUpDownTotalAbsGSE23586611213>1<0.05GSE1613431084286>1<0.05GSE1033424777279>1<0.05GSE7970512448>0<0.05Figure 10
The crosstalk gene-differentially expressed pathway interaction network.


## 3.10. The Identification of Core Crosstalk Genes

Among the 48 feature selection-obtained crosstalk genes and periodontitis-related genes in the DisGeNET database, 12 common genes were identified as core crosstalk genes including C3, C4A, C4B, CXCL12, FCGR3A, FCGR3B, HSPB1, IL1B, MME, MMP3, PLAT, and VEGFA (Table9). Among these 12 genes, 6 genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, were found associated with chronic periodontitis.Table 9
The 12 overlapping genes between 48 feature selection-obtained crosstalk genes and periodontitis-related genes in the DisGeNET database. Among the 12 genes, 6 genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, were associated with chronic periodontitis.
Gene symbolGene IDDisease IDDisease nameC3718C0031099PeriodontitisC4A720C0266929Chronic periodontitisC4B721C0266929Chronic periodontitisCXCL126387C0266929Chronic periodontitisFCGR3A2214C0266929Chronic periodontitisFCGR3A2214C0031099PeriodontitisFCGR3A2214C0031106Periodontitis, juvenileFCGR3A2214C0399447Early-onset periodontitisFCGR3A2214C0031030Periapical periodontitisFCGR3B2215C0399447Early-onset periodontitisFCGR3B2215C0031106Periodontitis, juvenileHSPB13315C0031030Periapical periodontitisIL1B3553C1719494Periodontitis, localized aggressiveIL1B3553C1719495Aggressive periodontitis, generalizedIL1B3553C0399447Early-onset periodontitisIL1B3553C0001342Acute periodontitisIL1B3553C0031106Periodontitis, juvenileIL1B3553C0031030Periapical periodontitisIL1B3553C4025886Severe periodontitisIL1B3553C0031099PeriodontitisIL1B3553C0266929Chronic periodontitisMME4311C0031099PeriodontitisMMP34314C0031106Periodontitis, juvenileMMP34314C0001342Acute periodontitisMMP34314C0266929Chronic periodontitisMMP34314C0031099PeriodontitisPLAT5327C1719495Aggressive periodontitis, generalizedVEGFA7422C0031099PeriodontitisFigure11 shows ROC curves for the 6 chronic periodontitis-related genes for each of the four datasets. An AUC value of more than 80% was presented by three genes, C4A, C4B, and CXCL12, for the dataset GSE10334, 5 genes, C4A, C4B, CXCL12, FCGR3A, and IL1B, for GSE16134, 3 genes, CXCL12, IL1B, and MMP3, for GSE23586, and all 6 genes for GSE79705. The 3 genes, C4A, C4B, and CXCL1, had the highest classification accuracy when the datasets with small samples (GSE23586 and GSE79705) were not considered.Figure 11
The ROC curves for 6 chronic periodontitis-related genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, in the four periodontitis-related datasets.
(a)
GSE10334(b)
GSE16134(c)
GSE23586(d)
GSE79705


## 4. Discussion

The present study addressed shared genetic mechanisms and molecular links between periodontitis and Alzheimer’s diseases by identifying gene expression, signaling pathways, and TFs that were most robustly associated with both these diseases. These findings are largely substantiated by preexisting experimental data.Six genes C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3 were identified as the most significant crosstalk genes linking chronic periodontitis and Alzheimer’s disease. C4B and C4A, respectively, encode the basic and acidic forms of the complement factor 4, and C4 gene deficiency has been noted to predispose the development of severe chronic periodontitis [38]. In AD, the expression levels of C4 mRNA were shown to be 3.27-fold increased in temporal cortex samples as compared to controls [39]. The CXCL12 (C-X-C motif chemokine ligand 12) gene is the ligand of the C-X-C motif chemokine receptor 4 (CXCR4). CXCL12 expression in the gingival crevicular fluid of periodontitis patients was shown to be significantly higher than that of healthy subjects, suggesting that it might play a role in enhancing neutrophil migration and further the progression of periodontitis [40]. A decreased level of CXCL12 in Alzheimer’s disease has been documented as affecting cognitive function, impairing learning and memory [41]. FCGR3A (Fc fragment of IgG receptor IIIa) encodes a receptor for the Fc portion of immunoglobulin G, and FCGR3A polymorphisms are shown to confer susceptibility to periodontitis in Caucasians [42]. In AD, the Fc gamma receptor (FcγR) was recently found to exacerbate neurodegeneration [43]. Cytokines are considered a primary link between chronic periodontitis and Alzheimer’s disease as they can enter systemic circulation through periodontal pockets [6]. The classical proinflammatory cytokine, IL1B, is elevated in periodontitis and can induce resorption of alveolar bone [44]. In AD, IL1B gene polymorphisms are linked to disease susceptibility [45]. MMP3 (matrix metalloproteinase 3) is implicated in the progression of chronic periodontitis and can degrade the periodontal tissue matrix [46]. Elevated brain levels of MMP3 have been associated with the duration of Alzheimer’s disease, and it has been found to increase the activity of MMP9, thereby indirectly promoting aggregation and cerebral accumulation of tau deposits [47].More interestingly, the six genes discussed in the last paragraph were also found to be the molecular crosstalks in linking the peripheral immune system and central nervous system (CNS). It has been well demonstrated that the periodontal disease-evoked peripheral systemic host immune response can aggravate the progression of neuroinflammation and neurodegeneration in Alzheimer’s disease by switching the microglia from the primed phenotype to an aggressive proinflammatory phenotype [48, 49]. The immune-inflammatory mediators (e.g., cytokines and chemokines) abundantly expressed during periodontal inflammation can circulate into the bloodstream and travel into the brain by crossing the blood-brain barrier (BBB) and impact the function of CNS [50]. Therefore, the crosstalk between the peripheral immune system and the CNS might be an important mechanism underlying periodontitis, increasing the risk of AD. This paragraph will provide a description regarding the potential role of the six crosstalk genes in linking periodontal disease and AD, especially by means of neuroimmune interaction. For example, the complement components C4A and C4B highly expressed in periodontal disease were found to modulate T cell immune response by stimulating the activation and migration of T cells [51–53]. The migration of T cells enhanced by C4A and C4B might allow T cells to traffic across the BBB and enter the brain. For another example, the chemokine and its receptor-composed system CXCL12/CXCR4-7 system were found to be a significant player of the neuroimmune interface [54]. On the one hand, the chemokine CXCL12 mediated the immune-inflammatory response in CNS by recruiting lymphocytes and macrophages [55]. On the other hand, CXCL12 can lead to neurotoxicity and neurodegenerescence by activating the neuronal survival-associated G protein-activated inward rectifier K(+) (GIRK) [54]. FCGR3A (also named CD16) is essential for the antibody-dependent cellular cytotoxicity (ADCC) mediated by natural killer (NK) cells [56]. The increased cytotoxic activity of NK cells was found to cause the dysregulation of protein kinase C and further led to the cognitive deficits in Alzheimer’s disease, indicating the contribution of immunological factors to the dysfunction of CNS [57]. IL1B, which was upregulated in periodontitis and transported through the vascular circulation into the brain, was found to play promoting roles in neuroinflammation by enhancing the expression of leukocyte chemotactic chemokines, cell surface adhesion molecules, cyclooxygenases, and MMPs within the brain parenchyma [58]. Likewise, MMP3 abundantly produced in periodontitis was also found to be associated with neuroinflammation via activating microglial cells, as well as participating in the BBB breakdown through the proteolysis of fibronectin and type IV collagen [59]. Taken together, the six crosstalk genes identified in the present research were well evidenced to be involved in periodontitis-triggered peripheral systemic host immune response caused CNS dysfunction in Alzheimer’s disease.Three transcription factors, FOS, MEF2C, and USF2, were identified as related to the regulation of the crosstalk genes and were also found to be dysregulated in chronic periodontitis. The proto-oncogene FOS (also named C-Fos) was found to be involved in the transcriptional regulation of collagenase and cell proliferation genes in periodontal gingival fibroblasts [60]. In AD, FOS is reported to initiate amyloid-β-mediated apoptosis and found to be increased in the hippocampal regions of AD patients [61]. Myocyte-specific enhancer factor 2C (MEF2C) was identified as a critical transcription factor involved in the coexpression network of chronic periodontitis [62]. Genome-wide association studies (GWAS) have shown the linkage between mutation of MEF2C and aging-associated late-onset Alzheimer’s disease [63, 64]. Experimentally, a lack of MEF2C expression was shown to exaggerate microglial response and negatively affect brain function [65]. The USF2 (upstream transcription factor 2, C-Fos interacting) transcription factor is reported to enhance osteogenic differentiation of periodontal ligament cells (PDLCs) [66]; however, its involvement in periodontal inflammation has not been reported. In the context of AD, the USF2 gene was shown to regulate the expression of genes Dhcr24, Aplp2, Tia1, Pdrx1, Vdac1, and Syn2, which drive the neuropathological mechanisms [67]. However, a study of Japanese participants found that the single nucleotide polymorphisms of the USF2 gene were not significantly related to the onset of AD [68].Differentially expressed pathways were identified from the crosstalk gene-pathway network, and several pathways including JAK-STAT, MAPK, NF-kappa B, and natural killer cell-mediated cytotoxicity were found as the most robust differentially expressed pathways in at least two periodontitis datasets. Overall, experimental evidence supports these as linkage mechanisms between periodontitis and AD. The activation of the JAK-STAT pathway induced by thePorphyromonas gingivalis lipopolysaccharide (LPS) and nicotine was shown to increase the expression of cyclooxygenase-2 (COX-2), prostaglandin E2 (PGE2), and proinflammatory cytokines in osteoblasts, thus further accelerating periodontitis progression [69, 70]. In AD, the inhibitor of the JAK-STAK pathway is reported as a therapeutic target, and blocking this pathway can protect against neuroinflammation and dopaminergic neurodegeneration [71]. The MAPK pathway is the upstream signaling intermediate to many inflammatory cytokines such as TNF-α, IL-1β, IL-6, and prostaglandin E2 [72], and the blockage of this pathway could be beneficial for treating inflammatory diseases like chronic periodontitis and AD [73, 74]. The activation of MAPK signaling is noted to promote the production of MMPs and RANKL, leading to osteoclastogenesis and the acceleration of alveolar bone loss [73]. MAPK signaling is also implicated in multiple aspects of the neuropathology of AD, such as promoting neuroinflammation, amyloid-beta toxicity and aggregation, autophagy, and apoptosis [75]. The overexpression of NF-κB signaling plays a pivotal role in periodontitis-associated bone destruction by promoting the differentiation and activation of osteoclasts [76]. The blockade of NF-κB signaling is found to trigger detrimental neural alterations including neuroinflammation, activation of microglia, oxidative stress-related complications, and apoptosis [77]. Natural killer (NK) cells are important regulators of innate and adaptive immunity and are closely linked to the regulation of cytotoxicity [78]. Experimental data shows NK cells can directly recognize the Fusobacterium nucleatum pathogen, leading to alveolar bone resorption and periodontitis [79]. The overactivity of NK is also purported to play a driving role in the progression of AD by producing a series of proinflammatory cytokines [80].The findings of thisin silico analysis must be considered in light of the strengths and limitations of this work. By using a machine learning-based feature selection method as the core technique, the most putatively robust crosstalk genes could be identified. Furthermore, functional molecular links were also analyzed in terms of differentially expressed pathways. The integrated analysis of multiple periodontitis-related GEO datasets enabled a larger sample size for improved accuracy of our computational prediction in the present study. The major limitation of the current approach is that no experimental validation of the identified pivotal genetic molecular linkage candidates was performed. This work has multiple implications for future research. Experimental and clinical studies focused on these candidates could be valuable from the perspectives of identification of shared susceptibility, exaggerating pathogenic mechanisms, biomarkers, and therapeutic targets relevant to precision medicine and drug development or repurposing.


## 5. Conclusion

Bioinformatic analysis integrating experimental transcriptomic data from Alzheimer’s disease and periodontitis revealed the most robust potentially shared molecular linkages. Six crosstalk genes, C4A, C4B, CXCL12, FCGR3A, IL1B, and MMP3, three transcription factors, FOS, MEF2C, and USF2, and several pathways, JAK-STAT, MAPK, NF-kappa B, and natural killer cell-mediated cytotoxicity, emerged as top candidate shared molecular linkage entities and merit future research in experimental and clinical studies.



---

*Source: 6633563-2021-04-01.xml*