Transcriptomic profiling of normal plasma cells and MM models using bulk RNA-Seq. Bulk RNA-Seq analysis was performed using 50 000 sorted CD138⁺B220⁻CD45.2⁺ BM cells of WT (n = 3), V (n = 3), recipients of tVk12653 (n = 2), VQ-D1 (n = 2), VQ-D4 (n = 1), and VQ-D5 (n = 2). (A) Normalized expression levels of Myc in all samples. *P < .05; **P < .01; ***P < .001. (B) Scatterplots of average (avg) log₁₀ FPKM values of all mapped genes of V, tVk12653, or VQ samples vs WT samples. Red and blue dots represent differentially upregulated and downregulated genes, respectively, for each comparison (log₂ fold change [FC] >1 or <−1, adjusted P [P.adj], <.05); r² value calculated from Pearson correlation coefficient. (C) Volcano plots illustrating the differentially expressed genes in tVk12653 or VQ MM cells vs WT plasma cells. GSEA was performed between tVk12653 or VQ and WT plasma cells against Hallmark (D) or UAMS-70 (E) gene sets. (F) Volcano plots illustrating the differentially expressed genes in VQ vs tVk12653 MM cells. GSEA analysis was performed between VQ and tVk12653 MM cells against Hallmark (G) or C2 curated (H) gene sets. Adjusted P value and normalized enrichment score (NES) are shown on each GSEA plot.