Background. The molecular mechanisms underlying non-union bone fractures largely remain elusive. Recently, spatial transcriptomics approaches for musculoskeletal tissue samples have been developed requiring direct placement of histology sections on barcoded slides. However, Formalin-Fixed-Paraffin-Embedded (FFPE) bone sections have been associated with limited RNA quality and read depth compared to soft tissue. Here, we test spatial transcriptomics workflows based on transcriptomic probe transfer to characterize molecular features discriminating non-union and union bone fractures in mice. Method. Histological sections (n=8) used for spatial transcriptomics (Visium CytAssist FFPE; 10x Genomics, n=4 on glass slides, n=4 on hydrogel-coated slides) were obtained from a fracture healing study in female 20-week-old C57BL/6J mice receiving either a femur osteotomy (0.7mm) or a segmental defect (2.4mm) (license 22/2022, Grisons CH). Sequence alignment and manual segmentation of different tissues (bone, defect region/callus, bone marrow, muscle) were performed using SpaceRanger and LoupeBrowser (10x Genomics). Differential gene expression was performed using DESeq2 (Seurat) followed by Gene-Set-Enrichment-Analysis (GSEA) of
Introduction. The most frequent diagnosis in young adults undergoing total hip arthroplasty (THA) is osteonecrosis of the femoral head (ONFH), an evolving and disabling condition with an increasing prevalence worldwide. Treatment of ONFH remains a challenge mainly because of a lack of understanding of the disease's pathophysiological basis. This study investigated the biological processes that could be affected by ONFH by comparing the microstructure, histological characteristics and transcriptomic profile of trabecular bone from the femoral head (FH) and the intertrochanteric region (IT) of patients suffering from this condition. Method. A total of 18 patients with idiopathic ONFH undergoing THA in our institution were included. Trabecular bone explants were taken intraoperatively from the FH and the IT of patients. Bone microstructure was examined by micro-computed tomography (micro-CT). After bone sectioning, histological features were studied by hematoxylin and eosin staining. Differential gene expression was investigated using a microarray platform. Result. Micro-CT imaging showed higher trabecular separation and lower trabecular thickness and bone volume in trabecular bone from the FH than from the IT. Histological staining revealed that the number of osteoblasts on the bone surface and the percentage of empty lacunae were higher in trabecular bone from the FH. Transcriptome analysis identified a differential signature in trabecular bone from the FH compared to the IT. The
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. 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.Aims
Methods
This study aimed to demonstrate the promoting effect of elastic fixation on fracture, and further explore its mechanism at the gene and protein expression levels. A closed tibial fracture model was established using 12 male Japanese white rabbits, and divided into elastic and stiff fixation groups based on different fixation methods. Two weeks after the operation, a radiograph and pathological examination of callus tissue were used to evaluate fracture healing. Then, the differentially expressed proteins (DEPs) were examined in the callus using proteomics. Finally, in vitro cell experiments were conducted to investigate hub proteins involved in this process.Aims
Methods
Rotator cuff tear (RCT) is the leading cause of shoulder pain, primarily associated with age-related tendon degeneration. This study aimed to elucidate the potential differential gene expressions in tendons across different age groups, and to investigate their roles in tendon degeneration. Linear regression and differential expression (DE) analyses were performed on two transcriptome profiling datasets of torn supraspinatus tendons to identify age-related genes. Subsequent functional analyses were conducted on these candidate genes to explore their potential roles in tendon ageing. Additionally, a secondary DE analysis was performed on candidate genes by comparing their expressions between lesioned and normal tendons to explore their correlations with RCTs.Aims
Methods
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. 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.Aims
Methods
This study aimed to explore the biological and clinical importance of dysregulated key genes in osteoarthritis (OA) patients at the cartilage level to find potential biomarkers and targets for diagnosing and treating OA. Six sets of gene expression profiles were obtained from the Gene Expression Omnibus database. Differential expression analysis, weighted gene coexpression network analysis (WGCNA), and multiple machine-learning algorithms were used to screen crucial genes in osteoarthritic cartilage, and genome enrichment and functional annotation analyses were used to decipher the related categories of gene function. Single-sample gene set enrichment analysis was performed to analyze immune cell infiltration. Correlation analysis was used to explore the relationship among the hub genes and immune cells, as well as markers related to articular cartilage degradation and bone mineralization.Aims
Methods
To investigate the effects of senescent osteocytes on bone homeostasis in the progress of age-related osteoporosis and explore the underlying mechanism. In a series of in vitro experiments, we used tert-Butyl hydroperoxide (TBHP) to induce senescence of MLO-Y4 cells successfully, and collected conditioned medium (CM) and senescent MLO-Y4 cell-derived exosomes, which were then applied to MC3T3-E1 cells, separately, to evaluate their effects on osteogenic differentiation. Furthermore, we identified differentially expressed microRNAs (miRNAs) between exosomes from senescent and normal MLO-Y4 cells by high-throughput RNA sequencing. Based on the key miRNAs that were discovered, the underlying mechanism by which senescent osteocytes regulate osteogenic differentiation was explored. Lastly, in the in vivo experiments, the effects of senescent MLO-Y4 cell-derived exosomes on age-related bone loss were evaluated in male SAMP6 mice, which excluded the effects of oestrogen, and the underlying mechanism was confirmed.Aims
Methods
The pathophysiological basis of alterations in trabecular bone of patients with osteonecrosis of the femoral head (ONFH) remains unclear. ONFH has classically been considered a vascular disease with secondary changes in the subchondral bone. However, there is increasing evidence suggesting that ONFH could be a bone disease, since alterations in the functionality of bone tissue distant from the necrotic lesion have been observed. We comparatively studied the transcriptomic profile of trabecular bone obtained from the intertrochanteric region of patients with ONFH without an obvious aetiological factor, and patients with osteoarthritis (OA) undergoing total hip replacement in our Institution. To explore the biological processes that could be affected by ONFH, we compared the transcriptomic profile of trabecular bone from the intertrochanteric region and the femoral head of patients affected by this condition. Differential gene expression was studied using an Affymetrix microarray platform. Transcriptome analysis showed a differential signature in trabecular bone from the intertrochanteric region between patients with ONFH and those with OA. The
Aims. Knee osteoarthritis (OA) involves a variety of tissues in the joint. Gene expression profiles in different tissues are of great importance in order to understand OA. Methods. First, we obtained gene expression profiles of cartilage, synovium, subchondral bone, and meniscus from the Gene Expression Omnibus (GEO). Several datasets were standardized by merging and removing batch effects. Then, we used unsupervised clustering to divide OA into three subtypes. The
Aims. This study aimed, through bioinformatics analysis and in vitro experiment validation, to identify the key extracellular proteins of intervertebral disc degeneration (IDD). Methods. The gene expression profile of GSE23130 was downloaded from the Gene Expression Omnibus (GEO) database. Extracellular protein-differentially expressed genes (EP-DEGs) were screened by protein annotation databases, and we used
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,
Aims. Lumbar spinal stenosis (LSS) is a common skeletal system disease that has been partly attributed to genetic variation. However, the correlation between genetic variation and pathological changes in LSS is insufficient, and it is difficult to provide a reference for the early diagnosis and treatment of the disease. Methods. We conducted a transcriptome-wide association study (TWAS) of spinal canal stenosis by integrating genome-wide association study summary statistics (including 661 cases and 178,065 controls) derived from Biobank Japan, and pre-computed gene expression weights of skeletal muscle and whole blood implemented in FUSION software. To verify the TWAS results, the candidate genes were furthered compared with messenger RNA (mRNA) expression profiles of LSS to screen for common genes. Finally, Metascape software was used to perform enrichment analysis of the candidate genes and common genes. Results. TWAS identified 295 genes with permutation p-values < 0.05 for skeletal muscle and 79 genes associated for the whole blood, such as RCHY1 (PTWAS = 0.001). Those genes were enriched in 112
This study aimed to investigate the role and mechanism of meniscal cell lysate (MCL) in fibroblast-like synoviocytes (FLSs) and osteoarthritis (OA). Meniscus and synovial tissue were collected from 14 patients with and without OA. MCL and FLS proteins were extracted and analyzed by liquid chromatography‒mass spectrometry (LC‒MS). The roles of MCL and adenine nucleotide translocase 3 (ANT3) in FLSs were examined by enzyme-linked immunosorbent assay (ELISA), flow cytometry, immunofluorescence, and transmission electron microscopy. Histological analysis was performed to determine ANT3 expression levels in a male mouse model.Aims
Methods
Osteoarthritis (OA) is a prevalent joint disorder with inflammatory response and cartilage deterioration as its main features. Dihydrocaffeic acid (DHCA), a bioactive component extracted from natural plant ( In vitro, interleukin-1 beta (IL-1β) was used to establish the mice OA chondrocytes. Cell counting kit-8 evaluated chondrocyte viability. Western blotting analyzed the expression levels of collagen II, aggrecan, SOX9, inducible nitric oxide synthase (iNOS), IL-6, matrix metalloproteinases (MMPs: MMP1, MMP3, and MMP13), and signalling molecules associated with nuclear factor-kappa B (NF-κB) and mitogen-activated protein kinase (MAPK) pathways. Immunofluorescence analysis assessed the expression of aggrecan, collagen II, MMP13, and p-P65. In vivo, a destabilized medial meniscus (DMM) surgery was used to induce mice OA knee joints. After injection of DHCA or a vehicle into the injured joints, histological staining gauged the severity of cartilage damage.Aims
Methods
Aims. Osteoporosis (OP) is a metabolic bone disease, characterized by a decrease in bone mineral density (BMD). However, the research of regulatory variants has been limited for BMD. In this study, we aimed to explore novel regulatory genetic variants associated with BMD. Methods. We conducted an integrative analysis of BMD genome-wide association study (GWAS) and regulatory single nucleotide polymorphism (rSNP) annotation information. Firstly, the discovery GWAS dataset and replication GWAS dataset were integrated with rSNP annotation database to obtain BMD associated SNP regulatory elements and SNP regulatory element-target gene (E-G) pairs, respectively. Then, the common genes were further subjected to HumanNet v2 to explore the biological effects. Results. Through discovery and replication integrative analysis for BMD GWAS and rSNP annotation database, we identified 36 common BMD-associated genes for BMD irrespective of regulatory elements, such as FAM3C (p. discovery GWAS. = 1.21 × 10. -25. , p. replication GWAS. = 1.80 × 10. -12. ), CCDC170 (p. discovery GWAS. = 1.23 × 10. -11. , p. replication GWAS. = 3.22 × 10. -9. ), and SOX6 (p. discovery GWAS. = 4.41 × 10. -15. , p. replication GWAS. = 6.57 × 10. -14. ). Then, for the 36 common target genes, multiple
Rheumatoid arthritis (RA) is a common chronic immune disease. Berberine, as its main active ingredient, was also contained in a variety of medicinal plants such as Berberaceae, Buttercup, and Rutaceae, which are widely used in digestive system diseases in traditional Chinese medicine with anti-inflammatory and antibacterial effects. The aims of this article were to explore the therapeutic effect and mechanism of berberine on rheumatoid arthritis. Cell Counting Kit-8 was used to evaluate the effect of berberine on the proliferation of RA fibroblast-like synoviocyte (RA-FLS) cells. The effect of berberine on matrix metalloproteinase (MMP)-1, MMP-3, receptor activator of nuclear factor kappa-Β ligand (RANKL), tumour necrosis factor alpha (TNF-α), and other factors was determined by enzyme-linked immunoassay (ELISA) kit. Transcriptome technology was used to screen related pathways and the potential targets after berberine treatment, which were verified by reverse transcription-polymerase chain reaction (RT-qPCR) and Western blot (WB) technology.Aims
Methods
Aims. Degenerative cervical spondylosis (DCS) is a common musculoskeletal disease that encompasses a wide range of progressive degenerative changes and affects all components of the cervical spine. DCS imposes very large social and economic burdens. However, its genetic basis remains elusive. Methods. Predicted whole-blood and skeletal muscle gene expression and genome-wide association study (GWAS) data from a DCS database were integrated, and functional summary-based imputation (FUSION) software was used on the integrated data. A transcriptome-wide association study (TWAS) was conducted using FUSION software to assess the association between predicted gene expression and DCS risk. The TWAS-identified genes were verified via comparison with differentially expressed genes (DEGs) in DCS RNA expression profiles in the Gene Expression Omnibus (GEO) (Accession Number: GSE153761). The Functional Mapping and Annotation (FUMA) tool for genome-wide association studies and Meta tools were used for gene functional enrichment and annotation analysis. Results. The TWAS detected 420 DCS genes with p < 0.05 in skeletal muscle, such as ribosomal protein S15A (RPS15A) (PTWAS = 0.001), and 110 genes in whole blood, such as selectin L (SELL) (PTWAS = 0.001). Comparison with the DCS RNA expression profile identified 12 common genes, including Apelin Receptor (APLNR) (PTWAS = 0.001, PDEG = 0.025). In total, 148 DCS-enriched
Aims. This study aimed to explore the role of small colony variants (SCVs) of Staphylococcus aureus in intraosseous invasion and colonization in patients with periprosthetic joint infection (PJI). Methods. A PJI diagnosis was made according to the MusculoSkeletal Infection Society (MSIS) for PJI. Bone and tissue samples were collected intraoperatively and the intracellular invasion and intraosseous colonization were detected. Transcriptomics of PJI samples were analyzed and verified by polymerase chain reaction (PCR). Results. SCVs can be isolated from samples collected from chronic PJIs intraoperatively. Transmission electron microscopy (TEM) and immunofluorescence (IF) showed that there was more S. aureus in bone samples collected from chronic PJIs, but much less in bone samples from acute PJIs, providing a potential mechanism of PJI. Immunofluorescence results showed that SCVs of S. aureus were more likely to invade osteoblasts in vitro. Furthermore, TEM and IF also demonstrated that SCVs of S. aureus were more likely to invade and colonize in vivo. Cluster analysis and principal component analysis (PCA) showed that there were substantial differences in gene expression profiles between chronic and acute PJI.
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.
