A Proof-of-Concept Study of Sequential Treatment with the HDAC Inhibitor Vorinostat following BRAF and MEK Inhibitors in BRAFV600-Mutated Melanoma.

PURPOSE: The development of resistance limits the clinical benefit of BRAF and MEK inhibitors (BRAFi/MEKi) in BRAFV600-mutated melanoma. It has been shown that short-term treatment (14 days) with vorinostat was able to initiate apoptosis of resistant tumor cells. We aimed to assess the antitumor activity of sequential treatment with vorinostat following BRAFi/MEKi in patients with BRAFV600-mutated melanoma who progressed after initial response to BRAFi/MEKi. PATIENTS AND METHODS: Patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma were treated with vorinostat 360 mg once daily for 14 days followed by BRAFi/MEKi. The primary endpoint was an objective response rate of progressive lesions of at least 30% according to Response Evaluation Criteria in Solid Tumors 1.1. Secondary endpoints included progression-free survival, overall survival, safety, pharmacokinetics of vorinostat, and translational molecular analyses using ctDNA and tumor biopsies. RESULTS: Of the 26 patients with progressive BRAFi/MEKi-resistant BRAFV600-mutated melanoma receiving treatment with vorinostat, 22 patients were evaluable for response. The objective response rate was 9%, with one complete response for 31.2 months and one partial response for 14.9 months. Median progression-free survival and overall survival were 1.4 and 5.4 months, respectively. Common adverse events were fatigue (23%) and nausea (19%). ctDNA analysis showed emerging secondary mutations in NRAS and MEK in eight patients at the time of BRAFi/MEKi resistance. Elimination of these mutations by vorinostat treatment was observed in three patients. CONCLUSIONS: Intermittent treatment with vorinostat in patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma is well tolerated. Although the primary endpoint of this study was not met, durable antitumor responses were observed in a minority of patients (9%).

Cancer cell genetics shaping of the tumor microenvironment reveals myeloid cell-centric exploitable vulnerabilities in hepatocellular carcinoma.

Myeloid cells are abundant and plastic immune cell subsets in the liver, to which pro-tumorigenic, inflammatory and immunosuppressive roles have been assigned in the course of tumorigenesis. Yet several aspects underlying their dynamic alterations in hepatocellular carcinoma (HCC) progression remain elusive, including the impact of distinct genetic mutations in shaping a cancer-permissive tumor microenvironment (TME). Here, in newly generated, clinically-relevant somatic female HCC mouse models, we identify cancer genetics' specific and stage-dependent alterations of the liver TME associated with distinct histopathological and malignant HCC features. Mitogen-activated protein kinase (MAPK)-activated, NrasG12D-driven tumors exhibit a mixed phenotype of prominent inflammation and immunosuppression in a T cell-excluded TME. Mechanistically, we report a NrasG12D cancer cell-driven, MEK-ERK1/2-SP1-dependent GM-CSF secretion enabling the accumulation of immunosuppressive and proinflammatory monocyte-derived Ly6Clow cells. GM-CSF blockade curbs the accumulation of these cells, reduces inflammation, induces cancer cell death and prolongs animal survival. Furthermore, GM-CSF neutralization synergizes with a vascular endothelial growth factor (VEGF) inhibitor to restrain HCC outgrowth. These findings underscore the profound alterations of the myeloid TME consequential to MAPK pathway activation intensity and the potential of GM-CSF inhibition as a myeloid-centric therapy tailored to subsets of HCC patients.

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.

Long-term outcomes of young, node-negative, chemotherapy-naïve, triple-negative breast cancer patients according to BRCA1 status.

