Introduction

Osteoarthritis (OA) is a degenerative disease characterized by the deterioration of articular carti-lage and affecting all components of the joint, including the synovium, subchondral bone, and meniscus1. Knee OA (KOA) is the leading cause of disability and increased living costs among elderly patients2. Currently, the treatment options for osteoarthritis are limited to symptomatic relief and joint replacement. However, the insufficient relief of symptoms, potential medication side effects, and the economic burden and complications associated with joint replacement sur-geries have prompted the need to identify new, effective, and safe treatment methods and explore unknown targets for individuals with KOA3.

Exosomes are extracellular vesicles that are distributed in body fluids such as serum and play a crucial biological role by delivering molecules such as MicroRNAs (miRNAs)4. They serve as vital carriers for intercellular communication and transfer of genetic information. MiRNAs are a class of non-coding RNAs ranging from 18 to 24 nucleotides in length, which exert inhibitory effects on the expression of target genes through interactions with mRNA. Studies have demon-strated that miRNAs play a crucial role as a pathogenic factor in osteoarthritis (OA). 5. For in-stance, miR-199b-5p was found to contribute to the osteogenic differentiation of bone marrow stromal cells 6; another miR-140 showed cartilage-specific expression, and its expression was significantly reduced in OA cartilage 7;8.recently, intra-articular injection of antisense oligonu-cleotides of miR-181a-5p produced chondroprotective effects in OA mice 9.

In this study, we initially screened miR-199b-5p as a potential key miRNA based on clinical data by detecting differentially expressed exosomal miRNAs. To elucidate the role and function of miR-199b-5p, we investigated its impact on chondrocytes, which are crucial pathogenic cells in KOA, and identified its molecular targets through in vitro experiments. Furthermore, we ex-plored its role in vivo. Hence, this article not only identifies miR-199b-5p as a potential micro-target for KOA but also provides a potential strategy for future identification of new molecular drugs.

Results

Identification and enrichment analysis of differentially expressed miRNAs in serum exo-somes

To investigate the dysregulated miRNAs in serum exosomes of KOA patients, we extracted exo-somal miRNAs from serum and performed sequencing. Serum samples from 15 patients with KOA and 10 healthy subjects were collected (Supplementary Table 2). After extraction, the se-rum exosomes were observed under transmission electron microscopy and nanoparticle tracking analysis, revealing a diameter ranging from approximately 70 to 150 nm (Fig. S1A, B). Further-more, the surface marker proteins CD9, CD63, and CD81 were found to be expressed on the ex-osomes (Fig. S1C). Next, we performed sequencing of miRNAs in serum exosomes.

The results showed that 88 miRNAs were up-regulated and 89 miRNAs were downregulated in KOA patients compared with the control group based on fold change > 1.5 and p < 0.05 (Fig. 1A, B). Afterwards, we performed bioinformatics Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses on these differentially expressed genes. Among the top 5 enriched results of up-regulated differentially expressed miRNAs, the CC (Cell, intracellular, and cytoplasm) GO terms were notable. The MF (Translation activator activity, binding to ubiq-uitin-conjugating enzymes, and possessing transferase activity) GO terms and BP (Cell differen-tiation and the stabilization of mRNA through 3’-UTR mediation) GO terms were also significant (Fig. 1C). The top 10 enriched KEGG pathways, PI3K-AKT signaling, Ras signaling, apoptosis, and other types of O-glycan biosynthesis are considered to be of significant importance (Fig. 1D).The analysis of the top 5 down-regulated differential miRNAs reveals interesting findings. the CC (Cell, intracellular) GO terms were particularly noteworthy. The MF (RNA polymerase II transcription factor activity and sequence-specific DNA binding) GO terms, as well as the BP (collagen-activated signaling pathway) GO terms, were also deemed significant (Fig. 1E). More-over, the analysis of the top 10 enriched KEGG pathways revealed a strong association with the Hedgehog signaling pathway, pluripotency of stem cells, and glycosaminoglycan biosynthesis-keratan sulfate (Fig. 1F).

MiRNA expression in KOA patients.

(A, B) Volcano plot and heatmap showing dif-ferential miRNAs between KOA and healthy groups (healthy, n = 10; KOA, n = 15, | log2 Fold Change | 0.585, P < 0.05). (C) GO and KEGG enrichment analyses of upregulated genes. (D) GO and KEGG enrichment analyses of downregulated genes.

