Epistasis and evolutionary dependencies in human cancers.

Cancer evolution is driven by the concerted action of multiple molecular alterations, which emerge and are selected during tumor progression. An alteration is selected when it provides an advantage to the tumor cell. However, the advantage provided by a specific alteration depends on the tumor lineage, cell epigenetic state, and presence of additional alterations. In this case, we say that an evolutionary dependency exists between an alteration and what influences its selection. Epistatic interactions between altered genes lead to evolutionary dependencies (EDs), by favoring or vetoing specific combinations of events. Large-scale cancer genomics studies have discovered examples of such dependencies, and showed that they influence tumor progression, disease phenotypes, and therapeutic response. In the past decade, several algorithmic approaches have been proposed to infer EDs from large-scale genomics datasets. These methods adopt diverse strategies to address common challenges and shed new light on cancer evolutionary trajectories. Here, we review these efforts starting from a simple conceptualization of the problem, presenting the tackled and still unmet needs in the field, and discussing the implications of EDs in cancer biology and precision oncology.

Systematic assessment of gene co-regulation within chromatin domains determines differentially active domains across human cancers.

BACKGROUND: Spatial interactions and insulation of chromatin regions are associated with transcriptional regulation. Domains of frequent chromatin contacts are proposed as functional units, favoring and delimiting gene regulatory interactions. However, contrasting evidence supports the association between chromatin domains and transcription. RESULT: Here, we assess gene co-regulation in chromatin domains across multiple human cancers, which exhibit great transcriptional heterogeneity. Across all datasets, gene co-regulation is observed only within a small yet significant number of chromatin domains. We design an algorithmic approach to identify differentially active domains (DADo) between two conditions and show that these provide complementary information to differentially expressed genes. Domains comprising co-regulated genes are enriched in the less active B sub-compartments and for genes with similar function. Notably, differential activation of chromatin domains is not associated with major changes of domain boundaries, but rather with changes of sub-compartments and intra-domain contacts. CONCLUSION: Overall, gene co-regulation is observed only in a minority of chromatin domains, whose systematic identification will help unravel the relationship between chromatin structure and transcription.

Cell-autonomous inflammation of BRCA1-deficient ovarian cancers drives both tumor-intrinsic immunoreactivity and immune resistance via STING.

In this study, we investigate mechanisms leading to inflammation and immunoreactivity in ovarian tumors with homologous recombination deficiency (HRD). BRCA1 loss is found to lead to transcriptional reprogramming in tumor cells and cell-intrinsic inflammation involving type I interferon (IFN) and stimulator of IFN genes (STING). BRCA1-mutated (BRCA1mut) tumors are thus T cell inflamed at baseline. Genetic deletion or methylation of DNA-sensing/IFN genes or CCL5 chemokine is identified as a potential mechanism to attenuate T cell inflammation. Alternatively, in BRCA1mut cancers retaining inflammation, STING upregulates VEGF-A, mediating immune resistance and tumor progression. Tumor-intrinsic STING elimination reduces neoangiogenesis, increases CD8+ T cell infiltration, and reverts therapeutic resistance to dual immune checkpoint blockade (ICB). VEGF-A blockade phenocopies genetic STING loss and synergizes with ICB and/or poly(ADP-ribose) polymerase (PARP) inhibitors to control the outgrowth of Trp53-/-Brca1-/- but not Brca1+/+ ovarian tumors in vivo, offering rational combinatorial therapies for HRD cancers.

Systematic inference and comparison of multi-scale chromatin sub-compartments connects spatial organization to cell phenotypes.

Chromatin compartmentalization reflects biological activity. However, inference of chromatin sub-compartments and compartment domains from chromosome conformation capture (Hi-C) experiments is limited by data resolution. As a result, these have been characterized only in a few cell types and systematic comparisons across multiple tissues and conditions are missing. Here, we present Calder, an algorithmic approach that enables the identification of multi-scale sub-compartments at variable data resolution. Calder allows to infer and compare chromatin sub-compartments and compartment domains in >100 cell lines. Our results reveal sub-compartments enriched for poised chromatin states and undergoing spatial repositioning during lineage differentiation and oncogenic transformation.

