RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage

A Damaging COL6A3 Variant Alters the MIR31HG-Regulated Response of Chondrocytes in Neocartilage Organoids to Hyperphysiologic Mechanical Loading

The pericellular matrix (PCM), with its hallmark proteins collagen type VI (COLVI) and fibronectin (FN), surrounds chondrocytes and is critical in transducing the biomechanical cues. To identify genetic variants that change protein function, exome sequencing is performed in a patient with symptomatic OA at multiple joint sites. A predicted damaging variant in COL6A3 is identified and introduced by CRISPR-Cas9 genome engineering in two established human induced pluripotent stem cell-derived in-vitro neocartilage organoid models. The downstream effects of the COL6A3 variant on the chondrocyte phenotypic state are studied by a multi-omics (mRNA and lncRNA) approach in interaction with hyper-physiological mechanical loading conditions. The damaging variant in COL6A3 results in significantly lower binding between the PCM proteins COLVI and FN and provokes an osteoarthritic chondrocyte state. By subsequently exposing the neocartilage organoids to hyperphysiological mechanical stress, it is demonstrated that the COL6A3 variant in chondrocytes abolishes the characteristic inflammatory signaling response after mechanical loading with PTGS2, PECAM1, and ADAMTS5, as central genes. Finally, by integrating epigenetic regulation, the lncRNA MIR31HG is identified as key regulator of the characteristic inflammatory signaling response to mechanical loading.

Evolution and advancements in genomics and epigenomics in OA research: How far we have come

OBJECTIVE

Osteoarthritis (OA) is the most prevalent musculoskeletal disease affecting articulating joint tissues, resulting in local and systemic changes that contribute to increased pain and reduced function. Diverse technological advancements have culminated in the advent of high throughput "omic" technologies, enabling identification of comprehensive changes in molecular mediators associated with the disease. Amongst these technologies, genomics and epigenomics - including methylomics and miRNomics, have emerged as important tools to aid our biological understanding of disease.

DESIGN

In this narrative review, we selected articles discussing advancements and applications of these technologies to OA biology and pathology. We discuss how genomics, deoxyribonucleic acid (DNA) methylomics, and miRNomics have uncovered disease-related molecular markers in the local and systemic tissues or fluids of OA patients.

RESULTS

Genomics investigations into the genetic links of OA, including using genome-wide association studies, have evolved to identify 100+ genetic susceptibility markers of OA. Epigenomic investigations of gene methylation status have identified the importance of methylation to OA-related catabolic gene expression. Furthermore, miRNomic studies have identified key microRNA signatures in various tissues and fluids related to OA disease.

CONCLUSIONS

Sharing of standardized, well-annotated omic datasets in curated repositories will be key to enhancing statistical power to detect smaller and targetable changes in the biological signatures underlying OA pathogenesis. Additionally, continued technological developments and analysis methods, including using computational molecular and regulatory networks, are likely to facilitate improved detection of disease-relevant targets, in-turn, supporting precision medicine approaches and new treatment strategies for OA.

miRNA-seq analysis of high glucose induced osteoblasts provides insight into the mechanism underlying diabetic osteoporosis

The present study aims to explore the etiology of Diabetic osteoporosis (DOP), a chronic complication associated with diabetes mellitus. Specifically, the research seeks to identify potential miRNA biomarkers of DOP and investigated role in regulating osteoblasts. To achieve this, an animal model of DOP was established through the administration of a high-sugar and high-fat diet, and then injection of streptozotocin. Bone microarchitecture and histopathology analysis were analyzed. Rat calvarial osteoblasts (ROBs) were stimulated with high glucose (HG). MiRNA profiles of the stimulated osteoblasts were compared to control osteoblasts using sequencing. Proliferation and mineralization abilities were assessed using MTT assay, alkaline phosphatase, and alizarin red staining. Expression levels of OGN, Runx2, and ALP were determined through qRT-PCR and Western blot. MiRNA-sequencing results revealed increased miRNA-702-5p levels. Luciferase reporter gene was utilized to study the correlation between miR-702-5p and OGN. High glucose impaired cell proliferation and mineralization in vitro by inhibiting OGN, Runx2, and ALP expressions. Interference with miR-702-5p decreased OGN, Runx2, and ALP levels, which were restored by OGN overexpression. Additionally, downregulation of OGN and Runx2 in DOP rat femurs was confirmed. Therefore, the miRNA-702-5p/OGN/Runx2 signaling axis may play a role in DOP, and could be diagnostic biomarker and therapeutic target for not only DOP but also other forms of osteoporosis.

Gene expression and immune infiltration analysis comparing lesioned and preserved subchondral bone in osteoarthritis

Background

Osteoarthritis (OA) is a degenerative disease requiring additional research. This study compared gene expression and immune infiltration between lesioned and preserved subchondral bone. The results were validated using multiple tissue datasets and experiments.

Methods

Differentially expressed genes (DEGs) between the lesioned and preserved tibial plateaus of OA patients were identified in the GSE51588 dataset. Moreover, functional annotation and protein-protein interaction (PPI) network analyses were performed on the lesioned and preserved sides to explore potential therapeutic targets in OA subchondral bones. In addition, multiple tissues were used to screen coexpressed genes, and the expression levels of identified candidate DEGs in OA were measured by quantitative real-time polymerase chain reaction. Finally, an immune infiltration analysis was conducted.

Results

A total of 1,010 DEGs were identified, 423 upregulated and 587 downregulated. The biological process (BP) terms enriched in the upregulated genes included "skeletal system development", "sister chromatid cohesion", and "ossification". Pathways were enriched in "Wnt signaling pathway" and "proteoglycans in cancer". The BP terms enriched in the downregulated genes included "inflammatory response", "xenobiotic metabolic process", and "positive regulation of inflammatory response". The enriched pathways included "neuroactive ligand-receptor interaction" and "AMP-activated protein kinase signaling". JUN, tumor necrosis factor α, and interleukin-1β were the hub genes in the PPI network. Collagen XI A1 and leucine-rich repeat-containing 15 were screened from multiple datasets and experimentally validated. Immune infiltration analyses showed fewer infiltrating adipocytes and endothelial cells in the lesioned versus preserved samples.

Conclusion

Our findings provide valuable information for future studies on the pathogenic mechanism of OA and potential therapeutic and diagnostic targets.

Hyper-physiologic mechanical cues, as an osteoarthritis disease-relevant environmental perturbation, cause a critical shift in set points of methylation at transcriptionally active CpG sites in neo-cartilage organoids

BACKGROUND

Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (Guilak in Best Pract Res Clin Rheumatol 25:815-823, 2011). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state.

RESULTS

We showed that hyper-physiological loading evokes consistent changes in CpGs associated with expression changes (ML-tCpGs) in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology (OA-tCpGs), we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity.

CONCLUSION

Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCpGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.

Thyroid hormone induces ossification and terminal maturation in a preserved OA cartilage biomimetic model

OBJECTIVE

To characterize aspects of triiodothyronine (T3) induced chondrocyte terminal maturation within the molecular osteoarthritis pathophysiology using the previously established T3 human ex vivo osteochondral explant model.

DESIGNS

RNA-sequencing was performed on explant cartilage obtained from OA patients (n = 8), that was cultured ex vivo with or without T3 (10 ng/ml), and main findings were validated using RT-qPCR in an independent sample set (n = 22). Enrichment analysis was used for functional clustering and comparisons with available OA patient RNA-sequencing and GWAS datasets were used to establish relevance for OA pathophysiology by linking to OA patient genomic profiles.

RESULTS

Besides the upregulation of known hypertrophic genes EPAS1 and ANKH, T3 treatment resulted in differential expression of 247 genes with main pathways linked to extracellular matrix and ossification. CCDC80, CDON, ANKH and ATOH8 were among the genes found to consistently mark early, ongoing and terminal maturational OA processes in patients. Furthermore, among the 37 OA risk genes that were significantly affected in cartilage by T3 were COL12A1, TNC, SPARC and PAPPA.

CONCLUSIONS

RNA-sequencing results show that metabolic activation and recuperation of growth plate morphology are induced by T3 in OA chondrocytes, indicating terminal maturation is accelerated. The molecular mechanisms involved in hypertrophy were linked to all stages of OA pathophysiology and will be used to validate disease models for drug testing.

Identified senescence endotypes in aged cartilage are reflected in the blood metabolome

