Identification of potential diagnostic gene biomarkers in patients with osteoarthritis

The Actin Cytoskeleton as a Regulator of Proteoglycan 4

OBJECTIVE

The superficial zone (SZ) of articular cartilage is responsible for distributing shear forces for optimal cartilage loading and contributes to joint lubrication through the production of proteoglycan 4 (PRG4). PRG4 plays a critical role in joint homeostasis and is chondroprotective. Normal PRG4 production is affected by inflammation and irregular mechanical loading in post-traumatic osteoarthritis (PTOA). THe SZ chondrocyte (SZC) phenotype, including PRG4 expression, is regulated by the actin cytoskeleton in vitro. There remains a limited understanding of the regulation of PRG4 by the actin cytoskeleton in native articular chondrocytes. The filamentous (F)-actin cytoskeleton is a potential node in crosstalk between mechanical stimulation and cytokine activation and the regulation of PRG4 in SZCs, therefore developing insights in the regulation of PRG4 by actin may identify molecular targets for novel PTOA therapies.

MATERIALS AND METHODS

A comprehensive literature search on PRG4 and the regulation of the SZC phenotype by actin organization was performed.

RESULTS

PRG4 is strongly regulated by the actin cytoskeleton in isolated SZCs in vitro. Biochemical and mechanical stimuli have been characterized to regulate PRG4 and may converge upon actin cytoskeleton signaling.

CONCLUSION

Actin-based regulation of PRG4 in native SZCs is not fully understood and requires further elucidation. Understanding the regulation of PRG4 by actin in SZCs requires an in vivo context to further potential of leveraging actin arrangement to arthritic therapeutics.

Identification and validation of FPR1, FPR2, IL17RA and TLR7 as immunogenic cell death related genes in osteoarthritis

Immunogenic cell death (ICDs) has gained increasing attention for its significant clinical efficacy in various diseases. Similarly, more and more attention has been paid in the role of immune factors in the pathological process of osteoarthritis (OA). The objective of this study is to reveal the relationship between ICD-related genes and the process of OA at the gene level through bioinformatics analysis. In this study, Limma R package was applied to identify differentially expressed genes (DEG), and OA related module genes were determined by weighted gene co-expression network analysis. The ICD-related genes were extracted from a previous study. The module genes related to DEGs and ICD were overlapped. Then, hub genes were identified by a series of analyses using the Least absolute shrinkage and selection operator and random forest algorithm, the expression level and diagnostic value of hub genes were evaluated by Logistic regression. In addition, we used Spearman rank correlation analysis to clarify the relationship between hub genes and infiltrating immune cells and immune pathways. The expression levels of FPR1, FPR2, IL17RA, and TLR7 was verified in SD rat knee joint model of OA by immunohistochemistry. The expression levels of FPR1, FPR2, IL17RA, and TLR7 mRNA were detected in the IL-1β induced rat chondrocytes in qPCR experiment in vitro. Four hub genes (FPR1, FPR2, IL17RA, and TLR7) were ultimately identified as OA biomarkers associated with ICD. And knockdown of TLR7 reversed collagen II and ADAMTS-5 degradation in IL-1β-stimulated chondrocytes. This research may provide new immune related biomarkers for the diagnosis of OA and serve as a reference for disease treatment monitoring.

Synovial Fluid Proteomics From Serial Aspirations of ACL-Injured Knees Identifies Candidate Biomarkers

BACKGROUND

Anterior cruciate ligament (ACL) tears often result in knee effusion and an increased risk for developing knee osteoarthritis (OA) in the long run. The molecular profile of these effusions could be informative regarding initial steps in the development of posttraumatic OA after an ACL tear.

HYPOTHESIS

The proteomics of knee synovial fluid changes over time after ACL injury.

STUDY DESIGN

Descriptive laboratory study.

METHODS