Nongenetic Evolution Drives Lung Adenocarcinoma Spatial Heterogeneity and Progression.

Cancer evolution determines molecular and morphologic intratumor heterogeneity and challenges the design of effective treatments. In lung adenocarcinoma, disease progression and prognosis are associated with the appearance of morphologically diverse tumor regions, termed histologic patterns. However, the link between molecular and histologic features remains elusive. Here, we generated multiomics and spatially resolved molecular profiles of histologic patterns from primary lung adenocarcinoma, which we integrated with molecular data from >2,000 patients. The transition from indolent to aggressive patterns was not driven by genetic alterations but by epigenetic and transcriptional reprogramming reshaping cancer cell identity. A signature quantifying this transition was an independent predictor of patient prognosis in multiple human cohorts. Within individual tumors, highly multiplexed protein spatial profiling revealed coexistence of immune desert, inflamed, and excluded regions, which matched histologic pattern composition. Our results provide a detailed molecular map of lung adenocarcinoma intratumor spatial heterogeneity, tracing nongenetic routes of cancer evolution. SIGNIFICANCE: Lung adenocarcinomas are classified based on histologic pattern prevalence. However, individual tumors exhibit multiple patterns with unknown molecular features. We characterized nongenetic mechanisms underlying intratumor patterns and molecular markers predicting patient prognosis. Intratumor patterns determined diverse immune microenvironments, warranting their study in the context of current immunotherapies.This article is highlighted in the In This Issue feature, p. 1307.

Discovering functional evolutionary dependencies in human cancers.

Cancer cells retain genomic alterations that provide a selective advantage. The prediction and validation of advantageous alterations are major challenges in cancer genomics. Moreover, it is crucial to understand how the coexistence of specific alterations alters response to genetic and therapeutic perturbations. In the present study, we inferred functional alterations and preferentially selected combinations of events in >9,000 human tumors. Using a Bayesian inference framework, we validated computational predictions with high-throughput readouts from genetic and pharmacological screenings on 2,000 cancer cell lines. Mutually exclusive and co-occurring cancer alterations reflected, respectively, functional redundancies able to rescue the phenotype of individual target inhibition, or synergistic interactions, increasing oncogene addiction. Among the top scoring dependencies, co-alteration of the phosphoinositide 3-kinase (PI3K) subunit PIK3CA and the nuclear factor NFE2L2 was a synergistic evolutionary trajectory in squamous cell carcinomas. By integrating computational, experimental and clinical evidence, we provide a framework to study the combinatorial functional effects of cancer genomic alterations.

Cathepsin S Regulates Antigen Processing and T Cell Activity in Non-Hodgkin Lymphoma.

Genomic alterations in cancer cells can influence the immune system to favor tumor growth. In non-Hodgkin lymphoma, physiological interactions between B cells and the germinal center microenvironment are coopted to sustain cancer cell proliferation. We found that follicular lymphoma patients harbor a recurrent hotspot mutation targeting tyrosine 132 (Y132D) in cathepsin S (CTSS) that enhances protein activity. CTSS regulates antigen processing and CD4+ and CD8+ T cell-mediated immune responses. Loss of CTSS activity reduces lymphoma growth by limiting communication with CD4+ T follicular helper cells while inducing antigen diversification and activation of CD8+ T cells. Overall, our results suggest that CTSS inhibition has non-redundant therapeutic potential to enhance anti-tumor immune responses in indolent and aggressive lymphomas.

Dynamic Emergence of Observed and Hidden Intra-tumor Heterogeneity.