It is noteworthy that in the aforementioned results, both the up-regulated and down-regulated dif-ferential miRNA enrichment analyses are involved in processes related to cartilage synthesis, including O-glycan biosynthesis and glycosaminoglycan biosynthesis-keratan sulfate10;11.

Dysregulated serum exosomal miRNAs were found to be expressed in both serum and car-tilage

Subsequently, we aimed to investigate whether these dysregulated miRNAs in serum exosomes show differential expression in other affected sites or tissues of KOA. Leveraging the KOA pa-tient cartilage and serum sequencing data available in the GEO database, we compared our dysregulated miRNA list with the datasets. Remarkably, we identified 169 miRNAs that exhibit-ed differential expression in the serum of KOA patients (Fig. 2A), suggesting their involvement in the disease. Moreover, these miRNAs were also found to be expressed in KOA patient carti-lage (Fig. 2B). This robust validation reinforces the reliability of our data. Consequently, we pro-ceeded with further screening of differentially expressed genes.

Verification with GEO data and miRNA screening.

(A, B) Venn plot showing GEO dataset and our result. (C) RT-qPCR in clinical samples to verify the expression (HC, n = 5, KOA, n = 5). (D) RT-qPCR results of mouse joint samples (n = 6). Data are shown as means ± SD. *P<0.05, **P<0.01.

miRNA-199b-5p has been identified as a potential target molecule in KOA

Based on the P-value and exosomal expression, five miRNAs(hsa-mir-3168, hsa-mir-1296-5p, hsa-mir-15b-3p, hsa-mir-338-3p and hsa-mir-199b-5p) were selected to further research and val-idated in independent human samples by RT-qPCR (hsa-mir-199b-5p, p < 0.01; hsa-mir-3168, p=0.33; hsa-mir-15b-3p, p < 0.01; hsa-mir-338-3p, p=0.20 and hsa-mir-1296-5p, p < 0.01) (Fig. 2C). To further explore the functional roles of these miRNAs, we established a mouse model of KOA to evaluate their expression in the joints. Due to species limitations, we examined the ex-pression of miR-199b-5p (p <0.01), miR-15b-3p (p <0.05) and miR-338-3p (p <0.05) in joint tissue samples (control vs M-0.5) (Fig. 2D). The results showed that only the expression trend of miR-199b-5p was consistent between the clinical samples and the mouse arthritis model. There-fore, we selected miR-199b-5p as the target for our subsequent research.

miR-199b-5p mimic inhibits the cell viability and extracellular matrix (ECM) of chondro-cytes and inhibitor restores LPS-induced chondrocytes damage

We first extracted mouse primary chondrocytes (Fig. 3A-C). In vitro experiments showed that overexpression miR-199b-5p inhibited the viability of chondrocytes(p<0.01) (Fig. 3D, Fig. S2). we also find miR-199b-5p mimic increased the mRNA expressions of MMP3 (p=0.09) and ADAMTS5 (p<0.05) and decreased the mRNA expression of COL2A1(p=0.05), AGGRECAN (p=0.20) and SOX9 (p=0.22), which are often used as the biomarkers of chondrocytes ECM met-abolic balance. In contrast, miR-199b-5p inhibitor decreased the mRNA expression of MMP3 (p<0.01) and ADAMTS5 (p=0.11) and increased the mRNA expression of COL2A1 (p=0.07), AGGRECAN (p<0.01) and SOX9 (p<0.01) (Fig. 3E–I). While some gene expression changes may not be significant in statistic, but the modulation of miR-199b-5p expression has been ob-served to exert an influence on the metabolic alterations of chondrocytes.

Chondrocyte proliferation and marker expression changes after miR-199b-5p mimic or inhibitor treatment.

(A) Second generation primary mouse chondrocytes. (B) Toluene blue staining. (C) Type Ⅱ collagen immunoassay. (D) CCK-8 assay for cell viability (n=6). (E, F) RT-qPCR detection of MMP-3 and ADAMTS5 mRNA expression (n=3). (G-I) RT-qPCR detec-tion of COL-2A1, AGGRECAN, and SOX9 mRNA expression (n=3). (J) CCK-8 cell viability assay after different doses of LPS induction (n=5). (K) CCK-8 cell viability assay after virus in-fection (n=6). (L, M) RT-qPCR detection of MMP-3 and ADAMTS5 mRNA expression (n=3). (N-P) RT-qPCR detection of COL2A1, AGGRECAN, and SOX9 mRNA expression (n=3). Data are shown as means ± SD. *P<0.05, **P<0.01.

