Apolipoprotein D as a Potential Biomarker and Construction of a Transcriptional Regulatory-Immune Network Associated with Osteoarthritis by Weighted Gene Coexpression Network Analysis

Chondroprotective effects of Apolipoprotein D in knee osteoarthritis mice through the PI3K/AKT/mTOR signaling pathway

BACKGROUND

Because the pathophysiology of osteoarthritis (OA) has not been fully elucidated, targeted treatments are lacking. In this study, we assessed the role and underlying mechanism apolipoprotein D (APOD) on the development of OA.

METHODS

To establish an in vitro OA model, we extracted primary chondrocytes from the cartilage of C57BL/6 mice and stimulated the chondrocytes with IL-1β. After APOD intervention or incubation with an overexpressing plasmid, we detected inflammatory-related markers using RT-qPCR, Western blotting, and ELISA. To detect apoptosis and autophagy-related markers, we used flow cytometry, immunofluorescence, and transmission electron microscopy (TEM). Finally, we measured the level of oxidative stress. We also used RNA-seq to identify the APOD-regulated downstream signaling pathways. We used an in vivo mice OA model of the anterior cruciate ligament transection (ACLT) and administered intra-articular adenovirus overexpressing APOD. To examine cartilage damage severity, we used immunohistochemical analysis (IHC), micro-CT, scanning electron microscopy (SEM), and Safranin O-fast green staining.

RESULTS

Our results showed that APOD inhibited chondrocyte inflammation, degeneration, and apoptosis induced by IL-1β. Additionally, APOD reversed autophagy inhibition and oxidative stress and also blocked activation of the PI3K/AKT/mTOR signaling pathway induced by IL-1β. Finally, overexpression of the APOD gene through adenovirus was sufficient to mitigate OA progression.

CONCLUSIONS

Our findings revealed that APOD had a chondroprotective role in OA progression by the PI3K/AKT/mTOR signaling pathway.

The potential molecular markers of inflammatory response in KOA with AD based on single-cell transcriptome sequencing analysis and identification of ligands by virtual screening

Alzheimer's disease (AD) and osteoarthritis (OA) are both senile degenerative diseases. Clinical studies have found that OA patients have a significantly increased risk of AD in their later life. This study hypothesized that chronic aseptic inflammation might lead to AD in KOA patients. However, current research has not yet clarified the potential mechanism between AD and KOA. Therefore, this study intends to use KOA transcriptional profiling and single-cell sequencing analysis technology to explore the molecular mechanism of KOA affecting AD development, and screen potential molecular biomarkers and drugs for the prediction, diagnosis, and prognosis of AD in KOA patients. It was found that the higher the expression of TXNIP, MMP3, and MMP13, the higher the risk coefficient of AD was. In addition, the AUC of TXNIP, MMP3, and MMP13 were all greater than 0.70, which had good diagnostic significance for AD. Finally, through the virtual screening of core proteins in FDA drugs and molecular dynamics simulation, it was found that compound Cobicistat could be targeted to TXNIP, Itc could be targeted to MMP3, and Isavuconazonium could be targeted to MMP13. To sum up, TXNIP, MMP3, and MMP13 are prospective molecular markers in KOA with AD, which could be used to predict, diagnose, and prognosis.

APOD acts on fibroblast-like synoviocyte and chondrocyte to alleviate the process of osteoarthritis in vitro

The pathogenesis of osteoarthritis (OA) is still unclear, leading to the lack of targeted treatment. We aimed to probe into the effect of apolipoprotein D (APOD), the key gene from our previous study through bioinformatics analysis, on fibroblast-like synoviocyte (FLS) and chondrocytes in vitro to confirm its potential roles on the delay of OA progression. Primary FLS and chondrocytes were extracted from synovium and cartilage of OA patients and stimulated with interleukin 1β (IL-1β) in vitro. After APOD intervention, viability and proliferation of FLS and chondrocytes were detected. Subsequently, the inflammatory factors of the two cells were detected by quantitative reverse-transcription polymerase chain reaction, enzyme-linked immunosorbent assay, and western blot, and the apoptosis and autophagy-related substances were determined at the same time. Finally, the oxidation level in FLS and chondrocytes were detected. APOD reversed the change of gene expression stimulated by IL-1β in FLS and chondrocytes. APOD alleviated the proliferation of FLS while promoted proliferation of chondrocytes, and reduced the expression of inflammatory factors. Moreover, APOD promoted apoptosis of FLS and autography of chondrocytes, while reduced apoptosis of chondrocytes. Finally, decrease level of reactive oxygen species (ROS) in both cells were observed after APOD intervention, as well as the increased expression of antioxidant-related genes. APOD had effects on the proliferation of FLS and chondrocytes through apoptosis and autography as well as the reduction of oxidative stress, delaying the progress of OA.