Intra-tumor heterogeneity is frequently observed in cancer patients, and it is associated with therapeutic resistance and disease relapse. However, its systematic assessment is still limited and often unfeasible. Here, we use a mathematical model of tumor progression to decipher how multiple clones emerge and organize into complex architectures. We found a trade-off between cancer cell alteration and proliferation rates that defines a transition between low and high heterogeneity, the latter characterized by branching tumor phylogenies. We predict the existence of observed and hidden intra-tumor heterogeneity, which challenges the correct estimation of intrinsic tumor complexity. Although the numbers of observed and hidden clones do not always correlate, we demonstrate that population frequencies of observed clones can be used to estimate the extent of hidden heterogeneity in both simulated and human tumors. The characterization of complex clonal architectures is a critical first step toward understanding their organizing principles and predicting their emergence.

The Oncogenic Action of NRF2 Depends on De-glycation by Fructosamine-3-Kinase.

The NRF2 transcription factor controls a cell stress program that is implicated in cancer and there is great interest in targeting NRF2 for therapy. We show that NRF2 activity depends on Fructosamine-3-kinase (FN3K)-a kinase that triggers protein de-glycation. In its absence, NRF2 is extensively glycated, unstable, and defective at binding to small MAF proteins and transcriptional activation. Moreover, the development of hepatocellular carcinoma triggered by MYC and Keap1 inactivation depends on FN3K in vivo. N-acetyl cysteine treatment partially rescues the effects of FN3K loss on NRF2 driven tumor phenotypes indicating a key role for NRF2-mediated redox balance. Mass spectrometry reveals that other proteins undergo FN3K-sensitive glycation, including translation factors, heat shock proteins, and histones. How glycation affects their functions remains to be defined. In summary, our study reveals a surprising role for the glycation of cellular proteins and implicates FN3K as targetable modulator of NRF2 activity in cancer.

Publisher Correction: Uncoupling protein 2 reprograms the tumor microenvironment to support the anti-tumor immune cycle.

In the version of this article initially published, the bars were not aligned with the data points or horizontal axis labels in Fig. 5d, and the labels along each horizontal axis of Fig. 5j-l indicating the presence (+) or absence (-) of doxycycline (Dox) were incorrectly included with the labels below that axis. Also, the right vertical bar above Fig. 7b linking 'P = 0.0001' to the key was incorrect; the correct comparison is αPD-1 versus Dox + αPD-1. Similarly, the right vertical bar above Fig. 7e linking 'P = 0.0002' to the key was incorrect; the correct comparison is αPD-1 versus Rosig + αPD-1. The errors have been corrected in the HTML and PDF versions of the article.

EZH2 oncogenic mutations drive epigenetic, transcriptional, and structural changes within chromatin domains.

Chromatin is organized into topologically associating domains (TADs) enriched in distinct histone marks. In cancer, gain-of-function mutations in the gene encoding the enhancer of zeste homolog 2 protein (EZH2) lead to a genome-wide increase in histone-3 Lys27 trimethylation (H3K27me3) associated with transcriptional repression. However, the effects of these epigenetic changes on the structure and function of chromatin domains have not been explored. Here, we found a functional interplay between TADs and epigenetic and transcriptional changes mediated by mutated EZH2. Altered EZH2 (p.Tyr646* (EZH2Y646X)) led to silencing of entire domains, synergistically inactivating multiple tumor suppressors. Intra-TAD gene silencing was coupled with changes of interactions between gene promoter regions. Notably, gene expression and chromatin interactions were restored by pharmacological inhibition of EZH2Y646X. Our results indicate that EZH2Y646X alters the topology and function of chromatin domains to promote synergistic oncogenic programs.

Uncoupling protein 2 reprograms the tumor microenvironment to support the anti-tumor immune cycle.

