Transcriptomic similarities and conserved genes between human periodontitis and mouse ligature-induced periodontitis: A secondary analysis of gene expression omnibus datasets

Article information

J Korean Dent Assoc. 2025;63(10):322-328
Publication date (electronic) : 2025 October 31
doi : https://doi.org/10.22974/jkda.2025.63.10.002
Department of Periodontics and Dental Research Institute, Pusan National University Dental Hospital, Yangsan, Korea
*Corresponding author: Dr. Shin-Kyu Lee Department of Periodontics and Dental Research Institute, Pusan National University Dental Hospital, Beomeo-ri, Mulgeum-eup, Yangsan 50612, Korea Tel: +82-55-360-5190, E-mail: pouiw2000@gmail.com
Received 2025 July 23; Revised 2025 September 24; Accepted 2025 October 13.

Abstract

Purpose

This study aimed to evaluate how well the ligature-induced periodontitis mouse model reflects the transcriptomic features of human periodontitis and to identify periodontitis-associated genes and functions conserved across species.

Materials and Methods

RNA sequencing data from human and mouse gingival tissues were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using DESeq2. Functional enrichment analysis was performed using Gene Ontology and KEGG pathway databases. Cross-species transcriptomic similarity was evaluated by comparing DEG overlap and enrichment similarity between human and mouse.

Results

A total of 223 human DEGs and 622 mouse DEGs were identified. Among these, 25 DEGs were shared between species, including 23 showing concordant regulation direction. The mouse-to-human overlap ratio was 3.70%. Functional enrichment analysis identified 189 significant terms or pathways in humans and 602 in mice, with 129 shared results. Specifically, 41.18% of mouse KEGG pathway results overlapped with human results. The shared concordant genes, including IL-1β, PTGS2 (also known as COX-2), and MMP13, were associated with immune and inflammatory functions.

Conclusion

The ligature-induced periodontitis mouse model reflects the transcriptomic features of human periodontitis in a limited manner, showing low similarity at the DEG level and moderate similarity at the enrichment level. Conserved DEGs such as IL-1β, PTGS2, and MMP13 may represent fundamental molecular mechanisms of periodontitis.

Introduction

Periodontitis is a chronic inflammatory disease characterized by the destruction of periodontal tissues, potentially leading to tooth loss [1]. Various animal models were utilized to investigate the pathophysiological mechanisms of human periodontitis [2]. Among these, the ligature-induced periodontitis mouse model was widely used due to its genetic similarity to humans, cost-effectiveness, and ease of handling [3]. Histological evaluation of ligature-induced periodontitis in mice showed prominent inflammatory responses from approximately 4 days after ligation, progressively increasing until day 8, after which no significant changes were observed [4]. In mice, periodontitis was also induced by injecting lipopolysaccharide or Porphyromonas gingivalis [5]. However, these methods required 4–6 weeks to produce significant periodontal tissue destruction, which was much longer than the ligation method [6]. Therefore, the ligature-induced periodontitis mouse model was effectively used for studying periodontal tissue destruction-related changes.

In addition, numerous studies used the ligature-induced mouse model to explore the molecular characteristics of periodontitis [6-8]. Differentially expressed genes (DEGs) in ligature-induced periodontitis were significantly enriched in functions such as neutrophil chemotaxis and inflammatory responses, and the expression of innate immune response-related genes S100a8 and S100a9 was markedly increased at both RNA and tissue levels in periodontitis [7]. Substance P played a key role in regulating inflammation and alveolar bone loss in ligature-induced periodontitis, whereas calcitonin gene-related peptide had a minimal effect [8]. Through such studies, the molecular characteristics of periodontitis were elucidated using the ligature-induced periodontitis mouse model. However, to the best of our knowledge, no study has directly evaluated how similar the molecular profiles observed in the ligature-induced periodontitis mouse model are to those in human periodontitis. In other words, it remains unclear to what extent this mouse model reflects the gene expression features of human periodontitis.

The aim of this study was to evaluate how well ligature-induced periodontitis in mice reflects the transcriptomic features of human periodontitis and furthermore to identify conserved periodontitis-related genes and their functions across species.

