The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

Round Table Discussion on Optimal Clinical Trial Design in Precursor Multiple Myeloma

SUMMARY

While the current approach to precursor hematologic conditions is to "watch and wait," this may change with the development of therapies that are safe and extend survival or delay the onset of symptomatic disease. The goal of future therapies in precursor hematologic conditions is to improve survival and prevent or delay the development of symptomatic disease while maximizing safety. Clinical trial considerations in this field include identifying an appropriate at-risk population, safety assessments, dose selection, primary and secondary trial endpoints including surrogate endpoints, control arms, and quality-of-life metrics, all of which may enable more precise benefit-risk assessment.

Genomic Profiling to Contextualize the Results of Intervention for Smoldering Multiple Myeloma

PURPOSE

Early intervention for High-Risk Smoldering Multiple Myeloma (HR-SMM) achieves deep and prolonged responses. It is unclear if beneficial outcomes are due to treatment of less complex, susceptible disease or inaccuracy in clinical definition of cases entered.

EXPERIMENTAL DESIGN

Here, we interrogated whole genome and whole exome sequencing for 54 patients across two HR-SMM interventional studies (NCT01572480, NCT02279394).

RESULTS

We reveal that the genomic landscape of treated HR-SMM is generally simple as compared to Newly Diagnosed (ND)MM counterparts with less inactivation of tumor suppressor genes, RAS pathway mutations, MYC disruption, and APOBEC contribution. The absence of these events parallels that of indolent precursor conditions, possibly explaining overall excellent outcomes. However, some patients harboring genomic complexity fail to sustain response and experience resistant, progressive disease. Overall, clinical risk scores do not effectively discriminate between genomically indolent and aggressive disease.

CONCLUSIONS

Genomic profiling can contextualize the advantage of early intervention in SMM and guide personalization of therapy.

Genomic Classification and Individualized Prognosis in Multiple Myeloma

PURPOSE

Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years.

METHODS

To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data.

RESULTS

Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited.

CONCLUSION

Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

Molecular Pathogenesis of Multiple Myeloma: Clinical Implications

Multiple myeloma is a malignancy of bone-marrow-localized, isotype-switched plasma cells that secrete a monoclonal immunoglobulin and cause hyperCalcemia, Anemia, Renal failure, and lytic Bone disease. It is preceded, often for decades, by a relatively stable monoclonal gammopathy lacking these clinical and malignant features. Both conditions are characterized by the presence of types of immunoglobulin heavy gene translocations that dysregulate a cyclin D family gene on 11q13 (CCND1), 6p21 (CCND3), or 12q11 (CCND2), a maf family gene on 16q23 (MAF), 20q11 (MAFB), or 8q24 (MAFA), or NSD2/FGFR3 on 4p16, or the presence of hyperdiploidy. Subsequent loss of function of tumor suppressor genes and mutations activating MYC, RAS, NFkB, and cell cycle pathways are associated with the progression to malignant disease.

Childhood cancer mutagenesis caused by transposase-derived PGBD5

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

Emerging Strategies for the Prevention of Immune Toxicities Associated with T cell-Engaging Cancer Therapies

SUMMARY

Immune-related toxicities including cytokine release syndrome (CRS) and immune effector cell-associated neurotoxicity syndrome (ICANS) are common side effects of bispecific antibody and chimeric antigen receptor (CAR) T-cell therapies of hematologic malignancies. As anti-inflammatory therapy (the standard of care) is variably effective in mitigating these toxicities after onset, here we discuss emerging evidence for shifting the strategy from mitigation to prevention.

Impact of Rare Structural Variant Events in Newly Diagnosed Multiple Myeloma

PURPOSE

Whole-genome sequencing (WGS) of patients with newly diagnosed multiple myeloma (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e., hotspots) and causing recurrent copy-number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs." More than half of SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined.

EXPERIMENTAL DESIGN

In this study, we utilize 752 WGS and 591 RNA sequencing data from patients with NDMM to determine the role of rare SVs in myeloma pathogenesis.

RESULTS

Ninety-four percent of patients harbored at least one rare SV event. Rare SVs showed an SV class-specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs.

CONCLUSIONS

Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a potential driver role in myeloma pathogenesis.

Genomic and immune determinants of resistance to anti-CD38 monoclonal antibody-based therapy in relapsed refractory multiple myeloma

Anti-CD38 antibody therapies have transformed multiple myeloma (MM) treatment. However, a large fraction of patients inevitably relapses. To understand this, we investigated 32 relapsed MM patients treated with daratumumab, lenalidomide, and dexamethasone (Dara-Rd; NCT03848676 ). Whole genome sequencing (WGS) before and after treatment pinpointed genomic drivers associated with early progression, including RPL5 loss and APOBEC mutagenesis. Flow cytometry on 202 blood samples, collected every three months until progression for 31 patients, revealed distinct immune changes significantly impacting clinical outcomes. Progressing patients exhibited significant depletion of CD38+ NK cells, persistence of T cell exhaustion, and reduced depletion of T-reg cells over time. These findings underscore the influence of immune composition and daratumumab-induced immune changes in promoting MM resistance. Integrating genomics and flow cytometry unveiled associations between adverse genomic features and immune patterns. Overall, this study sheds light on the intricate interplay between genomic complexity and the immune microenvironment driving resistance to Dara-Rd.

Genomic and immune signatures predict clinical outcome in newly diagnosed multiple myeloma treated with immunotherapy regimens

Despite improving outcomes, 40% of patients with newly diagnosed multiple myeloma treated with regimens containing daratumumab, a CD38-targeted monoclonal antibody, progress prematurely. By integrating tumor whole-genome and microenvironment single-cell RNA sequencing from upfront phase 2 trials using carfilzomib, lenalidomide and dexamethasone with daratumumab ( NCT03290950 ), we show how distinct genomic drivers including high APOBEC mutational activity, IKZF3 and RPL5 deletions and 8q gain affect clinical outcomes. Furthermore, evaluation of paired bone marrow profiles, taken before and after eight cycles of carfilzomib, lenalidomide and dexamethasone with daratumumab, shows that numbers of natural killer cells before treatment, high T cell receptor diversity before treatment, the disappearance of sustained immune activation (that is, B cells and T cells) and monocyte expansion over time are all predictive of sustained minimal residual disease negativity. Overall, this study provides strong evidence of a complex interplay between tumor cells and the immune microenvironment that is predictive of clinical outcome and depth of treatment response in patients with newly diagnosed multiple myeloma treated with highly effective combinations containing anti-CD38 antibodies.

Scientific techniques in adolescent and young adult classic Hodgkin lymphoma