Heterogeneous accumulation of senescent cells expressing the senescence-associated secretory phenotype (SASP) affects tissue homeostasis which leads to diseases, such as osteoarthritis (OA). In this study, we set out to characterize heterogeneity of cellular senescence within aged articular cartilage and explored the presence of corresponding metabolic profiles in blood that could function as representative biomarkers. Hereto, we set out to perform cluster analyses, using a gene-set of 131 senescence genes (N = 57) in a previously established RNA sequencing dataset of aged articular cartilage and a generated metabolic dataset in overlapping blood samples. Using unsupervised hierarchical clustering and pathway analysis, we identified two robust cellular senescent endotypes. Endotype-1 was enriched for cell proliferating pathways, expressing forkhead box protein O4 (FOXO4), RB transcriptional corepressor like 2 (RBL2), and cyclin-dependent kinase inhibitor 1B (CDKN1B); the FOXO mediated cell cycle was identified as possible target for endotype-1 patients. Endotype-2 showed enriched inflammation-associated pathways, expressed by interleukin 6 (IL6), matrix metallopeptidase (MMP)1/3, and vascular endothelial growth factor (VEGF)C and SASP pathways were identified as possible targets for endotype-2 patients. Notably, plasma-based metabolic profiles in overlapping blood samples (N = 21) showed two corresponding metabolic clusters in blood. These non-invasive metabolic profiles could function as biomarkers for patient-tailored targeting of senescence in OA.

MiR-29a-3p mediates phosphatase and tensin homolog and inhibits osteoarthritis progression

Despite substantial progress in clinical trials of osteoarthritis (OA) gene therapy, the prevalence of OA is still on the rise. MiRNAs have a potential biomarker and therapeutic target for OA. OA cartilage and chondrosarcoma cells were studied to determine the role of miR-29a-3p and PTEN. OA cartilage and human chondrosarcoma cells (SW1353) were obtained. miR-29a-3p and PTEN signature expression was determined by RT-qPCR. The binding relationship between miR-29a-3p and PTEN was investigated by dual-luciferase reporter gene and western blot assay. TUNEL, immunohistochemistry, CCK-8, and flow cytometry were utilized to determine the proliferation and apoptosis of SW1353 cells. This study indicated downregulation of miR-29a-3p expression and upregulation of PTEN expression in human OA primary chondrocytes or OA tissue samples, compared with the normal cartilage cells or tissues. PTEN expression was negatively correlated with miR-29a-3p expression, and miR-29a-3p targeted PTEN mechanistically. miR-29a-3p reduced SW1353 cell activity and proliferation and promoted cell apoptosis. However, the aforementioned effects could be reversed by downregulating PTEN. miR-29a-3p can stimulate chondrocyte proliferation and inhibit apoptosis by inhibiting PTEN expression.

Articular cartilage repair biomaterials: strategies and applications

Articular cartilage injury is a frequent worldwide disease, while effective treatment is urgently needed. Due to lack of blood vessels and nerves, the ability of cartilage to self-repair is limited. Despite the availability of various clinical treatments, unfavorable prognoses and complications remain prevalent. However, the advent of tissue engineering and regenerative medicine has generated considerable interests in using biomaterials for articular cartilage repair. Nevertheless, there remains a notable scarcity of comprehensive reviews that provide an in-depth exploration of the various strategies and applications. Herein, we present an overview of the primary biomaterials and bioactive substances from the tissue engineering perspective to repair articular cartilage. The strategies include regeneration, substitution, and immunization. We comprehensively delineate the influence of mechanically supportive scaffolds on cellular behavior, shedding light on emerging scaffold technologies, including stimuli-responsive smart scaffolds, 3D-printed scaffolds, and cartilage bionic scaffolds. Biologically active substances, including bioactive factors, stem cells, extracellular vesicles (EVs), and cartilage organoids, are elucidated for their roles in regulating the activity of chondrocytes. Furthermore, the composite bioactive scaffolds produced industrially to put into clinical use, are also explicitly presented. This review offers innovative solutions for treating articular cartilage ailments and emphasizes the potential of biomaterials for articular cartilage repair in clinical translation.

Mechanical stress and inflammation have opposite effects on Wnt signaling in human chondrocytes

Dysregulation of Wingless and Int-1 (Wnt) signaling has been strongly associated with development and progression of osteoarthritis (OA). Here, we set out to investigate the independent effects of either mechanical stress (MS) or inflammation on Wnt signaling in human neocartilage pellets, and to relate this Wnt signaling to OA pathophysiology. OA synovium-conditioned media (OAS-CM) was collected after incubating synovium from human end-stage OA joints for 24 h in medium. Cytokine levels in the OAS-CM were determined with a multiplex immunoassay (Luminex). Human neocartilage pellets were exposed to 20% MS, 2% OAS-CM or 1 ng/mL Interleukin-1β (IL-1β). Effects on expression levels of Wnt signaling members were determined by reverse transcription-quantitative polymerase chain reaction. Additionally, the expression of these members in articular cartilage from human OA joints was analyzed in association with joint space narrowing (JSN) and osteophyte scores. Protein levels of IL-1β, IL-6, IL-8, IL-10, tumor necrosis factor α, and granulocyte-macrophage colony-stimulating factor positively correlated with each other. MS increased noncanonical WNT5A and FOS expression. In contrast, these genes were downregulated upon stimulation with OAS-CM or IL-1β. Furthermore, Wnt inhibitors DKK1 and FRZB decreased in response to OAS-CM or IL-1β exposure. Finally, expression of WNT5A in OA articular cartilage was associated with increased JSN scores, but not osteophyte scores. Our results demonstrate that MS and inflammatory stimuli have opposite effects on canonical and noncanonical Wnt signaling in human neocartilage. Considering the extent to which MS and inflammation contribute to OA in individual patients, we hypothesize that targeting specific Wnt pathways offers a more effective, individualized approach.

Exploring miRNA-target gene pair detection in disease with coRmiT

A wide range of approaches can be used to detect micro RNA (miRNA)-target gene pairs (mTPs) from expression data, differing in the ways the gene and miRNA expression profiles are calculated, combined and correlated. However, there is no clear consensus on which is the best approach across all datasets. Here, we have implemented multiple strategies and applied them to three distinct rare disease datasets that comprise smallRNA-Seq and RNA-Seq data obtained from the same samples, obtaining mTPs related to the disease pathology. All datasets were preprocessed using a standardized, freely available computational workflow, DEG_workflow. This workflow includes coRmiT, a method to compare multiple strategies for mTP detection. We used it to investigate the overlap of the detected mTPs with predicted and validated mTPs from 11 different databases. Results show that there is no clear best strategy for mTP detection applicable to all situations. We therefore propose the integration of the results of the different strategies by selecting the one with the highest odds ratio for each miRNA, as the optimal way to integrate the results. We applied this selection-integration method to the datasets and showed it to be robust to changes in the predicted and validated mTP databases. Our findings have important implications for miRNA analysis. coRmiT is implemented as part of the ExpHunterSuite Bioconductor package available from https://bioconductor.org/packages/ExpHunterSuite.

MicroRNA-34a-5p promotes the progression of osteoarthritis secondary to developmental dysplasia of the hip by restraining SESN2-induced autophagy

Osteoarthritis (OA), a late-stage complication of developmental dysplasia of the hip (DDH), is a key factor leading to further degeneration of joint function. Studies have shown that Sestrin2 (SESN2) is a positive regulator in protecting articular cartilage from degradation. However, the regulatory effects of SESN2 on DDH-OA and its upstream regulators remain obscure. Here, we first identified that the expression of SESN2 significantly decreased in the cartilage of DDH-OA samples, with an expression trend negatively correlated with OA severity. Using RNA sequencing, we identified that the upregulation of miR-34a-5p may be an important factor for the decrease in SESN2 expression. Further exploring the regulation mechanism of miR-34a-5p/SESN2 is of great significance for understanding the mechanism of DDH occurrence and development. Mechanistically, we showed that miR-34a-5p could significantly inhibit the expression of SESN2, thereby promoting the activity of the mTOR signaling pathway. We also found that miR-34a-5p significantly inhibited SESN2-induced autophagy, thereby suppressing the proliferation and migration of chondrocytes. We further validated that knocking down miR-34a-5p in vivo resulted in a significant increase in SESN2 expression and autophagy activity in DDH-OA cartilage. Our study suggests that miR-34a-5p is a negative regulator of DDH-OA, and may provide a new target for the prevention of DDH-OA.

Identification and Validation of PTGS2 Gene as an Oxidative Stress-Related Biomarker for Arteriovenous Fistula Failure

(1) Background: Arteriovenous fistulas (AVFs) are the preferred site for hemodialysis. Unfortunately, approximately 60% of patients suffer from AVF failure within one year. Oxidative stress plays an important role in the occurrence and development of AVF. However, the underlying mechanisms remain unclear. Therefore, specific oxidative stress-related biomarkers are urgently needed for the diagnosis and treatment of AVF failure. (2) Methods: Bioinformatics analysis was carried out on dataset GSE119296 to screen for PTGS2 as a candidate gene related to oxidative stress and to verify the expression level and diagnostic efficacy of PTGS2 in clinical patients. The effects of NS398, a PTGS2 inhibitor, on hemodynamics, smooth muscle cell proliferation, migration, and oxidative stress were evaluated in a mouse AVF model. (3) Results: Based on 83 oxidative stress-related differentially expressed genes, we identified the important pathways related to oxidative stress. PTGS2 may have diagnostic and therapeutic efficacy for AVF failure. We further confirmed this finding using clinical specimens and validation datasets. The animal experiments illustrated that NS398 administration could reduce neointimal area (average decrease: 49%) and improve peak velocity (average increase: 53%). (4) Conclusions: Our study identified PTGS2 as an important oxidative stress-related biomarker for AVF failure. Targeting PTGS2 reduced oxidative stress and improved hemodynamics in an AVF mouse model.