Materials and Methods

The overall workflow of this study is presented in Figure 1.

Figure 1.

Workflow of the study. DEGs: differentially expressed genes

Data collection and preprocessing

RNA sequencing data from both human and mouse gingival tissues were obtained from the publicly accessible Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). Human data were retrieved from GSE223924, including 10 healthy gingiva samples and 10 periodontitis gingiva samples. All human subjects were systemically healthy. Periodontitis cases corresponded to stage II and III disease, classified according to the 2017 World Workshop criteria, while patients with Grade C (rapid progression or early onset) were excluded. Inflamed gingival tissues were obtained during open-flap surgery at diseased sites, and healthy tissues were collected during crown-lengthening or orthodontic extractions. Mouse data were retrieved from GSE184556, including 3 healthy gingiva samples and 3 ligature-induced periodontitis samples. In this dataset, periodontitis was induced by placing a ligature around the second maxillary molar, and gingival tissues were collected at 8 days post-induction. Three age- and sex-matched mice were used per group. In both groups, gingival tissues were harvested from the maxillary right molar region. Detailed characteristics of each sample can be found in the series matrix files available in the corresponding GEO entries. Raw count data were downloaded, and ortholog mapping between species was performed using Ensembl BioMart (Ensembl, European Molecular Biology Laboratory–European Bioinformatics Institute [EMBL–EBI], Hinxton, Cambridge, UK) in R (version 4.4.3; R Foundation for Statistical Computing, Vienna, Austria).

Differential gene expression analysis

Differential expression analysis was performed independently for the human and mouse datasets using the DESeq2 package (version 1.46.0) in R. Genes with an absolute log2 fold change ≥ 1 and a Benjamini-Hochberg adjusted p-value < 0.05 were defined as DEGs.

Functional enrichment analysis