Synovial fluid was collected from patients with an acute traumatic ACL tear presenting to the office for evaluation (18.31 ± 19.07 days from injury) (aspiration 1) and again at the time of surgery (35.41 ± 58.15 days after aspiration 1 (aspiration 2). High-resolution liquid chromatography mass spectrometry was used to assess the quantitative protein profile of synovial fluid, and differences in protein profile between the 2 aspirations were determined computationally.

RESULTS

A total of 58 synovial fluid samples collected from 29 patients (12 male, 17 female; 12 isolated ACL tear, 17 combined ACL and meniscal tear) with a mean age and body mass index of 27.01 ± 12.78 years and 26.30 ± 4.93, respectively, underwent unbiased proteomics analysis. The levels of 130 proteins in the synovial fluid changed over time (87 high, 43 low). Proteins of interest that were significantly higher in aspiration 2 included CRIP1, S100A11, PLS3, POSTN, and VIM, which represent catabolic/inflammatory activities in the joint. Proteins with a known role in chondroprotection and joint homeostasis such as CHI3L2 (YKL-39), TNFAIP6/TSG6, DEFA1, SPP1, and CILP were lower in aspiration 2.

CONCLUSION

Synovial fluid from knees with ACL tears exhibits an increased burden of inflammatory (catabolic) proteins relevant to OA with reduced levels of chondroprotective (anabolic) proteins.

CLINICAL RELEVANCE

This study identified a set of novel proteins that provide new biological insights into the aftermath of ACL tears. Elevated inflammation and decreased chondroprotection could represent initial disruption of homeostasis, potentially initiating the development of OA. Longitudinal follow-up and mechanistic studies are necessary to assess the functional role of these proteins in the joint. Ultimately, these investigations could lead to better approaches to predict and possibly improve patient outcomes.

Genome-wide DNA methylation and transcription analysis in psoriatic epidermis

Aim: To identify DNA methylation and transcription biomarkers in the psoriatic epidermis. Materials & methods: Gene transcription and DNA methylation datasets of psoriatic epidermal tissue were obtained from the Gene Expression Omnibus. Machine learning algorithm analysis and weighted gene coexpression network analysis were carried out to screen hub genes. Results: Differentially methylated and expressed genes were identified in the psoriatic epidermis. Six hub genes were selected - GZMB, CRIP1, S100A12, ISG15, CRABP2 and VNN1 - whose transcript levels showed a significant correlation with Psoriasis Area and Severity Index scores and immune infiltration. Conclusion: Psoriatic epidermis is primarily in a hypermethylated status. Epidermis-specific hub differentially methylated and expressed genes are potential biomarkers to help judge the condition of psoriasis.

Development and validation of a gene signature predicting the risk of postmenopausal osteoporosis

AIMS

We aimed to develop a gene signature that predicts the occurrence of postmenopausal osteoporosis (PMOP) by studying its genetic mechanism.

METHODS

Five datasets were obtained from the Gene Expression Omnibus database. Unsupervised consensus cluster analysis was used to determine new PMOP subtypes. To determine the central genes and the core modules related to PMOP, the weighted gene co-expression network analysis (WCGNA) was applied. Gene Ontology enrichment analysis was used to explore the biological processes underlying key genes. Logistic regression univariate analysis was used to screen for statistically significant variables. Two algorithms were used to select important PMOP-related genes. A logistic regression model was used to construct the PMOP-related gene profile. The receiver operating characteristic area under the curve, Harrell's concordance index, a calibration chart, and decision curve analysis were used to characterize PMOP-related genes. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the expression of the PMOP-related genes in the gene signature.

RESULTS

We identified three PMOP-related subtypes and four core modules. The muscle system process, muscle contraction, and actin filament-based movement were more active in the hub genes. We obtained five feature genes related to PMOP. Our analysis verified that the gene signature had good predictive power and applicability. The outcomes of the GSE56815 cohort were found to be consistent with the results of the earlier studies. qRT-PCR results showed that RAB2A and FYCO1 were amplified in clinical samples.

CONCLUSION

The PMOP-related gene signature we developed and verified can accurately predict the risk of PMOP in patients. These results can elucidate the molecular mechanism of RAB2A and FYCO1 underlying PMOP, and yield new and improved treatment strategies, ultimately helping PMOP monitoring.Cite this article: Bone Joint Res 2022;11(8):548-560.

Integrated Analysis of miRNAs and Gene Expression Profiles Reveals Potential Biomarkers for Osteoarthritis

Purpose: Currently, the early diagnosis and treatment of osteoarthritis (OA) remain a challenge. In the present study, we attempted to explore potential biomarkers for the diagnosis and treatment of OA.

Methods: The differentially expressed genes (DEGs) were identified based on three mRNA datasets of synovial tissues for OA patients and normal controls downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used for evaluating gene function related categories. Then, miRNA sequencing was performed for differentially expressed miRNAs’ identification. Finally, weighted gene co-expression network analysis (WGCNA) was performed for genes detected by the three mRNA datasets and a competing endogenous RNA (ceRNA) network with DEGs and differentially expressed microRNAs (miRNAs) was constructed for central genes identification. In addition, the relationship between central gene expression and immune infiltration was analyzed, and the candidate agents for OA were predicted based on the Connectivity Map database. Quantitative RT-PCR (qRT-PCR), Western blotting analysis, and immunofluorescent staining were performed to validate the expression levels of differentially expressed miRNAs and differentially expressed target genes in normal and OA tissues and chondrocytes. MiRNA–mRNA network was also validated in chondrocytes in vitro.

Results: A total of 259 DEGs and 26 differentially expressed miRNAs were identified, among which 94 miRNA–mRNA interactions were predicted. The brown module in WGCNA was most closely correlated with the clinical traits of OA. After overlapping the brown module genes with miRNA–mRNA pairs, 27 miRNA–mRNA pairs were obtained. A ceRNA network was constructed with 5505 lncRNA–miRNA–mRNA interactions. B-cell translocation gene 2(BTG2), Abelson-related gene (ABL2), and vascular endothelial growth factor A (VEGFA) were identified to be the central genes with good predictive performance, which were significantly correlated with immune cell infiltration in OA, reflected by declined activated dendritic cells (aDCs), and elevated contents of B cells, macrophages, neutrophils, and T helper cells. Anisomycin, MG-132, thapsigargin, and lycorine were predicted to be the potential candidate agents for OA intervention. In vitro, the expression levels of differentially expressed miRNAs and biomarkers identified in the present study were consistent with the results obtained in normal or OA knee cartilage tissues and chondrocytes. Furthermore, BTG2 was identified to be negatively regulated by miR-125a-5p.

Conclusion: BTG2, ABL2, and VEGFA can be regarded as potential predictive and treatment biomarkers for OA, which might guide the clinical therapy of OA.

Association of X-linked TLR-7 gene polymorphism with the risk of knee osteoarthritis: a case-control study

Knee osteoarthritis (OA) is the most prevalent type of OA, and Toll-like receptor 7 (TLR7) may lead to the pathogenesis of OA. Recently, X-linked TLR7 polymorphism has been confirmed to be associated with arthritis. However, there is a lack of studies on TLR7 gene polymorphism associated with knee OA susceptibility. The current study aimed to determine whether TLR7 gene polymorphism is associated with the risk of knee OA. Genotyping of two polymorphic sites (rs3853839 and rs179010) in the TLR7 gene was performed in 252 OA patients, and 265 healthy controls using the SNaPshot sequencing technique. Data were analyzed statistically by Chi-square tests and logistic regression. Rs3853839-C allele showed frequencies of 28% and 27% in the healthy control and female knee OA groups, respectively. The differences were not statistically significant (P > 0.05). The rs3853839-CG genotype frequency was significantly lower in the female knee OA group as compared to the healthy control group (OR 0.60; 95%CI 0.36–0.99; P = 0.044). In the male hemizygote population, the rs3853839-CC showed significantly lower frequencies in the male knee OA group as compared to the healthy control group (OR 0.35; 95%CI 0.17–0.71; P = 0.0025). Regarding rs179010, there were no differences in the genotype distribution and allele frequencies between OA patients and healthy subjects under any models (P > 0.05). Stratified analysis showed that the frequency of the rs3853839-CG genotypes was lower in high Kellgren-Lawrence grades (KLG) (OR 0.48; 95%CI 0.21–1.08; P = 0.066), and significantly lower in OA patients with effusion synovitis (OR 0.38; 95%CI 0.17–0.88; P = 0.013). TLR7 rs3853839 polymorphism may play a role in the susceptibility of knee OA in the Chinese Han Population and may be associated with OA severity and the risk of effusion synovitis in Knee OA.

Single-Cell Profiles of Age-Related Osteoarthritis Uncover Underlying Heterogeneity Associated With Disease Progression

Objective: Osteoarthritis (OA) is the most common chronic degenerative joint disease, which represents the leading cause of age-related disability. Here, this study aimed to depict the intercellular heterogeneity of OA synovial tissues.

Methods: Single-cell RNA sequencing (scRNA-seq) data were preprocessed and quality controlled by the Seurat package. Cell cluster was presented and cell types were annotated based on the mRNA expression of corresponding marker genes by the SingleR package. Cell-cell communication was assessed among different cell types. After integrating the GSE55235 and GSE55457 datasets, differentially expressed genes were identified between OA and normal synovial tissues. Then, differentially expressed marker genes were overlapped and their biological functions were analyzed.

Results: Totally, five immune cell subpopulations were annotated in OA synovial tissues including macrophages, dendritic cells, T cells, monocytes and B cells. Pseudo-time analysis revealed the underlying evolution process in the inflammatory microenvironment of OA synovial tissue. There was close crosstalk between five cell types according to the ligand-receptor network. The genetic heterogeneity was investigated between OA and normal synovial tissues. Furthermore, functional annotation analysis showed the intercellular heterogeneity across immune cells in OA synovial tissues.

Conclusion: This study offered insights into the heterogeneity of OA, which provided in-depth understanding of the transcriptomic diversities within synovial tissue. This transcriptional heterogeneity may improve our understanding on OA pathogenesis and provide potential molecular therapeutic targets for OA.

Silencing of TLR7 protects against lipopolysaccharide-induced chondrocyte apoptosis and injury by blocking the p21-mediated JAK2/STAT3 pathway

Involvement of toll-like receptor 7 (TLR7) in the immune response has been reported in diverse inflammatory diseases. However, the role of TLR7 in the pathogenesis of osteoarthritis (OA) is poorly understood. In this study, we sought to investigate the contribution of TLR7 in regulating chondrocyte apoptosis, inflammation, and degradation of the extracellula matrix (ECM), and its underlying mechanisms. We found that TLR7 expression was increased in cartilage tissues of OA patients and in lipopolysaccharide (LPS)-induced chondrocytes. Silencing of TLR7 alleviated LPS-induced chondrocyte apoptosis, inflammation, and ECM degradation. Mechanistically, TLR7 silencing inhibited the JAK2/STAT3 signaling pathway by inducing p21 expression. Moreover, p21 knockdown and colivein (an activator of JAK2/STAT3 signaling) partially rescued the suppressive role of TLR7 silencing on chondrocyte apoptosis, the inflammatory response, and ECM underproduction. Taken together, our data revealed that knockdown of TLR7 attenuated chondrocyte apoptosis and injury by blocking the p21-mediated JAK2/STAT3 pathway, suggesting that TLR7 may be a therapeutic target in OA.

Recent Advances in Pharmacological Intervention of Osteoarthritis: A Biological Aspect

Osteoarthritis (OA) is a degenerative joint disease in the musculoskeletal system with a relatively high incidence and disability rate in the elderly. It is characterized by the degradation of articular cartilage, inflammation of the synovial membrane, and abnormal structure in the periarticular and subchondral bones. Although progress has been made in uncovering the molecular mechanism, the etiology of OA is still complicated and unclear. Nevertheless, there is no treatment method that can effectively prevent or reverse the deterioration of cartilage and bone structure. In recent years, in the field of pharmacology, research focus has shifted to disease prevention and early treatment rather than disease modification in OA. Biologic agents become more and more attractive as their direct or indirect intervention effects on the initiation or development of OA. In this review, we will discuss a wide spectrum of biologic agents ranging from DNA, noncoding RNA, exosome, platelet-rich plasma (PRP), to protein. We searched for key words such as OA, DNA, gene, RNA, exosome, PRP, protein, and so on. From the pharmacological aspect, stem cell therapy is a very special technique, which is not included in this review. The literatures ranging from January 2016 to August 2021 were included and summarized. In this review, we aim to help readers have a complete and precise understanding of the current pharmacological research progress in the intervention of OA from the biological aspect and provide an indication for the future translational studies.

Maternal caffeine intake and DNA methylation in newborn cord blood

BACKGROUND

Epigenetic mechanisms may underlie associations between maternal caffeine consumption and adverse childhood metabolic outcomes. However, limited studies have examined neonate DNA methylation (DNAm) patterns in the context of preconception or prenatal exposure to caffeine metabolites.

OBJECTIVES

We examined preconception and pregnancy caffeine exposure with DNAm alterations in neonate cord blood (n = 378).

METHODS

In a secondary analysis of the Effects of Aspirin in Gestation and Reproduction Trial (EAGeR), we measured maternal caffeine, paraxanthine, and theobromine concentrations from stored serum collected preconception (on average 2 months before pregnancy) and at 8 weeks of gestation. In parallel, self-reported caffeinated beverage intake was captured via administration of questionnaires and daily diaries. We profiled DNAm from the cord blood buffy coat of singletons using the MethylationEPIC BeadChip. We assessed associations of maternal caffeine exposure and methylation β values using multivariable robust linear regression. A false discovery rate (FDR) correction was applied using the Benjamini-Hochberg method.

RESULTS

In preconception, the majority of women reported consuming 1 or fewer servings/day of caffeine on average, and caffeine and paraxanthine metabolite levels were 88 and 36 µmol/L, respectively. Preconception serum caffeine metabolites were not associated with individual cytosine-guanine (CpG) sites (FDR >5%), though pregnancy theobromine was associated with DNAm at cg09460369 near RAB2A (β = 0.028; SE = 0.005; FDR P = 0.012). Preconception self-reported caffeinated beverage intake compared to no intake was associated with DNAm at cg09002832 near GLIS3 (β = -0.013; SE = 0.002; FDR P = 0.036). No associations with self-reported intake during pregnancy were found.

CONCLUSIONS

Few effects of maternal caffeine exposure on neonate methylation differences in leukocytes were identified in this population with relatively low caffeine consumption.

Differentially expressed genes in the femur cartilage transcriptome clarify the understanding of femoral head separation in chickens

Locomotor problems are among one of the main concerns in the current poultry industry, causing major economic losses and affecting animal welfare. The most common bone anomalies in the femur are dyschondroplasia, femoral head separation (FHS), and bacterial chondronecrosis with osteomyelitis (BCO), also known as femoral head necrosis (FHN). The present study aimed to identify differentially expressed (DE) genes in the articular cartilage (AC) of normal and FHS-affected broilers by RNA-Seq analysis. In the transcriptome analysis, 12,169 genes were expressed in the femur AC. Of those, 107 genes were DE (FDR < 0.05) between normal and affected chickens, of which 9 were downregulated and 98 were upregulated in the affected broilers. In the gene-set enrichment analysis using the DE genes, 79 biological processes (BP) were identified and were grouped into 12 superclusters. The main BP found were involved in the response to biotic stimulus, gas transport, cellular activation, carbohydrate-derived catabolism, multi-organism regulation, immune system, muscle contraction, multi-organism process, cytolysis, leukocytes and cell adhesion. In this study, the first transcriptome analysis of the broilers femur articular cartilage was performed, and a set of candidate genes (AvBD1, AvBD2, ANK1, EPX, ADA, RHAG) that could trigger changes in the broiler´s femoral growth plate was identified. Moreover, these results could be helpful to better understand FHN in chickens and possibly in humans.

Males and Females Have Distinct Molecular Events in the Articular Cartilage during Knee Osteoarthritis

Osteoarthritis (OA) is a major public health challenge that imposes a remarkable burden on the affected individuals and the healthcare system. Based on the clinical observation, males and females have different prevalence rates and severity levels of OA. Thus, sex-based differences may play essential roles in OA’s prognosis and treatment outcomes. To date, the comprehensive understanding of the relationship between sex and OA is still largely lacking. In the current study, we analyzed a published transcriptome dataset of knee articular cartilage (GSE114007) from 18 healthy (five females, 13 males) and 20 OA (11 females, nine males) donors to provide a slight insight into this important but complex issue. First, comparing female healthy cartilage samples with those of males revealed 36 differential expression genes (DEGs), indicating the fundamental sex-related differences at the molecular level. Meanwhile, 923 DEGs were distinguished between OA and healthy female cartilage, which can be enriched to 15 Reactome pathways. On the other hand, when comparing OA and healthy male cartilage, there are only 419 DEGs were identified, and only six pathways were enriched against the Reactome database. The different signaling response to OA in the male and female cartilage was further enforced by recognizing 50 genes with significantly different OA-responsive expression fold changes in males and females. Particularly, 14 Reactome pathways, such as “Extracellular matrix organization”, “Collagen biosynthesis and modifying enzymes”, “Dissolution of fibrin clot”, and “Platelet Aggregation (Plug formation)”, can be noted from these 50 sex-dependent OA-responsive genes. Overall, the current study explores the Sex as a Biological Variable (SABV) at the transcriptomic level in the knee articular cartilage in both healthy status and OA event, which could help predict the differential OA prognosis and treatment outcome of males and female patients.