Hyper-physiologic mechanical cues, as an osteoarthritis disease relevant environmental perturbation, cause a critical shift in set-points of methylation at transcriptionally active CpG sites in neo-cartilage organoids

Background

Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (1). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state.

Results

We showed that hyper-physiological loading evokes consistent changes in ML-tCpGs associated with expression changes in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology, we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity.

Conclusion

Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCPGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.

Development of microRNA-based therapeutics for central nervous system diseases

MicroRNA (miRNA)-mediated gene silencing is a method of RNA interference in which a miRNA binds to messenger RNA sequences and regulates target gene expression. MiRNA-based therapeutics have shown promise in treating a variety of central nervous system diseases, as verified by results from diverse preclinical model organisms. Over the last decade, several miRNA-based therapeutics have entered clinical trials for various kinds of diseases, such as tumors, infections, and inherited diseases. However, such clinical trials for central nervous system diseases are scarce, and many central nervous system diseases, including hemorrhagic stroke, ischemic stroke, traumatic brain injury, intractable epilepsy, and Alzheimer's disease, lack effective treatment. Considering its effectiveness for central nervous system diseases in preclinical experiments, microRNA-based intervention may serve as a promising treatment for these kinds of diseases. This paper reviews basic principles and recent progress of miRNA-based therapeutics and summarizes general procedures to develop such therapeutics for treating central nervous system diseases. Then, the current obstacles in drug development are discussed. This review also provides a new perspective on possible solutions to these obstacles in the future.

PSD95 as a New Potential Therapeutic Target of Osteoarthritis: A Study of the Identification of Hub Genes through Self-Contrast Model

Osteoarthritis (OA) is a worldwide joint disease. However, the precise mechanism causing OA remains unclear. Our primary aim was to identify vital biomarkers associated with the mechano-inflammatory aspect of OA, providing potential diagnostic and therapeutic targets for OA. Thirty OA patients who underwent total knee arthroplasty were recruited, and cartilage samples were obtained from both the lateral tibial plateau (LTP) and medial tibial plateau (MTP). GO and KEGG enrichment analyses were performed, and the protein-protein interaction (PPI) assessment was conducted for hub genes. The effect of PSD95 inhibition on cartilage degeneration was also conducted and analyzed. A total of 1247 upregulated and 244 downregulated DEGs were identified. Significant differences were observed between MTP and LTP in mechanical stress-related genes and activated sensory neurons based on a self-contrast model of human knee OA. Cluster analysis identified DLG4 as the hub gene. Cyclic loading stress increased PSD95 (encoded by DLG4) expression in LTP cartilage, and PSD95 inhibitors could alleviate OA progression. This study suggests that inhibiting PSD95 could be a potential therapeutic strategy for preventing articular cartilage degradation.

Microstructural and cellular characterisation of the subchondral trabecular bone in human knee and hip osteoarthritis using synchrotron tomography

OBJECTIVE

It is unclear if different factors influence osteoarthritis (OA) progression and degenerative changes characterising OA disease in hip and knee. We investigated the difference between hip OA and knee OA at the subchondral bone (SCB) tissue and cellular level, relative to the degree of cartilage degeneration.

DESIGN

Bone samples were collected from 11 patients (aged 70.4 ± 10.7years) undergoing knee arthroplasty and 8 patients (aged 62.3 ± 13.4years) undergoing hip arthroplasty surgery. Trabecular bone microstructure, osteocyte-lacunar network, and bone matrix vascularity were evaluated using synchrotron micro-CT imaging. Additionally, osteocyte density, viability, and connectivity were determined histologically.

RESULTS