Risk of metabolic abnormalities in osteoarthritis: a new perspective to understand its pathological mechanisms

Although aging has traditionally been viewed as the most important risk factor for osteoarthritis (OA), an increasing amount of epidemiological evidence has highlighted the association between metabolic abnormalities and OA, particularly in younger individuals. Metabolic abnormalities, such as obesity and type II diabetes, are strongly linked to OA, and they affect both weight-bearing and non-weight-bearing joints, thus suggesting that the pathogenesis of OA is more complicated than the mechanical stress induced by overweight. This review aims to explore the recent advances in research on the relationship between metabolic abnormalities and OA risk, including the impact of abnormal glucose and lipid metabolism, the potential pathogenesis and targeted therapeutic strategies.

Proteomic profiling of human menisci from mild joint degeneration and end-stage osteoarthritis versus healthy controls

Objective

To gain new insight into the molecular changes of the meniscus by comparing the proteome profiles of healthy controls with mild degeneration and end-stage osteoarthritis (OA).

Method

We obtained tissue plugs from lateral and medial menisci of 37 individuals (central part of the posterior horn) classified as healthy (n ​= ​12), mild signs of joint damage (n ​= ​13) and end-stage OA (n ​= ​12). The protein profile was analysed by nano-liquid chromatography-mass spectrometry using data-independent acquisition and quantified by Spectronaut. Linear-mixed effects modelling was applied to extract the between-group comparisons.

Results

A similar protein profile was observed for the mild group as compared to healthy controls while the most different group was end-stage OA mainly for the medial compartment. When a pattern of gradual change in protein levels from healthy to end-stage OA was required, a 42-proteins panel was identified, suggesting a potential role in OA development. The levels of QSOX1 were lower and G6PD higher in the mild group following the proposed protein abundance pattern. Qualitative protein changes suggest lower levels of CYTL1 as a potential biomarker of early joint degradation.

Conclusion

For future targeted proteomic approaches, we propose a candidate panel of 42 proteins based on gradually altered meniscal posterior horn protein abundance patterns associated with joint degradation.

Identification of osteoarthritis-characteristic genes and immunological micro-environment features through bioinformatics and machine learning-based approaches

BACKGROUND

Osteoarthritis (OA) is a multifaceted chronic joint disease characterized by complex mechanisms. It has a detrimental impact on the quality of life for individuals in the middle-aged and elderly population while also imposing a significant socioeconomic burden. At present, there remains a lack of comprehensive understanding regarding the pathophysiology of OA. The objective of this study was to examine the genes, functional pathways, and immune infiltration characteristics associated with the development and advancement of OA.

METHODS

The Gene Expression Omnibus (GEO) database was utilized to acquire gene expression profiles. The R software was employed to conduct the screening of differentially expressed genes (DEGs) and perform enrichment analysis on these genes. The OA-characteristic genes were identified using the Weighted Gene Co-expression Network Analysis (WGCNA) and the Lasso algorithm. In addition, the infiltration levels of immune cells in cartilage were assessed using single-sample gene set enrichment analysis (ssGSEA). Subsequently, a correlation analysis was conducted to examine the relationship between immune cells and the OA-characteristic genes.

RESULTS

A total of 80 DEGs were identified. As determined by functional enrichment, these DEGs were associated with chondrocyte metabolism, apoptosis, and inflammation. Three OA-characteristic genes were identified using WGCNA and the lasso algorithm, and their expression levels were then validated using the verification set. Finally, the analysis of immune cell infiltration revealed that T cells and B cells were primarily associated with OA. In addition, Tspan2, HtrA1 demonstrated a correlation with some of the infiltrating immune cells.

