Complex Regulatory Interactions at GDF5 Shape Joint Morphology and Osteoarthritis Disease Risk

OBJECTIVE: The objective of this study was to reveal causal-level osteoarthritis (OA) disease biology by targeting regulatory interactions at GDF5. METHODS: By investigating different GDF5 regulatory regions (R2, R3-R5, R7-R9, R18-R20, GROW1), we explored their functional impacts on gene expression and joint morphology in vivo and in vitro. We additionally modeled OA variants in said enhancers in in vitro and in vivo mouse models for expression and disease effects. RESULTS: For all regulatory regions, we found evidence of activation and repression between or within said regions that impacted patterns of joint-specific expression. Examples are as follows: (1) the R4 enhancer, although considered to be activating, has dual roles repressing expression in adjacent tissues and sites, and (2) growth plate-specific expression patterns by the GROW1 regulatory region are confined by adjacent sequences to restrict its expression to the perichondrium. We next targeted different regions and variants in vivo. Testing the R2de region resulted in ~40% reduction in Gdf5 expression and joint morphology changes but no increase in OA risk; likewise, modeling the most cited OA risk variant (rs143384) in mice had no impact on expression, joint morphology, or disease. However, we identified epistatic interactions between this rs143384 risk variant and downstream disease risk variants lying within regulatory regions subject to repression, which compound to impact expression. CONCLUSION: These findings, at the best studied OA locus to date, serve as lessons on the nature of how gene regulatory interactions and local epistasis work in the etiology of OA disease risk, and that assessment of individual variants of high genome-wide association study significance need not alone be considered causal.