Efficacy of pembrolizumab and biomarker analysis in patients with WGS-based intermediate to high tumor mutational load: results from the Drug Rediscovery Protocol

PURPOSE

To evaluate efficacy of pembrolizumab across multiple cancer types harboring different levels of Whole-Genome Sequencing (WGS)-based tumor mutational load (TML; total of non-synonymous mutations across the genome) in patients included in the Drug Rediscovery Protocol (NCT02925234).

PATIENTS AND METHODS

Patients with solid, treatment-refractory, microsatellite-stable tumors were enrolled in cohort A: breast cancer TML 140-290, cohort B: tumor-agnostic cohort TML 140-290, and cohort C: tumor-agnostic cohort TML >290. Patients received pembrolizumab 200 mg every three weeks. Primary endpoint was clinical benefit (CB: objective response or stable disease (SD) ≥16 weeks). Pre-treatment tumor biopsies were obtained for WGS and RNA-sequencing.

RESULTS

Seventy-two evaluable patients with 26 different histotypes were enrolled. CB rate was 13% in cohort A (3/24 with partial response (PR)), 21% in cohort B (3/24 with SD, 2/24 with PR), and 42% in cohort C (4/24 with SD, 6/24 with PR). In cohort C, neoantigen burden estimates and expression of inflammation and innate immune biomarkers were significantly associated with CB. Similar associations were not identified in cohort A and B. In cohort A, CB was significantly associated with mutations in the chromatin remodeling gene PBRM1, while in cohort B, CB was significantly associated with expression of MICA/MICB and butyrophilins. CB and clonal TML were not significantly associated.

CONCLUSION

While in cohort A pembrolizumab lacked activity, cohort B and cohort C met the study's primary endpoint. Further research is warranted to refine selection of patients with tumors harboring lower TMLs and may benefit from a focus on innate immunity.

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.

Evaluating the effectiveness of pre-operative diagnosis of ovarian cancer using minimally invasive liquid biopsies by combining serum human epididymis protein 4 and cell-free DNA in patients with an ovarian mass

OBJECTIVE

To assess the feasibility of scalable, objective, and minimally invasive liquid biopsy-derived biomarkers such as cell-free DNA copy number profiles, human epididymis protein 4 (HE4), and cancer antigen 125 (CA125) for pre-operative risk assessment of early-stage ovarian cancer in a clinically representative and diagnostically challenging population and to compare the performance of these biomarkers with the Risk of Malignancy Index (RMI).

METHODS

In this case-control study, we included 100 patients with an ovarian mass clinically suspected to be early-stage ovarian cancer. Of these 100 patients, 50 were confirmed to have a malignant mass (cases) and 50 had a benign mass (controls). Using WisecondorX, an algorithm used extensively in non-invasive prenatal testing, we calculated the benign-calibrated copy number profile abnormality score. This score represents how different a sample is from benign controls based on copy number profiles. We combined this score with HE4 serum concentration to separate cases and controls.

RESULTS

Combining the benign-calibrated copy number profile abnormality score with HE4, we obtained a model with a significantly higher sensitivity (42% vs 0%; p<0.002) at 99% specificity as compared with the RMI that is currently employed in clinical practice. Investigating performance in subgroups, we observed especially large differences in the advanced stage and non-high-grade serous ovarian cancer groups.

CONCLUSION

This study demonstrates that cell-free DNA can be successfully employed to perform pre-operative risk of malignancy assessment for ovarian masses; however, results warrant validation in a more extensive clinical study.

Distinct spatiotemporal dynamics of CD8+ T cell-derived cytokines in the tumor microenvironment

Cells in the tumor microenvironment (TME) influence each other through secretion and sensing of soluble mediators, such as cytokines and chemokines. While signaling of interferon γ (IFNγ) and tumor necrosis factor α (TNFα) is integral to anti-tumor immune responses, our understanding of the spatiotemporal behavior of these cytokines is limited. Here, we describe a single cell transcriptome-based approach to infer which signal(s) an individual cell has received. We demonstrate that, contrary to expectations, CD8+ T cell-derived IFNγ is the dominant modifier of the TME relative to TNFα. Furthermore, we demonstrate that cell pools that show abundant IFNγ sensing are characterized by decreased expression of transforming growth factor β (TGFβ)-induced genes, consistent with IFNγ-mediated TME remodeling. Collectively, these data provide evidence that CD8+ T cell-secreted cytokines should be categorized into local and global tissue modifiers, and describe a broadly applicable approach to dissect cytokine and chemokine modulation of the TME.

Multi-omics analysis reveals distinct non-reversion mechanisms of PARPi resistance in BRCA1- versus BRCA2-deficient mammary tumors

BRCA1 and BRCA2 both function in DNA double-strand break repair by homologous recombination (HR). Due to their HR defect, BRCA1/2-deficient cancers are sensitive to poly(ADP-ribose) polymerase inhibitors (PARPis), but they eventually acquire resistance. Preclinical studies yielded several PARPi resistance mechanisms that do not involve BRCA1/2 reactivation, but their relevance in the clinic remains elusive. To investigate which BRCA1/2-independent mechanisms drive spontaneous resistance in vivo, we combine molecular profiling with functional analysis of HR of matched PARPi-naive and PARPi-resistant mouse mammary tumors harboring large intragenic deletions that prevent reactivation of BRCA1/2. We observe restoration of HR in 62% of PARPi-resistant BRCA1-deficient tumors but none in the PARPi-resistant BRCA2-deficient tumors. Moreover, we find that 53BP1 loss is the prevalent resistance mechanism in HR-proficient BRCA1-deficient tumors, whereas resistance in BRCA2-deficient tumors is mainly induced by PARG loss. Furthermore, combined multi-omics analysis identifies additional genes and pathways potentially involved in modulating PARPi response.

A living biobank of patient-derived ductal carcinoma in situ mouse-intraductal xenografts identifies risk factors for invasive progression

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer (IBC). Due to a lack of biomarkers able to distinguish high- from low-risk cases, DCIS is treated similar to early IBC even though the minority of untreated cases eventually become invasive. Here, we characterized 115 patient-derived mouse-intraductal (MIND) DCIS models reflecting the full spectrum of DCIS observed in patients. Utilizing the possibility to follow the natural progression of DCIS combined with omics and imaging data, we reveal multiple prognostic factors for high-risk DCIS including high grade, HER2 amplification, expansive 3D growth, and high burden of copy number aberrations. In addition, sequential transplantation of xenografts showed minimal phenotypic and genotypic changes over time, indicating that invasive behavior is an intrinsic phenotype of DCIS and supporting a multiclonal evolution model. Moreover, this study provides a collection of 19 distributable DCIS-MIND models spanning all molecular subtypes.

