Increased DNA Methylation and Reduced Expression of Transcription Factors in Human Osteoarthritis Cartilage

WTAP-mediated m6A modification of FRZB triggers the inflammatory response via the Wnt signaling pathway in osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis. However, the exact pathogenesis remains unclear. Emerging evidence shows that N6-methyladenosine (m6A) modification may have an important role in OA pathogenesis. This study aimed to investigate the role of m6A writers and the underlying mechanisms in osteoarthritic cartilage. Among m6A methyltransferases, Wilms tumor 1-associated protein (WTAP) expression most significantly differed in clinical osteoarthritic cartilage. WTAP regulated extracellular matrix (ECM) degradation, inflammation and antioxidation in human chondrocytes. Mechanistically, the m6A modification and relative downstream targets in osteoarthritic cartilage were assessed by methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing, which indicated that the expression of frizzled-related protein (FRZB), a secreted Wnt antagonist, was abnormally decreased and accompanied by high m6A modification in osteoarthritic cartilage. In vitro dysregulated WTAP had positive effects on β-catenin expression by targeting FRZB, which finally contributed to the cartilage injury phenotype in chondrocytes. Intra-articular injection of adeno-associated virus-WTAP alleviated OA progression in a mouse model, while this protective effect could be reversed by the application of a Wnt/β-catenin activator. In summary, this study revealed that WTAP-dependent RNA m6A modification contributed to Wnt/β-catenin pathway activation and OA progression through post-transcriptional regulation of FRZB mRNA, thus providing a potentially effective therapeutic strategy for OA treatment.

Integrated Analysis of DNA Methylation and Gene Expression Profiles in a Rat Model of Osteoarthritis

Osteoarthritis (OA) is common and affected by several factors, such as age, weight, sex, and genetics. The pathogenesis of OA remains unclear. Therefore, using a rat model of monosodium iodoacetate (MIA)-induced OA, we examined genomic-wide DNA methylation using methyl-seq and characterized the transcriptome using RNA-seq in the articular cartilage tissue from a negative control (NC) and MIA-induced rats. We identified 170 genes (100 hypomethylated and upregulated genes and 70 hypermethylated and downregulated genes) regulated by DNA methylation in OA. DNA methylation-regulated genes were enriched in functions related to focal adhesion, extracellular matrix (ECM)-receptor interaction and the PI3K-Akt and Hippo signaling pathways. Functions related to extracellular matrix organization, extracellular matrix proteoglycans, and collagen formation were involved in OA. A molecular and protein-protein network was constructed using methylated expression-correlated genes. Erk1/2 was a downstream target of OA-induced changes in DNA methylation and RNA expression. We found that the integrin subunit alpha 2 (ITGA2) gene is important in focal adhesion, alpha6-beta4 integrin signaling, and the inflammatory response pathway in OA. Overall, gene expression changes because DNA methylation influences OA pathogenesis. ITGA2, whose gene expression changes are regulated by DNA methylation during OA onset, is a candidate gene. Our findings provide insights into the epigenetic targets of OA processes in rats.

ETS2 overexpression ameliorates cartilage injury in osteoarthritis by the ETS2/miR-155/STAT1/DNMT1 feedback loop pathway

Osteoarthritis (OA) is the most common irreversible chronic joint dysfunction disease, which is pathologically characterized by disturbance of articular cartilage homeostasis leading to subsequent inflammatory response and cartilage extracellular matrix (ECM) degradation. Increasing evidence has demonstrated the dysregulation of transcription factors play crucial roles in the occurrence and development of osteoarthritis (OA), but the potential functions and mechanism of most transcription factors in OA has not been completely illuminated. In this study, we identified that transcription factor V-ets erythroblastosis virus E26 oncogene homolog 2 (ETS2) was significantly down-regulated in OA cartilage and IL-1β-induced OA chondrocytes. Functional experiments in vitro demonstrated that the overexpressed ETS2 strikingly enhanced proliferation, outstandingly suppressed apoptosis, and dramatically reduced inflammation and ECM degradation in IL-1β-induced OA chondrocytes, whereas the knockdown of ETS2 led to the opposite effects. Further in vivo studies have shown that up-regulated ETS2 dramatically ameliorates cartilage injury in DMM-induced OA mice. Mechanical studies have disclosed that DNMT1-mediated downregulation of ETS2 dramatically promotes STAT1 by inhibiting miR-155 transcription, and increased STAT1 initiates a feedback loop that may enhance DNMT1-mediated hypermethylation of ETS2 to inhibit ETS2 expression, thus forming a DNMT1/ETS2/miR-155/STAT1 feedback loop that inhibits MAPK signaling pathways and aggravates OA cartilage injury. In all, our results revealed that overexpression of ETS2 markedly ameliorated OA cartilage injury through the ETS2/miR-155/STAT1/DNMT1 feedback loop, providing a new perspective on the pathogenesis and therapeutic strategies for OA.

Identification of Diagnostic Markers in Synovial Tissue of Osteoarthritis by Weighted Gene Coexpression Network

Osteoarthritis (OA) is a serious threat to human health. However, the etiology and pathogenesis of the disease are not fully understood. Most researchers believe that the degeneration and imbalance of articular cartilage, extracellular matrix, and subchondral bone are the fundamental causes of osteoarthritis. However, recent studies have shown that synovial lesions may precede cartilage, which may be an important precipitating factor in the early stage of OA and the whole course of the disease. This study aimed to conduct an analysis based on sequence data from the Gene Expression Omnibus (GEO) database to investigate the presence of effective biomarkers in the synovial tissue of osteoarthritis for the diagnosis and control of OA progression. In this study, the differentially expressed OA-related genes (DE-OARGs) in osteoarthritis synovial tissues were extracted in the GSE55235 and GSE55457 datasets using the Weighted Gene Co-expression Network Analysis (WGCNA) and limma. Least-Absolute Shrinkage and Selection Operator (LASSO) algorithm was used to select the diagnostic genes based on the DE-OARGs by glmnet package. 7 genes were selected as diagnostic genes including SAT1, RLF, MAFF, SIK1, RORA, ZNF529, and EBF2. Subsequently, the diagnostic model was constructed and the results of the Area Under the Curve (AUC) demonstrated that the diagnostic model had high diagnostic performance for OA. Additionally, among the 22 immune cells of the Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) and the 24 immune cells of the single sample Gene Set Enrichment Analysis (ssGSEA), 3 immune cells and 5 immune cells were different between the OA and normal samples, respectively. The expression trends of the 7 diagnostic genes were consistent in the GEO datasets and the results of the real-time reverse transcription PCR (qRT-PCR). The results of this study demonstrate that these diagnostic markers have important significance in the diagnosis and treatment of OA, and will provide further evidence for the clinical and functional studies of OA.

Identification of four-gene signature to diagnose osteoarthritis through bioinformatics and machine learning methods

BACKGROUND

Although osteoarthritis (OA) is one of the most prevalent joint disorders, effective biomarkers to diagnose OA are still unavailable. This study aimed to acquire some key synovial biomarkers (hub genes) and analyze their correlation with immune infiltration in OA.

METHODS

Gene expression profiles and clinical characteristics of OA and healthy synovial samples were retrieved from the Gene Expression Omnibus (GEO) database. Hub genes for OA were mined based on a combination of weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) algorithms. A diagnostic nomogram model for OA prediction was developed based on the hub genes. Receiver operating characteristic curves (ROC) were performed to confirm the abnormal expression of hub genes in the experimemtal and validation datasets. qRT-PCR using patients' samples were conducted as well. In addition, the infiltration level of 28 immune cells in the expression profile and their relationship with hub genes were analyzed using single-sample GSEA (ssGSEA).

RESULTS

4 hub genes (ZBTB16, TNFSF11, SCRG1 and KDELR3) were obtained by WGCNA, lasso, SVM-RFE, RF algorithms as potential biomarkers for OA. The immune infiltration analyses revealed that hub genes were most correlated with regulatory T cell and natural killer cell.

CONCLUSION

