Genome-wide profiling of bone reveals differentially methylated regions in osteoporosis and osteoarthritis

Objective To determine genome-wide methylation profiles of bone from patients with hip osteoarthritis (OA) and those with osteoporotic (OP) hip fractures. Methods Trabecular bone pieces were obtained from the central part of the femoral head of 27 patients with hip fractures and 26 patients with hip OA. DNA was isolated, and methylation was explored with Illumina methylation arrays. RNA was extracted, pooled, and deep-sequenced to obtain the whole transcriptome. Differentially methylated regions were identified, and connections between genes with differentially methylated regions were explored by pathway and text-mining analyses. Results After quality control, methylation of 23,367 CpG sites (13,463 genes) was analyzed. There was a genome-wide inverse relationship between methylation and gene expression in both patient groups. Comparison of OP and OA bones revealed 241 CpG sites, located in 228 genes, with significant differences in methylation (false discovery rate <0.05). Of them, 217 were less methylated in OP than in OA. The absolute methylation differences were >5% in 128 CpG sites and >10% in 45 CpG sites. The differentially methylated genes were enriched for association with bone traits in the genome-wide association study catalog. Pathway analysis and text-mining analysis with Gene Relationships Across Implicated Loci software revealed enrichment in genes participating in glycoprotein metabolism or cell differentiation, and particularly in the homeobox superfamily of transcription factors. Conclusion Genome-wide methylation profiling of bone samples revealed differentially methylated regions in OP and OA. These regions were enriched in genes associated with cell differentiation and skeletal embryogenesis, such as those in the homeobox superfamily, suggesting the existence of a developmental component in the predisposition to these disorders.