The HVEM (TNFRSF14) receptor gene is among the most frequently mutated genes in germinal center lymphomas. We report that loss of HVEM leads to cell-autonomous activation of B cell proliferation and drives the development of GC lymphomas in vivo. HVEM-deficient lymphoma B cells also induce a tumor-supportive microenvironment marked by exacerbated lymphoid stroma activation and increased recruitment of T follicular helper (TFH) cells. These changes result from the disruption of inhibitory cell-cell interactions between the HVEM and BTLA (B and T lymphocyte attenuator) receptors. Accordingly, administration of the HVEM ectodomain protein (solHVEM(P37-V202)) binds BTLA and restores tumor suppression. To deliver solHVEM to lymphomas in vivo, we engineered CD19-targeted chimeric antigen receptor (CAR) T cells that produce solHVEM locally and continuously. These modified CAR-T cells show enhanced therapeutic activity against xenografted lymphomas. Hence, the HVEM-BTLA axis opposes lymphoma development, and our study illustrates the use of CAR-T cells as "micro-pharmacies" able to deliver an anti-cancer protein.

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.

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.

Chromatin folding gives rise to structural elements among which are clusters of densely interacting DNA regions termed topologically associating domains (TADs). TADs have been characterized across multiple species, tissue types, and differentiation stages, sometimes in association with regulation of biological functions. The reliability and reproducibility of these findings are intrinsically related with the correct identification of these domains from high-throughput chromatin conformation capture (Hi-C) experiments.

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.

NOTCH1 is a well-established lineage specifier for T cells and among the most frequently mutated genes throughout all subclasses of T cell acute lymphoblastic leukemia (T-ALL). How oncogenic NOTCH1 signaling launches a leukemia-prone chromatin landscape during T-ALL initiation is unknown. Here we demonstrate an essential role for the high-mobility-group transcription factor Tcf1 in orchestrating chromatin accessibility and topology, allowing aberrant Notch1 signaling to convey its oncogenic function. Although essential, Tcf1 is not sufficient to initiate leukemia. The formation of a leukemia-prone epigenetic landscape at the distal Notch1-regulated Myc enhancer, which is fundamental to this disease, is Tcf1-dependent and occurs within the earliest progenitor stage even before cells adopt a T lymphocyte or leukemic fate. Moreover, we discovered a unique evolutionarily conserved Tcf1-regulated enhancer element in the distal Myc-enhancer, which is important for the transition of preleukemic cells to full-blown disease.

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.

Uveal melanoma (UM) is a rare intraocular tumor that, similar to cutaneous melanoma, originates from melanocytes. To gain insights into its genetics, we performed whole-genome sequencing at very deep coverage of tumor-control pairs in 33 samples (24 primary and 9 metastases). Genome-wide, the number of coding mutations was rather low (only 17 variants per tumor on average; range 7-28), thus radically different from cutaneous melanoma, where hundreds of exonic DNA insults are usually detected. Furthermore, no UV light-induced mutational signature was identified. Recurrent coding mutations were found in the known UM drivers GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. Other genes, i.e., TP53BP1, CSMD1, TTC28, DLK2, and KTN1, were also found to harbor somatic mutations in more than one individual, possibly indicating a previously undescribed association with UM pathogenesis. De novo assembly of unmatched reads from non-coding DNA revealed peculiar copy-number variations defining specific UM subtypes, which in turn could be associated with metastatic transformation. Mutational-driven comparison with other tumor types showed that UM is very similar to pediatric tumors, characterized by very few somatic insults and, possibly, important epigenetic changes. Through the analysis of whole-genome sequencing data, our findings shed new light on the molecular genetics of uveal melanoma, delineating it as an atypical tumor of the adult for which somatic events other than mutations in exonic DNA shape its genetic landscape and define its metastatic potential.

