15 June 2024: Database Analysis
Comprehensive Analysis of Sphingolipid Metabolism-Related Genes in Osteoarthritic Diagnosis and Synovial Immune Dysregulation
Zheng Zhu1DE, Bizhi Tu1B, Run Fang1C, Jun Tong1D, Yulong Liu1F, Rende Ning1A*DOI: 10.12659/MSM.943369
Med Sci Monit 2024; 30:e943369
Abstract
BACKGROUND: Osteoarthritis (OA) is a chronic degenerative disease characterized by synovitis and has been implicated in sphingolipid metabolism disorder. However, the role of sphingolipid metabolism pathway (SMP)-related genes in the occurrence of OA and synovial immune dysregulation remains unclear.
MATERIAL AND METHODS: In this study, we obtained synovium-related databases from GEO (n=40 for both healthy controls and OA) and analyzed the expression levels of SMP-related genes. Using 2 algorithms, we identified hub genes and developed a diagnostic model incorporating these hub genes to predict the occurrence of OA. Subsequently, the hub genes were further validated in peripheral blood samples from OA patients. Additionally, CIBERSORT and MCP-counter analyses were employed to explore the correlation between hub genes and immune dysregulation in OA synovium. WGCNA was used to determine enriched modules in different clusters.
RESULTS: Overall, the expression levels of SMP genes were upregulated in OA synovium. We identified 6 hub genes of SMP and constructed an excellent diagnostic model (AUC=0.976). The expression of re-confirmed hub genes showed associations with immune-related cell infiltration and levels of inflammatory cytokines. Furthermore, we observed heterogeneity in the expression patterns of hub genes across different clusters of OA. Notably, older patients displayed increased susceptibility to elevated levels of pain-related inflammatory cytokines and infiltration of immune cells.
CONCLUSIONS: The SMP-related hub genes have the potential to serve as diagnostic markers for OA patients. Moreover, the 4 hub genes of SMP demonstrate wide participation in immune dysregulation in OA synovium. The activation of different pathways is observed among different populations of patients with OA.
Keywords: biomarkers, Immune System Phenomena, lipid metabolism disorders, Osteoarthritis, Knee, Sphingolipids
Introduction
Osteoarthritis (OA) is a chronic degenerative disease that affects joint function. Its global prevalence rate is estimated to be 6–32%, with the elderly population being most affected [1]. Various risk factors contribute to the development of OA, including aging, excessive physical labor, dyslipidemia, and metabolic syndrome [2,3]. OA is characterized by progressive cartilage degeneration, subchondral sclerosis, synovial inflammation, meniscus degeneration, and osteophyte formation [4–6]. Clinical manifestations of OA include knee dysfunction, stiffness, and pain, which impose physical, psychological, and socio-economic burdens, particularly on the elderly [7]. Previous studies have highlighted the significant role of the synovium in OA development, as joint degeneration occurs following apparent pathological alterations in synovial inflammatory responses [8,9]. Osteoarthritic synovial fibroblasts help support immune cell activation [10] and secrete joint-damaged cytokines such as IL-6, TNF, and MMP3 [11]. The inflammatory response aroused by synovitis triggers local infiltration of macrophages, T cells, and B cells [12], which in turn lead to articular cartilage destruction and accelerated joint degeneration through inflammatory cascades [13,14]. Reducing synovial inflammation can help alleviate joint pain, or even reverse articular structural remodeling [15]. There is an urgent need to explore the underlying cellular mechanism of synovitis and methods of protecting the OA synovium from damaged cytokines in people with OA [16].
There are currently several therapies for OA, including limitation of activity, alleviation of pain, and mitigation of cartilage injury in the early stage, for delaying OA progression [17]. As joint degeneration progresses, many patients eventually receive knee replacement [18]. Thus, adequate diagnosis of OA becomes very important. However, the confirmation of OA depends on imaging data and patient symptoms [19,20], and timely diagnosis of OA is difficult due to lack of specific clinical manifestations. The effective recognition of OA is key to preventing the loss of articular cartilage and OA progression in joints [17, 21], and numerous studies have unveiled early diagnostic biomarkers for OA [22,23], but the clinical utilization of these biomarkers remains a challenge, underscoring the need for further exploration of potential candidate markers for clinical application.
Lipid metabolism disorder has been demonstrated to be important in OA pathogenesis, with significant differences in lipid metabolites observed between the synovial fluid of people with OA versus those without OA [24,25]. Specifically, sphingolipids have been found to be upregulated in OA and play a protective role in inhibiting immune infiltration and attenuating synovitis [26,27]. Therefore, sphingolipid metabolism pathways (SMP) may be involved in regulating synovitis in OA. However, the exact impact of SMP gene expression on synovitis has not been reported, and it remains unclear whether SMP genes can be applied to diagnosis of OA. In this study, we identified 4 hub genes associated with SMP in OA synovium. These hub genes have the potential to serve as diagnostic markers for OA and play a crucial role in immune cell infiltration and inflammatory factors. A deeper understanding of the molecular function of SMP genes in osteoarthritic synovium holds promise for exploring novel diagnostic biomarkers for OA and guiding future research endeavors.
Material and Methods
PREPROCESSING OF RAW DATA SOURCE:
The workflow chart of this study is shown in Figure 1. Gene sets expression (GSE) profiles containing the gene expression information of healthy (CT) and OA synovium (GSE29746, GSE12021, GSE55235, GSE55457) were downloaded from the GEO online database with search term “osteoarthritis AND (synovitis OR synovium)”. Then, the batch effects were removed with the “sva” package [28], and we obtained an integrated dataset including 40 OA and 40 healthy samples. Gene annotation file profiles, which used for converting gene IDs within included transcriptome data into gene symbols, were downloaded from the GPL96 (GSE55235, GSE55457, and GSE12021) and GPL4133 (GSE29746) platforms. The genes of SMP were acquired from Reactome (https://reactome.org/). The demographic information of these patients is shown in Table 1.
IDENTIFICATION OF SMP-RELATED DIFFERENTIAL EXPRESSION GENES (DEGS) IN OA SYNOVIUM:
Using the “limma” package in R software, the differential expression analysis was used to explore the DEGs between OA and healthy synovium in the integrated dataset. Using Venn analysis, the SMP-related DEGs were identified by overlapping the genes of SMP and DEGs. The “heatmap” and “ggplot2” packages in R software were used to draw the heatmap and volcano plots, and the “ggpubr” package was used to create boxplots.
IDENTIFICATION OF SMP-RELATED HUB GENES IN OA SYNOVIUM:
Using the “randomForest” package, we conducted random forest classifiers consisting of 1000 decision trees and then validated the 10-fold cross. The SMP-related genes were analyzed with the random forest function in the “randomForest” package and the “forestplot” package was used to generate a forest plot. The significance of each SMP-related gene was evaluated by importance using random forest analysis as previously described [29], and the top 10 were reserved. SVM recursive feature elimination (SVM-RFE) was employed to discard the relatively low predictive feature in each iteration [30]. Hence, SMP genes were ranked from most to least important, and the top 10 were acquired for further analysis. Then, the top 10 SMP-related genes that were acquired separately by random forest and SVM-RFE analysis were applied to Venn analysis, and the SMP-related hub genes were further identified in the peripheral blood samples of OA (GSE63359).
CONSTRUCTION OF OA PREDICTION MODEL:
Multi-factor logistic regression analysis was used to evaluate the diagnostic model. We calculated the area under the receiver operator characteristic (ROC) curve to predict whether the identified hub genes have diagnostic value in OA. According to the results of multi-factor analysis and clinical information, a predictive nomogram was constructed using the “regplot” and “rms” packages in R software [31]. Calibration plots and decision curve analysis (DCA) were used to assess the applicability of the constructed nomogram.
CORRELATION ANALYSIS BETWEEN THE EXPRESSION OF SMP-RELATED HUB GENES AND THE IMMUNE MICROENVIRONMENT:
Using the “relative” and “absolute” methods available on the CIBERSORT (
CONSENSUS CLUSTERING AND WGCNA ANALYSIS:
“ConsensusClusterPlus” package was utilized to conduct the consensus cluster analysis. The median absolute deviation was used to measure the top 5000 most variable genes applied for samples clustered. The cluster of OA samples and identification of the optimal number of clusters was conducted through unsupervised consensus clustering. Among each 1000 resampling iterations in consensus clustering, we applied ward linkage as an agglomeration method and Euclidean distance as a distance measure to perform hierarchical clustering. The gene sets with similar mRNA transcription models across OA synovium samples were weighted by WGCN analysis [32]. We first normalized the expression of OA synovium genes. Using a power adjacency function, we obtained connection strengths between any 2 genes by converting this matrix into an adjacency matrix. The scale-free topology (SFT) criterion was used to choose the parameter (soft threshold) for adjacency, where the scale-free topology is necessary to construct the network. Based on the SFT criterion recommendation, the optimal threshold parameter value was accepted due to its model-fit saturation of more than 0.85. Then, the optimal topology as soft-threshold power was identified by the “pickSoftThreshold” function. The module was determined using the “tree” method with “deepSplit” to identify modules with no less than 10 genes. Each modular feature was represented by an eigengene generated using the DynamicTreeCut algorithm in the context of WGCNA, as previously described [33]. Default values of WGCNA were consistently applied throughout the process unless stated otherwise. The co-expression modules that enriched for “Biological Process” in Gene Ontology (GO) enrichment analysis were described using Gorilla [34], which is based on pre-inputted gene sets into WGCNA.
GO/KEGG ENRICHMENT ANALYSIS:
The genes that were most obviously different in the WGCNA model were utilized to perform GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis based on the DAVID database (
STATISTICAL ANALYSIS:
All data analysis and processing were performed in R (version 4.2.1). SPSS software (version 23.0) was used for statistical analysis, which began by thoroughly evaluating normality and homoscedasticity to determine the appropriate test method, whether parametric or non-parametric. Normally distributed data underwent the two-tailed unpaired
Results
IDENTIFICATION OF 6 SMP-RELATED HUB GENES IN OA AYNOVIUM:
After removing the batch effects, we performed the differential expression analysis in the integrated database (Figure 2A). We found that the SMP-related genes were expressed at higher levels in the OA group than in the CT group (Figure 2B, P<0.05). The 11 314 DEGs were detected by differential expression analysis between OA and the healthy group. Then, 48 SMP-related DEGs were obtained through overlapping the 11 314 DEGs and 93 SMP genes. A volcano and a box plot presented the level of 48 SMP-related DEGs (Figure 2C). The heatmap presents the expression heterogeneity of 48 SMP-related DEGs between healthy and OA synovium (Figure 2C). The 48 SMP-related DEGs were ranked by SVM-RFE and random forest analysis (Figure 2D). Then, 6 candidate hub SMP-related genes in OA synovium were obtained by intersecting the top 10 genes ranked by SVM-RFE and random forest analysis (Figure 3A).
THE 6 SMP-RELATED HUB GENES CAN SERVE AS DIAGNOSTIC MARKERS FOR OA:
To verify whether the 6 candidate hub genes can be used as clinical diagnostic markers for OA, we performed a binary univariate/multivariable logistic regression analysis for the 6 SMP-related hub genes, and statistically significant variables were included to develop a nomogram line chart for predicting the onset of OA. The nomogram shows that the possibility of OA was determined by the total scores corresponding to the sum of 6 included hub genes (Figure 3B, nomogram for 1 of the random samples). Then, the ROC curve was drawn to verify the reliability of the constructed diagnostic model. The model has a good diagnostic value for OA (Figure 3C, AUC=0.976). The random bootstrap procedure drawn with replacement from the OA samples (Figure 3C, n=100 bootstraps) was applied to assess the reliability of the constructed model. The calibration curves revealed that the model predictions of the nomogram plots were close to those of the ideal model (Figure 3C). The decision curve analysis revealed that the diagnostic model that included the 6 hub genes worked better for OA patients than the model that included only 1 hub gene (Figure 3D). Figure 4 shows the range of AUC, sensitivity, and specificity in the diagnostic model.
FOUR HUB GENES WERE IDENTIFIED IN PERIPHERAL BLOOD SAMPLES IN OA:
To determine whether the diagnostic marker we identified have clinical practicality, we assessed the expression levels of the 6 hub genes. We acquired a dataset comprising transcriptome data from 26 healthy individuals and 46 OA patients’ peripheral blood samples. Upon comparing the expression levels of 4 key genes in the peripheral blood of healthy subjects and those with OA, we observed significant disparities in the expression of these genes (Figure 5). Thus, the following 4 nominated hub genes were considered hub genes of SMP in OA synovium: beta-1,3-N-acetylgalactosaminyltransferase1 (B3GALNT1), sphingosine-1-phosphate phosphatase 1 (SGPP1), sphingomyelin synthase 1 (SGMS1), and sphingosine kinase 1 (SPHK1).
HUB GENES REGULATED THE IMMUNE MICROENVIRONMENT IN OA SYNOVIUM:
CIBERSORT analysis [35] was utilized to deconvolve the gene expression profiles of hub genes. Subsequently, adhering to the operational guidelines outlined on the CIBERSORT website (https://cibersortx.stanford.edu/), an mRNA-based CIBERSORT algorithm incorporating 100 permutations was created to investigate the relationship between the expression levels of the 4 SMP hub genes, immune cell infiltration, and inflammatory factors (Figure 6). Our findings revealed significant associations for B3GALNT1 with 9 out of 22 immune cell infiltrations, 5 out of 9 immune cell populations, and 13 out of 24 inflammatory cytokine levels. Additionally, we observed clear correlations between the expression level of SGMS1 and 7 out of 22 immune cell infiltrations, 6 out of 24 inflammatory cytokine levels, and a positive relationship with B lineage abundance. In the case of SPHK1, we found significant correlations with 5 out of 22 immune cell infiltrations, 9 out of 24 inflammatory cytokine levels, and the abundance of 2 out of 9 immune cell populations. Furthermore, the expression level of SGPP1 was associated with 4 out of 22 immune cell infiltrations, and 4 out of 24 inflammatory cytokine levels, and showed a negative correlation with endothelial cell abundance.
COMPARISON OF CLINIC PARAMETERS AND IMMUNE CHARACTERISTICS AMONG CLUSTERS:
The connection between identified hub genes and clinical factors was explored in R software with “ggpubr” and “ConsensusClusterPlus” packages [36]. We identified the 3 subtypes of OA patents based on the optimal K value (Figure 7A). We found the expression of 4 hub genes among 3 clusters was significantly heterogeneous, as were the clinical factors containing age and sex in OA samples (Figure 7B). There was no significant association between numbers and age, while there was a significant difference in the probability of different age periods among the 3 clusters (Figure 7C, 7D). The expressions of B3GALNT1, SGPP1, CSF3, and IL2 were upregulated, while SPHK1, PDGFA, and IL1A were downregulated in cluster 2 (Figure 8A, 8B). In cluster 2, there were fewer monocytes, but there were more T cells, cytotoxic lymphocytes, myeloid dendritic cells, and endothelial cells (Figure 8C).
CONSTRUCTION OF CO-EXPRESSION PATTERNS:
No samples were removed after performing quality checks among the 40 OA samples. By employing a power function (beta=10; Figure 9A), an adjacency matrix was created from a pairwise gene correlation matrix, initially computed from individual gene residual values. Subsequently, this matrix was converted into a topological overlap matrix, facilitating module assignments through the utilization of a dynamic tree cutting technique. Using the WGCNA algorithm, we identified all co-expression modules with statistical significance using optimal dynamic tree cut (Figure 9B) and hierarchical clustering (Figure 9C). The MEred modules and cluster 2 had the highest absolute correlation values, and the most relevant molecular changes of cluster 2 were detected in MEred modules. Thus, enrichment analysis was performed using genes from the MEred modules (Table 2). We observed a negative association between the phospholipase C-activating G protein-coupled receptor signaling pathway, signaling receptor activator activity, receptor ligand activity, and cluster 1. Conversely, these factors showed a positive association with cluster 2. No statistically significant difference was observed in cluster 3 (Figure 9D).
Discussion
OA is the most common degenerative osteoarthritic disease; it causes irreversible cartilage destruction and synovial inflammation [1]. In recent years, there have been increasing numbers of studies showing that synovitis plays a vital role in the pathogenic process of OA [8]. In this study, a total of 11 314 DEGs were identified between osteoarthritis (OA) and healthy synovium tissue. Among these, 19 genes were found to overlap (15 upregulated, 4 downregulated) and were associated with OA and SMP. Notably, B3GALNT1, SGPP1, SGMS1, and SPHK1 were identified as pivotal hub genes. Utilizing these hub genes, an OA risk prediction model was constructed, demonstrating their potential utility as diagnostic markers for OA. Additionally, the expression levels of these hub genes showed correlations with immune infiltration and inflammatory factors. Through unsupervised clustering and WGCNA analysis, heterogeneity in GO/KEGG enrichment was observed across different clusters.
Previous studies have shown that the triacylceramide 3-β-N-acetylgalactosamine transferase encoded by B3GALNT1 can encode glycosides and P antigens elicit immune responses [37]. Qin et al found that B3GALNT1 gene expression can regulate immune response [38]. In the present study, we found that B3GALNT1 expression shows a strong positive correlation with B lineage, memory B cells, and plasma cells, while it was negatively correlated with the abundance of naive B cells in OA. The result of B3GALNT1 expression may indicate the increased consumption of naive B cells and the activation of B cells in OA synovium. This observation suggests that the upregulation of B3GALNT1 in OA synovium can induce synovitis by enhancing the transformation of naive B cells into memory B cells and plasma cells. Furthermore, B3GALNT1 upregulates the level of IL15 and aggravates bone destruction [39], which is consistent with our results. Thus, downregulated B3GALNT1 expression may improve OA synovitis, but the specific mechanisms need further exploration. For SGPP1, previous studies reported that SGPP1 expression promotes the infiltration of immune cells [40] and the release of pro-inflammation factors [41]. Some pro-inflammation factors, such as IL-6, IL-1A, and IL-1B, can aggravate the development of local inflammation response of OA synovium [42]. In the present study, we found that the level of SGPP1 was negatively correlated with the level of IL1A but was positively correlated with TGFB1, TGFB2, and CSF3. Thus, SGPP1 expression may regulate synovitis development due to releasing different cytokines. Interestingly, SGPP1 increased the proliferation of rheumatoid arthritis (RA) synovial endothelial cells [43], but decreased the number of OA synovial endothelial cells. This inconsistency may indicate the different pathogeneses of OA and RA. The level of SGMS1 involves the release of inflammatory cytokines [44]. In the present study, SGMS1 expression was negatively correlated with IL6. In addition, it has been reported that IL6 accelerates cartilage matrix destruction [42]. Thus, the level of SGMS1 may be protective against OA injury. Previous studies found that SPHK1 regulates various immune cell functions and inflammatory mediator levels, such as IL6 and TGFB2 [45, 46]. In present study, SPHK1 level was positively correlated with IL-6 level and endothelial cell amount, but was negatively correlated with TGFB2. These results may indicate that SPHK1 can decrease the IL6 level and reduce the number of synovial endothelial cells. Meanwhile, consistent with a previous report, TGFB2 can inhibit the inflammatory response and alleviate OA progression [46]. Notable, the large confidence interval observed in the multifactorial analysis of B3GALNT1, SGMS1, and SGPL1 may be attributed to the relatively small sample size included in this study and confounding factors or unmeasured variables that were not accounted for in our analysis.
Early diagnosis plays a crucial role in timely intervention for osteoarthritis. Utilizing transcriptome data from various pathological tissues within joints, numerous studies have identified a range of biomarkers (hub genes) for early OA detection. For example, Xinyue et al [47] employed diverse bioinformatics screening methods to pinpoint key genes showing distinct expression in OA synovial tissue, suggesting their potential as diagnostic markers for OA. Similarly, Hannah et al [48] used single-cell analysis to highlight the significant role of the gene ZEB1, selected from single-cell transcriptome data of joint cartilage and meniscus, in the pathogenesis of OA. Previous research from our team has also indicated the involvement of genes within the arachidonic acid metabolism pathway in OA diagnosis and local inflammation within synovial tissues [23]. In this study, we found that the expression of SMP-related hub genes in synovium can serve as robust early markers for OA and are correlated with local inflammatory cell infiltration and levels of inflammatory factors, suggesting the tight correlation between the SMP-related hub genes and immune disorder in OA synovium. This provides a treatment direction of synovial immune disorders for OA in the future.
Four hub genes involved in sphingolipid metabolism showed heterogeneity among these 3 clusters. We observed that B3GALNT1 and SGPP1 were specifically upregulated in cluster 2, and both of them can be checked in OA patients’ peripheral blood. Thus, the B3GALNT1 and SGPP1 genes may serve as accessible potential biomarkers for OA diagnosis, especially in the elderly population. Utilizing peripheral blood testing for diagnosing OA presents the advantages of convenience and cost-effectiveness, highlighting the importance of gathering additional peripheral blood samples for validation in future studies. The relative score of T cells, endothelial cells, myeloid dendritic cells, and cytotoxic lymphocytes were upregulated, while monocytic lineage was downregulated in cluster 2. This indicates that elderly patients with OA may show a higher level of synovial vascularization and immune cell disorder [49]. We also found that CSF3 and IL2 were upregulated in cluster 2 and PDGFA was downregulated. Downregulation of CSF3 and IL2 significantly reduces the production of inflammatory markers [50,51]. Low levels of CSF3 and IL2 can relieve intractable chronic pain [52–54], and high levels of PDGFA can improve the development of OA synovitis [55]. These results may indicate that downregulation of CSF3 and IL2 and upregulation of PDGFA are potential therapeutic targets for cluster 2 (elderly OA patients). This study offers a direction for personalized treatment of osteoarthritis patients in the future, emphasizing the importance of addressing local immune dysregulation in elderly individuals with OA.
WGCNA analysis showed that MEred module eigengenes were negatively correlated with cluster 1 and cluster 2. MEred module eigengenes mainly involve the “phospholipase C-activating G protein-coupled receptor signaling pathway” in the biological process (BP). The opposite process may indicate higher activation of phospholipase C in cluster 2 compared with cluster 1. Interestingly, the activation of phospholipase C was proved to be harmful in OA, and inhibition of phospholipase C-mediated autophagy may benefit OA cartilage protection [56]. Therefore, targeting the activation of phospholipase C may be a sensitive intervention in cluster 2 patients. Meanwhile, the treatment strategies for inhibiting the activity of phospholipase C in cluster 1 patients may be ineffective. Therefore, the “phospholipase C-activating G protein-coupled receptor signaling pathway” should receive more attention in future research on treatment of elderly OA patients.
The study has limitations that need to be acknowledged. The influence of additional variables affecting the transcription level of inflammatory genes in OA synovium, such as corticosteroids, NSAIDs, and disease duration, remains unexplored. While the study proposed that SMP genes could potentially function as diagnostic markers for OA and interact with immune-mediated inflammation in the synovium, further validation through animal and clinical studies is crucial.
Conclusions
We conducted an investigation into the impact of sphingolipid metabolism in osteoarthritis (OA) synovium and observed significant upregulation of SMP genes. Among these, 4 hub genes (B3GALNT1, SGPP1, SGMS1, and SPHK1) identified from peripheral blood exhibit potential as novel diagnostic markers in OA patients. These hub genes showed significant correlations with inflammatory factors and immune-related cells. Additionally, our findings revealed differential expression of SMP hub genes among different clusters of OA patients, along with variations in the infiltration of immune-related cells and the level of inflammatory cytokines across the 3 clusters. This study suggests that SMP may play a role in regulating the immune microenvironment in OA. The participation of SMP-related hub genes in local synovial tissue immune dysregulation could be a pathogenic mechanism of OA, and targeting these hub genes may offer a potential solution to alleviate joint inflammation. However, further research is needed to fully elucidate the underlying mechanisms involved.
Data Sharing Statement
The datasets used and/or analyzed during the current study are available from the corresponding authors upon reasonable request.
Figures
Figure 1. Work flow chart. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation). Figure 2. Differential expression of genes of the sphingolipid metabolism pathways (SMP). (A) Principal component analysis (PCA) plot shows the batch effect was removed among the 4 obtained datasets. (B) Upregulated expression of sphingolipid metabolism in OA synovium. (C) Volcano plot, box plot, and heat map of differentially expressed SMP genes between OA and healthy (CT) synovium. (D) The rank of SMP-related genes by the random forest analysis and SVM-RFE algorithm. * P<0.05. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 3. Prediction model for OA was constructed by identified hub genes related to sphingolipid metabolism. (A) Six candidate hub genes were identified by random forest and SVM-RFE analysis. (B) Nomogram of 6 candidate hub genes in the diagnosis of OA patients. (C) Receiver operating characteristic (ROC) curve and its verification of repeated sampling, as well as calibration curve indicates the excellence and robustness of model in predicting the OA occurrence. (D) Model evaluation curves: A diagnostic model containing multiple nominated hub genes makes patients more profitable than a single gene. * P<0.05, ** P<0.01. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 4. Distribution of AUC (A), sensitivity (B), and specificity (C) in repeating sample procedure. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 5. Four hub genes identified by peripheral blood-related database (GSE63359). Differentially expression between OA and healthy synovium for B3GALNT1, SGPP1, SGMS1, SPHK1, SGPL1, and HEXB, respectively. * P<0.05, ** P<0.01. Drawn by R (version 4.2.1, JJ Allaire Company), AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and GraphPad Prism (Version 8.0, GraphPad Corporation). Figure 6. The correlation between 4 hub genes expression and immune model. (A) Correlation between hub genes expression and 22 immune-related cell infiltration (CIBERSORT). (B) Correlation between hub genes expression and the abundance of 9 immune-related cells (MCP-counter). (C) Correlation between hub genes expression and the level of 24 inflammatory cytokines. * P<0.05, ** P<0.01, *** P<0.001. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 7. Differences among 3 clusters. (A) Identified the 3 subtypes of OA patents based on the most optimal K value. (B) The heterogeneity of gene expression is related to age and gender. (C) Correlation between age and 3 clusters. (D) Differences in the age-related possibility of OA among 3 clusters. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 8. Correlation between different clustering and hub gene expression and immune module. (A) Differential expression of hub genes in 3 clusters. (B) Differential express level of inflammatory cytokines in 3 clusters. (C) The abundance of 8 immune-related cells differs among their clusters. * P<0.05, ** P<0.01, *** P<0.001. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company). Figure 9. Identification of key modules correlated with hub genes through WGCNA. (A) Analysis of the scale-free fit index (β) and the mean connectivity for various soft-thresholding powers. (B) Clustering dendrogram of genes based on topological overlapping. (C) Heatmap of the correlation between module eigengenes and 3 clusters. (D) Functional enrichment analysis for 3 clusters that most likely enriched in MEred model. Drawn by AI (Adobe Illustrator CC2021, v25.0.0.60 for Mac version, Adobe Corporation) and R (version 4.2.1, JJ Allaire Company).References
1. McDougall C, Hurd K, Barnabe C, Systematic review of rheumatic disease epidemiology in the indigenous populations of Canada, the United States, Australia, and New Zealand: Semin Arthritis Rheum, 2017; 46(5); 675-86
2. Bijlsma JWJ, Berenbaum F, Lafeber FPJG, Osteoarthritis: An update with relevance for clinical practice: Lancet, 2011; 377(9783); 2115-26
3. Vincent TL, Targeting mechanotransduction pathways in osteoarthritis: A focus on the pericellular matrix: Curr Opin Pharmacol, 2013; 13(3); 449-54
4. Chen D, Shen J, Zhao W, Osteoarthritis: Toward a comprehensive understanding of pathological mechanism: Bone Research, 2017; 5; 16044
5. Rousseau J, Garnero P, Biological markers in osteoarthritis: Bone, 2012; 51(2); 265-77
6. Sulzbacher I, Osteoarthritis: Histology and pathogenesis: Wiener Medizinische Wochenschrift, 2012; 163(9–10); 212-19
7. Wu ZY, Du G, Lin YC, Identifying hub genes and immune infiltration of osteoarthritis using comprehensive bioinformatics analysis: J Orthop Surg Res, 2021; 16(1); 630
8. Mathiessen A, Conaghan PG, Synovitis in osteoarthritis: Current understanding with therapeutic implications: Arthritis Res Ther, 2017; 19(1); 18
9. Prieto-Potin I, Largo R, Roman-Blas JA, Herrero-Beaumont G, Walsh DA, Characterization of multinucleated giant cells in synovium and subchondral bone in knee osteoarthritis and rheumatoid arthritis: BMC Musculoskeletal Disorders, 2015; 16(1); 226
10. Xie X, Doody GM, Shuweihdi F, B-cell capacity for expansion and differentiation into plasma cells are altered in osteoarthritis: Osteoarthritis and Cartilage, 2023; 31(9); 1176-88
11. Sanchez-Lopez E, Coras R, Torres A, Synovial inflammation in osteoarthritis progression: Nat Rev Rheumatol, 2022; 18(5); 258-75
12. Wang T, Wang L, Zhang L, Single-cell RNA sequencing in orthopedic research: Bone Res, 2023; 11(1); 10
13. Li N, Xu Q, Liu Q, Leonurine attenuates fibroblast-like synoviocyte-mediated synovial inflammation and joint destruction in rheumatoid arthritis: Rheumatology (Oxford), 2017; 56(8); 1417-27
14. Miao CG, Yang YY, He X, Wnt signaling pathway in rheumatoid arthritis, with special emphasis on the different roles in synovial inflammation and bone remodeling: Cell Signal, 2013; 25(10); 2069-78
15. Mathiessen A, Conaghan PG, Synovitis in osteoarthritis: Current understanding with therapeutic implications: Arthritis Res Ther, 2017; 19(1); 18
16. Conaghan PG, Cook AD, Hamilton JA, Tak PP, Therapeutic options for targeting inflammatory osteoarthritis pain: Nature Rev Rheumatol, 2019; 15(6); 355-63
17. Nelson AE, Allen KD, Golightly YM, A systematic review of recommendations and guidelines for the management of osteoarthritis: The chronic osteoarthritis management initiative of the U.S. bone and joint initiative: Semin Arthritis Rheum, 2014; 43(6); 701-12
18. Roos EM, Arden NK, Strategies for the prevention of knee osteoarthritis: Nat Rev Rheumatol, 2016; 12(2); 92-101
19. Bradke F, Fawcett JW, Spira ME, Assembly of a new growth cone after axotomy: The precursor to axon regeneration: Nat Rev Neurosci, 2012; 13(3); 183-93
20. Zhang Y, Guan Z, Reader B, Autonomic dysreflexia causes chronic immune suppression after spinal cord injury: J Neurosci, 2013; 33(32); 12970-81
21. Liu S, Deng Z, Chen K, Cartilage tissue engineering: From proinflammatory and anti-inflammatory cytokines to osteoarthritis treatments (review): Mol Med Rep, 2022; 25(3); 99
22. Chen Z, Wang W, Zhang Y, Identification of four-gene signature to diagnose osteoarthritis through bioinformatics and machine learning methods: Cytokine, 2023; 169; 156300
23. Tu B, Fang R, Zhu Z, Comprehensive analysis of arachidonic acid metabolism-related genes in diagnosis and synovial immune in osteoarthritis: based on bulk and single-cell RNA sequencing data: Inflamm Res, 2023; 72(5); 955-70
24. Van de Vyver A, Clockaerts S, van de Lest CHA, Synovial fluid fatty acid profiles differ between osteoarthritis and healthy patients: Cartilage, 2020; 11(4); 473-78
25. Gonzalez-Perilli L, Prolo C, Alvarez MN, Arachidonic acid and nitroarachidonic: Effects on NADPH oxidase activity: Adv Exp Med Biol, 2019; 1127; 85-95
26. Kosinska MK, Liebisch G, Lochnit G, Sphingolipids in human synovial fluid – a lipidomic study: PLoS One, 2014; 9(3); e91769
27. Liu SC, Chiu CP, Tsai CH, Soya-cerebroside, an extract of Cordyceps militaris, suppresses monocyte migration and prevents cartilage degradation in inflammatory animal models: Sci Rep, 2017; 7; 43205
28. Chakraborty S, Datta S, Datta S, Surrogate variable analysis using partial least squares (SVA-PLS) in gene expression studies: Bioinformatics, 2012; 28(6); 799-806
29. Wang H, Zhou L, Random survival forest with space extensions for censored data: Artif Intell Med, 2017; 79; 52-61
30. Zhou X, Tuck DP, MSVM-RFE: Extensions of SVM-RFE for multiclass gene selection on DNA microarray data: Bioinformatics, 2007; 23(9); 1106-14
31. Tu B, Jia Y, Qian J, Bioinformatics analysis identified five widely expressed genes associated with prognosis in sarcoma: Int J Gen Med, 2022; 15; 3711-25
32. Zhang B, Horvath S, A general framework for weighted gene co-expression network analysis: Stat Appl Genet Mol Biol, 2005; 4; Article17
33. Tang J, Kong D, Cui Q, Prognostic genes of breast cancer identified by gene co-expression network analysis: Front Oncol, 2018; 8; 374
34. Eden E, Navon R, Steinfeld I, GOrilla: A tool for discovery and visualization of enriched GO terms in ranked gene lists: BMC Bioinformatics, 2009; 10; 48
35. Newman AM, Liu CL, Green MR, Robust enumeration of cell subsets from tissue expression profiles: Nat Methods, 2015; 12(5); 453-57
36. Wilkerson MD, Hayes DN, ConsensusClusterPlus: A class discovery tool with confidence assessments and item tracking: Bioinformatics, 2010; 26(12); 1572-73
37. Westman JS, Benktander J, Storry JR, Identification of the molecular and genetic basis of PX2, a glycosphingolipid blood group antigen lacking on globoside-deficient erythrocytes: J Biol Chem, 2015; 290(30); 18505-18
38. Qin WY, Gan LN, Xia RW, Promoter identification and analysis of key glycosphingolipid biosynthesis-globo series pathway genes in piglets: Genet Mol Res, 2017; 16(2); 16029574
39. Okabe I, Kikuchi T, Mogi M, IL-15 and RANKL play a synergistically important role in osteoclastogenesis: J Cell Biochem, 2017; 118(4); 739-47
40. Nema R, Kumar A, Sphingosine-1-phosphate catabolizing enzymes predict better prognosis in triple-negative breast cancer patients and correlates with tumor-infiltrating immune cells: Front Mol Biosci, 2021; 8; 697922
41. Huang WC, Liang J, Nagahashi M, Sphingosine-1-phosphate phosphatase 2 promotes disruption of mucosal integrity, and contributes to ulcerative colitis in mice and humans: FASEB J, 2016; 30(8); 2945-58
42. Wang T, He C, Pro-inflammatory cytokines: The link between obesity and osteoarthritis: Cytokine Growth Factor Rev, 2018; 44; 38-50
43. Turathum B, Gao EM, Grataitong K, Dysregulated sphingolipid metabolism and autophagy in granulosa cells of women with endometriosis: Front Endocrinol (Lausanne), 2022; 13; 906570
44. Koh EH, Yoon JE, Ko MS, Sphingomyelin synthase 1 mediates hepatocyte pyroptosis to trigger non-alcoholic steatohepatitis: Gut, 2021; 70(10); 1954-64
45. Bu Y, Wu H, Deng R, Wang Y, Therapeutic potential of SphK1 inhibitors based on abnormal expression of SphK1 in inflammatory immune related-diseases: Front Pharmacol, 2021; 12; 733387
46. Tchetina EV, Antoniou J, Tanzer M, Transforming growth factor-beta2 suppresses collagen cleavage in cultured human osteoarthritic cartilage, reduces expression of genes associated with chondrocyte hypertrophy and degradation, and increases prostaglandin E(2) production: Am J Pathol, 2006; 168(1); 131-40
47. Hu X, Ni S, Zhao K, Bioinformatics-led discovery of osteoarthritis biomarkers and inflammatory infiltrates: Front Immunol, 2022; 13; 871008
48. Swahn H, Li K, Duffy T, Senescent cell population with ZEB1 transcription factor as its main regulator promotes osteoarthritis in cartilage and meniscus: Ann Rheum Dis, 2023; 82(3); 403-15
49. Pessler F, Ogdie AR, Mayer CT, Amyloid arthropathy associated with multiple myeloma: Polyarthritis without synovial infiltration of CD20+ or CD38+ cells: Amyloid, 2014; 21(1); 28-34
50. Son ES, Ko UW, Jeong HY, miRNA-6515-5p regulates particulate matter-induced inflammatory responses by targeting CSF3 in human bronchial epithelial cells: Toxicol In Vitro, 2022; 84; 105428
51. Corrigall VM, Arastu M, Khan S, Functional IL-2 receptor beta (CD122) and gamma (CD132) chains are expressed by fibroblast-like synoviocytes: Activation by IL-2 stimulates monocyte chemoattractant protein-1 production: J Immunol, 2001; 166(6); 4141-47
52. Cook AD, Pobjoy J, Steidl S, Granulocyte-macrophage colony-stimulating factor is a key mediator in experimental osteoarthritis pain and disease development: Arthritis Res Ther, 2012; 14(5); R199
53. Eller OC, Glidden N, Knight B, A role for global DNA methylation level and IL2 expression in the transition from acute to chronic low back pain: Front Pain Res (Lausanne), 2021; 2; 744148
54. Fan Z, Fan C, Qi B, Sympathetic nerve-mediated fellow eye pain during sequential cataract surgery by regulating granulocyte colony stimulating factor CSF3: Front Cell Neurosci, 2022; 16; 841733
55. Almasry SM, Soliman HM, El-Tarhouny SA, Platelet rich plasma enhances the immunohistochemical expression of platelet derived growth factor and vascular endothelial growth factor in the synovium of the meniscectomized rat models of osteoarthritis: Ann Anat, 2015; 197; 38-49
56. Chen X, Wang Y, Qu N, PLCgamma1 inhibition-driven autophagy of IL-1beta-treated chondrocyte confers cartilage protection against osteoarthritis, involving AMPK, Erk and Akt: J Cell Mol Med, 2021; 25(3); 1531-45
Figures
Tables
In Press
Clinical Research
High-Resolution Anorectal Manometry and Balloon Expulsion Test Outcomes in Functional Constipation: A Compa...Med Sci Monit In Press; DOI: 10.12659/MSM.944599
Clinical Research
Plaque Removal and Gingival Bleeding Using Biodegradable Toothbrushes: Salvadora persica, Bamboo, and Nylo...Med Sci Monit In Press; DOI: 10.12659/MSM.944469
Clinical Research
Comparative Analysis of Extended Curettage with Plate Fixation and Extended Curettage with Intramedullary N...Med Sci Monit In Press; DOI: 10.12659/MSM.945157
Clinical Research
Clinical Outcomes of Talus Bone Marrow Edema: Comparing Plaster Cast, Crutches, and Walking Boot TreatmentsMed Sci Monit In Press; DOI: 10.12659/MSM.946072
Most Viewed Current Articles
17 Jan 2024 : Review article 6,227,150
Vaccination Guidelines for Pregnant Women: Addressing COVID-19 and the Omicron VariantDOI :10.12659/MSM.942799
Med Sci Monit 2024; 30:e942799
14 Dec 2022 : Clinical Research 1,908,702
Prevalence and Variability of Allergen-Specific Immunoglobulin E in Patients with Elevated Tryptase LevelsDOI :10.12659/MSM.937990
Med Sci Monit 2022; 28:e937990
16 May 2023 : Clinical Research 695,232
Electrophysiological Testing for an Auditory Processing Disorder and Reading Performance in 54 School Stude...DOI :10.12659/MSM.940387
Med Sci Monit 2023; 29:e940387
07 Jan 2022 : Meta-Analysis 260,405
Efficacy and Safety of Light Therapy as a Home Treatment for Motor and Non-Motor Symptoms of Parkinson Dise...DOI :10.12659/MSM.935074
Med Sci Monit 2022; 28:e935074