Super enhancers targeting ZBTB16 in osteogenesis protect against osteoporosis

Abstract

As the major cell precursors in osteogenesis, mesenchymal stem cells (MSCs) are indispensable for bone homeostasis and development. However, the primary mechanisms regulating osteogenic differentiation are controversial. Composed of multiple constituent enhancers, super enhancers (SEs) are powerful cis-regulatory elements that identify genes that ensure sequential differentiation. The present study demonstrated that SEs were indispensable for MSC osteogenesis and involved in osteoporosis development. Through integrated analysis, we identified the most common SE-targeted and osteoporosis-related osteogenic gene, ZBTB16. ZBTB16, positively regulated by SEs, promoted MSC osteogenesis but was expressed at lower levels in osteoporosis. Mechanistically, SEs recruited bromodomain containing 4 (BRD4) at the site of ZBTB16, which then bound to RNA polymerase II-associated protein 2 (RPAP2) that transported RNA polymerase II (POL II) into the nucleus. The subsequent synergistic regulation of POL II carboxyterminal domain (CTD) phosphorylation by BRD4 and RPAP2 initiated ZBTB16 transcriptional elongation, which facilitated MSC osteogenesis via the key osteogenic transcription factor SP7. Bone-targeting ZBTB16 overexpression had a therapeutic effect on the decreased bone density and remodeling capacity of Brd4 ^fl/fl Prx1-cre mice and osteoporosis (OP) models. Therefore, our study shows that SEs orchestrate the osteogenesis of MSCs by targeting ZBTB16 expression, which provides an attractive focus and therapeutic target for osteoporosis.

Without SEs located on osteogenic genes, BRD4 is not able to bind to osteogenic identity genes due to its closed structure before osteogenesis. During osteogenesis, histones on osteogenic identity genes are acetylated, and OB-gain SEs appear, enabling the binding of BRD4 to the osteogenic identity gene ZBTB16. RPAP2 transports RNA Pol II from the cytoplasm to the nucleus and guides Pol II to target ZBTB16 via recognition of the navigator BRD4 on SEs. After the binding of the RPAP2-Pol II complex to BRD4 on SEs, RPAP2 dephosphorylates Ser5 at the Pol II CTD to terminate the transcriptional pause, and BRD4 phosphorylates Ser2 at the Pol II CTD to initiate transcriptional elongation, which synergistically drives efficient transcription of ZBTB16, ensuring proper osteogenesis. Dysregulation of SE-mediated ZBTB16 expression leads to osteoporosis, and bone-targeting ZBTB16 overexpression is efficient in accelerating bone repair and treating osteoporosis.

[graphic not available: see fulltext]

References

[1] Salhotra A, Shah HN, Levi B, Longaker MT. Mechanisms of bone development and repair. Nat. Rev. Mol. Cell Biol., 2020, 21: 696-711

[2] Wang R, Wang Y, Zhu L, Liu Y, Li W. Epigenetic regulation in mesenchymal stem cell aging and differentiation and osteoporosis. Stem Cells Int., 2020, 2020: 8836258

[3] Whyte WA et al. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell, 2013, 153: 307-319

[4] Brown JD et al. BET bromodomain proteins regulate enhancer function during adipogenesis. Proc. Natl Acad. Sci. USA, 2018, 115: 2144-2149

[5] Zhao Y et al. MyoD induced enhancer RNA interacts with hnRNPL to activate target gene transcription during myogenic differentiation. Nat. Commun., 2019, 10

[6] Chen Z et al. Fusion between a novel Krüppel-like zinc finger gene and the retinoic acid receptor-alpha locus due to a variant t(11;17) translocation associated with acute promyelocytic leukaemia. Embo J., 1993, 12: 1161-1167

[7] Vincent-Fabert C et al. PLZF mutation alters mouse hematopoietic stem cell function and cell cycle progression. Blood, 2016, 127: 1881-1885

[8] Hosokawa H et al. Bcl11b sets pro-T cell fate by site-specific cofactor recruitment and by repressing Id2 and Zbtb16. Nat. Immunol., 2018, 19: 1427-1440