Melanoma susceptibility differs significantly in male versus female populations. Low levels of androgen receptor (AR) in melanocytes of the two sexes are accompanied by heterogeneous expression at various stages of the disease. Irrespective of expression levels, genetic and pharmacological suppression of AR activity in melanoma cells blunts proliferation and induces senescence, while increased AR expression or activation exert opposite effects. AR down-modulation elicits a shared gene expression signature associated with better patient survival, related to interferon and cytokine signaling and DNA damage/repair. AR loss leads to dsDNA breakage, cytoplasmic leakage, and STING activation, with AR anchoring the DNA repair proteins Ku70/Ku80 to RNA Pol II and preventing RNA Pol II-associated DNA damage. AR down-modulation or pharmacological inhibition suppresses melanomagenesis, with increased intratumoral infiltration of macrophages and, in an immune-competent mouse model, cytotoxic T cells. AR provides an attractive target for improved management of melanoma independent of patient sex.

Whole-genome doubling (WGD) is a recurrent event in human cancers and it promotes chromosomal instability and acquisition of aneuploidies1-8. However, the three-dimensional organization of chromatin in WGD cells and its contribution to oncogenic phenotypes are currently unknown. Here we show that in p53-deficient cells, WGD induces loss of chromatin segregation (LCS). This event is characterized by reduced segregation between short and long chromosomes, A and B subcompartments and adjacent chromatin domains. LCS is driven by the downregulation of CTCF and H3K9me3 in cells that bypassed activation of the tetraploid checkpoint. Longitudinal analyses revealed that LCS primes genomic regions for subcompartment repositioning in WGD cells. This results in chromatin and epigenetic changes associated with oncogene activation in tumours ensuing from WGD cells. Notably, subcompartment repositioning events were largely independent of chromosomal alterations, which indicates that these were complementary mechanisms contributing to tumour development and progression. Overall, LCS initiates chromatin conformation changes that ultimately result in oncogenic epigenetic and transcriptional modifications, which suggests that chromatin evolution is a hallmark of WGD-driven cancer.

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.

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.

The eIF4E translation initiation factor has oncogenic properties and concordantly, the inhibitory eIF4E-binding protein (4EBP1) is considered a tumor suppressor. The exact molecular effects of 4EBP1 activation in cancer are still unknown. Surprisingly, 4EBP1 is a target of genomic copy number gains (Chr. 8p11) in breast and lung cancer. We noticed that 4EBP1 gains are genetically linked to gains in neighboring genes, including WHSC1L1 and FGFR1. Our results show that FGFR1 gains act to attenuate the function of 4EBP1 via PI3K-mediated phosphorylation at Thr37/46, Ser65, and Thr70 sites. This implies that not 4EBP1 but instead FGFR1 is the genetic target of Chr. 8p11 gains in breast and lung cancer. Accordingly, these tumors show increased sensitivity to FGFR1 and PI3K inhibition, and this is a therapeutic vulnerability through restoring the tumor-suppressive function of 4EBP1. Ribosome profiling reveals genes involved in insulin signaling, glucose metabolism, and the inositol pathway to be the relevant translational targets of 4EBP1. These mRNAs are among the top 200 translation targets and are highly enriched for structure and sequence motifs in their 5'UTR, which depends on the 4EBP1-EIF4E activity. In summary, we identified the translational targets of 4EBP1-EIF4E that facilitate the tumor suppressor function of 4EBP1 in cancer.

Uveal melanoma (UM) is a rare cancer resulting from the transformation of melanocytes in the uveal tract. Integrative analysis has identified four molecular and clinical subsets of UM. To improve our molecular understanding of UM, we performed extensive multi-omics characterization comparing two aggressive UM patient-derived xenograft models with normal choroidal melanocytes, including DNA optical mapping, specific histone modifications, and DNA topology analysis using Hi-C. Our gene expression and cytogenetic analyses suggest that genomic instability is a hallmark of UM. We also identified a recurrent deletion in the BAP1 promoter resulting in loss of expression and associated with high risk of metastases in UM patients. Hi-C revealed chromatin topology changes associated with the upregulation of PRAME, an independent prognostic biomarker in UM, and a potential therapeutic target. Our findings illustrate how multi-omics approaches can improve our understanding of tumorigenesis and reveal two distinct mechanisms of gene expression dysregulation in UM.