ecpc: an R-package for generic co-data models for high-dimensional prediction

BACKGROUND

High-dimensional prediction considers data with more variables than samples. Generic research goals are to find the best predictor or to select variables. Results may be improved by exploiting prior information in the form of co-data, providing complementary data not on the samples, but on the variables. We consider adaptive ridge penalised generalised linear and Cox models, in which the variable-specific ridge penalties are adapted to the co-data to give a priori more weight to more important variables. The R-package ecpc originally accommodated various and possibly multiple co-data sources, including categorical co-data, i.e. groups of variables, and continuous co-data. Continuous co-data, however, were handled by adaptive discretisation, potentially inefficiently modelling and losing information. As continuous co-data such as external p values or correlations often arise in practice, more generic co-data models are needed.

RESULTS

Here, we present an extension to the method and software for generic co-data models, particularly for continuous co-data. At the basis lies a classical linear regression model, regressing prior variance weights on the co-data. Co-data variables are then estimated with empirical Bayes moment estimation. After placing the estimation procedure in the classical regression framework, extension to generalised additive and shape constrained co-data models is straightforward. Besides, we show how ridge penalties may be transformed to elastic net penalties. In simulation studies we first compare various co-data models for continuous co-data from the extension to the original method. Secondly, we compare variable selection performance to other variable selection methods. The extension is faster than the original method and shows improved prediction and variable selection performance for non-linear co-data relations. Moreover, we demonstrate use of the package in several genomics examples throughout the paper.

CONCLUSIONS

