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.