Enhanced osteogenic potential of iPSC-derived mesenchymal progenitor cells following genome editing of GWAS variants in the RUNX1 gene

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.

References

[1.] CefaluCA. Is bone mineral density predictive of fracture risk reduction?. Curr. Med. Res. Opin., 2004, 20: 341-349

[2.] ClarkEM, NessAR, BishopNJ, TobiasJH. Association between bone mass and fractures in children: a prospective cohort study. J. Bone Miner. Res. Off. J. Am. Soc. Bone Miner. Res., 2006, 21: 1489-1495

[3.] GouldingA, JonesIE, TaylorRW, WilliamsSM, ManningPJ. Bone mineral density and body composition in boys with distal forearm fractures: a dual-energy x-ray absorptiometry study. J. Pediatr., 2001, 139: 509-515

[4.] KalkwarfHJ, LaorT, BeanJA. Fracture risk in children with a forearm injury is associated with volumetric bone density and cortical area (by peripheral QCT) and areal bone density (by DXA). Osteoporos. Int., 2011, 22: 607-616

[5.] HeaneyRP, et al. . Peak bone mass. Osteoporos. Int., 2000, 11: 985-1009

[6.] MorrisJA, et al. . An atlas of genetic influences on osteoporosis in humans and mice. Nat. Genet., 2019, 51: 258-266

[7.] YergesLM, et al. . Candidate gene analysis of femoral neck trabecular and cortical volumetric bone mineral density in older men. J. Bone Miner. Res. Off. J. Am. Soc. Bone Miner. Res., 2010, 25: 330-338

[8.] ZhengHF, et al. . Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture. Nature, 2015, 526: 112-117

[9.] Díaz-HernándezME, et al. . Sexually dimorphic increases in bone mass following tissue-specific overexpression of Runx1 in Osteoclast Precursors. Endocrinology, 2022, 163: bqac113

[10.] PagliaDN, et al. . Deletion of Runx1 in osteoclasts impairs murine fracture healing through progressive woven bone loss and delayed cartilage remodeling. J. Orthop. Res. Off. Publ. Orthop. Res. Soc., 2020, 38: 1007-1015

[11.] PagliaDN, et al. . Runx1 regulates myeloid precursor differentiation into osteoclasts without affecting differentiation into antigen presenting or phagocytic cells in both males and females. Endocrinology, 2016, 157: 3058-3069

[12.] Soung doY, et al. . Runx1-mediated regulation of osteoclast differentiation and function. Mol. Endocrinol., 2014, 28: 546-553

[13.] SollisE, et al. . The NHGRI-EBI GWAS Catalog: knowledgebase and deposition resource. Nucleic acids Res., 2023, 51: 977-985

[14.] MauranoMT, et al. . Systematic localization of common disease-associated variation in regulatory DNA. Science, 2012, 337: 1190-1195

[15.] MusunuruK, et al. . From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus. Nature, 2010, 466: 714-719

[16.] SanyalA, LajoieBR, JainG, DekkerJ. The long-range interaction landscape of gene promoters. Nature, 2012, 489: 109-113

[17.] ClaussnitzerM, et al. . FTO obesity variant circuitry and adipocyte browning in humans. N. Engl. J. Med., 2015, 373: 895-907

[18.] SmemoS, et al. . Obesity-associated variants within FTO form long-range functional connections with IRX3. Nature, 2014, 507: 371-375

[19.] GeevenG, TeunissenH, de LaatW, de WitE. peakC: a flexible, non-parametric peak calling package for 4C and Capture-C data. Nucleic Acids Res., 2018, 46: e91

[20.] AutonA, et al. . A global reference for human genetic variation. Nature, 2015, 526: 68-74

[21.] KundajeA, et al. . Integrative analysis of 111 reference human epigenomes. Nature, 2015, 518: 317-330

[22.] Diaz-HernandezME, et al. . Derivation of notochordal cells from human embryonic stem cells reveals unique regulatory networks by single cell-transcriptomics. J. Cell. Physiol., 2020, 235: 5241-5255

[23.] BilousovaG, et al. . Osteoblasts derived from induced pluripotent stem cells form calcified structures in scaffolds both in vitro and in vivo. Stem cells, 2011, 29: 206-216

[24.] ZujurD, et al. . Stepwise strategy for generating osteoblasts from human pluripotent stem cells under fully defined xeno-free conditions with small-molecule inducers. Regen. Ther., 2020, 14: 19-31

[25.] GuzzoRM, DrissiH. Differentiation of human induced pluripotent stem cells to chondrocytes. Methods Mol. Biol., 2015, 1340: 79-95

[26.] Khan, N. M. et al. Differential chondrogenic differentiation between iPSC derived from healthy and OA cartilage is associated with changes in epigenetic regulation and metabolic transcriptomic signatures. eLife 12, e83138 (2023).

[27.] KhanNM, Diaz-HernandezME, DrissiH. Differentiation of human induced pluripotent stem cells (iPSCs)-derived mesenchymal progenitors into chondrocytes. Bio-Protoc., 2023, 13: e4874

[28.] LiuX, FangT, ShiT, WangY, LiuG. Hydrogels provide microenvironments to mesenchymal stem cells for craniofacial bone regeneration: review. J. Biomater. Appl., 2023, 38: 3-24

