Abstract: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 and . Conversely, the largest hypomethylated dDMRs mapped to genes encoding extracellular matrix proteins including and 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 and ), yet the majority (73%) were novel (including those mapping to and ). 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.

Primary cartilage transcriptional signatures reflect cell-type-specific molecular pathways underpinning osteoarthritis

Comprehensive profiling of transcriptional regulation in cartilage reveals pathogenesis of osteoarthritis

Transcriptome regulation analysis in osteoarthritis primary tissues identifies high-confidence effector genes

Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis

Single-cell RNA-seq analysis reveals the progression of human osteoarthritis.

Genome-wide expression profiling reveals new candidate genes associated with osteoarthritis

Integrated Analysis of Transcriptome Changes in Osteoarthritis: Gene Expression, Pathways and Alternative Splicing

Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel candidate effector genes

Insights into osteoarthritis progression revealed by analyses of both knee tibiofemoral compartments

Single-Cell RNA-Seq Reveals Transcriptomic Heterogeneity and Post-Traumatic Osteoarthritis-Associated Early Molecular Changes in Mouse Articular Chondrocytes

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk

OA30 Primary osteoarthritis chondrocyte map of chromatin conformation reveals novel effector genes

Epigenomic differences between osteoarthritis grades in primary cartilage

Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration

Knee Osteoarthritis Molecular Subgroups Revealed by Pathogenic Joint Gene Expression

Epigenetic mechanisms of osteoarthritis risk in human skeletal development

RNA Microarray Analysis of Macroscopically Normal Articular Cartilage from Knees Undergoing Partial Medial Meniscectomy: Potential Prediction of the Risk for Developing Osteoarthritis

Investigation of Candidate Genes for Osteoarthritis Based on Gene Expression Profiles.

Hyper-physiologic mechanical cues, as an osteoarthritis disease-relevant environmental perturbation, cause a critical shift in set points of methylation at transcriptionally active CpG sites in neo-cartilage organoids

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

Single-cell RNA sequencing analysis of human chondrocytes reveals cell–cell communication alterations mediated by interactive signaling pathways in osteoarthritis