Abstract
Objective
Modic changes (MC), a form of intervertebral disc degeneration visible as subchondral and vertebral bone marrow changes on spine magnetic resonance (MR), are known to be associated with low back pain. This study aimed to identify genes contributing to the development of MC using genome-wide association study.
Methods
Presence of MC was evaluated in lumbar MR images in the Northern Finland Birth Cohort 1966 (NFBC1966, N=1182) and TwinsUK (N=647). Genome-wide association analyses were carried out in the cohorts separately using a linear regression model fitted to test for additive effects of SNPs and adjusting for age, sex, BMI, and either family relatedness via a kinship matrix (TwinsUK) or population stratification using principal components (NFBC1966). Meta-analysis of the two studies was carried out using the inverse-variance weighting approach.
Results
A locus associated with MC reaching genome-wide significance (p<5e-8) was found on chromosome 9 with the lead SNP rs1934268 in intron 6 of the PTPRD gene. The SNP is located in the region of binding for a number of transcription factors which are involved in the development of the musculoskeletal system and spine cord.
Conclusions
The first GWAS of MC has identified a likely functional intronic locus in PTPRD on chromosome 9 implicating musculoskeletal development. This work sheds light on the genesis of MC and paves the way for further studies on the shared genetic factors underlying the various features of spine degeneration.
No conflicts of interest
Sources of Funding: The study was supported by EU FP7 project PainOMICs (grant agreement #602736), University of Oulu (grant #24000692), Oulu University Hospital (grant #24301140), and the European Regional Development Fund (grant # 539/2010 A31592). MBF, MK, and SS contributed equally to this study.
