Integrated epigenetic and transcriptional single-cell analysis of t(11;14) multiple myeloma and its BCL2 dependency.

ABSTRACT: The translocation t(11;14) occurs in 20% of patients with multiple myeloma (MM) and results in the upregulation of CCND1. Nearly two-thirds of t(11;14) MM cells are BCL2 primed and highly responsive to the oral BCL2 inhibitor venetoclax. Although it is evident that this unique sensitivity to venetoclax depends on the Bcl-2 homology domain 3- proapoptotic protein priming of BCL2, the biology underlying t(11;14) MM dependency on BCL2 is poorly defined. Importantly, the epigenetic regulation of t(11;14) transcriptomes and its impact on gene regulation and clinical response to venetoclax remain elusive. In this study, by integrating assay for transposase-accessible chromatin by sequencing (ATAC-seq) and RNA-seq at the single-cell level in primary MM samples, we have defined the epigenetic regulome and transcriptome associated with t(11;14) MM. A B-cell-like epigenetic signature was enriched in t(11;14) MM, confirming its phylogeny link to B-cell rather than plasma cell biology. Of note, a loss of a B-cell-like epigenetic signature with a gain of canonical plasma cell transcription factors was observed at the time of resistance to venetoclax. In addition, MCL1 and BCL2L1 copy number gains and structural rearrangements were linked to venetoclax resistance in patients with t(11;14) MM. To date, this is the first study in which both single-cell (sc) ATAC-seq and scRNA-seq analysis are integrated into primary MM cells to obtain a deeper resolution of the epigenetic regulome and transcriptome associated with t(11;14) MM biology and venetoclax resistance.

ETV4-Dependent Transcriptional Plasticity Maintains MYC Expression and Results in IMiD Resistance in Multiple Myeloma.

Immunomodulatory drugs (IMiD) are a backbone therapy for multiple myeloma (MM). Despite their efficacy, most patients develop resistance, and the mechanisms are not fully defined. Here, we show that IMiD responses are directed by IMiD-dependent degradation of IKZF1 and IKZF3 that bind to enhancers necessary to sustain the expression of MYC and other myeloma oncogenes. IMiD treatment universally depleted chromatin-bound IKZF1, but eviction of P300 and BRD4 coactivators only occurred in IMiD-sensitive cells. IKZF1-bound enhancers overlapped other transcription factor binding motifs, including ETV4. Chromatin immunoprecipitation sequencing showed that ETV4 bound to the same enhancers as IKZF1, and ETV4 CRISPR/Cas9-mediated ablation resulted in sensitization of IMiD-resistant MM. ETV4 expression is associated with IMiD resistance in cell lines, poor prognosis in patients, and is upregulated at relapse. These data indicate that ETV4 alleviates IKZF1 and IKZF3 dependency in MM by maintaining oncogenic enhancer activity and identify transcriptional plasticity as a previously unrecognized mechanism of IMiD resistance. SIGNIFICANCE: We show that IKZF1-bound enhancers are critical for IMiD efficacy and that the factor ETV4 can bind the same enhancers and substitute for IKZF1 and mediate IMiD resistance by maintaining MYC and other oncogenes. These data implicate transcription factor redundancy as a previously unrecognized mode of IMiD resistance in MM. See related article by Welsh, Barwick, et al., p. 34. See related commentary by Yun and Cleveland, p. 5. This article is featured in Selected Articles from This Issue, p. 4.

Mechanisms of antigen escape from BCMA- or GPRC5D-targeted immunotherapies in multiple myeloma.

