Alteration of NK cell transcriptomic landscape in patients with MM. A scRNA-seq analysis of fresh blood and BM CD3⁻CD56⁺ NK cells from 10 patients with NDMM and 10 age/sex–matched HDs was performed. (A) Uniform manifold approximation and projection (UMAP) embeddings of NK cells colored by clusters. (B) Variations of NK cluster frequency across disease and tissue status (box-and-whisker plots, in which whiskers represent 1.5∗ interquartile range and a box represents the 25th, 50th, and 75th percentiles; false discovery rate [FDR] derived from generalized linear model; see supplemental Methods; FDR < 0.05∗, FDR < 0.01∗∗, and FDR < 0.001∗∗∗). (C) MELD (manifold enhancement of latent dimensions) likelihood (see supplemental Methods) highlighting MM vs HD enrichment computed separately on BM or PB. Likelihood values are overlaid on UMAP embeddings and per-cluster likelihood distributions displayed as violin plots. (D) Heatmap of marker genes standardized expression across clusters. (E) Heat map of standardized gene module scores from published data sets (Crinier et al²⁵ and Yang et al²⁶).