Epigenetic mechanisms of osteoarthritis risk in human skeletal development.

The epigenome, including the methylation of cytosine bases at CG dinucleotides, is intrinsically linked to transcriptional regulation. The tight regulation of gene expression during skeletal development is essential, with ~1/500 individuals born with skeletal abnormalities. Furthermore, increasing evidence is emerging to link age-associated complex genetic musculoskeletal diseases, including osteoarthritis (OA), to developmental factors including joint shape. Multiple studies have shown a functional role for DNA methylation in the genetic mechanisms of OA risk using articular cartilage samples taken from aged patients. Despite this, our knowledge of temporal changes to the methylome during human cartilage development has been limited. We quantified DNA methylation at ~700,000 individual CpGs across the epigenome of developing human articular cartilage in 72 samples ranging from 7-21 post-conception weeks, a time period that includes cavitation of the developing knee joint. We identified significant changes in 8% of all CpGs, and >9400 developmental differentially methylated regions (dDMRs). The largest hypermethylated dDMRs mapped to transcriptional regulators of early skeletal patterning including MEIS1 and IRX1. Conversely, the largest hypomethylated dDMRs mapped to genes encoding extracellular matrix proteins including SPON2 and TNXB and were enriched in chondrocyte enhancers. Significant correlations were identified between the expression of these genes and methylation within the hypomethylated dDMRs. We further identified 811 CpGs at which significant dimorphism was present between the male and female samples, with the majority (68%) being hypermethylated in female samples. Following imputation, we captured the genotype of these samples at >5 million variants and performed epigenome-wide methylation quantitative trait locus (mQTL) analysis. Colocalization analysis identified 26 loci at which genetic variants exhibited shared impacts upon methylation and OA genetic risk. This included loci which have been previously reported to harbour OA-mQTLs (including GDF5 and ALDH1A2), yet the majority (73%) were novel (including those mapping to CHST3, FGF1 and TEAD1). To our knowledge, this is the first extensive study of DNA methylation across human articular cartilage development. We identify considerable methylomic plasticity within the development of knee cartilage and report active epigenomic mediators of OA risk operating in prenatal joint tissues.

Genetic risk of osteoarthritis operates during human skeletogenesis.

Osteoarthritis (OA) is a polygenic disease of older people resulting in the breakdown of cartilage within articular joints. Although it is a leading cause of disability, there are no disease-modifying therapies. Evidence is emerging to support the origins of OA in skeletogenesis. Whereas methylation quantitative trait loci (mQTLs) co-localizing with OA genome-wide association study signals have been identified in aged human cartilage and used to identify effector genes and variants, such analyses have never been conducted during human development. Here, for the first time, we have investigated the developmental origins of OA genetic risk at seven well-characterized OA risk loci, comprising 39 OA-mQTL CpGs, in human fetal limb (FL) and cartilage (FC) tissues using a range of molecular genetic techniques. We identified significant OA-mQTLs at 14 and 29 CpGs in FL and FC tissues, respectively, and compared our results with aged cartilage samples (AC). Differential methylation was observed at 26 sites between FC and AC, with the majority becoming actively hypermethylated in old age. Notably, 6/9 OA effector genes showed allelic expression imbalances during fetal development. Finally, we conducted ATAC-sequencing in cartilage from the developing and aged hip and knee to identify accessible chromatin regions and found enrichment for transcription factor binding motifs including SOX9 and FOS/JUN. For the first time, we have demonstrated the activity of OA-mQTLs and expression imbalance of OA effector genes during human skeletogenesis. We show striking differences in the spatiotemporal function of these loci, contributing to our understanding of OA aetiology, with implications for the timing and strategy of pharmacological interventions.

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.

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.