[29.] KamalakarA, et al. . JAGGED1 stimulates cranial neural crest cell osteoblast commitment pathways and bone regeneration independent of canonical NOTCH signaling. Bone, 2021, 143: 115657

[30.] GarcíaAJ. PEG-maleimide hydrogels for protein and cell delivery in regenerative medicine. Ann. Biomed. Eng., 2014, 42: 312-322

[31.] YeS, et al. . SET domain-containing protein 4 epigenetically controls breast cancer stem cell quiescence. Cancer Res., 2019, 79: 4729-4743

[32.] AgredoA, KasinskiAL. Histone 4 lysine 20 tri-methylation: a key epigenetic regulator in chromatin structure and disease. Front. Genet., 2023, 14: 1243395

[33.] LiaoX, et al. . SETD4 in the proliferation, migration, angiogenesis, myogenic differentiation and genomic methylation of bone marrow mesenchymal stem cells. Stem Cell Rev. Rep., 2021, 17: 1374-1389

[34.] DixonJR, et al. . Chromatin architecture reorganization during stem cell differentiation. Nature, 2015, 518: 331-336

[35.] NaikNG, et al. . Epigenetic factor siRNA screen during primary KSHV infection identifies novel host restriction factors for the lytic cycle of KSHV. PLoS Pathog., 2020, 16: e1008268

[36.] ImJY, et al. . Bone regeneration of mouse critical-sized calvarial defects with human mesenchymal stem cells in scaffold. Lab Anim. Res, 2013, 29: 196-203,

[37.] LeviB, et al. . Human adipose derived stromal cells heal critical size mouse calvarial defects. PLoS ONE, 2010, 5: e11177

[38.] AmbrosiTH, et al. . Aged skeletal stem cells generate an inflammatory degenerative niche. Nature, 2021, 597: 256-262

[39.] ChanCK, et al. . Identification and specification of the mouse skeletal stem cell. Cell, 2015, 160: 285-298

[40.] DebnathS, et al. . Discovery of a periosteal stem cell mediating intramembranous bone formation. Nature, 2018, 562: 133-139

[41.] KhanNM, et al. . Wogonin, a plant derived small molecule, exerts potent anti-inflammatory and chondroprotective effects through the activation of ROS/ERK/Nrf2 signaling pathways in human Osteoarthritis chondrocytes. Free Radic. Biol. Med., 2017, 106: 288-301

[42.] KhanNM, CliftonKB, LorenzoJ, HansenMF, DrissiH. Comparative transcriptomic analysis identifies distinct molecular signatures and regulatory networks of chondroclasts and osteoclasts. Arthritis Res. Ther., 2020, 22: 168

[43.] KhanNM, Diaz-HernandezME, PresciuttiSA-O, DrissiH. Network analysis identifies gene regulatory network indicating the role of RUNX1 in human intervertebral disc degeneration. Genes (Basel), 2020, 11: 771

[44.] GuzzoRM, GibsonJ, XuRH, LeeFY, DrissiH. Efficient differentiation of human iPSC-derived mesenchymal stem cells to chondroprogenitor cells. J. Cell. Biochem., 2013, 114: 480-490

[45.] DominiciM, et al. . Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy, 2006, 8: 315-317

[46.] KhanNM, PoduvalTB. Immunomodulatory and immunotoxic effects of bilirubin: molecular mechanisms. J. Leukoc. Biol., 2011, 90: 997-1015

[47.] Diaz-HernandezME, KhanNM, DrissiH. Efficient differentiation of human induced pluripotent stem cell (hiPSC)-derived mesenchymal progenitors into adipocytes and osteoblasts. Bio-Protoc., 2023, 13: e4885

[48.] KhanNM, AhmadI, HaqqiTM. Nrf2/ARE pathway attenuates oxidative and apoptotic response in human osteoarthritis chondrocytes by activating ERK1/2/ELK1-P70S6K-P90RSK signaling axis. Free Radic. Biol. Med., 2018, 116: 159-171

[49.] BouxseinML, et al. . Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J. Bone Miner. Res. Off. J. Am. Soc. Bone Miner. Res., 2010, 25: 1468-1486

[50.] KhanNM, et al. . pH-sensing G protein-coupled orphan receptor GPR68 is expressed in human cartilage and correlates with degradation of extracellular matrix during OA progression. PeerJ, 2023, 11: e16553

[51.] van de WerkenHJ, et al. . 4C technology: protocols and data analysis. Methods Enzymol., 2012, 513: 89-112

[52.] CotneyJ, et al. . Chromatin state signatures associated with tissue-specific gene expression and enhancer activity in the embryonic limb. Genome Res., 2012, 22: 1069-1080

[53.] BolgerAM, LohseM, UsadelB. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30: 2114-2120

[54.] LangmeadB, SalzbergSL. Fast gapped-read alignment with Bowtie 2. Nat. Methods, 2012, 9: 357-359

[55.] ServantN, et al. . HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol., 2015, 16 259

[56.] RamírezF, et al. . High-resolution TADs reveal DNA sequences underlying genome organization in flies. Nat. Commun., 2018, 9 189

Funding

U.S. Department of Health & Human Services | NIH | National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS)(AR071536)