To explore the effect of miR-199b-5p under pathological conditions, an LPS-stimulated inflam-mation chondrocyte cell model was established. We examined the effect of LPS at 5, 10, and 15 μg/ml on cell viability and found that cell viability was significantly decreased at 15 μg/ml(p<0.001) (Fig. 2J). Next, we chose 15 μg/ml of LPS to establish a chondrocyte injury model. The CCK-8 assay showed that miR-199b-5p inhibitor reversed the decrease of cell viabil-ity caused by LPS (p<0.01) (Fig. 3K). Also, we revealed that LPS elevated MMP3 (p<0.01) and ADAMTS5 (p<0.01) and decreased COL2A1(p=0.06), AGGRECAN (p<0.01), and SOX9 (p=0.13) expression. In the presence of miR-199b-5p inhibitor, the changes in mRNA levels were re-versed (Fig. 3L-P). These results suggest that miR-199b-5p overexpression reduces cell viability and miR-199b-5p inhibition partly restores LPS-induced cell damage and ECM degradation.

miR-199b-5p mimic induces inflammation in normal mice, and miR-199b-5p inhibitor alle-viates symptoms in KOA mice

Now, we want to know the vivo role of miRNA-199b-5p. Firstly, we screened Adneovirus(AD) and utilized High-AD as a vector to either overexpress or inhibit the expression of miR-199b- 5p(Fig. S3). Then, miR-199b-5p mimic was injected to the joint of mice (Fig. 4A) and a decrease in pain threshold was found (p<0.01) (Fig. 4B). Four weeks later, serum IFN-γ (p<0.01) and TNF-α (p<0.01) were also significantly increased (Fig. 4C, D). The cartilage of mice injected with the miR-199b-5p mimic was slightly degraded (p=0.02) (Fig. 4E, F). Additionally, the ar-ticular surface showed insect erosion (Fig. 4G). interestingly, the level of serum inflammation in the miR-199b-5p inhibitor injection group was significantly decreased compared with the mimic group (IFN-γ, p<0.01; TNF-α, p<0.01). These results indicated that intra-articular injection of miR-199b-5p mimic induced inflammation response in mice.

Injection of adenovirus expressing miR-199b-5p mimic results in inflammation and pain threshold sensitivity in mice.

(A) Animal experiment schematic. (B) Behavioral detection of animal thermal pain threshold. (C, D) Detection of serum levels of IFN-γ and TNF-α in controls by ELISA. (E, F) Safranin-fast green staining and semiquantitative scoring of articular cartilage. (G) 3D reconstruction and 2D images of joints from μCT scans. Data are shown as mean ± SD. *P<0.05, **P<0.01, n=6.

Followingly, we established an MIA-induced KOA model to further investigate the role of miR- 199b-5p (Fig. S4). We observed a decrease in pain threshold in the model group, and recovery was observed on the 10th day in the inhibitor group (p<0.01) (Fig. 5A, B). Joint tissues were tak-en at the fourth week, and revealed that the expression of IFN-γ (p<0.01) and TNF-α (p<0.01) decreased after inhibiting miRNA-199b-5p (Fig. 5C, D). Safranin-fast green staining of joints showed recovery of articular cartilage degradation (p<0.05) (Fig. 5E, F), and the moth erosion of the cartilage μCT was partly improved (Fig. 5G). These results proved that intra-articular injec-tion of the miR-199b-5p inhibitor partly recovered pain, inflammation, and cartilage degenera-tion in KOA mice.

Injection of adenovirus expressing miR-199b-5p inhibitor partly recovered pathologi-cal changes in KOA mice.

(A) Animal experiment schematic. (B) Behavioral detection of animal thermal pain threshold. (C, D) Detection of serum levels of IFN-γ and TNF-α in controls by ELISA. (E, F) Safranin-fast green staining and semiquantitative scoring of articular cartilage. (G) 3D reconstruction and 2D images of joints from μCT scans. Data are shown as mean ± SD. *P<0.05, **P<0.01, n=6.

Gcnt2 and Fzd6 are two target genes of miR-199b-5p

