Fig. 9

Dynamic changes of RUNX2 and CREB3L1 regulon activity. A RUNX2 regulatory network. Blue dots represent skeletal system development-related target genes, other target genes are shown as gray dots. B CREB3L1 regulatory network. Blue dots represent skeletal system development-related target genes, other target genes are shown as gray dots. C RUNX2 and CREB3L1 regulon activity, embedded on SCENIC clustering layout and colored by RUNX2 or CREB3L1 regulon active levels. Gray dots represent completely inactive cells with zero active scores. Other dot colors represent active score in each cell from black (low) to blue (high). D Violin plots showing expression levels of RUNX2 (left) and CREB3L1 (right). *Adjusted p value < 0.05 (Kruskal–Wallis test). E Skeletal system development-related target gene expression in RUNX2 and CREB3L1 regulon among four cell clusters. Dot color indicates relative expression level and the dot size shows the proportion of cells expressing the indicated genes. F Violin plots showing average expression levels of all target genes in RUNX2 (left) and CREB3L1 (right) regulons. *Adjusted p value < 0.05 (Kruskal–Wallis test). G Continuum of dynamic RUNX2 regulon target gene expression in pseudotime of osteoblasts. Pixel color indicates expression level. Early-stage cells are listed at left. H Continuum of dynamic RUNX2 regulon’s target gene expression in pseudotime of osteoblasts. Pixel color indicates the expression levels. Early-stage cells are listed at left. I Motif swissregulon_hs_RUNX1..3.p2 of RUNX2. J Motif hocomoco_CR3L1_HUMAN.H11MO.0.D of CREB3L1. The relative sizes of the letters indicate the frequency of four bases in the sequences