Understanding the tumor microenvironment and genomic landscape is crucial for better prediction of treatment outcomes and developing novel therapies in Hodgkin lymphoma (HL). Recent advancements in genomics have enabled researchers to gain deeper insights into the genomic characteristics of HL at both single-cell resolution and the whole genome level. The use of noninvasive methods such as liquid biopsies and formalin-fixed paraffin-embedded-based imaging techniques has expanded the possibilities of applying cutting-edge analyses to routine clinically available samples. Collaborative efforts between adult and pediatric group are imperative to translate novel findings into routine patient care.

Multiomic Mapping of Acquired Chromosome 1 Copy-Number and Structural Variants to Identify Therapeutic Vulnerabilities in Multiple Myeloma

PURPOSE

Chromosome 1 (chr1) copy-number abnormalities (CNA) and structural variants (SV) are frequent in newly diagnosed multiple myeloma (NDMM) and are associated with a heterogeneous impact on outcomes, the drivers of which are largely unknown.

EXPERIMENTAL DESIGN

A multiomic approach comprising CRISPR, gene mapping of CNAs and SVs, methylation, expression, and mutational analysis was used to document the extent of chr1 molecular variants and their impact on pathway utilization.

RESULTS

We identified two distinct groups of gain(1q): focal gains associated with limited gene-expression changes and a neutral prognosis, and whole-arm gains, which are associated with substantial gene-expression changes, complex genetics, and an adverse prognosis. CRISPR identified a number of dependencies on chr1 but only limited variants associated with acquired CNAs. We identified seven regions of deletion, nine of gain, three of chromothripsis (CT), and two of templated insertion (TI), which contain a number of potential drivers. An additional mechanism involving hypomethylation of genes at 1q may contribute to the aberrant gene expression of a number of genes. Expression changes associated with whole-arm gains were substantial and gene set enrichment analysis identified metabolic processes, apoptotic resistance, signaling via the MAPK pathway, and upregulation of transcription factors as being key drivers of the adverse prognosis associated with these variants.

CONCLUSIONS

Multiple layers of genetic complexity impact the phenotype associated with CNAs on chr1 to generate its associated clinical phenotype. Whole-arm gains of 1q are the critically important prognostic group that deregulate multiple pathways, which may offer therapeutic vulnerabilities.

Immunophenotypic correlates of sustained MRD negativity in patients with multiple myeloma

The role of the immune microenvironment in maintaining disease remission in patients with multiple myeloma (MM) is not well understood. In this study, we comprehensively profile the immune system in patients with newly diagnosed MM receiving continuous lenalidomide maintenance therapy with the aim of discovering correlates of long-term treatment response. Leveraging single-cell RNA sequencing and T cell receptor β sequencing of the peripheral blood and CyTOF mass cytometry of the bone marrow, we longitudinally characterize the immune landscape in 23 patients before and one year after lenalidomide exposure. We compare patients achieving sustained minimal residual disease (MRD) negativity to patients who never achieved or were unable to maintain MRD negativity. We observe that the composition of the immune microenvironment in both the blood and the marrow varied substantially according to both MRD negative status and history of autologous stem cell transplant, supporting the hypothesis that the immune microenvironment influences the depth and duration of treatment response.

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 Vk*MYC Mouse Model recapitulates human multiple myeloma evolution and genomic diversity

Despite advancements in profiling multiple myeloma (MM) and its precursor conditions, there is limited information on mechanisms underlying disease progression. Clincal efforts designed to deconvolute such mechanisms are challenged by the long lead time between monoclonal gammopathy and its transformation to MM. MM mouse models represent an opportunity to overcome this temporal limitation. Here, we profile the genomic landscape of 118 genetically engineered Vk*MYC MM and reveal that it recapitulates the genomic heterogenenity and life history of human MM. We observed recurrent copy number alterations, structural variations, chromothripsis, driver mutations, APOBEC mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identified frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC expression, that drives the progression of monoclonal gammopathy to MM.

Tracking the evolution of therapy-related myeloid neoplasms using chemotherapy signatures

Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Preleukemic clones (ie, clonal hematopoiesis [CH]) are detectable years before the development of these aggressive malignancies, although the genomic events leading to transformation and expansion are not well defined. Here, by leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of prechemotherapy samples, we reconstructed the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy was revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan were hypermutated and enriched for complex structural variants (ie, chromothripsis), whereas neoplasms with nonmutagenic chemotherapy exposures were genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to discrete clinical exposure in each patient's life, we estimated that several complex events and genomic drivers were acquired after chemotherapy was administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes melphalan exposure and is selected after reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA damage. Overall, we revealed a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select preexisting CH but also promote the acquisition of recurrent genomic drivers.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis.

SIGNIFICANCE

Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

Abstract 6067: Characterizing the landscape of rare structural variant events in newly diagnosed multiple myeloma

Background: Structural variants (SV) are known to play a critical role in the pathogenesis of multiple cancer types. Using whole genome sequencing (WGS), we recently characterized the SV landscape of 752 newly diagnosed multiple myeloma (MM) patients, identifying 68 SV hotspots and 152 recurrent copy number aberrations (CNA; Rustad et al. Blood Cancer Discovery 2020). Despite comprehensive annotation, more than half of SVs were not linked to any known MM genomic driver. The biological impact of these SV events, here defined as rare SV, occurring in 93% (702/752) of patients, is unknown.

Methods: To study the biological impact of rare SVs, we interrogated WGS (n=752) and RNAseq (n=591) in the CoMMpass trial. Recurrent SVs identified by involvement in canonical Ig translocations, recurrent CNAs, or SV hotspots were excluded. All SVs within an event must not involve a recurrent region to be defined rare. To determine SV class-specific gene relationships, breakpoint enrichment was compared against a permuted background model for each SV class and gene expression direction, up to 1 Mb. Genes were considered affected if expression was above a gene specific outlier Z-score of +/- 2. Lastly, we modeled breakpoint density to the nearest MM superenhancer up to 10 MB, and compared to permuted background rates.

Results: Of the total 8,942 SVs, 4,959 (55%) were identified as rare. 201 (34%) patients had at least 1 enriched rare SV event associated with gene expression outliers. Amongst over-expressed gene outliers, rare templated insertions and duplications were enriched within the gene body and up to 1 MB away. Rare inversions were enriched in genes 100kb and 1MB away, and rare translocations were associated with outliers 1 MB away. Amongst under-expressed gene outliers, rare complex SVs were enriched within the gene body, while deletion SVs were enriched in the gene body and up to 1 MB away. Rare duplications, translocations and templated insertions were enriched up to 1 MB of superenhancers. Rare templated insertions were significantly enriched against the background model (p < 0.001, Fisher Exact). Overall, 82% (104/126) of gene outliers affected by rare templated insertions were associated with superenhancers, (95 over-expressed, e.g. IRF6 and 9 under-expressed), 54% (130/237) by rare translocations, (105 overexpressed, e.g. FAM46A and 25 under-expressed), and 55% (96/172) by rare duplication events, (93 overexpressed, e.g. MAPK13, and 7 under-expressed). In addition, among the 853 outlier genes affected by enriched rare SVs, at least 15 are involved with B-cell development, suggesting a potential driver role in myeloma pathogenesis.