BACKGROUND: Due to the abundant usage of chemotherapy in young triple-negative breast cancer (TNBC) patients, the unbiased prognostic value of BRCA1-related biomarkers in this population remains unclear. In addition, whether BRCA1-related biomarkers modify the well-established prognostic value of stromal tumor-infiltrating lymphocytes (sTILs) is unknown. This study aimed to compare the outcomes of young, node-negative, chemotherapy-naïve TNBC patients according to BRCA1 status, taking sTILs into account. METHODS: We included 485 Dutch women diagnosed with node-negative TNBC under age 40 between 1989 and 2000. During this period, these women were considered low-risk and did not receive chemotherapy. BRCA1 status, including pathogenic germline BRCA1 mutation (gBRCA1m), somatic BRCA1 mutation (sBRCA1m), and tumor BRCA1 promoter methylation (BRCA1-PM), was assessed using DNA from formalin-fixed paraffin-embedded tissue. sTILs were assessed according to the international guideline. Patients' outcomes were compared using Cox regression and competing risk models. RESULTS: Among the 399 patients with BRCA1 status, 26.3% had a gBRCA1m, 5.3% had a sBRCA1m, 36.6% had tumor BRCA1-PM, and 31.8% had BRCA1-non-altered tumors. Compared to BRCA1-non-alteration, gBRCA1m was associated with worse overall survival (OS) from the fourth year after diagnosis (adjusted HR, 2.11; 95% CI, 1.18-3.75), and this association attenuated after adjustment for second primary tumors. Every 10% sTIL increment was associated with 16% higher OS (adjusted HR, 0.84; 95% CI, 0.78-0.90) in gBRCA1m, sBRCA1m, or BRCA1-non-altered patients and 31% higher OS in tumor BRCA1-PM patients. Among the 66 patients with tumor BRCA1-PM and ≥ 50% sTILs, we observed excellent 15-year OS (97.0%; 95% CI, 92.9-100%). Conversely, among the 61 patients with gBRCA1m and < 50% sTILs, we observed poor 15-year OS (50.8%; 95% CI, 39.7-65.0%). Furthermore, gBRCA1m was associated with higher (adjusted subdistribution HR, 4.04; 95% CI, 2.29-7.13) and tumor BRCA1-PM with lower (adjusted subdistribution HR, 0.42; 95% CI, 0.19-0.95) incidence of second primary tumors, compared to BRCA1-non-alteration. CONCLUSIONS: Although both gBRCA1m and tumor BRCA1-PM alter BRCA1 gene transcription, they are associated with different outcomes in young, node-negative, chemotherapy-naïve TNBC patients. By combining sTILs and BRCA1 status for risk classification, we were able to identify potential subgroups in this population to intensify and optimize adjuvant treatment.

Locus-specific proteome decoding reveals Fpt1 as a chromatin-associated negative regulator of RNA polymerase III assembly.

Transcription of tRNA genes by RNA polymerase III (RNAPIII) is tuned by signaling cascades. The emerging notion of differential tRNA gene regulation implies the existence of additional regulatory mechanisms. However, tRNA gene-specific regulators have not been described. Decoding the local chromatin proteome of a native tRNA gene in yeast revealed reprogramming of the RNAPIII transcription machinery upon nutrient perturbation. Among the dynamic proteins, we identified Fpt1, a protein of unknown function that uniquely occupied RNAPIII-regulated genes. Fpt1 binding at tRNA genes correlated with the efficiency of RNAPIII eviction upon nutrient perturbation and required the transcription factors TFIIIB and TFIIIC but not RNAPIII. In the absence of Fpt1, eviction of RNAPIII was reduced, and the shutdown of ribosome biogenesis genes was impaired upon nutrient perturbation. Our findings provide support for a chromatin-associated mechanism required for RNAPIII eviction from tRNA genes and tuning the physiological response to changing metabolic demands.

Comparison of human recombinant protein coatings and fibroblast-ECM to Matrigel for induced pluripotent stem cell culture and renal podocyte differentiation.

Human induced pluripotent stem cells (hiPSCs) offer great opportunities within the 3R framework. In the field of toxicology, they may contribute greatly to the reduction and eventually replacement of animal models. However, culturing hiPSCs as well as differentiation of hiPSCs into target cells that are used for toxicity testing depend on the presence of extracellular matrix (ECM) coating the growth surface. The most widely used ECM is MatrigelR, an animal product that is derived from mouse sarcoma. Drawbacks of Matrigel are widely recognized and include batch-to batch variations, use of animal rather than human material, and ethical concerns about its production. While alternative coatings exist, higher cost and limited characterizations may hinder their broader uptake by the scientific community. Here, we report an extensive comparison of three commercially available human ECM coatings, vitronectin, laminin-511, and laminin-521, to Matrigel in three different hiPSC lines in long-term culture (≥ 9 passages). Characterization included expression of pluripotent markers in a genome-wide transcriptomics study (TempO-Seq), capacity to differentiate into embryoid bodies, and karyotype stability assessed by analyzing copy number variations by shallow DNA sequencing. Furthermore, a low-cost, decellularized ECM produced by human neonatal dermal fibroblasts was tested. In addition, all alternative coatings were tested for hiPSC differentiation into renal podocyte-like cells in a genome-wide transcriptomics screen. Our results show that all tested coatings were highly comparable to animal-derived Matrigel for both hiPSC maintenance and differentiation into renal podocyte-like cells. Furthermore, decellularized fibroblast-ECM could be a novel, attractive low-cost coating material.

