Abstract
Genome-wide-association studies (GWASs) have discovered genetic signals robustly associated with BMD, but typically not the precise localization of effector genes. By intersecting genome-wide promoter-focused Capture C and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data generated in human mesenchymal progenitor cell (hMSC)-derived osteoblasts, consistent contacts were previously predicted between the EPDR1 promoter and multiple BMD-associated candidate causal variants at the 'STARD3NL' locus. RNAi knockdown of EPDR1 expression in hMSC-derived osteoblasts was shown to lead to inhibition of osteoblastogenesis. To fully characterize the physical connection between these putative noncoding causal variants at this locus and the EPDR1 gene, clustered regularly interspaced short-palindromic repeat Cas9 endonuclease (CRISPR-Cas9) genome editing was conducted in hFOB1.19 cells across the single open-chromatin region harboring candidates for the underlying causal variant, rs1524068, rs6975644, and rs940347, all in close proximity to each other. RT-qPCR and immunoblotting revealed dramatic and consistent downregulation of EPDR1 specifically in the edited differentiated osteoblast cells. Consistent with EPDR1 expression changes, alkaline phosphatase staining was also markedly reduced in the edited differentiated cells. Collectively, CRISPR-Cas9 genome editing in the hFOB1.19 cell model supports previous observations, where this regulatory region harboring GWAS-implicated variation operates through direct long-distance physical contact, further implicating a key role for EPDR1 in osteoblastogenesis and BMD determination. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.