A Landscape of Pharmacogenomic Interactions in Cancer.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

Effective drug combinations in breast, colon and pancreatic cancer cells.

Combinations of anti-cancer drugs can overcome resistance and provide new treatments1,2. The number of possible drug combinations vastly exceeds what could be tested clinically. Efforts to systematically identify active combinations and the tissues and molecular contexts in which they are most effective could accelerate the development of combination treatments. Here we evaluate the potency and efficacy of 2,025 clinically relevant two-drug combinations, generating a dataset encompassing 125 molecularly characterized breast, colorectal and pancreatic cancer cell lines. We show that synergy between drugs is rare and highly context-dependent, and that combinations of targeted agents are most likely to be synergistic. We incorporate multi-omic molecular features to identify combination biomarkers and specify synergistic drug combinations and their active contexts, including in basal-like breast cancer, and microsatellite-stable or KRAS-mutant colon cancer. Our results show that irinotecan and CHEK1 inhibition have synergistic effects in microsatellite-stable or KRAS-TP53 double-mutant colon cancer cells, leading to apoptosis and suppression of tumour xenograft growth. This study identifies clinically relevant effective drug combinations in distinct molecular subpopulations and is a resource to guide rational efforts to develop combinatorial drug treatments.

Predicting patient response with models trained on cell lines and patient-derived xenografts by nonlinear transfer learning.

Preclinical models have been the workhorse of cancer research, producing massive amounts of drug response data. Unfortunately, translating response biomarkers derived from these datasets to human tumors has proven to be particularly challenging. To address this challenge, we developed TRANSACT, a computational framework that builds a consensus space to capture biological processes common to preclinical models and human tumors and exploits this space to construct drug response predictors that robustly transfer from preclinical models to human tumors. TRANSACT performs favorably compared to four competing approaches, including two deep learning approaches, on a set of 23 drug prediction challenges on The Cancer Genome Atlas and 226 metastatic tumors from the Hartwig Medical Foundation. We demonstrate that response predictions deliver a robust performance for a number of therapies of high clinical importance: platinum-based chemotherapies, gemcitabine, and paclitaxel. In contrast to other approaches, we demonstrate the interpretability of the TRANSACT predictors by correctly identifying known biomarkers of targeted therapies, and we propose potential mechanisms that mediate the resistance to two chemotherapeutic agents.

Subunit composition of VRAC channels determines substrate specificity and cellular resistance to Pt-based anti-cancer drugs.

Although platinum-based drugs are widely used chemotherapeutics for cancer treatment, the determinants of tumor cell responsiveness remain poorly understood. We show that the loss of subunits LRRC8A and LRRC8D of the heteromeric LRRC8 volume-regulated anion channels (VRACs) increased resistance to clinically relevant cisplatin/carboplatin concentrations. Under isotonic conditions, about 50% of cisplatin uptake depended on LRRC8A and LRRC8D, but neither on LRRC8C nor on LRRC8E. Cell swelling strongly enhanced LRRC8-dependent cisplatin uptake, bolstering the notion that cisplatin enters cells through VRAC. LRRC8A disruption also suppressed drug-induced apoptosis independently from drug uptake, possibly by impairing VRAC-dependent apoptotic cell volume decrease. Hence, by mediating cisplatin uptake and facilitating apoptosis, VRAC plays a dual role in the cellular drug response. Incorporation of the LRRC8D subunit into VRAC substantially increased its permeability for cisplatin and the cellular osmolyte taurine, indicating that LRRC8 proteins form the channel pore. Our work suggests that LRRC8D-containing VRACs are crucial for cell volume regulation by an important organic osmolyte and may influence cisplatin/carboplatin responsiveness of tumors.

Pathogenic and targetable genetic alterations in 70 urachal adenocarcinomas.

