Chromosome 7 to the rescue: overcoming chromosome 10 loss in gliomas.

The co-occurrence of chromosome 10 loss and chromosome 7 gain in gliomas is the most frequent loss-gain co-aneuploidy pair in human cancers, a phenomenon that has been investigated without resolution since the late 1980s. Expanding beyond previous gene-centric studies, we investigate the co-occurrence in a genome-wide manner taking an evolutionary perspective. First, by mining large tumor aneuploidy data, we predict that the more likely order is 10 loss followed by 7 gain. Second, by analyzing extensive genomic and transcriptomic data from both patients and cell lines, we find that this co-occurrence can be explained by functional rescue interactions that are highly enriched on 7, which can possibly compensate for any detrimental consequences arising from the loss of 10. Finally, by analyzing transcriptomic data from normal, non-cancerous, human brain tissues, we provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

CRISPR screens decode cancer cell pathways that trigger γδ T cell detection.

γδ T cells are potent anticancer effectors with the potential to target tumours broadly, independent of patient-specific neoantigens or human leukocyte antigen background1-5. γδ T cells can sense conserved cell stress signals prevalent in transformed cells2,3, although the mechanisms behind the targeting of stressed target cells remain poorly characterized. Vγ9Vδ2 T cells-the most abundant subset of human γδ T cells4-recognize a protein complex containing butyrophilin 2A1 (BTN2A1) and BTN3A1 (refs. 6-8), a widely expressed cell surface protein that is activated by phosphoantigens abundantly produced by tumour cells. Here we combined genome-wide CRISPR screens in target cancer cells to identify pathways that regulate γδ T cell killing and BTN3A cell surface expression. The screens showed previously unappreciated multilayered regulation of BTN3A abundance on the cell surface and triggering of γδ T cells through transcription, post-translational modifications and membrane trafficking. In addition, diverse genetic perturbations and inhibitors disrupting metabolic pathways in the cancer cells, particularly ATP-producing processes, were found to alter BTN3A levels. This induction of both BTN3A and BTN2A1 during metabolic crises is dependent on AMP-activated protein kinase (AMPK). Finally, small-molecule activation of AMPK in a cell line model and in patient-derived tumour organoids led to increased expression of the BTN2A1-BTN3A complex and increased Vγ9Vδ2 T cell receptor-mediated killing. This AMPK-dependent mechanism of metabolic stress-induced ligand upregulation deepens our understanding of γδ T cell stress surveillance and suggests new avenues available to enhance γδ T cell anticancer activity.

A landscape of response to drug combinations in non-small cell lung cancer.

Combination of anti-cancer drugs is broadly seen as way to overcome the often-limited efficacy of single agents. The design and testing of combinations are however very challenging. Here we present a uniquely large dataset screening over 5000 targeted agent combinations across 81 non-small cell lung cancer cell lines. Our analysis reveals a profound heterogeneity of response across the tumor models. Notably, combinations very rarely result in a strong gain in efficacy over the range of response observable with single agents. Importantly, gain of activity over single agents is more often seen when co-targeting functionally proximal genes, offering a strategy for designing more efficient combinations. Because combinatorial effect is strongly context specific, tumor specificity should be achievable. The resource provided, together with an additional validation screen sheds light on major challenges and opportunities in building efficacious combinations against cancer and provides an opportunity for training computational models for synergy prediction.

INCISOR: An Algorithm to Identify Synthetic Rescue Mediators of Resistance to Targeted and Immunotherapy.

Despite the success of targeted therapies including immunotherapies in cancer treatments, tumor resistance to targeted therapies remains a fundamental challenge. Tumors can evolve resistance to a therapy that targets one gene by acquiring compensatory alterations in another gene, such compensatory interaction between two genes is referred to as synthetic rescue (SR) interactions. To identify SRs, here we describe an algorithm, INCISOR, that leverages tumor transcriptomics and clinical information from 10,000 patients as well as data from experimental screens. INCISOR can identify SRs that are common across several cancer-types in genome-wide fashion by sifting through half a billion possible gene-gene combinations and provide a framework to design therapies to tackle resistance.