Drug-Induced Epigenomic Plasticity Reprograms Circadian Rhythm Regulation to Drive Prostate Cancer toward 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 multiomics analyses on tissues isolated before and after 3 months of AR-targeting enzalutamide monotherapy from patients with high-risk prostate cancer enrolled in a neoadjuvant clinical trial. Transcriptomic analyses demonstrated that AR inhibition drove tumors toward a neuroendocrine-like disease state. Additionally, epigenomic profiling revealed massive enzalutamide-induced reprogramming of pioneer factor FOXA1 from inactive chromatin sites toward active cis-regulatory elements that dictate prosurvival signals. Notably, treatment-induced FOXA1 sites were enriched for the circadian clock component ARNTL. Posttreatment ARNTL levels were associated with patients' clinical outcomes, 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 the circadian regulator ARNTL, a novel candidate therapeutic target. SIGNIFICANCE: Understanding how prostate cancers adapt to AR-targeted interventions is critical for identifying novel drug targets to improve the clinical management of treatment-resistant disease. Our study revealed an enzalutamide-induced epigenomic plasticity toward prosurvival signaling and uncovered the circadian regulator ARNTL as an acquired vulnerability after AR inhibition, presenting a novel lead for therapeutic development. See related commentary by Zhang et al., p. 2017. This article is highlighted in the In This Issue feature, p. 2007.

Adjuvant capecitabine-containing chemotherapy benefit and homologous recombination deficiency in early-stage triple-negative breast cancer patients.

BACKGROUND: The addition of adjuvant capecitabine to standard chemotherapy of early-stage triple-negative breast cancer (TNBC) patients has improved survival in a few randomised trials and in meta-analyses. However, many patients did not benefit. We evaluated the BRCA1-like DNA copy number signature, indicative of homologous recombination deficiency, as a predictive biomarker for capecitabine benefit in the TNBC subgroup of the FinXX trial. METHODS: Early-stage TNBC patients were randomised between adjuvant capecitabine-containing (TX + CEX: capecitabine-docetaxel, followed by cyclophosphamide-epirubicin-capecitabine) and conventional chemotherapy (T + CEF: docetaxel, followed by cyclophosphamide-epirubicin-fluorouracil). Tumour BRCA1-like status was determined on low-coverage, whole genome next-generation sequencing data using an established DNA comparative genomic hybridisation algorithm. RESULTS: For 129/202 (63.9%) patients the BRCA1-like status could be determined, mostly due to lack of tissue. During a median follow-up of 10.7 years, 35 recurrences and 32 deaths occurred. Addition of capecitabine appears to improve recurrence-free survival more among 61 (47.3%) patients with non-BRCA1-like tumours (HR 0.23, 95% CI 0.08-0.70) compared to 68 (52.7%) patients with BRCA1-like tumours (HR 0.66, 95% CI 0.24-1.81) (P-interaction = 0.17). CONCLUSION: Based on our data, patients with non-BRCA1-like TNBC appear to benefit from the addition of capecitabine to adjuvant chemotherapy. Patients with BRCA1-like TNBC may also benefit. Additional research is needed to define the subgroup within BRCA1-like TNBC patients who may not benefit from adjuvant capecitabine.

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.

Predictive Immune-Checkpoint Blockade Classifiers Identify Tumors Responding to Inhibition of PD-1 and/or CTLA-4.

PURPOSE: Combining anti-PD-1 + anti-CTLA-4 immune-checkpoint blockade (ICB) shows improved patient benefit, but it is associated with severe immune-related adverse events and exceedingly high cost. Therefore, there is a dire need to predict which patients respond to monotherapy and which require combination ICB treatment. EXPERIMENTAL DESIGN: In patient-derived melanoma xenografts (PDX), human tumor microenvironment (TME) cells were swiftly replaced by murine cells upon transplantation. Using our XenofilteR deconvolution algorithm we curated human tumor cell RNA reads, which were subsequently subtracted in silico from bulk (tumor cell + TME) patients' melanoma RNA. This produced a purely tumor cell-intrinsic signature ("InTumor") and a signature comprising tumor cell-extrinsic RNA reads ("ExTumor"). RESULTS: We show that whereas the InTumor signature predicts response to anti-PD-1, the ExTumor predicts anti-CTLA-4 benefit. In PDX, InTumorLO, but not InTumorHI, tumors are effectively eliminated by cytotoxic T cells. When used in conjunction, the InTumor and ExTumor signatures identify not only patients who have a substantially higher chance of responding to combination treatment than to either monotherapy, but also those who are likely to benefit little from anti-CTLA-4 on top of anti-PD-1. CONCLUSIONS: These signatures may be exploited to distinguish melanoma patients who need combination ICB blockade from those who likely benefit from either monotherapy.