A machine learning model to diagnose OA based on ZBTB16, TNFSF11, SCRG1 and KDELR3 using synovial tissue was constructed, providing theoretical foundation and guideline for diagnostic and treatment targets in OA.

Whole-genome DNA methylation and DNA methylation-based biomarkers in lung squamous cell carcinoma

Exploring early detection methods through comprehensive evaluation of DNA methylation for lung squamous cell carcinoma (LUSC) patients is of great significance. By using different machine learning algorithms for feature selection and model construction based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, five methylation biomarkers in LUSC (along with mapped genes) were identified including cg14823851 (TBX4), cg02772121 (TRIM15), cg10424681 (C6orf201), cg12910906 (ARHGEF4), and cg20181079 (OR4D11), achieving extremely high sensitivity and specificity in distinguishing LUSC from normal samples in independent cohorts. Pyrosequencing assay verified DNA methylation levels, meanwhile qRT-PCR and immunohistochemistry results presented their accordant methylation-related gene expression statuses in paired LUSC and normal lung tissues. The five methylation-based biomarkers proposed in this study have great potential for the diagnosis of LUSC and could guide studies in methylation-regulated tumor development and progression.

The Role of Genetics and Epigenetic Regulation in the Pathogenesis of Osteoarthritis

Osteoarthritis (OA) is progressive disease characterised by cartilage degradation, subchondral bone remodelling and inflammation of the synovium. The disease is associated with obesity, mechanical load and age. However, multiple pro-inflammatory immune mediators regulate the expression of metalloproteinases, which take part in cartilage degradation. Furthermore, genetic factors also contribute to OA susceptibility. Recent studies have highlighted that epigenetic mechanisms may regulate the expression of OA-associated genes. This review aims to present the mechanisms of OA pathogenesis and summarise current evidence regarding the role of genetics and epigenetics in this process.

ZBTB16 eases lipopolysaccharide‑elicited inflammation, apoptosis and degradation of extracellular matrix in chondrocytes during osteoarthritis by suppressing GRK2 transcription

Osteoarthritis (OA) is a chronic degenerative disease of the bone that is a major contributor of disability in the elderly population. Zinc finger and BTB domain-containing 16 (ZBTB16) is a transcription factor that has been previously revealed to be impaired in human OA tissues. The present study was designed to elaborate the potential impact of ZBTB16 on OA and to possibly assess any latent regulatory mechanism. ZBTB16 expression in human OA tissues was examined using the Gene Expression Series (GSE) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE169077) whereas ZBTB16 expression in chondrocytes was examined using reverse transcription-quantitative PCR (RT-qPCR) and western blotting. Cell viability was examined using a Cell Counting Kit-8 assay. A TUNEL assay and western blotting were used to assess cell apoptosis and apoptosis-related markers, including Bcl-2, Bax and cleaved caspase-3. The levels and expression of inflammatory factors, including TNF-α, IL-1β and IL-6, were determined by ELISA and western blotting. RT-qPCR and western blotting were also used to analyze the expression levels of extracellular matrix (ECM)-degrading enzymes, including MMP-13, a disintegrin-like and metalloproteinase with thrombospondin type-1 motifs-5, aggrecan and collagen type II α1. After the potential binding of ZBTB16 with the G protein coupled receptor kinase type 2 (GRK2) promoter was predicted using the Cistrome DB database, GRK2 expression was confirmed by RT-qPCR and western blotting. Chromatin immunoprecipitation and luciferase reporter assays were then used to determine the potential interaction between ZBTB16 and the GRK2 promoter. Following GRK2 overexpression in ZBTB16-overexpressing chondrocytes by co-transfection of GRK2 and ZBTB16 overexpression plasmids, the aforementioned functional experiments were performed again. ZBTB16 expression was found to be reduced in human OA tissues compared with in normal cartilage tissues and lipopolysaccharide (LPS)-stimulated chondrocytes. ZBTB16 overexpression increased cell viability whilst decreasing apoptosis, inflammation and ECM degradation by LPS-treated chondrocytes. In addition, GRK2 expression was found to be increased in LPS-stimulated chondrocytes. ZBTB16 successfully bound to the GRK2 promoter, which negatively modulated GRK2 expression. GRK2 upregulation reversed the effects of ZBTB16 overexpression on the viability, apoptosis, inflammation and ECM degradation by LPS-challenged chondrocytes. In conclusion, these data suggest that ZBTB16 may inhibit the development of OA through the transcriptional inactivation of GRK2.

Local diversity of drought resistance and resilience in Populus tomentosa correlates with the variation of DNA methylation

Little information is known about DNA methylation variation in shaping environment-specific drought resistance and resilience for tree adaptation. In this study, we leveraged RNA sequencing and whole-genome bisulfite sequencing data to dissect the distinction of epigenetic regulation under drought stress and rewater condition of Populus tomentosa accessions from three geographical regions. We demonstrated low resistance and high resilience for accessions from South. Non-CG methylation levels in promoter regions of Southern accessions were lower than accessions from higher latitudes and negatively regulated gene expression. CHH context methylation was more sensitive to drought stress, and the geographical-specific differentially methylated regions were scarcely changed by environmental fluctuation. We identified 60 conserved hub genes within the co-expression networks that correlate with photosynthetic and stomatal morphological traits. Epigenome-wide association studies and genome-wide association studies of these 60 hub genes revealed the interdependency between genetic and epigenetic variation in GATA9 and LECRK-VIII.2, which was associated with stomatal morphology and chlorophyll content. The natural epigenetic variation in GATA9 was also faithfully transmitted to progenies in two family-based F1 populations. This study indicates a functional relationship of DNA methylation diversity with drought resistance and resilience which offers new insights into plants' local adaptation to a stressful environment.

The epigenetic players and the chromatin marks involved in the articular cartilage during osteoarthritis

Epigenetics defines the modifications of the genome that do not involve a change in the nucleotide sequence of DNA. These modifications constitute a mechanism of gene regulation poorly explored in the context of cartilage physiology. They are now intensively studied by the scientific community working on articular cartilage and its related pathology such as osteoarthritis. Indeed, epigenetic regulations can control the expression of crucial gene in the chondrocytes, the only resident cells of cartilage. Some epigenetic changes are considered as a possible cause of the abnormal gene expression and the subsequent alteration of the chondrocyte phenotype (hypertrophy, proliferation, senescence…) as observed in osteoarthritic cartilage. Osteoarthritis is a joint pathology, which results in impaired extracellular matrix homeostasis and leads ultimately to the progressive destruction of cartilage. To date, there is no pharmacological treatment and the exact causes have yet to be defined. Given that the epigenetic modifying enzymes can be controlled by pharmacological inhibitors, it is thus crucial to describe the epigenetic marks that enable the normal expression of extracellular matrix encoding genes, and those associated with the abnormal gene expression such as degradative enzyme or inflammatory cytokines encoding genes. In this review, only the DNA methylation and histone modifications will be detailed with regard to normal and osteoarthritic cartilage. Although frequently referred as epigenetic mechanisms, the regulatory mechanisms involving microRNAs will not be discussed. Altogether, this review will show how this nascent field influences our understanding of the pathogenesis of OA in terms of diagnosis and how controlling the epigenetic marks can help defining epigenetic therapies.

Unbiased transcriptome mapping and modeling identify candidate genes and compounds of osteoarthritis

Osteoarthritis (OA) is a chronic degenerative joint disease characterized by progressive cartilage loss, subchondral bone remodeling, and synovial inflammation. Given that the current therapies for advanced OA patients are limited, the understanding of mechanisms and novel therapies are urgently needed. In this study, we employed the weighted gene co-expression network (WGCNA) method and the connectivity map (CMap) database to identify the candidate target genes and potential compounds. Four groups of co-expressing genes were identified as the OA-related modules. The biological annotations of these modules indicated some critical hallmarks of OA and aging, such as mitochondrial dysfunctions and abnormal energy metabolism, and the signaling pathways, such as MAPK, TNF, and PI3K/Akt signaling pathways. Some genes, such as RELA and GADD45B, were predicted to extensively involve these critical pathways, indicating their potential functions in OA mechanisms. Moreover, we constructed the co-expressing networks of modules and identified the hub genes based on network topology. GADD45B, MAFF, and MYC were identified and validated as the hub genes. Finally, anisomycin and MG-262 were predicted to target these OA-related modules, which may be the potential drugs for OA therapy. In conclusion, this study identified the significant modules, signaling pathways, and hub genes relevant to OA and highlighted the potential clinical value of anisomycin and MG-262 as novel therapies in OA management.