Supporting the Next Generation of Scientists to Lead Cancer Immunology Research.

Recent success in the use of immunotherapy for a broad range of cancers has propelled the field of cancer immunology to the forefront of cancer research. As more and more young investigators join the community of cancer immunologists, the Arthur L. Irving Family Foundation Cancer Immunology Symposium provided a platform to bring this expanding and vibrant community together and support the development of the future leaders in the field. This commentary outlines the lessons that emerged from the inaugural symposium highlighting the areas of scientific and career development that are essential for professional growth in the field of cancer immunology and beyond. Leading scientists and clinicians in the field provided their experience on the topics of scientific trajectory, career trajectory, publishing, fundraising, leadership, mentoring, and collaboration. Herein, we provide a conceptual and practical framework for career development to the broader scientific community.

Therapeutically Increasing MHC-I Expression Potentiates Immune Checkpoint Blockade.

Immune checkpoint blockade (ICB) therapy revolutionized cancer treatment, but many patients with impaired MHC-I expression remain refractory. Here, we combined FACS-based genome-wide CRISPR screens with a data-mining approach to identify drugs that can upregulate MHC-I without inducing PD-L1. CRISPR screening identified TRAF3, a suppressor of the NFκB pathway, as a negative regulator of MHC-I but not PD-L1. The Traf3-knockout gene expression signature is associated with better survival in ICB-naïve patients with cancer and better ICB response. We then screened for drugs with similar transcriptional effects as this signature and identified Second Mitochondria-derived Activator of Caspase (SMAC) mimetics. We experimentally validated that the SMAC mimetic birinapant upregulates MHC-I, sensitizes cancer cells to T cell-dependent killing, and adds to ICB efficacy. Our findings provide preclinical rationale for treating tumors expressing low MHC-I expression with SMAC mimetics to enhance sensitivity to immunotherapy. The approach used in this study can be generalized to identify other drugs that enhance immunotherapy efficacy. SIGNIFICANCE: MHC-I loss or downregulation in cancer cells is a major mechanism of resistance to T cell-based immunotherapies. Our study reveals that birinapant may be used for patients with low baseline MHC-I to enhance ICB response. This represents promising immunotherapy opportunities given the biosafety profile of birinapant from multiple clinical trials.This article is highlighted in the In This Issue feature, p. 1307.

Clonal tracing reveals diverse patterns of response to immune checkpoint blockade.

BACKGROUND: Immune checkpoint blockade (ICB) therapy has improved patient survival in a variety of cancers, but only a minority of cancer patients respond. Multiple studies have sought to identify general biomarkers of ICB response, but elucidating the molecular and cellular drivers of resistance for individual tumors remains challenging. We sought to determine whether a tumor with defined genetic background exhibits a stereotypic or heterogeneous response to ICB treatment. RESULTS: We establish a unique mouse system that utilizes clonal tracing and mathematical modeling to monitor the growth of each cancer clone, as well as the bulk tumor, in response to ICB. We find that tumors derived from the same clonal populations showed heterogeneous ICB response and diverse response patterns. Primary response is associated with higher immune infiltration and leads to enrichment of pre-existing ICB-resistant cancer clones. We further identify several cancer cell-intrinsic gene expression signatures associated with ICB resistance, including increased interferon response genes and glucocorticoid response genes. These findings are supported by clinical data from ICB treatment cohorts. CONCLUSIONS: Our study demonstrates diverse response patterns from the same ancestor cancer cells in response to ICB. This suggests the value of monitoring clonal constitution and tumor microenvironment over time to optimize ICB response and to design new combination therapies. Furthermore, as ICB response may enrich for cancer cell-intrinsic resistance signatures, this can affect interpretations of tumor RNA-seq data for response-signature association studies.

Immune receptor repertoires in pediatric and adult acute myeloid leukemia.

