Fig. 4

Dynamic changes of CREM regulon activity. A CREM regulatory network. Red dots represent immunity-related target genes, blue dots represent cell proliferation/differentiation-related target genes, yellow dots represent both immunity and cell proliferation/differentiation-related target genes; other target genes are shown as gray dots. B CREM regulon activity, embedded on SCENIC clustering layout and colored by CREM 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). C Violin plots showing expression levels of CREM. *Adjusted p value < 0.05 (Kruskal–Wallis test). D Violin plots showing average expression levels of all target genes in CREM regulon. *Adjusted p value < 0.05 (Kruskal–Wallis test). E Immunity and cell proliferation/differentiation-related target gene expression in CREM regulon among four cell clusters. Dot color indicates the relative expression levels and the dot size shows the proportion of cells expressing the indicated genes. F Continuum of dynamic target gene expression in pseudotime of osteoblasts. Pixel color indicates the expression levels. Early-stage cells are list in on the left. G Motif taipale_cyt_meth_ATF4_GGATGACGTCATCC_eDBD of CREM. The relative sizes of the letters indicate the frequency of four bases in the sequences