[9] Wasim M et al. PLZF/ZBTB16, a glucocorticoid response gene in acute lymphoblastic leukemia, interferes with glucocorticoid-induced apoptosis. J. Steroid Biochem. Mol. Biol., 2010, 120: 218-227

[10] Sharma, M. et al. Identification of EOMES-expressing spermatogonial stem cells and their regulation by PLZF. Elife 8, e43352 (2019).

[11] Barna M, Hawe N, Niswander L, Pandolfi PP. Plzf regulates limb and axial skeletal patterning. Nat. Genet., 2000, 25: 166-172

[12] Onizuka S et al. ZBTB16 as a downstream target gene of osterix regulates osteoblastogenesis of human multipotent mesenchymal stromal cells. J. Cell Biochem., 2016, 117: 2423-2434

[13] Felthaus O, Gosau M, Morsczeck C. ZBTB16 induces osteogenic differentiation marker genes in dental follicle cells independent from RUNX2. J. Periodontol., 2014, 85: e144-151

[14] Rauch A et al. Osteogenesis depends on commissioning of a network of stem cell transcription factors that act as repressors of adipogenesis. Nat. Genet., 2019, 51: 716-727

[15] Najafova Z et al. BRD4 localization to lineage-specific enhancers is associated with a distinct transcription factor repertoire. Nucleic Acids Res., 2017, 45: 127-141

[16] Sabari, B. R. et al. Coactivator condensation at super-enhancers links phase separation and gene control. Science 361, eaar3958 (2018).

[17] Loven J et al. Selective inhibition of tumor oncogenes by disruption of super-enhancers. Cell, 2013, 153: 320-334

[18] Alghamdi S et al. BET protein inhibitor JQ1 inhibits growth and modulates WNT signaling in mesenchymal stem cells. Stem Cell Res. Ther., 2016, 7: 22

[19] Geng Y et al. Systematic analysis of mRNAs and ncRNAs in BMSCs of senile osteoporosis patients. Front. Genet., 2021, 12: 776984

[20] Forget D et al. Nuclear import of RNA polymerase II is coupled with nucleocytoplasmic shuttling of the RNA polymerase II-associated protein 2. Nucleic Acids Res., 2013, 41: 6881-6891

[21] Egloff S, Murphy S. Cracking the RNA polymerase II CTD code. Trends Genet., 2008, 24: 280-288

[22] Ni Z et al. RPRD1A and RPRD1B are human RNA polymerase II C-terminal domain scaffolds for Ser5 dephosphorylation. Nat. Struct. Mol. Biol., 2014, 21: 686-695

[23] Devaiah BN et al. BRD4 is an atypical kinase that phosphorylates serine2 of the RNA polymerase II carboxy-terminal domain. Proc. Natl Acad. Sci. USA, 2012, 109: 6927-6932

[24] Yang YS et al. Bone-targeting AAV-mediated silencing of Schnurri-3 prevents bone loss in osteoporosis. Nat. Commun., 2019, 10

[25] Gao J et al. SIRT3/SOD2 maintains osteoblast differentiation and bone formation by regulating mitochondrial stress. Cell Death Differ., 2018, 25: 229-240

[26] Pal S, Porwal K, Rajak S, Sinha RA, Chattopadhyay N. Selective dietary polyphenols induce differentiation of human osteoblasts by adiponectin receptor 1-mediated reprogramming of mitochondrial energy metabolism. Biomed. Pharmacother., 2020, 127: 110207

[27] Chen X et al. Regulatory role of RNA N(6)-methyladenosine modification in bone biology and osteoporosis. Front. Endocrinol., 2019, 10: 911

[28] Zhang W et al. Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells. Stem Cell Res. Ther., 2017, 8: 30

[29] Liu, Z. et al. Myeloma cells shift osteoblastogenesis to adipogenesis by inhibiting the ubiquitin ligase MURF1 in mesenchymal stem cells. Sci. Signal 13, eaay8203 (2020).

[30] Pott S, Lieb JD. What are super-enhancers? Nat. Genet., 2015, 47: 8-12