Urachal cancer (UrC) is a rare but aggressive malignancy often diagnosed in advanced stages requiring systemic treatment. Although cytotoxic chemotherapy is of limited effectiveness, prospective clinical studies can hardly be conducted. Targeted therapeutic treatment approaches and potentially immunotherapy based on a biological rationale may provide an alternative strategy. We therefore subjected 70 urachal adenocarcinomas to targeted next-generation sequencing, conducted in situ and immunohistochemical analyses (including PD-L1 and DNA mismatch repair proteins [MMR]) and evaluated the microsatellite instability (MSI) status. The analytical findings were correlated with clinicopathological and outcome data and Kaplan-Meier and univariable/multivariable Cox regression analyses were performed. The patients had a mean age of 50 years, 66% were male and a 5-year overall survival (OS) of 58% and recurrence-free survival (RFS) of 45% was detected. Sequence variations were observed in TP53 (66%), KRAS (21%), BRAF (4%), PIK3CA (4%), FGFR1 (1%), MET (1%), NRAS (1%), and PDGFRA (1%). Gene amplifications were found in EGFR (5%), ERBB2 (2%), and MET (2%). We detected no evidence of MMR-deficiency (MMR-d)/MSI-high (MSI-h), whereas 10 of 63 cases (16%) expressed PD-L1. Therefore, anti-PD-1/PD-L1 immunotherapy approaches might be tested in UrC. Importantly, we found aberrations in intracellular signal transduction pathways (RAS/RAF/PI3K) in 31% of UrCs with potential implications for anti-EGFR therapy. Less frequent potentially actionable genetic alterations were additionally detected in ERBB2 (HER2), MET, FGFR1, and PDGFRA. The molecular profile strengthens the notion that UrC is a distinct entity on the genomic level with closer resemblance to colorectal than to bladder cancer.

TANDEM: a two-stage approach to maximize interpretability of drug response models based on multiple molecular data types.