Conclusion: In summary, leveraging WGS and RNA-seq of clinical samples, we demonstrate that rare SVs are frequently associated with aberrant gene expression, expanding our understanding of their potential role in heterogenous clinical response in patients diagnosed with MM.

Citation Format: Monika Chojnacka, Benjamin Diamond, Bachisio Ziccheddu, Even Rustad, Kylee Maclachlan, Marios Papadimitriou, Eileen Boyle, Patrick Blaney, Saad Usmani, Gareth Morgan, Ola Landgren, Francesco Maura. Characterizing the landscape of rare structural variant events in newly diagnosed multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6067.

Waldenström macroglobulinemia whole genome reveals prolonged germinal center activity and late copy number aberrations

The genomic landscape of Waldenström macroglobulinemia (WM) is characterized by somatic mutations in MYD88, present from the precursor stages. Using the comprehensive resolution of whole genome sequencing (WGS) in 14 CD19-selected primary WM samples; comparing clonal and subclonal mutations revealed that germinal center (GC) mutational signatures SBS9 (poly-eta) and SBS84 (AID) have sustained activity, suggesting that the interaction between WM and the GC continues over time. Expanding our cohort size with 33 targeted sequencing samples, we interrogated the WM copy number aberration (CNA) landscape and chronology. Of interest, CNA prevalence progressively increased in symptomatic WM and relapsed disease when compared with stable precursor stages, with stable precursors lacking genomic complexity. Two MYD88 wild-type WGS contained a clonal gain affecting chromosome 12, which is typically an early event in chronic lymphocytic leukemia. Molecular time analysis demonstrated that both chromosomal 12 gain events occurred early in cancer development whereas other CNA changes tend to occur later in the disease course and are often subclonal. In summary, WGS analysis in WM allows the demonstration of sustained GC activity over time and allows the reconstruction of the temporal evolution of specific genomic features. In addition, our data suggest that, although MYD88-mutations are central to WM clone establishment and can be observed in precursor disease, CNA may contribute to later phases, and may be used as a biomarker for progression risk from precursor conditions to symptomatic disease.

Impact of rare structural variant events in newly diagnosed multiple myeloma

Whole genome sequencing (WGS) of newly diagnosed multiple myeloma patients (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e. hotspots) and causing recurrent copy number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs". More than half SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined. In this study, we utilize 752 WGS and 591 RNA-seq data from NDMM patients to determine the role of rare SVs in myeloma pathogenesis. 94% of patients harbored at least one rare SV event. Rare SVs showed an SV-class specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs. Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a driver role in myeloma pathogenesis.

SIGNIFICANCE

Characterization of multiple myeloma genome revealed that more than half structural variants are not involved in recurrent events. Here, we demonstrate that these rare SVs hold potential for myeloma pathogenesis through their gene expression impact. Rare SVs contribute to MM heterogeneity and have implications for development of individualized treatment.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)^high–B cell receptor (BCR)^low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Single-cell genomic and transcriptomic analyses of longitudinal samples of patients with Richter syndrome reveal the presence and dynamics of clones driving transformation from chronic lymphocytic leukemia years before clinical manifestation

Whole-genome sequencing reveals complex genomic features underlying anti-CD19 CAR T-cell treatment failures in lymphoma

CD19-directed chimeric antigen receptor (CAR-19) T cells are groundbreaking immunotherapies approved for use against large B-cell lymphomas. Although host inflammatory and tumor microenvironmental markers associate with efficacy and resistance, the tumor-intrinsic alterations underlying these phenomena remain undefined. CD19 mutations associate with resistance but are uncommon, and most patients with relapsed disease retain expression of the wild-type receptor, implicating other genomic mechanisms. We therefore leveraged the comprehensive resolution of whole-genome sequencing to assess 51 tumor samples from 49 patients with CAR-19-treated large B-cell lymphoma. We found that the pretreatment presence of complex structural variants, APOBEC mutational signatures, and genomic damage from reactive oxygen species predict CAR-19 resistance. In addition, the recurrent 3p21.31 chromosomal deletion containing the RHOA tumor suppressor was strongly enriched in patients for whom CAR T-cell therapy failed. Pretreatment reduced expression or monoallelic loss of CD19 did not affect responses, suggesting CAR-19 therapy success and resistance are related to multiple mechanisms. Our study showed that tumor-intrinsic genomic alterations are key among the complex interplay of factors that underlie CAR-19 efficacy and resistance for large B-cell lymphomas.

The spatio-temporal evolution of multiple myeloma from baseline to relapse-refractory states

Deciphering Multiple Myeloma evolution in the whole bone marrow is key to inform curative strategies. Here, we perform spatial-longitudinal whole-exome sequencing, including 140 samples collected from 24 Multiple Myeloma patients during up to 14 years. Applying imaging-guided sampling we observe three evolutionary patterns, including relapse driven by a single-cell expansion, competing/co-existing sub-clones, and unique sub-clones at distinct locations. While we do not find the unique relapse sub-clone in the baseline focal lesion(s), we show a close phylogenetic relationship between baseline focal lesions and relapse disease, highlighting focal lesions as hotspots of tumor evolution. In patients with ≥3 focal lesions on positron-emission-tomography at diagnosis, relapse is driven by multiple distinct sub-clones, whereas in other patients, a single-cell expansion is typically seen (p < 0.01). Notably, we observe resistant sub-clones that can be hidden over years, suggesting that a prerequisite for curative therapies would be to overcome not only tumor heterogeneity but also dormancy.

Diverse mutational landscapes in human lymphocytes

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Abstract 4067: Combination venetoclax and selinexor effective in relapsed refractory multiple myeloma with translocation t(11;14)

Translocation t(11;14) positive multiple myeloma (MM) is sensitive to venetoclax, yet the drug lacks FDA approval in MM. Selinexor is an inhibitor of nuclear exporter XPO1 and is approved in relapsed refractory MM. We recently analyzed RNA-seq from 859 MM patient samples (MMRF CoMMpass study) and compared them to healthy bone marrow samples (GSE114922). After correcting for batch effects, differential gene expression analysis found apoptosis and cell cycle regulation amongst the top upregulated gene pathways in MM. We therefore hypothesized that MM would be sensitive to the combination of the pro-apoptotic drug, venetoclax, and selinexor, which causes cell cycle arrest.