The associations between severe cartilage degeneration and increase of bone volume fraction (%) [- 8.7, 95% CI (-14.1, -3.4)], trabecular number (#/mm) [- 1.5, 95% CI (-0.8, -2.3)], osteocyte lacunar density (#/mm3) [4714.9; 95% CI (2079.1, 7350.6)] and decrease of trabecular separation (mm) [- 0.07, 95% CI (0.02, 0.1)] were found in both knee and hip OA. When compared to knee OA, hip OA was characterised by larger (µm3) but less spheric osteocyte lacunae [47.3; 95% CI (11.2, 83.4), - 0.04; 95% CI (-0.06, -0.02), respectively], lower vascular canal density (#/mm3) [- 22.8; 95% CI (-35.4, -10.3)], lower osteocyte cell density (#/mm2) [- 84.2; 95% CI (-102.5, -67.4)], and less senescent (#/mm2) but more apoptotic osteocytes (%) [- 2.4; 95% CI (-3.6, -1.2), 24.9; 95% CI (17.7, 32.1)], respectively.

CONCLUSION

SCB from hip OA and knee OA exhibits different characteristics at the tissue and cellular levels, suggesting different mechanisms of OA progression in different joints.

Long Non-coding RNA SNHG7 Suppresses Inflammation and Apoptosis of Chondrocytes Through Inactivating of p38 MAPK Signaling Pathway in Osteoarthritis

This study aims to explore the molecular mechanism of LncRNA SNHG7 in Osteoarthritis (OA). Cartilage tissues of OA patients or patients with trauma or amputation were collected. Compared to normal cartilage tissues, SNHG7 was downregulated while miR-324-3p was upregulated in cartilage tissues of OA patients. IL-1β was used to induce damage to chondrocytes and treatment with IL-1β reduced SNHG7 expression in OA chondrocytes. In IL-1β-treated OA chondrocytes, SNHG7 overexpression reduced the levels of TNF-α and IL-6, inhibited cell apoptosis, and increased cell viability. Additionally, the luciferase reporter assay proved that SNHG7 upregulated dual-specificity phosphatase 1 (DUSP1) by sponging miR-324-3p, thereby inactivating the p38 MAPK signaling pathway by regulating the miR-324-3p/DUSP1 axis. Anisomycin (a p38 MAPK activator) enhanced OA chondrocytes inflammation, promoted cell apoptosis, and reduced cell viability; however, this was reversed by SNHG7 overexpression. This study demonstrates that the SNHG7/miR-324-3p/DUSP1 axis suppresses OA chondrocytes inflammation and apoptosis by inhibiting the p38 MAPK signaling pathway. Thus, this study indicates that SNHG7 is a novel target for OA treatment.

Osteoclast-Derived Exosomal miR-5134-5p Interferes with Alveolar Bone Homeostasis by Targeting the JAK2/STAT3 Axis

Background

In chronic periodontitis, exosomes transport various informative substances between osteoclasts and osteoblasts in alveolar bone. Herein, we aimed to investigate the effect of exosomal micro-ribonucleic acid (miRNA/miR)-5134-5p derived from osteoclasts on osteoblastic proliferation and differentiation and the development of periodontitis in vivo and in vitro.

Methods

The effects of OC-Exos on the proliferation and differentiation of osteoblasts were identified by Real-time quantitative reverse polymerase chain reaction (qRT-PCR), Western blot(WB), alkaline phosphatase(ALP) staining, etc. Exosomal miRNA expression was analyzed by sequencing. The sites of miRNA action were predicted through TargetScan and tested by double luciferase assay. After transfecting miR-5134-5p mimic/inhibitor into osteoblasts, we measured the proliferation and differentiation of osteoblasts by ALP staining and WB, etc. Furthermore, OC-Exos were injected into the gingival sulcus at the ligation site. Inflammation was observed by Hematoxylin-eosin (H&E) staining, the expression of inflammatory factors were detected by qRT-PCR, the resorption of alveolar bone was observed by Micro CT.

Results

Osteoblastic proliferation and differentiation were negatively regulated by OC-Exos in vitro. miRNA sequencing analysis revealed that miR-5134-5p expression was significantly elevated in OC-Exos, which also increased in osteoblasts following OC-Exo intervention. The dual-luciferase assay revealed that miR-5134-5p and Janus kinase 2 (JAK2) had binding sites. miR-5134-5p-mimics could upregulate miR-5134-5p expression in osteoblasts while downregulating Runt-related transcription factor 2(Runx2), phosphorylated-JAK2 (p-JAK2), and phosphorylated-signal transducer and activator of transcription 3 (p-STAT3) expression and inhibited osteogenic differentiation. However, miR-5134-5p-inhibitor had the opposite effect. In vivo, the OC-Exo group demonstrated morphological disruption of periodontal tissue, massive inflammatory cell infiltration, upregulation of inflammatory factors mRNA expression, a significant decrease in BV/TV, and an increase in the cementoenamel junction and alveolar bone crest distance.

Conclusion

Osteoclast-derived exosomal miR-5134-5p inhibits osteoblastic proliferation and differentiation via the JAK2/STAT3 pathway. OC-Exos exacerbate periodontal tissue inflammation and accelerate alveolar bone resorption in mice with experimental periodontitis.

Integrated Analysis of Transcriptome Changes in Osteoarthritis: Gene Expression, Pathways and Alternative Splicing

OBJECTIVE

Osteoarthritis (OA) is the most prevalent joint disease characterized by the degeneration of articular cartilage and the remodeling of its underlying bones, resulting in pain and loss of function in the knees and hips. As far as we know, no curative treatments are available except for the joint replacement. The precise molecular mechanisms which are involved in the degradation of cartilage matrix and development of osteoarthritis are still unclear.

DESIGN

By analyzing RNA-seq data, we found the molecular changes at the transcriptome level such as alternative splicing, gene expression, and molecular pathways in OA knees cartilage.

RESULTS

Expression analysis have identified 457 differential expressed genes including 266 up-regulated genes such as TNFSF15, ST6GALNAC5, TGFBI, ASPM, and TYM, and 191 down-regulated genes such as ADM, JUN, IRE2, PIGA, and MAFF. Gene set enrichment analysis (GSEA) analysis identified down-regulated pathways related to translation, transcription, immunity, PI3K/AKT, and circadian as well as disturbed pathways related to extracellular matrix and collagen. Splicing analysis identified 442 differential alternative splicing events within 284 genes in osteoarthritis, including genes involved in extracellular matrix (ECM) and alternative splicing, and TIA1 was identified as a key regulator of these splicing events.

CONCLUSIONS

These findings provide insights into disease etiology, and offer favorable information to support the development of more effective interventions in response to the global clinical challenge of osteoarthritis.

Analysis of the regulatory role of miR-34a-5p/PLCD3 in the progression of osteoarthritis

Osteoarthritis is a heterogeneous disease with a complex etiology. However, there is no effective treatment strategy at present. The purpose of this study was to explore the miRNA‒mRNA regulatory network and molecular mechanism that regulate the progression of osteoarthritis. In this article, we downloaded datasets (GSE55457, GSE82107, GSE143514 and GSE55235) from Gene Expression Omnibus (GEO) to screen differentially expressed mRNAs in osteoarthritis. Then, through weighted gene coexpression network (WGCNA), functional enrichment, protein‒protein interaction (PPI) network, miRNA‒mRNA coexpression network, ROC curve, and immune infiltration analyses and qPCR, the mRNA PLCD3, which was highly expressed in osteoarthritis and had clinical predictive value, was screened. We found that PLCD3 directly targets miR-34a-5p through DIANA and dual-luciferase experiments. The expression levels of PLCD3 and miR-34a-5p were negatively correlated. In addition, CCK-8 and wound healing assays showed that the miR-34a-5p mimic inhibited hFLS-OA cell proliferation and promoted hFLS-OA cell migration. PLCD3 overexpression showed the opposite trend. Western blotting further found that overexpression of miR-34a-5p reduced the protein expression levels of p-PI3K and p-AKT, while overexpression of PLCD3 showed the opposite trend. In addition, combined with the effect of the PI3K/AKT pathway inhibitor BIO (IC50 = 5.95 μM), the results showed that overexpression of miR-34a-5p increased the inhibitory effects of BIO on p-PI3K and p-AKT protein expression, while overexpression of PLCD3 significantly reversed these inhibitory effects. Overall, the miR-34a-5p/PLCD3 axis may mediate the PI3K/AKT pathway in regulating cartilage homeostasis in synovial osteoarthritis. These data indicate that miR-34a-5p/PLCD3 may be a new prognostic factor in the pathology of synovial osteoarthritis.

MicroRNA Signatures in Cartilage Ageing and Osteoarthritis

Osteoarthritis is the most common degenerative joint disorder. MicroRNAs are gene expression regulators that act post-transcriptionally to control tissue homeostasis. Microarray analysis was undertaken in osteoarthritic intact, lesioned and young intact cartilage. Principal component analysis showed that young intact cartilage samples were clustered together; osteoarthritic samples had a wider distribution; and osteoarthritic intact samples were separated into two subgroups, osteoarthritic-Intact-1 and osteoarthritic-Intact-2. We identified 318 differentially expressed microRNAs between young intact and osteoarthritic lesioned cartilage, 477 between young intact and osteoarthritic-Intact-1 cartilage and 332 between young intact and osteoarthritic-Intact-2 cartilage samples. For a selected list of differentially expressed microRNAs, results were verified in additional cartilage samples using qPCR. Of the validated DE microRNAs, four-miR-107, miR-143-3p, miR-361-5p and miR-379-5p-were selected for further experiments in human primary chondrocytes treated with IL-1β. Expression of these microRNAs decreased in human primary chondrocytes treated with IL-1β. For miR-107 and miR-143-3p, gain- and loss-of-function approaches were undertaken and associated target genes and molecular pathways were investigated using qPCR and mass spectrometry proteomics. Analyses showed that WNT4 and IHH, predicted targets of miR-107, had increased expression in osteoarthritic cartilage compared to young intact cartilage and in primary chondrocytes treated with miR-107 inhibitor, and decreased expression in primary chondrocytes treated with miR-107 mimic, suggesting a role of miR-107 in chondrocyte survival and proliferation. In addition, we identified an association between miR-143-3p and EIF2 signalling and cell survival. Our work supports the role of miR-107 and miR-143-3p in important chondrocyte mechanisms regulating proliferation, hypertrophy and protein translation.

MicroRNA-155 suppressed cholesterol-induced matrix degradation, pyroptosis and apoptosis by targeting RORα in nucleus pulposus cells

Intervertebral disc degeneration (IDD) is associated with low back pain, yet its inherent mechanism remains obscure. Hypercholesteremia was regarded as a risk factor for IDD, and our previous study showed that cholesterol accumulation could elicit matrix degradation in the nucleus pulposus (NP). MicroRNA-155 (miR-155) was substantiated as protective in IDD, but its role in cholesterol-induced IDD was unclear. The present study investigated whether miR-155 could mediate cholesterol-related IDD and its internal mechanisms. In vivo experiments revealed high-fat diet-induced hypercholesteremia in wild-type (WT) mice along with the occurrence of IDD, whereas Rm155LG transgenic mice showed milder NP degeneration, as evidenced by Safranin O-fast green (SF) staining and immunohistochemistry (IHC). Meanwhile, IHC showed that NLRP3 and Bax expression was also suppressed in Rm155LG mice. In vitro studies using Western blotting (WB) and immunofluorescence (IF) confirmed that the miR-155 mimic could alleviate cholesterol-induced matrix degradation, apoptosis and pyroptosis in NP. Moreover, RORα was upregulated in severely degenerated NP compared to mild IDD. It was also noted that RORα was suppressed in Rm155LG mice. In this study, we demonstrated that miR-155 could target RORα and that inhibition of RORα could prevent cholesterol-induced matrix degradation, apoptosis, and pyroptosis in NP, indicating the protective effect of miR-155 in cholesterol-induced IDD by targeting RORα.

Inhibition of miR-182-5p Targets FGF9 to Alleviate Osteoarthritis

Background. The pathogenesis of osteoarthritis (OA) is complex and there is no specific drug for treatment. The aim of this study was to identify the molecular targets of OA therapy, focusing on the expression and biological functions of miR-182-5p and its target genes in OA. Methods. miR-182-5p and fibroblast growth factor 9 (FGF9) were overexpressed or knocked down in IL-1β-induced chondrocytes. An OA knee model was performed by surgically destroying the medial meniscus. The gene expression of miR-182-5p and FGF9 was calculated. The protein FGF9 was tested by western blotting. Cell counting kit-8 (CCK8), plate cloning assay, and flow cytometry were conducted to evaluate cell proliferation and apoptosis. The expression of inflammatory factors, tumor necrosis factor-alpha (TNF-α), interleukin (IL)-6, and interleukin (IL)-8, was evaluated using enzyme-linked immunosorbent assay (ELISA). Dual-luciferase reporter assays validated the targeting relationship between miR-182-5p and FGF9. Hematoxylin–eosin (HE) and safranin O-fast Green (S–O) staining were utilized to access cartilage damage. Ki67 expression in cartilage was detected using immunohistochemistry (IHC). TdT-mediated dUTP nick-end labeling (TUNEL) assays were used to calculate the apoptosis rate of cartilage. Results. The expression of miR-182-5p was upregulated, and FGF9 was downregulated in the IL-1β-induced chondrocytes. OA chondrocytes proliferation ability in the miR-182-5p mimics group was decreased, and the apoptosis rate and inflammatory factor were increased. Transfection with miR-182-5p inhibitor increased the proliferative ability and decreased the apoptosis rate in the IL-1β-induced chondrocytes. Transfection with miR-182-5p inhibitor reversed IL-1β-induced inflammatory factor release in chondrocytes. Targeted binding sites existed between miR-182-5p and FGF9. After overexpression of FGF9, the miR-182-5p effect on OA chondrocytes was reversed. The hyaline cartilage thickness and proteoglycan content decreased in OA rats, and this was reversed by miR-182-5p inhibitor treatment. Conclusions. miR-182-5p expression levels were increased in OA chondrocytes and regulated chondrocyte proliferation, apoptosis, and inflammation by targeting FGF9. miR-182-5p is a potential gene for OA treatment.

MiR-760 targets HBEGF to control cartilage extracellular matrix degradation in osteoarthritis

The present study was developed to explore whether microRNA (miR)-760 targets heparin-binding EGF-like growth factor (HBEGF) to control cartilage extracellular matrix degradation in osteoarthritis. Both miR-760 and HBEGF expression levels were analysed in human degenerative cartilage tissues and in interleukin (IL)-1β/tumour necrosis factor (TNF)-α-treated chondrocytes in vitro. A series of knockdown and overexpression assays were then used to gauge the functional importance of miR-760 and HBEGF in OA, with qPCR and western immunoblotting analyses. Bioinformatics assays were used to identify putative miR-760 target genes, with these predictions then being validated through RNA pulldown and luciferase reporter assays. A murine anterior cruciate ligament transection model of OA was then established to prove the in vivo relevance of these findings. These experiments revealed that human degenerative cartilage tissues exhibited significant increases in miR-760 expression with a concomitant drop in HBEGF levels. IL-1β/TNF-α-treated chondrocytes also exhibited significant increases in miR-760 expression with a concomitant drop in HBEGF expression. When chondrocytes were transfected with either miR-760 inhibitor or HBEGF overexpression constructs, this was sufficient to interfere with degradation of the extracellular matrix (ECM). Moreover, miR-760 was confirmed to control chondrocyte matrix homeostasis by targeting HBEGF, and the overexpression of HBEGF partially reversed the effects of miR-760 mimic treatment on the degradation of the cartilage ECM. When OA model mice were administered an intra-articular knee injection of an adenoviral vector encoding a miR-760 mimic construct, cartilage ECM degradation was aggravated. Conversely, the overexpression of HBEGF in OA model mice partially reversed the effects of miR-760 overexpression, restoring appropriate ECM homeostasis. In summary, these data indicated that the miR-760/HBEGF axis plays a central role in orchestrating the pathogenesis of OA, making it a candidate target for therapeutic efforts in OA.

CCN4/WISP1 Promotes Migration of Human Primary Osteoarthritic Chondrocytes

OBJECTIVES

Previously, we have shown the involvement of cellular communication network factor 4/Wnt-activated protein Wnt-1-induced signaling protein 1 (CCN4/WISP1) in osteoarthritic (OA) cartilage and its detrimental effects on cartilage. Here, we investigated characteristics of CCN4 in chondrocyte biology by exploring correlations of CCN4 with genes expressed in human OA cartilage with functional follow-up.

DESIGN

Spearman correlation analysis was performed for genes correlating with CCN4 using our previously established RNA sequencing dataset of human preserved OA cartilage of the RAAK study, followed by a pathway enrichment analysis for genes with ρ ≥|0.6.| Chondrocyte migration in the absence or presence of CCN4 was determined in a scratch assay, measuring scratch size using a live cell imager for up to 36 h. Changes in expression levels of 12 genes, correlating with CCN4 and involved in migratory processes, were determined with reverse transcription-quantitative polymerase chain reaction (RT-qPCR).

RESULTS

Correlation of CCN4 with ρ ≥|0.6| was found for 58 genes in preserved human OA cartilage. Pathway analysis revealed "neural crest cell migration" as most significant enriched pathway, containing among others CORO1C, SEMA3C, and SMO. Addition of CCN4 to primary chondrocytes significantly enhance chondrocyte migration as demonstrated by reduced scratch size over the course of 36 h, but at the timepoints measured no effect was observed on mRNA expression of the 12 genes.

CONCLUSION

CCN4 increases cell migration of human primary OA chondrocytes. Since WISP1 expression is known to be increased in OA cartilage, this may serve to direct chondrocytes toward cartilage defects and orchestrate repair.

Osteoarthritis: pathogenic signaling pathways and therapeutic targets

Osteoarthritis (OA) is a chronic degenerative joint disorder that leads to disability and affects more than 500 million population worldwide. OA was believed to be caused by the wearing and tearing of articular cartilage, but it is now more commonly referred to as a chronic whole-joint disorder that is initiated with biochemical and cellular alterations in the synovial joint tissues, which leads to the histological and structural changes of the joint and ends up with the whole tissue dysfunction. Currently, there is no cure for OA, partly due to a lack of comprehensive understanding of the pathological mechanism of the initiation and progression of the disease. Therefore, a better understanding of pathological signaling pathways and key molecules involved in OA pathogenesis is crucial for therapeutic target design and drug development. In this review, we first summarize the epidemiology of OA, including its prevalence, incidence and burdens, and OA risk factors. We then focus on the roles and regulation of the pathological signaling pathways, such as Wnt/β-catenin, NF-κB, focal adhesion, HIFs, TGFβ/ΒΜP and FGF signaling pathways, and key regulators AMPK, mTOR, and RUNX2 in the onset and development of OA. In addition, the roles of factors associated with OA, including MMPs, ADAMTS/ADAMs, and PRG4, are discussed in detail. Finally, we provide updates on the current clinical therapies and clinical trials of biological treatments and drugs for OA. Research advances in basic knowledge of articular cartilage biology and OA pathogenesis will have a significant impact and translational value in developing OA therapeutic strategies.

Protective effects of Tripterygium glycoside on IL-1β-induced inflammation and apoptosis of rat chondrocytes via microRNA-216a-5p/TLR4/NF-κB axis

BACKGROUND

This study is designed to fill the research gap concerning the efficacy of Tripterygium glycoside (TG) on Interleukin-1β (IL-1β)-induced inflammation and injury in chondrocytes.

METHODS

Chondrocytes were isolated from Sprague-Dawley rats. After the treatment with IL-1β and TG and transfection, the viability and apoptosis of chondrocytes were determined via Cell Counting Kit-8 (CCK-8) assay and flow cytometry. The levels of inflammatory cytokines tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), and IL-8 were determined by enzyme-linked immunosorbent assay (ELISA). Relative expression levels of potential microRNAs (miRNAs, miRs) that may target toll-like receptor 4 (TLR4), as well as apoptosis- and TLR4/nuclear factor-κB (TLR4/NF-κB) pathway-associated factors were quantified using quantitative real-time (qRT) PCR and western blot. The targeting relationship between miR-216a-5p and TLR4 was predicted by TargetScan and further confirmed by dual-luciferase reporter assay.

