Objectives. Previously, we reported the improved transfection efficiency of a plasmid DNA-chitosan (pDNA-CS) complex using a phosphorylatable nuclear localization signal-linked nucleic kinase substrate short peptide (pNNS) conjugated to chitosan (pNNS-CS). This study investigated the effects of pNNS-CS-mediated miR-140 and interleukin-1 receptor antagonist protein (IL-1Ra) gene transfection both in rabbit chondrocytes and a cartilage defect model. Methods. The pBudCE4.1-miR-140, pBudCE4.1-IL-1Ra, and negative control pBudCE4.1 plasmids were constructed and combined with pNNS-CS to form pDNA/pNNS-CS complexes. These complexes were transfected into chondrocytes or injected into the knee joint cavity. Results. High IL-1Ra and miR-140 expression levels were detected both in vitro and in vivo. In vitro, compared with the pBudCE4.1 group, the transgenic group presented with significantly increased chondrocyte proliferation and glycosaminoglycan (GAG) synthesis, as well as increased collagen type II alpha 1 chain (COL2A1), aggrecan (ACAN), and TIMP metallopeptidase inhibitor 1 (TIMP-1) levels. Nitric oxide (NO) synthesis was reduced, as were a disintegrin and metalloproteinase with thrombospondin type 1 motif 5 (ADAMTS-5) and matrix metalloproteinase (MMP)-13 levels. In vivo, the exogenous genes reduced the synovial fluid GAG and NO concentrations and the ADAMTS-5 and MMP-13 levels in cartilage. In contrast, COL2A1, ACAN, and TIMP-1 levels were increased, and the cartilage Mankin score was decreased in the transgenic group compared with the pBudCE4.1 group. Double gene combination produced greater efficacies than each single gene, both in vitro and in vivo. Conclusion. This study suggests that pNNS-CS is a good candidate for treating cartilage defects via gene therapy, and that IL-1Ra in combination with miR-140 produces promising biological effects on cartilage defects. Cite this article: R. Zhao, S. Wang, L. Jia, Q. Li, J. Qiao, X. Peng. Interleukin-1 receptor antagonist protein (IL-1Ra) and miR-140 overexpression via pNNS-conjugated chitosan-mediated gene transfer enhances the repair of full-thickness cartilage defects in a rabbit model. Bone Joint Res 2019;8:165–178. DOI: 10.1302/2046-3758.83.BJR-2018-0222.R1