Using a panel of 8 MM cell lines (U266-B1, KMS12BM, SKMM2, RPMI-8226, LP1, OPM2, KMS12PE, and NCIH929), we performed cell viability assays in increasing concentrations of selinexor, venetoclax, and a combination of the two drugs at 72 hours. Overall, the combination showed synergy in MM cell lines but strikingly, cell lines that possessed t(11;14) were more sensitive to the drug combination and showed enhanced synergy compared to the non-t(11;14). Average Bliss model synergy score were 0.5 in non-t(11;14) and 10.2 in t(11;14) MM cells (&gt;10 indicates synergistic effects and &lt;-10 indicates antagonistic drug effects).

To better understand molecular mechanisms underlying the observed synergistic effect in t(11;14), we performed western blot analysis in these cell lines treated with selinexor, venetoclax, the combination, or DMSO control for 24 hours. We measured protein levels of Cyclin D1, which is overexpressed in t(11;14) and a cargo of XPO1. We confirmed Cyclin D1 overexpression in t(11;14) cells lines but not in non-t(11;14) cells. Cyclin D1 levels decreased with selinexor, and the reduction was enhanced by adding venetoclax. Additionally, we measured levels of XPO1, p53, MCL-1, and p65, which we have previously shown to be altered by selinexor treatment. While venetoclax monotherapy caused increased MCL-1 expression, this was abrogated by combining with selinexor and the effects were statistically more significant in t(11;14) cell lines.

Concurrently, two patients with relapsed, refractory t(11;14) MM and disease progression after multiple lines of therapy were treated with venetoclax; however, after initial response, they developed resistance and progressive disease. Based off our preclinical results, selinexor was added to venetoclax and recaptured responses (VGPR and MR, respectively) suggesting a beneficial effect of the combination over single agent venetoclax. To our knowledge, this is the first report of patients with MM treated with the combination of selinexor and venetoclax. Further clinical evaluation of this combination in t(11;14) MM is planned and samples will be collected to better understand the mechanistic correlation between XPO1, BCL2/MCL-1, and Cyclin D1 levels and response.

Citation Format: Sana Chaudhry, Nina Nguyen, Rimpi Khurana, Tulasigeri M. Totiger, Skye Montoya, Jumana Afaghani, Terrence Bradley, Francesco Maura, Dickran Kazandjian, Jennifer Chapman, Stephan Schürer, Ola Landgren, James Hoffman, Justin Taylor. Combination venetoclax and selinexor effective in relapsed refractory multiple myeloma with translocation t(11;14) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 4067.

Abstract 3795: Early seeding of Richter transformation in chronic lymphocytic leukemia

Introduction: Clonal evolution drives cancer development due to the emergence and/or selection of proliferatively advantageous subclones. Its understanding may facilitate the design of anticipation-based management strategies. Richter transformation (RT) is a paradigmatic tumor evolution in which chronic lymphocytic leukemia (CLL), an indolent neoplasia of mature B-cells, transforms into a high-grade lymphoma, usually diffuse large B-cell lymphoma (DLBCL), conferring a dismal prognosis. The evolutionary trajectories of RT and its driving (epi)genomic mechanisms remain largely unknown.

Aims: To reconstruct the evolutionary history of RT and to reveal the molecular processes underlying this transformation.

Methods: We characterized the whole genome (WGS), epigenome (DNA methylation, H3K27ac, ATAC-seq), and transcriptome (RNA-seq), combined with single-cell DNA and RNA sequencing analyses, of 19 CLL patients developing RT before (n=3) or after treatment with chemoimmunotherapy (n=6) and targeted therapies (BCR or BCL2 inhibitors, n=10). We analyzed 54 longitudinal samples covering up to 19 years of disease course.

Results: Our WGS analyses uncovered that RT is characterized by a remarkable structural complexity. We also identified a novel treatment-independent RT-specific mutational process, which we named SBS-RT. The genetic driver landscape of RT is a compendium of alterations in genes involved in cell cycle, MYC, and NF-κB pathways, frequently targeted in single catastrophic events including chromothripsis and chromoplexy. The WGS-based phylogenic reconstruction and single-cell DNA/RNA-seq analyses identified a very early diversification of CLL leading to emergence of RT-cells carrying specific genetic drivers and transcriptomic profiles of RT already at CLL diagnosis. These small subclones were dormant for 6-19 years until rapid expansion associated with the clinical transformation. While the DNA methylome kept track of the cell of origin and proliferative history of RT cells, their chromatin configuration and transcriptional program converged into the overexpression of cell cycle regulators, Toll-like receptors, MYC, MTORC1, and OXPHOS related transcripts, as well as downregulation of BCR pathway. This phenotypic shift was related to de novo activation of key transcription factors. In vitro experiments confirmed that RT cells have a 4-fold higher oxygen consumption at routine respiration and electron transfer system capacity compared to CLL. The resistance of RT to BCR inhibition is consistent with its high OXPHOS and low BCR signaling, which mimics de novo DLBCL-OXPHOS insensitive to BCR inhibition. This OXPHOShigh-BCRlow transcriptional axis of RT can be exploited therapeutically.

Conclusions: These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover therapeutic targets for the, once expanded, lethal Richter transformation.

Citation Format: Ferran Nadeu, Romina Royo, Ramon Massoni-Badosa, Beatriz Garcia-Torre, Martí Duran-Ferrer, Kevin J. Dawson, Marta Kulis, Ander Diaz-Navarro, Neus Villamor, Juan L. Melero, Vicente Chapaprieta, Ana Dueso-Barroso, Julio Delgado, Riccardo Moia, Sara Ruiz-Gil, Domenica Marchese, Núria Verdaguer-Dot, Mónica Romo, Maria Rozman, Gerard Frigola, Alfredo Rivas-Delgado, Tycho Baumann, Miguel Alcoceba, Marcos González, Fina Climent, Pau Abrisqueta, Josep Castellví, Francesc Bosch, Marta Aymerich, Anna Enjuanes, Sílvia Ruiz-Gaspà, Armando López-Guillermo, Pedro Jares, Sílvia Beà, Dolors Colomer, Núria López-Bigas, Josep LlGelpí, David Torrents, Peter J. Campbell, Ivo Gut, Pablo M. Garcia-Roves, Davide Rossi, Gianluca Gaidano, Xose S. Puente, Holger Heyn, Francesco Maura, José I. Martín-Subero, Elías Campo. Early seeding of Richter transformation in chronic lymphocytic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3795.

