Fig. 4

Identification of hub genes and functional annotation. A Differential analysis among three H4M modification patterns. B Heat map of the correlations between the module eigengenes and H4M modification patterns. The Pearson correlation coefficient and corresponding p value were displayed in each cell. C Scatterplots of correlation between MEturquoise membership and gene significance for H4Mcluster-B and H4Mcluster-A. D KEGG enrichment analysis for 133 hub genes identified by WGCNA analysis. E Enriched pathway regulatory network for 20 pathways. F KEGG enrichment analysis for 39 feature genes identified by LASSO regression analysis. G Enriched pathway regulatory network for 13 pathways. H The gene regulatory network for the cell cycle signaling pathway