BACKGROUND: Acute myeloid leukemia (AML), caused by the abnormal proliferation of immature myeloid cells in the blood or bone marrow, is one of the most common hematologic malignancies. Currently, the interactions between malignant myeloid cells and the immune microenvironment, especially T cells and B cells, remain poorly characterized. METHODS: In this study, we systematically analyzed the T cell receptor and B cell receptor (TCR and BCR) repertoires from the RNA-seq data of 145 pediatric and 151 adult AML samples as well as 73 non-tumor peripheral blood samples. RESULTS: We inferred over 225,000 complementarity-determining region 3 (CDR3) sequences in TCR α, ß, γ, and δ chains and 1,210,000 CDR3 sequences in B cell immunoglobulin (Ig) heavy and light chains. We found higher clonal expansion of both T cells and B cells in the AML microenvironment and observed many differences between pediatric and adult AML. Most notably, adult AML samples have significantly higher level of B cell activation and more secondary Ig class switch events than pediatric AML or non-tumor samples. Furthermore, adult AML with highly expanded IgA2 B cells, which might represent an immunosuppressive microenvironment, are associated with regulatory T cells and worse overall survival. CONCLUSIONS: Our comprehensive characterization of the AML immune receptor repertoires improved our understanding of T cell and B cell immunity in AML, which may provide insights into immunotherapies in hematological malignancies.

Allele-specific enhancers mediate associations between LCAT and ABCA1 polymorphisms and HDL metabolism.

For most complex traits, the majority of SNPs identified through genome-wide association studies (GWAS) reside within noncoding regions that have no known function. However, these regions are enriched for the regulatory enhancers specific to the cells relevant to the specific trait. Indeed, many of the GWAS loci that have been functionally characterized lie within enhancers that regulate expression levels of key genes. In order to identify polymorphisms with potential allele-specific regulatory effects, we developed a bioinformatics pipeline that harnesses epigenetic signatures as well as transcription factor (TF) binding motifs to identify putative enhancers containing a SNP with potential allele-specific TF binding in linkage disequilibrium (LD) with a GWAS-identified SNP. We applied the approach to GWAS findings for blood lipids, revealing 7 putative enhancers harboring associated SNPs, 3 of which lie within the introns of LCAT and ABCA1, genes that play crucial roles in cholesterol biogenesis and lipoprotein metabolism. All 3 enhancers demonstrated allele-specific in vitro regulatory activity in liver-derived cell lines. We demonstrated that these putative enhancers are in close physical proximity to the promoters of their respective genes, in situ, likely through chromatin looping. In addition, the associated alleles altered the likelihood of transcription activator STAT3 binding. Our results demonstrate that through our approach, the LD blocks that contain GWAS signals, often hundreds of kilobases in size with multiple SNPs serving as statistical proxies to the true functional site, can provide an experimentally testable hypothesis for the underlying regulatory mechanism linking genetic variants to complex traits.

Genome-wide prediction of synthetic rescue mediators of resistance to targeted and immunotherapy.

Most patients with advanced cancer eventually acquire resistance to targeted therapies, spurring extensive efforts to identify molecular events mediating therapy resistance. Many of these events involve synthetic rescue (SR) interactions, where the reduction in cancer cell viability caused by targeted gene inactivation is rescued by an adaptive alteration of another gene (the rescuer). Here, we perform a genome-wide in silico prediction of SR rescuer genes by analyzing tumor transcriptomics and survival data of 10,000 TCGA cancer patients. Predicted SR interactions are validated in new experimental screens. We show that SR interactions can successfully predict cancer patients' response and emerging resistance. Inhibiting predicted rescuer genes sensitizes resistant cancer cells to therapies synergistically, providing initial leads for developing combinatorial approaches to overcome resistance proactively. Finally, we show that the SR analysis of melanoma patients successfully identifies known mediators of resistance to immunotherapy and predicts novel rescuers.

Study of cell differentiation by phylogenetic analysis using histone modification data.