Abstract 5747: Chemotherapy-related mutational signatures reveal the origins of therapy-related myeloid neoplasms

Patients treated with chemotherapy (CT) and/or autologous stem-cell transplantation (ASCT) are at risk for therapy-related myeloid neoplasms (tMN). Certain cytotoxic agents introduce mutations within distinct trinucleotide contexts resulting in a unique barcode for each exposed cell. We leveraged mutational signatures to investigate the role of CT in the genomic landscape of tMN with respect to antecedent clonal hematopoiesis (CH).

We analyzed 32 tMN and 2 tALL from 33 patients and interrogated for copy number abnormalities (CNA), structural variants (SV), single nucleotide variants (SNV), and mutational signatures. For 7 patients with tMN post-melphalan/ASCT, we investigated antecedent CH using targeted sequencing on pre-melphalan samples, including autograft products.

CH variants that became clonal in tumor were seen in 5/7 pre-melphalan/ASCT samples (TP53, RUNX1, NCOR1, NF1, CREBBP, DNMT3A, and PPM1D). Complex SV were seen in 7 tMNs; including chromothripsis in 6 (19.4%). In 4 cases, chromothripsis involved chromosome 19 with hyper-amplification of the SMARCA4 locus (≥5 copies). Mutational signature analysis revealed 6 known single base substitution (SBS) signatures in tMN including melphalan (SBS-MM1) and platinum signatures (SBS31, SBS35, and E-SBS37). TMNs with CT signatures had higher mutation burden than those without (p = 0.004). 17 patients with exposure to agents other than melphalan/platinum did not have increased mutational burden with respect to de novo AML (TCGA; NEJM, 2013).

All patients with prior platinum exposure (including tALL, n=9) had platinum SBS signatures while only 2 of 7 patients with prior melphalan/ASCT had a melphalan signature (SBS-MM1). Detection of CT signatures in bulk sequencing relies on one cell, with its barcode of mutations, to expand to clonal dominance. Given pre-existent CH, including in 3/3 autograft products, absence of a CT signature despite melphalan exposure implies progression by a clone that escaped CT exposure with stem-cell collection and reinfusion. Conversely, all platinum-exposed tAML had signature evidence of exposure confirming existence of CH prior to exposure and supporting post-CT single-cell expansion. TMNs from 3 patients exposed to sequential platinum and melphalan/ASCT had platinum but not melphalan signatures confirming single-cell expansion of the pre-tMN CH clone post-platinum but with escape from exposure to melphalan via leukapheresis. Chromothripsis events bore only non-duplicated CT-induced mutations, indicative of acquisition prior to, and not directly caused by, CT exposure.

These disparities suggest that ASCT provides a mechanism for CH clones to escape CT and re-engraft with transplant. Coupled with driver events accrued prior to CT, this suggest that CT-induced mutagenesis may be less important than other factors, such as CT-induced immunosuppression, in the expansion of pre-TMN CH clones.