In order to investigate the underlying mechanism of miRNA-199b-5p, we utilized five widely used miRNA target gene prediction tools, namely miRWalk, miRDB, TarBase, starbase, and TargetScan, to identify potential target genes. Consequently, we identified six putative target genes of miRNA-199b-5p (Fig. 6A). Bioinformatics analysis of the six possible target genes showed that BP is in posttranscriptional regulation of gene expression and angiogenesis; CC is in cytoplasmic vesicle and cell leading edge; and MF is in ubiquitin protein ligase binding and pro-tein domain specific binding (Fig. 6B).

Validation of the miR-199b-5p target gene.

(A) Prediction of target genes of miR- 199b-5p using search sites. (B) Target gene GO analysis. (C, D) Detection of the expression of target genes under different conditions (n=3). (E, F) Predicting the binding site of miR-199b-5p and target genes. (G, H) Validation by luciferase reporter gene assay (n=3). Data are show as mean ± SD.*P<0.05, **P<0.01, n=6. (I-L) The expression of FZD6 and GCNT2 in the synovial membrane and chondrocytes of GEO Profiles KOA.

After Mimic infected cells, we found decreases expression in Fzd6 (P<0.01), Gcnt2 (P<0.01), and Caprin1 (P=0.024) in the chondrocytes (Fig. 6C). Conversely, after inhibitor infection, Hif-1α (P=0.048), Fzd6 (P<0.01) and Gcnt2 (P<0.01) were increased and Atg14 (P=0.042) increased (Fig. 6D). Notably, we observed corresponding changes in the expression of Fzd6 and Gcnt2 up-on miR-199b-5p overexpression and under-expression. Moreover, we predicted potential binding sites for miRNA-199b-5p within the 3’-untranslated region (UTR) of these two target genes (Fig. 6E, F) and luciferase reporter assays confirmed that miR-199b-5p can bind to and suppress the expression of both Gcnt2 (P < 0.01) and Fzd6 (P < 0.01) via their complementary sequences (Fig. 6G, H). Furthermore, we found differential expression of Fzd6 in both synovial tissue data (GDS5401) and chondrocyte data (GDS3758) from KOA patients in the GEO profile (Fig. 6I, J). Similarly, differential expression tendency of Gcnt2 was observed in both the synovial tissue da-ta (GDS5403) and GDS3758 from the same KOA patients (Fig. 6K, L). These findings provide further validation to our results. Therefore, Fzd6 and Gcnt2 was confirmed as an important downstream target for further research.

Discussion

In this study, we initially performed sequencing of serum exosomal miRNAs from clinical pa-tients and identified 177 dysregulated miRNAs. Subsequently, through comparison with GEO data, we found that 169 miRNAs were expressed in both KOA serum and cartilage. Following a screening process, miR-199b-5p was selected for further experiments. In cell-based assays, we discovered that miR-199b-5p can influence the viability of chondrocytes and cytokine-mediated extracellular matrix metabolism. Moreover, in vitro experiments demonstrated that it can induce inflammation and abnormal pain threshold in normal mice. Importantly, inhibition of miR-199b- 5p alleviated the pathological symptoms of KOA. Finally, these effects were achieved by target-ing Gcnt2 and Fzd6 (Fig. 7). Thus, our findings demonstrated that miR-199b-5p might be a novel potential therapeutic target for OA prevention and treatment.

miR-199b-5p exerts its effects on in vitro cells and in vivo mice by targeting Fzd6 and Gcnt2.

miRNAs function as regulators of gene expression in biological processes by regulate mRNA translation by specifically binding to the 3′UTRs of target mRNAs 12. Exosomes are rich in miRNAs, and various cells can secrete exosomes to target cells under physiological and patho-logical conditions, which function as delivery vehicles for miRNAs 13;14. There is evidence sug-gesting that miRNAs can participate in various cellular processes (inflammation, cell viability, ECM dysregulation) and signaling pathways (Hedgehog signaling, PI3K-AKT signaling) rele-vant to OA15;16. Our differential miRNA enrichment analysis also strongly supports the correla-tion with these findings.