Functional enrichment analysis was conducted using g:Profiler (version e113_eg59_p19_f6a03c19, https://biit.cs.ut.ee/gprofiler/). Analyses included Gene Ontology (GO) categories — Molecular Function (MF), Biological Process (BP), and Cellular Component (CC) — and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Significance was determined using the g:SCS-adjusted p-value provided by g:Profiler, with values < 0.05 considered significant.

Cross-species similarity analysis

To evaluate similarity between human and mouse DEGs, the mouse-to-human overlap ratio and the Jaccard similarity index were calculated. The mouse-to-human overlap ratio was defined as the proportion of mouse DEGs that were also identified as human DEGs. The Jaccard similarity index was calculated as the size of the intersection divided by the size of the union between the human and mouse DEG sets. In addition, similarity between the functional enrichment results derived from human and mouse DEGs was assessed using the same indices.

Results

Cross-species similarity of periodontitis-associated DEGs between human and mouse

Periodontitis-associated DEGs totaled 223 in humans and 622 in mice. Of these, 25 DEGs were shared between the two species, including 23 shared concordant DEGs showing the same regulation direction. Shared concordant DEGs accounted for 3.70% of mouse DEGs, and the Jaccard similarity was 2.80%.

Cross-species functional similarity based on periodontitis-associated DEGs

Functional enrichment analysis based on DEGs identified 189 significant GO terms or KEGG pathways in humans and 602 in mice (Table 1). Among these, 129 terms or pathways were shared between the two species. Overall, 21.43% of the mouse results were shared with human results, and the Jaccard similarity was 19.49%. Specifically, 41.18% of the KEGG pathway results from mice were shared with humans, while 9.52% of the GO:MF terms overlapped.

Similarity of functional enrichment results between human and mouse DEGs

Functional and network analysis of shared concordant DEGs between human and mouse

A total of 23 shared concordant DEGs between human and mouse were identified (Table 2). Protein–protein interaction (PPI) relationships among these genes in human were assessed and visualized (Fig. 2). IL-1β showed the highest node degree, followed by PTGS2, LOX, and MMP13. ZFP57 and CCER2 exhibited the largest absolute log₂ fold change, while both had a node degree of 0. Functional enrichment analysis of the shared concordant DEGs identified 79 GO terms (3 MF, 68 BP, 8 CC) and 3 KEGG pathways as significantly enriched (Fig. 3; only the driver terms are shown for GO).

List of shared concordant DEGs between human and mouse with fold change

Figure 2.

Protein–protein interaction networks. The five genes with the highest node degrees are IL1B (11), LOX, MMP13, and PTGS2 (each 7), and KDR (6).

Figure 3.

Functional enrichment results for shared concordant DEGs between human and mMouse. For GO terms, only the driver terms—representative key functions identified through grouping related terms in g:Profiler analysis—are shown

Discussion

To our knowledge, this is the first study to analyze DEGs in both human and mouse periodontitis. This study provides evidence for evaluating whether the ligature-induced periodontitis mouse model is suitable for exploring transcriptomic features of human periodontitis. We also identified periodontitis-associated DEGs conserved between human and mouse and characterized their functions.

In the present study, the mouse-to-human overlap was 3.70% at the DEG level and 21.43% at the functional enrichment level, indicating that the ligature-induced periodontitis mouse model only partially reflected human periodontitis transcriptomic features. Meanwhile, 41.18% of the KEGG pathway results from the mouse overlapped with those from humans. This suggested that the ligature-induced periodontitis mouse model might have more effectively reproduced pathway-level functions than individual DEGs. A previous study in peri-implantitis conducted a similar analysis, performing RNA sequencing using a ligature-induced peri-implantitis mouse model and comparing the results with those from humans [9]. It reported a mouse-to-human overlap of shared concordant DEGs of 28.63% and a Pearson correlation coefficient of gene expression patterns of 0.6935. Taken together with our findings, this suggested that despite inducing inflammation in the periodontal tissue using a similar method, the transcriptomic reproducibility of the mouse model was higher in peri-implantitis than in periodontitis. This might have been attributed to the greater differences in tooth structure and morphology between humans and mice [10], compared with those in implant fixtures produced through standardized manufacturing processes.

A total of 23 shared concordant DEGs were identified between human and mouse periodontitis, showing consistent regulation direction across species. These shared concordant DEGs might have represented conserved molecular responses and underlying mechanisms of periodontitis, potentially playing a more fundamental role in the disease than other DEGs. Among them, IL-1β and PTGS2 were well-established periodontitis-associated genes. IL-1β, which showed the highest node degree in this study, is known as a key gene in the onset and progression of periodontitis [11]. PTGS2, also known as COX-2, was reported to be highly expressed in inflamed periodontal tissues [12] and played a critical role in periodontal inflammatory responses [13]. Our findings supported the notion that these genes were central regulators of periodontitis and suggested that they were involved in fundamental periodontitis-related biological responses conserved across species. In addition to IL-1β and PTGS2, among the remaining 21 shared concordant DEGs, ZFP57 [14], ADAMTS5 [15], OSM [16], MMP13 [17], APLNR [18], and TREM1 [19] were reported as potentially associated with periodontitis. The other 13 genes had not been clearly linked to periodontitis in previous studies. PPI analysis of the 23 shared concordant DEGs revealed significantly greater interactions compared to random protein sets of the same size. Furthermore, these genes were significantly enriched in GO terms related to immune responses and inflammation, such as metalloendopeptidase activity, positive regulation of acute inflammatory response, and regulation of MAPK cascade. They were also significantly enriched in the IL-17 signaling pathway, which is known to play a key role in periodontitis [20]. These findings suggested that the shared concordant DEGs might have regulated periodontitis through these biological functions, and further studies were needed to clarify their specific roles individually or collectively.

This study conducted a comparative analysis of DEGs identified by comparing healthy and periodontitis samples in both ligature-induced periodontitis mouse models and human subjects. Our findings are significant as one of the few attempts to evaluate how well the ligature-induced periodontitis mouse model reflected the transcriptomic features of human periodontitis. Furthermore, this study identified periodontitis-related DEGs that were conserved between mice and humans.

However, there were limitations in that this study utilized publicly available datasets from GEO, and the sample sizes for both mouse and human data were relatively small. In addition, because the datasets were derived from independent studies, the number of samples differed between species. Future studies with larger and more balanced sample sets are warranted to validate the cross-species comparisons. Moreover, it is important to assess whether this varies with the method and duration of periodontitis induction, to develop more relevant mouse models. In particular, comparative studies with other induction methods, such as lipopolysaccharide injection or models induced by periodontal pathogens like Porphyromonas gingivalis, are needed to further clarify their translational relevance. Also, there is a need to verify whether the conserved periodontitis-related DEGs identified in this study are consistently observed across different datasets and to clarify their specific roles in the pathogenesis of periodontitis. This study would serve as a starting point for such efforts and contributes to a fundamental understanding of the molecular mechanisms involved in the onset and progression of periodontitis.

Notes

Conflicts of Interest

None

Acknowledgment

Sincere thanks to Professor Hyun-Joo Kim for her kind advice and encouragement.

References

1. Slots J. Periodontitis: facts, fallacies and the future. Periodontol 2000 2017;75:7–23.
2. Oz HS, Puleo DA. Animal models for periodontal disease. J Biomed Biotechnol 2011;2011:754857.
3. Abe T, Hajishengallis G. Optimization of the ligature-induced periodontitis model in mice. J Immunol Methods 2013;394:49–54.
4. Yoon H, Jung BH, Yoo KY, Lee JB, Um HS, Chang BS, et al. Temporal changes of periodontal tissue pathology in a periodontitis animal model. J Periodontal Implant Sci 2023;53:248–58.
5. Rojas C, García MP, Polanco AF, González-Osuna L, Sierra-Cristancho A, Melgar-Rodríguez S, et al. Humanized mouse models for the study of periodontitis: an opportunity to elucidate unresolved aspects of its immunopathogenesis and analyze new immunotherapeutic strategies. Front Immunol 2021;12:663328.
6. Lin P, Niimi H, Ohsugi Y, Tsuchiya Y, Shimohira T, Komatsu K, et al. Application of ligature-induced periodontitis in mice to explore the molecular mechanism of periodontal disease. Int J Mol Sci 2021;22:8900.
7. Maekawa S, Onizuka S, Katagiri S, Hatasa M, Ohsugi Y, Sasaki N, et al. RNA sequencing for ligature induced periodontitis in mice revealed important role of S100A8 and S100A9 for periodontal destruction. Sci Rep 2019;9:14663.
8. Siddiqui YD, Nie X, Wang S, Abbasi Y, Park L, Fan X, et al. Substance P aggravates ligature-induced periodontitis in mice. Front Immunol 2023;14:1099017.
9. Tong ZA, Chen L, Shen L, Lu YF, Zhang JW, Qi YD, et al. A modified ligature-induced peri-implantitis murine model and RNA sequencing analysis compared with human subjects. J Clin Periodontol 2025;52:787–98.
10. Beniash E, Stifler CA, Sun CY, Jung GS, Qin Z, Buehler MJ, et al. The hidden structure of human enamel. Nat Commun 2019;10:4383.
11. Cheng R, Wu Z, Li M, Shao M, Hu T. Interleukin-1β is a potential therapeutic target for periodontitis: a narrative review. Int J Oral Sci 2020;12:2.
12. Morton RS, Dongari-Bagtzoglou AI. Cyclooxygenase-2 is upregulated in inflamed gingival tissues. J Periodontol 2001;72:461–9.
13. Ebersole JL, Dawson D 3rd, Emecen-Huja P, Nagarajan R, Howard K, Grady ME, et al. The periodontal war: microbes and immunity. Periodontol 2000 2017;75:52–115.
14. Hernández HG, Hernández-Castañeda AA, Pieschacón MP, Arboleda H. ZNF718, HOXA4, and ZFP57 are differentially methylated in periodontitis in comparison with periodontal health: epigenome-wide DNA methylation pilot study. J Periodontal Res 2021;56:710–25.
15. Tayman MA, Koyuncu İ, Köklü NÖ. Expression levels of A disintegrin-like metalloproteinase with thrombospondin motifs-4 and -5 (ADAMTS-4 and ADAMTS-5) in inflamed and healthy gingival tissues. Comb Chem High Throughput Screen 2020;23:168–76.
16. Thorat M, Pradeep AR, Garg G. Correlation of levels of oncostatin M cytokine in crevicular fluid and serum in periodontal disease. Int J Oral Sci 2010;2:198–207.
17. Hernandez M, Valenzuela MA, Lopez-Otin C, Alvarez J, Lopez JM, Vernal R, et al. Matrix metalloproteinase-13 is highly expressed in destructive periodontal disease activity. J Periodontol 2006;77:1863–70.
18. Cores Ziskoven P, Nogueira AV, Yoldaş O, Buduneli N, Wild PS, Koeck T, et al. Apelin - a new kid on the block in periodontology. Oral Health Prev Dent 2024;22:417–24.
19. Wu D, Weng Y, Feng Y, Liang B, Wang H, Li L, et al. Trem1 induces periodontal inflammation via regulating M1 polarization. J Dent Res 2022;101:437–47.
20. Feng Y, Chen Z, Tu SQ, Wei JM, Hou YL, Kuang ZL, et al. Role of interleukin-17A in the pathomechanisms of periodontitis and related systemic chronic inflammatory diseases. Front Immunol 2022;13:862415.

Article information Continued

Figure 1.

Workflow of the study. DEGs: differentially expressed genes

Figure 2.

Protein–protein interaction networks. The five genes with the highest node degrees are IL1B (11), LOX, MMP13, and PTGS2 (each 7), and KDR (6).

Figure 3.

Functional enrichment results for shared concordant DEGs between human and mMouse. For GO terms, only the driver terms—representative key functions identified through grouping related terms in g:Profiler analysis—are shown

Table 1.

Similarity of functional enrichment results between human and mouse DEGs

Source Counts Similarity
Human Mouse Shared Mouse-to-human overlap Jaccard similarity
GO: BP 151 470 108 22.98% 21.05%
GO: CC 15 52 8 15.38% 13.56%
GO: MF 13 63 6 9.52% 8.57%
KEGG 10 17 7 41.18% 35.00%
Total 189 602 129 21.43% 19.49%

Table 2.

List of shared concordant DEGs between human and mouse with fold change

Gene symbol Fold change
Human Mouse Human Mouse
Log2 (fold change) Adjusted p-value Log2 (fold change) Adjusted p-value
ZFP57 Zfp57 8.1520 0.0026 1.4868 0.0000
SFRP4 Sfrp4 4.9785 0.0228 2.1974 0.0000
ADAMTS5 Adamts5 3.7677 0.0128 1.4896 0.0000
ADAMTS6 Adamts6 3.6203 0.0006 1.1164 0.0417
PDZRN4 Pdzrn4 3.5159 0.0447 4.3179 0.0001
OSM Osm 3.4963 0.0049 2.4800 0.0011
MMP13 Mmp13 3.4531 0.0361 3.5774 0.0060
APLNR Aplnr 2.9780 0.0121 1.6531 0.0000
INHBA Inhba 2.9521 0.0049 2.0642 0.0000
TREM1 Trem1 2.9244 0.0351 1.7101 0.0279
B3GALT5 B3galt5 2.9140 0.0137 2.7406 0.0001
IL1B Il1b 2.8609 0.0031 1.8198 0.0004
TGFB3 Tgfb3 2.8176 0.0056 1.0572 0.0009
MMP16 Mmp16 2.8053 0.0133 1.3566 0.0399
PLEK Plek 2.5744 0.0292 1.4520 0.0032
CDH5 Cdh5 2.5088 0.0239 1.2211 0.0000
IL1RL1 Il1rl1 2.2247 0.0425 1.8857 0.0014
KDR Kdr 2.1794 0.0466 1.0046 0.0019
PTGS2 Ptgs2 2.0291 0.0037 1.6771 0.0405
RGS2 Rgs2 1.6455 0.0000 1.5615 0.0028
LOX Lox 1.0785 0.0440 2.5971 0.0000
TPPP3 Tppp3 -2.8378 0.0247 -1.1882 0.0044
CCER2 Ccer2 -5.1263 0.0005 -1.4253 0.0000