Citation Format: Benjamin Diamond, Bachisio Ziccheddu, Eileen M. Boyle, Kylee Maclachlan, Justin Taylor, Justin M. Watts, Sydney X. Lu, David G. Coffey, Niccolo Bolli, Elli Papaemmanuil, Kelly Bolton, Jae H. Park, Heather Landau, Karuna Ganesh, Mikkael A. Sekeres, Stephen Nimer, David J. Chung, Caleb H. Ho, Mikhail Roshal, Alexander Lesokhin, Gareth Morgan, Ola Landgren, Francesco Maura. Chemotherapy-related mutational signatures reveal the origins of therapy-related myeloid neoplasms [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5747.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

Bayesian networks elucidate complex genomic landscapes in cancer

Bayesian networks (BNs) are disciplined, explainable Artificial Intelligence models that can describe structured joint probability spaces. In the context of understanding complex relations between a number of variables in biological settings, they can be constructed from observed data and can provide a guiding, graphical tool in exploring such relations. Here we propose BNs for elucidating the relations between driver events in large cancer genomic datasets. We present a methodology that is specifically tailored to biologists and clinicians as they are the main producers of such datasets. We achieve this by using an optimal BN learning algorithm based on well established likelihood functions and by utilising just two tuning parameters, both of which are easy to set and have intuitive readings. To enhance value to clinicians, we introduce (a) the use of heatmaps for families in each network, and (b) visualising pairwise co-occurrence statistics on the network. For binary data, an optional step of fitting logic gates can be employed. We show how our methodology enhances pairwise testing and how biologists and clinicians can use BNs for discussing the main relations among driver events in large genomic cohorts. We demonstrate the utility of our methodology by applying it to 5 cancer datasets revealing complex genomic landscapes. Our networks identify central patterns in all datasets including a central 4-way mutual exclusivity between HDR, t(4,14), t(11,14) and t(14,16) in myeloma, and a 3-way mutual exclusivity of three major players: CALR, JAK2 and MPL, in myeloproliferative neoplasms. These analyses demonstrate that our methodology can play a central role in the study of large genomic cancer datasets.

Bayesian network inference on several blood and solid cancer genomic datasets provides more accessible ways to explore driver events in cancer.

Defining genomic events involved in the evolutionary trajectories of myeloma and its precursor conditions

All patients with a diagnosis of multiple myeloma (MM) have a preceding, asymptomatic expansion of clonal plasma cells, clinically recognized as monoclonal gammopathy of undetermined significance or smoldering multiple myeloma (SMM). While most patients with monoclonal gammopathy of undetermined significance have a very small rate of progression, SMM is a widely heterogeneous condition where a fraction of patients will progress to symptomatic MM rather quickly, while others will experience an indolent clinical course. The differentiation between progressive and stable precursor condition thus represents one of the most important unmet clinical needs in the MM community. The ability to identify patients at high-risk of progression before major clonal expansion and onset of end-organ damage would enable strategies for early prevention and perhaps more effective intervention. All proposed criteria to predict the progression of myeloma precursor conditions are built around indirect markers of disease burden and, therefore, are generally able to accurately identify only a small fraction of patients in whom progression to MM is already occurring. Leveraging whole genome and exome sequencing, it has been shown that patients with stable myeloma precursor conditions are characterized by either absence or lower prevalence of distinct genomic events that are detectable in progressive precursor condition years before the progression. In this review, we discuss evolving genomic concepts and tools; and their ability to differentiate myeloma precursor conditions into two distinct entities: one benign (monoclonal gammopathy of benign significance) and another malignant (asymptomatic multiple myeloma).

Functional Impact of Genomic Complexity on the Transcriptome of Multiple Myeloma

PURPOSE

Multiple myeloma is a biologically heterogenous plasma-cell disorder. In this study, we aimed at dissecting the functional impact on transcriptome of gene mutations, copy-number abnormalities (CNA), and chromosomal rearrangements (CR). Moreover, we applied a geno-transcriptomic approach to identify specific biomarkers for personalized treatments.

EXPERIMENTAL DESIGN

We analyzed 514 newly diagnosed patients from the IA12 release of the CoMMpass study, accounting for mutations in multiple myeloma driver genes, structural variants, copy-number segments, and raw-transcript counts. We performed an in silico drug sensitivity screen (DSS), interrogating the Cancer Dependency Map (DepMap) dataset after anchoring cell lines to primary tumor samples using the Celligner algorithm.

RESULTS

Immunoglobulin translocations, hyperdiploidy and chr(1q)gain/amps were associated with the highest number of deregulated genes. Other CNAs and specific gene mutations had a lower but very distinct impact affecting specific pathways. Many recurrent genes showed a hotspot (HS)-specific effect. The clinical relevance of double-hit multiple myeloma found strong biological bases in our analysis. Biallelic deletions of tumor suppressors and chr(1q)-amplifications showed the greatest impact on gene expression, deregulating pathways related to cell cycle, proliferation, and expression of immunotherapy targets. Moreover, our in silico DSS showed that not only t(11;14) but also chr(1q)gain/amps and CYLD inactivation predicted differential expression of transcripts of the BCL2 axis and response to venetoclax.

CONCLUSIONS

The multiple myeloma genomic architecture and transcriptome have a strict connection, led by CNAs and CRs. Gene mutations impacted especially with HS-mutations of oncogenes and biallelic tumor suppressor gene inactivation. Finally, a comprehensive geno-transcriptomic analysis allows the identification of specific deregulated pathways and candidate biomarkers for personalized treatments in multiple myeloma.

675 Genomic drivers of large B-cell lymphoma resistance to CD19 CAR-T therapy

Background

CD19-directed chimeric antigen receptor-reprogrammed autologous T cells are breakthrough immunotherapies for heavily pretreated patients with diffuse large B-cell lymphoma (DLBCL), but across CAR-19 products, ~60% of patients fail to respond or relapse. Inflammatory markers and clinical factors associate with impaired responses, but tumor-intrinsic resistance drivers are largely undefined.

Methods

To characterize the genomic mechanisms involved resistance to CAR-19, we interrogated whole genome sequencing (WGS) from 28 relapsed/refractory (r/r) aggressive lymphoma patients uniformly treated with axicabtagene ciloleucel (axi-cel).

Results

Because prognostic factors defined in the frontline treatment setting are largely inapplicable to CAR-19, we leveraged the WGS data, including comparative analyses with untreated DLBCL cases in the Pan-Cancer Analysis of Whole Genomes (PCAWG) (figure 1). In analyses of individual mutated genes, TP53 was significantly enriched (p=0.002) in CAR-19 patients, but did not predict outcome. However, mutations in either NFKBIA or MYC associated with worse PFS after CAR-19 (p=0.04, p=0.025 respectively). We next identified 12 single base substitution (SBS) mutational signatures in our cohort and found presence of APOBEC (SBS2 and SBS13) signatures associated with worse PFS, with 4/5 patients progressing (p=0.03). Copy number analysis by GISTIC2.0 revealed focal deletions of RHOA and RB1 to be significantly enriched in our cohort and independently predicted poor outcome (p=0.0007, p=0.05 respectively). WGS identifies structural variants and complex events. We found chromothripsis, a catastrophic shattering and reassembly of chromosomes, in 39.3% of r/r DLBCL, which was strongly associated with poor CAR-19 outcome, with 9/11 affected cases progressing (p=0.041). Finally, reduced expression (n=3) or genomic alteration (n=3) of CD19 did not associate with poor outcome. One case with durable response contained a sub-clonal CD19 mutation (L174V) previously reported as associated with CAR-19 resistance. These findings demonstrate predominance of CD19-independent resistance and indicate antigen-mediated tumor killing is not the only mechanism of tumor eradication. Genomic complexity appears to promote an immunosuppressive tumor microenvironment (TME), limiting CAR-19 efficacy.

Conclusions

Leveraging the resolution of WGS, we observed that markers of genomic complexity (chromothripsis and APOBEC) and specific genomic alterations (RHOA and RB1 deletions) associate with resistance to CAR-19 immunotherapy for aggressive B-cell lymphomas (figure 1). 93.8% of CAR-19 relapsed patients contained at least one or these genomic alterations. Recent patient data demonstrate that an immunosuppressed TME leads to CAR-19 failure. Combining these findings with our genomics findings, successful CAR-19 therapy must overcome the immune-exhausted TME to mobilize the host immune system and eliminate the tumor.

Abstract 675 Figure 1

Genomic alterations associated with disease progression. (a) The heatmap shows the significant genomic alteration present in at least 4 patients associated with progression after CD19 CAR-T cell therapy. (b) Kaplan-Meier curve of progression free survival with the combination of statistically significant genomic anomalies

CDKN2A deletion is a frequent event associated with poor outcome in patients with peripheral T-cell lymphoma not otherwise specified (PTCL-NOS)

Nodal peripheral T-cell lymphoma not otherwise specified (PTCL-NOS) remains a diagnosis encompassing a heterogenous group of PTCL cases not fitting criteria for more homogeneous subtypes. They are characterized by a poor clinical outcome when treated with anthracycline-containing regimens. A better understanding of their biology could improve prognostic stratification and foster the development of novel therapeutic approaches. Recent targeted and whole exome sequencing studies have shown recurrent copy number abnormalities (CNAs) with prognostic significance. Here, investigating 5 formalin-fixed, paraffin embedded cases of PTCL-NOS by whole genome sequencing (WGS), we found a high prevalence of structural variants and complex events, such as chromothripsis likely responsible for the observed CNAs. Among them, CDKN2A and PTEN deletions emerged as the most frequent aberration, as confirmed in a final cohort of 143 patients with nodal PTCL. The incidence of CDKN2A and PTEN deletions among PTCL-NOS was 46% and 26%, respectively. Furthermore, we found that co-occurrence of CDKN2A and PTEN deletions is an event associated with PTCL-NOS with absolute specificity. In contrast, these deletions were rare and never co-occurred in angioimmunoblastic and anaplastic lymphomas. CDKN2A deletion was associated with shorter overall survival in multivariate analysis corrected by age, IPI, transplant eligibility and GATA3 expression (adjusted HR =2.53; 95% CI 1.006-6.3; p=0.048). These data suggest that CDKN2A deletions may be relevant for refining the prognosis of PTCL-NOS and their significance should be evaluated in prospective trials.

Copy number signatures predict chromothripsis and clinical outcomes in newly diagnosed multiple myeloma

Chromothripsis is detectable in 20–30% of newly diagnosed multiple myeloma (NDMM) patients and is emerging as a new independent adverse prognostic factor. In this study we interrogate 752 NDMM patients using whole genome sequencing (WGS) to investigate the relationship of copy number (CN) signatures to chromothripsis and show they are highly associated. CN signatures are highly predictive of the presence of chromothripsis (AUC = 0.90) and can be used identify its adverse prognostic impact. The ability of CN signatures to predict the presence of chromothripsis is confirmed in a validation series of WGS comprised of 235 hematological cancers (AUC = 0.97) and an independent series of 34 NDMM (AUC = 0.87). We show that CN signatures can also be derived from whole exome data (WES) and using 677 cases from the same series of NDMM, we are able to predict both the presence of chromothripsis (AUC = 0.82) and its adverse prognostic impact. CN signatures constitute a flexible tool to identify the presence of chromothripsis and is applicable to WES and WGS data.

Chromothripsis is associated with unfavourable outcomes in multiple myeloma (MM), but its detection usually requires whole genome sequencing. Here the authors develop an approach to detect chromothripsis in MM based on copy-number signatures that also works with whole exome sequencing data.

The mutagenic impact of melphalan in multiple myeloma

The introduction of whole genome and exome sequencing partnered with advanced bioinformatic pipelines has allowed the comprehensive characterization of mutational processes (i.e., mutational signatures) in individual cancer patients. Studies focusing on multiple myeloma have defined several mutational processes, including a recently identified mutational signature (called "SBS-MM1") directly caused by exposure to high-dose melphalan (i.e., autologous stem cell transplant). High-dose melphalan exposure increases both the overall and nonsynonymous mutational burden detected between diagnosis and relapse by ~10-20%. Nevertheless, most of these mutations are acquired within the heterochromatin and late-replicating regions, rarely involving key myeloma driver genes. In this review, we summarize key studies that made this discovery possible, and we discuss potential clinical implications.

Dynamics of minimal residual disease in patients with multiple myeloma on continuous lenalidomide maintenance: a single-arm, single-centre, phase 2 trial

Background Lenalidomide maintenance improves progression-free survival for patients with multiple myeloma, although its optimal duration is unknown. Clearance of minimal residual disease (MRD) in the bone marrow results in superior outcomes, although its attainment or sustainment does not alter clinical decision-making. Studies that have evaluated MRD serially are limited in length. We therefore aimed to evaluate longitudinal changes in MRD-status (dynamics) and their association with progression-free survival in patients with multiple myeloma.

METHODS

In this single-centre, single-arm, phase 2 study, we enrolled patients aged 18 years and older from the Memorial Sloan Kettering Cancer Center (New York, NY, USA) who had newly diagnosed multiple myeloma following unrestricted frontline therapy and an Eastern Cooperative Oncology Group Performance Status of 2 or lower, including patients who started maintenance before study enrolment. All participants received lenalidomide maintenance at 10 mg for 21 days of 28-day cycles until progression or unacceptable toxic effects for up to 5 years on protocol. The primary endpoint was progression-free survival at 60 months per protocol and key secondary endpoints were MRD rates after completion of the 12th, 24th, and 36th cycle of maintenance and the association between progression-free survival and annual measurement of MRD status. MRD was assessed from first-pull bone marrow aspirates at baseline and annually by flow cytometry per International Myeloma Working Group criteria, (limit of detection of at least 1 × 10^-5) up to a maximum of 5 years. Patients who completed at least four cycles of treatment were included in the analysis of the primary endpoint, and patients who had completed at least one dose of treatment on protocol were assessable for secondary endpoints. The study was registered at ClinicalTrials.gov, NCT02538198, and is now closed to accrual.

FINDINGS

Between Sept 8, 2015, and Jan 25, 2019, 108 patients (100 evaluable for the primary endpoint) were enrolled. Median follow-up was 40·7 months (95% CI 38·7-45·0). At 60 months, progression-free survival was 64% (95% CI 52-79). Median progression-free survival was unreached (95% CI unreached-unreached). MRD dynamics were assessed using 340 MRD assessments done over 5 years for 103 evaluable patients. Patients who sustained MRD negativity for 2 years (n=34) had no recorded disease progression at median 19·8 months (95% CI 15·8-22·3) past the 2-year maintenance landmark. By contrast, patients who lost their MRD-negative responses (n=10) were more likely to progress than those with sustained MRD negativity (HR infinite; p<0·0001) and those with persistent MRD positivity (HR 5·88, 95% CI 1·18-33·33; p=0·015) at the 2-year landmark. Haematological and non-haematological serious adverse events occurred in 19 patients (18%). The most common adverse events of grade 3 or worse were decreased lymphocyte count in 48 (44%) patients and decreased neutrophil count in 47 (44%) patients. One death occurred on study due to sepsis and heart failure and was considered unrelated to the study drug.

INTERPRETATION

Serial measurements of MRD allow for dynamic assessment of risk for disease progression. Early intervention should be investigated for patients with loss of MRD negativity. Sustained MRD positivity is not categorically an unfavourable outcome and might portend prolonged stability of low-level disease.

FUNDING

Memorial Sloan Kettering and Celgene.

Safety and Effectiveness of Weekly Carfilzomib, Lenalidomide, Dexamethasone, and Daratumumab Combination Therapy for Patients With Newly Diagnosed Multiple Myeloma: The MANHATTAN Nonrandomized Clinical Trial

Importance

Recently, the benefit of adding daratumumab to the proteasome inhibitor-based, 3-drug combination of bortezomib, lenalidomide, and dexamethasone for patients with newly diagnosed multiple myeloma who underwent high-dose melphalan chemotherapy and autologous hemopoietic cell transplant was assessed. Here, we examine the addition of daratumumab to the second-generation proteasome inhibitor-based, 3-drug combination of carfilzomib, lenalidomide, and dexamethasone.

Objective

To assess the safety and effectiveness of carfilzomib-lenalidomide-dexamethasone-daratumumab combination therapy for patients with newly diagnosed multiple myeloma, in the absence of high-dose melphalan chemotherapy and autologous hemopoietic cell transplant.

Design, Setting, and Participants

Clinical and correlative pilot study at the Memorial Sloan Kettering Cancer Center in New York, New York. Patients with newly diagnosed multiple myeloma were enrolled between October 1, 2018, and November 15, 2019. The median follow-up from start of treatment was 20.3 months (95% CI, 19.2-21.9 months).

Interventions

Eight 28-day cycles with intravenous carfilzomib, 20/56 mg/m2 (days 1, 8, and 15); oral lenalidomide, 25 mg, (days 1-21); dexamethasone, 40 mg weekly, orally or intravenously (cycles 1-4), and 20 mg after cycle 4; and intravenous daratumumab, 16 mg/kg (days 1, 8, 15, and 22 [cycles 1-2]; days 1 and 15 [cycles 3-6]; and day 1 [cycles 7 and 8]).

Main Outcomes and Measures

The primary end point was the minimal residual disease (MRD) rate, in the absence of high-dose melphalan chemotherapy and autologous hemopoietic cell transplant. Secondary end points included determining safety and tolerability, evaluating rates of clinical response per the International Myeloma Working Group, and estimating progression-free survival (PFS) and overall survival (OS) rates.

Results

Forty-one evaluable patients were enrolled (median age, 59 years; range, 30-70 years); 25 (61%) were female, and 20 (49%) had high-risk multiple myeloma. The primary end point (MRD negativity in the bone marrow; 10-5 sensitivity) was achieved in 29 of 41 patients (71%; 95% CI, 54%-83%), and therefore the trial was deemed successful. Median time to MRD negativity was 6 cycles (range, 1-8 cycles). Secondary end points of the overall response rate and the very good partial response or complete response rate were 100% (41 of 41 patients) and 95% (39 of 41 patients), respectively. At 11 months of the median follow-up, the 1-year PFS rate and the OS rate were 98% (95% CI, 93%-100%) and 100%, respectively. Most common (≥2 patients) grade 3 or 4 adverse events were neutropenia (12 patients [27%]), rash (4 patients [9%]), lung infection (3 patients [7%]), and increased alanine aminotransferase level (2 patients [4%]). There were no deaths.

Conclusions and Relevance

In this nonrandomized clinical trial, carfilzomib-lenalidomide-dexamethasone-daratumumab combination therapy was associated with high rates of MRD negativity in patients with newly diagnosed multiple myeloma and high rates of PFS.

DIS3 mutations in multiple myeloma impact the transcriptional signature and clinical outcome

DIS3 gene mutations occur in roughly 10% of patients with multiple myeloma (MM); furthermore, DIS3 expression can be affected by monosomy 13 and del(13q), which occur in approximately 40% of MM cases. Despite several reports on the prevalence of DIS3 mutations, their contribution to the pathobiology of MM remains largely unknown. We took advantage of the large public CoMMpass dataset to investigate the spectrum of DIS3 mutations in MM and its impact on the transcriptome and clinical outcome. We found that the clinical relevance of DIS3 mutations strictly depended on the co-occurrence of del(13q). In particular, bi-allelic DIS3 lesions significantly affected progression-free survival, independently of other predictors of poor clinical outcome, while mono-allelic events mostly affected overall survival. As expected, DIS3 mutations affect the MM transcriptome involving cellular processes and signaling pathways associated with RNA metabolism, and the deregulation of a large number of long non-coding RNA, among which we identified five distinct transcripts as independent predictors of poorer overall survival and nine of worse progression-free survival, with two (AC015982.2 and AL445228.3) predicting both unfavorable outcomes. These findings strongly prompt further studies investigating the relevance of these long non-coding RNA in MM.

Chromothripsis as a pathogenic driver of multiple myeloma

Analysis of the genetic basis for multiple myeloma (MM) has informed many of our current concepts of the biology that underlies disease initiation and progression. Studying these events in further detail is predicted to deliver important insights into its pathogenesis, prognosis and treatment. Information from whole genome sequencing of structural variation is revealing the role of these events as drivers of MM. In particular, we discuss how the insights we have gained from studying chromothripsis suggest that it can be used to provide information on disease initiation and that, as a consequence, it can be used for the clinical classification of myeloma precursor diseases allowing for early intervention and prognostic determination. For newly diagnosed MM, the integration of information on the presence of chromothripsis has the potential to significantly enhance current risk prediction strategies and to better characterize patients with high-risk disease biology. In this article we summarize the genetic basis for MM and the role played by chromothripsis as a critical pathogenic factor active at early disease phases.

Initial Whole-Genome Sequencing of Plasma Cell Neoplasms in First Responders and Recovery Workers Exposed to the World Trade Center Attack of September 11, 2001

PURPOSE

The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole-genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack.

EXPERIMENTAL DESIGN

We performed WGS of nine CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster.

RESULTS

No significant differences were observed in comparing the post-WTC driver and mutational signature landscapes with 110 previously published WGSs from 56 patients with multiple myeloma and the CoMMpass WGS cohort (n = 752). Leveraging constant activity of the single-base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure.

CONCLUSIONS

Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

The molecular make up of smoldering myeloma highlights the evolutionary pathways leading to multiple myeloma

Smoldering myeloma (SMM) is associated with a high-risk of progression to myeloma (MM). We report the results of a study of 82 patients with both targeted sequencing that included a capture of the immunoglobulin and MYC regions. By comparing these results to newly diagnosed myeloma (MM) we show fewer NRAS and FAM46C mutations together with fewer adverse translocations, del(1p), del(14q), del(16q), and del(17p) in SMM consistent with their role as drivers of the transition to MM. KRAS mutations are associated with a shorter time to progression (HR 3.5 (1.5–8.1), p = 0.001). In an analysis of change in clonal structure over time we studied 53 samples from nine patients at multiple time points. Branching evolutionary patterns, novel mutations, biallelic hits in crucial tumour suppressor genes, and segmental copy number changes are key mechanisms underlying the transition to MM, which can precede progression and be used to guide early intervention strategies.

Progression from asymptomatic smoldering multiple myeloma (SMM) to symptomatic Multiple Myeloma occurs at different rates in different patients. Here, the authors report fewer NRAS and FAM46C mutations and adverse translocations in SMM compared to MM, while KRAS mutations are associated with a shorter time to progression.

Designing Evolutionary-based Interception Strategies to Block the Transition from Precursor Phases to Multiple Myeloma

The development of next-generation sequencing technology has dramatically improved our understanding of the genetic landscape of multiple myeloma. Several new drivers and recurrent events have been reported and linked to a potential driver role. This complex landscape is enhanced by intraclonal mutational heterogeneity and variability introduced through the dimensions of time and space. The evolutionary history of multiple myeloma is driven by both the accumulation of different genomic drivers and by the activity of different mutational processes active overtime. In this review, we describe how these new findings and sequencing technologies have been progressively allowed to understand and reshape our knowledge of the complexity of multiple myeloma at each of its developmental stages: premalignant, at diagnosis, and in relapsed/refractory states. We discuss how these evolutionary concepts can be utilized in the clinic to alter evolutionary trajectories providing a framework for therapeutic intervention at early-disease stages.