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

Epigenome editing revealed the role of DNA methylation of T-DMR/CpG island shore on Runx2 transcription

RUNX2 is a transcription factor crucial for bone formation. Mutant mice with varying levels of Runx2 expression display dosage-dependent skeletal abnormalities, underscoring the importance of Runx2 dosage control in skeletal formation. RUNX2 activity is regulated by several molecular mechanisms, including epigenetic modification such as DNA methylation. In this study, we investigated whether targeted repressive epigenome editing including hypermethylation to the Runx2-DMR/CpG island shore could influence Runx2 expression using Cas9-based epigenome-editing tools. Through the transient introduction of CRISPRoff-v2.1 and gRNAs targeting Runx2-DMR into MC3T3-E1 cells, we successfully induced hypermethylation of the region and concurrently reduced Runx2 expression during osteoblast differentiation. Although the epigenome editing of Runx2-DMR did not impact the expression of RUNX2 downstream target genes, these results indicate a causal relationship between the epigenetic status of the Runx2-DMR and Runx2 transcription. Additionally, we observed that hypermethylation of the Runx2-DMR persisted for at least 24 days under growth conditions but decreased during osteogenic differentiation, highlighting an endogenous DNA demethylation activity targeting the Runx2-DMR during the differentiation process. In summary, our study underscore the usefulness of the epigenome editing technology to evaluate the function of endogenous genetic elements and revealed that the Runx2-DMR methylation is actively regulated during osteoblast differentiation, subsequently could influence Runx2 expression.

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Potential therapeutic strategies for osteoarthritis via CRISPR/Cas9 mediated gene editing

Osteoarthritis (OA) is an incapacitating and one of the most common physically degenerative conditions with an assorted etiology and a highly complicated molecular mechanism that to date lacks an efficient treatment. The capacity to design biological networks and accurately modify existing genomic sites holds an apt potential for applications across medical and biotechnological sciences. One of these highly specific genomes editing technologies is the CRISPR/Cas9 mechanism, referred to as the clustered regularly interspaced short palindromic repeats, which is a defense mechanism constituted by CRISPR associated protein 9 (Cas9) directed by small non-coding RNAs (sncRNA) that bind to target DNA through Watson-Crick base pairing rules where subsequent repair of the target DNA is initiated. Up-to-date research has established the effectiveness of the CRISPR/Cas9 mechanism in targeting the genetic and epigenetic alterations in OA by suppressing or deleting gene expressions and eventually distributing distinctive anti-arthritic properties in both in vitro and in vivo osteoarthritic models. This review aims to epitomize the role of this high-throughput and multiplexed gene editing method as an analogous therapeutic strategy that could greatly facilitate the clinical development of OA-related treatments since it's reportedly an easy, minimally invasive technique, and a comparatively less painful method for osteoarthritic patients.

Differences of RUNX2 gene promoter methylation and transcription level in ankylosing spondylitis

BACKGROUND AND OBJECTIVE

Ankylosing spondylitis is a refractory immune disease that seriously affects the life and work of patients. Epigenetic modifications, especially DNA methylation, have become a research hotspot in complex diseases. We aim to explore the changes in runt-related transcription factor 2 (RUNX2) gene promoter methylation and transcription level in AS.

METHOD

We detected the RUNX2 gene promoter methylation in 83 AS patients and 83 healthy controls (HCs), then inspected the mRNA difference of RUNX2 between 30 AS patients and 30 HCs by the quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR).

RESULTS

The RUNX2 gene promoter was hypomethylated in AS patients compared to HCs (p < .001). The research involved 4 CpG regions and 74 CpG sites of RUNX2, of which CpG-2, CpG-4 regions, and 18 CpG sites have been differentially methylated. The CpG-4 island methylation was negatively correlated with C-reactive protein (p < .05) in AS patients. In the qRT-PCR validation phase, the mRNA level of RUNX2 in AS patients was significantly higher than HCs (p < .05), and in AS patients who were treated with biologics, the methylation level of CpG-2 island showed a negative correlation to mRNA (p < .05). ROC results indicated that RUNX2 methylation and its transcription level have good potential to distinguish AS patients from HCs.

CONCLUSION

The RUNX2 gene promoter was hypomethylated in AS patients. Meanwhile, the qRT-PCR verified the up-regulated expression on the transcription level, suggesting the abnormal methylation of RUNX2 contributes to the pathogenesis of AS.

Role of epigenetics and the transcription factor Sp1 in the expression of the D prostanoid receptor 1 in human cartilage