B cell maturation antigen (BCMA) target loss is considered to be a rare event that mediates multiple myeloma (MM) resistance to anti-BCMA chimeric antigen receptor T cell (CAR T) or bispecific T cell engager (TCE) therapies. Emerging data report that downregulation of G-protein-coupled receptor family C group 5 member D (GPRC5D) protein often occurs at relapse after anti-GPRC5D CAR T therapy. To examine the tumor-intrinsic factors that promote MM antigen escape, we performed combined bulk and single-cell whole-genome sequencing and copy number variation analysis of 30 patients treated with anti-BCMA and/or anti-GPRC5D CAR T/TCE therapy. In two cases, MM relapse post-TCE/CAR T therapy was driven by BCMA-negative clones harboring focal biallelic deletions at the TNFRSF17 locus at relapse or by selective expansion of pre-existing subclones with biallelic TNFRSF17 loss. In another five cases of relapse, newly detected, nontruncating, missense mutations or in-frame deletions in the extracellular domain of BCMA negated the efficacies of anti-BCMA TCE therapies, despite detectable surface BCMA protein expression. In the present study, we also report four cases of MM relapse with biallelic mutations of GPRC5D after anti-GPRC5D TCE therapy, including two cases with convergent evolution where multiple subclones lost GPRC5D through somatic events. Immunoselection of BCMA- or GPRC5D-negative or mutant clones is an important tumor-intrinsic driver of relapse post-targeted therapies. Mutational events on BCMA confer distinct sensitivities toward different anti-BCMA therapies, underscoring the importance of considering the tumor antigen landscape for optimal design and selection of targeted immunotherapies in MM.

The pre-existing T cell landscape determines the response to bispecific T cell engagers in multiple myeloma patients.

Bispecific T cell engagers (TCEs) have shown promise in the treatment of various cancers, but the immunological mechanism and molecular determinants of primary and acquired resistance to TCEs remain poorly understood. Here, we identify conserved behaviors of bone marrow-residing T cells in multiple myeloma patients undergoing BCMAxCD3 TCE therapy. We show that the immune repertoire reacts to TCE therapy with cell state-dependent clonal expansion and find evidence supporting the coupling of tumor recognition via major histocompatibility complex class I (MHC class I), exhaustion, and clinical response. We find the abundance of exhausted-like CD8+ T cell clones to be associated with clinical response failure, and we describe loss of target epitope and MHC class I as tumor-intrinsic adaptations to TCEs. These findings advance our understanding of the in vivo mechanism of TCE treatment in humans and provide the rationale for predictive immune-monitoring and conditioning of the immune repertoire to guide future immunotherapy in hematological malignancies.

Needlestack: an ultra-sensitive variant caller for multi-sample next generation sequencing data.

The emergence of next-generation sequencing (NGS) has revolutionized the way of reaching a genome sequence, with the promise of potentially providing a comprehensive characterization of DNA variations. Nevertheless, detecting somatic mutations is still a difficult problem, in particular when trying to identify low abundance mutations, such as subclonal mutations, tumour-derived alterations in body fluids or somatic mutations from histological normal tissue. The main challenge is to precisely distinguish between sequencing artefacts and true mutations, particularly when the latter are so rare they reach similar abundance levels as artefacts. Here, we present needlestack, a highly sensitive variant caller, which directly learns from the data the level of systematic sequencing errors to accurately call mutations. Needlestack is based on the idea that the sequencing error rate can be dynamically estimated from analysing multiple samples together. We show that the sequencing error rate varies across alterations, illustrating the need to precisely estimate it. We evaluate the performance of needlestack for various types of variations, and we show that needlestack is robust among positions and outperforms existing state-of-the-art method for low abundance mutations. Needlestack, along with its source code is freely available on the GitHub platform: https://github.com/IARCbioinfo/needlestack.

Deregulation of Adaptive T Cell Immunity in Multiple Myeloma: Insights Into Mechanisms and Therapeutic Opportunities.

Immunotherapy has recently emerged as a promising treatment option for multiple myeloma (MM) patients. Profound immune dysfunction and evasion of immune surveillance are known to characterize MM evolution and disease progression. Along with genomic changes observed in malignant plasma cells, the bone marrow (BM) milieu creates a protective environment sustained by the complex interaction of BM stromal cells (BMSCs) and malignant cells that using bidirectional connections and cytokines released stimulate disease progression, drug resistance and enable immune escape. Local immune suppression and T-cell exhaustion are important mediating factors of clinical outcomes and responses to immune-based approaches. Thus, further characterization of the defects present in the immune system of MM patients is essential to develop novel therapies and to repurpose the existing ones. This review seeks to provide insights into the mechanisms that promote tumor escape, cause inadequate T-cell stimulation and impaired cytotoxicity in MM. Furthermore, it highlights current immunotherapies being used to restore adaptive T-cell immune responses in MM and describes strategies created to escape these multiple immune evasion mechanisms.