Kinome capture sequencing of high-grade serous ovarian carcinoma reveals novel mutations in the JAK3 gene.

High-grade serous ovarian carcinoma (HGSOC) remains the deadliest form of epithelial ovarian cancer and despite major efforts little improvement in overall survival has been achieved. Identification of recurring "driver" genetic lesions has the potential to enable design of novel therapies for cancer. Here, we report on a study to find such new therapeutic targets for HGSOC using exome-capture sequencing approach targeting all kinase genes in 127 patient samples. Consistent with previous reports, the most frequently mutated gene was TP53 (97% mutation frequency) followed by BRCA1 (10% mutation frequency). The average mutation frequency of the kinase genes mutated from our panel was 1.5%. Intriguingly, after BRCA1, JAK3 was the most frequently mutated gene (4% mutation frequency). We tested the transforming properties of JAK3 mutants using the Ba/F3 cell-based in vitro functional assay and identified a novel gain-of-function mutation in the kinase domain of JAK3 (p.T1022I). Importantly, p.T1022I JAK3 mutants displayed higher sensitivity to the JAK3-selective inhibitor Tofacitinib compared to controls. For independent validation, we re-sequenced the entire JAK3 coding sequence using tagged amplicon sequencing (TAm-Seq) in 463 HGSOCs resulting in an overall somatic mutation frequency of 1%. TAm-Seq screening of CDK12 in the same population revealed a 7% mutation frequency. Our data confirms that the frequency of mutations in kinase genes in HGSOC is low and provides accurate estimates for the frequency of JAK3 and CDK12 mutations in a large well characterized cohort. Although p.T1022I JAK3 mutations are rare, our functional validation shows that if detected they should be considered as potentially actionable for therapy. The observation of CDK12 mutations in 7% of HGSOC cases provides a strong rationale for routine somatic testing, although more functional and clinical characterization is required to understand which nonsynonymous mutations alterations are associated with homologous recombination deficiency.

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.

Publisher Correction: Immune induction strategies in metastatic triple-negative breast cancer to enhance the sensitivity to PD-1 blockade: the TONIC trial.

In the version of this article originally published, there was an error in Fig. 3j. A label on the heatmap read "TGF-α signaling via NF-κB". It should have read "TNF-α signaling via NF-κB". The error has been corrected in the HTML and PDF versions of this article.

Immune induction strategies in metastatic triple-negative breast cancer to enhance the sensitivity to PD-1 blockade: the TONIC trial.

The efficacy of programmed cell death protein 1 (PD-1) blockade in metastatic triple-negative breast cancer (TNBC) is low1-5, highlighting a need for strategies that render the tumor microenvironment more sensitive to PD-1 blockade. Preclinical research has suggested immunomodulatory properties for chemotherapy and irradiation6-13. In the first stage of this adaptive, non-comparative phase 2 trial, 67 patients with metastatic TNBC were randomized to nivolumab (1) without induction or with 2-week low-dose induction, or with (2) irradiation (3 × 8 Gy), (3) cyclophosphamide, (4) cisplatin or (5) doxorubicin, all followed by nivolumab. In the overall cohort, the objective response rate (ORR; iRECIST14) was 20%. The majority of responses were observed in the cisplatin (ORR 23%) and doxorubicin (ORR 35%) cohorts. After doxorubicin and cisplatin induction, we detected an upregulation of immune-related genes involved in PD-1-PD-L1 (programmed death ligand 1) and T cell cytotoxicity pathways. This was further supported by enrichment among upregulated genes related to inflammation, JAK-STAT and TNF-α signaling after doxorubicin. Together, the clinical and translational data of this study indicate that short-term doxorubicin and cisplatin may induce a more favorable tumor microenvironment and increase the likelihood of response to PD-1 blockade in TNBC. These data warrant confirmation in TNBC and exploration of induction treatments prior to PD-1 blockade in other cancer types.

