Patterns of Oncogene Coexpression at Single-Cell Resolution Influence Survival in Lymphoma.

Cancers often overexpress multiple clinically relevant oncogenes, but it is not known if combinations of oncogenes in cellular subpopulations within a cancer influence clinical outcomes. Using quantitative multispectral imaging of the prognostically relevant oncogenes MYC, BCL2, and BCL6 in diffuse large B-cell lymphoma (DLBCL), we show that the percentage of cells with a unique combination MYC+BCL2+BCL6- (M+2+6-) consistently predicts survival across four independent cohorts (n = 449), an effect not observed with other combinations including M+2+6+. We show that the M+2+6- percentage can be mathematically derived from quantitative measurements of the individual oncogenes and correlates with survival in IHC (n = 316) and gene expression (n = 2,521) datasets. Comparative bulk/single-cell transcriptomic analyses of DLBCL samples and MYC/BCL2/BCL6-transformed primary B cells identify molecular features, including cyclin D2 and PI3K/AKT as candidate regulators of M+2+6- unfavorable biology. Similar analyses evaluating oncogenic combinations at single-cell resolution in other cancers may facilitate an understanding of cancer evolution and therapy resistance. SIGNIFICANCE: Using single-cell-resolved multiplexed imaging, we show that selected subpopulations of cells expressing specific combinations of oncogenes influence clinical outcomes in lymphoma. We describe a probabilistic metric for the estimation of cellular oncogenic coexpression from IHC or bulk transcriptomes, with possible implications for prognostication and therapeutic target discovery in cancer. This article is highlighted in the In This Issue feature, p. 1027.

An ex vivo platform to guide drug combination treatment in relapsed/refractory lymphoma.

Although combination therapy is the standard of care for relapsed/refractory non-Hodgkin's lymphoma (RR-NHL), combination treatment chosen for an individual patient is empirical, and response rates remain poor in individuals with chemotherapy-resistant disease. Here, we evaluate an experimental-analytic method, quadratic phenotypic optimization platform (QPOP), for prediction of patient-specific drug combination efficacy from a limited quantity of biopsied tumor samples. In this prospective study, we enrolled 71 patients with RR-NHL (39 B cell NHL and 32 NK/T cell NHL) with a median of two prior lines of treatment, at two academic hospitals in Singapore from November 2017 to August 2021. Fresh biopsies underwent ex vivo testing using a panel of 12 drugs with known efficacy against NHL to identify effective single and combination treatments. Individualized QPOP reports were generated for 67 of 75 patient samples, with a median turnaround time of 6 days from sample collection to report generation. Doublet drug combinations containing copanlisib or romidepsin were most effective against B cell NHL and NK/T cell NHL samples, respectively. Off-label QPOP-guided therapy offered at physician discretion in the absence of standard options (n = 17) resulted in five complete responses. Among patients with more than two prior lines of therapy, the rates of progressive disease were lower with QPOP-guided treatments than with conventional chemotherapy. Overall, this study shows that the identification of patient-specific drug combinations through ex vivo analysis was achievable for RR-NHL in a clinically applicable time frame. These data provide the basis for a prospective clinical trial evaluating ex vivo-guided combination therapy in RR-NHL.