Transcriptional changes following Kdm1a depletion in Eμ-TCL1A mice. (A) PCA of RNA-seq data reveals distinct groups by genotype. Each data point represents 1 mouse. (B) Left: volcano plot showing differential gene expression between splenic MCs from Eμ-TCL1A vs from iKdm1a^KD;Eμ-TCL1A mice (light blue: upregulated in iKdm1a^KD;Eμ-TCL1A, dark blue: downregulated in iKdm1a^KD;Eμ-TCL1A). Right: pie chart indicating the percentage of up- and downregulated genes. One thousand thirteen genes were differentially expressed between splenic MCs from Eμ-TCL1A vs from iKdm1a^KD;Eμ-TCL1A mice, FCh > 1.5, q ≤ 0.01, with most of them (81%) upregulated in the iKdm1a^KD;Eμ-TCL1A model. (C) GSEA of DEGs between the Eμ-TCL1A and iKdm1a^KD;Eμ-TCL1A systems. The bars represent the percentage of genes from the data set that mapped to the corresponding KEGG pathways. Light blue: upregulated genes and dark blue: downregulated genes in iKdm1a^KD;Eμ-TCL1A splenocytes, FCh > 1.5, P ≤ .01. (D) Enrichment plots of GSEA of apoptosis (left) and p53 (right) pathways. NES and FDR q values are indicated. (E) Heatmap of DEGs associated with p53- and proapoptotic pathways. (F) Heatmap of DEGs associated with the pathway clusters of cell adhesion and migration. (G) Quantitative reverse transcriptase polymerase chain reactionanalysis demonstrating reduced expression of Itgb2l, Rapgef3, and Cxcr4 in splenic MCs and in PB leukocytes of iKdm1a^KD;Eμ-TCL1A vs Eμ-TCL1A mice (all P < .05, Student t test, mean ± SEM).