Multi-omics molecular biomarkers and database of osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis in the adult population and is a leading cause of disability. OA-related genetic loci may play an important role in clinical diagnosis and disease progression. With the rapid development of diverse technologies and omics methods, many OA-related public data sets have been accumulated. Here, we retrieved a diverse set of omics experimental results from 159 publications, including genome-wide association study, differentially expressed genes and differential methylation regions, and 2405 classified OA-related gene markers. Meanwhile, based on recent single-cell RNA-seq data from different joints, 5459 cell-type gene markers of joints were collected. The information has been integrated into an online database named OAomics and molecular biomarkers (OAOB). The database (http://ibi.zju.edu.cn/oaobdb/) provides a web server for OA marker genes, omics features and so on. To our knowledge, this is the first database of molecular biomarkers for OA.

The crucial mechanism and therapeutic implication of RNA methylation in bone pathophysiology

Methylation is the most common posttranscriptional modification in cellular RNAs, which has been reported to modulate the alteration of RNA structure for initiating relevant functions such as nuclear translocation and RNA degradation. Recent studies found that RNA methylation especially N6-methyladenosine (m⁶A) regulates the dynamic balance of bone matrix and forms a complicated network in bone metabolism. The modulation disorder of RNA methylation contributes to several pathological bone diseases including osteoporosis (OP), osteoarthritis (OA), rheumatoid arthritis (RA), and so on. In the review, we will discuss advanced technologies for detecting RNA methylation, summarize RNA methylation-related biological impacts on regulating bone homeostasis and pathological bone diseases. In addition, we focus on the promising roles of RNA methylation in early diagnosis and therapeutic implications for bone-related diseases. Then, we aim to establish a theoretical basis for further investigation in this meaningful field.

DNA methylation and noncoding RNA in OA: Recent findings and methodological advances

Introduction:

Osteoarthritis (OA) is a chronic musculoskeletal disease characterized by progressive loss of joint function. Historically, it has been characterized as a disease caused by mechanical trauma, so-called ‘wear and tear’. Over the past two decades, it has come to be understood as a complex systemic disorder involving gene-environmental interactions. Epigenetic changes have been increasingly implicated. Recent improvements in microarray and next-generation sequencing (NGS) technologies have allowed for ever more complex evaluations of epigenetic aberrations associated with the development and progression of OA.

Methods:

A systematic review was conducted in the Pubmed database. We curated studies that presented the results of DNA methylation and noncoding RNA research in human OA and OA animal models since 1985.

Results:

Herein, we discuss recent findings and methodological advancements in OA epigenetics, including a discussion of DNA methylation, including microarray and NGS studies, and noncoding RNAs. Beyond cartilage, we also highlight studies in subchondral bone and peripheral blood mononuclear cells, which highlight widespread and potentially clinically important alterations in epigenetic patterns seen in OA patients. Finally, we discuss epigenetic editing approaches in the context of OA.

Conclusions:

Although a substantial body of literature has already been published in OA, much is still unknown. Future OA epigenetics studies will no doubt continue to broaden our understanding of underlying pathophysiology and perhaps offer novel diagnostics and/or treatments for human OA.

Single-cell transcriptomics of popliteal lymphatic vessels and peripheral veins reveals altered lymphatic muscle and immune cell populations in the TNF-Tg arthritis model

Background

Lymphatic dysfunction exists in tumor necrosis factor transgenic (TNF-Tg) mice and rheumatoid arthritis (RA) patients. While joint-draining TNF-Tg popliteal lymphatic vessels (PLVs) have deficits in contractility during end-stage arthritis, the nature of lymphatic muscle cells (LMCs) and their TNF-altered transcriptome remain unknown. Thus, we performed single-cell RNA-sequencing (scRNAseq) on TNF-Tg LMCs in PLVs efferent to inflamed joints versus wild-type (WT) controls.

Methods

Single-cell suspensions of PLVs were sorted for smooth muscle cells (SMCs), which was validated by Cspg4-Cre;tdTomato reporter gene expression. Single-cell RNA-seq was performed on a 10x Genomics platform and analyzed using the Seurat R package. Uniform Manifold Approximation and Projections (UMAPs) and Ingenuity Pathway Analysis software were used to assess cell clusters and functional genomics in WT vs. TNF-Tg populations.

Results

Fluorescent imaging of Cspg4-Cre;tdTomato vessels demonstrated dim PLVs and strong reporter gene expression in the adjacent superficial saphenous vein, which was corroborated by flow cytometry of LMCs and vascular smooth muscle cells (VSMCs) from these vessels. Due to their unique morphology, these populations could also be readily detected by scatter analysis of cells from non-fluorescent mice. Bioinformatics analysis of flow sorted WT and TNF-Tg cells identified 20 unique cell clusters that together were 22.4% LMCs, 15.0% VSMCs, and 62.6% non-muscle cells of 8879 total cells. LMCs and M2-macrophages were decreased, while inflammatory monocytes were increased in TNF-Tg lower limb vasculature. SMC populations were defined by Cald1, Tpm1, and Pdgfrb expression and were enriched in myofibroblast-like gene expression. TNF-Tg LMCs exhibited enhanced functional genomics associated with cell death, phagocyte recruitment, and joint inflammation. Among the most prominent TNF-induced genes in SMCs were Mmp3, Cxcl12, and Ccl19, and the most downregulated genes were Zbtb16, Galnt15, and Apod.

Conclusions

Single-cell RNA-seq can be used to investigate functional genomics of lower limb vasculature in mice. Our findings confirm the inflammatory transcriptome of TNF-Tg vessels and altered gene expression in SMC populations. This study further supports a potential role of mesenchymal stromal cells in inflammatory-erosive arthritis pathogenesis, and warrants future studies to define the effects of this TNF-altered transcriptome on PLV function and joint homeostasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-022-02730-z.

NCOR2 is a novel candidate gene for migraine-epilepsy phenotype

HYPOTHESIS

To identify genetic factors predisposing to migraine-epilepsy phenotype utilizing a multi-generational family with known linkage to chr12q24.2-q24.3.

METHODS

We used single nucleotide polymorphism (SNP) genotyping and next-generation sequencing technologies to perform linkage, haplotype, and variant analyses in an extended Finnish migraine-epilepsy family (n = 120). In addition, we used a large genome-wide association study (GWAS) dataset of migraine and two biobank studies, UK Biobank and FinnGen, to test whether variants within the susceptibility region associate with migraine or epilepsy related phenotypes in a population setting.

RESULTS

The family showed the highest evidence of linkage (LOD 3.42) between rs7966411 and epilepsy. The haplotype shared among 12 out of 13 epilepsy patients in the family covers almost the entire NCOR2 and co-localizes with one of the risk loci of the recent GWAS on migraine. The haplotype harbors nine low-frequency variants with potential regulatory functions. Three of them, in addition to two common variants, show nominal associations with neurological disorders in either UK Biobank or FinnGen.

CONCLUSION

We provide several independent lines of evidence supporting association between migraine-epilepsy phenotype and NCOR2. Our study suggests that NCOR2 may have a role in both migraine and epilepsy and thus would provide evidence for shared pathophysiology underlying these two diseases.

General Caspase Inhibition in Primary Chondrogenic Cultures Impacts Their Transcription Profile Including Osteoarthritis-Related Factors

OBJECTIVE

The knowledge about functions of caspases, usually associated with cell death and inflammation, keeps expanding also regarding cartilage. Active caspases are present in the growth plate, and caspase inhibition in limb-derived chondroblasts altered the expression of osteogenesis-related genes. Caspase inhibitors were reported to reduce the severity of cartilage lesions in osteoarthritis (OA), and caspase-3 might represent a promising biomarker for OA prognosis. The objective of this investigation was to decipher the transcriptomic regulation of caspase inhibition in chondrogenic cells.

DESIGN

Limb-derived chondroblasts were cultured in the presence of 2 different inhibitors: Z-VAD-FMK (FMK) and Q-VD-OPH (OPH). A whole transcriptome RNA sequencing was performed as the key analysis.

RESULTS

The analysis revealed a statistically significant increase in the expression of 252 genes in the FMK samples and 163 genes in the OPH samples compared with controls. Conversely, there was a significant decrease in the expression of 290 genes in the FMK group and 188 in the OPH group. Among the top up- and downregulated genes (more than 10 times changed), almost half of them were associated with OA. Both inhibitors displayed the highest upregulation of the inflammatory chemokine Ccl5, the most downregulated gene was the one for mannose receptors Mrc1.

CONCLUSIONS

The obtained datasets pointed to a significant impact of caspase inhibition on the expression of several chondro-/osteogenesis-related markers in an in vitro model of endochondral ossification. Notably, the list of these genes included some encoding for factors associated with cartilage/bone pathologies such as OA.

DNA methylation at a nutritionally sensitive region of the PAX8 gene is associated with thyroid volume and function in Gambian children

PAX8 is a key thyroid transcription factor implicated in thyroid gland differentiation and function, and PAX8 gene methylation is reported to be sensitive to the periconceptional environment. Using a novel recall-by-epigenotype study in Gambian children, we found that PAX8 hypomethylation at age 2 years is associated with a 21% increase in thyroid volume and an increase in free thyroxine (T4) at 5 to 8 years, the latter equivalent to 8.4% of the normal range. Free T4 was associated with a decrease in DXA-derived body fat and bone mineral density. Furthermore, offspring PAX8 methylation was associated with periconceptional maternal nutrition, and methylation variability was influenced by genotype, suggesting that sensitivity to environmental exposures may be under partial genetic control. Together, our results demonstrate a possible link between early environment, PAX8 gene methylation and thyroid gland development and function, with potential implications for early embryonic programming of thyroid-related health and disease.

An evolutionary perspective of DNA methylation patterns in skeletal tissues using a baboon model of osteoarthritis

Epigenetic factors, such as DNA methylation, play an influential role in the development of the degenerative joint disease osteoarthritis (OA). These molecular mechanisms have been heavily studied in humans, and although OA affects several other animals in addition to humans, few efforts have taken an evolutionary perspective. This study explores the evolution of OA epigenetics by assessing the relationship between DNA methylation variation and knee OA development in baboons (Papio spp.) and by comparing these findings to human OA epigenetic associations. Genome-wide DNA methylation patterns were identified in bone and cartilage of the right distal femora from 56 pedigreed, adult baboons (28 with and 28 without knee OA) using the Illumina Infinium MethylationEPIC BeadChip. Several significantly differentially methylated positions (DMPs) and regions were found between tissue types. Substantial OA-related differential methylation was also identified in cartilage, but not in bone, suggesting that cartilage epigenetics may be more influential in OA than bone epigenetics. Additionally, some genes containing OA-related DMPs overlap with and display methylation patterns similar to those previously identified in human OA, revealing a mixture of evolutionarily conserved and divergent OA-related methylation patterns in primates. Overall, these findings reinforce the current etiological perspectives of OA and enhance our evolutionary understanding of epigenetic mechanisms associated with OA. This study further establishes baboons as a valuable nonhuman primate model of OA, and continued investigations in baboons will help to disentangle the molecular mechanisms contributing to OA and their evolutionary histories.

Mohawk is a transcription factor that promotes meniscus cell phenotype and tissue repair and reduces osteoarthritis severity

Meniscus tears are common knee injuries and a major osteoarthritis (OA) risk factor. Knowledge gaps that limit the development of therapies for meniscus injury and degeneration concern transcription factors that control the meniscus cell phenotype. Analysis of RNA sequencing data from 37 human tissues in the Genotype-Tissue Expression database and RNA sequencing data from meniscus and articular cartilage showed that transcription factor Mohawk (MKX) is highly enriched in meniscus. In human meniscus cells, MKX regulates the expression of meniscus marker genes, OA-related genes, and other transcription factors, including Scleraxis (SCX), SRY Box 5 (SOX5), and Runt domain-related transcription factor 2 (RUNX2). In mesenchymal stem cells (MSCs), the combination of adenoviral MKX (Ad-MKX) and transforming growth factor-β3 (TGF-β3) induced a meniscus cell phenotype. When Ad-MKX-transduced MSCs were seeded on TGF-β3-conjugated decellularized meniscus scaffold (DMS) and inserted into experimental tears in meniscus explants, they increased glycosaminoglycan content, extracellular matrix interconnectivity, cell infiltration into the DMS, and improved biomechanical properties. Ad-MKX injection into mouse knee joints with experimental OA induced by surgical destabilization of the meniscus suppressed meniscus and cartilage damage, reducing OA severity. Ad-MKX injection into human OA meniscus tissue explants corrected pathogenic gene expression. These results identify MKX as a previously unidentified key transcription factor that regulates the meniscus cell phenotype. The combination of Ad-MKX with TGF-β3 is effective for differentiation of MSCs to a meniscus cell phenotype and useful for meniscus repair. MKX is a promising therapeutic target for meniscus tissue engineering, repair, and prevention of OA.

Effect of DNA methylation on gene transcription is associated with the distribution of methylation sites across the genome in osteoarthritis

Genetics and epigenetics are important subjects in the field of osteoarthritis (OA) research. DNA methylation may affect gene transcription, but the specific mechanisms have remained to be fully elucidated. In the present study, the ChAMP methylation analysis package was used to identify differentially methylated genes (DMGs) from the dataset GSE63695 from the Gene Expression Omnibus (GEO) database. The distribution of differentially methylated sites (DMS) and the total array sites across the genome were analyzed by enrichment analysis. Subsequently, two mRNA expression profiling datasets, GSE114007 and GSE113825, were obtained from the GEO database and common differentially expressed genes (DEGs) were identified using the Limma package. Key genes were screened by analyzing the distribution of DMS across the genome consisting of DEGs and DMGs. A total of 1,662 and 1,986 DEGs were identified between OA and normal human cartilage from the GSE113825 and GSE114007 dataset, respectively. A further screening revealed 292 genes with common differences between the two datasets. A total of 574 DMS containing 394 DMGs were observed between OA and normal cartilage. Integrative analysis revealed a corresponding subset of 15 genes. Of these, 6 genes were verified by reverse transcription-quantitative PCR, confirming that the mRNA expression of 5 genes (MAP1B, FNDC1, ANLN, SCNN1A and STC2) in OA cartilage was consistent with the mRNA expression from the analysis of the datasets. Upon treatment with the DNA methylation inhibitor 5-aza-2'-deoxycytidine, the mRNA levels of FNDC1 and SCNN1A were decreased, and no significant alteration in the mRNA levels of MAP1B, ANLN, KCNN4 and STC2 was observed. The incidence of differential methylation varied in subregions of the genome and the effects on transcription were associated with the distribution of DEGs across the genome. The regulation of this appears more complex than initially postulated. Combining the data on epigenetic differences of OA with the genome or transcriptome data for analysis may improve the understanding of the pathophysiological processes of OA. FNDC1 and SCNN1A may potentially be valuable biomarkers for OA.

Cathelicidin antimicrobial peptide (CAMP) gene promoter methylation induces chondrocyte apoptosis

Objective

The occurrence of osteoarthritis is related to genetic and environmental factors. Among them, the change of chondrocyte gene expression pattern regulated by epigenetic modification is an important participant. This study analyzed the effect of CAMP gene methylation on the level of oxidative stress and inflammation of chondrocytes.

Methods

We analyzed the changes of the transcriptome in the articular cartilage tissue of osteoarthritis (OA) patients from the GSE117999 dataset. The GSE48422 dataset was used to analyze the changes in the methylation level of osteoarthritis cells. Cell Counting Kit-8 (CCK-8) and flow cytometry analysis of short hairpin RNA (shRNA) silencing CAMP gene and 5-μM 5-Aza-2’-Deoxycytidine (AZA) treatment on the proliferation and apoptosis of Human chondrocytes osteoarthritis (HC-OA) cells. The Dichloro-dihydro-fluorescein diacetate (DCFH-DA) assay was used to detect the level of reactive oxygen species (ROS), and the expression level of inflammatory factors was analyzed by Western Blot.

Results

The expression of CAMP in cartilage tissue of OA patients was upregulated, and the level of methylation was downregulated. CAMP was highly expressed in osteoarthritis articular cartilage cells. Silencing CAMP inhibited the proliferation of HC-OA cells and promoted their apoptosis. CAMP gene methylation inhibited ROS levels and tumor necrosis factor-α (TNF-α) expression levels in HC-OA cells, and promoted transforming growth factor beta (TGF-β) expression. CAMP gene methylation inhibited the proliferation of HC-OA cells and promoted their apoptosis.

Conclusion

CAMP gene promoter methylation inhibits ROS levels and inflammation and induces chondrocyte apoptosis.

DNA Methylation Signatures of Bone Metabolism in Osteoporosis and Osteoarthritis Aging-Related Diseases: An Updated Review

DNA methylation is one of the most studied epigenetic mechanisms that play a pivotal role in regulating gene expression. The epigenetic component is strongly involved in aging-bone diseases, such as osteoporosis and osteoarthritis. Both are complex multi-factorial late-onset disorders that represent a globally widespread health problem, highlighting a crucial point of investigations in many scientific studies. In recent years, new findings on the role of DNA methylation in the pathogenesis of aging-bone diseases have emerged. The aim of this systematic review is to update knowledge in the field of DNA methylation associated with osteoporosis and osteoarthritis, focusing on the specific tissues involved in both pathological conditions.

Progress in the use of mesenchymal stromal cells for osteoarthritis treatment

LITERATURE REVIEW OF MSCS IN THE TREATMENT OF OSTEOARTHRITIS IN THE PAST FIVE YEARS: Osteoarthritis (OA) is one of the most common chronic joint diseases, with prominent symptoms caused by many factors. However, current medical interventions for OA have resulted in poor clinical outcomes, demonstrating that there are huge unmet medical needs in this area. Cell therapy has opened new avenues of OA treatment. Different sources of mesenchymal stromal cells (MSCs) may have different phenotypes and cellular functions. Pre-clinical and clinical studies have demonstrated the feasibility, safety and efficacy of MSC therapy. Mitogen-activated protein kinase, Wnt and Notch signaling pathways are involved in the chondrogenesis of MSC-mediated treatments. MSCs may also exert effective immunoregulatory and paracrine effects to stimulate tissue repair. Therapy with extracellular vesicles containing cytokines, which are secreted by MSCs, might be a potential treatment for OA.

Long non-coding RNA CIR inhibits chondrogenic differentiation of mesenchymal stem cells by epigenetically suppressing ATOH8 via methyltransferase EZH2

Background

Osteoarthritis (OA) is the most common articular disorder, leading to joint malfunction and disability. Although the incidence of OA is increasing globally, the treatment of OA is very limited. LncRNA CIR has been implicated in OA through unclear mechanisms. Here, we investigated the role of lncRNA CIR in chondrogenic differentiation.

Methods

Human umbilical-cord-derived mesenchymal stem cells (hUC-MSCs) were obtained from human umbilical cords. Flow cytometry was used to analyze the surface markers of hUC-MSCs. Various culture conditions and corresponding staining assays were employed to assess the differentiation abilities of hUC-MSC. qRT-PCR, western blot, and immunostaining were used to measure expression levels of related genes and proteins such as lncRNA CIR, ATOH8, EZH2, and H3K27me3. RNA immunoprecipitation assay, biotin pull-down, and chromatin immunoprecipitaion assay were performed to analyze the interactions of lncRNA CIR, EZH2, H3K27me3 and ATOH8 promoter.

Results

hUC-MSCs exhibited MSCs features and could differentiate into chondrocytes under specific conditions. LncRNA CIR was downregulated while ATOH8 was upregulated during the chondrogenic differentiation of hUC-MSCs. Knockdown lncRNA CIR or overexpression of ATOH8 promoted chondrogenic differentiation. Further, lncRNA CIR bound to EZH2 and repressed ATOH8 expression via EZH2-mediated H3K27me3, which promotes the methylation of ATOH8. Inhibition of ATOH8 reversed the effects of knockdown lncRNA CIR on chondrogenic differentiation.

Conclusion

LncRNA CIR suppresses chondrogenic differentiation of hUC-MSCs. Mechanistically, lncRNA CIR could inhibit ATOH8 expression that functions to promote chondrogenic differentiation through EZH2-mediated epigenetic modifications.

HNRNPD interacts with ZHX2 regulating the vasculogenic mimicry formation of glioma cells via linc00707/miR-651-3p/SP2 axis

Studies have found that RNA-binding proteins (RBPs) are dysfunctional and play a significant regulatory role in the development of glioma. Based on The Cancer Genome Atlas database and the previous studies, we selected heterogeneous nuclear ribonucleoprotein (HNRNPD) as the research candidate and sought its downstream targeted genes. In the present study, HNRNPD, linc00707, and specific protein 2 (SP2) were highly expressed, while zinc fingers and homeboxes 2 (ZHX2) and miR-651-3p were remarkedly downregulated in glioma tissues and cells. HNRNPD, linc00707, and SP2 knockdown or ZHX2 and miR-651-3p overexpression suppressed glioma cells proliferation, migration, and invasion and vasculogenic mimicry (VM) formation. Knockdown of HNRNPD increased the stability of ZHX2 mRNA. ZHX2 bound to the promoter region of linc00707 and negatively regulate its expression. Linc00707 could bind with miR-651-3p, while miR-651-3p bound to the 3' untranslated region (3'UTR) of SP2 mRNA to negatively regulate its expression. The transcription factor SP2 directly bound to the promoter regions of the VM formation-related proteins MMP2, MMP9, and VE-cadherin, playing a role in promoting transcription in order to regulate the VM formation ability of glioma cells.