D prostanoid receptor 1 (DP1), a prostaglandin D2 receptor, plays a central role in the modulation of inflammation and cartilage metabolism. We have previously shown that activation of DP1 signaling downregulated catabolic responses in cultured chondrocytes and was protective in mouse osteoarthritis (OA). However, the mechanisms underlying its transcriptional regulation in cartilage remained poorly understood. In the present study, we aimed to characterize the human DP1 promoter and the role of DNA methylation in DP1 expression in chondrocytes. In addition, we analyzed the expression level and methylation status of the DP1 gene promoter in normal and OA cartilage. Deletion and site-directed mutagenesis analyses identified a minimal promoter region (-250/-120) containing three binding sites for specificity protein 1 (Sp1). Binding of Sp1 to the DP1 promoter was confirmed using electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP) assays. Treatment with the Sp1 inhibitor mithramycin A reduced DP1 promoter activity and DP1 mRNA expression. Inhibition of DNA methylation by 5-Aza-2'-deoxycytidine upregulated DP1 expression, and in vitro methylation reduced the DP1 promoter activity. Neither the methylation status of the DP1 promoter nor the DP1 expression level were different between normal and OA cartilage. In conclusion, our results suggest that the transcription factor Sp1 and DNA methylation are important determinants of DP1 transcription regulation. They also suggest that the methylation status and expression level of DP1 are not altered in OA cartilage. These findings will improve our understanding of the regulatory mechanisms of DP1 transcription and may facilitate the development of intervention strategies involving DP1.

Mediation of the Same Epigenetic and Transcriptional Effect by Independent Osteoarthritis Risk-Conferring Alleles on a Shared Target Gene, COLGALT2

OBJECTIVE

Over 100 DNA variants have been associated with osteoarthritis (OA), including rs1046934, located within a linkage disequilibrium block encompassing part of COLGALT2 and TSEN15. The present study was undertaken to determine the target gene(s) and the mechanism of action of the OA locus using human fetal cartilage, cartilage from OA and femoral neck fracture arthroplasty patients, and a chondrocyte cell model.

METHODS

Genotyping and methylation array data of DNA from human OA cartilage samples (n = 87) were used to determine whether the rs1046934 genotype is associated with differential DNA methylation at proximal CpGs. Results were replicated in DNA from human arthroplasty (n = 132) and fetal (n = 77) cartilage samples using pyrosequencing. Allelic expression imbalance (AEI) measured the effects of genotype on COLGALT2 and TSEN15 expression. Reporter gene assays and epigenetic editing determined the functional role of regions harboring differentially methylated CpGs. In silico analyses complemented these experiments.

RESULTS

Three differentially methylated CpGs residing within regulatory regions were detected in the human OA cartilage array data, and 2 of these were replicated in human arthroplasty and fetal cartilage. AEI was detected for COLGALT2 and TSEN15, with associations between expression and methylation for COLGALT2. Reporter gene assays confirmed that the CpGs are in chondrocyte enhancers, with epigenetic editing results directly linking methylation with COLGALT2 expression.

CONCLUSION

COLGALT2 is a target of this OA locus. We previously characterized another OA locus, marked by rs11583641, that independently targets COLGALT2. The genotype of rs1046934, like rs11583641, mediates its effect by modulating expression of COLGALT2 via methylation changes to CpGs located in enhancers. Although the single-nucleotide polymorphisms, CpGs, and enhancers are distinct between the 2 independent OA risk loci, their effect on COLGALT2 is the same. COLGALT2 is the target of independent OA risk loci sharing a common mechanism of action.

Osteoarthritis genetic risk acting on the galactosyltransferase gene COLGALT2 has opposing functional effects in articulating joint tissues

BACKGROUND

Investigation of cartilage and chondrocytes has revealed that the osteoarthritis risk marked by the independent DNA variants rs11583641 and rs1046934 mediate  their effects by decreasing the methylation status of CpG dinucleotides in enhancers and increasing the expression of shared target gene COLGALT2. We set out to investigate if these functional effects operate in a non-cartilaginous joint tissue.

METHODS

Nucleic acids were extracted from the synovium of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at CpGs within the COLGALT2 enhancers. CpGs were tested for enhancer effects using a synovial cell line and a reporter gene assay. DNA methylation was altered using epigenetic editing, with the impact on gene expression determined using quantitative polymerase chain reaction. In silico analysis complemented laboratory experiments.

RESULTS

The rs1046934 genotype did not associate with DNA methylation or COLGALT2 expression in the synovium, whereas the rs11583641 genotype did. Surprisingly, the effects for rs11583641 were opposite to those previously observed in cartilage. Epigenetic editing in synovial cells revealed that enhancer methylation is causally linked to COLGALT2 expression.

