Gene expression profiles of Ezh1-deficient HSPCs. (A) Venn diagram showing the overlap between genes upregulated (at least twofold) and downregulated (at least twofold) in Ezh1^−/− and Ezh2^Δ/Δ LSK cells compared with WT LSK cells. RNA-sequence analysis was performed using LSK cells from WT, Ezh1^−/−, and Ezh2^Δ/Δ mice at 1 month postinjection of tamoxifen. (B) GSEA plots demonstrating enrichment levels of indicated gene sets in Ezh1^−/− and Ezh2^Δ/Δ LSK cells compared with WT LSK cells. NES, NOM P value, and FDR are indicated. (C) Box-and-whisker plots showing the expression levels of H3K27me3-only genes and bivalent genes in Ezh1^−/− and Ezh2^Δ/Δ LSK cells relative to WT LSK cells. Boxes represent 25 to 75 percentile ranges. Vertical lines represent 10 to 90 percentile ranges. Horizontal bars represent median. ***P < .001 (Student t test.). (D) Graphic presentation of RNA-sequence data. RPKM of Ezh1 in Ezh1^−/− and Ezh2^Δ/Δ LSK cells were shown relative to that of WT LSK cells. (E) Expression of EZH1 and EZH2 in CD34⁺ MDS cells. Relative expression of EZH1 and EZH2 in CD34⁺ cells from healthy controls, MDS patients, and MDS patients with 7q deletion. The data are presented as scatter diagrams with median values (bars). ***P < .001 (Student t test). The data from MDS patients were retrieved from published database.³⁸