RESULTS

The viability was reduced yet the apoptosis and inflammation were promoted in IL-1β-treated chondrocytes, where upregulation of Bax, Cleaved caspase 3, TLR4, Myeloid differentiation factor 88 (MyD88), phosphorylation of P65 and IκBα yet downregulation of Bcl-2 and IκBα were evidenced. Strikingly, the above changes were reversed by TG. TG also offset the effects of IL-1β on repressing the expression of miR-216a-5p, the miRNA targeting TLR4. Additionally, TLR4 overexpression neutralized the impacts of TG upon viability, apoptosis, and TLR4 expression in IL-1β-treated chondrocytes, while all these effects induced by TLR4 overexpression could be restored by miR-216a-5p.

CONCLUSIONS

TG protects chondrocytes against IL-1β-induced inflammation and apoptosis via miR-216a-5p/TLR4/NF-κB axis.

Identification and functional characterization of imbalanced osteoarthritis-associated fibronectin splice variants

OBJECTIVE

To identify FN1 transcripts associated with OA pathophysiology and investigate the downstream effects of modulating FN1 expression and relative transcript ratio.

METHODS

FN1 transcriptomic data was obtained from our previously assessed RNA-seq dataset of lesioned and preserved OA cartilage samples from the Research osteoArthritis Articular Cartilage (RAAK) study. Differential transcript expression analysis was performed on all 27 FN1 transcripts annotated in the Ensembl database. Human primary chondrocytes were transduced with lentiviral particles containing short hairpin RNA (shRNA) targeting full-length FN1 transcripts or non-targeting shRNA. Subsequently, matrix deposition was induced in our 3D in vitro neo-cartilage model. Effects of changes in the FN1 transcript ratio on sulphated glycosaminoglycan (sGAG) deposition were investigated by Alcian blue staining and dimethylmethylene blue assay. Moreover, gene expression levels of 17 cartilage-relevant markers were determined by reverse transcription quantitative polymerase chain reaction.

RESULTS

We identified 16 FN1 transcripts differentially expressed between lesioned and preserved cartilage. FN1-208, encoding migration-stimulating factor, was the most significantly differentially expressed protein coding transcript. Downregulation of full-length FN1 and a concomitant increased FN1-208 ratio resulted in decreased sGAG deposition as well as decreased ACAN and COL2A1 and increased ADAMTS-5, ITGB1 and ITGB5 gene expression levels.

CONCLUSION