Immune checkpoint blockade therapy has shifted the paradigm for cancer treatment. However, the majority of patients lack effective responses due to insufficient T cell infiltration in tumors. Here we show that expression of mitochondrial uncoupling protein 2 (UCP2) in tumor cells determines the immunostimulatory feature of the tumor microenvironment (TME) and is positively associated with prolonged survival. UCP2 reprograms the immune state of the TME by altering its cytokine milieu in an interferon regulatory factor 5-dependent manner. Consequently, UCP2 boosts the conventional type 1 dendritic cell- and CD8+ T cell-dependent anti-tumor immune cycle and normalizes the tumor vasculature. Finally we show, using either a genetic or pharmacological approach, that induction of UCP2 sensitizes melanomas to programmed cell death protein-1 blockade treatment and elicits effective anti-tumor responses. Together, this study demonstrates that targeting the UCP2 pathway is a potent strategy for alleviating the immunosuppressive TME and overcoming the primary resistance of programmed cell death protein-1 blockade.

Pan-Cancer Landscape of Aberrant DNA Methylation across Human Tumors.

The discovery of cancer-associated alterations has primarily focused on genetic variants. Nonetheless, altered epigenomes contribute to deregulate transcription and promote oncogenic pathways. Here, we designed an algorithmic approach (RESET) to identify aberrant DNA methylation and associated cis-transcriptional changes across >6,000 human tumors. Tumors exhibiting mutations of chromatin remodeling factors and Wnt signaling displayed DNA methylation instability, characterized by numerous hyper- and hypo-methylated loci. Most silenced and enhanced genes coalesced in specific pathways including apoptosis, DNA repair, and cell metabolism. Cancer-germline antigens (CG) were frequently epigenomically enhanced and their expression correlated with response to anti-PD-1, but not anti-CTLA4, in skin melanoma. Finally, we demonstrated the potential of our approach to explore DNA methylation changes in pediatric tumors, which frequently lack genetic drivers and exhibit epigenomic modifications. Our results provide a pan-cancer map of aberrant DNA methylation to inform functional and therapeutic studies.

Pan-cancer inference of intra-tumor heterogeneity reveals associations with different forms of genomic instability.

Genomic instability is a major driver of intra-tumor heterogeneity. However, unstable genomes often exhibit different molecular and clinical phenotypes, which are associated with distinct mutational processes. Here, we algorithmically inferred the clonal phylogenies of ~6,000 human tumors from 32 tumor types to explore how intra-tumor heterogeneity depends on different implementations of genomic instability. We found that extremely unstable tumors associated with DNA repair deficiencies or high chromosomal instability are not the most intrinsically heterogeneous. Conversely, intra-tumor heterogeneity is greatest in tumors exhibiting relatively high numbers of both mutations and copy number alterations, a feature often observed in cancers associated with exogenous mutagens. Independently of the type of instability, tumors with high number of clones invariably evolved through branching phylogenies that could be stratified based on the extent of clonal (early) and subclonal (late) instability. Interestingly, tumors with high number of subclonal mutations frequently exhibited chromosomal instability, TP53 mutations, and APOBEC-related mutational signatures. Vice versa, mutations of chromatin remodeling genes often characterized tumors with few subclonal but multiple clonal mutations. Understanding how intra-tumor heterogeneity depends on genomic instability is critical to identify markers predictive of the tumor complexity and envision therapeutic strategies able to exploit this association.

Oncogenic Signaling Pathways in The Cancer Genome Atlas.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFß signaling, p53 and ß-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

Machine Learning Detects Pan-cancer Ras Pathway Activation in The Cancer Genome Atlas.

Precision oncology uses genomic evidence to match patients with treatment but often fails to identify all patients who may respond. The transcriptome of these "hidden responders" may reveal responsive molecular states. We describe and evaluate a machine-learning approach to classify aberrant pathway activity in tumors, which may aid in hidden responder identification. The algorithm integrates RNA-seq, copy number, and mutations from 33 different cancer types across The Cancer Genome Atlas (TCGA) PanCanAtlas project to predict aberrant molecular states in tumors. Applied to the Ras pathway, the method detects Ras activation across cancer types and identifies phenocopying variants. The model, trained on human tumors, can predict response to MEK inhibitors in wild-type Ras cell lines. We also present data that suggest that multiple hits in the Ras pathway confer increased Ras activity. The transcriptome is underused in precision oncology and, combined with machine learning, can aid in the identification of hidden responders.