TRAP1 Regulates Wnt/β-Catenin Pathway through LRP5/6 Receptors Expression Modulation

Wnt/β-Catenin signaling is involved in embryonic development, regeneration, and cellular differentiation and is responsible for cancer stemness maintenance. The HSP90 molecular chaperone TRAP1 is upregulated in 60-70% of human colorectal carcinomas (CRCs) and favors stem cells maintenance, modulating the Wnt/β-Catenin pathway and preventing β-Catenin phosphorylation/degradation. The role of TRAP1 in the regulation of Wnt/β-Catenin signaling was further investigated in human CRC cell lines, patient-derived spheroids, and CRC specimens. TRAP1 relevance in the activation of Wnt/β-Catenin signaling was highlighted by a TCF/LEF Cignal Reporter Assay in Wnt-off HEK293T and CRC HCT116 cell lines. Of note, this regulation occurs through the modulation of Wnt ligand receptors LRP5 and LRP6 that are both downregulated in TRAP1-silenced cell lines. However, while LRP5 mRNA is significantly downregulated upon TRAP1 silencing, LRP6 mRNA is unchanged, suggesting independent mechanisms of regulation by TRAP1. Indeed, LRP5 is regulated upon promoter methylation in CRC cell lines and human CRCs, whereas LRP6 is controlled at post-translational level by protein ubiquitination/degradation. Consistently, human CRCs with high TRAP1 expression are characterized by the co-upregulation of active β-Catenin, LRP5 and LRP6. Altogether, these data suggest that Wnt/β-Catenin signaling is modulated at multiple levels by TRAP1.