Local network alignment is an important component of the analysis of protein-protein interaction networks that may lead to the identification of evolutionary related complexes. We present AlignNemo, a new algorithm that, given the networks of two organisms, uncovers subnetworks of proteins that relate in biological function and topology of interactions. The discovered conserved subnetworks have a general topology and need not to correspond to specific interaction patterns, so that they more closely fit the models of functional complexes proposed in the literature. The algorithm is able to handle sparse interaction data with an expansion process that at each step explores the local topology of the networks beyond the proteins directly interacting with the current solution. To assess the performance of AlignNemo, we ran a series of benchmarks using statistical measures as well as biological knowledge. Based on reference datasets of protein complexes, AlignNemo shows better performance than other methods in terms of both precision and recall. We show our solutions to be biologically sound using the concept of semantic similarity applied to Gene Ontology vocabularies. The binaries of AlignNemo and supplementary details about the algorithms and the experiments are available at: sourceforge.net/p/alignnemo.

In cancer genomics, recurrence of mutations in independent tumor samples is a strong indicator of functional impact. However, rare functional mutations can escape detection by recurrence analysis owing to lack of statistical power. We enhance statistical power by extending the notion of recurrence of mutations from single genes to gene families that share homologous protein domains. Domain mutation analysis also sharpens the functional interpretation of the impact of mutations, as domains more succinctly embody function than entire genes. By mapping mutations in 22 different tumor types to equivalent positions in multiple sequence alignments of domains, we confirm well-known functional mutation hotspots, identify uncharacterized rare variants in one gene that are equivalent to well-characterized mutations in another gene, detect previously unknown mutation hotspots, and provide hypotheses about molecular mechanisms and downstream effects of domain mutations. With the rapid expansion of cancer genomics projects, protein domain hotspot analysis will likely provide many more leads linking mutations in proteins to the cancer phenotype.

Pancreatic neuroendocrine tumors (PanNET) comprise two molecular subtypes, relatively benign islet tumors (IT) and invasive, metastasis-like primary (MLP) tumors. Until now, the origin of aggressive MLP tumors has been obscure. Herein, using multi-omics approaches, we revealed that MLP tumors arise from IT via dedifferentiation following a reverse trajectory along the developmental pathway of islet β cells, which results in the acquisition of a progenitor-like molecular phenotype. Functionally, the miR-181cd cluster induces the IT-to-MLP transition by suppressing expression of the Meis2 transcription factor, leading to upregulation of a developmental transcription factor, Hmgb3. Notably, the IT-to-MLP transition constitutes a distinct step of tumorigenesis and is separable from the classic proliferation-associated hallmark, temporally preceding accelerated proliferation of cancer cells. Furthermore, patients with PanNET with elevated HMGB3 expression and an MLP transcriptional signature are associated with higher-grade tumors and worse survival. Overall, our results unveil a new mechanism that modulates cancer cell plasticity to enable malignant progression. SIGNIFICANCE: Dedifferentiation has long been observed as a histopathologic characteristic of many cancers, albeit inseparable from concurrent increases in cell proliferation. Herein, we demonstrate that dedifferentiation is a mechanistically and temporally separable step in the multistage tumorigenesis of pancreatic islet cells, retracing the developmental lineage of islet β cells.This article is highlighted in the In This Issue feature, p. 2355.

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.

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.

Cancer therapy is challenged by the diversity of molecular implementations of oncogenic processes and by the resulting variation in therapeutic responses. Projects such as The Cancer Genome Atlas (TCGA) provide molecular tumor maps in unprecedented detail. The interpretation of these maps remains a major challenge. Here we distilled thousands of genetic and epigenetic features altered in cancers to ∼500 selected functional events (SFEs). Using this simplified description, we derived a hierarchical classification of 3,299 TCGA tumors from 12 cancer types. The top classes are dominated by either mutations (M class) or copy number changes (C class). This distinction is clearest at the extremes of genomic instability, indicating the presence of different oncogenic processes. The full hierarchy shows functional event patterns characteristic of multiple cross-tissue groups of tumors, termed oncogenic signature classes. Targetable functional events in a tumor class are suggestive of class-specific combination therapy. These results may assist in the definition of clinical trials to match actionable oncogenic signatures with personalized therapies.