Heterogeneity of Molecular Mechanisms Determining Blood Cancer Cell Lines Resistance to Natural Killer Cells in the Context of Tumor-Stromal Interactions: Insights from Studies of Pooled "DNA-Barcoded" Cell Line Panels

Immunotherapy has remarkably changed the treatment paradigm in hematologic malignancies and natural killer (NK) cell therapy represents an attractive option, as it has been feasible and safe in early clinical trials, without graft-versus-host effects. Nevertheless, the molecular markers determining cancer cell sensitivity or resistance to NK cells, especially in the context of tumor cell interaction with the bone marrow (BM) stromal microenvironment remain incompletely understood, but have major translational relevance since these tumor-stromal interactions have been known to attenuate the response of blood cancer cells to diverse classes of pharmacological agents.

To address these questions, we performed NK cell treatment of a series of cell lines from hematologic malignancies by applying a pooled "DNA-barcoded" format of these cell lines (PRISM system). We specifically quantified the dose-dependent responses to primary NK cells for 70 molecularly-annotated blood cancer cell lines, including myeloid and lymphoblastic leukemia, diffuse large B cell lymphoma and 15 multiple myeloma (MM) lines, either in presence or in absence of BM stromal cells (BMSCs) and interferon gamma (IFNg), followed by integrated computational analyses to identify candidate molecular markers correlating with tumor cell sensitivity or resistance to NK cells. NK cell cytotoxicity, quantified by the relative abundance of barcodes in treated cells compared to controls, was correlated with the transcriptional, mutational and other molecular features of each of the 70 cell lines from publicly available databases. Furthermore, data from MM cell lines were compared to our genome-wide loss of function (LOF) and gain of function (GOF) CRISPR screen data in the MM cell line MM.1S.

Two distinct clusters of cell lines, sensitive and resistant to NK cell treatment, were identified. Such clusters retained distinct pattern of in vitro resistance in the presence of stroma, while showing an overall markedly decreased NK cell responsiveness, which underscores the protective effect of stromal microenvironment in blood malignancies, regardless of the addition of IFNg. RNA-seq data showed no differences in dependencies between the two clusters and no distinct gene expression patterns at baseline that clearly allows to predict NK cell response, which underscores the heterogeneity of resistance patterns at single gene level, across different hematologic malignancies. However, when comparing baseline RNAseq data to data obtained from previous GOF and LOF CRISPR screens in MM.1S, surface antigens such as PVR, ULBP1, ULBP3 were more frequently downregulated, whereas MUC1 was upregulated in resistant cells clusters. An important observation is that gene lesions such as TP53, PTEN, MMSET, commonly associated with high-risk diseases, do not affect NK cell responses in the cell lines tested. Interestingly, a gene set enrichment analysis (GSEA) showed that the cluster of resistant cells displays upregulation of class I MHC complex, class II MHC complex binding, IL7 pathway and a downregulation of transmembrane receptor protein serine/threonine kinase signaling pathway. GSEA also showed that baseline state of IFN-JAK-STAT signaling correlates with BMSCs-induced NK cell resistance, a result further confirmed by addition of IFNgto tumor-NK cocultures in the absence of BMSCs. No significant differences in NK cell response were observed when comparing cell lines of different hematologic neoplasms, suggesting that candidate markers from these studies may be relevant across different hematologic malignancies.

In conclusion, this is the first study of this size correlating the molecular annotation of different concurrently-treated hematologic cell lines with their response to a NK-based treatment in the context of BMSC interaction. This study of a large panel of pooled "DNA-barcoded" cell lines provided complementary and orthogonal information to our LOF and GOF screens, expanding our potential to identify and validate molecular markers for individualized use of NK cell-based therapies in hematologic malignancies.