In cancer research, overexpression of miR-199b-5p can inhibit the proliferation, migration and invasion of prostate cancer cells in vitro and tumor growth and metastasis in vivo by targeting DDR117. In addition, exogenous miR-199b-5p inhibited the growth of bone marrow mesenchy-mal stem cells (MSCs) and promoted the differentiation of bone MSCs into chondrocytes by tar-geting the JAG1 pathway 18;19. Our functional experiments showed that overexpression of miR- 199b-5p reduced chondrocyte viability and the expression of anabolic factors such as COL2A1, AGGREGN, and SOX9. It also increased inflammation levels and decreased pain thresholds in control mice. Although no pathological changes were observed in the articular cartilage of con-trol mice, a similar study demonstrated that pathological cartilage changes only occurred after six months in mice with miR-211 and miR-204 knockout20. Therefore, we speculate that the time of overexpression of miR-199b-5p in our experiment was too short, so it only caused inflammation and pain threshold response. To best of our knowledge, this is the initial investigation document-ing the involvement of miR-199b-5p in KOA.

Fzd6 is known to be up-regulated during the osteogenic differentiation of MSCs and can be regu-lated by miR-194-5p to activate the WNT signaling pathway, thereby promoting osteogenic dif-ferentiation of MSCs21. Gcnt2 has been found to induce epithelial-mesenchymal transition and enhance migration and invasion of esophageal squamous cell carcinoma cells22. Our findings in-dicate that miR-199b-5p plays a crucial role in KOA by targeting Fzd6 and Gcnt2. Future re-search should further explore the roles and mechanisms of Fzd6 and Gcnt2 in the context of KOA.

This study also has some limitations. We initially detected serum exosomes miRNA and later examined the miRNA through in vitro experiment and animal study. However, whether the miRNAs targeting the joints have the same role as the serum exosome miRNAs or blood miRNA has not been clear. We also performed direct intra-articular injection of adenoviral vectors. The intra-articular injection can reduce the exposure of extra-articular tissue, thereby minimizing the associated side effects and improving targeting. Besides, more in-vivo study such as gene knockout mice study can be used in the future study.

In conclusion, we found that miR-199b-5p is elevated in osteoarthritis and may affect cell viabil-ity and related cytokines by targeting Gcnt2 and Fzd6. Overexpression of miR-199b-5p induced OA-like pathological changes in normal mice and inhibiting miR-199b-5p alleviated symptoms in KOA mice, suggesting it may be a target for treatment of OA. Through human clinical trials, cell experiments, and animal models, we not only identified a new OA-related miR-199b-5p but also examined the biological function of miR-199b-5p in OA.

Materials and methods

Human samples

Study participants were recruited from the Hospital of Chengdu University of Traditional Chi-nese Medicine and the surrounding communities. The study was approved by the Ethics Review Committee of the Hospital at Chengdu University of Traditional Chinese Medicine (2016KL-017) and conformed to the ethical guidelines of the 1975 Declaration of Helsinki. Serum was collected from all participants and serum samples were stored at -80°C.

Human subjects

Patients were enrolled if they fulfilled the following criteria: (1) diagnosed with KOA according to the ACR; (2) age between 40 and 65 years; (3) agreed to cooperate with researchers in all re-search procedures after enrollment; (4) provided with written informed consent. Patients with any of the following conditions were excluded: (1) accompanied with other diseases such as rheumatoid arthritis, bone tumors, and bone tuberculosis; (2) treated with intra-articular gluco-corticoid or viscoelastic supplementation in the last six months; (3) knee replacement history; (4) had complicated cardiovascular disease, diabetes, skin disease, and liver or kidney impairment; (5) pregnant or breastfeeding women; (6) accompanied with mental and intellectual disabilities; or (7) undergoing other clinical trials.

Extraction and sequencing of serum exosome miRNAs

Serum samples were filtered using 0.22 µM filters. Exosomes were isolated from the serum sample using ExoQuick™ Exosome Precipitation Solution (System Biosciences) following the manufacturer’s instructions. Briefly, serum was thawed on ice and centrifuged at 4000 rpm for 15 min to remove any cells or cellular debris. Next, 50 μL ExoQuick Solution was added into the 200 μL serum sample and mixed thoroughly. The exosomes were suspended in PBS.