CONCLUSIONS

This is the first direct demonstration for osteoarthritis genetic risk of a functional link between DNA methylation and gene expression operating in opposite directions between articular joint tissues. It highlights pleiotropy in the action of osteoarthritis risk and provides a cautionary note in the application of future genetically based osteoarthritis therapies: an intervention that decreases the detrimental effect of a risk allele in one joint tissue may inadvertently increase its detrimental effect in another joint tissue.

SARS-CoV-2 airway infection results in the development of somatosensory abnormalities in a hamster model

Although largely confined to the airways, SARS-CoV-2 infection has been associated with sensory abnormalities that manifest in both acute and chronic phenotypes. To gain insight on the molecular basis of these sensory abnormalities, we used the golden hamster model to characterize and compare the effects of infection with SARS-CoV-2 and influenza A virus (IAV) on the sensory nervous system. We detected SARS-CoV-2 transcripts but no infectious material in the cervical and thoracic spinal cord and dorsal root ganglia (DRGs) within the first 24 hours of intranasal virus infection. SARS-CoV-2-infected hamsters exhibited mechanical hypersensitivity that was milder but prolonged compared with that observed in IAV-infected hamsters. RNA sequencing analysis of thoracic DRGs 1 to 4 days after infection suggested perturbations in predominantly neuronal signaling in SARS-CoV-2-infected animals as opposed to type I interferon signaling in IAV-infected animals. Later, 31 days after infection, a neuropathic transcriptome emerged in thoracic DRGs from SARS-CoV-2-infected animals, which coincided with SARS-CoV-2-specific mechanical hypersensitivity. These data revealed potential targets for pain management, including the RNA binding protein ILF3, which was validated in murine pain models. This work elucidates transcriptomic signatures in the DRGs triggered by SARS-CoV-2 that may underlie both short- and long-term sensory abnormalities.

Targeted Resequencing of Otosclerosis Patients from Different Populations Replicates Results from a Previous Genome-Wide Association Study

Otosclerosis is one of the most common causes of hearing loss in young adults. It has a prevalence of 0.3-0.4% in the European population. Clinical symptoms usually occur between the second and fifth decade of life. Different studies have been performed to unravel the genetic architecture of the disease. Recently, a genome-wide association study (GWAS) identified 15 novel risk loci and replicated the regions of three previously reported candidate genes. In this study, seven candidate genes from the GWAS were resequenced using single molecule molecular inversion probes (smMIPs). smMIPs were used to capture the exonic regions and the 3' and 5' untranslated regions (UTR). Discovered variants were tested for association with the disease using single variant and gene-based association analysis. The single variant results showed that 13 significant variants were associated with otosclerosis. Associated variants were found in five of the seven genes studied here, including AHSG, LINC01482, MARK3, SUPT3H and RELN. Conversely, burden testing did not show a major role of rare variants in the disease. In conclusion, this study was able to replicate five out of seven candidate genes reported in the previous GWAS. This association is likely mainly driven by common variants.

SARS-CoV-2 Airway Infection Results in Time-dependent Sensory Abnormalities in a Hamster Model

Despite being largely confined to the airways, SARS-CoV-2 infection has been associated with sensory abnormalities that manifest in both acute and long-lasting phenotypes. To gain insight on the molecular basis of these sensory abnormalities, we used the golden hamster infection model to characterize the effects of SARS-CoV-2 versus Influenza A virus (IAV) infection on the sensory nervous system. Efforts to detect the presence of virus in the cervical/thoracic spinal cord and dorsal root ganglia (DRGs) demonstrated detectable levels of SARS-CoV-2 by quantitative PCR and RNAscope uniquely within the first 24 hours of infection. SARS-CoV-2-infected hamsters demonstrated mechanical hypersensitivity during acute infection; intriguingly, this hypersensitivity was milder, but prolonged when compared to IAV-infected hamsters. RNA sequencing (RNA-seq) of thoracic DRGs from acute infection revealed predominantly neuron-biased signaling perturbations in SARS-CoV-2-infected animals as opposed to type I interferon signaling in tissue derived from IAV-infected animals. RNA-seq of 31dpi thoracic DRGs from SARS-CoV-2-infected animals highlighted a uniquely neuropathic transcriptomic landscape, which was consistent with substantial SARS-CoV-2-specific mechanical hypersensitivity at 28dpi. Ontology analysis of 1, 4, and 30dpi RNA-seq revealed novel targets for pain management, such as ILF3. Meta-analysis of all SARS-CoV-2 RNA-seq timepoints against preclinical pain model datasets highlighted both conserved and unique pro-nociceptive gene expression changes following infection. Overall, this work elucidates novel transcriptomic signatures triggered by SARS-CoV-2 that may underlie both short- and long-term sensory abnormalities while also highlighting several therapeutic targets for alleviation of infection-induced hypersensitivity.