Changes in DNA methylation accompany changes in gene expression during chondrocyte hypertrophic differentiation in vitro

During osteoarthritis (OA), articular chondrocytes undergo phenotypic changes that resemble developmental patterns characteristic of growth plate chondrocytes. These phenotypic alterations lead to a hypertrophy-like phenotype characterized by altered production of extracellular matrix constituents and increased collagenase activity, which, in turn, results in cartilage destruction in OA disease. Recent studies have shown that the phenotypic instability and dysregulated gene expression in OA are associated with changes in DNA methylation patterns. Subsequent efforts have aimed to identify changes in DNA methylation with functional impact in OA disease, to potentially uncover therapeutic targets. Here, we paired an in vitro 3D/pellet culture system that mimics chondrocyte hypertrophy with RNA sequencing (RNA-Seq) and enhanced reduced representation of bisulfite sequencing (ERRBS) to identify transcriptomic and epigenomic changes in murine primary articular chondrocytes undergoing hypertrophy-like differentiation. We identified hypertrophy-associated changes in DNA methylation patterns in vitro. Integration of RNA-Seq and ERRBS datasets identified associations between changes in methylation and gene expression. Our integrative analyses showed that hypertrophic differentiation of articular chondrocytes is accompanied by transcriptomic and epigenomic changes in vitro. We believe that our integrative approaches have the potential to uncover new targets for therapeutic intervention.