Low and variable tumor reactivity of the intratumoral TCR repertoire in human cancers.

Infiltration of human cancers by T cells is generally interpreted as a sign of immune recognition, and there is a growing effort to reactivate dysfunctional T cells at such tumor sites1. However, these efforts only have value if the intratumoral T cell receptor (TCR) repertoire of such cells is intrinsically tumor reactive, and this has not been established in an unbiased manner for most human cancers. To address this issue, we analyzed the intrinsic tumor reactivity of the intratumoral TCR repertoire of CD8+ T cells in ovarian and colorectal cancer-two tumor types for which T cell infiltrates form a positive prognostic marker2,3. Data obtained demonstrate that a capacity to recognize autologous tumor is limited to approximately 10% of intratumoral CD8+ T cells. Furthermore, in two of four patient samples tested, no tumor-reactive TCRs were identified, despite infiltration of their tumors by T cells. These data indicate that the intrinsic capacity of intratumoral T cells to recognize adjacent tumor tissue can be rare and variable, and suggest that clinical efforts to reactivate intratumoral T cells will benefit from approaches that simultaneously increase the quality of the intratumoral TCR repertoire.

Role of variant allele fraction and rare SNP filtering to improve cellular DNA repair endpoint association.

BACKGROUND: Large cancer genome studies continue to reveal new players in treatment response and tumorigenesis. The discrimination of functional alterations from the abundance of passenger genetic alterations still poses challenges and determines DNA sequence variant selection procedures. Here we evaluate variant selection strategies that select homozygous variants and rare SNPs and assess its value in detecting tumor cells with DNA repair defects. METHODS: To this end we employed a panel of 29 patient-derived head and neck squamous cell carcinoma (HNSCC) cell lines, of which a subset harbors DNA repair defects. Mitomycin C (MMC) sensitivity was used as functional endpoint of DNA crosslink repair deficiency. 556 genes including the Fanconi anemia (FA) and homologous recombination (HR) genes, whose products strongly determine MMC response, were capture-sequenced. RESULTS: We show a strong association between MMC sensitivity, thus loss of DNA repair function, and the presence of homozygous and rare SNPs in the relevant FA/HR genes. Excluding such selection criteria impedes the discrimination of crosslink repair status by mutation analysis. Applied to all KEGG pathways, we find that the association with MMC sensitivity is strongest in the KEGG FA pathway, therefore also demonstrating the value of such selection strategies for exploratory analyses. Variant analyses in 56 clinical samples demonstrate that homozygous variants occur more frequently in tumor suppressor genes than oncogenes further supporting the role of a homozygosity criterion to improve gene function association or tumor suppressor gene identification studies. CONCLUSION: Together our data show that the detection of relevant genes or of repair pathway defected tumor cells can be improved by the consideration of allele zygosity and SNP allele frequencies.

XenofilteR: computational deconvolution of mouse and human reads in tumor xenograft sequence data.

BACKGROUND: Mouse xenografts from (patient-derived) tumors (PDX) or tumor cell lines are widely used as models to study various biological and preclinical aspects of cancer. However, analyses of their RNA and DNA profiles are challenging, because they comprise reads not only from the grafted human cancer but also from the murine host. The reads of murine origin result in false positives in mutation analysis of DNA samples and obscure gene expression levels when sequencing RNA. However, currently available algorithms are limited and improvements in accuracy and ease of use are necessary. RESULTS: We developed the R-package XenofilteR, which separates mouse from human sequence reads based on the edit-distance between a sequence read and reference genome. To assess the accuracy of XenofilteR, we generated sequence data by in silico mixing of mouse and human DNA sequence data. These analyses revealed that XenofilteR removes > 99.9% of sequence reads of mouse origin while retaining human sequences. This allowed for mutation analysis of xenograft samples with accurate variant allele frequencies, and retrieved all non-synonymous somatic tumor mutations. CONCLUSIONS: XenofilteR accurately dissects RNA and DNA sequences from mouse and human origin, thereby outperforming currently available tools. XenofilteR is open source and available at https://github.com/PeeperLab/XenofilteR .

Integrative Modeling Identifies Key Determinants of Inhibitor Sensitivity in Breast Cancer Cell Lines.