One Sentence Summary

SARS-CoV-2 infection results in an interferon-associated transcriptional response in sensory tissues underlying time-dependent hypersensitivity.

Identification of TMEM129, encoding a ubiquitin-protein ligase, as an effector gene of osteoarthritis genetic risk

BACKGROUND

Osteoarthritis is highly heritable and genome-wide studies have identified single nucleotide polymorphisms (SNPs) associated with the disease. One such locus is marked by SNP rs11732213 (T > C). Genotype at rs11732213 correlates with the methylation levels of nearby CpG dinucleotides (CpGs), forming a methylation quantitative trait locus (mQTL). This study investigated the regulatory activity of the CpGs to identify a target gene of the locus.

METHODS

Nucleic acids were extracted from the articular cartilage of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at 14 CpGs within a 259-bp interval. CpGs were tested for enhancer effects in immortalised chondrocytes using a reporter gene assay. DNA methylation at the locus was altered using targeted epigenome editing, with the impact on gene expression determined using quantitative polymerase chain reaction.

RESULTS

rs11732213 genotype correlated with DNA methylation at nine CpGs, which formed a differentially methylated region (DMR), with the osteoarthritis risk allele T corresponding to reduced levels of methylation. The DMR acted as an enhancer and demethylation of the CpGs altered expression of TMEM129. Allelic imbalance in TMEM129 expression was identified in cartilage, with under-expression of the risk allele.

CONCLUSIONS

TMEM129 is a target of osteoarthritis genetic risk at this locus. Genotype at rs11732213 impacts DNA methylation at the enhancer, which, in turn, modulates TMEM129 expression. TMEM129 encodes an enzyme involved in protein degradation within the endoplasmic reticulum, a process previously implicated in osteoarthritis. TMEM129 is a compelling osteoarthritis susceptibility target.

Epigenetic Regulation in Knee Osteoarthritis

Osteoarthritis (OA) is a complicated disease with both hereditary and environmental causes. Despite an increase in reports of possible OA risk loci, it has become clear that genetics is not the sole cause of osteoarthritis. Epigenetics, which can be triggered by environmental influences and result in transcriptional alterations, may have a role in OA pathogenesis. The majority of recent research on the epigenetics of OA has been focused on DNA methylation, histone modification, and non-coding RNAs. However, this study will explore epigenetic regulation in OA at the present stage. How genetics, environmental variables, and epigenetics interact will be researched, shedding light for future studies. Their possible interaction and control processes open up new avenues for the development of innovative osteoarthritis treatment and diagnostic techniques.

An epigenome-wide view of osteoarthritis in primary tissues

Osteoarthritis is a complex degenerative joint disease. Here, we investigate matched genotype and methylation profiles of primary chondrocytes from macroscopically intact (low-grade) and degraded (high-grade) osteoarthritis cartilage and from synoviocytes collected from 98 osteoarthritis-affected individuals undergoing knee replacement surgery. We perform an epigenome-wide association study of knee cartilage degeneration and report robustly replicating methylation markers, which reveal an etiologic mechanism linked to the migration of epithelial cells. Using machine learning, we derive methylation models of cartilage degeneration, which we validate with 82% accuracy in independent data. We report a genome-wide methylation quantitative trait locus (mQTL) map of articular cartilage and synovium and identify 18 disease-grade-specific mQTLs in osteoarthritis cartilage. We resolve osteoarthritis GWAS loci through causal inference and colocalization analyses and decipher the epigenetic mechanisms that mediate the effect of genotype on disease risk. Together, our findings provide enhanced insights into epigenetic mechanisms underlying osteoarthritis in primary tissues.

Genetics of osteoarthritis