Effects of 5-aza-2´-deoxycytidine on primary human chondrocytes from osteoarthritic patients

Chondrocytes, comparable to many cells from the connective tissue, dedifferentiate and end up in a similar fibroblastoid cell type, marked by the loss of the specific expression pattern. Here, chondrocytes isolated from osteoarthritic (OA) patients were investigated. The OA chondrocytes used in this work were not affected by the loss of specific gene expression upon cell culture. The mRNA levels of known cartilage markers, such as SOX5, SOX6, SOX9, aggrecan and proteoglycan 4, remained unchanged. Since chondrocytes from OA and healthy tissue show different DNA methylation patterns, the underlying mechanisms of cartilage marker maintenance were investigated with a focus on the epigenetic modification by DNA methylation. The treatment of dedifferentiated chondrocytes with the DNA methyltransferase inhibitor 5-aza-2´-deoxycytidine (5-aza-dC) displayed no considerable impact on the maintenance of marker gene expression observed in the dedifferentiated state, while the chondrogenic differentiation capacity was compromised. On the other hand, the pre-cultivation with 5-aza-dC improved the osteogenesis and adipogenesis of OA chondrocytes. Contradictory to these effects, the DNA methylation levels were not reduced after treatment for four weeks with 1 μM 5-aza-dC. In conclusion, 5-aza-dC affects the differentiation capacity of OA chondrocytes, while the global DNA methylation level remains stable. Furthermore, dedifferentiated chondrocytes isolated from late-stage OA patients represent a reliable cell source for in vitro studies and disease models without the need for additional alterations.

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next‐generation sequencing (ChIP‐seq) was used to assess a range of N‐terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP‐seq and cis‐regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis‐related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de‐methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell‐types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell‐types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross‐talk between different epigenetic mechanisms during chondrocyte differentiation.

Deep-reaching acceleration of global mean ocean circulation over the past two decades

Ocean circulation redistributes Earth's energy and water masses and influences global climate. Under historical greenhouse warming, regional ocean currents show diverse tendencies, but whether there is an emerging trend of the global mean ocean circulation system is not yet clear. Here, we show a statistically significant increasing trend in the globally integrated oceanic kinetic energy since the early 1990s, indicating a substantial acceleration of global mean ocean circulation. The increasing trend in kinetic energy is particularly prominent in the global tropical oceans, reaching depths of thousands of meters. The deep-reaching acceleration of the ocean circulation is mainly induced by a planetary intensification of surface winds since the early 1990s. Although possibly influenced by wind changes associated with the onset of a negative Pacific decadal oscillation since the late 1990s, the recent acceleration is far larger than that associated with natural variability, suggesting that it is principally part of a long-term trend.

DNA hypomethylation during MSC chondrogenesis occurs predominantly at enhancer regions

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DNA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated CpGs in osteoarthritis. Taken together we have associated DNA methylation levels with the chondrocyte phenotype. The consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.

Long non-coding XIST raises methylation of TIMP-3 promoter to regulate collagen degradation in osteoarthritic chondrocytes after tibial plateau fracture

Background

Hypermethylation of gene promoters has been regarded as an epigenetic regulator for gene inactivation in the development of several diseases. In the current study, we aimed to explore how long noncoding RNA X-inactive specific transcript (lncRNA XIST) function in collagen degradation in chondrocytes of osteoarthritis (OA) after tibial plateau fracture by regulating tissue inhibitor of metalloproteinase-3 (TIMP-3) promoter methylation.

Methods

In silico analysis was used to screen differentially expressed lncRNAs in cartilage tissues of OA. Chondrocytes were then successfully isolated from normal and OA cartilage tissues and identified, with the expressions of lncRNA XIST and TIMP-3 examined. The methylation levels of TIMP-3 promoter were determined by MS-PCR. The binding of lncRNA XIST to DNA methyltransferase and the binding of TIMP-3 promoter to DNA methyltransferase were determined by a series of experiments, including RIP, RNA pull-down, and ChIP assays.

Results

The differentially expressed lncRNA XIST was determined in OA. In addition, cartilage tissues of OA showed upregulation of lncRNA XIST and downregulation of TIMP-3. LncRNA XIST was primarily localized in the nucleus and was capable of binding to the promoter of TIMP-3. The silencing of lncRNA XIST decreased the methylation levels of TIMP-3 promoter and increased the expressions of TIMP-3, which consequently inhibited collagen degradation in OA chondrocytes. Furthermore, TIMP-3 over-expression reversed the effect of lncRNA XIST on collagen degradation in OA chondrocytes.

Conclusion

Collectively, lncRNA XIST raises collagen degradation in OA chondrocytes after tibial plateau fracture by accelerating the methylation of TIMP-3 promoter by recruiting DNA methyltransferase.

Genome-wide analysis of DNA methylation profile identifies differentially methylated loci associated with human intervertebral disc degeneration

BACKGROUND

Environmental and endogenous factors under genetic predisposition are considered to initiate the human intervertebral disc (IVD) degeneration. DNA methylation is an essential mechanism to ensure cell-specific gene expression for normal development and tissue stability. Aberrant epigenetic alterations play a pivotal role in several diseases, including osteoarthritis. However, epigenetic alternations, including DNA methylation, in IVD degeneration have not been evaluated. The purpose of this study was to comprehensively compare the genome-wide DNA methylation profiles of human IVD tissues, specifically nucleus pulpous (NP) tissues, with early and advanced stages of disc degeneration.

METHODS