BACKGROUND: In cell differentiation, a cell of a less specialized type becomes one of a more specialized type, even though all cells have the same genome. Transcription factors and epigenetic marks like histone modifications can play a significant role in the differentiation process. RESULTS: In this paper, we present a simple analysis of cell types and differentiation paths using phylogenetic inference based on ChIP-Seq histone modification data. We precisely defined the notion of cell-type trees and provided a procedure of building such trees. We propose new data representation techniques and distance measures for ChIP-Seq data and use these together with standard phylogenetic inference methods to build biologically meaningful cell-type trees that indicate how diverse types of cells are related. We demonstrate our approach on various kinds of histone modifications for various cell types, also using the datasets to explore various issues surrounding replicate data, variability between cells of the same type, and robustness. We use the results to get some interesting biological findings like important patterns of histone modification changes during cell differentiation process. CONCLUSIONS: We introduced and studied the novel problem of inferring cell type trees from histone modification data. The promising results we obtain point the way to a new approach to the study of cell differentiation. We also discuss how cell-type trees can be used to study the evolution of cell types.

Epigenomic model of cardiac enhancers with application to genome wide association studies.

Mammalian gene regulation is often mediated by distal enhancer elements, in particular, for tissue specific and developmental genes. Computational identification of enhancers is difficult because they do not exhibit clear location preference relative to their target gene and also because they lack clearly distinguishing genomic features. This represents a major challenge in deciphering transcriptional regulation. Recent ChIP-seq based genome-wide investigation of epigenomic modifications have revealed that enhancers are often enriched for certain epigenomic marks. Here we utilize the epigenomic data in human heart tissue along with validated human heart enhancers to develop a Support Vector Machine (SVM) model of cardiac enhancers. Cross-validation classification accuracy of our model was 84% and 92% on positive and negative sets respectively with ROC AUC = 0.92. More importantly, while P300 binding has been used as gold standard for enhancers, our model can distinguish P300-bound validated enhancers from other P300-bound regions that failed to exhibit enhancer activity in transgenic mouse. While GWAS studies reveal polymorphic regions associated with certain phenotypes, they do not immediately provide causality. Next, we hypothesized that genomic regions containing a GWAS SNP associated with a cardiac phenotype might contain another SNP in a cardiac enhancer, which presumably mediates the phenotype. Starting with a comprehensive set of SNPs associated with cardiac phenotypes in GWAS studies, we scored other SNPs in LD with the GWAS SNP according to its probability of being an enhancer and choose one with best score in the LD as enhancer. We found that our predicted enhancers are enriched for known cardiac transcriptional regulator motifs and are likely to regulate the nearby gene. Importantly, these tendencies are more favorable for the predicted enhancers compared with an approach that uses P300 binding as a marker of enhancer activity.

ChIPnorm: a statistical method for normalizing and identifying differential regions in histone modification ChIP-seq libraries.

The advent of high-throughput technologies such as ChIP-seq has made possible the study of histone modifications. A problem of particular interest is the identification of regions of the genome where different cell types from the same organism exhibit different patterns of histone enrichment. This problem turns out to be surprisingly difficult, even in simple pairwise comparisons, because of the significant level of noise in ChIP-seq data. In this paper we propose a two-stage statistical method, called ChIPnorm, to normalize ChIP-seq data, and to find differential regions in the genome, given two libraries of histone modifications of different cell types. We show that the ChIPnorm method removes most of the noise and bias in the data and outperforms other normalization methods. We correlate the histone marks with gene expression data and confirm that histone modifications H3K27me3 and H3K4me3 act as respectively a repressor and an activator of genes. Compared to what was previously reported in the literature, we find that a substantially higher fraction of bivalent marks in ES cells for H3K27me3 and H3K4me3 move into a K27-only state. We find that most of the promoter regions in protein-coding genes have differential histone-modification sites. The software for this work can be downloaded from http://lcbb.epfl.ch/software.html.