CONCLUSIONS

The findings of an extensive bioinformatics analysis revealed that OA is correlated with a variety of distinct genes, functional pathways, and processes involving immune cell infiltration. The present study has successfully identified characteristic genes and functional pathways that hold potential as biomarkers for guiding drug treatment and facilitating molecular-level research on OA.

Proteomic profiling of extracellular vesicles in synovial fluid and plasma from Oligoarticular Juvenile Idiopathic Arthritis patients reveals novel immunopathogenic biomarkers

Introduction

New early low-invasive biomarkers are demanded for the management of Oligoarticular Juvenile Idiopathic Arthritis (OJIA), the most common chronic pediatric rheumatic disease in Western countries and a leading cause of disability. A deeper understanding of the molecular basis of OJIA pathophysiology is essential for identifying new biomarkers for earlier disease diagnosis and patient stratification and to guide targeted therapeutic intervention. Proteomic profiling of extracellular vesicles (EVs) released in biological fluids has recently emerged as a minimally invasive approach to elucidate adult arthritis pathogenic mechanisms and identify new biomarkers. However, EV-prot expression and potential as biomarkers in OJIA have not been explored. This study represents the first detailed longitudinal characterization of the EV-proteome in OJIA patients.

Methods

Fourty-five OJIA patients were recruited at disease onset and followed up for 24 months, and protein expression profiling was carried out by liquid chromatography-tandem mass spectrometry in EVs isolated from plasma (PL) and synovial fluid (SF) samples.

Results

We first compared the EV-proteome of SF vs paired PL and identified a panel of EV-prots whose expression was significantly deregulated in SF. Interaction network and GO enrichment analyses performed on deregulated EV-prots through STRING database and ShinyGO webserver revealed enrichment in processes related to cartilage/bone metabolism and inflammation, suggesting their role in OJIA pathogenesis and potential value as early molecular indicators of OJIA development. Comparative analysis of the EV-proteome in PL and SF from OJIA patients vs PL from age/gender-matched control children was then carried out. We detected altered expression of a panel of EV-prots able to differentiate new-onset OJIA patients from control children, potentially representing a disease-associated signature measurable at both the systemic and local levels with diagnostic potential. Deregulated EV-prots were significantly associated with biological processes related to innate immunity, antigen processing and presentation, and cytoskeleton organization. Finally, we ran WGCNA on the SF- and PL-derived EV-prot datasets and identified a few EV-prot modules associated with different clinical parameters stratifying OJIA patients in distinct subgroups.

Discussion

These data provide novel mechanistic insights into OJIA pathophysiology and an important contribution in the search of new candidate molecular biomarkers for the disease.

Identification of susceptibility modules and hub genes of osteoarthritis by WGCNA analysis

Osteoarthritis (OA) is a major cause of pain, disability, and social burden in the elderly throughout the world. Although many studies focused on the molecular mechanism of OA, its etiology remains unclear. Therefore, more biomarkers need to be explored to help early diagnosis, clinical outcome measurement, and new therapeutic target development. Our study aimed to retrieve the potential hub genes of osteoarthritis (OA) by weighted gene co-expression network analysis (WGCNA) and assess their clinical utility for predicting OA. Here, we integrated WGCNA to identify novel OA susceptibility modules and hub genes. In this study, we first selected 477 and 834 DEGs in the GSE1919 and the GSE55235 databases, respectively, from the Gene Expression Omnibus (GEO) website. Genes with p-value<0.05 and | log2FC | > 1 were included in our analysis. Then, WGCNA was conducted to build a gene co-expression network, which filtered out the most relevant modules and screened out 23 overlapping WGCNA-derived hub genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses elucidated that these hub genes were associated with cell adhesion molecules pathway, leukocyte activation, and inflammatory response. In addition, we conducted the protein–protein interaction (PPI) network in 23 hub genes, and the top four upregulated hub genes were sorted out (CD4, SELL, ITGB2, and CD52). Moreover, our nomogram model showed good performance in predicting the risk of OA (C-index = 0.76), and this model proved to be efficient in diagnosis by ROC curves (AUC = 0.789). After that, a single-sample gene set enrichment (ssGSEA) analysis was performed to discover immune cell infiltration in OA. Finally, human primary synoviocytes and immunohistochemistry study of synovial tissues confirmed that those candidate genes were significantly upregulated in the OA groups compared with normal groups. We successfully constructed a co-expression network based on WGCNA and found out that OA-associated susceptibility modules and hub genes, which may provide further insight into the development of pre-symptomatic diagnosis, may contribute to understanding the molecular mechanism study of OA risk genes.