Osteoarthritis genetics has been transformed in the past decade through the application of large-scale genome-wide association scans. So far, over 100 polymorphic DNA variants have been associated with this common and complex disease. These genetic risk variants account for over 20% of osteoarthritis heritability and the vast majority map to non-protein coding regions of the genome where they are presumed to act by regulating the expression of target genes. Statistical fine mapping, in silico analyses of genomics data, and laboratory-based functional studies have enabled the identification of some of these targets, which encode proteins with diverse roles, including extracellular signaling molecules, intracellular enzymes, transcription factors, and cytoskeletal proteins. A large number of the risk variants correlate with epigenetic factors, in particular cartilage DNA methylation changes in cis, implying that epigenetics may be a conduit through which genetic effects on gene expression are mediated. Some of the variants also appear to have been selected as humans adapted to bipedalism, suggesting that a proportion of osteoarthritis genetic susceptibility results from antagonistic pleiotropy, with risk variants having a positive role in joint formation but a negative role in the long-term health of the joint. Although data from an osteoarthritis genetic study has not yet directly led to a novel treatment, some of the osteoarthritis associated genes code for proteins that have available therapeutics. Genetic investigations are therefore revealing fascinating fundamental insights into osteoarthritis and can expose options for translational intervention.

Kindlin-2 preserves integrity of the articular cartilage to protect against osteoarthritis

Osteoarthritis (OA) is an aging-related degenerative joint disease with a poorly defined mechanism. Here we report that kindlin-2 is highly expressed in articular chondrocytes and downregulated in the degenerated cartilage of aged mice and patients with OA. Kindlin-2 deletion in articular chondrocytes leads to spontaneous OA and exacerbates instability-induced OA lesions in adult mice. Kindlin-2 deficiency promotes mitochondrial oxidative stress and activates Stat3, leading to Runx2-mediated chondrocyte catabolism. Pharmacological inhibition of Stat3 activation or genetic ablation of Stat3 in chondrocytes reverses aberrant accumulation of Runx2 and extracellular-matrix-degrading enzymes and limits OA deteriorations caused by kindlin-2 deficiency. Deleting Runx2 in chondrocytes reverses structural changes and OA lesions caused by kindlin-2 deletion without downregulating p-Stat3. Intra-articular injection of AAV5-kindlin-2 decelerates progression of aging- and instability-induced knee joint OA in mice. Collectively, we identify a pathway consisting of kindlin-2, Stat3 and Runx2 in articular chondrocytes that is responsible for maintaining articular cartilage integrity and define a potential therapeutic target for OA.

Translating osteoarthritis genetics research: challenging times ahead

The ultimate goal of molecular genetic studies of human diseases is to translate the discoveries for patient benefit. For diseases that lack licensed disease-modifying therapeutics, such as osteoarthritis (OA), the need is acute. OA is polygenic and affects older individuals, with a recent genome-wide study of over 800 000 individuals adding 52 novel association signals to those already reported on for this common arthritis. Many of the predicted effector genes of these signals encode proteins that are targets of drugs for other indications, highlighting repurposing opportunities. Here, the potential for OA genetic data to translate is discussed, including whether the developmental origin of OA will limit the application of genetic risk data for disease-modification purposes.

Genetic and Epigenetic Interplay Within a COLGALT2 Enhancer Associated With Osteoarthritis

OBJECTIVE

The osteoarthritis (OA)-associated single-nucleotide polymorphism (SNP) rs11583641 is located in COLGALT2, encoding a posttranslational modifier of collagen. In cartilage, the SNP genotype correlates with DNA methylation in a putative enhancer. This study was undertaken to characterize the mechanistic relationship between rs11583641, the putative enhancer, and COLGALT2 expression using cartilage samples from human patients and a chondrocyte cell model.

METHODS

Nucleic acids were extracted from articular cartilage samples obtained from patients with OA (n = 137). Samples were genotyped, and DNA methylation was quantified at 12 CpGs using pyrosequencing. The putative enhancer was deleted in Tc28a2 chondrocytes using clustered regularly interspaced short palindromic repeat/Cas9, and the impact on nearby gene expression was determined using real-time quantitative polymerase chain reaction. Targeted modulation of the epigenome using catalytically dead Cas9 (dCas9) constructs fused to DNA methyltransferase 3a or ten-eleven translocase 1 allowed for the investigation of a causal relationship between DNA methylation and enhancer activity.

RESULTS

The genotype at rs11583641 correlated with DNA methylation at 3 CpGs, and the presence of the OA risk allele, C, corresponded to reduced levels of methylation. Deletion of the enhancer resulted in a 2.7-fold reduction in COLGALT2 expression. Targeted methylation and demethylation of the CpGs had antagonistic effects on COLGALT2 expression. An allelic imbalance in the expression of COLGALT2 was identified in the cartilage from patients with OA, with relative overexpression of the OA risk allele. Allelic expression ratios correlated with DNA methylation at 4 CpGs.

CONCLUSION

COLGALT2 is a target of OA genetic risk at this locus. The genotype at rs11583641 impacts DNA methylation in a gene enhancer, which, in turn, modulates COLGALT2 expression. COLGALT2 encodes an enzyme that initiates posttranslational glycosylation of collagens and is therefore a compelling OA susceptibility target.