The R-package ecpc accommodates linear, generalised additive and shape constrained additive co-data models for the purpose of improved high-dimensional prediction and variable selection. The extended version of the package as presented here (version number 3.1.1 and higher) is available on ( https://cran.r-project.org/web/packages/ecpc/ ).

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.

PD-L1 blockade in combination with carboplatin as immune induction in metastatic lobular breast cancer: the GELATO trial

Invasive lobular breast cancer (ILC) is the second most common histological breast cancer subtype, but ILC-specific trials are lacking. Translational research revealed an immune-related ILC subset, and in mouse ILC models, synergy between immune checkpoint blockade and platinum was observed. In the phase II GELATO trial ( NCT03147040 ), patients with metastatic ILC were treated with weekly carboplatin (area under the curve 1.5 mg ml-1 min-1) as immune induction for 12 weeks and atezolizumab (PD-L1 blockade; triweekly) from the third week until progression. Four of 23 evaluable patients had a partial response (17%), and 2 had stable disease, resulting in a clinical benefit rate of 26%. From these six patients, four had triple-negative ILC (TN-ILC). We observed higher CD8+ T cell infiltration, immune checkpoint expression and exhausted T cells after treatment. With this GELATO trial, we show that ILC-specific clinical trials are feasible and demonstrate promising antitumor activity of atezolizumab with carboplatin, particularly for TN-ILC, and provide insights for the design of highly needed ILC-specific trials.

Codon-specific KRAS mutations predict survival benefit of trifluridine/tipiracil in metastatic colorectal cancer

Genomics has greatly improved how patients with cancer are being treated; however, clinical-grade genomic biomarkers for chemotherapies are currently lacking. Using whole-genome analysis of 37 patients with metastatic colorectal cancer (mCRC) treated with the chemotherapy trifluridine/tipiracil (FTD/TPI), we identified KRAS codon G12 (KRAS^G12) mutations as a potential biomarker of resistance. Next, we collected real-world data of 960 patients with mCRC receiving FTD/TPI and validated that KRAS^G12 mutations were significantly associated with poor survival, also in analyses restricted to the RAS/RAF mutant subgroup. We next analyzed the data of the global, double-blind, placebo-controlled, phase 3 RECOURSE trial (n = 800 patients) and found that KRAS^G12 mutations (n = 279) were predictive biomarkers for reduced overall survival (OS) benefit of FTD/TPI versus placebo (unadjusted interaction P = 0.0031, adjusted interaction P = 0.015). For patients with KRAS^G12 mutations in the RECOURSE trial, OS was not prolonged with FTD/TPI versus placebo (n = 279; hazard ratio (HR) = 0.97; 95% confidence interval (CI) = 0.73-1.20; P = 0.85). In contrast, patients with KRAS^G13 mutant tumors showed significantly improved OS with FTD/TPI versus placebo (n = 60; HR = 0.29; 95% CI = 0.15-0.55; P < 0.001). In isogenic cell lines and patient-derived organoids, KRAS^G12 mutations were associated with increased resistance to FTD-based genotoxicity. In conclusion, these data show that KRAS^G12 mutations are biomarkers for reduced OS benefit of FTD/TPI treatment, with potential implications for approximately 28% of patients with mCRC under consideration for treatment with FTD/TPI. Furthermore, our data suggest that genomics-based precision medicine may be possible for a subset of chemotherapies.

Tumour-educated platelets for breast cancer detection: biological and technical insights

BACKGROUND

Studies have shown that blood platelets contain tumour-specific mRNA profiles tumour-educated platelets (TEPs). Here, we aim to train a TEP-based breast cancer detection classifier.

METHODS

Platelet mRNA was sequenced from 266 women with stage I-IV breast cancer and 212 female controls from 6 hospitals. A particle swarm optimised support vector machine (PSO-SVM) and an elastic net-based classifier (EN) were trained on 71% of the study population. Classifier performance was evaluated in the remainder (29%) of the population, followed by validation in an independent set (37 cases and 36 controls). Potential confounding was assessed in post hoc analyses.

RESULTS

Both classifiers reached an area under the curve (AUC) of 0.85 upon internal validation. Reproducibility in the independent validation set was poor with an AUC of 0.55 and 0.54 for the PSO-SVM and EN classifier, respectively. Post hoc analyses indicated that 19% of the variance in gene expression was associated with hospital. Genes related to platelet activity were differentially expressed between hospitals.

CONCLUSIONS

We could not validate two TEP-based breast cancer classifiers in an independent validation cohort. The TEP protocol is sensitive to within-protocol variation and revision might be necessary before TEPs can be reconsidered for breast cancer detection.

IL-5-producing CD4+ T cells and eosinophils cooperate to enhance response to immune checkpoint blockade in breast cancer

Immune checkpoint blockade (ICB) has heralded a new era in cancer therapy. Research into the mechanisms underlying response to ICB has predominantly focused on T cells; however, effective immune responses require tightly regulated crosstalk between innate and adaptive immune cells. Here, we combine unbiased analysis of blood and tumors from metastatic breast cancer patients treated with ICB with mechanistic studies in mouse models of breast cancer. We observe an increase in systemic and intratumoral eosinophils in patients and mice responding to ICB treatment. Mechanistically, ICB increased IL-5 production by CD4⁺ T cells, stimulating elevated eosinophil production from the bone marrow, leading to systemic eosinophil expansion. Additional induction of IL-33 by ICB-cisplatin combination or recombinant IL-33 promotes intratumoral eosinophil infiltration and eosinophil-dependent CD8⁺ T cell activation to enhance ICB response. This work demonstrates the critical role of eosinophils in ICB response and provides proof-of-principle for eosinophil engagement to enhance ICB efficacy.

Homologous Recombination Deficiency and Cyclin E1 Amplification Are Correlated with Immune Cell Infiltration and Survival in High-Grade Serous Ovarian Cancer

BACKGROUND

How molecular profiles are associated with tumor microenvironment (TME) in high-grade serous ovarian cancer (HGSOC) is incompletely understood. Therefore, we analyzed the TME and molecular profiles of HGSOC and assessed their associations with overall survival (OS).

METHODS

Patients with advanced-stage HGSOC treated in three Dutch hospitals between 2008-2015 were included. Patient data were collected from medical records. BRCA1/2 mutation, BRCA1 promotor methylation analyses, and copy number variations were used to define molecular profiles. Immune cells were assessed with immunohistochemical staining.

RESULTS

348 patients were categorized as BRCA mutation (BRCAm) (BRCAm or promotor methylation) (30%), non-BRCA mutated HRD (19%), Cyclin E1 (CCNE1)-amplification (13%), non-BRCAmut HRD and CCNE1-amplification (double classifier) (20%), and no specific molecular profile (NSMP) (18%). BRCAm showed highest immune cell densities and CCNE1-amplification lowest. BRCAm showed the most favorable OS (52.5 months), compared to non-BRCAmut HRD (41.0 months), CCNE1-amplification (28.0 months), double classifier (27.8 months), and NSMP (35.4 months). Higher immune cell densities showed a favorable OS compared to lower, also within the profiles. CD8+, CD20+, and CD103+ cells remained associated with OS in multivariable analysis.

CONCLUSIONS

Molecular profiles and TME are associated with OS. TME differs per profile, with higher immune cell densities showing a favorable OS, even within the profiles. HGSOC does not reflect one entity but comprises different entities based on molecular profiles and TME.

Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients' outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.

Truncated FGFR2 is a clinically actionable oncogene in multiple cancers

Somatic hotspot mutations and structural amplifications and fusions that affect fibroblast growth factor receptor 2 (encoded by FGFR2) occur in multiple types of cancer¹. However, clinical responses to FGFR inhibitors have remained variable^1–9, emphasizing the need to better understand which FGFR2 alterations are oncogenic and therapeutically targetable. Here we apply transposon-based screening^10,11 and tumour modelling in mice^12,13, and find that the truncation of exon 18 (E18) of Fgfr2 is a potent driver mutation. Human oncogenomic datasets revealed a diverse set of FGFR2 alterations, including rearrangements, E1–E17 partial amplifications, and E18 nonsense and frameshift mutations, each causing the transcription of E18-truncated FGFR2 (FGFR2^ΔE18). Functional in vitro and in vivo examination of a compendium of FGFR2^ΔE18 and full-length variants pinpointed FGFR2-E18 truncation as single-driver alteration in cancer. By contrast, the oncogenic competence of FGFR2 full-length amplifications depended on a distinct landscape of cooperating driver genes. This suggests that genomic alterations that generate stable FGFR2^ΔE18 variants are actionable therapeutic targets, which we confirmed in preclinical mouse and human tumour models, and in a clinical trial. We propose that cancers containing any FGFR2 variant with a truncated E18 should be considered for FGFR-targeted therapies.

Truncation of exon 18 of FGFR2 (FGFR2^ΔE18) is a potent driver mutation in mice and humans, and FGFR-targeted therapy should be considered for patients with cancer expressing stable FGFR2^ΔE18 variants.

Drug-induced epigenomic plasticity reprograms circadian rhythm regulation to drive prostate cancer towards androgen-independence

In prostate cancer, androgen receptor (AR)-targeting agents are very effective in various disease stages. However, therapy resistance inevitably occurs and little is known about how tumor cells adapt to bypass AR suppression. Here, we performed integrative multi-omics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from high-risk prostate cancer patients enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors towards a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 - from inactive chromatin sites towards active cis-regulatory elements that dictate pro-survival signals. Notably, treatment-induced FOXA1 sites were enriched for circadian clock component ARNTL. Post-treatment ARNTL levels associated with poor outcome, and ARNTL knockout strongly decreased prostate cancer cell growth. Our data highlight a remarkable cistromic plasticity of FOXA1 following AR-targeted therapy, and revealed an acquired dependency on circadian regulator ARNTL, a novel candidate therapeutic target.

Genomic analysis defines clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancer

Ductal carcinoma in situ (DCIS) is the most common form of preinvasive breast cancer and, despite treatment, a small fraction (5-10%) of DCIS patients develop subsequent invasive disease. A fundamental biologic question is whether the invasive disease arises from tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial DCIS lesion and paired invasive recurrent tumors in 95 patients together with single-cell DNA sequencing in a subset of cases. Our data show that in 75% of cases the invasive recurrence was clonally related to the initial DCIS, suggesting that tumor cells were not eliminated during the initial treatment. Surprisingly, however, 18% were clonally unrelated to the DCIS, representing new independent lineages and 7% of cases were ambiguous. This knowledge is essential for accurate risk evaluation of DCIS, treatment de-escalation strategies and the identification of predictive biomarkers.

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.

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.

Flexible co-data learning for high-dimensional prediction

Clinical research often focuses on complex traits in which many variables play a role in mechanisms driving, or curing, diseases. Clinical prediction is hard when data is high-dimensional, but additional information, like domain knowledge and previously published studies, may be helpful to improve predictions. Such complementary data, or co-data, provide information on the covariates, such as genomic location or P-values from external studies. We use multiple and various co-data to define possibly overlapping or hierarchically structured groups of covariates. These are then used to estimate adaptive multi-group ridge penalties for generalized linear and Cox models. Available group adaptive methods primarily target for settings with few groups, and therefore likely overfit for non-informative, correlated or many groups, and do not account for known structure on group level. To handle these issues, our method combines empirical Bayes estimation of the hyperparameters with an extra level of flexible shrinkage. This renders a uniquely flexible framework as any type of shrinkage can be used on the group level. We describe various types of co-data and propose suitable forms of hypershrinkage. The method is very versatile, as it allows for integration and weighting of multiple co-data sets, inclusion of unpenalized covariates and posterior variable selection. For three cancer genomics applications we demonstrate improvements compared to other models in terms of performance, variable selection stability and validation.

Contralateral parenchymal enhancement on MRI is associated with tumor proteasome pathway gene expression and overall survival of early ER+/HER2-breast cancer patients

Purpose

To assess whether contralateral parenchymal enhancement (CPE) on MRI is associated with gene expression pathways in ER+/HER2-breast cancer, and if so, whether such pathways are related to survival.

Methods

Preoperative breast MRIs were analyzed of early ER+/HER2-breast cancer patients eligible for breast-conserving surgery included in a prospective observational cohort study (MARGINS). The contralateral parenchyma was segmented and CPE was calculated as the average of the top-10% delayed enhancement. Total tumor RNA sequencing was performed and gene set enrichment analysis was used to reveal gene expression pathways associated with CPE (N = 226) and related to overall survival (OS) and invasive disease-free survival (IDFS) in multivariable survival analysis. The latter was also done for the METABRIC cohort (N = 1355).

Results

CPE was most strongly correlated with proteasome pathways (normalized enrichment statistic = 2.04, false discovery rate = .11). Patients with high CPE showed lower tumor proteasome gene expression. Proteasome gene expression had a hazard ratio (HR) of 1.40 (95% CI = 0.89, 2.16; P = .143) for OS in the MARGINS cohort and 1.53 (95% CI = 1.08, 2.14; P = .017) for IDFS, in METABRIC proteasome gene expression had an HR of 1.09 (95% CI = 1.01, 1.18; P = .020) for OS and 1.10 (95% CI = 1.02, 1.18; P = .012) for IDFS.

Conclusion

CPE was negatively correlated with tumor proteasome gene expression in early ER+/HER2-breast cancer patients. Low tumor proteasome gene expression was associated with improved survival in the METABRIC data.

•

Contralateral parenchymal enhancement on MRI was associated with tumor proteasome gene expression in ER+/HER2-breast cancer.

•

A high contralateral parenchymal enhancement was associated with a low proteasome gene expression in the breast cancer.

•

Low proteasome tumor gene expression was associated with improved survival in an independent patient cohort.

The Prognostic Potential of Human Prostate Cancer-Associated Macrophage Subtypes as Revealed by Single-Cell Transcriptomics

Macrophages in the tumor microenvironment are causally linked with prostate cancer development and progression, yet little is known about their composition in neoplastic human tissue. By performing single cell transcriptomic analysis of human prostate cancer resident macrophages, three distinct populations were identified in the diseased prostate. Unexpectedly, no differences were observed between macrophages isolated from the tumorous and nontumorous portions of the prostatectomy specimens. Markers associated with canonical M1 and M2 macrophage phenotypes were identifiable, however these were not the main factors defining unique subtypes. The genes selectively associated with each macrophage cluster were used to develop a gene signature which was highly associated with both recurrence-free and metastasis-free survival. These results highlight the relevance of tissue-specific macrophage subtypes in the tumor microenvironment for prostate cancer progression and demonstrates the utility of profiling single-cell transcriptomics in human tumor samples as a strategy to design gene classifiers for patient prognostication. IMPLICATIONS: The specific macrophage subtypes present in a diseased human prostate have prognostic value, suggesting that the relative proportions of these populations are related to patient outcome. Understanding the relative contributions of these subtypes will not only inform patient prognostication, but will enable personalized immunotherapeutic strategies to increase beneficial populations or reduce detrimental populations.

Limited evolution of the actionable metastatic cancer genome under therapeutic pressure

Genomic profiling is critical for the identification of treatment options for patients with metastatic cancer, but it remains unclear how frequently this procedure should be repeated during the course of the disease. To address this, we analyzed whole-genome sequencing (WGS) data of 250 biopsy pairs, longitudinally collected over the treatment course of 231 adult patients with a representative variety of metastatic solid malignancies. Within the biopsy interval (median, 6.4 months), patients received one or multiple lines of (mostly) standard-of-care (SOC) treatments, with all major treatment modalities being broadly represented. SOC biomarkers and biomarkers for clinical trial enrollment could be identified in 23% and 72% of biopsies, respectively. For SOC genomic biomarkers, we observed full concordance between the first and the second biopsy in 99% of pairs. Of the 219 biomarkers for clinical trial enrollment that were identified in the first biopsies, we recovered 94% in the follow-up biopsies. Furthermore, a second WGS analysis did not identify additional biomarkers for clinical trial enrollment in 91% of patients. More-frequent genomic evolution was observed when considering specific genes targeted by small-molecule inhibitors or hormonal therapies (21% and 22% of cases, respectively). Together, our data demonstrate that there is limited evolution of the actionable genome of treated metastases. A single WGS analysis of a metastatic biopsy is generally sufficient to identify SOC genomic biomarkers and to identify investigational treatment opportunities.

The Cancer SENESCopedia: A delineation of cancer cell senescence

Cellular senescence is characterized as a stable proliferation arrest that can be triggered by multiple stresses. Most knowledge about senescent cells is obtained from studies in primary cells. However, senescence features may be different in cancer cells, since the pathways that are involved in senescence induction are often deregulated in cancer. We report here a comprehensive analysis of the transcriptome and senolytic responses in a panel of 13 cancer cell lines rendered senescent by two distinct compounds. We show that in cancer cells, the response to senolytic agents and the composition of the senescence-associated secretory phenotype are more influenced by the cell of origin than by the senescence trigger. Using machine learning, we establish the SENCAN gene expression classifier for the detection of senescence in cancer cell samples. The expression profiles and senescence classifier are available as an interactive online Cancer SENESCopedia.

•

Senescent cancer cells respond differently to senolytic ABT-263

•

SASP expression in cancer is heterogeneous and influenced by cell origin

•

The SENCAN classifier detects cancer cell senescence in vitro

•

The Cancer SENESCopedia contains transcriptome data from 37 senescence models

Breast adipocyte size associates with ipsilateral invasive breast cancer risk after ductal carcinoma in situ

Although ductal carcinoma in situ (DCIS) is a non-obligate precursor to ipsilateral invasive breast cancer (iIBC), most DCIS lesions remain indolent. Hence, overdiagnosis and overtreatment of DCIS is a major concern. There is an urgent need for prognostic markers that can distinguish harmless from potentially hazardous DCIS. We hypothesised that features of the breast adipose tissue may be associated with risk of subsequent iIBC. We performed a case-control study nested in a population-based DCIS cohort, consisting of 2658 women diagnosed with primary DCIS between 1989 and 2005, uniformly treated with breast conserving surgery (BCS) alone. We assessed breast adipose features with digital pathology (HALO®, Indica Labs) and related these to iIBC risk in 108 women that developed subsequent iIBC (cases) and 168 women who did not (controls) by conditional logistic regression, accounting for clinicopathological and immunohistochemistry variables. Large breast adipocyte size was significantly associated with iIBC risk (odds ratio (OR) 2.75, 95% confidence interval (95% CI) = 1.25-6.05). High cyclooxygenase (COX)-2 protein expression in the DCIS cells was also associated with subsequent iIBC (OR 3.70 (95% CI = 1.59-8.64). DCIS with both high COX-2 expression and large breast adipocytes was associated with a 12-fold higher risk (OR 12.0, 95% CI = 3.10-46.3, P < 0.001) for subsequent iIBC compared with women with smaller adipocyte size and low COX-2 expression. Large breast adipocytes combined with high COX-2 expression in DCIS is associated with a high risk of subsequent iIBC. Besides COX-2, adipocyte size has the potential to improve clinical management in patients diagnosed with primary DCIS.

Correction to: A role for the unfolded protein response stress sensor ERN1 in regulating the response to MEK inhibitors in KRAS mutant colon cancers

An amendment to this paper has been published and can be accessed via the original article.

Epigenetic profiling demarcates molecular subtypes of muscle-invasive bladder cancer

Muscle-invasive bladder cancer (MIBC) is a heterogeneous disease that often recurs despite aggressive treatment with neoadjuvant chemotherapy and (radical) cystectomy. Basal and luminal molecular subtypes have been identified that are linked to clinical characteristics and have differential sensitivities to chemotherapy. While it has been suggested that epigenetic mechanisms play a role in defining these subtypes, a thorough understanding of the biological mechanisms is lacking. This report details the first genome-wide analysis of histone methylation patterns of human primary bladder tumours by chromatin immunoprecipitations and next-generation sequencing (ChIP-seq). We profiled multiple histone marks: H3K27me3, a marker for repressed genes, and H3K4me1 and H3K4me3, which are indicators of active enhancers and active promoters. Integrated analysis of ChIP-seq data and RNA sequencing revealed that H3K4 mono-methylation demarcates MIBC subtypes, while no association was found for the other two histone modifications in relation to basal and luminal subtypes. Additionally, we identified differentially methylated H3K4me1 peaks in basal and luminal tumour samples, suggesting that active enhancers play a role in defining subtypes. Our study is the first analysis of histone modifications in primary bladder cancer tissue and provides an important resource for the bladder cancer community.

Functional Categorization of BRCA1 Variants of Uncertain Clinical Significance in Homologous Recombination Repair Complementation Assays

PURPOSE

Because BRCA1 is a high-risk breast/ovarian cancer susceptibility gene, BRCA1 sequence variants of uncertain clinical significance (VUS) complicate genetic counseling. As most VUS are rare, reliable classification based on clinical and genetic data is often impossible. However, all pathogenic BRCA1 variants analyzed result in defective homologous recombination DNA repair (HRR). Thus, BRCA1 VUS may be categorized based on their functional impact on this pathway.

EXPERIMENTAL DESIGN

Two hundred thirty-eight BRCA1 VUS-comprising most BRCA1 VUS known in the Netherlands and Belgium-were tested for their ability to complement Brca1-deficient mouse embryonic stem cells in HRR, using cisplatin and olaparib sensitivity assays and a direct repeat GFP (DR-GFP) HRR assay. Assays were validated using 25 known benign and 25 known pathogenic BRCA1 variants. For assessment of pathogenicity by a multifactorial likelihood analysis method, we collected clinical and genetic data for functionally deleterious VUS and VUS occurring in three or more families.

RESULTS

All three assays showed 100% sensitivity and specificity (95% confidence interval, 83%-100%). Out of 238 VUS, 45 showed functional defects, 26 of which were deleterious in all three assays. For 13 of these 26 variants, we could calculate the probability of pathogenicity using clinical and genetic data, resulting in the identification of 7 (likely) pathogenic variants.

CONCLUSIONS

We have functionally categorized 238 BRCA1 VUS using three different HRR-related assays. Classification based on clinical and genetic data alone for a subset of these variants confirmed the high sensitivity and specificity of our functional assays.

PRECISE: a domain adaptation approach to transfer predictors of drug response from pre-clinical models to tumors

Motivation

Cell lines and patient-derived xenografts (PDXs) have been used extensively to understand the molecular underpinnings of cancer. While core biological processes are typically conserved, these models also show important differences compared to human tumors, hampering the translation of findings from pre-clinical models to the human setting. In particular, employing drug response predictors generated on data derived from pre-clinical models to predict patient response remains a challenging task. As very large drug response datasets have been collected for pre-clinical models, and patient drug response data are often lacking, there is an urgent need for methods that efficiently transfer drug response predictors from pre-clinical models to the human setting.

Results

We show that cell lines and PDXs share common characteristics and processes with human tumors. We quantify this similarity and show that a regression model cannot simply be trained on cell lines or PDXs and then applied on tumors. We developed PRECISE, a novel methodology based on domain adaptation that captures the common information shared amongst pre-clinical models and human tumors in a consensus representation. Employing this representation, we train predictors of drug response on pre-clinical data and apply these predictors to stratify human tumors. We show that the resulting domain-invariant predictors show a small reduction in predictive performance in the pre-clinical domain but, importantly, reliably recover known associations between independent biomarkers and their companion drugs on human tumors.

Availability and implementation

PRECISE and the scripts for running our experiments are available on our GitHub page (https://github.com/NKI-CCB/PRECISE).

Supplementary information

Supplementary data are available at Bioinformatics online.

Radiogenomic Analysis of Breast Cancer by Linking MRI Phenotypes with Tumor Gene Expression

Background Better understanding of the molecular biology associated with MRI phenotypes may aid in the diagnosis and treatment of breast cancer. Purpose To discover the associations between MRI phenotypes of breast cancer and their underlying molecular biology derived from gene expression data. Materials and Methods This is a secondary analysis of the Multimodality Analysis and Radiologic Guidance in Breast-Conserving Therapy, or MARGINS, study. MARGINS included patients eligible for breast-conserving therapy between November 2000 and December 2008 for preoperative breast MRI. Tumor RNA was collected for sequencing from surgical specimen. Twenty-one computer-generated MRI features of tumors were condensed into seven MRI factors related to tumor size, shape, initial enhancement, late enhancement, smoothness of enhancement, sharpness, and sharpness variation. These factors were associated with gene expression levels from RNA sequencing by using gene set enrichment analysis. Statistical significance of these associations was evaluated by using a sample permutation test and the false discovery rate. Results Gene expression and MRI data were obtained for 295 patients (mean age, 56 years ± 10.3 [standard deviation]). Larger and more irregular tumors showed increased expression of cell cycle and DNA damage checkpoint genes (false discovery rate <0.25; normalized enrichment statistic [NES], 2.15). Enhancement and sharpness of the tumor margin were associated with expression of ribosomal proteins (false discovery rate <0.25; NES, 1.95). Smoothness of enhancement, tumor size, and tumor shape were associated with expression of genes involved in the extracellular matrix (false discovery rate <0.25; NES, 2.25). Conclusion Breast cancer MRI phenotypes were related to their underlying molecular biology revealed by using RNA sequencing. The association between enhancements and sharpness of the tumor margin with the ribosome suggests that these MRI features may be imaging biomarkers for drugs targeting the ribosome. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Cho in this issue.

Spatial proximity of CD8 T cells to tumor cells as an independent biomarker for response to anti-PD-1 therapy

10038

Background: Only a subset of advanced melanoma patients respond to anti-PD-1 (aPD1) monotherapy. Upfront identification of (non-)responsiveness would help guide first-line treatment decisions, prevent overtreatment and unnecessary risk for toxicities. T cell density and expression of T cell related genes have been associated with response to aPD1, but are imperfect predictors. We investigated whether spatial proximity of CD8 T cells to tumor cells improves upon the predictive value of T cell density alone. Methods: Pretreatment tumor specimens from melanoma patients treated with aPD1 in the Netherlands Cancer Institute were stained for DAPI, SOX10/Melan-A, CD4, CD8, FOXP3 and PD-1 by multiplex immunofluorescence. Sections were imaged on Vectra and analyzed using HALO to optimize marker thresholds and demarcate tumor and stroma. T cell proximity to tumor cells was evaluated as difference in area under the curve between i) a spatial G-function quantifying T cell density around tumor cells in tumor areas and ii) analogous null distributions obtained by random permutation of cell labels. This assessment of co-clustering is independent of cell density and heterogeneity therein and does not reflect repulsion of T cells to stromal/marginal areas. Clinical characteristics, RECIST response and survival were collected from patient records. Associations between T cell density, T cell proximity to Sox10/Melan-A+ tumor cells, other clinical biomarkers (LDH, M stage and WHO performance status) and response were examined in a Bayesian hierarchical logistic regression. Results: Tumor specimens of 98 patients were included, of whom 45 were treated with aPD1 as first-line therapy and 33 had an objective response. CD8 T cell proximity to tumor cells was associated with response in an independent, comparatively strong, and tissue dependent manner (cutaneous tissue: 2.78 [2.45, 3.17], visceral: 2.30 [1.95, 2.72], lymphoid: 2.12 [1.88, 2.40], format: maximal posteriori odds ratio [89% equal-tailed credibility interval]), in a multivariate model correcting for CD8 T cell density (1.74 [1.62, 1.88]), LDH (1.93 [1.72, 2.16]), M stage (0.92 [0.87, 0.98]) and WHO performance status (0.79 [0.72, 0.88]). Our model achieved an area under the ROC curve of 77.7%, whereas an analogous model omitting the proximity variable achieved 73.1%. Conclusions: Our analyses show that spatial proximity of CD8 T cells to tumor cells functions as an independent biomarker for response to aPD1 and suggests that preexisting CD8 T cell tumor reactivity is reflected by this spatial proximity.

Identifying Epistasis in Cancer Genomes: A Delicate Affair

Recent studies of the tumor genome seek to identify cancer pathways as groups of genes in which mutations are epistatic with one another or, specifically, "mutually exclusive." Here, we show that most mutations are mutually exclusive not due to pathway structure but to interactions with disease subtype and tumor mutation load. In particular, many cancer driver genes are mutated preferentially in tumors with few mutations overall, causing mutations in these cancer genes to appear mutually exclusive with numerous others. Researchers should view current epistasis maps with caution until we better understand the multiple cause-and-effect relationships among factors such as tumor subtype, positive selection for mutations, and gross tumor characteristics including mutational signatures and load.

Endonuclease FEN1 Coregulates ERα Activity and Provides a Novel Drug Interface in Tamoxifen-Resistant Breast Cancer

Estrogen receptor α (ERα) is a key transcriptional regulator in the majority of breast cancers. ERα-positive patients are frequently treated with tamoxifen, but resistance is common. In this study, we refined a previously identified 111-gene outcome prediction-classifier, revealing FEN1 as the strongest determining factor in ERα-positive patient prognostication. FEN1 levels were predictive of outcome in tamoxifen-treated patients, and FEN1 played a causal role in ERα-driven cell growth. FEN1 impacted the transcriptional activity of ERα by facilitating coactivator recruitment to the ERα transcriptional complex. FEN1 blockade induced proteasome-mediated degradation of activated ERα, resulting in loss of ERα-driven gene expression and eradicated tumor cell proliferation. Finally, a high-throughput 465,195 compound screen identified a novel FEN1 inhibitor, which effectively blocked ERα function and inhibited proliferation of tamoxifen-resistant cell lines as well as ex vivo-cultured ERα-positive breast tumors. Collectively, these results provide therapeutic proof of principle for FEN1 blockade in tamoxifen-resistant breast cancer. SIGNIFICANCE: These findings show that pharmacologic inhibition of FEN1, which is predictive of outcome in tamoxifen-treated patients, effectively blocks ERα function and inhibits proliferation of tamoxifen-resistant tumor cells.

Response of metastatic mouse invasive lobular carcinoma to mTOR inhibition is partly mediated by the adaptive immune system

Effective treatment of invasive lobular carcinoma (ILC) of the breast is hampered by late detection, invasive growth, distant metastasis, and poor response to chemotherapy. Phosphoinositide 3-kinase (PI3K) signaling, one of the major druggable oncogenic signaling networks, is frequently activated in ILC. We investigated treatment response and resistance to AZD8055, an inhibitor of mammalian target of rapamycin (mTOR), in the K14-cre;Cdh1^Flox/Flox;Trp53^Flox/Flox (KEP) mouse model of metastatic ILC. Inhibition of mTOR signaling blocked the growth of primary KEP tumors as well as the progression of metastatic disease. However, primary tumors and distant metastases eventually acquired resistance after long-term AZD8055 treatment, despite continued effective suppression of mTOR signaling in cancer cells. Interestingly, therapeutic responses were associated with increased expression of genes related to antigen presentation. Consistent with this observation, increased numbers of tumor-infiltrating major histocompatibility complex class II-positive (MHCII+) immune cells were observed in treatment-responsive KEP tumors. Acquisition of treatment resistance was associated with loss of MHCII+ cells and reduced expression of genes related to the adaptive immune system. The therapeutic efficacy of mTOR inhibition was reduced in Rag1^−/- mice lacking mature T and B lymphocytes, compared to immunocompetent mice. Furthermore, therapy responsiveness could be partially rescued by transplanting AZD8055-resistant KEP tumors into treatment-naïve immunocompetent hosts. Collectively, these data indicate that the PI3K signaling pathway is an attractive therapeutic target in invasive lobular carcinoma, and that part of the therapeutic effect of mTOR inhibition is mediated by the adaptive immune system.

Androgen modulation of XBP1 is functionally driving part of the AR transcriptional program

Prostate cancer development and progression is largely dependent on androgen receptor (AR) signaling. AR is a hormone-dependent transcription factor, which binds to thousands of sites throughout the human genome to regulate expression of directly responsive genes, including pro-survival genes that enable tumor cells to cope with increased cellular stress. ERN1 and XBP1 - two key players of the unfolded protein response (UPR) - are among such stress-associated genes. Here, we show that XBP1 levels in primary prostate cancer are associated with biochemical recurrence in five independent cohorts. Patients who received AR-targeted therapies had significantly lower XBP1 expression, whereas expression of the active form of XBP1 (XBP1s) was elevated. In vitro results show that AR-induced ERN1 expression led to increased XBP1s mRNA and protein levels. Furthermore, ChIP-seq analysis revealed that XBP1s binds enhancers upon stress stimuli regulating genes involved in UPR processes, eIF2 signaling and protein ubiquitination. We further demonstrate genomic overlap of AR- and XBP1s-binding sites, suggesting genomic conversion of the two signaling cascades. Transcriptomic effects of XBP1 were further studied by knockdown experiments, which lead to decreased expression of androgen-responsive genes and UPR genes. These results suggest a two-step mechanism of gene regulation, which involves androgen-induced expression of ERN1, thereby enhancing XBP1 splicing and transcriptional activity. This signaling cascade may prepare the cells for the increased protein folding, mRNA decay and translation that accompanies AR-regulated tumor cell proliferation.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

RUNX2/CBFB modulates the response to MEK inhibitors through activation of receptor tyrosine kinases in KRAS-mutant colorectal cancer

Intrinsic and acquired resistances are major hurdles preventing the effective use of MEK inhibitors for treatment of colorectal cancer (CRC). Some 35-45% of colorectal cancers are KRAS-mutant and their treatment remains challenging as these cancers are refractory to MEK inhibitor treatment, because of feedback activation of receptor tyrosine kinases (RTKs). We reported previously that loss of ERN1 sensitizes a subset of KRAS-mutant colon cancer cells to MEK inhibition. Here we show that the loss of RUNX2 or its cofactor CBFB can confer MEK inhibitor resistance in CRC cells. Mechanistically, we find that cells with genetically ablated RUNX2 or CBFB activate multiple RTKs, which coincides with high SHP2 phosphatase activity, a phosphatase that relays signals from the cell membrane to downstream pathways governing growth and proliferation. Moreover, we show that high activity of SHP2 is causal to loss of RUNX2-induced MEK inhibitor resistance, as a small molecule SHP2 inhibitor reinstates sensitivity to MEK inhibitor in RUNX2 knockout cells. Our results reveal an unexpected role for loss of RUNX2/CBFB in regulating RTK activity in colon cancer, resulting in reduced sensitivity to MEK inhibitors.

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.

Modelling signalling networks from perturbation data

Motivation

Intracellular signalling is realized by complex signalling networks, which are almost impossible to understand without network models, especially if feedbacks are involved. Modular Response Analysis (MRA) is a convenient modelling method to study signalling networks in various contexts.

Results

We developed the software package STASNet (STeady-STate Analysis of Signalling Networks) that provides an augmented and extended version of MRA suited to model signalling networks from incomplete perturbation schemes and multi-perturbation data. Using data from the Dialogue on Reverse Engineering Assessment and Methods challenge, we show that predictions from STASNet models are among the top-performing methods. We applied the method to study the effect of SHP2, a protein that has been implicated in resistance to targeted therapy in colon cancer, using a novel dataset from the colon cancer cell line Widr and a SHP2-depleted derivative. We find that SHP2 is required for mitogen-activated protein kinase signalling, whereas AKT signalling only partially depends on SHP2.

Availability and implementation

An R-package is available at https://github.com/molsysbio/STASNet.

Supplementary information

Supplementary data are available at Bioinformatics online.

iTOP: inferring the topology of omics data

Motivation

In biology, we are often faced with multiple datasets recorded on the same set of objects, such as multi-omics and phenotypic data of the same tumors. These datasets are typically not independent from each other. For example, methylation may influence gene expression, which may, in turn, influence drug response. Such relationships can strongly affect analyses performed on the data, as we have previously shown for the identification of biomarkers of drug response. Therefore, it is important to be able to chart the relationships between datasets.

Results

We present iTOP, a methodology to infer a topology of relationships between datasets. We base this methodology on the RV coefficient, a measure of matrix correlation, which can be used to determine how much information is shared between two datasets. We extended the RV coefficient for partial matrix correlations, which allows the use of graph reconstruction algorithms, such as the PC algorithm, to infer the topologies. In addition, since multi-omics data often contain binary data (e.g. mutations), we also extended the RV coefficient for binary data. Applying iTOP to pharmacogenomics data, we found that gene expression acts as a mediator between most other datasets and drug response: only proteomics clearly shares information with drug response that is not present in gene expression. Based on this result, we used TANDEM, a method for drug response prediction, to identify which variables predictive of drug response were distinct to either gene expression or proteomics.

Availability and implementation

An implementation of our methodology is available in the R package iTOP on CRAN. Additionally, an R Markdown document with code to reproduce all figures is provided as Supplementary Material.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comparative Network Reconstruction using mixed integer programming

Motivation

Signal-transduction networks are often aberrated in cancer cells, and new anti-cancer drugs that specifically target oncogenes involved in signaling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. A quantitative understanding of these signaling networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells can help in addressing this problem.

Results

Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on possibly incomplete perturbation data, and to identify which edges differ quantitatively between two or more signaling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant, PTPN11 KO cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re-establishing wild-type MAPK signaling, possibly through an alternative RAF-isoform.

Availability and implementation

CNR is available as a python module at https://github.com/NKI-CCB/cnr. Additionally, code to reproduce all figures is available at https://github.com/NKI-CCB/CNR-analyses.

Supplementary information

Supplementary data are available at Bioinformatics online.

Loss of p53 triggers WNT-dependent systemic inflammation to drive breast cancer metastasis

Cancer-associated systemic inflammation is strongly linked to poor disease outcome in patients with cancer^1,2. For most human epithelial tumour types, high systemic neutrophil-to-lymphocyte ratios are associated with poor overall survival³, and experimental studies have demonstrated a causal relationship between neutrophils and metastasis^4,5. However, the cancer-cell-intrinsic mechanisms that dictate the substantial heterogeneity in systemic neutrophilic inflammation between tumour-bearing hosts are largely unresolved. Here, using a panel of 16 distinct genetically engineered mouse models for breast cancer, we uncover a role for cancer-cell-intrinsic p53 as a key regulator of pro-metastatic neutrophils. Mechanistically, loss of p53 in cancer cells induced the secretion of WNT ligands that stimulate tumour-associated macrophages to produce IL-1β, thus driving systemic inflammation. Pharmacological and genetic blockade of WNT secretion in p53-null cancer cells reverses macrophage production of IL-1β and subsequent neutrophilic inflammation, resulting in reduced metastasis formation. Collectively, we demonstrate a mechanistic link between the loss of p53 in cancer cells, secretion of WNT ligands and systemic neutrophilia that potentiates metastatic progression. These insights illustrate the importance of the genetic makeup of breast tumours in dictating pro-metastatic systemic inflammation, and set the stage for personalized immune intervention strategies for patients with cancer.

Principal component analysis of binary genomics data

Motivation

Genome-wide measurements of genetic and epigenetic alterations are generating more and more high-dimensional binary data. The special mathematical characteristics of binary data make the direct use of the classical principal component analysis (PCA) model to explore low-dimensional structures less obvious. Although there are several PCA alternatives for binary data in the psychometric, data analysis and machine learning literature, they are not well known to the bioinformatics community. Results: In this article, we introduce the motivation and rationale of some parametric and nonparametric versions of PCA specifically geared for binary data. Using both realistic simulations of binary data as well as mutation, CNA and methylation data of the Genomic Determinants of Sensitivity in Cancer 1000 (GDSC1000), the methods were explored for their performance with respect to finding the correct number of components, overfit, finding back the correct low-dimensional structure, variable importance, etc. The results show that if a low-dimensional structure exists in the data, that most of the methods can find it. When assuming a probabilistic generating process is underlying the data, we recommend to use the parametric logistic PCA model, while when such an assumption is not valid and the data are considered as given, the nonparametric Gifi model is recommended.

Availability

The codes to reproduce the results in this article are available at the homepage of the Biosystems Data Analysis group (www.bdagroup.nl).

Basket of baskets (BoB): A modular, open label, phase II, multicenter study to evaluate targeted agents in molecularly selected populations with advanced solid tumors

TPS3151

Background: Basket trials with targeted agents can show high response rates for tumors with specific molecular profiles, granting extension of the label of some drugs. In other cases, study results were disappointing, likely due to the rarity of molecular alterations, limits in trial design and the difficulties in applying molecular tumor profiling in the clinical setting. Methods: Basket of Baskets (BoB), NCT03767075, aims to bridge the gap between Academic Genomics and clinical applications (ready-to market multi-marker Companion Diagnostics) by providing a sustainable and adaptable (to new technologies, markers, and therapeutic agents) platform for co-development of drug/companion diagnostic. BoB is a novel platform trial from Cancer Core Europe, a recently established sustainable European network for innovative cancer research. This protocol has two parts: (1) Part A includes a molecular profiling program for subjects with advanced solid tumors (iPROFILER), a variant annotation tool, and a molecular tumor board to select the most appropriate treatment. It also enables testing/developing companion diagnostics linked with the therapeutic part (part B). (2) Part B includes iBASKET, a modular multi-arm basket trial for subjects with tumors harboing selected molecular alterations. Each module is focused on a certain molecular pathway or on certain molecular alterations that may confer sensitivity to the study drug or study drug combination evaluated in that module/arm. The current version of iBASKET (Module 1- Atezolizumab in genomically-selected patients) is open for enrollment for patients with advanced neoplasms bearing one of the following alterations: Arm 1A: BRCA1 or BRCA2; Arm 1B: MLH1, MSH2, MSH6, PMS2; Arm 1C: POLE, POLD1 mutations; Arm 1D: hypermutated tumors; Arm 1E: other mutations in DNA-repair genes; Arm 1F: PDL1 gene amplification. All patients enrolled in Module 1 will receive single-agent atezolizumab. New Modules for genomically selected populations can be added through amendments. Our final aim is to achieve drug repurposing of treatments, co-develop multi-marker companion diagnostics and a large database of knowledge in Precision Medicine. Clinical trial information: NCT03767075.

Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency

Understanding the impact of guide RNA (gRNA) and genomic locus on CRISPR-Cas9 activity is crucial to design effective gene editing assays. However, it is challenging to profile Cas9 activity in the endogenous cellular environment. Here we leverage our TRIP technology to integrate ~ 1k barcoded reporter genes in the genomes of mouse embryonic stem cells. We target the integrated reporters (IRs) using RNA-guided Cas9 and characterize induced mutations by sequencing. We report that gRNA-sequence and IR locus explain most variation in mutation efficiency. Predominant insertions of a gRNA-specific nucleotide are consistent with template-dependent repair of staggered DNA ends with 1-bp 5' overhangs. We confirm that such staggered ends are induced by Cas9 in mouse pre-B cells. To explain observed insertions, we propose a model generating primarily blunt and occasionally staggered DNA ends. Mutation patterns indicate that gRNA-sequence controls the fraction of staggered ends, which could be used to optimize Cas9-based insertion efficiency.

Inhibition of the Replication Stress Response Is a Synthetic Vulnerability in SCLC That Acts Synergistically in Combination with Cisplatin

Small cell lung cancer (SCLC) is generally regarded as very difficult to treat, mostly due to the development of metastases early in the disease and a quick relapse with resistant disease. SCLC patients initially show a good response to treatment with the DNA damaging agents cisplatin and etoposide. This is, however, quickly followed by the development of resistant disease, which urges the development of novel therapies for this type of cancer. In this study, we set out to compile a comprehensive overview of the vulnerabilities of SCLC. A functional genome-wide screen where all individual genes were knocked out was performed to identify novel vulnerabilities of SCLC. By analysis of the knockouts that were lethal to these cancer cells, we identified several processes to be synthetic vulnerabilities in SCLC. We were able to validate the vulnerability to inhibition of the replication stress response machinery by use of Chk1 and ATR inhibitors. Strikingly, SCLC cells were more sensitive to these inhibitors than nontransformed cells. In addition, these inhibitors work synergistically with either etoposide and cisplatin, where the interaction is largest with the latter. ATR inhibition by VE-822 treatment in combination with cisplatin also outperforms the combination of cisplatin with etoposide in vivo Altogether, our study uncovered a critical dependence of SCLC on the replication stress response and urges the validation of ATR inhibitors in combination with cisplatin in a clinical setting.