Transcription factors (TFs) regulate the expression of target genes, inducing changes in cell morphology or activities needed for cell fate determination and differentiation. The BMP signaling pathway is widely regarded as one of the most important pathways in vertebrate skeletal biology, of which BMP2 is a potent inducer, governing the osteoblast differentiation of bone marrow stromal cells (BMSCs). However, the mechanism by which BMP2 initiates its downstream transcription factor cascade and determines the direction of differentiation remains largely unknown. In this study, we used RNA-seq, ATAC-seq, and animal models to characterize the BMP2-dependent gene regulatory network governing osteoblast lineage commitment. Sp7-Cre; Bmp2fx/fx mice (BMP2-cKO) were generated and exhibited decreased bone density and lower osteoblast number (n > 6). In vitro experiments showed that BMP2-cKO mouse bone marrow stromal cells (mBMSCs) had an impact on osteoblast differentiation and deficient cell proliferation. Osteogenic medium induced mBMSCs from BMP2-cKO mice and control were subjected to RNA-seq and ATAC-seq analysis to reveal differentially expressed TFs, along with their target open chromatin regions. Combined with H3K27Ac CUT&Tag during osteoblast differentiation, we identified 2338 BMP2-dependent osteoblast-specific active enhancers. Motif enrichment assay revealed that over 80% of these elements were directly targeted by RUNX2, DLX5, MEF2C, OASIS, and KLF4. We deactivated Klf4 in the Sp7 + lineage to validate the role of KLF4 in osteoblast differentiation of mBMSCs. Compared to the wild-type, Sp7-Cre; Klf4fx/+ mice (KLF4-Het) were smaller in size and had abnormal incisors resembling BMP2-cKO mice. Additionally, KLF4-Het mice had fewer osteoblasts and decreased osteogenic ability. RNA-seq and ATAC-seq revealed that KLF4 mainly “co-bound” with RUNX2 to regulate downstream genes. Given the significant overlap between KLF4- and BMP2-dependent NFRs and enriched motifs, our findings outline a comprehensive BMP2-dependent gene regulatory network specifically governing osteoblast differentiation of the Sp7 + lineage, in which Klf4 is a novel transcription factor.

Cellular differentiation is caused by highly controlled modifications in the gene expression but rarely involves a change in the DNA sequence itself. Histone acetylation is a major epigenetic factor that adds an acetyl group to histone proteins, thus altering their interaction with DNA and nuclear proteins.Illumination of the histone acetylation during dentinogenesis is important for odontoblast differentiation and dentinogenesis. In the current study, we aimed to discover the roles and regulation of acetylation at histone 3 lysine 9 (H3K9ac) and H3K27ac during dentinogenesis. We first found that both of these modifications were enhanced during odontoblast differentiation and dentinogenesis. These modifications are dynamically catalyzed by histone acetyltransferases (HATs) and deacetylases (HDACs), among which HDAC3 was decreased while p300 increased during odontoblast differentiation. Moreover, overexpression of HDAC3 or knockdown p300 inhibited odontoblast differentiation in vitro, and inhibition of HDAC3 and p300 with trichostatin A or C646 regulated odontoblast differentiation. Taken together, the results of our present study suggest that histone acetylation is involved in dentinogenesis and coordinated expression of p300-and HDAC3-regulated odontoblast differentiation through upregulating histone acetylation.

Runt‐related transcription factor 2 (Runx2) has been shown to regulate osteoblast differentiation by directly or indirectly regulating numerous osteoblast‐related genes. However, our understanding of the transcriptional mechanisms of RUNX2 is mainly restricted to its transactivation, while the mechanism underlying its inhibitory effect during osteoblast differentiation remains largely unknown. Here, we incorporated the anti‐RUNX2 chromatin immunoprecipitation (ChIP) sequencing in MC3T3‐E1 cells and RNA‐sequencing of parietal bone from Runx2 heterozygous mutant mice, to identify the putative genes negatively regulated by RUNX2. We identified HtrA serine peptidase 1 (Htra1) as a target gene and found ten candidate Htra1 enhancers potentially regulated by RUNX2, among which seven were verified by dual‐luciferase assays. Furthermore, we investigated the motifs in the vicinity of RUNX2‐binding sites and identified early growth response 1 (EGR1) as a potential partner transcription factor (TF) potentially regulating Htra1 expression, which was subsequently confirmed by Re‐ChIP assays. RUNX2 and EGR1 co‐repressed Htra1 and increased the expression levels of other osteoblast marker genes, such as osterix, osteocalcin, and osteoprotegerin at the messenger RNA and protein level. Moreover, Alizarin red staining combined with alkaline phosphatase (ALP) staining showed decreased calcified nodules and ALP activity in the siRUNX2+siEGR1 group compared with siRUNX2 group. Our findings revealed the detailed mechanism of the inhibitory function of RUNX2 towards its downstream genes, along with its partner TFs, to promote osteoblast differentiation.

