Chujun Peng1
Jinhang Huang1Mingyue Li2Guanru Liu2Lingxian Liu2Jiechun Lin2Weijun Sun2Hongtao Liu2Yonghui Huang2*
Xin Chen2*- 1School of Physics and Optoelectronic Engineering, Guangdong University of Technology, Guangzhou, China
- 2School of Automation, Guangdong University of Technology, Guangzhou, China
Introduction: Periodontitis, a common chronic inflammatory disease, significantly impacted oral health. To provide novel biological indicators for the diagnosis and treatment of periodontitis, we analyzed public microarray datasets to identify biomarkers associated with periodontitis.
Method: The Gene Expression Omnibus (GEO) datasets GSE16134 and GSE106090 were downloaded. We performed differential analysis and robust rank aggregation (RRA) to obtain a list of differential genes. To obtain the core modules and core genes related to periodontitis, we evaluated differential genes through enrichment analysis, correlation analysis, protein-protein interaction (PPI) network and competing endogenous RNA (ceRNA) network analysis. Potential biomarkers for periodontitis were identified through comparative analysis of dual networks (PPI network and ceRNA network). PPI network analysis was performed in STRING. The ceRNA network consisted of RRA differentially expressed messenger RNAs (RRA_DEmRNAs) and RRA differentially expressed long non-coding RNAs (RRA_DElncRNAs), which regulated each other’s expression by sharing microRNA (miRNA) target sites.
Results: RRA_DEmRNAs were significantly enriched in inflammation-related biological processes, osteoblast differentiation, inflammatory response pathways and immunomodulatory pathways. Comparing the core ceRNA module and the core PPI module, C1QA, CENPK, CENPU and BST2 were found to be the common genes of the two core modules, and C1QA was highly correlated with inflammatory functionality. C1QA and BST2 were significantly enriched in immune-regulatory pathways. Meanwhile, LINC01133 played a significant role in regulating the expression of the core genes during the pathogenesis of periodontitis.
Conclusion: The identified biomarkers C1QA, CENPK, CENPU, BST2 and LINC01133 provided valuable insight into periodontitis pathology.
1 Significance
The development of periodontitis involves multiple biological processes. If not treated promptly, it will have a significant impact on an individual’s oral health, causing aggressive periodontitis, which leads to the destruction of periodontal tissues (Carvalho et al., 2018). Hence, we present a bioinformatics-based methodology for the identification of markers associated with periodontitis. This approach entails integrating diverse transcriptomic datasets and conducting a comparative network analysis to delineate functional pathways related to potential markers and differentially expressed genes in periodontitis. Through this methodology, we aim to elucidate the underlying mechanisms of periodontitis and unearth novel biomarkers. These findings are anticipated to contribute to the development of more precise strategies for early diagnosis, treatment, and prevention of periodontitis.
2 Introduction
Periodontitis is a common oral disease characterized by an inflammatory response of the periodontal tissues. The development and progression of periodontitis involve several biological processes, including bacterial infection, inflammatory response, immune regulation and tissue repair. If periodontitis is left untreated, it can ultimately result in tooth loss, significantly impacting an individual’s oral health and overall quality of life. Currently, diagnosing periodontitis relies on methods such as oral examinations, gingival probing and X-rays. Nevertheless, these methods have some limitations. In recent years, transcriptomics has become an essential tool for studying complex diseases (Supplitt et al., 2021). By analyzing the expression of genes, transcriptomics can reveal the molecular mechanisms and biological processes of diseases. For the study of periodontitis, transcriptomics can provide insight into the molecular basis of disease onset and progression (Liu et al., 2022a). However, transcriptomics research also faces some challenges (Jiang et al., 2015). Initially, periodontitis is a complex disease involving multiple biological processes (Cecoro et al., 2020). Transcriptomics data from a single source may not fully reflect the complexity of periodontitis. Secondly, transcriptomics data have high-dimensional characteristics, which poses challenges for data processing and analysis. Therefore, when studying the gene expression patterns related to periodontitis, appropriate methods must be employed to tackle and analyze these high-dimensional data. For instance, one may employ a differential analysis to discern variations in gene expression levels between healthy tissues and those afflicted by periodontitis.
Recently, scientists have initiated utilizing bioinformatics methodologies to discern and characterize biomarkers, striving to enhance their proficiency in recognizing disease-specific biomarkers. For instance, Ji et al. reveal the progression of human osteoarthritis with the help of single-cell RNA-seq analysis, integrating transcriptomics data from multiple biological specimens, such as blood, tissues and cell lines (Ji et al., 2019). Liu et al. identify osteosarcoma metastasis-associated signaling pathways with the help of logistic regression analysis (Liu et al., 2018). Yang et al. explored the relationship between MitoEVs and the immune microenvironment in periodontitis by using machine learning and bioinformatics methods (Yang et al., 2024). Cai et al. found a common pattern of gene expression between obesity and periodontitis by analyzing transcriptomic data and identified five important biomarkers (Cai et al., 2023). Huang et al. used machine learning combined with L1 regularisation and the LIME model interpreter to identify genes associated with periodontitis (Huang et al., 2023). He et al. analyzed the data for differential expression by bioinformatics techniques and screened for lncRNAs associated with periodontitis (He et al., 2023). Liu et al. used quantitative TMT proteomics and transcriptomics analyses to determine the protein expression profiles of patients with periodontitis and constructed nine representative biomarkers using machine learning models (Liu et al., 2022b). This approach can obtain more comprehensive and consistent information, improving the diagnosis and prediction of diseases. However, few studies have reported cases of using integrated multi-source transcriptomics data to identify markers related to oral diseases. For example, Wang et al. identify key markers of gingival tissue and immune cell infiltration studies in periodontitis (Wang Zihui et al., 2021). However, there may be batch effects and inconsistencies in data from different sources. It is necessary to consider some algorithms to remove the batch effects. Therefore, the introduction of RRA analysis in our study contributed to the further removal of batch effects Liu et al. use existing methods to analyze transcriptomics data from multiple sources, so it is possible to identify differential genes associated with pulpitis (Liu et al., 2021). Their approach took into account the batch effects from different sources of transcriptomic data. Therefore, We proposed methods for integrating multi-source transcriptome data to identify periodontitis-related markers by referring to the approach of Liu et al. Meanwhile, we refined the core gene identification process through a network comparative analysis involving both the core PPI network and the ceRNA network. The amalgamation of information from diverse networks contributed to the heightened accuracy and robustness of the identified biomarkers. Besides, we supplemented ceRNA network analysis based on their experiments. Functional enrichment analysis was performed in the core ceRNA network, which helped to uncover the major biological processes and pathways involved in the core ceRNA network. We performed correlation analysis of adjacent genes of core genes and identified lncRNAs with important regulatory roles.
3 Materials and methods
3.1 Workflow of the study
This study summarized the workflow diagram in Figure 1.
Figure 1. Flowchart of the study. The black part represents data preprocessing, the red part represents the integration and identification of differential genes, the brown part represents the processing of post-polymerization differential mRNAs, the green part represents the identification of inflammatory and odontogenic genes, the blue part represents the construction and analysis of the ceRNA network, the purple part represents the identification and analysis of the PPI network, and the yellow part is the final identification of the potential genes related to the mechanisms of periodontitis genesis.
Transcriptome datasets (GSE16134 and GSE106090) were collected from GEO, including proximal gingiva, distal gingiva and medial gingiva. The preprocessing included steps such as batch effect detection and probe remapping, which ensured data standardization between different datasets to eliminate differences due to batch effects. Then, the preprocessed data were subjected to differential analysis. The ranked list of differential genes (p < 0.05) was aggregated using RobustRankAggreg (version 0.6.1) to mitigate potential effects from smaller datasets. By employing an aggregation p-value threshold of less than 0.05 and aggregation ranking, we identified a list of differently expressed genes. Additionally, the differential mRNAs were subjected to functional enrichment analysis to understand their functions. To validate the correlation between the periodontitis genes and specific genes, we extracted inflammation-related genes and odontogenesis-related genes from geneCards. Subsequently, utilizing pre-existing miRNA targeting information, we constructed a ceRNA network comprising differential mRNA and lncRNA. MCODE was utilized to identify the core ceRNA module. Enrichment analysis was performed on the core ceRNA module to determine its associated functional pathways. Then we constructed a PPI network utilizing STRING and identified the core module of this PPI network by employing the MCODE plug-in. The core ceRNA module was compared with the core PPI module to identify potential genes implicated in periodontitis. Subsequently, we undertook a more detailed analysis of the correlation between the potential genes and associated genes within the core ceRNA network to pinpoint lncRNAs that played a crucial regulatory role in periodontitis.
3.2 Data processing
Unless specified otherwise, all data processing and statistical analyses were conducted using the R programming environment (R version 4.2.2). Gene expression profiling of GSE106090 involved screening the expression of lncRNAs and mRNAs in gingival tissues from six patients with peri-implantitis, six patients with periodontitis and six healthy individuals. GSE16134 was obtained from 120 systemically healthy periodontitis patients. The entire samples above were divided into diseased and healthy periodontal inclusion treatments. The collected transcriptomics datasets were tested for batch effects with the help of principal component analysis. The possible batch effects of GSE16134 and GSE106090 were removed from the batch effects with the help of the combat package (version 3.42.0).
3.3 Integration and identification of differential genes
Differential analysis was performed with the help of the limma package (version 3.48.3) and differential mRNAs and differential lncRNAs of the two datasets were extracted. Differential mRNAs and differential lncRNAs were obtained by threshold (p < 0.05,
3.4 Functional enrichment analysis
The differential gene sets were functionally annotated by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The analysis of GO and KEGG was performed by using the clusterProfiler (version 3.18.0) package. The results of GO and KEGG enrichment analysis were subsequently visualized by the ggplot2 package (version 3.3.5) by plotting the bubble plot and the bar plot to show the enrichment of differential genes.
3.5 Correlation analysis
The expression data of differential genes after RRA analysis were extracted from the gene expression matrices of GSE16134 and GSE106090, and then the correlation analyses were performed respectively. The pearson method was chosen here to calculate the correlation coefficients. The correlation matrix was subsequently used to cluster genes with similar expression patterns together. The correlation coefficient matrix was visualized using heatmaps to show the correlation between different genes. Additionally, Inflammation-associated regulatory genes and odontogenesis-associated fundamental genes were discerned in geneCards, scrutinizing the correlation between these genes in both expression patterns.
3.6 Construction and scrutiny of ceRNA networks
The lists of differential mRNAs and differential lncRNAs were extracted from the results of differential analysis. Targets of miRNAs were derived from starBase (https://rnasysu.com/encori/). Based on the results of miRNA-mRNA prediction, relationships between differential mRNAs and miRNAs were screened. Additionally, in light of the miRNA-lncRNA target prediction outcomes, the interplay between the differently expressed lncRNAs and miRNAs was scrutinized. We ascertained mRNA-lncRNA pairs featuring identical miRNAs between differentially expressed mRNAs and differentially expressed lncRNAs. Furthermore, the mRNA-lncRNA pairs possessing overlapping miRNAs were deemed as nodes to construct a triangular relationship. Ultimately, we interlinked the triangular relationship between differentially expressed mRNAs and differentially expressed lncRNAs to form a comprehensive ceRNA network diagram. Subsequently, the ceRNA network identified in the preceding phase was imported into Cytoscape for comprehensive visualization of the ceRNA network. Afterward, the built-in plug-in NetworkAnalyzer of Cytoscape was used to calculate the metrics such as degree (number of connections) and centrality of each node (mRNA, lncRNA), and MCODE was used to screen the core ceRNA module. Subsequently, enrichment analysis was conducted on the pinpointed core ceRNA module to scrutinize their potential functional pathways.
3.7 Construction of PPI network
PPI network was constructed for differential mRNAs with the help of STRING’s built-in interaction relationships. Each gene represented a node, and their interaction relationship represented an edge. Additionally, the edges with a composite score ≥0.4 were selected for PPI network construction. For the constructed network, core module was selected by MCODE (parameters: degree cutoff: two, node score cutoff: 0.2, K-core: two) to identify tightly connected subgraphs in the network as gene modules with functional relevance. In PPI network analysis, Cytoscape’s built-in plug-in NetworkAnalyzer was used to calculate metrics such as degree (number of connections), centrality, etc., for each node (gene) to assess its importance in the network.
3.8 Correlation analysis of mRNAs and lncRNAs directly associated with core genes in the ceRNA core network
In the core ceRNA module, we selected those differential mRNAs and differential lncRNAs. Subsequently, core genes and their direct neighbor genes were extracted from gene expression matrices of GSE16134 and GSE106090. Subsequently, a correlation analysis was conducted. The pearson method was chosen to compute the correlation coefficient. Heat maps were utilized to project the correlation coefficient matrix, providing a visual display of the correlation between core genes and their direct neighbor genes.
4 Result
4.1 Differential mRNAs enriched in cell differentiation and regulation of toxic substances
The detailed depiction of distinct gene integration and identification was illustrated in Table 1. The differential mRNAs and lncRNAs have been systematically identified based on the defined threshold (p < 0.05,
Table 1. The number of genes in differential analysis and RRA analysis.
Figure 2. Results of GO/KEGG functional enrichment analysis of post-polymerisation differential mRNAs. Blueprint representation of (A) GO evaluation and (C) KEGG analysis via its network mapping. The distinctive mRNA GO/KEGG enrichment appraisal results are illustrated utilizing bubble charts and network maps, with the radius of the circles denoting the number of enriched genes. The tint of the circles mirrors the p-value, with all enriched terms deemed significant (p < 0.05). Network maps depict the distribution of genes residing within highly reinforced pathway clusters. (B) Enriched GO terms identified across the categories of biological processes, molecular functions, and cellular components. (D) Bubble chart of enriched pathway clusters identified through the KEGG pathway.
4.2 Periodontitis-related genes were highly correlated with inflammation and odontogenesis
Correlation analysis was a further validation and exploration of correlation between periodontitis genes. Based on the results of correlation analysis, it was known that in GSE16134 (Figure 3A), ICAM2, NCF4, CD79A, FKBP11, FCGR2B, LY96, DSPP, FOXC1, HES1, RASGRP2, PSMB9, TYROBP, C1QA, RUNX3, ITGA3, CARM1, SH3GL1, SPI1, A4GALT, MSC, MGMT and BAX had high correlation, among which DSPP, FOXC1 and BAX was associated with odontogenesis. In GSE106090 (Figure 3B), there were two highly correlated blocks, region1 for C1QA, FOXC1, A4GALT, SPI1, MGMT, RUNX3, CXCL12, COL9A3, TRYOBP, CD79A, ICAM2 and NCF4, region2 for MGST1, LBR, GSTO1, ITGA3, DMBT1, HSPA1A, CRYAB, DSC2, EFNA1, ABCG2, HES1, MX1, SH3GL1, KRT8, ESRP1, PKP1, SLPI, IVL, CARM1, IL36RN and KRT18. FOXC1, CXCL12, KRT8 and PKP1 were correlated with odontogenesis. Comparing the high correlation between the two data sets revealed that ABCG2, LBR, ITGA3, CARM1 and SH3GL1 had high correlation in both GSE16134 high correlation region and GSE106090 region2. C1QA, FOXC1, SPI1, MGMT, CXCL12, COL9A3, TYROBP, CD79A, ICAM2 and NCF4 were highly correlated in both GSE16134 and GSE106090 region1.
Figure 3. Results of correlation analysis of differential mRNAs for inflammation and odontogenic functions. We signifyed the extent of the positive or negative correlation through the color intensity (p < 0.05). Correlation coefficients exhibiting red or blue color respectively signify a positive or negative correlation. (A) The correlation analysis was conducted to examine the relationship between genes associated with odontogenesis and inflammation within the mRNAs that were differentially expressed in GSE16134. Among these, boxed in red are significantly associated with region. (B) The correlation analysis was conducted to examine the relationship between genes associated with odontogenesis and inflammation within the mRNAs that were differentially expressed in GSE106090. The top left corner red triangular area is designated as region 1. The central red large triangular area is designated as region 2. Red areas represent strong correlations, blue areas represent weak correlations, genes named in blackness are inflammation-related genes, and genes named in blue are odontogenesis-related genes.
4.3 Inflammatory, odontogenic genes played hub gene roles in ceRNA networks
The correlation between RRA_DEmRNAs and RRA_DElncRNAs was selected by identical miRNAs to establish a ceRNA network, of which ten mRNAs were associated with inflammation, and two mRNAs were associated with odontogenesis (Figure 4A). The PPI network was extracted by the STRING platform (Supplementary Figure S1A), and core PPI module was identified by using the MCODE plug-in (Supplementary Figure S1B). They revealed a nucleus composed of C1QA, CENPK, CENPU and BST2 as central elements within the complex network. Additionally, C1QA exhibited a correlation with inflammation in geneCards, which implied that CENPK, CENPU and BST2 could feasibly manifest as factors linked to periodontitis. They potentially contributed to regulating the physiological processes associated with periodontitis. The core ceRNA module comprised 28 mRNAs and 21 lnRNAs (Figure 4B). Among these, seven mRNAs were found to be associated with inflammation, while two mRNAs were associated with odontogenesis. Among them, C1QA, CENPK, CENPU and BST2 were also core genes identified in the PPI network. Supplementary Figure S2 demonstrated that 11 lncRNAs were intrinsically linked to C1QA within the core ceRNA module. Additionally, 13 lncRNAs were directly linked to CENPU, and 17 lncRNAs were directly associated with CENPK. BST2 was directly related to 18 lncRNAs. Additionally, among the intersection of the lncRNAs linked to these core genes, there were six lncRNAs (LINC00943, LINC00174, DSCAM-AS1, MAGI1-IT1, MIR4458HG and LINC01133) that were universally linked to them. Notably, LINC00943, LINC00174, DSCAM-AS1 and MAGI1-IT1 were upregulated lncRNAs. Therefore, it was postulated that these six lncRNAs might play key roles in the regulation of gene expression in periodontitis.
Figure 4. Construction of the ceRNA network and core network function mining results. (A) The ceRNA network: circles represent mRNAs, triangles represent lncRNAs, the size of the nodes indicates the degree of sparseness of connection with other nodes, node borders in purple indicate inflammation-related, blue indicates odontogenesis-related, and orange indicates that no correlation has been detected with both genes for the time being. (B) The ceRNA core module: node size and border settings refer to the ceRNA network. For the core ceRNA module, GO function enrichment analysis was performed for (C), (D), (E), and (F). Box plots and mesh plots were generated to visualize the result. GO functional enrichment analysis of core ceRNA modules, box plots, bubble plots, and mesh plots were drawn, and the settings of bubble plots and mesh plots were referred to in Figure 2. In the box plots, the longer the bands are, the more significant the functional enrichment is.
4.4 Core ceRNA module enriched in immunoregulatory and complement activation
GO analysis of the core ceRNA module was performed to further reveal the biological significance of genes within the core ceRNA network. In the GO analysis, activation of immune response (Figures 4C, D) [GO:0002253, p = 1.06e-6, involving genes BAX/BTN3A3/C1QA/C1QB/CFD/CFH/CFI (Figure 4E)] was the most significantly enriched BP in the core ceRNA module. Complement activation (Figures 4C, D) [GO:0006956, p = 1.10e-6, involving genes C1QA/C1QB/CFD/CFH/CFI (Figure 4E)] and humoral immune response (Figures 4C, D) [GO:0006959, p = 8.01e-5, involving genes C1QA/C1QB/CFD/CFH/CFI (Figure 4E)] also significantly enriched. It can be seen that C1QA was significantly enriched in immunomodulation-related pathways. The top ten terms in the BP were summarized in Figure 4E.
4.5 LINC01133 played a role in regulating the expression of core genes in the pathogenesis of periodontitis
The results of differential mRNA and lncRNA correlation analysis of direct association of core genes were performed in Figure 5. In GSE16134 expression profile analysis (Figure 5A), two regions appeared to have a high correlation of expression effect at either end of the diagonal with the upper right expression highly correlated. In the GSE106090 expression profile analysis (Figure 5B), most of the core-associated genes showed significant positive correlations. Then we compared with the higher correlation region in the GSE16134 expression profile, which revealed that C1QA and A4GALT, BST2; BST2 and A4GALT, ABCB7, C1QA; CENPU and ABCA13, ABCB7, ABI1, CENPK, LINC01133; CENPK and ABCA13, ABCB7, ABI1, CENPU, LINC01133 were positively correlated in both expression profiles. This observation strongly suggested that LINC01133 might serve as regulators in regulating the expressions of these core genes implicated in the development of periodontitis. Meanwhile, C1QA-BST2 and CENPU-CENPK were highly correlated with each other in the core genes. It can be hypothesized that CENPK, CENPU, BST2 and LINC01133 were potential periodontitis-associated genes.
Figure 5. Results of correlation analyses of core gene-associated genes. Red hue signifies a positive correlation, whereas blue signals a negative one (p < 0.05). (A) Shows the correlation analysis of mRNAs and lncRNAs associated with the four core genes in the ceRNA core network in the GSE16134 expression profile, and (B) shows the correlation analysis of mRNAs and lncRNAs associated with the four core genes in the ceRNA core network in the GSE106090 expression profile. Among them, region 1 and region 2 within (A) have a high correlation.
5 Discussion
Transcriptomics data analysis has emerged as a promising approach for periodontitis research. It enables the identification of markers associated with periodontitis, providing valuable support for diagnosis and treatment. However, transcriptomics data analysis is challenging due to the high dimension, complexity and diversity of the data. Li et al. point out that although a variety of algorithms have been designed to integrate spatial and single-cell transcriptome data, there are significant differences in how these algorithms work and their scope of application (Li et al., 2022). Spatial transcriptome data are highly non-ideal, including features such as complex data structure, low signal-to-noise ratio, high sparsity, and uneven coverage, which pose challenges for in-depth analysis of the data and parsing of biological information (Li et al., 2024). Therefore, efficient algorithms are necessary for identification and analysis. Additionally, integrating and normalizing different transcriptomics data sources is crucial for subsequent research (Hauser et al., 2017).
Therefore, the paper proposes an integrated approach to identify periodontitis-related markers using transcriptomics data from multiple sources. This method efficiently processes large-scale transcriptomics data and employs advanced analysis algorithms to accurately identify these markers. It also enables the integration of periodontitis-related transcriptomics data from different sources. Furthermore, concerns regarding diagnostic techniques and their impact on analysis quality are addressed by checking for batch effects and ensuring sample homogeneity and between-group variability. Remapped probes retain updated annotation information from various platforms. Therefore, this paper presents a comprehensive approach to identify periodontitis-related markers using transcriptomic data from multiple sources. Meanwhile, considering batch effects for both datasets, we performed batch effect removal and introduce RRA analysis to enhance the robustness of results. In addition, core genes and potential markers were previously determined by network analysis for experimental single network determination. In this study, PPI network and ceRNA network are used to identify core genes and potential markers. This allows for more comprehensive biological information. Moreover, the accuracy and robustness of the identified biomarkers can be enhanced by integrating the information from different networks.
From the GO and KEGG analysis, it can be seen that differential mRNAs were involved in several processes related to inflammation and immunity. Among the GO enrichment pathways, immunomodulation, B-cell activation, cell detoxification, epidermal cell differentiation, developmental regulation and apoptosis showed a more significant enrichment. The medical literature provides a wealth of information on the molecular and immunological mechanisms by which T cells and B cells are involved in the pathogenesis of inflammatory diseases (Gonzales, 2000). Cell detoxification has been shown to be significantly enriched in periodontitis (Lehmann, 2020). Suzuki et al. indicate that periodontitis-related genes are significantly enriched in epidermal cell differentiation (Suzuki et al., 2019). Salmon et al. have suggested that developmental regulation is involved in the pathogenesis of periodontal disease (Salmon et al., 2017). Recent studies have demonstrated the involvement of ER stress in periodontal disease (Jiang et al., 2022). Therefore, we assessed the relevance and differential expression of the expressions of ER proteins and activation of immune response pathway in the GO-enriched pathway. The expressions of ER proteins and activation of immune response pathway have a strong correlation (Supplementary Figure S3A). As can be seen from the box plots of differential expression, the expressions of ER proteins and activation of immune response pathway were significantly different in both inflammatory and normal samples (Supplementary Figure S3B, Supplementary Figure S3C, Supplementary Figure S4D). In KEGG enrichment pathways, osteoblastic cell differentiation, glutathione metabolism, the tonicity and coagulation cascades and ECM-receptor interactions presented more enriched pathways. Comparison of the core PPI module and the core ceRNA module revealed that C1QA, BST2 and TYROBP in the PPI core module were significantly enriched in the immunomodulatory and B-cell activation pathways. Willems et al. clearly indicate that the C1Q family (C1QA) is associated with immunoregulatory pathways and autoimmune diseases (Willems, 2021). Alvarez et al. found that BST2 genes associated with antiviral defense, interferon signaling and Toll-like receptor signaling were significantly upregulated in the OPM of VEH/SIV (Alvarez et al., 2020). Huo et al. found that TYROBP rich in complement, inflammatory response, interferon γ response, and TNF-α signaling via NF-κB (Huo and Wang, 2021). The core ceRNA module genes C1QA, BAX, CFD, CFH and CFI were significantly enriched in the immune regulation pathway. Ma et al. noted that CFD + fibroblasts show high expression of chemokines similar to iCAF in some types of tumors (Ma et al., 2023). Duan et al. found that CFH and CFI were associated with immunity and characteristic reflection of periodontitis (Duan et al., 2023). BST2 was significantly enriched in the B-cell activation pathway. Moreover, according to geneCards, C1QA, TYROBP and CFI were found to be associated with inflammation, while BAX was associated with the process of odontogenesis. Additionally, the results of the correlation analysis indicated a strong correlation between C1QA and TYROBP within the two datasets. In summary, it can be inferred that changes in differential mRNA expression were closely associated with inflammation and immune-related biological processes. These findings may contribute to an in-depth understanding of the molecular mechanisms of inflammation and immune regulation.
LncRNAs play a pivotal role in periodontitis and there is growing evidence that lncRNAs have diagnostic value. Wu et al. found that inflammation in diabetes-associated periodontitis can be attenuated by activating the CTBP1-AS2/miR-155/SIRT1 axis (Ng et al., 2024). Xia et al. found that long-chain non-coding RNA PVT1 can be involved in pulpitis pathogenesis by regulating miR-128-3p (Xia et al., 2022). We constructed a ceRNA network aiming to discover the regulatory relationship between differential mRNAs and differential lncRNAs. Meanwhile, we found ten inflammation-associated genes and two odontogenesis-associated genes within the geneCards. Additionally, C1QA, CENPU, CENPK and BST2 were also present in the core PPI module, and C1QA was identified to be highly correlated with inflammatory genes within geneCards from the core ceRNA network. We also found 11 lncRNAs related to the regulation of C1QA, 13 lncRNAs linked to CENPU, 13 lncRNAs linked to CENPK, 17 lncRNAs linked to CENPK and 18 lncRNAs linked to BST2. Six lncRNAs were found to be co-linked to the above four genes, namely, LINC00943, LINC00174, DSCAM-AS1, MAGI1-IT1, MIR4458HG and LINC01133, suggesting that these six lncRNAs may play a role in the regulation of gene expression in periodontitis. Meng et al. found that LINC00943 could attenuate MPP + -induced neuronal injury through the RAB3IP axis in SK-N-SH cells (Meng et al., 2021). Su et al. found that LINC00174 could attenuate cardiac muscle injury through p53-mediated autophagy and apoptosis (Su et al., 2021). Maimeti et al. found that LINC00174 as an immune regulator may have a regulatory role in low-grade gliomas (Maimaiti et al., 2021). Ning et al. found that DSCAM-AS1 can accelerate cell value-addition and migration in osteosarcoma through GPRC5A signaling (Ning and Bai, 2021). Wang et al. found that MAGI1-IT1 has a regulatory role in controlling the value-addition of gastric cancer (Wang et al., 2021). Sun et al. found that LINC01133 also has a regulatory effect on the value-added of gastric cancer (Sun et al., 2022). Zeng et al. found that CENPK has the potential to serve as a predictive marker gene for clinical prognosis and personalized immunotherapy in cancer patients (Zeng et al., 2021). Zhou et al. found that CENPU was a key gene in the development of LUAD, closely associated with the infiltration of various immune cells (Zhou et al., 2021). Shan et al. found that BST2 contributes to the promotion of metastasis, invasion and proliferation of oral squamous cell carcinoma (Shan et al., 2023). Through the above core gene-associated lncRNAs, it has been shown that the above lncRNAs were related to the value-added and differentiation regulation of cancer cells, and LINC00174 is related to immune regulation. Meanwhile, in the enrichment analysis of the core ceRNA network, it was found that the immune-regulatory pathway was the most significant enriched pathway, and C1QA was one of the genes of the pathway. Therefore, we suspected that LINC00943, LINC00174, DSCAM-AS1, MAGI1-IT1, MIR4458HG and LINC01133 may have similar regulatory roles in the proliferation and differentiation of stromal and neutrophil cells as those in the proliferation and differentiation of cancer cells. They may also have a role in the process of immune regulation.
In the enrichment analysis of core ceRNA network, C1QA was significantly enriched in immune activation, complement activation, humoral immune response and leukocyte-mediated immune pathway. Hajishengallis et al. collation of an exposition of how to intervene in periodontal disease mechanisms using complement dependence (Hajishengallis et al., 2019). In this study (Supplementary Figure S4A), ficolin-1 and C1QA have some correlation. Supplementary Figure S4B, Supplementary Figure S4C and Supplementary Figure S4D pointed out that the two have significant differences. BST2 was significantly enriched in leukocyte-mediated immune pathway. The enriched pathways in the PPI network were also immune-regulation related, which suggested to a certain extent that C1QA played a certain role in immune regulation in periodontitis, and BST2 was a more likely potential periodontitis regulator. Meanwhile, CENPU and CENPK also existed in the core PPI network and were associated with cancer cell value-added and differentiation-related processes. Thus, it was speculated that they were also potential regulators in diagnosing periodontitis.
In the core ceRNA network, LINC01133 were identified as positively correlated with coregulators. It implied that LINC01133 may play a role in regulating coregulators' expression in cells. The positive correlation indicated that the expression levels of LINC01133 were consistent with the trend of change in the expression levels of regulatory genes. This may imply that LINC01133 may be involved in the regulatory network of core genes by interacting with core genes, which in turn affects the expression levels of core genes. This positive correlation may provide a new explanatory mechanism that LINC01133 may regulate the expression of core genes by sharing miRNA binding sites with core genes as ceRNAs. Differential lncRNAs increase the expression of core genes by adsorbing miRNAs and blocking the inhibitory effect of miRNAs on core genes. Therefore, the identification of some differential lncRNAs positively correlated with core factor correlations may help to further understand the role of these differential lncRNAs in regulating core factor expression and cellular functions.
In the PPI network, we chose to polymerize the post-differential mRNAs for the study of the interactions network, once again reducing the problem of heterogeneity present in the data from different platforms. The core PPI network was identified within geneCards as containing seven genes related to inflammation, in which C1QA, BST2, CENPU and CENPK were simultaneously present in the core ceRNA module. So it can be inferred that C1QA played an essential regulatory role in periodontitis. Jahanimoghadam et al. describe the interactions of common DEGs by constructing a protein interaction network, in which C1QA is one of the core ceRNA regulators. They find that IF135, MX1, SPI1 and IF144L associated with it are in the periodontitis core PPI network, and C1QB in the periodontitis core ceRNA network (Jahanimoghadam et al., 2022). Wang et al. construct the PPI network of differential genes when studying the changes in gene expression profiles in the dorsal horn of the spinal cord after sciatic nerve injury, and C1QA plays a more central role in this network. Meanwhile, C1QA has a strong interaction with TYROBP and C1QB, while TYROBP also plays a certain regulatory role (Wang et al., 2017). C1QA and C1QB are found to have a high interplay relationship in the bioinformatics analysis studies on the regulatory role of inflammatory genes in dwarfism diseases by Yuan et al. and vascular dementia molecules by Shu et al. (Yuan et al., 2021; Shu et al., 2022). In our study, C1QB exists in the core ceRNA network in periodontitis differential mRNAs and differential lncRNAs, which can be seen in the existing studies indicating that C1QB and C1QA have a strong reciprocal relationship. It can be speculated that there is also a potential reciprocal relationship between C1QB and C1QA in the regulation of periodontitis genes. BST2, CENPU and CENPK are the potential undiscovered genes that are related to the regulation of periodontitis.
Integrating transcriptomics data from multiple sources can reduce bias and error, improving data reliability. It enables the identification of more differentially expressed genes related to periodontitis, enhancing our understanding of its pathogenesis. It allows for comprehensive bioinformatics analysis, including gene function annotation, pathway analysis and protein interaction networks, deepening our understanding of periodontitis-related markers. At the same time, the deletion of the six dental implant samples did not affect the final conclusion. However, limitations exist due to potential differences in sample processing and sequencing platforms, which may affect data consistency.
In comparison with other studies, our research highlights four major points of distinction. Firstly, we consider that multi-source transcriptomic data may introduce batch effects due to factors such as laboratory and sequencing platform. Therefore, the detection and removal of batch effects are carried out before the differential analysis, in order to increase the consistency of the data and reduce the possibility of false positives, and to improve the accuracy of the differential analysis. Secondly, two-net control analysis is used to identify potential markers of periodontitis molecules to help obtain a more comprehensive understanding of gene regulatory networks. The accuracy and robustness of identifying biomarkers can be enhanced by integrating the information of the PPI network and the ceRNA network. Moreover, the functional enrichment analysis of the core ceRNA network is added, which compares with the enrichment analysis of differential genes, and find that the two are interlinked in immune regulation. It can also help us to further determine whether the genes in the core ceRNA network are involved in specific signaling pathways. This can help to further reveal the key pathways that may be affected in the development of periodontitis. Finally, we performed a correlation analysis of the core adjacent genes in the core ceRNA network and identify some of the genes highly positively associated with the core genes. Potential cooperative relationships and therapeutic targets may exist for highly positively correlated genes. This helps to provide a deeper understanding of the molecular mechanisms and disease development of periodontitis. Although there are limitations due to potential differences in sample processing and sequencing platforms, overall it does not affect the identification and identification of potential markers of periodontitis.
In conclusion, this study identified potential markers in the diagnostic process of periodontitis, and analyzed the functional pathways and interactions of the core modules. These results provided candidates for molecular diagnosis.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
CP: Data curation, Formal Analysis, Software, Visualization, Writing – original draft, Writing – review and editing. JH: Software, Writing – review and editing. ML: Writing – original draft. GL: Validation, Writing – review and editing. LL: Formal Analysis, Writing – review and editing. JL: Visualization, Writing – review and editing. WS: Supervision, Writing – review and editing. HL: Supervision, Writing – review and editing. YH: Supervision, Writing – review and editing. XC: Methodology, Supervision, Validation, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by grants from National Natural Science Foundation of China (Grant Nos. 62003094 and 82003615), and Science and Technology Projects in Guangzhou (Grant No. 202201010147).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1398582/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | PPI network construction and MCODE mining of core networks. (A) The PPI network: the PPI network was constructed for the mRNAs with differences after polymerization, the node size indicates the degree of the sparseness of the association with other nodes, the node border color is purple for inflammation, the border color is blue for odontogenesis, and the border color is orange for the time being, no association was found with the two types of physiological processes. (B) The PPI core network: Significant PPI modules delineated by MCODE. Node attributes such as color, dimension, and bordering styles are altered by the parameters stipulated for the PPI network.
SUPPLEMENTARY FIGURE S2 | Results of the identification of genes directly associated with core genes. (A) Represents an inseparable gene cluster directly linked to C1QA, (B) signifies another such gene cluster linked to CENPU, while (C) denotes yet another about CENPK. Lastly, (D) indicates a distinct gene cluster directly related to BST2. Within this figure, circles signify mRNA molecules whereas triangles denote lncRNAs.
SUPPLEMENTARY FIGURE S3 | The correlation assessment of the expressions of ER proteins and activation of immune response pathway. The correlation between the expressions of ER proteins and activation of immune response pathway in this study was assessed using normalized data of GSE16134 and GSE106090. (A) The scatter plot of correlation, the abscissa is the mean of expression of genes within the activation of immune response pathway in GSE16134 and GSE106090. The ordinate is the mean expression of genes within the ER pathway. R is the correlation coefficient while the P-values reflect the significance of the correlation. (B) “ai_ average” is the mean of expression of genes within the activation of immune response pathway in GSE16134 and GSE106090. “er_ average” is 16 the mean expression of genes within the ER pathway in GSE16134 and GSE106090. The expression of two pathways in the (C) normal samples and (D) inflammatory samples. The P-values were calculated by t-test.
SUPPLEMENTARY FIGURE S4 | The correlation assessment of C1QA and ficolin-1. The assessment method was consistent with the Supplementary Figure S3. (A) The scatter plot of correlation with the abscissa representing the mean of expression of genes within the ficolin-1 in GSE16134 and GSE106090. The ordinate is the mean expression of C1QA in GSE16134 and GSE106090. (B) “f1_ average” is the mean of expression of genes within ficolin-1 in GSE16134 and GSE106090. “C1QA” is the expression of C1QA in GSE16134 and GSE106090. The expression of C1QA and ficolin-1 in (C) normal samples and (D) inflammatory samples.
References
Alvarez, X., Sestak, K., Byrareddy, S. N., and Mohan, M. (2020). Long term Delta-9-tetrahydrocannabinol administration inhibits proinflammatory responses in minor salivary glands of chronically simian immunodeficieny virus infected rhesus macaques. Viruses 12 (7), 713. doi:10.3390/v12070713
Cai, Y., Zuo, X., Zuo, Y., Wu, S., Pang, W., Ma, K., et al. (2023). Transcriptomic analysis reveals shared gene signatures and molecular mechanisms between obesity and periodontitis. Front. Immunol. 14, 1101854. doi:10.3389/fimmu.2023.1101854
Carvalho, C. V., Saraiva, L., Bauer, F. P. F., Kimura, R. Y., Souto, M. L. S., Bernardo, C. C., et al. (2018). Orthodontic treatment in patients with aggressive periodontitis. Am. J. Orthod. Dentofac. Orthop. 153 (4), 550–557. doi:10.1016/j.ajodo.2017.08.018
Cecoro, G., Annunziata, M., Iuorio, M. T., Nastri, L., and Guida, L. (2020). Periodontitis, low-grade inflammation and systemic health: a scoping review. Med. Kaunas. 56 (6), 272. doi:10.3390/medicina56060272
Duan, A., Zhang, Y., and Yuan, G. (2023). Screening of feature genes related to immune and inflammatory responses in periodontitis. BMC Oral Health 23 (1), 234. doi:10.1186/s12903-023-02925-z
Gonzales, J. R. (2000). T and B-cell subsets in periodontitis. Periodontology 69 (1), 181–200. doi:10.1111/prd.12090
Hajishengallis, G., Kajikawa, T., Hajishengallis, E., Maekawa, T., Reis, E. S., Mastellos, D. C., et al. (2019). Complement-Dependent mechanisms and interventions in periodontal disease. Front. Immunol. 12 (10), 406. doi:10.3389/fimmu.2019.00406
Hauser, A. S., Attwood, M. M., Rask-Andersen, M., Schiöth, H. B., and Gloriam, D. E. (2017). Trends in GPCR drug discovery: new agents, targets and indications. Nat. Rev. Drug Discov. 16 (12), 829–842. doi:10.1038/nrd.2017.178
He, J., Zheng, Z., Li, S., Liao, C., and Li, Y. (2023). Identification and assessment of differentially expressed necroptosis long non-coding RNAs associated with periodontitis in human. BMC Oral Health 23 (1), 632. doi:10.1186/s12903-023-03308-0
Huang, Q., Wang, Z., and Han, J. (2023). Identification of key gene associated with periodontitis and prediction of therapeutic drugs using machine learning in combination with LIME model explainer. MEDS Clin. Med. 4 (5), 133–142. doi:10.23977/medsc.2023.040518
Huo, T., and Wang, Z. (2021). Comprehensive analysis to identify key genes involved in advanced atherosclerosis. Dis. Markers 2021, 2021–2025. doi:10.1155/2021/4026604
Jahanimoghadam, A., Abdolahzadeh, H., Rad, N. K., and Zahiri, J. (2022). Discovering common pathogenic mechanisms of COVID-19 and Parkinson disease: an integrated bioinformatics analysis. J. Mol. Neurosci. 72 (11), 2326–2337. doi:10.1007/s12031-022-02068-w
Ji, Q., Zheng, Y., Zhang, G., Hu, Y., Fan, X., Hou, Y., et al. (2019). Single-cell RNA-seq analysis reveals the progression of human osteoarthritis. Ann. Rheum. Dis. 78 (1), 100–110. doi:10.1136/annrheumdis-2017-212863
Jiang, Z., Zhou, X., Li, R., Michal, J. J., Zhang, S., Dodson, M. V., et al. (2015). Whole transcriptome analysis with sequencing: methods, challenges and potential solutions. Cell Mol. Life Sci. 72 (18), 3425–3439. doi:10.1007/s00018-015-1934-y
Jiang, , Li, Z., and Zhu, G. (2022). The role of endoplasmic reticulum stress in the pathophysiology of periodontal disease. J. Periodontal Res. 57 (5), 915–932. doi:10.1111/jre.13031
Lehmann, K. (2020). Oral porphyromonas gingivalis infection induces epigenetic changes that promote persistence of cardiovascular disease risk.
Li, B., Zhang, W., Guo, C., Xu, H., Li, L., Fang, M., et al. (2022). Benchmarking spatial and single-cell transcriptomics integration methods for transcript distribution prediction and cell type deconvolution. Nat. methods 19 (6), 662–670. doi:10.1038/s41592-022-01480-9
Li, R., Chen, X., and Yang, X. (2024). Navigating the landscapes of spatial transcriptomics: how computational methods guide the way. Wiley Interdiscip. Rev. RNA 15 (2), e1839. doi:10.1002/wrna.1839
Liu, L., Wang, T., Huang, D., and Song, D. (2021). Comprehensive analysis of differentially expressed genes in clinically diagnosed irreversible pulpitis by multiplatform data integration using a robust rank aggregation approach. J. Endod. 47 (9), 1365–1375. doi:10.1016/j.joen.2021.07.007
Liu, W., Qiu, W., Huang, Z., Zhang, K., Wu, K., Deng, K., et al. (2022a). Identification of nine signature proteins involved in periodontitis by integrated analysis of TMT proteomics and transcriptomics. Front. Immunol. 13, 963123. doi:10.3389/fimmu.2022.963123
Liu, W., Qiu, W., Huang, Z., Zhang, K., Wu, K., Deng, K., et al. (2022b). Identification of nine signature proteins involved in periodontitis by integrated analysis of TMT proteomics and transcriptomics. Front. Immunol. 13, 963123. doi:10.3389/fimmu.2022.963123
Liu, Y., Sun, W., Ma, X., Hao, Y., Liu, G., Hu, X., et al. (2018). Logistic regression analysis for the identification of the metastasis-associated signaling pathways of osteosarcoma. Int. J. Mol. Med. 41 (3), 1233–1244. doi:10.3892/ijmm.2018.3360
Ma, C., Yang, C., Peng, A., Sun, T., Ji, X., Mi, J., et al. (2023). Pan-cancer spatially resolved single-cell analysis reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment. Mol. Cancer 22 (1), 170. doi:10.1186/s12943-023-01876-x
Maimaiti, A., Jiang, L., Wang, X., Shi, X., Pei, Y., Hao, Y., et al. (2021). Identification and validation of an individualized prognostic signature of lower-grade glioma based on nine immune related long non-coding RNA. Clin. Neurol. Neurosurg. 201, 106464. doi:10.1016/j.clineuro.2020.106464
Meng, C., Gao, J., Ma, Q., Sun, Q., and Qiao, T. (2021). LINC00943 knockdown attenuates MPP+-induced neuronal damage via miR-15b-5p/RAB3IP axis in SK-N-SH cells. Neurol. Res. 43 (3), 181–190. doi:10.1080/01616412.2020.1834290
Ng, M. Y., Yu, C. C., Chen, S. H., Liao, Y. W., and Lin, T. (2024). Er: YAG laser alleviates inflammaging in diabetes-associated periodontitis via activation CTBP1-AS2/miR-155/SIRT1 Axis. Int. J. Mol. Sci. 25 (4), 2116. doi:10.3390/ijms25042116
Ning, Y., and Bai, Z. (2021). DSCAM-AS1 accelerates cell proliferation and migration in osteosarcoma through miR-186-5p/GPRC5A signaling. Cancer Biomark. 30 (1), 29–39. doi:10.3233/CBM-190703
Salmon, C. R., Giorgetti, A. P. O., Leme, A. F. P., Domingues, R. R., Kolli, T. N., Foster, B. L., et al. (2017). Microproteome of dentoalveolar tissues. Bone 101, 219–229. doi:10.1016/j.bone.2017.05.014
Shan, F., Shen, S., Wang, X., and Chen, G. (2023). BST2 regulated by the transcription factor STAT1 can promote metastasis, invasion and proliferation of oral squamous cell carcinoma via the AKT/ERK1/2 signaling pathway. Int. J. Oncol. 62 (4), 54. doi:10.3892/ijo.2023.5502
Shu, J., Wei, W., and Zhang, L. (2022). Identification of molecular signatures and candidate drugs in vascular dementia by bioinformatics analyses. Front. Mol. Neurosci. 15, 751044. doi:10.3389/fnmol.2022.751044
Su, Q., Lv, X. W., Xu, Y. L., Cai, R. P., Dai, R. X., Yang, X. H., et al. (2021). Exosomal LINC00174 derived from vascular endothelial cells attenuates myocardial I/R injury via p53-mediated autophagy and apoptosis. Mol. Ther. Nucleic Acids 23, 1304–1322. doi:10.1016/j.omtn.2021.02.005
Sun, Y., Tian, Y., He, J., Zhang, G., Zhao, R., et al. (2022). Linc01133 contributes to gastric cancer growth by enhancing YES1-dependent YAP1 nuclear translocation via sponging miR-145-5p. Cell Death Dis. 13 (1), 51. doi:10.1038/s41419-022-04500-w
Supplitt, S., Karpinski, P., Sasiadek, M., and Laczmanska, I. (2021). Current achievements and applications of transcriptomics in personalized cancer medicine. Int. J. Mol. Sci. 22 (3), 1422. doi:10.3390/ijms22031422
Suzuki, A., Horie, T., and Numabe, Y. (2019). Investigation of molecular biomarker candidates for diagnosis and prognosis of chronic periodontitis by bioinformatics analysis of pooled microarray gene expression datasets in Gene Expression Omnibus (GEO). BMC Oral Health 19, 52–15. doi:10.1186/s12903-019-0738-0
Wang, J., Ma, S. H., Tao, R., Xia, L. J., Liu, L., and Jiang, Y. H. (2017). Gene expression profile changes in rat dorsal horn after sciatic nerve injury. Neurol. Res. 39 (2), 176–182. doi:10.1080/01616412.2016.1273590
Wang, Q., Gu, M., Zhuang, Y., and Chen, J. (2021b). The long noncoding RNA MAGI1-IT1 regulates the miR-302d-3p/IGF1 Axis to control gastric cancer cell proliferation. Cancer Manag. Res. 13, 2959–2967. doi:10.2147/CMAR.S305108
Wang, Z., Gao, H., Feifei, Mo, Guangjie, T., and Wang, Y. (2021a). Identification of key biomarkers and immune infiltration in human gingival tissue of periodontitis by bioinformatics analysis. J. Oral Sci. Res. 37 (4), 304–309. doi:10.13701/j.cnki.kqyxyj.2021.04.007
Willems, E. (2021). Clinical (glyco) proteomics: from pre-clinical discovery to translational diagnostics in infectious disease and complement deficiency.
Xia, L., Wang, J., Qi, Y., Fei, Y., and Wang, D. (2022). Long non-coding RNA PVT1 is involved in the pathological mechanism of pulpitis by regulating miR-128-3p. Oral Health Prev. Dent. 20 (1), 263–270. doi:10.3290/j.ohpd.b3147193
Yang, H., Zhao, A., Chen, Y., Cheng, T., Zhou, J., and Li, Z. (2024). Exploring the potential link between MitoEVs and the immune microenvironment of periodontitis based on machine learning and bioinformatics methods. BMC Oral Health 24 (1), 169. doi:10.1186/s12903-024-03912-8
Yuan, J., Du, Z., Wu, Z., Yang, Y., Cheng, X., Liu, X., et al. (2021). A novel diagnostic predictive model for idiopathic short stature in children. Front. Endocrinol. (Lausanne) 12, 721812. doi:10.3389/fendo.2021.721812
Zeng, H., Shen, Y., Hirachan, S., Bhandari, A., and Zhang, X. (2021). Pan-cancer investigation of CENPK gene: clinical significance and oncogenic immunology. Am. J. Transl. Res. 13 (12), 13336–13355.
Keywords: periodontitis, biomarkers, network analysis, integration, public microarrary datasets
Citation: Peng C, Huang J, Li M, Liu G, Liu L, Lin J, Sun W, Liu H, Huang Y and Chen X (2024) Uncovering periodontitis-associated markers through the aggregation of transcriptomics information from diverse sources. Front. Genet. 15:1398582. doi: 10.3389/fgene.2024.1398582
Received: 19 March 2024; Accepted: 10 May 2024;
Published: 11 June 2024.
Edited by:
Junichi Iwata, University of Texas Health Science Center at Houston, United StatesReviewed by:
Anshul Tiwari, Vanderbilt University, United StatesTakao Fukuda, Kyushu University, Japan
Copyright © 2024 Peng, Huang, Li, Liu, Liu, Lin, Sun, Liu, Huang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin Chen, xinchen@gdut.edu.cn; Yonghui Huang, msh_huizi@gdut.edu.cn