Fig. 6

Analysis of the major ZGA genes affected by HDACi at the 2-cell stage. (a) The Sankey diagrams display the proportions of different expression patterns of up- and downregulated DEGs in the three groups. The upper part contains the up- and downregulated DEGs set from the three groups affected by HDACi in this experiment, whereas the lower part represents the gene sets classified according to programmed waves during maternal-to-zygotic gene activation. (b) The combined plot represents the presence of the top 20 major ZGA genes with the highest effect under HDACi treatment in the enriched GO terms. The main heatmap shows the presence of these 20 major ZGA genes in the GO terms that were calculated from all HDACi-affected major ZGA genes. If a gene is present in a GO term, the corresponding rectangle in the heatmap is colored based on the evidence code, which indicates how the gene is annotated in that GO term. The lollipop plot depicted on the left side of the heatmap shows the fold change in expression of these 20 genes, where|LFC| denotes the absolute value of the log2 fold change and the upregulated genes are represented in red, whereas the downregulated genes are indicated in blue. The bar chart at the bottom of the heatmap shows enriched GO terms and was divided into four clusters. The height of each bar represents the proportion of an HDACi-affected major ZGA gene in the corresponding GO term, with the color representing the negative logarithm of the P-value obtained from the enrichment analysis