[31] Siersbaek R et al. Transcription factor cooperativity in early adipogenic hotspots and super-enhancers. Cell Rep., 2014, 7: 1443-1455

[32] Lee BK et al. Super-enhancer-guided mapping of regulatory networks controlling mouse trophoblast stem cells. Nat. Commun., 2019, 10

[33] Paradise CR et al. The epigenetic reader Brd4 is required for osteoblast differentiation. J. Cell Physiol., 2020, 235: 5293-5304

[34] Paradise CR et al. Brd4 is required for chondrocyte differentiation and endochondral ossification. Bone, 2022, 154: 116234

[35] Lin L et al. Super-enhancer-associated MEIS1 promotes transcriptional dysregulation in Ewing sarcoma in co-operation with EWS-FLI1. Nucleic Acids Res., 2019, 47: 1255-1267

[36] Shin HY et al. Hierarchy within the mammary STAT5-driven Wap super-enhancer. Nat. Genet., 2016, 48: 904-911

[37] Marofi F et al. Gene expression of TWIST1 and ZBTB16 is regulated by methylation modifications during the osteoblastic differentiation of mesenchymal stem cells. J. Cell Physiol., 2019, 234: 6230-6243

[38] Hall DD, Spitler KM, Grueter CE. Disruption of cardiac Med1 inhibits RNA polymerase II promoter occupancy and promotes chromatin remodeling. Am. J. Physiol. Heart Circ. Physiol., 2019, 316: H314-h325

[39] Chen FX, Smith ER, Shilatifard A. Born to run: control of transcription elongation by RNA polymerase II. Nat. Rev. Mol. Cell Biol., 2018, 19: 464-478

[40] Adelman K, Lis JT. Promoter-proximal pausing of RNA polymerase II: emerging roles in metazoans. Nat. Rev. Genet., 2012, 13: 720-731

[41] Harlen KM, Churchman LS. The code and beyond: transcription regulation by the RNA polymerase II carboxy-terminal domain. Nat. Rev. Mol. Cell Biol., 2017, 18: 263-273

[42] Egloff S, Zaborowska J, Laitem C, Kiss T, Murphy S. Ser7 phosphorylation of the CTD recruits the RPAP2 Ser5 phosphatase to snRNA genes. Mol. Cell, 2012, 45: 111-122

[43] Reid IR, Billington EO. Drug therapy for osteoporosis in older adults. Lancet, 2022, 399: 1080-1092

[44] Cheng C, Wentworth K, Shoback DM. New frontiers in osteoporosis therapy. Annu. Rev. Med., 2020, 71: 277-288

[45] Ma, Y. et al. Autophagy controls mesenchymal stem cell properties and senescence during bone aging. Aging cell 17, e12709 (2018).

[46] Guo Y et al. Sirt3-mediated mitophagy regulates AGEs-induced BMSCs senescence and senile osteoporosis. Redox Biol., 2021, 41: 101915

[47] Colella P, Ronzitti G, Mingozzi F. Emerging Issues in AAV-Mediated In Vivo Gene Therapy. Mol. Ther. Methods Clin. Dev., 2018, 8: 87-104

[48] Hnisz D et al. Super-enhancers in the control of cell identity and disease. Cell, 2013, 155: 934-947

[49] Heinz S et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell, 2010, 38: 576-589

[50] Li R, Li Y, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics, 2008, 24: 713-714

[51] Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods, 2015, 12: 357-360

[52] Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat. Methods, 2012, 9: 357-359

[53] Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform., 2011, 12: 323

[54] Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol., 2014, 15

[55] Zhang Y et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol., 2008, 9

[56] Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res., 2014, 42: W187-191

Funding

National Natural Science Foundation of China (National Science Foundation of China)(82172385)

The Key-Area Research and Development Program of Guangdong Province [2019B020236001] Shenzhen Key Medical Discipline Construction Fund [ZDSYS20190902092851024]

The Natural Science Foundation of Guangdong Province [2020A1515010097]. Shenzhen Outstanding Science and Technology Innovation Talents - Outstanding Youth Fund project [RCYX20210706092106042]