Integrative genomic profiling of large-cell neuroendocrine carcinomas reveals distinct subtypes of high-grade neuroendocrine lung tumors.

Pulmonary large-cell neuroendocrine carcinomas (LCNECs) have similarities with other lung cancers, but their precise relationship has remained unclear. Here we perform a comprehensive genomic (n = 60) and transcriptomic (n = 69) analysis of 75 LCNECs and identify two molecular subgroups: "type I LCNECs" with bi-allelic TP53 and STK11/KEAP1 alterations (37%), and "type II LCNECs" enriched for bi-allelic inactivation of TP53 and RB1 (42%). Despite sharing genomic alterations with adenocarcinomas and squamous cell carcinomas, no transcriptional relationship was found; instead LCNECs form distinct transcriptional subgroups with closest similarity to SCLC. While type I LCNECs and SCLCs exhibit a neuroendocrine profile with ASCL1high/DLL3high/NOTCHlow, type II LCNECs bear TP53 and RB1 alterations and differ from most SCLC tumors with reduced neuroendocrine markers, a pattern of ASCL1low/DLL3low/NOTCHhigh, and an upregulation of immune-related pathways. In conclusion, LCNECs comprise two molecularly defined subgroups, and distinguishing them from SCLC may allow stratified targeted treatment of high-grade neuroendocrine lung tumors.

New Insights into the Molecular Characteristics of Pulmonary Carcinoids and Large Cell Neuroendocrine Carcinomas, and the Impact on Their Clinical Management.

Carcinoids and large cell neuroendocrine carcinomas (LCNECs) are rare neuroendocrine lung tumors. Here we provide an overview of the most updated data on the molecular characteristics of these diseases. Recent genomic studies showed that carcinoids generally contain a low mutational burden and few recurrently mutated genes. Most of the reported mutations occur in chromatin-remodeling genes (e.g., menin 1 gene [MEN1]), and few affect genes of the phosphoinositide 3-kinase (PI3K)-AKT-mechanistic target of rapamycin gene pathway. Aggressive disease has been related to chromothripsis, DNA-repair gene mutations, loss of orthopedia homeobox/CD44, and upregulation of ret proto-oncogene gene (RET) gene expression. In the case of LCNECs, which present with a high mutation burden, two major molecular subtypes have been identified: one with biallelic inactivation of tumor protein p53 gene (TP53) and retinoblastoma gene (RB1), a hallmark of SCLC; and the other one with biallelic inactivation of TP53 and serine/threonine kinase 11 gene (STK11)/kelch like ECH associated protein 1 gene (KEAP1), genes that are frequently mutated in NSCLC. These data, together with the identification of common mutations in the different components of combined LCNEC tumors, provide further evidence of the close molecular relation of LCNEC with other lung tumor types. In terms of therapeutic options, future studies should explore the association between mechanistic target of rapamycin pathway mutations and response to mechanistic target of rapamycin inhibitors in carcinoids. For LCNEC, preliminary data suggest that the two molecular subtypes might have a predictive value for chemotherapy response, but this observation needs to be validated in randomized prospective clinical trials. Finally, delta like Notch canonical ligand 3 inhibitors and immunotherapy may provide alternative options for patient-tailored therapy in LCNEC.

Molecular Subtypes of Pulmonary Large-cell Neuroendocrine Carcinoma Predict Chemotherapy Treatment Outcome.