Genetic and Epigenetic Fine-Tuning of TGFB1 Expression Within the Human Osteoarthritic Joint

OBJECTIVE

Osteoarthritis (OA) is an age-related disease characterized by articular cartilage degeneration. It is largely heritable, and genetic screening has identified single-nucleotide polymorphisms (SNPs) marking genomic risk loci. One such locus is marked by the G>A SNP rs75621460, downstream of TGFB1. This gene encodes transforming growth factor β1, the correct expression of which is essential for cartilage maintenance. This study investigated the regulatory activity of rs75621460 to characterize its impact on TGFB1 expression in disease-relevant patient samples (n = 319) and in Tc28a2 immortalized chondrocytes.

METHODS

Articular cartilage samples from human patients were genotyped, and DNA methylation levels were quantified using pyrosequencing. Gene reporter and electrophoretic mobility shift assays were used to determine differential nuclear protein binding to the region. The functional impact of DNA methylation on TGFB1 expression was tested using targeted epigenome editing.

RESULTS

The analyses showed that SNP rs75621460 was located within a TGFB1 enhancer region, and the OA risk allele A altered transcription factor binding, with decreased enhancer activity. Protein complexes binding to A (but not G) induced DNA methylation at flanking CG dinucleotides. Strong correlations between patient DNA methylation levels and TGFB1 expression were observed, with directly opposing effects in the cartilage and the synovium at this locus. This demonstrated biologic pleiotropy in the impact of the SNP within different tissues of the articulating joint.

CONCLUSION

The OA risk SNP rs75621460 impacts TGFB1 expression by modulating the function of a gene enhancer. We propose a mechanism by which the SNP impacts enhancer function, providing novel biologic insight into one mechanism of OA genetic risk, which may facilitate the development of future pharmacologic therapies.

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.

Transcription Factors in Cartilage Homeostasis and Osteoarthritis

Osteoarthritis (OA) is the most common degenerative joint disease, and it is characterized by articular cartilage loss. In part, OA is caused by aberrant anabolic and catabolic activities of the chondrocyte, the only cell type present in cartilage. These chondrocyte activities depend on the intra- and extracellular signals that the cell receives and integrates into gene expression. The key proteins for this integration are transcription factors. A large number of transcription factors exist, and a better understanding of the transcription factors activated by the various signaling pathways active during OA can help us to better understand the complex etiology of OA. In addition, establishing such a profile can help to stratify patients in different subtypes, which can be a very useful approach towards personalized therapy. In this review, we discuss crucial transcription factors for extracellular matrix metabolism, chondrocyte hypertrophy, chondrocyte senescence, and autophagy in chondrocytes. In addition, we discuss how insight into these factors can be used for treatment purposes.

Emerging Gene-Editing Modalities for Osteoarthritis

Osteoarthritis (OA) is a pathological degenerative condition of the joints that is widely prevalent worldwide, resulting in significant pain, disability, and impaired quality of life. The diverse etiology and pathogenesis of OA can explain the paucity of viable preventive and disease-modifying strategies to counter it. Advances in genome-editing techniques may improve disease-modifying solutions by addressing inherited predisposing risk factors and the activity of inflammatory modulators. Recent progress on technologies such as CRISPR/Cas9 and cell-based genome-editing therapies targeting the genetic and epigenetic alternations in OA offer promising avenues for early diagnosis and the development of personalized therapies. The purpose of this literature review was to concisely summarize the genome-editing options against chronic degenerative joint conditions such as OA with a focus on the more recently emerging modalities, especially CRISPR/Cas9. Future advancements in novel genome-editing therapies may improve the efficacy of such targeted treatments.

Osteoarthritis treatment with a novel nutraceutical acetylated ligstroside aglycone, a chemically modified extra-virgin olive oil polyphenol

Recent studies have shown that dietary patterns confer protection from certain chronic diseases related to oxidative stress, the immune system and chronic low-grade inflammatory diseases. The aim of this study was to evaluate the anti-inflammatory potential and the capacity to attenuate cartilage degradation using extra-virgin olive oil–derived polyphenols for the treatment of osteoarthritis. Results show that both nutraceuticals ligstroside aglycone and acetylated ligstroside aglycone showed an anti-inflammatory profile. Acetylated ligstroside aglycone significantly reduced the expression of pro-inflammatory genes including NOS2 and MMP13 at both RNA and protein levels; decreased nitric oxide release; and, importantly, reduced proteoglycan loss in human osteoarthritis cartilage explants. Our study demonstrated that a new synthetic acetylated ligstroside aglycone derivative offers enhanced anti-inflammatory profile than the natural nutraceutical compound in osteoarthritis. These results substantiate the role of nutraceuticals in osteoarthritis with implications for therapeutic intervention and our understanding of osteoarthritis pathophysiology.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Comprehensive Analysis of the Genetic and Epigenetic Mechanisms of Osteoporosis and Bone Mineral Density