Cancer cell lines differ greatly in their sensitivity to anticancer drugs as a result of different oncogenic drivers and drug resistance mechanisms operating in each cell line. Although many of these mechanisms have been discovered, it remains a challenge to understand how they interact to render an individual cell line sensitive or resistant to a particular drug. To better understand this variability, we profiled a panel of 30 breast cancer cell lines in the absence of drugs for their mutations, copy number aberrations, mRNA, protein expression and protein phosphorylation, and for response to seven different kinase inhibitors. We then constructed a knowledge-based, Bayesian computational model that integrates these data types and estimates the relative contribution of various drug sensitivity mechanisms. The resulting model of regulatory signaling explained the majority of the variability observed in drug response. The model also identified cell lines with an unexplained response, and for these we searched for novel explanatory factors. Among others, we found that 4E-BP1 protein expression, and not just the extent of phosphorylation, was a determinant of mTOR inhibitor sensitivity. We validated this finding experimentally and found that overexpression of 4E-BP1 in cell lines that normally possess low levels of this protein is sufficient to increase mTOR inhibitor sensitivity. Taken together, our work demonstrates that combining experimental characterization with integrative modeling can be used to systematically test and extend our understanding of the variability in anticancer drug response.Significance: By estimating how different oncogenic mutations and drug resistance mechanisms affect the response of cancer cells to kinase inhibitors, we can better understand and ultimately predict response to these anticancer drugs.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/78/15/4396/F1.large.jpg Cancer Res; 78(15); 4396-410. ©2018 AACR.

Integrative Kinome Profiling Identifies mTORC1/2 Inhibition as Treatment Strategy in Ovarian Clear Cell Carcinoma.

Purpose: Advanced-stage ovarian clear cell carcinoma (OCCC) is unresponsive to conventional platinum-based chemotherapy. Frequent alterations in OCCC include deleterious mutations in the tumor suppressor ARID1A and activating mutations in the PI3K subunit PIK3CA In this study, we aimed to identify currently unknown mutated kinases in patients with OCCC and test druggability of downstream affected pathways in OCCC models.Experimental Design: In a large set of patients with OCCC (n = 124), the human kinome (518 kinases) and additional cancer-related genes were sequenced, and copy-number alterations were determined. Genetically characterized OCCC cell lines (n = 17) and OCCC patient-derived xenografts (n = 3) were used for drug testing of ERBB tyrosine kinase inhibitors erlotinib and lapatinib, the PARP inhibitor olaparib, and the mTORC1/2 inhibitor AZD8055.Results: We identified several putative driver mutations in kinases at low frequency that were not previously annotated in OCCC. Combining mutations and copy-number alterations, 91% of all tumors are affected in the PI3K/AKT/mTOR pathway, the MAPK pathway, or the ERBB family of receptor tyrosine kinases, and 82% in the DNA repair pathway. Strong p-S6 staining in patients with OCCC suggests high mTORC1/2 activity. We consistently found that the majority of OCCC cell lines are especially sensitive to mTORC1/2 inhibition by AZD8055 and not toward drugs targeting ERBB family of receptor tyrosine kinases or DNA repair signaling. We subsequently demonstrated the efficacy of mTORC1/2 inhibition in all our unique OCCC patient-derived xenograft models.Conclusions: These results propose mTORC1/2 inhibition as an effective treatment strategy in OCCC. Clin Cancer Res; 24(16); 3928-40. ©2018 AACR.

Integration of genomic, transcriptomic and proteomic data identifies two biologically distinct subtypes of invasive lobular breast cancer.

Invasive lobular carcinoma (ILC) is the second most frequently occurring histological breast cancer subtype after invasive ductal carcinoma (IDC), accounting for around 10% of all breast cancers. The molecular processes that drive the development of ILC are still largely unknown. We have performed a comprehensive genomic, transcriptomic and proteomic analysis of a large ILC patient cohort and present here an integrated molecular portrait of ILC. Mutations in CDH1 and in the PI3K pathway are the most frequent molecular alterations in ILC. We identified two main subtypes of ILCs: (i) an immune related subtype with mRNA up-regulation of PD-L1, PD-1 and CTLA-4 and greater sensitivity to DNA-damaging agents in representative cell line models; (ii) a hormone related subtype, associated with Epithelial to Mesenchymal Transition (EMT), and gain of chromosomes 1q and 8q and loss of chromosome 11q. Using the somatic mutation rate and eIF4B protein level, we identified three groups with different clinical outcomes, including a group with extremely good prognosis. We provide a comprehensive overview of the molecular alterations driving ILC and have explored links with therapy response. This molecular characterization may help to tailor treatment of ILC through the application of specific targeted, chemo- and/or immune-therapies.