We show that full-length FN1 downregulation and concomitant relative FN1-208 upregulation was unbeneficial for deposition of cartilage matrix, likely due to decreased availability of the classical RGD (Arg-Gly-Asp) integrin-binding site of fibronectin.

WWP2 confers risk to osteoarthritis by affecting cartilage matrix deposition via hypoxia associated genes

OBJECTIVE

To explore the co-expression network of the osteoarthritis (OA) risk gene WWP2 in articular cartilage and study cartilage characteristics when mimicking the effect of OA risk allele rs1052429-A on WWP2 expression in a human 3D in vitro model of cartilage.

METHOD

Co-expression behavior of WWP2 with genes expressed in lesioned OA articular cartilage (N = 35 samples) was explored. By applying lentiviral particle mediated WWP2 upregulation in 3D in vitro pellet cultures of human primary chondrocytes (N = 8 donors) the effects of upregulation on cartilage matrix deposition was evaluated. Finally, we transfected primary chondrocytes with miR-140 mimics to evaluate whether miR-140 and WWP2 are involved in similar pathways.

RESULTS

Upon performing Spearman correlations in lesioned OA cartilage, 98 highly correlating genes (|ρ| > 0.7) were identified. Among these genes, we identified GJA1, GDF10, STC2, WDR1, and WNK4. Subsequent upregulation of WWP2 on 3D chondrocyte pellet cultures resulted in a decreased expression of COL2A1 and ACAN and an increase in EPAS1 expression. Additionally, we observed a decreased expression of GDF10, STC2, and GJA1. Proteomics analysis identified 42 proteins being differentially expressed with WWP2 upregulation, which were enriched for ubiquitin conjugating enzyme activity. Finally, upregulation of miR-140 in 2D chondrocytes resulted in significant upregulation of WWP2 and WDR1.

CONCLUSIONS

Mimicking the effect of OA risk allele rs1052429-A on WWP2 expression initiates detrimental processes in the cartilage shown by a response in hypoxia associated genes EPAS1, GDF10, and GJA1 and a decrease in anabolic markers, COL2A1 and ACAN.

Insights from multi-omics integration in complex disease primary tissues

Genome-wide association studies (GWAS) have provided insights into the genetic basis of complex diseases. In the next step, integrative multi-omics approaches can characterize molecular profiles in relevant primary tissues to reveal the mechanisms that underlie disease development. Here, we highlight recent progress in four examples of complex diseases generated by integrative studies: type 2 diabetes (T2D), osteoarthritis, Alzheimer's disease (AD), and systemic lupus erythematosus (SLE). High-resolution methodologies such as single-cell and spatial omics techniques will become even more important in the future. Furthermore, we emphasize the urgent need to include as yet understudied cell types and increase the diversity of studied populations.

Extracellular vesicles in osteoarthritis of peripheral joint and temporomandibular joint

Osteoarthritis (OA) is a disabling disease with significant morbidity worldwide. OA attacks the large synovial joint, including the peripheral joints and temporomandibular joint (TMJ). As a representative of peripheral joint OA, knee OA shares similar symptoms with TMJ OA. However, these two joints also display differences based on their distinct development, anatomy, and physiology. Extracellular vesicles (EVs) are phospholipid bilayer nanoparticles, including exosomes, microvesicles, and apoptotic bodies. EVs contain proteins, lipids, DNA, micro-RNA, and mRNA that regulate tissue homeostasis and cell-to-cell communication, which play an essential role in the progression and treatment of OA. They are likely to partake in mechanical response, extracellular matrix degradation, and inflammatory regulation during OA. More evidence has shown that synovial fluid and synovium-derived EVs may serve as OA biomarkers. More importantly, mesenchymal stem cell-derived EV shows a therapeutic effect on OA. However, the different function of EVs in these two joints is largely unknown based on their distinct biological characteristic. Here, we reviewed the effects of EVs in OA progression and compared the difference between the knee joint and TMJ, and summarized their potential therapeutic role in the treatment of OA.

Identification and validation of immune and oxidative stress-related diagnostic markers for diabetic nephropathy by WGCNA and machine learning

Background

Diabetic nephropathy (DN) is the primary cause of end-stage renal disease, but existing therapeutics are limited. Therefore, novel molecular pathways that contribute to DN therapy and diagnostics are urgently needed.

Methods

Based on the Gene Expression Omnibus (GEO) database and Limma R package, we identified differentially expressed genes of DN and downloaded oxidative stress-related genes based on the Genecard database. Then, immune and oxidative stress-related hub genes were screened by combined WGCNA, machine learning, and protein-protein interaction (PPI) networks and validated by external validation sets. We conducted ROC analysis to assess the diagnostic efficacy of hub genes. The correlation of hub genes with clinical characteristics was analyzed by the Nephroseq v5 database. To understand the cellular clustering of hub genes in DN, we performed single nucleus RNA sequencing through the KIT database.

Results

Ultimately, we screened three hub genes, namely CD36, ITGB2, and SLC1A3, which were all up-regulated. According to ROC analysis, all three demonstrated excellent diagnostic efficacy. Correlation analysis revealed that the expression of hub genes was significantly correlated with the deterioration of renal function, and the results of single nucleus RNA sequencing showed that hub genes were mainly clustered in endothelial cells and leukocyte clusters.

Conclusion

By combining three machine learning algorithms with WGCNA analysis, this research identified three hub genes that could serve as novel targets for the diagnosis and therapy of DN.

MicroRNAs and long non-coding RNAs in cartilage homeostasis and osteoarthritis

During the last decade, osteoarthritis (OA) has become one of the most prevalent musculoskeletal diseases worldwide. OA is characterized by progressive loss of articular cartilage, abnormal remodeling of subchondral bone, hyperplasia of synovial cells, and growth of osteophytes, which lead to chronic pain and disability. The pathological mechanisms underlying OA initiation and progression are still poorly understood. Non-coding RNAs (ncRNAs) constitute a large portion of the transcriptome that do not encode proteins but function in numerous biological processes. Cumulating evidence has revealed a strong association between the changes in expression levels of ncRNA and the disease progression of OA. Moreover, loss- and gain-of-function studies utilizing transgenic animal models have demonstrated that ncRNAs exert vital functions in regulating cartilage homeostasis, degeneration, and regeneration, and changes in ncRNA expression can promote or decelerate the progression of OA through distinct molecular mechanisms. Recent studies highlighted the potential of ncRNAs to serve as diagnostic biomarkers, prognostic indicators, and therapeutic targets for OA. MiRNAs and lncRNAs are two major classes of ncRNAs that have been the most widely studied in cartilage tissues. In this review, we focused on miRNAs and lncRNAs and provided a comprehensive understanding of their functional roles as well as molecular mechanisms in cartilage homeostasis and OA pathogenesis.

Whole-genome RNA sequencing identifies distinct transcriptomic profiles in impingement cartilage between patients with femoroacetabular impingement and hip osteoarthritis

Femoroacetabular impingement (FAI) has a strong clinical association with the development of hip osteoarthritis (OA); however, the pathobiological mechanisms underlying the transition from focal impingement to global joint degeneration remain poorly understood. The purpose of this study is to use whole-genome RNA sequencing to identify and subsequently validate differentially expressed genes (DEGs) in femoral head articular cartilage samples from patients with FAI and hip OA secondary to FAI. Thirty-seven patients were included in the study with whole-genome RNA sequencing performed on 10 gender-matched patients in the FAI and OA cohorts and the remaining specimens were used for validation analyses. We identified a total of 3531 DEGs between the FAI and OA cohorts with multiple targets for genes implicated in canonical OA pathways. Quantitative reverse transcription-polymerase chain reaction validation confirmed increased expression of FGF18 and WNT16 in the FAI samples, while there was increased expression of MMP13 and ADAMTS4 in the OA samples. Expression levels of FGF18 and WNT16 were also higher in FAI samples with mild cartilage damage compared to FAI samples with severe cartilage damage or OA cartilage. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. We independently validated the results of the sequencing analysis and found increased expression of anabolic markers in patients with FAI and minimal histologic cartilage damage, suggesting that anabolic signaling may be increased in early FAI with a transition to catabolic and inflammatory gene expression as FAI progresses towards more severe hip OA. Clinical significance:Cam-type FAI has a strong clinical association with hip OA; however, the cellular pathophysiology of disease progression remains poorly understood. Several previous studies have demonstrated increased expression of inflammatory markers in FAI cartilage samples, suggesting the involvement of these inflammatory pathways in the disease progression. Our study further expands the knowledge regarding distinct genetic reprogramming in the cartilage between FAI and hip OA patients. In addition to differences in inflammatory gene expression, we also identified differential expression in multiple pathways involved in hip OA progression.

Phosphoproteomics reveals the BRAF-ERK1/2 axis as an important pathogenic signaling node in cartilage degeneration

OBJECTIVE