Objectives. The objective of this study was to investigate the association of four single-nucleotide polymorphisms (SNPs) of the cannabinoid receptor 2 (CNR2) gene, gene-obesity interaction, and haplotype combination with osteoporosis (OP) susceptibility. Methods. Chinese patients with OP were recruited between March 2011 and December 2015 from our hospital. In this study, a total of 1267 post-menopausal female patients (631 OP patients and 636 control patients) were selected. The mean age of all subjects was 69.2 years (sd 15.8). A generalized multifactor dimensionality reduction (GMDR) model and logistic regression model were used to examine the interaction between SNP and obesity on OP. For OP patient-control haplotype analyses, the SHEsis online haplotype analysis software (. http://analysis.bio-x.cn/. ) was employed. Results. The logistic regression model revealed that the C allele of rs2501431 and the G allele of rs3003336 were associated with increased OP risk, compared with those with wild genotype. However, no significant correlations were found when analyzing the association of rs4237 and rs2229579 with OP risk. The GMDR analysis suggested that the interaction model composed of two factors, rs3003336 and abdominal obesity (AO), was the best model with statistical significance (p-value from sign test (P. sign. ) = 0.012), indicating a potential gene-environment interaction between rs3003336 and AO. Overall, the two-locus models had a cross-validation consistency of 10/10 and had a testing accuracy of 0.641. Abdominally obese subjects with the AG or GG genotype have the highest OP risk, compared with subjects with the AA genotype and normal waist circumference (WC) (odds ratio (OR) 2.23, 95% confidence interval (CI) 1.54 to 3.51). Haplotype analysis also indicated that the haplotype containing the rs3003336-G and rs2501431-C alleles was associated with a statistically increased OP risk. Conclusion. Our results suggested that the C allele of rs2501431 and the G allele of rs3003336 of the CNR2 gene, interaction between rs3003336 and AO, and the haplotype containing the rs3003336-G and rs2501431-C alleles were all associated with increased OP risk. Cite this article: Bone Joint Res 2019;8:544–549

Objectives. Previous genome-wide association studies (GWAS) have reported significant association of the single nucleotide polymorphism (SNP) rs8044769 in the fat mass and obesity-associated gene (FTO) with osteoarthritis (OA) risk in European populations. However, these findings have not been confirmed in Chinese populations. Methods. We systematically genotyped rs8044769 and evaluated the association between the genetic variants and OA risk in a case-controlled study including 196 OA cases and 442 controls in a northern Chinese population. Genotyping was performed using the Sequenom MassARRAY iPLEX platform. Results. We found that the variant T allele of rs8044769 showed no significant association of OA risk (p = 0.791), or association with body mass index (BMI) (pmeta = 0.786) in an additive genetic model. However, we detected a significant interaction between rs8044769 genotypes and BMI on OA risk (p = 0.037), as well as a borderline interaction between rs8044769 genotypes and age on OA risk (p = 0.062). Conclusions. Our findings indicate that rs8044769 in the FTO gene may not modify individual susceptibility to OA or increased BMI in the Chinese population. Further studies are warranted to validate and extend our findings. Cite this article: Prof L. Jiang. No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: A population-based study in China. Bone Joint Res 2016;5:169–174. DOI: 10.1302/2046-3758.55.2000589

Low bone mass and osteopenia have been described in the axial and peripheral skeleton of patients with adolescent idiopathic scoliosis (AIS). Recently, many studies have shown that gene polymorphism is related to osteoporosis. However, no studies have linked the association between IL6 gene polymorphism and bone mass in AIS. This study examined the association between bone mass and IL6 gene polymorphism in 198 girls with AIS. The polymorphisms of IL6-597 G→A, IL6-572 G→C and IL6-174 G→A and the bone mineral density in the lumbar spine and femoral neck were analysed and compared with their levels in healthy controls. The mean bone mineral density at both sites in patients with AIS was decreased compared with controls (p = 0.0022 and p = 0.0013, respectively). Comparison of genotype frequencies between AIS and healthy controls revealed a statistically significant difference in IL6-572 G→C polymorphism (p = 0.0305). There was a significant association between the IL6-572 G→C polymorphism and bone mineral density in the lumbar spine, with the CC genotype significantly higher with the GC (p = 0.0124) or GG (p = 0.0066) genotypes. These results suggest that the IL6-572 G→C polymorphism is associated with bone mineral density in the lumbar spine in Korean girls with AIS

Objectives. Effects of insulin-like growth factor 1 (IGF1), fibroblast growth factor 2 (FGF2) and bone morphogenetic protein 2 (BMP2) on the expression of genes involved in the proliferation and differentiation of osteoblasts in culture were analysed. The best sequence of growth factor addition that induces expansion of cells before their differentiation was sought. Methods. Primary human osteoblasts in in vitro culture were treated with IGF1, BMP2 or FGF2 (10 ng/ml) for 24 hours (IGF1) or 48 hours (BMP2 and FGF2). Experiments were performed during the exponential growth phase with approximately 1e7 cells per 75 cm. 2. flask. mRNA was reverse transcribed directly and analysed using RT-PCR Taqman assays. Expression levels of key genes involved in cell growth and differentiation (CDH11, TNFRSF11B, RUNX2, POSTN, ALP, WNT5A, LEF1, HSPA5, FOS, p21) were monitored using RT-PCR with gene-specific Taqman probes. . Results. Autocrine expression of BMP2 is stimulated by FGF2 and BMP2 itself. BMP2 and FGF2 act as proliferative factors as indicated by reduced expression of ALP and POSTN, whereas IGF1 exhibits a more subtle picture: the Wingless und Int-1 (Wnt) signalling pathway and the Smad pathway, but not p38 mitogen-activated protein (MAP) kinase signalling, were shown to be activated by IGF1, leading to proliferation and differentiation of the cells. . Conclusions. For future use of autologous bone cells in the management of bony defects, new treatment options take advantage of growth factors and differentiation factors. Thus, our results might help to guide the timely application of these factors for the expansion and subsequent differentiation of osteoblastic cells in culture. Cite this article: Bone Joint Res 2014;3:236–40

Objectives. This study aimed to investigate time-dependent gene expression of injured human anterior cruciate ligament (ACL), and to evaluate the histological changes of the ACL remnant in terms of cellular characterisation. Methods. Injured human ACL tissues were harvested from 105 patients undergoing primary ACL reconstruction and divided into four phases based on the period from injury to surgery. Phase I was < three weeks, phase II was three to eight weeks, phase III was eight to 20 weeks, and phase IV was ≥ 21 weeks. Gene expressions of these tissues were analysed in each phase by quantitative real-time polymerase chain reaction using selected markers (collagen types 1 and 3, biglycan, decorin, α-smooth muscle actin, IL-6, TGF-β1, MMP-1, MMP-2 and TIMP-1). Immunohistochemical staining was also performed using primary antibodies against CD68, CD55, Stat3 and phosphorylated-Stat3 (P-Stat3). . Results. Expression of IL-6 was mainly seen in phases I, II and III, collagen type 1 in phase II, MMP-1, 2 in phase III, and decorin, TGF-β1 and α-smooth muscle actin in phase IV. Histologically, degradation and scar formation were seen in the ACL remnant after phase III. The numbers of CD55 and P-Stat3 positive cells were elevated from phase II to phase III. . Conclusions. Elevated cell numbers including P-Stat3 positive cells were not related to collagens but to MMPs’ expressions

We stably transfected early passage chondrocytes with an anti-apoptotic Bcl-2 gene in vitro using a retrovirus vector. Samples of articular cartilage were obtained from 11 patients with a mean age of 69 years (61 to 75) who were undergoing total knee replacement for osteoarthritis. The Bcl-2-gene-transfected chondrocytes were compared with non-transfected and lac-Z-gene-transfected chondrocytes, both of which were used as controls. All three groups of cultured chondrocytes were incubated with nitric oxide (NO) for ten days. Using the Trypan Blue exclusion assay, an enzyme-linked immunosorbent assay and flow cytometric analysis, we found that the number of apoptotic chondrocytes was significantly higher in the non-transfected and lac-Z-transfected groups than in the Bcl-2-transfected group (p < 0.05). The Bcl-2-transfected chondrocytes were protected from NO-induced impairment of proteoglycan synthesis. We conclude that NO-induced chondrocyte death involves a mechanism which appears to be subject to regulation by an anti-apoptotic Bcl-2 gene. Therefore, Bcl-2 gene therapy may prove to be of therapeutic value in protecting human articular chondrocytes

Objectives. Degenerative disc disease (DDD) and osteoarthritis (OA) are relatively frequent causes of disability amongst the elderly; they constitute serious socioeconomic costs and significantly impair quality of life. Previous studies to date have found that aggrecan variable number of tandem repeats (VNTR) contributes both to DDD and OA. However, current data are not consistent across studies. The purpose of this study was to evaluate systematically the relationship between aggrecan VNTR, and DDD and/or OA. Methods. This study used a highly sensitive search strategy to identify all published studies related to the relationship between aggrecan VNTR and both DDD and OA in multiple databases from January 1996 to December 2016. All identified studies were systematically evaluated using specific inclusion and exclusion criteria. Cochrane methodology was also applied to the results of this study. Results. The final selection of seven studies was comprehensively evaluated and includes results for 2928 alleles. The most frequent allele among all the studies was allele 27. After comparing the distributions of each allele with others, statistically significant differences have been found in the distribution of the alleles by the two groups, with an over-representation of allele (A)21 (disease: 3.22%, control: 0.44%). Thus, carrying A21 increased the risk of DDD. Such an association was not found to be statistically significant when considering the risk of OA. Conclusions. The findings suggest that VNTR A21 seems to be associated with higher risk to DDD, however, such an association may not be statistically significant regarding the risk of OA. Cite this article: L. Cong, G. Tu, D. Liang. A systematic review of the relationship between the distributions of aggrecan gene VNTR polymorphism and degenerative disc disease/osteoarthritis. Bone Joint Res 2018;7:308–317. DOI: 10.1302/2046-3758.74.BJR-2017-0207.R1

We have examined whether primary human muscle-derived cells can be used in ex vivo gene therapy to deliver BMP-2 and to produce bone in vivo. Two in vitro experiments and one in vivo experiment were used to determine the osteocompetence and BMP-2 secretion capacity of cells isolated from human skeletal muscle. We isolated five different populations of primary muscle cells from human skeletal muscle in three patients. In the first in vitro experiment, production of alkaline phosphatase by the cells in response to stimulation by rhBMP-2 was measured and used as an indicator of cellular osteocompetence. In the second, secretion of BMP-2 was measured after the cell populations had been transduced by an adenovirus encoding for BMP-2. In the in vivo experiment, the cells were cotransduced with a retrovirus encoding for a nuclear localised β-galactosidase gene and an adenovirus encoding for BMP-2. The cotransduced cells were then injected into the hind limbs of severe combined immune-deficient (SCID) mice and analysed radiographically and histologically. The nuclear localised β-galactosidase gene allowed identification of the injected cells in histological specimens. In the first in vitro experiment, the five different cell populations all responded to in vitro stimulation of rhBMP-2 by producing higher levels of alkaline phosphatase when compared with non-stimulated cells. In the second, the five different cell populations were all successfully transduced by an adenovirus to express and secrete BMP-2. The cells secreted between 444 and 2551 ng of BMP-2 over three days. In the in vivo experiment, injection of the transduced cells into the hind-limb musculature of SCID mice resulted in the formation of ectopic bone at 1, 2, 3 and 4 weeks after injection. Retroviral labelling of the cell nuclei showed labelled human muscle-derived cells occupying locations of osteoblasts in the ectopic bone, further supporting their osteocompetence. Cells from human skeletal muscle, because of their availability to orthopaedic surgeons, their osteocompetence, and their ability to express BMP-2 after genetic engineering, are an attractive cell population for use in BMP-2 gene therapy approaches

We have examined the process of fusion of the intertransverse processes and bone graft in the rabbit by in situ hybridisation and evaluated the spatial and temporal expression of genes encoding pro-α1 (I) collagen (COL1A1), pro-α1 (II) collagen (COL2A1) and pro-α1 (X) collagen (COL10A1). Beginning at two weeks after operation, osteogenesis and chondrogenesis occurred around the transverse process and the grafted bone at the central portion of the area of the fusion mass. Osteoblasts and osteocytes at the newly-formed woven bone expressed COL1A1. At the cartilage, most chondrocytes expressed COL2A1 and some hypertrophic chondrocytes COL10A1. In some regions, co-expression of COL1A1 and COL2A1 was observed. At four weeks, such expressions for COL1A1, COL2A1 and COL10A1 became prominent at the area of the fusion mass. From four to six weeks, bone remodelling progressed from the area of the transverse processes towards the central zone. Osteoblasts lining the trabeculae expressed a strong signal for COL1A1. At the central portion of the area of the fusion mass, endochondral ossification progressed and chondrocytes expressed COL2A1 and COL10A1. Our findings show that the fusion process begins with the synthesis of collagens around the transverse processes and around the grafted bone independently. Various spatial and temporal osteogenic and chondrogenic responses, including intramembranous, endochondral and transchondroid bone formation, progress after bone grafting at the intertransverse processes. Bone formation through cartilage may play an important role in posterolateral spinal fusion

It has been suggested that matrix metalloproteinase-3 (MMP-3, stromelysin-1) has an important role in the degeneration of intervertebral discs (IVDs). A human MMP-3 promoter 5A/6A polymorphism was reported to be involved in the regulation of MMP-3 gene expression. We suggest that IVD degeneration is associated with 5A/6A polymorphism. We studied 54 young and 49 elderly Japanese subjects. Degeneration of the lumbar discs was graded using MRI in the younger group and by radiography in the elderly. 5A/6A polymorphism was determined by polymerase-chain reaction-based assays. We found that the 5A5A and 5A6A genotype in the elderly was associated with a significantly larger number of degenerative IVDs than the 6A6A (p < 0.05), but there was no significant difference in the young. In the elderly, the IVD degenerative scores were also distributed more highly in the 5A5A and 5A6A genotypes (p = 0.0029). Our findings indicate that the 5A allele is a possible risk factor for the acceleration of degenerative changes in the lumbar disc in the elderly

Aims. This study examined the relationship between obesity (OB) and osteoporosis (OP), aiming to identify shared genetic markers and molecular mechanisms to facilitate the development of therapies that target both conditions simultaneously. Methods. Using weighted gene co-expression network analysis (WGCNA), we analyzed datasets from the Gene Expression Omnibus (GEO) database to identify co-expressed gene modules in OB and OP. These modules underwent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and protein-protein interaction analysis to discover Hub genes. Machine learning refined the gene selection, with further validation using additional datasets. Single-cell analysis emphasized specific cell subpopulations, and enzyme-linked immunosorbent assay (ELISA), protein blotting, and cellular staining were used to investigate key genes. Results. WGCNA revealed critical gene modules for OB and OP, identifying the Toll-like receptor (TLR) signalling pathway as a common factor. TLR2 was the most significant gene, with a pronounced expression in macrophages. Elevated TLR2 expression correlated with increased adipose accumulation, inflammation, and osteoclast differentiation, linking it to OP development. Conclusion. Our study underscores the pivotal role of TLR2 in connecting OP and OB. It highlights the influence of TLR2 in macrophages, driving both diseases through a pro-inflammatory mechanism. These insights propose TLR2 as a potential dual therapeutic target for treating OP and OB. Cite this article: Bone Joint Res 2024;13(10):573–587

Aims. Osteoarthritis (OA) is a common degenerative joint disease. The osteocyte transcriptome is highly relevant to osteocyte biology. This study aimed to explore the osteocyte transcriptome in subchondral bone affected by OA. Methods. Gene expression profiles of OA subchondral bone were used to identify disease-relevant genes and signalling pathways. RNA-sequencing data of a bone loading model were used to identify the loading-responsive gene set. Weighted gene co-expression network analysis (WGCNA) was employed to develop the osteocyte mechanics-responsive gene signature. Results. A group of 77 persistent genes that are highly relevant to extracellular matrix (ECM) biology and bone remodelling signalling were identified in OA subchondral lesions. A loading responsive gene set, including 446 principal genes, was highly enriched in OA medial tibial plateaus compared to lateral tibial plateaus. Of this gene set, a total of 223 genes were identified as the main contributors that were strongly associated with osteocyte functions and signalling pathways, such as ECM modelling, axon guidance, Hippo, Wnt, and transforming growth factor beta (TGF-β) signalling pathways. We limited the loading-responsive genes obtained via the osteocyte transcriptome signature to identify a subgroup of genes that are highly relevant to osteocytes, as the mechanics-responsive osteocyte signature in OA. Based on WGCNA, we found that this signature was highly co-expressed and identified three clusters, including early, late, and persistently responsive genes. Conclusion. In this study, we identified the mechanics-responsive osteocyte signature in OA-lesioned subchondral bone. Cite this article: Bone Joint Res 2022;11(6):362–370

Aims. To explore the novel molecular mechanisms of histone deacetylase 4 (HDAC4) in chondrocytes via RNA sequencing (RNA-seq) analysis. Methods. Empty adenovirus (EP) and a HDAC4 overexpression adenovirus were transfected into cultured human chondrocytes. The cell survival rate was examined by real-time cell analysis (RTCA) and EdU and flow cytometry assays. Cell biofunction was detected by Western blotting. The expression profiles of messenger RNAs (mRNAs) in the EP and HDAC4 transfection groups were assessed using whole-transcriptome sequencing (RNA-seq). Volcano plot, Gene Ontology, and pathway analyses were performed to identify differentially expressed genes (DEGs). For verification of the results, the A289E/S246/467/632 A sites of HDAC4 were mutated to enhance the function of HDAC4 by increasing HDAC4 expression in the nucleus. RNA-seq was performed to identify the molecular mechanism of HDAC4 in chondrocytes. Finally, the top ten DEGs associated with ribosomes were verified by quantitative polymerase chain reaction (QPCR) in chondrocytes, and the top gene was verified both in vitro and in vivo. Results. HDAC4 markedly improved the survival rate and biofunction of chondrocytes. RNA-seq analysis of the EP and HDAC4 groups showed that HDAC4 induced 2,668 significant gene expression changes in chondrocytes (1,483 genes upregulated and 1,185 genes downregulated, p < 0.05), and ribosomes exhibited especially large increases. The results were confirmed by RNA-seq of the EP versus mutated HDAC4 groups and the validations in vitro and in vivo. Conclusion. The enhanced ribosome pathway plays a key role in the mechanism by which HDAC4 improves the survival rate and biofunction of chondrocytes. Cite this article: Bone Joint Res 2023;12(7):433–446

Aims. The antidiabetic agent metformin inhibits fibrosis in various organs. This study aims to elucidate the effects of hyperglycaemia and metformin on knee joint capsule fibrosis in mice. Methods. Eight-week-old wild-type (WT) and type 2 diabetic (db/db) mice were divided into four groups without or with metformin treatment (WT met(-/+), Db met(-/+)). Mice received daily intraperitoneal administration of metformin and were killed at 12 and 14 weeks of age. Fibrosis morphology and its related genes and proteins were evaluated. Fibroblasts were extracted from the capsules of 14-week-old mice, and the expression of fibrosis-related genes in response to glucose and metformin was evaluated in vitro. Results. The expression of all fibrosis-related genes was higher in Db met(-) than in WT met(-) and was suppressed by metformin. Increased levels of fibrosis-related genes, posterior capsule thickness, and collagen density were observed in the capsules of db/db mice compared with those in WT mice; these effects were suppressed by metformin. Glucose addition increased fibrosis-related gene expression in both groups of mice in vitro. When glucose was added, metformin inhibited the expression of fibrosis-related genes other than cellular communication network factor 2 (Ccn2) in WT mouse cells. Conclusion. Hyperglycaemia promotes fibrosis in the mouse knee joint capsule, which is inhibited by metformin. These findings can help inform the development of novel strategies for treating knee joint capsule fibrosis. Cite this article: Bone Joint Res 2024;13(7):321–331

Aims. Currently, the effect of drug treatment for osteoporosis is relatively poor, and the side effects are numerous and serious. Melatonin is a potential drug to improve bone mass in postmenopausal women. Unfortunately, the mechanism by which melatonin improves bone metabolism remains unclear. The aim of this study was to further investigate the potential mechanism of melatonin in the treatment of osteoporosis. Methods. The effects of melatonin on mitochondrial apoptosis protein, bmal1 gene, and related pathway proteins of RAW264.7 (mouse mononuclear macrophage leukaemia cells) were analyzed by western blot. Cell Counting Kit-8 was used to evaluate the effect of melatonin on cell viability. Flow cytometry was used to evaluate the effect of melatonin on the apoptosis of RAW264.7 cells and mitochondrial membrane potential. A reactive oxygen species (ROS) detection kit was used to evaluate the level of ROS in osteoclast precursors. We used bmal1-small interfering RNAs (siRNAs) to downregulate the Bmal1 gene. We established a postmenopausal mouse model and verified the effect of melatonin on the bone mass of postmenopausal osteoporosis in mice via micro-CT. Bmal1 lentiviral activation particles were used to establish an in vitro model of overexpression of the bmal1 gene. Results. Melatonin promoted apoptosis of RAW264.7 cells and increased the expression of BMAL1 to inhibit the activation of ROS and phosphorylation of mitogen-activated protein kinase (MAPK)-p38. Silencing the bmal1 gene weakened the above effects of melatonin. After that, we used dehydrocorydaline (DHC) to enhance the activation of MAPK-p38, and the effects of melatonin on reducing ROS levels and promoting apoptosis of RAW264.7 cells were also blocked. Then, we constructed a mouse model of postmenopausal osteoporosis and administered melatonin. The results showed that melatonin improves bone loss in ovariectomized mice. Finally, we established a model of overexpression of the bmal1 gene, and these results suggest that the bmal1 gene can regulate ROS activity and change the level of the MAPK-p38 signalling pathway. Conclusion. Our study confirmed that melatonin promotes the apoptosis of RAW264.7 cells through BMAL1/ROS/MAPK-p38, and revealed the therapeutic effect and mechanism of melatonin in postmenopausal osteoporosis. This finding enriches BMAL1 as a potential target for the treatment of osteoporosis and the pathogenesis of postmenopausal osteoporosis. Cite this article: Bone Joint Res 2023;12(11):677–690

Aims. Osteoarthritis (OA) is a common degenerative joint disease worldwide, which is characterized by articular cartilage lesions. With more understanding of the disease, OA is considered to be a disorder of the whole joint. However, molecular communication within and between tissues during the disease process is still unclear. In this study, we used transcriptome data to reveal crosstalk between different tissues in OA. Methods. We used four groups of transcription profiles acquired from the Gene Expression Omnibus database, including articular cartilage, meniscus, synovium, and subchondral bone, to screen differentially expressed genes during OA. Potential crosstalk between tissues was depicted by ligand-receptor pairs. Results. During OA, there were 626, 97, 1,060, and 2,330 differentially expressed genes in articular cartilage, meniscus, synovium, and subchondral bone, respectively. Gene Ontology enrichment revealed that these genes were enriched in extracellular matrix and structure organization, ossification, neutrophil degranulation, and activation at different degrees. Through ligand-receptor pairing and proteome of OA synovial fluid, we predicted ligand-receptor interactions and constructed a crosstalk atlas of the whole joint. Several interactions were reproduced by transwell experiment in chondrocytes and synovial cells, including TNC-NT5E, TNC-SDC4, FN1-ITGA5, and FN1-NT5E. After lipopolysaccharide (LPS) or interleukin (IL)-1β stimulation, the ligand expression of chondrocytes and synovial cells was upregulated, and corresponding receptors of co-culture cells were also upregulated. Conclusion. Each tissue displayed a different expression pattern in transcriptome, demonstrating their specific roles in OA. We highlighted tissue molecular crosstalk through ligand-receptor pairs in OA pathophysiology, and generated a crosstalk atlas. Strategies to interfere with these candidate ligands and receptors may help to discover molecular targets for future OA therapy. Cite this article: Bone Joint Res 2022;11(12):862–872

Aims. To assess the alterations in cell-specific DNA methylation associated with chondroitin sulphate response using peripheral blood collected from Kashin-Beck disease (KBD) patients before initiation of chondroitin sulphate treatment. Methods. Peripheral blood samples were collected from KBD patients at baseline of chondroitin sulphate treatment. Methylation profiles were generated using reduced representation bisulphite sequencing (RRBS) from peripheral blood. Differentially methylated regions (DMRs) were identified using MethylKit, while DMR-related genes were defined as those annotated to the gene body or 2.2-kilobase upstream regions of DMRs. Selected DMR-related genes were further validated by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) to assess expression levels. Tensor composition analysis was performed to identify cell-specific differential DNA methylation from bulk tissue. Results. This study revealed 21,060 hypermethylated and 44,472 hypomethylated DMRs, and 13,194 hypermethylated and 22,448 hypomethylated CpG islands for differential global methylation for chondroitin sulphate treatment response. A total of 12,666 DMR-related genes containing DMRs were identified in their promoter regions, such as CHL1 (false discovery rate (FDR) = 2.11 × 10. -11. ), RIC8A (FDR = 7.05 × 10. -4. ), and SOX12 (FDR = 1.43 × 10. -3. ). Additionally, RIC8A and CHL1 were hypermethylated in responders, while SOX12 was hypomethylated in responders, all showing decreased gene expression. The patterns of cell-specific differential global methylation associated with chondroitin sulphate response were observed. Specifically, we found that DMRs located in TESPA1 and ATP11A exhibited differential DNA methylation between responders and non-responders in granulocytes, monocytes, and B cells. Conclusion. Our study identified cell-specific changes in DNA methylation associated with chondroitin sulphate response in KBD patients. Cite this article: Bone Joint Res 2024;13(5):237–246