Osteoporosis is a skeletal disorder characterized by a systemic impairment of bone mineral density (BMD). Genome-wide association studies (GWAS) have identified hundreds of susceptibility loci for osteoporosis and BMD. However, the vast majority of susceptibility loci are located in non-coding regions of the genome and provide limited information about the genetic mechanisms of osteoporosis. Herein we performed a comprehensive functional analysis to investigate the genetic and epigenetic mechanisms of osteoporosis and BMD. BMD and osteoporosis are found to share many common susceptibility loci, and the corresponding susceptibility genes are significantly enriched in bone-related biological pathways. The regulatory element enrichment analysis indicated that BMD and osteoporosis susceptibility loci are significantly enriched in 5′UTR and DNase I hypersensitive sites (DHSs) of peripheral blood immune cells. By integrating GWAS and expression Quantitative Trait Locus (eQTL) data, we found that 15 protein-coding genes are regulated by the osteoporosis and BMD susceptibility loci. Our analysis provides new clues for a better understanding of the pathogenic mechanisms and offers potential therapeutic targets for osteoporosis.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

Role of Signal Transduction Pathways and Transcription Factors in Cartilage and Joint Diseases

Osteoarthritis and rheumatoid arthritis are common cartilage and joint diseases that globally affect more than 200 million and 20 million people, respectively. Several transcription factors have been implicated in the onset and progression of osteoarthritis, including Runx2, C/EBPβ, HIF2α, Sox4, and Sox11. Interleukin-1 β (IL-1β) leads to osteoarthritis through NF-ĸB, IκBζ, and the Zn²⁺-ZIP8-MTF1 axis. IL-1, IL-6, and tumor necrosis factor α (TNFα) play a major pathological role in rheumatoid arthritis through NF-ĸB and JAK/STAT pathways. Indeed, inhibitory reagents for IL-1, IL-6, and TNFα provide clinical benefits for rheumatoid arthritis patients. Several growth factors, such as bone morphogenetic protein (BMP), fibroblast growth factor (FGF), parathyroid hormone-related protein (PTHrP), and Indian hedgehog, play roles in regulating chondrocyte proliferation and differentiation. Disruption and excess of these signaling pathways cause genetic disorders in cartilage and skeletal tissues. Fibrodysplasia ossificans progressive, an autosomal genetic disorder characterized by ectopic ossification, is induced by mutant ACVR1. Mechanistic target of rapamycin kinase (mTOR) inhibitors can prevent ectopic ossification induced by ACVR1 mutations. C-type natriuretic peptide is currently the most promising therapy for achondroplasia and related autosomal genetic diseases that manifest severe dwarfism. In these ways, investigation of cartilage and chondrocyte diseases at molecular and cellular levels has enlightened the development of effective therapies. Thus, identification of signaling pathways and transcription factors implicated in these diseases is important.

Runx2 plays a central role in Osteoarthritis development

Osteoarthritis (OA) is the most common form of arthritis, is the leading cause of impaired mobility in the elderly, and accounts for more than a third of chronic moderate to severe pain. As a degenerative joint disorder, OA affects the whole joint and results in synovial hyperplasia, degradation of articular cartilage, subchondral sclerosis, osteophyte formation, and chronic pain. Currently, there is no effective drug to decelerate OA progression and molecular targets for drug development have been insufficiently investigated. Anti-OA drug development can benefit from more and precise knowledge of molecular targets for drug development. Runt-related transcription factor 2 (Runx2) is a key transcription factor controlling osteoblast and chondrocyte differentiation and is among the most promising potential therapeutic targets. Notably, Runx2 expression is upregulated in several murine OA models, suggesting a role in disease pathogenesis. In this review article, we summarized recent findings on Runx2 related to OA development and evaluated its potential as a therapeutic target.

The translational potential of this article

A better understanding of the role of Runx2 in osteoarthritis pathogenesis will contribute to the development of novel intervention of osteoarthritis disease.

Discovery and analysis of methylation quantitative trait loci (mQTLs) mapping to novel osteoarthritis genetic risk signals

OBJECTIVE