A novel serological biomarker are associated with disease severity in patients with osteoarthritis

BACKGROUND AND AIM

Apolipoprotein D (ApoD) is a 25-30 kDa glycoprotein, as a distinct component of the human plasma lipoprotein system. Its known biological functions are mainly related to lipid metabolism. The purpose of this study was to explore the potential role of ApoD concentration in knee osteoarthritis (KOA).

METHODS

This study was performed in a population of 113 KOA subjects and 97 healthy controls. Serum ApoD was detected using the commercial ELISA kit (Cusabio, Wu Han, China, Cat No. CSB-EL001935HU). Radiographic progression was evaluated using Kellgren-Lawrence (KL) classification. The clinical severity of KOA was assessed by visual analog score (VAS), Hospital for special surgery (HSS) knee score disease duration and TNF-α. Receiver-operating characteristic (ROC) curve was performed to test the potential diagnostic value of ApoD in radiographic progression.

RESULTS

The serum ApoD level of patients with KOA was significantly lower than that of healthy controls. The serum ApoD level was negatively correlated with KL grades, VAS score, disease duration, TNF-α concentrations. On the contrary, it was positively correlated with HSS score. However, there was no correlation between this serological indicator and which side was affected. ROC curve analysis indicated that attenuated serum ApoD could serve as an appropriate biomarker concerning the diagnosis of KOA.

CONCLUSIONS

Serum ApoD concentration was correlated with the presence and severity of KOA.

Protein N-glycosylation aberrations and glycoproteomic network alterations in osteoarthritis and osteoarthritis with type 2 diabetes

Whether the relationship between type 2 diabetes mellitus (T2DM) and osteoarthritis (OA) can be solely attributed to the shared risk factors, such as obesity, remains controversial. Several studies have revealed the critical role of abnormal glycosylation in the pathogenesis of OA and T2DM. Therefore, we speculate that T2DM may contribute to the pathogenesis of OA through the intrinsic mechanisms of N-glycosylation aberrations. Using N-glycoproteomics, we compared the changes in N-glycosylated protein abundance in cartilage samples from patients with OA without and with T2DM (DM-OA), and from patients with traumatic joint injury (NC) as controls. We identified 847 N-glycosylation sites corresponding to 729 peptides fragments from 374 proteins. The number of N-glycosylated proteins in the DM-OA group tended to decrease compared with that in the OA and NC groups. We identified 22 upregulated and 1 down-regulated N-glycosylated peptides in the OA group compared to the NC group, while only fibronectin 1 (FN1) at position N1007, cartilage intermediate layer protein 1 (CILP) at N346, and collagen type VI alpha 1 chain (COL6A1) at N804, were also identified in the DM-OA group. Compared to the OA group, the downregulation of secreted protein acidic and rich in cysteine (SPARC) at N116, collagen type VI alpha 1 chain (COL6A2) at N785, and asporin (ASPN) at N282, and the upregulation of complement component C8 alpha chain (C8α) at N437, were the most remarkable alterations in the DM-OA group. The differentially expressed N-glycosylated proteins between the OA and DM-OA groups were mainly located extracellularly and enriched in the KEGG pathways involving PI3K/Akt signaling, focal adhesion, and ECM-receptor interaction. Their predicted protein–protein interactions were also depicted. We were thus able to show the general characteristics of N-glycosylation aberrations in OA and DM-OA. Moreover, the upregulated glycosylated complement C8α in the DM-OA group might augment membrane attack complex activity, thereby exacerbating cartilage destruction. Although further confirmation is required, our hypothesis proposes a possible explanation for the deduction that T2DM is an independent risk factor for OA.