Exosomes were characterized by electron microscopy (Tecnai G2 Spirit 120KV, FEI), nanopar-ticle tracking analysis (NTA) (NTA 3.2 Dev Build 3.2.16), and western blot analysis (for CD9, CD63 and CD81). Sequencing was performed by single-end sequencing (1×150 bp) on Illumina NextSeq 500. The libraries were sequenced on an Agilent 2100 Bioanalyzer platform. The mirdeep2 software (https://www.mdc-berlin.de/content/mirdeep2-documentation) was used to analyze the miRNA sequences and quantification. The heatmap was plotted based on the log2 (fold change), using Heatmap Illustrator software (Heml 1.0).

Total RNA was isolated from human serum, cell and mice

Chondrocytes using a kit (Yeasen,Shanghai, China) according to the manufacturer’s instructions. Reverse transcription was performed using 1000 ng total RNA and a Prime Script RT Reagent Kit(Yeasen, Shanghai, China) or Prime Script RT Master Mix (Yeasen, Shanghai, China), which were used for miRNA and mRNA, respectively. For miRNA, the reactions were incubated at 42°C for 15 min followed by inactivation at 85°C for 5 s. qRT-PCR amplification was assessed in a CFX Connection Real-Time System (Bio-Rad) using the SYBR Premix Ex Taq II kit (Yeasen, Shanghai,China). The following cycling conditions were used: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. All reactions were performed in duplicate and normalized to the internal reference U6 for miRNA and GAPDH mRNA for mRNAs. The 2-ΔΔCT CT method was used to evaluate the relative mRNA/miRNA expression levels.

RNA primer

Primers are listed in Supplementary Table 1.

Bioinformatics analysis

GO analysis was performed, involving the biological process (BP), cellular component (CC), and molecular function (MF), using DAVID23 (https://david.ncifcrf.gov/home.jsp).We used the miRNA target gene prediction websites miRDB24(http://mirdb.org/), miR-Walk25(http://mirwalk.umm.uni-heidelberg.de/), Starbase26 (https://starbase. sysu.edu.cn/), DI-ANA-TarBase27 (http://diana.imis.athena-innovation.gr/DianaTools/ index.php?r=tarbase/index), Targetscan28 (https://www.targetscan.org/vert_72/). The target genes of mmu-miR-199b-5p were predicted, and the intersection of the target gene prediction list results of the five sites was taken.

The dysregulated miRNAs were compared to relevant published miRNA data from human KOA patients. The dataset GSE105027 was obtained from serum samples of KOA patients, while GSE175961 comprised sequencing data from KOA patient cartilage. Finally, we validated the expression of the target genes Gcnt2 and Fzd6 in KOA patients using synovial data (GDS5401, GDS5403) and cartilage data (GDS3758).

Primary mouse chondrocyte culture

Primary mouse cartilage was extracted following the protocol of Gosset et al.29. The cells were cultured in F12 medium containing 10% fetal bovine serum and 1% penicillin-streptomycin at 37°C and a 5% CO2. The medium was changed every three days, and cells were cultured to the second passage for experiments. Cells were stimulated with LPS for 6 h before subsequent inves-tigation.

Animal model of KOA

Animal experiments were performed using total of 108 8-week-old male C57BL/6 mice. The protocol was approved by the Committee on the Ethics of Animal Experiments of Chengdu Uni-versity of Traditional Chinese medicine (2019-04). After randomization, the KOA mouse model was established by orthotopic injection of MIA (sigma) into the knee joint of mice. After the an-imals were anesthetized with isoflurane in a small animal anesthesia machine, the right knee joint of the mouse was shaved, the knee joint was flexed 90°, and 10 μL of MIA solution was injected into the knee joint cavity using a 28G microsyringe. The syringe needle was slowly withdrawn and the knee joint was gently moved. Before model establishment and at 3, 7, 10, 14, 21, and 28 days after model establishment. After the 3rd day and 14th day after the start of the experiment, the experimental mice were injected with adenovirus (total volume 10μL) in the knee joint, and the control mice were injected with empty adenovirus as a control (Hanbio tech, Shanghai, China)30.

Thermal pain threshold detection

A pain threshold detection instrument was used to measure the latency of the right plantar leg raising reflex of mice under heat radiation. Three days before the experiment, mice were cut off from water for half a day and adapted to the environment for 1 h. A certain intensity of pyrogen was used to irradiate the plantar position of the mouse’s right limb, and a machine was used to automatically record the time when the mouse moved the limb from thermal pain. The resting intensity of the thermal pain stimulator was set to 10%, and the maximum duration was set to 20 s to avoid prolonged heating and burning of skin. After the light source is aimed at the sole of the mouse, we observed the mouse in real time. If paw raising, paw licking, or paw retraction was observed, the irradiation was stopped and the irradiation time of the instrument was recorded. The right limb of each mouse was tested five times, with an interval of 10 min each time; outliers were eliminated, and the average value was calculated and included in the statistics31.

ELISA and μCT scanning

Mice were sacrificed after four weeks, and serum was collected for ELISA assay (Jiangsu Jing-mei Biotechnology Co., Ltd., Yancheng, China). The knee joint specimens of the right hindlimb of mice were scanned using a high-resolution micro-CT skysan1267 instrument (Bruker, Germa-ny). The sample was removed from fixative and dried. The scanning parameters were as follows: voltage 55 KV, current 200 μA, and filter 0.25AL. After scanning, three-dimensional reconstruc-tion was performed using NRecon1.7.4.2 software.

Safranin Fast green staining

Mouse knee joint staining was performed using the Safranin Fast Green Staining Kit (Servicebio Biological Technology, Wuhan, China). The sample was examined under a microscope, and the cartilage integrity was scored according to the OARSI grading system, in which the score ranged from 0 to 6 points32.

Luciferase assay

The wild-type (WT) and mutant-type (MUT) sequences (according to the predicted binding site) were inserted into the pmiRGLO plasmid. HEK-293T cells were seeded in 6-well plates at 24 h before transfection. GCNT2 and FZD6 3′UTR-wt and gcnt2 and fzd6 3′UTR-mut plasmids (500 ng) and 20 nmol miR-199b-3p and NC were co-transfected with Lipofectamine 3000 (Hanbio Biotechnology, China) following the manufacturer’s instructions. After 48h, firefly and Renilla luciferase activities were calculated using the Promega Dual-Luciferase system following the manufacturer’s instructions (Hanbio tech, Shanghai, China). Firefly/Renilla luciferase was meas-ured to evaluate relative luciferase activity.

Statistical analysis

Data are reported as the means ± SD. Normal distribution and homogeneity of variance of data were first tested. Shapiro-Wilk test was used to verify data normality, while Levene’s test or Browne-Forsythe test was adopted for assessment of variance equality. Statistical analysis was performed by unpaired two t-test and Tukey-corrected one-way analysis of variance (ANOVA) for comparisons between groups. Tukey-corrected two-way ANOVA for the comparison of mice thermal pain threshold data. P < 0.05 was considered statistically significant for all statistical calculations. PRISM 8.0 (GraphPad Software, San Diego, CA, USA) was used for data analysis.

Contributions

F.T., Z.Q., and Q.-F.W. designed the study.

F.T. performed most in vitro and in vivo experiments.

F.T. and J.-M.W. analyzed data.

Z.Q., Y.J. and S.-H.L. performed the miRNA sequencing and interpreted the data.

Z.Q., and Y.J collected and inspected human patient samples.

F.T. and Q.-F.W. wrote the manuscript.

All authors read and edited the manuscript.

F.-R.L., S.-G.Y., and Q.-F.W. supervised the study.

Competing interests

The authors declare that they have no known competing financial interests or personal relation-ships that could have appeared to influence the work reported in this paper

Acknowledgements

This work was supported by the National Key R&D Program of China (No.2019YFC1709001); Innovation Team and Talents Cultivation Program of National Administration of Traditional Chinese Medicine (No. ZYYCXTD-D-202003); Fund of Science and Technology Department of Sichuan Province, China (No. 2021ZYD0081); and National Natural Science Foundation of Chi-na (No. 82174512, 81873383).

Data Sharing Statement

All data are available in the main text or the supplementary materials.

List of Supplementary Materials

Present a list of the Supplementary Materials in the following format.

Materials and Methods

Fig S1. Identification of exosome. (A) Transmission electron microscope scanning of isolated exosomes from the serum of participants. (B) The Nano sight particle analysis of isolated exo-somes from the serum of participants. (C) The Western blot analysis of symbolic surface markers of isolated exosomes from the serum of participants.

Fig S2. Fluorescence of Adenovirus-Infected Chondrocytes.

Fig S3. Expression of adenovirus in mouse knee joint.

Fig S4. Establishment of KOA model in mice by injection of MIA. (A) Behavioral detection of animal thermal pain threshold. (B) Microscopic observation of the surface of the mouse knee joint. (C) Safranin-fast green staining and semiquantitative scoring of articular cartilage. (D) 3D reconstruction images of joints from μCT scans. *P<0.05, **P<0.01, n=6.

Primer List

Basic information for recruiting patients