Purpose: Previous genomic studies have identified two mutually exclusive molecular subtypes of large-cell neuroendocrine carcinoma (LCNEC): the RB1 mutated (mostly comutated with TP53) and the RB1 wild-type groups. We assessed whether these subtypes have a predictive value on chemotherapy outcome.Experimental Design: Clinical data and tumor specimens were retrospectively obtained from the Netherlands Cancer Registry and Pathology Registry. Panel-consensus pathology revision confirmed the diagnosis of LCNEC in 148 of 232 cases. Next-generation sequencing (NGS) for TP53, RB1, STK11, and KEAP1 genes, as well as IHC for RB1 and P16 was performed on 79 and 109 cases, respectively, and correlated with overall survival (OS) and progression-free survival (PFS), stratifying for non-small cell lung cancer type chemotherapy including platinum + gemcitabine or taxanes (NSCLC-GEM/TAX) and platinum-etoposide (SCLC-PE).Results:RB1 mutation and protein loss were detected in 47% (n = 37) and 72% (n = 78) of the cases, respectively. Patients with RB1 wild-type LCNEC treated with NSCLC-GEM/TAX had a significantly longer OS [9.6; 95% confidence interval (CI), 7.7-11.6 months] than those treated with SCLC-PE [5.8 (5.5-6.1); P = 0.026]. Similar results were obtained for patients expressing RB1 in their tumors (P = 0.001). RB1 staining or P16 loss showed similar results. The same outcome for chemotherapy treatment was observed in LCNEC tumors harboring an RB1 mutation or lost RB1 protein.Conclusions: Patients with LCNEC tumors that carry a wild-type RB1 gene or express the RB1 protein do better with NSCLC-GEM/TAX treatment than with SCLC-PE chemotherapy. However, no difference was observed for RB1 mutated or with lost protein expression. Clin Cancer Res; 24(1); 33-42. ©2017 AACR.

BAP1 Is Altered by Copy Number Loss, Mutation, and/or Loss of Protein Expression in More Than 70% of Malignant Peritoneal Mesotheliomas.

INTRODUCTION: Malignant mesothelioma is a deadly disease that is strongly associated with asbestos exposure. Peritoneal mesotheliomas account for 10% of all the cases. BRCA1 associated protein 1 (BAP1) is a deubiquitinating hydrolase that plays a key role in various cellular processes. Germline and somatic inactivation of BRCA1 associated protein 1 gene (BAP1) is frequent in pleural mesothelioma; however, little is known about its status in peritoneal mesothelioma. METHODS: Taking advantage of the extensive French National Network for the Diagnosis of Malignant Pleural Mesothelioma and Rare Peritoneal Tumors and the French National Network for the Treatment of Rare Peritoneal Surface Malignancies, we collected biological material and clinical and epidemiological data for 46 patients with peritoneal mesothelioma. The status of BAP1 was evaluated at the mutational and protein expression levels and combined with our previous data on copy number alterations assessed in the same samples. RESULTS: We detected mutations in 32% of the malignant peritoneal mesotheliomas analyzed. In addition, we have previously reported that copy number losses occurred in 42% of the samples included in this series. Overall, 73% of the malignant peritoneal mesotheliomas analyzed carried at least one inactivated BAP1 allele, but only 57% had a complete loss of its protein nuclear expression. Better overall survival was observed for patients with BAP1 mutations (p = 0.04), protein expression loss (p = 0.016), or at least one of these alterations (p = 0.007) independently of tumor histological subtype, age, and sex. CONCLUSIONS: As in pleural mesothelioma, inactivation of BAP1 is frequent in peritoneal mesotheliomas. We found that BAP1 protein nuclear expression is a good prognostic factor and a more reliable marker for the complete loss of BAP1 activity than mutation or copy number loss.

Identification of Circulating Tumor DNA for the Early Detection of Small-cell Lung Cancer.

Circulating tumor DNA (ctDNA) is emerging as a key potential biomarker for post-diagnosis surveillance but it may also play a crucial role in the detection of pre-clinical cancer. Small-cell lung cancer (SCLC) is an excellent candidate for early detection given there are no successful therapeutic options for late-stage disease, and it displays almost universal inactivation of TP53. We assessed the presence of TP53 mutations in the cell-free DNA (cfDNA) extracted from the plasma of 51 SCLC cases and 123 non-cancer controls. We identified mutations using a pipeline specifically designed to accurately detect variants at very low fractions. We detected TP53 mutations in the cfDNA of 49% SCLC patients and 11.4% of non-cancer controls. When stratifying the 51 initial SCLC cases by stage, TP53 mutations were detected in the cfDNA of 35.7% early-stage and 54.1% late-stage SCLC patients. The results in the controls were further replicated in 10.8% of an independent series of 102 non-cancer controls. The detection of TP53 mutations in 11% of the 225 non-cancer controls suggests that somatic mutations in cfDNA among individuals without any cancer diagnosis is a common occurrence, and poses serious challenges for the development of ctDNA screening tests.