MOTIVATION: Clinical response to anti-cancer drugs varies between patients. A large portion of this variation can be explained by differences in molecular features, such as mutation status, copy number alterations, methylation and gene expression profiles. We show that the classic approach for combining these molecular features (Elastic Net regression on all molecular features simultaneously) results in models that are almost exclusively based on gene expression. The gene expression features selected by the classic approach are difficult to interpret as they often represent poorly studied combinations of genes, activated by aberrations in upstream signaling pathways. RESULTS: To utilize all data types in a more balanced way, we developed TANDEM, a two-stage approach in which the first stage explains response using upstream features (mutations, copy number, methylation and cancer type) and the second stage explains the remainder using downstream features (gene expression). Applying TANDEM to 934 cell lines profiled across 265 drugs (GDSC1000), we show that the resulting models are more interpretable, while retaining the same predictive performance as the classic approach. Using the more balanced contributions per data type as determined with TANDEM, we find that response to MAPK pathway inhibitors is largely predicted by mutation data, while predicting response to DNA damaging agents requires gene expression data, in particular SLFN11 expression. AVAILABILITY AND IMPLEMENTATION: TANDEM is available as an R package on CRAN (for more information, see http://ccb.nki.nl/software/tandem). CONTACT: m.michaut@nki.nl or l.wessels@nki.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Emergence of De Novo Conditions Following Remission of Cushing Syndrome: A Case Report and Scoping Review.

OBJECTIVE: Onset and exacerbation of autoimmune, inflammatory or steroid-responsive conditions have been reported following the remission of Cushing syndrome, leading to challenges in distinguishing a new condition versus expected symptomatology following remission. We describe a case of a 42-year-old man presenting with new-onset sarcoidosis diagnosed 12 months following the surgical cure of Cushing syndrome and synthesise existing literature reporting on de novo conditions presenting after Cushing syndrome remission. METHODS: A scoping review was conducted in Medline, Epub, Ovid and PubMed. Case reports and case series detailing adult patients presenting with new-onset conditions following Cushing syndrome remission were included. RESULTS: In total, 1641 articles were screened, 138 full-text studies were assessed for eligibility, and 43 studies were included, of which 84 cases (including our case) were identified. Most patients were female (85.7%), and the median reported age was 39.5 years old (IQR = 13). Thyroid diseases were the most commonly reported conditions (48.8%), followed by sarcoidosis (15.5%). Psoriasis, lymphocytic hypophysitis, idiopathic intracranial hypertension, multiple sclerosis, rheumatoid arthritis, lupus and seronegative arthritis were reported in more than one case. The median duration between Cushing remission and de novo condition diagnosis was 4.1 months (IQR = 3.75). Of those patients, 59.5% were receiving corticosteroid therapy at the time of onset. CONCLUSION: Our scoping review identified several cases of de novo conditions emerging following the remission of Cushing syndrome. They occurred mostly in women and within the year following remission. Clinicians should remain aware that new symptoms, particularly in the first year following the treatment of Cushing syndrome, may be manifestations of a wide range of conditions aside from adrenal insufficiency or glucocorticoid withdrawal syndrome.

A pan-cancer screen identifies drug combination benefit in cancer cell lines at the individual and population level.

Combining drugs can enhance their clinical efficacy, but the number of possible combinations and inter-tumor heterogeneity make identifying effective combinations challenging, while existing approaches often overlook clinically relevant activity. We screen one of the largest cell line panels (N = 757) with 51 clinically relevant combinations and identify responses at the level of individual cell lines and tissue populations. We establish three response classes to model cellular effects beyond monotherapy: synergy, Bliss additivity, and independent drug action (IDA). Synergy is rare (11% of responses) and frequently efficacious (>50% viability reduction), whereas Bliss and IDA are more frequent but less frequently efficacious. We introduce "efficacious combination benefit" (ECB) to describe high-efficacy responses classified as either synergy, Bliss, or IDA. We identify ECB biomarkers in vitro and show that ECB predicts response in patient-derived xenografts better than synergy alone. Our work here provides a valuable resource and framework for preclinical evaluation and the development of combination treatments.

Spatial relationships in the urothelial and head and neck tumor microenvironment predict response to combination immune checkpoint inhibitors.

Immune checkpoint inhibitors (ICI) can achieve remarkable responses in urothelial cancer (UC), which may depend on tumor microenvironment (TME) characteristics. However, the relationship between the TME, usually characterized by immune cell density, and response to ICI is unclear. Here, we quantify the TME immune cell densities and spatial relationships (SRs) of 24 baseline UC samples, obtained before pre-operative combination ICI treatment, using multiplex immunofluorescence. We describe SRs by approximating the first nearest-neighbor distance distribution with a Weibull distribution and evaluate the association between TME metrics and ipilimumab+nivolumab response. Immune cell density does not discriminate between response groups. However, the Weibull SR metrics of CD8+ T cells or macrophages to their closest cancer cell positively associate with response. CD8+ T cells close to B cells are characteristic of non-response. We validate our SR response associations in a combination ICI cohort of head and neck tumors. Our data confirm that SRs, in contrast to density metrics, are strong biomarkers of response to pre-operative combination ICIs.

Platinum-Based Chemotherapy Induces Opposing Effects on Immunotherapy Response-Related Spatial and Stromal Biomarkers in the Bladder Cancer Microenvironment.

PURPOSE: Platinum-based chemotherapy and immune checkpoint inhibitors are key components of systemic treatment for muscle-invasive and advanced urothelial cancer. The ideal integration of these two treatment modalities remains unclear as clinical trials have led to inconsistent results. Modulation of the tumor-immune microenvironment by chemotherapy is poorly characterized. We aimed to investigate this modulation, focusing on potential clinical implications for immune checkpoint inhibitor response. EXPERIMENTAL DESIGN: We assessed immune cell densities, spatial relations, and tumor/stromal components from 116 patients with urothelial bladder cancer (paired data for 95 patients) before and after platinum-based chemotherapy. RESULTS: Several published biomarkers for immunotherapy response changed upon chemotherapy treatment. The intratumoral CD8+ T-cell percentage increased after treatment and was associated with increased TNFα-via-NF-κB signaling. The percentage of PDL1+ immune cells was higher after chemotherapy. An increase in chemo-induced changes that potentially inhibit an antitumor immune response was also observed, including increased fibroblast-based TGFß signaling and distances from immune cells to the nearest cancer cell. The latter two parameters correlated significantly in posttreatment samples, suggesting that TGFß signaling in fibroblasts may play a role in spatially separating immune cells from cancer cells. We examined specific chemotherapy regimens and found that treatment with methotrexate, vinblastine, doxorubicin, and cisplatin was associated with an increase in the macrophage cell percentage. Gemcitabine-containing chemotherapy was associated with upregulation of fibroblast TGFß signaling. CONCLUSIONS: The opposing effects of platinum-based chemotherapy on the immune cell composition and stromal context of the tumor-immune microenvironment may explain the inconsistent results of clinical trials investigating chemotherapy and immune checkpoint inhibitor combinations in bladder cancer.

Genomic and transcriptomic features of androgen receptor signaling inhibitor resistance in metastatic castration-resistant prostate cancer.

BACKGROUNDAndrogen receptor signaling inhibitors (ARSIs) have improved outcomes for patients with metastatic castration-resistant prostate cancer (mCRPC), but their clinical benefit is limited by treatment resistance.METHODSTo investigate the mechanisms of ARSI resistance, we analyzed the whole-genome (n = 45) and transcriptome (n = 31) sequencing data generated from paired metastatic biopsies obtained before initiation of first-line ARSI therapy for mCRPC and after radiographic disease progression. We investigated the effects of genetic and pharmacologic modulation of SSTR1 in 22Rv1 cells, a representative mCRPC cell line.RESULTSWe confirmed the predominant role of tumor genetic alterations converging on augmenting androgen receptor (AR) signaling and the increased transcriptional heterogeneity and lineage plasticity during the emergence of ARSI resistance. We further identified amplifications involving a putative enhancer downstream of the AR and transcriptional downregulation of SSTR1, encoding somatostatin receptor 1, in ARSI-resistant tumors. We found that patients with SSTR1-low mCRPC tumors derived less benefit from subsequent ARSI therapy in a retrospective cohort. We showed that SSTR1 was antiproliferative in 22Rv1 cells and that the FDA-approved drug pasireotide suppressed 22Rv1 cell proliferation.CONCLUSIONOur findings expand the knowledge of ARSI resistance and point out actionable next steps, exemplified by potentially targeting SSTR1, to improve patient outcomes.FUNDINGNational Cancer Institute (NCI), NIH; Prostate Cancer Foundation; Conquer Cancer, American Society of Clinical Oncology Foundation; UCSF Benioff Initiative for Prostate Cancer Research; Netherlands Cancer Institute.

Assessment of Predictive Genomic Biomarkers for Response to Cisplatin-based Neoadjuvant Chemotherapy in Bladder Cancer.

Cisplatin-based neoadjuvant chemotherapy (NAC) followed by radical cystectomy is recommended for patients with muscle-invasive bladder cancer (MIBC). It has been shown that somatic deleterious mutations in ERCC2, gain-of-function mutations in ERBB2, and alterations in ATM, RB1, and FANCC are correlated with pathological response to NAC in MIBC. The objective of this study was to validate these genomic biomarkers in pretreatment transurethral resection material from an independent retrospective cohort of 165 patients with MIBC who subsequently underwent NAC and radical surgery. Patients with ypT0/Tis/Ta/T1N0 disease after surgery were defined as responders. Somatic deleterious mutations in ERCC2 were found in nine of 68 (13%) evaluable responders and two of 95 (2%) evaluable nonresponders (p = 0.009; FDR = 0.03). No correlation was observed between response and alterations in ERBB2 or in ATM, RB1, or FANCC alone or in combination. In an exploratory analysis, no additional genomic alterations discriminated between responders and nonresponders to NAC. No further associations were identified between the aforementioned biomarkers and pathological complete response (ypT0N0) after surgery. In conclusion, we observed a positive association between deleterious mutations in ERCC2 and pathological response to NAC, but not overall survival or recurrence-free survival. Other previously reported genomic biomarkers were not validated. PATIENT SUMMARY: It is currently unknown which patients will respond to chemotherapy before definitive surgery for bladder cancer. Previous studies described several gene mutations in bladder cancer that correlated with chemotherapy response. This study confirmed that patients with bladder cancer with a mutation in the ERCC2 gene often respond to chemotherapy.

A kinome-centered CRISPR-Cas9 screen identifies activated BRAF to modulate enzalutamide resistance with potential therapeutic implications in BRAF-mutated prostate cancer.

Resistance to drugs targeting the androgen receptor (AR) signaling axis remains an important challenge in the treatment of prostate cancer patients. Activation of alternative growth pathways is one mechanism used by cancer cells to proliferate despite treatment, conferring drug resistance. Through a kinome-centered CRISPR-Cas9 screen in CWR-R1 prostate cancer cells, we identified activated BRAF signaling as a determinant for enzalutamide resistance. Combined pharmaceutical targeting of AR and MAPK signaling resulted in strong synergistic inhibition of cell proliferation. The association between BRAF activation and enzalutamide resistance was confirmed in two metastatic prostate cancer patients harboring activating mutations in the BRAF gene, as both patients were unresponsive to enzalutamide. Our findings suggest that co-targeting of the MAPK and AR pathways may be effective in patients with an activated MAPK pathway, particularly in patients harboring oncogenic BRAF mutations. These results warrant further investigation of the response to AR inhibitors in BRAF-mutated prostate tumors in clinical settings.

The Tumor Immune Landscape and Architecture of Tertiary Lymphoid Structures in Urothelial Cancer.

Candidate immune biomarkers have been proposed for predicting response to immunotherapy in urothelial cancer (UC). Yet, these biomarkers are imperfect and lack predictive power. A comprehensive overview of the tumor immune contexture, including Tertiary Lymphoid structures (TLS), is needed to better understand the immunotherapy response in UC. We analyzed tumor sections by quantitative multiplex immunofluorescence to characterize immune cell subsets in various tumor compartments in tumors without pretreatment and tumors exposed to preoperative anti-PD1/CTLA-4 checkpoint inhibitors (NABUCCO trial). Pronounced immune cell presence was found in UC invasive margins compared to tumor and stroma regions. CD8+PD1+ T-cells were present in UC, particularly following immunotherapy. The cellular composition of TLS was assessed by multiplex immunofluorescence (CD3, CD8, FoxP3, CD68, CD20, PanCK, DAPI) to explore specific TLS clusters based on varying immune subset densities. Using a k-means clustering algorithm, we found five distinct cellular composition clusters. Tumors unresponsive to anti-PD-1/CTLA-4 immunotherapy showed enrichment of a FoxP3+ T-cell-low TLS cluster after treatment. Additionally, cluster 5 (macrophage low) TLS were significantly higher after pre-operative immunotherapy, compared to untreated tumors. We also compared the immune cell composition and maturation stages between superficial (submucosal) and deeper TLS, revealing that superficial TLS had more pronounced T-helper cells and enrichment of early TLS than TLS located in deeper tissue. Furthermore, superficial TLS displayed a lower fraction of secondary follicle like TLS than deeper TLS. Taken together, our results provide a detailed quantitative overview of the tumor immune landscape in UC, which can provide a basis for further studies.

Ovarian Cancer-Specific BRCA-like Copy-Number Aberration Classifiers Detect Mutations Associated with Homologous Recombination Deficiency in the AGO-TR1 Trial.

PURPOSE: Previously, we developed breast cancer BRCA1-like and BRCA2-like copy-number profile shrunken centroid classifiers predictive for mutation status and response to therapy, targeting homologous recombination deficiency (HRD). Therefore, we investigated BRCA1- and BRCA2-like classification in ovarian cancer, aiming to acquire classifiers with similar properties as those in breast cancer.Experimental Design: We analyzed DNA copy-number profiles of germline BRCA1- and BRCA2-mutant ovarian cancers and control tumors and observed that existing breast cancer classifiers did not sufficiently predict mutation status. Hence, we trained new shrunken centroid classifiers on this set and validated them in the independent The Cancer Genome Atlas dataset. Subsequently, we assessed BRCA1/2-like classification and obtained germline and tumor mutation and methylation status of cancer predisposition genes, among them several involved in HR repair, of 300 ovarian cancer samples derived from the consecutive cohort trial AGO-TR1 (NCT02222883). RESULTS: The detection rate of the BRCA1-like classifier for BRCA1 mutations and promoter hypermethylation was 95.6%. The BRCA2-like classifier performed less accurately, likely due to a smaller training set. Furthermore, three quarters of the BRCA1/2-like tumors could be explained by (epi)genetic alterations in BRCA1/2, germline RAD51C mutations and alterations in other genes involved in HR. Around half of the non-BRCA-mutated ovarian cancer cases displayed a BRCA-like phenotype. CONCLUSIONS: The newly trained classifiers detected most BRCA-mutated and methylated cancers and all tumors harboring a RAD51C germline mutations. Beyond that, we found an additional substantial proportion of ovarian cancers to be BRCA-like.

Endocrine pulse identification using penalized methods and a minimum set of assumptions.

The detection of hormone secretion episodes is important for understanding normal and abnormal endocrine functioning, but pulse identification from hormones measured with short interval sampling is challenging. Furthermore, to obtain useable results, the model underlying hormone secretion and clearance must be augmented with restrictions based on biologically acceptable assumptions. Here, using the assumption that there are only a few time points at which a hormone is secreted, we used a modern penalized nonlinear least-squares setup to select the number of secretion events. We did not assume a particular shape or frequency distribution for the secretion pulses. Our pulse identfication method, VisPulse, worked well with luteinizing hormone (LH), cortisol, growth hormone, or testosterone. In particular, applying our modeling strategy to previous LH data revealed a good correlation between the modeled and measured LH hormone concentrations, the estimated secretion pattern was sparse, and the small and structureless residuals indicated a proper model with a good fit. We benchmarked our method to AutoDecon, a commonly used hormone secretion model, and performed releasing hormone infusion experiments. The results of these experiments confirmed that our method is accurate and outperforms AutoDecon, especially for detecting silent periods and small secretion events, suggesting a high-secretion event resolution. Method validation using (releasing hormone) infusion data revealed sensitivities and selectivities of 0.88 and 0.95 and of 0.69 and 0.91 for VisPulse and AutoDecon, respectively.

Preoperative ipilimumab plus nivolumab in locoregionally advanced urothelial cancer: the NABUCCO trial.

Preoperative immunotherapy with anti-PD1 plus anti-CTLA4 antibodies has shown remarkable pathological responses in melanoma1 and colorectal cancer2. In NABUCCO (ClinicalTrials.gov: NCT03387761 ), a single-arm feasibility trial, 24 patients with stage III urothelial cancer (UC) received two doses of ipilimumab and two doses of nivolumab, followed by resection. The primary endpoint was feasibility to resect within 12 weeks from treatment start. All patients were evaluable for the study endpoints and underwent resection, 23 (96%) within 12 weeks. Grade 3-4 immune-related adverse events occurred in 55% of patients and in 41% of patients when excluding clinically insignificant laboratory abnormalities. Eleven patients (46%) had a pathological complete response (pCR), meeting the secondary efficacy endpoint. Fourteen patients (58%) had no remaining invasive disease (pCR or pTisN0/pTaN0). In contrast to studies with anti-PD1/PD-L1 monotherapy, complete response to ipilimumab plus nivolumab was independent of baseline CD8+ presence or T-effector signatures. Induction of tertiary lymphoid structures upon treatment was observed in responding patients. Our data indicate that combined CTLA-4 plus PD-1 blockade might provide an effective preoperative treatment strategy in locoregionally advanced UC, irrespective of pre-existing CD8+ T cell activity.

Genomic data integration by WON-PARAFAC identifies interpretable factors for predicting drug-sensitivity in vivo.

Integrative analyses that summarize and link molecular data to treatment sensitivity are crucial to capture the biological complexity which is essential to further precision medicine. We introduce Weighted Orthogonal Nonnegative parallel factor analysis (WON-PARAFAC), a data integration method that identifies sparse and interpretable factors. WON-PARAFAC summarizes the GDSC1000 cell line compendium in 130 factors. We interpret the factors based on their association with recurrent molecular alterations, pathway enrichment, cancer type, and drug-response. Crucially, the cell line derived factors capture the majority of the relevant biological variation in Patient-Derived Xenograft (PDX) models, strongly suggesting our factors capture invariant and generalizable aspects of cancer biology. Furthermore, drug response in cell lines is better and more consistently translated to PDXs using factor-based predictors as compared to raw feature-based predictors. WON-PARAFAC efficiently summarizes and integrates multiway high-dimensional genomic data and enhances translatability of drug response prediction from cell lines to patient-derived xenografts.