Osteoarthritis (OA) is polygenic with over 90 independent genome-wide association loci so far reported. A key next step is the identification of target genes and the molecular mechanisms through which this genetic risk operates. The majority of OA risk-conferring alleles are predicted to act by modulating gene expression. DNA methylation at CpG dinucleotides may be a functional conduit through which this occurs and is detectable by mapping methylation quantitative trait loci, or mQTLs. This approach can therefore provide functional insight into OA risk and will prioritize genes for subsequent investigation. That was our goal, with a focus on the largest set of OA loci yet to be reported.

METHOD

We investigated DNA methylation, genotype and RNA sequencing data derived from the cartilage of patients who had undergone arthroplasty and combined this with in silico analyses of expression quantitative trait loci, epigenomes and chromatin interactions.

RESULTS

We investigated 42 OA risk loci and in ten of these we identified 24 CpGs in which methylation correlated with genotype (false discovery rate (FDR) P-values ranging from 0.049 to 1.73x10^-25). In silico analyses of these mQTLs prioritised genes and regulatory elements at the majority of the ten loci, with COLGALT2 (encoding a collagen galactosyltransferase), COL11A2 (encoding a polypeptide chain of type XI collagen) and WWP2 (encoding a ubiquitin ligase active during chondrogenesis) emerging as particularly compelling target genes.

CONCLUSION

We have highlighted the pivotal role of DNA methylation as a link between genetic risk and OA and prioritized genes for further investigation.

Effect of Saliva Collection Methods on the Detection of Periodontium-Related Genetic and Epigenetic Biomarkers-A Pilot Study

Different collection methods may influence the ability to detect and quantify biomarker levels in saliva, particularly in the expression of DNA/RNA methylation regulators of several inflammations and tissue turnover markers. This pilot study recruited five participants and unstimulated saliva were collected by either spitting or drooling, and the relative preference for each method was evaluated using a visual analogue scale. Subsequently, total RNA, gDNA and proteins were isolated using the Trizol method. Thereafter, a systematic evaluation was carried out on the potential effects of different saliva collection methods on periodontium-associated genes, DNA/RNA epigenetic factors and periodontium-related DNA methylation levels. The quantity and quality of DNA and RNA were comparable from different collection methods. Periodontium-related genes, DNA/RNA methylation epigenetic factors and periodontium-associated DNA methylation could be detected in the saliva sample, with a similar expression for both methods. The methylation of tumour necrosis factor-alpha gene promoter from drooling method showed a significant positive correlation (TNF α, r = 0.9) with clinical parameter (bleeding on probing-BOP). In conclusion, the method of saliva collection has a minimal impact on detecting periodontium-related genetic and epigenetic regulators in saliva. The pilot data shows that TNF α methylation may be correlated with clinical parameters.

Expression analysis of the osteoarthritis genetic susceptibility mapping to the matrix Gla protein gene MGP

BACKGROUND

Osteoarthritis (OA) is a common disease of older individuals that impacts detrimentally on the quality and the length of life. It is characterised by the painful loss of articular cartilage and is polygenic and multifactorial. Genome-wide association scans have highlighted over 90 osteoarthritis genetic signals, some of which reside within or close to highly plausible candidate genes. An example is an association to polymorphisms within and adjacent to the matrix Gla protein gene MGP. We set out to undertake a functional study of this gene.

METHODS

Nucleic acid was extracted from cartilage, infrapatellar fat pad, synovium, trabecular bone, trapezium and peripheral whole blood from OA patients and also from mesenchymal stem cells (MSCs) subjected to chondrogenesis. Expression of MGP was measured by quantitative PCR (qPCR), RNA-sequencing and allelic expression imbalance (AEI) analysis. Matrix Gla protein was depleted in chondrocytes by knocking down MGP expression using RNA interference (RNAi) and the effect on a range of genes assessed by qPCR.

RESULTS

MGP is expressed in joint tissues, blood and chondrocytes cultured from MSCs. There is a higher expression in diseased versus non-diseased cartilage. Polymorphisms that are associated with OA also correlate with the expression of MGP, with the OA risk-conferring allele showing significantly reduced expression in cartilage, fat pad and synovium but increased expression in blood. Depletion of Matrix Gla protein had a significant effect on the majority of genes tested, with an increased expression of catabolic genes that encode enzymes that degrade cartilage.

CONCLUSIONS

MGP expression is subject to cis-acting regulators that correlate with the OA association signal. These are active in a range of joint tissues but have effects which are particularly strong in cartilage. An opposite effect is observed in blood, highlighting the context-specific nature of the regulation of this gene's expression. Recapitulation of the genetic deficit in cartilage chondrocytes is pro-catabolic.

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.