Osteoarthritis (OA) causes gradual cellular alterations, structural anomalies and joint dysfunction. Progressive decline of chondrocyte function plays a vital role on OA pathogenesis. Although protein phosphorylation controls cartilage metabolism, its regulation mechanism in OA remains unclear. Thus, proteomic methods were used to investigate phosphorylation changes in preserved and OA articular cartilage samples, and to explore the intervention targets of phosphorylated kinase.

METHODS

Preserved (control) and lesioned (OA) cartilage samples from OA cases were assessed by phosphoproteomics. Immobilized metal affinity chromatography was performed for phosphopeptide enrichment. Quantitated phosphosites were comparatively assessed in the cartilage sample pair. Kinase-substrate enrichment analyses were carried out for identifying OA-related kinases. BRAF expression in cartilage tissues was assessed by immunohistochemical staining. The effects of BRAF inhibitor on cartilage degeneration were examined in mouse chondrocytes and OA mouse model.

RESULTS

High-sensitivity mass spectrometry-based proteomics revealed 7,471 peptides and 4,375 phosphorylated peptides differing between preserved and lesioned cartilage samples, which represented the significant alteration of kinase hubs and transduction pathways. Phosphoproteomics identified BRAF may be involved in developing OA. BRAF regulated the downstream ERK signaling pathway. In addition, BRAF was upregulated in human OA cartilage as shown by immunohistochemistry. Remarkably, BRAF inhibition alleviated cartilage degradation in a mouse model of OA through its downstream of ERK pathway activation.

CONCLUSIONS

Jointly, these findings provide an overview of phosphoproteomic alterations occurring during cartilage degeneration, identifying the BRAF-ERK1/2 Axis signaling as a potential signaling pathway involved in OA.

miR-210-3p protects against osteoarthritis through inhibiting subchondral angiogenesis by targeting the expression of TGFBR1 and ID4

Excessive subchondral angiogenesis is a key pathological feature of osteoarthritis (OA), as it alters the balance of subchondral bone remodeling and causes progressive cartilage degradation. We previously found that miR-210-3p correlates negatively with angiogenesis, though the specific mechanism of miR-210-3p-related angiogenesis in subchondral bone during OA progression remains unclear. This study was conducted to identify the miR-210-3p-modulating subchondral angiogenesis mechanism in OA and investigate its therapeutic effect. We found that miR-210-3p expression correlated negatively with subchondral endomucin positive (Emcn+) vasculature in the knee joints of OA mice. miR-210-3p overexpression regulated the angiogenic ability of endothelial cells (ECs) under hypoxic conditions in vitro. Mechanistically, miR-210-3p inhibited ECs angiogenesis by suppressing transforming growth factor beta receptor 1 (TGFBR1) mRNA translation and degrading DNA-binding inhibitor 4 (ID4) mRNA. In addition, TGFBR1 downregulated the expression of ID4. Reduced ID4 levels led to a negative feedback regulation of TGFBR1, enhancing the inhibitory effect of miR-210-3p on angiogenesis. In OA mice, miR-210-3p overexpression in ECs via adeno-associated virus (AAV) alleviated cartilage degradation, suppressed the type 17 immune response and relieved symptoms by attenuating subchondral Emcn+ vasculature and subchondral bone remodeling. In conclusion, we identified a miR-210-3p/TGFBR1/ID4 axis in subchondral ECs that modulates OA progression via subchondral angiogenesis, representing a potential OA therapy target.

Current understanding of osteoarthritis pathogenesis and relevant new approaches

Osteoarthritis (OA) is the most common degenerative joint disease that causes painful swelling and permanent damage to the joints in the body. The molecular mechanisms of OA are currently unknown. OA is a heterogeneous disease that affects the entire joint, and multiple tissues are altered during OA development. To better understand the pathological mechanisms of OA, new approaches, methods, and techniques need to be used to understand OA pathogenesis. In this review, we first focus on the epigenetic regulation of OA, with a particular focus on DNA methylation, histone modification, and microRNA regulation, followed by a summary of several key mediators in OA-associated pain. We then introduce several innovative techniques that have been and will continue to be used in the fields of OA and OA-associated pain, such as CRISPR, scRNA sequencing, and lineage tracing. Next, we discuss the timely updates concerning cell death regulation in OA pathology, including pyroptosis, ferroptosis, and autophagy, as well as their individual roles in OA and potential molecular targets in treating OA. Finally, our review highlights new directions on the role of the synovial lymphatic system in OA. An improved understanding of OA pathogenesis will aid in the development of more specific and effective therapeutic interventions for OA.

MiR-99a alleviates apoptosis and extracellular matrix degradation in experimentally induced spine osteoarthritis by targeting FZD8

Background

Our previous study identified miR-99a as a negative regulator of early chondrogenic differentiation. However, the functional role of miR-99a in the pathogenesis of osteoarthritis (OA) remains unclear.

Methods

We examined the levels of miR-99a and Frizzled 8 (FZD8) expression in tissue specimens. Human SW1353 chondrosarcoma cells were stimulated with IL-6 and TNF-α to construct an in vitro OA environment. A luciferase reporter assay was performed to analyze the relationship between miR-99a and FZD8. CCK-8 assays, flow cytometry, and ELISA assays were used to assess cell viability, apoptosis, and inflammatory molecule expression, respectively. Percutaneous intra-spinal injections of papain mixed solution were performed to create an OA Sprague–Dawley rat model. Alcian Blue staining, Safranin O Fast Green staining, and Toluidine Blue O staining were performed to detect the degrees of cartilage injury.

Results

MiR-99a expression was downregulated in the severe spine OA patients when compared with the mild spine OA patients, and was also decreased in the experimentally induced in vitro OA environment when compared with the control environment. Functionally, overexpression of miR-99a significantly suppressed cell apoptosis and extracellular matrix degradation stimulated by IL-6 and TNF-α. FZD8 was identified as a target gene of miR-99a. Furthermore, the suppressive effects of miR-99a on cell injury induced by IL-6 and TNF-α were reversed by FZD8 overexpression. Moreover, the levels of miR-99a expression were also reduced in the induced OA model rats, and miR-99a agomir injection relieved the cartilage damage. At the molecular level, miR-99a overexpression downregulated the levels of MMP13, β-catenin, Bax, and caspase-3 protein expression and upregulated the levels of COL2A1 and Bcl-2 protein expression in the in vitro OA-like chondrocyte model and also in the experimental OA model rats.

Conclusions

Our data showed that miR-99a alleviated apoptosis and extracellular matrix degradation by targeting FZD8, and thereby suppressed the development and progression of experimentally induced spine osteoarthritis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05822-8.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis

OBJECTIVES

To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients.

METHODS

RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the GWAS catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk-SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP.

RESULTS

We identified 30 transcript SNPs (28 in cartilage and 2 in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7x10-4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75x10-5). Moreover, 12 drugs were interacting with 7 genes displaying AEI, of which 7 drugs have been already approved.

CONCLUSIONS

By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harboring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further, translation towards underlying biological mechanisms.

Circular RNA hsa_circ_0002360 Promotes Proliferation and Invasion and Inhibits Oxidative Stress in Gastric Cancer by Sponging miR-629-3p and Regulating the PDLIM4 Expression

Many studies have found that circRNA hsa_0002360 (circ0002360) plays an important role in cancer onset and progression. However, its role in gastric cancer (GC) remains uncertain. Circ0002360 was found to be upregulated in GC cells using QRT-PCR. Furthermore, miR-629-3p, a target miRNA of circ0002360, was the most suppressed miRNA following circ0002360 overexpression. RNA immunoprecipitation (RIP), dual-luciferase reporter analyses, clone formation, transwell, DCFH-DA, and ELISA assays demonstrated that circ0002360-targeted miR-629-3p promotes cell proliferation and migration while inhibiting oxidative stress. GC-related mRNA microarrays from the GEO and TCGA databases, including GSE103236, GSE79973, GSE33429, GSE22804, GSE84437, and TCGA-STAD datasets, were used to find hub biomarkers between normal and gastric cancer samples. WGCNA and uni-Cox analysis were used to identify 27 survival-related risk genes, which were then used to build a risk model for prognosis prediction. Following that, all patients from the GSE84437 and TCGA-STAD datasets with 27 survival-related genes and enough data on survival status and time were randomly assigned to train (n = 433) and test (n = 375) cohorts. Furthermore, ROC and Kaplan-Meier (KM) analyses were used to validate the risk model for both cohorts. randomForest analysis indicated that PDLIM4 was the target gene of miR-629-3p, whose level was increased by circ0002360 but reversed by miR-629-3p mimics. Finally, this study confirmed that circ0002360 sponged miR-629-3p and then upregulated PDLIM4 expression. As a result, circ0002360 may be a useful marker for predicting GC prognosis and an anti-GC treatment target.

miR-132-3p participates in the pathological mechanism of temporomandibular joint osteoarthritis by targeting PTEN

OBJECTIVE

This study aimed to investigate the role of miR-132-3p in the progression of temporomandibular joint osteoarthritis (TMJOA) and its potential pathological mechanism.

DESIGN