A comprehensive study of odontoblastic differentiation is essential to understand the process of tooth development and to achieve the ability of tooth regeneration in the future. Zinc finger E-box-binding homeobox 1 ( Zeb1) is a transcription factor expressed in various neural crest–derived tissues, including the mesenchyme of the tooth germ. However, its role in odontoblastic differentiation remains unknown. In this study, we found the expression of Zeb1 gradually increased during odontoblast differentiation in vivo, as well as during induced differentiation of cultured primary murine dental papilla cells (mDPCs) in vitro. In addition, the differentiation of mDPCs was repressed in Zeb1-silenced cells. We used RNA sequencing (RNA -seq) to identify the transcriptome-wide targets of Zeb1 and used assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to explore the direct targets of Zeb1 in both the early stage (embryonic day 16.5; E16.5) and the late stage (postnatal day 0; PN0) of tooth development. We identified the motifs of transcription factors enriched in Zeb1-dependent accessible chromatin regions and observed that only in the early stage of mDPCs could Zeb1 significantly change the accessibility of chromatin regions. In vivo and in vitro experiments confirmed that silencing of Zeb1 at E16.5 inhibited dentinogenesis. Analysis of RNA-seq and ATAC-seq resulted in the identification of Runx2, a gene directly regulated by Zeb1 during early odontoblast differentiation. Zeb1 enhances the expression of Runx2 by binding to its cis-elements, and ZEB1 interacts with RUNX2. In the late stage of tooth development, we found that ZEB1 could directly bind to and increase the enhancer activity of an element upstream of Dspp and promote dentinogenesis. In this study, for the first time, we revealed that ZEB1 promoted odontoblast differentiation in the early stage by altering chromatin accessibility of cis-elements near genes such as Runx2, while in the late stage, it directly enhanced Dspp transcription, thereby performing a dual role.

Embryonic development and stem cell differentiation are orchestrated by changes in sequential binding of regulatory transcriptional factors to their motifs. These processes are invariably accompanied by the alternations in chromatin accessibility, conformation, and histone modification. Odontoblast lineage originates from cranial neural crest cells and is crucial in dentinogenesis. Our previous work revealed several transcription factors (TFs) that promote odontoblast differentiation. However, it remains elusive as to whether chromatin accessibility affects odontoblast terminal differentiation. Herein, integration of single-cell RNA-seq and bulk RNA-seq revealed that in vitro odontoblast differentiation using dental papilla cells at E18.5 was comparable to the crown odontoblast differentiation trajectory of OC (osteocalcin)-positive odontogenic lineage. Before in vitro odontoblast differentiation, ATAC-seq and H3K27Ac CUT and Tag experiments demonstrated high accessibility of chromatin regions adjacent to genes associated with odontogenic potential. However, following odontoblastic induction, regions near mineralization-related genes became accessible. Integration of RNA-seq and ATAC-seq results further revealed that the expression levels of these genes were correlated with the accessibility of nearby chromatin. Time-course ATAC-seq experiments further demonstrated that odontoblast terminal differentiation was correlated with the occupation of the basic region/leucine zipper motif (bZIP) TF family, whereby we validated the positive role of ATF5 in vitro. Collectively, this study reports a global mapping of open chromatin regulatory elements during dentinogenesis and illustrates how these regions are regulated via dynamic binding of different TF families, resulting in odontoblast terminal differentiation. The findings also shed light on understanding the genetic regulation of dentin regeneration using dental mesenchymal stem cells.