Conditional Selection of Genomic Alterations Dictates Cancer Evolution and Oncogenic Dependencies.

Cancer evolves through the emergence and selection of molecular alterations. Cancer genome profiling has revealed that specific events are more or less likely to be co-selected, suggesting that the selection of one event depends on the others. However, the nature of these evolutionary dependencies and their impact remain unclear. Here, we designed SELECT, an algorithmic approach to systematically identify evolutionary dependencies from alteration patterns. By analyzing 6,456 genomes from multiple tumor types, we constructed a map of oncogenic dependencies associated with cellular pathways, transcriptional readouts, and therapeutic response. Finally, modeling of cancer evolution shows that alteration dependencies emerge only under conditional selection. These results provide a framework for the design of strategies to predict cancer progression and therapeutic response.

PD-L1 Is a Therapeutic Target of the Bromodomain Inhibitor JQ1 and, Combined with HLA Class I, a Promising Prognostic Biomarker in Neuroblastoma.

Purpose: This study sought to evaluate the expression of programmed cell death-ligand-1 (PD-L1) and HLA class I on neuroblastoma cells and programmed cell death-1 (PD-1) and lymphocyte activation gene 3 (LAG3) on tumor-infiltrating lymphocytes to better define patient risk stratification and understand whether this tumor may benefit from therapies targeting immune checkpoint molecules.Experimental Design:In situ IHC staining for PD-L1, HLA class I, PD-1, and LAG3 was assessed in 77 neuroblastoma specimens, previously characterized for tumor-infiltrating T-cell density and correlated with clinical outcome. Surface expression of PD-L1 was evaluated by flow cytometry and IHC in neuroblastoma cell lines and tumors genetically and/or pharmacologically inhibited for MYC and MYCN. A dataset of 477 human primary neuroblastomas from GEO and ArrayExpress databases was explored for PD-L1, MYC, and MYCN correlation.Results: Multivariate Cox regression analysis demonstrated that the combination of PD-L1 and HLA class I tumor cell density is a prognostic biomarker for predicting overall survival in neuroblastoma patients (P = 0.0448). MYC and MYCN control the expression of PD-L1 in neuroblastoma cells both in vitro and in vivo Consistently, abundance of PD-L1 transcript correlates with MYC expression in primary neuroblastoma.Conclusions: The combination of PD-L1 and HLA class I represents a novel prognostic biomarker for neuroblastoma. Pharmacologic inhibition of MYCN and MYC may be exploited to target PD-L1 and restore an efficient antitumor immunity in high-risk neuroblastoma. Clin Cancer Res; 23(15); 4462-72. ©2017 AACR.

Tumor-infiltrating T lymphocytes improve clinical outcome of therapy-resistant neuroblastoma.

Neuroblastoma grows within an intricate network of different cell types including epithelial, stromal and immune cells. The presence of tumor-infiltrating T cells is considered an important prognostic indicator in many cancers, but the role of these cells in neuroblastoma remains to be elucidated. Herein, we examined the relationship between the type, density and organization of infiltrating T cells and clinical outcome within a large collection of neuroblastoma samples by quantitative analysis of immunohistochemical staining. We found that infiltrating T cells have a prognostic value greater than, and independent of, the criteria currently used to stage neuroblastoma. A variable in situ structural organization and different concurrent infiltration of T-cell subsets were detected in tumors with various outcomes. Low-risk neuroblastomas were characterized by a higher number of proliferating T cells and a more structured T-cell organization, which was gradually lost in tumors with poor prognosis. We defined an immunoscore based on the presence of CD3+, CD4+ and CD8+ infiltrating T cells that associates with favorable clinical outcome in MYCN-amplified tumors, improving patient survival when combined with the v-myc avian myelocytomatosis viral oncogene neuroblastoma derived homolog (MYCN) status. These findings support the hypothesis that infiltrating T cells influence the behavior of neuroblastoma and might be of clinical importance for the treatment of patients.