A TMJOA model was established using six rats via the unilateral anterior crossbite method. The differential expression of miR-132-3p in the TMJOA (n = 6) and control groups (n = 6) was detected via miRNA sequencing and verified via PCR. The chondrocytes in the condylar cartilage of the temporomandibular joint were cultured and stimulated with IL-1β to simulate TMJOA in vitro. The changes in the proliferation, apoptosis, inflammation and extracellular matrix of these chondrocytes were detected after the upregulation of miR-132-3p expression. The targeted relationship of miR-132-3p and PTEN in TMJOA was verified, and rescue experiments were conducted via co-upregulation of the expression of both miR-132-3p and PTEN.

RESULTS

Compared with that in the control group, miR-132-3p expression was lower in the cartilage tissues of TMJOA rats and IL-1β-induced TMJ chondrocytes. After upregulating the expression of miR-132-3p, the cell proliferation activity and expression levels of aggrecan and type II collagen of IL-1β-induced TMJ chondrocytes were increased, and the apoptosis rate and levels of inflammatory factors were decreased. miR-132-3p can regulate PTEN expression in a targeted manner, and upregulating PTEN expression could reverse the influences of the upregulation of miR-132-3p expression on TMJOA cells.

CONCLUSION

miR-132-3p is less expressed in TMJOA, and it regulates the proliferation, extracellular matrix, and inflammatory response of TMJOA chondrocytes and participates in TMJOA progression by targeting PTEN.

The Involvement of Neutrophils in the Pathophysiology and Treatment of Osteoarthritis

Osteoarthritis (OA) is a chronic disability that significantly impairs quality of life. OA is one of the most prevalent joint pathologies in the world, characterized by joint pain and stiffness due to the degeneration of articular cartilage and the remodeling of subchondral bone. OA pathogenesis is unique in that it involves simultaneous reparative and degradative mechanisms. Low-grade inflammation as opposed to high-grade allows for this coexistence. Previously, macrophages and T cells have been identified as playing major roles in the inflammation and destruction of OA joints, but recent studies have demonstrated that neutrophils also contribute to the pathogenesis. Neutrophils are the first immune cells to enter the synovium after joint injury, and neutrophilic activity is indispensably a requisite for the progression of OA. Neutrophils act through multiple mechanisms including tissue degeneration via neutrophil elastase (NE), osteophyte development, and the release of inflammatory cytokines and chemokines. As the actions of neutrophils in OA are discovered, the potential for novel therapeutic targets as well as diagnostic methods are revealed. The use of chondrogenic progenitor cells (CPCs), microRNAs, and exosomes are among the newest therapeutic advances in OA treatment, and this review reveals how they can be used to mitigate destructive neutrophil activity.

Otto Aufranc Award: Identification of Key Molecular Players in the Progression of Hip Osteoarthritis Through Transcriptomes and Epigenetics

BACKGROUND

This study aimed: (1) to compare the transcriptome profile of articular cartilage in cam-FAI (early stage) to advanced OA secondary to cam-FAI (late stage) and (2) to investigate epigenetic changes through the expression of DNA methylation enzymes DNMT3B, DNMT1, and DNMT3A and peroxisome proliferator-activated receptor gamma (PPARγ) in human cartilage samples during the progression of hip OA.

METHODS

Full-thickness cartilage samples were collected from the anterolateral head-neck junction (impingement zone) of 22 patients (9 early-FAI and 13 late-FAI). RNA sequencing and in vitro cartilage cultures with histological analysis and immunohistochemistry staining for PPARγ and DNMT3B were performed. Target gene validation was confirmed with RT-PCR.

RESULTS

Fifty genes and 42 pathways were identified differentially between early and late-FAI (fold change <-1.5 or >1.5, P < .01). PPARγ and DNMT3B were gradually suppressed with disease progression. Contrarily, disease progression induced expression of DNMT1/3A.

CONCLUSION

By comparing comprehensive gene expression in early and late stage hip degeneration at the whole-genome level, distinct transcriptome profiles for early and late stage disease were identified along with key molecular contributors to the progression of hip OA. Preservation of endogenous PPARγ may have therapeutic potential to delay or prevent hip OA.

Potential Methods of Targeting Cellular Aging Hallmarks to Reverse Osteoarthritic Phenotype of Chondrocytes

Osteoarthritis (OA) is a chronic degenerative joint disease that causes pain, physical disability, and life quality impairment. The pathophysiology of OA remains largely unclear, and currently no FDA-approved disease-modifying OA drugs (DMOADs) are available. As has been acknowledged, aging is the primary independent risk factor for OA, but the mechanisms underlying such a connection are not fully understood. In this review, we first revisit the changes in OA chondrocytes from the perspective of cellular hallmarks of aging. It is concluded that OA chondrocytes share many alterations similar to cellular aging. Next, based on the findings from studies on other cell types and diseases, we propose methods that can potentially reverse osteoarthritic phenotype of chondrocytes back to a healthier state. Lastly, current challenges and future perspectives are summarized.

Pentraxin 3 regulated by miR-224-5p modulates macrophage reprogramming and exacerbates osteoarthritis associated synovitis by targeting CD32

Emerging evidence has shown an imbalance in M1/M2 macrophage polarization to play an essential role in osteoarthritis (OA) progression. However, the underlying mechanistic basis for this polarization is unknown. RNA sequencing of OA M1-polarized macrophages found highly expressed levels of pentraxin 3 (PTX3), suggesting a role for PTX3 in OA occurrence and development. Herein, PTX3 was found to be increased in the synovium and articular cartilage of OA patients and OA mice. Intra-articular injection of PTX3 aggravated, while PTX3 neutralization reversed synovitis and cartilage degeneration. No metabolic disorder or proteoglycan loss were observed in cartilage explants when treated with PTX3 alone. However, cartilage explants exhibited an OA phenotype when treated with culture supernatants of macrophages stimulated with PTX3, suggesting that PTX3 did not have a direct effect on chondrocytes. Therefore, the OA anti-chondrogenic effects of PTX3 are primarily mediated through macrophages. Mechanistically, PTX3 was upregulated by miR-224-5p deficiency, which activated the p65/NF-κB pathway to promote M1 macrophage polarization by targeting CD32. CD32 was expressed by macrophages, that when stimulated with PTX3, secreted abundant pro-inflammation cytokines that induced severe articular cartilage damage. The paracrine interaction between macrophages and chondrocytes produced a feedback loop that enhanced synovitis and cartilage damage. The findings of this study identified a functional pathway important to OA development. Blockade of this pathway and PTX3 may prevent and treat OA.

miRNA-432 and SLC38A1 as Predictors of Hepatocellular Carcinoma Complicated with Alcoholic Steatohepatitis

Alcoholic steatohepatitis (ASH) is asymptomatic in the early stages and is typically advanced at the time of diagnosis. With the global rise in alcohol abuse, ASH is currently among the most detrimental diseases around the world. Hepatocellular carcinoma (HCC) is one of the final outcomes of numerous liver diseases. However, at present, HCC screening is mostly focused on liver cancer development. Moreover, there is no effective biomarker to predict the prognosis and recurrence of liver cancer. Meanwhile, there are limited studies on the prognosis and recurrence of HCC patients complicated with ASH. In this study, using bioinformatic analysis as well as cellular and animal models, we screened the differentially expressed (DE) miRNA-432 and SLC38A1 gene in ASH. Based on our analysis, miRNA-432 targeted SLC38A1, and the levels of miRNA-432 and SLC38A1 could accurately predict the overall survival (OS) and relapse free survival (RFS) in patients with liver cancer. Hence, these two genetic elements have the potential to synergistically predict the prognosis and recurrence of HCC complicated with ASH.

Identification of Key Genes and Pathways Associated with PIEZO1 in Bone-Related Disease Based on Bioinformatics

PIEZO1 is a mechano-sensitive ion channel that can sense various forms of mechanical stimuli and convert them into biological signals, affecting bone-related diseases. The present study aimed to identify key genes and signaling pathways in Piezo1-regulated bone-related diseases and to explain the potential mechanisms using bioinformatic analysis. The differentially expressed genes (DEGs) in tendon, femur, and humerus bone tissue; cortical bone; and bone-marrow-derived macrophages were identified with the criteria of |log2FC| > 1 and adjusted p-value < 0.05 analysis based on a dataset from GSE169261, GSE139121, GSE135282, and GSE133069, respectively, and visualized in a volcano plot. Venn diagram analyses were performed to identify the overlapping DEGs expressed in the above-mentioned tissues. Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, protein−protein interaction (PPI) analysis, and module analysis were also conducted. Furthermore, qRT-PCR was performed to validate the above results using primary chondrocytes. As a result, a total of 222 overlapping DEGs and 12 mostly overlapping DEGs were identified. Key Piezo1-related genes, such as Lcn2, Dkk3, Obscn, and Tnnt1, were identified, and pathways, such as Wnt/β-catenin and PI3k-Akt, were also identified. The present informatic study provides insight, for the first time, into the potential therapeutic targets of Piezo1-regulated bone-related diseases