Human NP tissues were used in this study. The samples were divided into two groups: early stage degeneration (n = 8, Pfirrmann's MRI grade: I-III) and advanced stage degeneration (n = 8, grade: IV). Genomic DNA was processed for genome-wide DNA methylation profiling using the Infinium MethylationEPIC BeadChip array. Extraction of raw methylation data, clustering and scatter plot of each group values of each sample were performed using a methylation module in GenomeStudio software. The identification of differentially methylated loci (DMLs) and the Gene Ontology (GO) analysis were performed using R software with the ChAMP package.

RESULTS

Unsupervised hierarchical clustering revealed that early and advanced stage degenerated IVD samples segregated into two main clusters by their DNA methylome. A total of 220 DMLs were identified between early and advanced disc degeneration stages. Among these, four loci were hypomethylated and 216 loci were hypermethylated in the advanced disc degeneration stage. The GO enrichment analysis of genes containing DMLs identified two significant GO terms for biological processes, hemophilic cell adhesion and cell-cell adhesion.

CONCLUSIONS

We conducted a genome-wide DNA methylation profile comparative study and observed significant differences in DNA methylation profiles between early and advanced stages of human IVD degeneration. These results implicate DNA methylation in the process of human IVD degeneration.

Racial differences in skeletal fragility but not osteoarthritis among women and men with cerebral palsy

Background

Adults with cerebral palsy (CP) have increased risk for skeletal fragility and osteoarthritis. However, racial differences in these outcomes have not been examined. Such knowledge could improve patient-specific clinical care for the prevention and management of these conditions. The purpose of this study was to determine if there were racial differences in the prevalence of osteoporosis, all-cause fracture, and osteoarthritis among young and middle-aged adults with CP.

Methods

Data from 2016 were extracted from a random 20% sample of the Medicare fee-for-service database. International Classification of Diseases, Tenth Revision, Clinical Modification codes were used to identify 18-64 year olds with CP, as well as osteoporosis, all-cause fracture, osteoarthritis, and neurodevelopmental and noncommunicable disease comorbidities.

Results

Of the 16,488 adults with CP, 13,334 were White, 2477 were Black, and 677 were Hispanic. The age-standardized prevalence of osteoporosis (women: 12.9%, 9.0%, 8.3%, respectively; men: 9.2%, 4.8%, 7.9%, respectively) and fracture (women: 7.4%, 4.2%, 9.9%; men: 6.0%, 2.3%, 6.0%) was lower for Black adults with CP compared to White adults with CP (all p < 0.005). No racial differences were observed for age-standardized prevalence of osteoarthritis (women: 13.6%, 14.4%, 9.6%; men: 10.7%, 10.4%, 12.7%). The racial differences between Black and White adults with CP remained even after adjusting for age, U.S. region, neurodevelopmental comorbidities, and several noncommunicable diseases for osteoporosis (women: OR = 0.66, 99.5% CI = 0.48-0.91; men: OR = 0.51, 99.5% CI = 0.35-0.75) and fracture (women: OR = 0.57, 99.5% CI = 0.37-0.89; men: OR = 0.39, 99.5% CI = 0.23-0.68).

Conclusions

Study findings suggest racial differences in skeletal fragility among young and middle-aged adults with CP, with White women and men with CP having greater risk compared to Black women and men with CP. This study found no evidence of racial differences in the prevalence of osteoarthritis.

Assessment of DNA Methylation Patterns in the Bone and Cartilage of a Nonhuman Primate Model of Osteoarthritis

OBJECTIVE

Osteoarthritis (OA) affects humans and several other animals. Thus, the mechanisms underlying this disorder, such as specific skeletal tissue DNA methylation patterns, may be evolutionary conserved. However, associations between methylation and OA have not been readily studied in nonhuman animals. Baboons serve as important models of disease and develop OA at rates similar to those in humans. Therefore, this study investigated the associations between methylation and OA in baboons to advance the evolutionary understanding of OA.

DESIGN

Trabecular bone and cartilage was collected from the medial condyles of adult female baboon femora, 5 with and 5 without knee OA. The Infinium HumanMethylation450 BeadChip (450K array) was used to identify DNA methylation patterns in these tissues.

RESULTS

Approximately 44% of the 450K array probes reliably align to the baboon genome, contain a CpG site of interest, and maintain a wide distribution throughout the genome. Of the 2 filtering methods tested, both identified significantly differentially methylated positions (DMPs) between healthy and OA individuals in cartilage tissues, and some of these patterns overlap with those previously identified in humans. Conversely, no DMPs were found between tissue types or between disease states in bone tissues.

CONCLUSIONS

Overall, the 450K array can be used to measure genome-wide DNA methylation in baboon tissues and identify significant associations with complex traits. The results of this study indicate that some DNA methylation patterns associated with OA are evolutionarily conserved, while others are not. This warrants further investigation in a larger and more phylogenetically diverse sample set.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Recent advances in understanding the regulation of metalloproteinases

Metalloproteinases remain important players in arthritic disease, in part because members of this large enzymatic family, namely matrix metalloproteinase-1 (MMP-1) and MMP-13, are responsible for the irreversible degradation of articular cartilage collagen. Although direct inhibition of MMPs fell out of vogue with the initial clinical disappointment of the first generation of compounds, interest in other mechanisms that control these important enzymes has always been maintained. Since these enzymes are critically important for tissue homeostasis, their expression and activity are tightly regulated at many levels, not just by direct inhibition by their endogenous inhibitors the tissue inhibitors of metalloproteinases (TIMPs). Focussing on MMP-13, we discuss recent work that highlights new discoveries in the transcriptional regulation of this enzyme, from defined promoter functional analysis to how more global technologies can provide insight into the enzyme’s regulation, especially by epigenetic mechanisms, including non-coding RNAs. In terms of protein regulation, we highlight recent findings into enzymatic cascades involved in MMP-13 regulation and activation. Importantly, we highlight a series of recent studies that describe how MMP-13 activity, and in fact that of other metalloproteinases, is in part controlled by receptor-mediated endocytosis. Together, these new discoveries provide a plethora of novel regulatory mechanisms, besides direct inhibition, which with renewed vigour could provide further therapeutic opportunities for regulating the activity of this class of important enzymes.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

Epigenetics and bone diseases

Owing to the development of new technologies, the epigenome, a second dimensional method for genome analysis has emerged. Epigenetic mechanisms, including DNA methylation, histone modifications and noncoding RNAs, regulate gene expression without changing the genetic sequence. These epigenetic mechanisms normally modulate gene expression, trans-generational effects and inherited expression states in various biological processes. Abnormal epigenetic patterns typically cause pathological conditions, including cancers, age-related diseases, and specific cartilage and bone diseases. Facing the rapidly developing epigenetic field, we reviewed epigenetic mechanisms and their involvement with the skeletal system and their role in skeletal development, homeostasis and degeneration. Finally, we discuss the prospects for the future of epigenetics.

Phenotypic instability of chondrocytes in osteoarthritis: on a path to hypertrophy

Articular chondrocytes are quiescent, fully differentiated cells responsible for the homeostasis of adult articular cartilage by maintaining cellular survival functions and the fine-tuned balance between anabolic and catabolic functions. This balance requires phenotypic stability that is lost in osteoarthritis (OA), a disease that affects and involves all joint tissues and especially impacts articular cartilage structural integrity. In OA, articular chondrocytes respond to the accumulation of injurious biochemical and biomechanical insults by shifting toward a degradative and hypertrophy-like state, involving abnormal matrix production and increased aggrecanase and collagenase activities. Hypertrophy is a necessary, transient developmental stage in growth plate chondrocytes that culminates in bone formation; in OA, however, chondrocyte hypertrophy is catastrophic and it is believed to initiate and perpetuate a cascade of events that ultimately result in permanent cartilage damage. Emphasizing changes in DNA methylation status and alterations in NF-κB signaling in OA, this review summarizes the data from the literature highlighting the loss of phenotypic stability and the hypertrophic differentiation of OA chondrocytes as central contributing factors to OA pathogenesis.

What did we learn from 'omics' studies in osteoarthritis

PURPOSE OF REVIEW

'Omics' technologies developed for the massive analysis of the major biologically relevant molecules (genes, proteins, metabolites) have been applied to the study of osteoarthritis (OA) for more than a decade.

RECENT FINDINGS

'Omics' studies have undoubtedly contributed to increase the knowledge on pathogenic processes related with OA and have provided hundreds to thousands of molecules that might have a putative biomarker utility for this disease.

SUMMARY

This review describes the most recent 'omics' studies in OA research, their conclusions, and discuss those remaining challenges. Still many validation studies must be performed in large and well-characterized cohorts for the translation of the findings from 'omics' strategies to clinical applications. The development of tools for the intelligent integration of 'omics' data with clinical and imaging information is also mandatory to take full profit of the work that has been already performed.

Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage

Aim

Osteoarthritis (OA) is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control the expression of genes and are likely to regulate the OA transcriptome. We performed integrative genomic analyses to define methylation-gene expression relationships in osteoarthritic cartilage.

Patients and Methods

Genome-wide DNA methylation profiling of articular cartilage from five patients with OA of the knee and five healthy controls was conducted using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, California). Other independent genome-wide mRNA expression profiles of articular cartilage from three patients with OA and three healthy controls were obtained from the Gene Expression Omnibus (GEO) database. Integrative pathway enrichment analysis of DNA methylation and mRNA expression profiles was performed using integrated analysis of cross-platform microarray and pathway software. Gene ontology (GO) analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID).

