Abstract
Recent genome-wide association studies (GWAS) identified 518 significant loci associated with bone mineral density (BMD), including variants at the RUNX1 locus (rs13046645, rs2834676, and rs2834694). However, their regulatory impact on RUNX1 expression and bone formation remained unclear. This study utilized human induced pluripotent stem cells (iPSCs) differentiated into osteoblasts to investigate these variants' regulatory roles. CRISPR/Cas9 was employed to generate mutant (Δ) iPSC lines lacking these loci at the RUNX1 locus. Deletion lines (Δ1 and Δ2) were created in iPSCs to assess the effects of removing regions containing these loci. Deletion lines exhibited enhanced osteogenic potential, with increased expression of osteogenic marker genes and Alizarin Red staining. Circularized chromosome conformation capture (4C-Seq) was utilized to analyze interactions between BMD-associated loci and the RUNX1 promoter during osteogenesis. Analysis revealed altered chromatin interactions with multiple gene promoters including RUNX1 isoform, as well as SETD4, a histone methyltransferase, indicating their regulatory influence. Interestingly, both deletion lines notably stimulated the expression of the long isoform of RUNX1, with more modest effects on the shorter isoform. Consistent upregulation of SETD4 and other predicted targets within the Δ2 deletion suggested its removal removed a regulatory hub constraining expression of multiple genes at this locus. In vivo experiments using a bone defect model in mice demonstrated increased bone regeneration with homozygous deletion of the Δ2 region. These findings indicate that BMD-associated variants within the RUNX1 locus regulate multiple effector genes involved in osteoblast commitment, providing valuable insights into genetic regulation of bone density and potential therapeutic targets.