Results

We identified 1265 differentially methylated genes, of which 145 are associated with significant changes in gene expression, such as DLX5, NCOR2 and AXIN2 (all p-values of both DNA methylation and mRNA expression < 0.05). Pathway enrichment analysis identified 26 OA-associated pathways, such as mitogen-activated protein kinase (MAPK) signalling pathway (p = 6.25 × 10-4), phosphatidylinositol (PI) signalling system (p = 4.38 × 10-3), hypoxia-inducible factor 1 (HIF-1) signalling pathway (p = 8.63 × 10-3 pantothenate and coenzyme A (CoA) biosynthesis (p = 0.017), ErbB signalling pathway (p = 0.024), inositol phosphate (IP) metabolism (p = 0.025), and calcium signalling pathway (p = 0.032).

Conclusion

We identified a group of genes and biological pathwayswhich were significantly different in both DNA methylation and mRNA expression profiles between patients with OA and controls. These results may provide new clues for clarifying the mechanisms involved in the development of OA.

Cite this article: A. He, Y. Ning, Y. Wen, Y. Cai, K. Xu, Y. Cai, J. Han, L. Liu, Y. Du, X. Liang, P. Li, Q. Fan, J. Hao, X. Wang, X. Guo, T. Ma, F. Zhang. Use of integrative epigenetic and mRNA expression analyses to identify significantly changed genes and functional pathways in osteoarthritic cartilage. Bone Joint Res 2018;7:343–350. DOI: 10.1302/2046-3758.75.BJR-2017-0284.R1.

Osteoarthritis year in review 2017: genetics and epigenetics

OBJECTIVE

The purpose of this review is to describe highlights from original research publications related to osteoarthritis (OA), epigenetics and genomics with the intention of recognising significant advances.

DESIGN

To identify relevant papers a Pubmed literature search was conducted for articles published between April 2016 and April 2017 using the search terms 'osteoarthritis' together with 'genetics', 'genomics', 'epigenetics', 'microRNA', 'lncRNA', 'DNA methylation' and 'histone modification'.

RESULTS

The search term OA generated almost 4000 references. Publications using the combination of descriptors OA and genetics provided the most references (82 references). However this was reduced compared to the same period in the previous year; 8.1-2.1% (expressed as a percentage of the total publications combining the terms OA and genetics). Publications combining the terms OA with genomics (29 references), epigenetics (16 references), long non-coding RNA (lncRNA) (11 references; including the identification of novel lncRNAs in OA), DNA methylation (21 references), histone modification (3 references) and microRNA (miR) (79 references) were reviewed. Potential OA therapeutics such as histone deacetylase (HDAC) inhibitors have been identified. A number of non-coding RNAs may also provide targets for future treatments.

CONCLUSION

There continues to be a year on year increase in publications researching miRs in OA (expressed as a percentage of the total publications), with a doubling over the last 4 years. An overview on the last year's progress within the fields of epigenetics and genomics with respect to OA will be given.

Transcriptional profiling of articular cartilage in a porcine model of early post-traumatic osteoarthritis

To identify the molecular pathophysiology present in early post-traumatic osteoarthritis (PTOA), the transcriptional profile of articular cartilage and its response to surgical PTOA induction were determined. Thirty six Yucatan minipigs underwent anterior cruciate ligament (ACL) transection and were randomly assigned in equal numbers to no further treatment, reconstruction or ligament repair. Cartilage was harvested at 1 and 4 weeks post-operatively and histology and RNA-sequencing were performed and compared to controls. Microscopic cartilage scores significantly worsened at 1 (p = 0.028) and 4 weeks (p = 0.001) post-surgery relative to controls, but did not differ between untreated, reconstruction or repair groups. Gene expression after ACL reconstruction and ACL transection were similar, with only 0.03% (including SERPINB7 and CR2) and 0.2% of transcripts (including INHBA) differentially expressed at 1 and 4 weeks respectively. COL2A1, COMP, SPARC, CHAD, and EF1ALPHA were the most highly expressed non ribosomal, non mitochondrial genes in the controls and remained abundant after surgery. A total of 1,275 genes were differentially expressed between 1 and 4 weeks post-surgery. With the treatment groups pooled, 682 genes were differentially expressed at both time-points, with the most significant changes observed in MMP1, COCH, POSTN, CYTL1, and PTGFR. This study confirmed the development of a microscopic PTOA stage after ACL surgery in the porcine model. Upregulation of multiple proteases (including MMP1 and ADAMTS4) were found; however, the level of expression remained orders of magnitude below that of extracellular matrix protein-coding genes (including COL2A1 and ACAN). In summary, genes with established roles in PTOA as well as novel targets for specific intervention were identified. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:318-329, 2018.

TWIST1 induces MMP3 expression through up-regulating DNA hydroxymethylation and promotes catabolic responses in human chondrocytes

The objective was to investigate the levels of TWIST1 in normal and OA cartilage and examine its role in regulating gene expression in chondrocytes. Human cartilage tissues and chondrocytes were obtained at autopsy from normal knee joints and from OA-affected joints at the time of total knee arthroplasty. TWIST1 expression was increased in human OA knee cartilage compared to normal knee cartilage. TWIST1 induced matrix metalloproteinase 3 (MMP3) expression without direct binding to MMP3 promoter and increased the 5-hydroxymethylcytosine (5hmC) level at the MMP3 promoter. The effect of TWIST1 on expression of TET family (TET1, 2 and 3) was measured in stable TWIST1 transfected TC28 cells, and TET1 expression was up-regulated. TWIST1 dependent upregulation of Mmp3 expression was suppressed in Tet triple KO fibroblast derived from mouse ES cells. Increased TWIST1 expression is a feature of OA-affected cartilage. We identified a novel mechanism of catabolic reaction where TWIST1 up-regulates MMP3 expression by enriching 5hmC levels at the MMP3 promoter via TET1 induction. These findings implicate TWIST1 as an important factor regulating OA related gene expression. Clarifying epigenetic mechanisms of 5hmC induced by TWIST1 is a critical molecule to understanding OA pathogenesis.

Osteoarthritis year in review 2016: genetics, genomics and epigenetics

The purpose of this narrative review is to provide an overview of last year's publications in the field of genetics, genomics and epigenetics in the osteoarthritis (OA) field. Major themes arising from a Pubmed search on (epi)genetics in OA were identified. In addition, general developments in the fast evolving field of (epi)genetics are reviewed and relevance for the OA field is summarized. In the last 5 years, a number of genome-wide association studies have identified a modest number of genetic loci associated to OA. Continued functional research into these DNA variants is showing putative biological mechanisms underlying these associations. Over the last year, no additional large genome-wide association studies were published, but there clearly remains much to be discovered in the OA genetic field. A lot of research has been done into the epigenetics of OA over the last year. Several genome-wide screens examining the methylome of osteoarthritic cartilage were done. Pathway analysis confirmed deregulation of developmental and extracellular pathways in OA cartilage. Over the last year many microRNAs (miRNAs) have been identified that potentially play important roles in cartilage homeostasis and/or OA process. Continued research will learn whether these identified miRNAs are truly causal and can be used in clinical applications. Many of the epigenetic findings need further confirmation, but they highlight potential novel pathways involved in cartilage biology and OA.

Analysis of genetics and DNA methylation in osteoarthritis: What have we learnt about the disease?

Osteoarthritis (OA) is a chronic musculoskeletal disease characterised by the destruction of articular cartilage, synovial inflammation and bone remodelling. Disease aetiology is complex and highly heritable, with genetic variation estimated to contribute to 50% of OA occurrence. Epigenetic alterations, including DNA methylation changes, have also been implicated in OA pathophysiology. This review examines what genetic and DNA methylation studies have taught us about the genes and pathways involved in OA pathology. The influence of DNA methylation on the molecular mechanisms underlying OA genetic risk and the consequence of this interaction on disease susceptibility and penetrance are also discussed.