Ortiz, María

Introduction: Diffuse large B cell lymphoma (DLBCL) patients experience favorable outcomes following frontline immunochemotherapy, however, 30-40% of patients will exhibit primary refractory disease or relapse (rrDLBCL). Recent studies have expanded upon well-defined characteristics of newly diagnosed DLBCL (ndDLBCL) to elucidate genomic landscapes underlying rrDLBCL; yet, high-throughput examination of resistance and relapse mechanisms within rrDLBCL subtypes, including the newly defined dark zone signature (DZsig) addition to cell of origin (COO), have yet to be described. To expand upon the current tumor etiology and define subtype-specific driving mechanisms, we describe a multi-omic analysis leveraging clinical and biologic variables in rrDLBCL. Methods: Tumor biopsies (N=228) collected at the time of relapse were obtained from n=85 participants in the Mayo Clinic/University of Iowa Lymphoma Molecular Epidemiology Resource (MER), n=73 in the Mayo Clinic Biobank, and n=38 from CC-122-DLBCL-001 or n=32 CC-122-ST-001 trials. Tumor samples were analyzed by WES (n=158) or WGS (n=70) and by RNAseq (n=145). Comparative analysis against ndDLBCL MER samples (N=444) utilized WES (n=404) and RNAseq (n=321). Chi Square or Fisher's exact test was used for enrichment analysis of categorical variables, while Wilcoxon or Kruskal-Wallis test was used for continuous variables. Results: rrDLBCL samples were enriched for ABC (46% vs 33% in ndDLBCL, p=0.01), and fewer GCB (27% vs 40%, p=0.006). Amongst rrDLBCL patients, variance in overall survival (from time of original diagnosis; log-rank p=0.001) was observed between COO subtypes, with poor outcomes in patients with DZsig (median 15.0 mo) and GCB (36.5 mo) compared to ABC tumors (74.4 mo). Patients with progression or relapse within 6 months featured a more even COO distribution (ABC - 22%, GCB - 34%, DZsig - 28%) than patients with early (6-24 mo; 51%, 22%, 22%) or late relapses (>24 mo; 62%, 19%, 4%), which were predominantly ABC. Comparing mutation frequency within subtypes, rrDLBCL GCB tumors were enriched for B2M (41% vs 15% in ndDLBCL, p=0.002), EP300 (19% vs 7%, p=0.059), CIITA (19% vs 5%, p=0.02), PRKCB (11% vs 2%, p=0.035), CD79B (11% vs 0%, p=0.004), and CD83 (8% vs 1%, p=0.049) mutations, and CNV events including gains at 8q24.21 (MYC locus; 48% vs 14%, p<0.001) and 18q21.33 (BCL2, TCF4 locus; 36% vs 15%, p=0.021), and losses at 17p13.1 (TP53 locus; 40% vs 19%, p=0.035). rrDLBCL DZsig tumors featured enrichment for BCL6 (22% vs 3%, p=0.025), and SPEN (17% vs 3%, p=0.062) mutations and 2q22.2 losses (LRP1B locus; 20% vs 0%, p=0.019). Despite differential mutation and CNV profiles, gene expression analysis revealed minimal variance within COO subtypes in rrDLBCL and ndDLBCL, however upregulation of lymphoma driving genes including MYBPC2 (LogFC=1.68, FDR<0.001), MYOM2 (LogFC=1.33, FDR=0.007), MYC (LogFC=0.61, FDR=0.04), and CELSR2 (LogFC=0.53, FDR=0.02) were observed in ABC rrDLBCL. Increased frequency of LymphGen classifier A53 was observed in ABC (43% vs 17%, p=0.005) and GCB (40% vs 14%, p=0.009) rrDLBCL. Using CIBERSORTx to infer lymphoma microenvironment (LME) cell type abundance, an LME-depleted class was evident in rrDLBCL tumors, including reduction of CD8+ T cells (effect size (r)=-0.78, p<0.001) and T follicular helper cells (r=-0.48, p=0.012) in GCB rrDLBCL and reduction of T follicular helper cells (r=-0.43, p=0.023) and Tregs (r=-0.32, p=0.065) in DZsig rrDLBCL. Significance: We present an update to the genomic landscape of rrDLBCL, highlighting GCB and DZsig tumors with differential biology and advanced clinical course. Our data supports increased implementation of sequencing efforts in rrDLBCL to more broadly evaluate relapse mechanisms, with optimism that a more detailed understanding of r/r disease may guide trials for targeted treatment strategies. For example, rrDLBCL patients may be strong candidates for MYC pathway targeting agents (BETi), or a Lenalidomide-based regimen (R2 or Tafa-Len), as Lenalidomide has been shown to reprogram the LME via PD-L1 expression in addition to downregulation of MYC and MYC target genes. Or, due to enrichment of variants affecting anti-apoptotic proteins including BCL6 or BCL2 in DZsig and GCB rrDLBCL, a venetoclax-based therapeutic may be considered. Further targeted strategies should be considered as development of a more personalized approach remains an unmet need.

Background Diffuse large B cell lymphoma (DLBCL) is the most common aggressive lymphoid malignancy in adults. Though standard immunochemotherapy regimens can result in clinical remission and cure in a majority of patients, approximately 30% of patients are primary treatment resistant or eventually relapse (rrDLBCL), and their prognosis is very poor. Recent molecular classification by large scale genomic and transcriptomic profiling of newly diagnosed DLBCL tumors has resulted in a deeper understanding of disease drivers, however, the molecular heterogeneity of rrDLBCL is yet to be fully explored. Therefore, in this study, we used an unsupervised approach to better define transcriptional and genetic programs in rrDLBCL samples hypothesizing that molecular stratification within this high-risk group may lead to a better understanding of this disease, which remains an unmet need, and a more personalized therapeutic approach. Methods Data used in this study included RNA (n=143) and whole exome (n=126) sequencing data from available FFPE tumor samples at the time of a relapse (any line of treatment, r1-r10 relapse timepoints included in analysis, one per patient), consented to the Molecular Epidemiology Resource (n=61), banked in the Mayo Lymphoma Biobank (n=50), or consented to the CC-122-ST-001 clinical trial (n=32, NCT01421524). Unsupervised clustering was performed on protein coding gene expression values using the non-negative matrix factorization (NMF) approach. NMF was carried out by performing 200 runs on cluster sizes from 2 to 6. Cell of origin (COO) was determined using the method described by Reddy et al. The tumor microenvironment was analyzed using CIBERSORTx and Lymphoma Microenvironment Classification (LME). Genetic classification was done using LymphGen. Results In order to identify unique transcriptional programs in rrDLBCL, we implemented NMF and consensus clustering to define an ideal quantity of stable patient clusters. We identified 4 rrDLBCL patient clusters (rrC1-rrC4) as the most stable solution, based on consensus clustering selection criteria of iterative NMF runs. While relapse timepoint was not associated with individual clusters, we did observe significant differences in tumor characteristics between the clusters. COO classification found that rrC1 and rrC4 were enriched for GCB-DLBCL while rrC3 was enriched for ABC-DLBCL classification (P < .001). To explore the biologic programs in each cluster, we used pathFindR, which leverages protein-protein interaction networks to identify disease related pathways, rrC1 and rrC4 were defined by pathways involved in cell cycle and mismatch repair; oxidative phosphorylation in rrC2; and MAPK and NF-kB in rrC3. Next, LME classification was done to define patterns in the tumor microenvironment and we found a significant difference between groups (P < .001), with rrC2 and rrC4 tumors enriched for a depleted TME, while rrC1 and rrC3 had a more immune rich or inflammatory TME. This was supported by CIBERSORTx analysis identifying significant differences (P > .001) in abundance of CD8 T cells, memory CD4 T cells, and M1 macrophages between clusters, with rrC4 having the lowest abundance. We then explored genetic patterns in each cluster. Gene level mutation enrichment analysis showed BCL7A mutation exclusively in C4 (6/6, P = .03). The clusters were not associated with LymphGen classification. Lastly, our recently published high risk signature classification that is associated with early clinical failure and is characterized by metabolic dysregulation, a depleted TME, T53 alterations, and poor outcome, was significantly different between the groups (P = .03) with enrichment for high-risk cases in rrC4. In summary, we show for the first time that rrDLBCL patients can be classified into four gene expression clusters that are associated with distinct pathway, TME, and genetic programs. These clusters should now be tested to learn if they can help select patients for newer therapies for rrDLBCL such as CAR-T and bispecific antibodies.

Background: DLBCL has traditionally been classified by cell of origin (COO) subcategories based on tumor gene expression profiles which include Activated B-Cell (ABC) and Germinal Center B-Cell (GCB). Recently, using tumor samples from patients treated with RCHOP, new classification models have focused on DNA alterations. However, a comprehensive integrative approach using a large transcriptomic data set across both newly diagnosed (nd) and relapsed/refractory (r/r) DLBCL is yet to be accomplished. A robust clustering of this type will allow for identification of biologically driven DLBCL patient subgroups and may predict patient outcome and inform treatment approaches. Methods: RNAseq was performed on a total of 882 DLBCL tumor FFPE biopsies from 2 ndDLBCL cohorts (cohort 1 and 2) and 2 r/r DLBCL cohorts (cohorts 3 and 4). Cohort 1 (N=267) was commercially sourced and served as the discovery cohort. Cohort 2 (N=340) was from the Mayo/Iowa Lymphoma SPORE Molecular Epidemiology Resource (MER) and served as the replication cohort. Cohort 3 (N=189) was from the CC-122-ST-001 and CC-122-DLBCL-001 clinical trials (NCT01421524 and NCT02031419), and cohort 4 (N=86) was from r/r patients from the MER. Clustering input consisted of gene expression data, gene set variation analysis (GSVA) scores computed from the hallmark gene sets of MSigDB gene sets, as well as immune cell abundance estimates from a DLBCL-specific deconvolution method. An integrative clustering method iClusterPlus was applied to the input data to identify patient subgroups. A multinomial generalized linear model classifier was trained on the discovery dataset and applied to cohorts 2, 3, and 4 to assess patterns of gene expression and clinical features among the subgroups. Results: Integrative clustering identified 8 subgroups of ndDLBCL patients (DLBCL1-8; D1-D8) in cohort 1. Classifiers trained on cohort 1 were applied to cohort 2 and the same 8 clusters were identified. Among RCHOP treated patients in cohort 2, subgroups D4 (p<0.01) and D8 (p<0.0001) had significantly worse survival outcomes than the rest of the population. D4 comprised 21% of the MER ndDLBCL replication cohort (cohort 2) with a median event-free survival (mEFS) of 38.2 months and a median overall survival (mOS) of 80.3 months. D8 comprised 5% of the cohort with a mEFS of 7.5 months and a mOS of 12.1 months. The remaining 6 subgroups were standard risk, with mEFS ranging from 82.1 months to not reached, and none reaching mOS. The subgroups were not uniquely defined by previously known molecular classification methods such as COO or double hit signature (DHITsig), nor by clinical risk factors such as age or international prognostic index (IPI). Within D4, 92% of patients were ABC, representing a high risk subset of ABC patients. The mEFS in D4 ABCs was 38.2 months, while mEFS of non-D4 ABCs was not reached (p<0.005). Transcriptomic analysis revealed a lower abundance of immune infiltration. D4 was associated with high IPI, with 49% of D4 having IPI>2, compared to 33% of non-D4 with IPI>2 (p<0.05). D8 represented a high-risk subset which was 73% GCB. The mEFS of D8 GCBs was 5.4 months, while mEFS of non-D8 GCBs was not reached (p<0.0001). Transcriptomic analysis revealed low expression of immune response and cytokine signaling pathways, consistent with the low abundance of immune cells in D8. This subgroup consisted of 63% DHITsig positive patients. Although only 20% of all DHITsig positive patients were in D8, these D8 DHITsig patients showed significantly worse survival than non-D8 DHITsig patients (mEFS 11.3 months vs. not reached, p<0.0001). In the r/r DLBCL setting, D1-D8 were all present, with an increased prevalence of D4 and D8 in Cohort 3 (30% and 17%, respectively) and Cohort 4 (30% and 14%) compared to the newly diagnosed setting. Mutational data for these cohorts has been collected and is being interpreted in the context of the discovered subgroups. Conclusion: A novel integrative clustering of transformed gene expression data revealed 8 biologically homogeneous groups, two of which had inferior outcomes when treated with RCHOP therapy. Furthermore, these two subgroups were more prevalent in r/r DLBCL. This classification allows for the transcriptomic identification of high-risk patients underserved by RCHOP therapy. *Ortiz, Wenzl and Stokes contributed equally **Gandhi and Novak contributed equally Disclosures Ortiz: Celgene Corporation: Employment, Equity Ownership. Stokes:Celgene Corporation: Employment, Equity Ownership. Huang:Celgene Corporation: Employment, Equity Ownership. Maurer:Celgene: Research Funding; Morphosys: Membership on an entity's Board of Directors or advisory committees; Nanostring: Research Funding. Towfic:Celgene Corporation: Employment, Equity Ownership. Hagner:Celgene Corporation: Employment, Equity Ownership, Patents & Royalties. Ratushny:Celgene Corporation: Employment, Equity Ownership. McConnell:Celgene Corporation: Employment, Equity Ownership. Danziger:Celgene Corporation: Employment, Equity Ownership. Stong:Celgene Corporation: Employment, Equity Ownership. Lata:Celgene Corporation: Employment, Equity Ownership. Kamalakaran:Celgene Corporation: Employment, Equity Ownership. Mavrommatis:Celgene Corporation: Employment, Equity Ownership. Trotter:Celgene Corporation: Employment, Equity Ownership. Czuczman:Celgene Corporation: Employment, Equity Ownership. Ansell:Seattle Genetics: Research Funding; Affimed: Research Funding; Regeneron: Research Funding; Trillium: Research Funding; Mayo Clinic Rochester: Employment; Seattle Genetics: Research Funding; Mayo Clinic Rochester: Employment; Regeneron: Research Funding; Trillium: Research Funding; Bristol-Myers Squibb: Research Funding; Bristol-Myers Squibb: Research Funding; LAM Therapeutics: Research Funding; LAM Therapeutics: Research Funding; Affimed: Research Funding. Cerhan:NanoString: Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees; Celgene: Research Funding. Nowakowski:Genentech, Inc.: Research Funding; F. Hoffmann-La Roche Ltd: Research Funding; Curis: Research Funding; Bayer: Consultancy, Research Funding; Celgene: Consultancy, Research Funding; Selvita: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding; MorphoSys: Consultancy, Research Funding. Gandhi:Celgene Corporation: Employment, Equity Ownership, Patents & Royalties. Novak:Celgene Coorperation: Research Funding.

Background: Diffuse large B cell lymphoma (DLBCL) is the most common aggressive lymphoid malignancy in adults. Though standard immunochemotherapy regimens can result in clinical remission and cure in a majority of patients, approximately 30% of patients are primary refractory or eventually relapse, and their prognosis is very poor. Recent progress by large scale genomic and transcriptomic profiling of DLBCL patients has resulted in a deep understanding of disease drivers in newly diagnosed DLBCL (ndDLBCL). However, previous genomic studies of relapsed and/or refractory DLBCL (rrDLBCL) are limited by small sample sizes and much less is known about the genomic landscape or the changes in clonal populations that occur within a patient undergoing R-CHOP therapy. Analysis of a large cohort of rrDLBCL is needed to uncover the importance of known drivers, define clonal shifts that occur in relapsing patients, and discover novel mutations that impact resistance to therapy. Methods: DNA sequencing (whole exome or genome sequencing) was performed on tumor and matching germline samples from two rrDLBCL patient cohorts (N=143). The first rrDLBCL cohort consisted of 68 formalin fixed paraffin embedded (FFPE) tumor samples with matched germline obtained from the Mayo/Iowa Lymphoma SPORE Molecular Epidemiology Resource (MER); including 18 diagnostic/relapse paired samples. The second cohort consisted of 75 FFPE samples from 2 rrDLBCL clinical trials CC-122-ST-001 and CC-122-DLBCL-001 (NCT01421524 and NCT02031419). In addition, a ndDLBCL cohort (n=364) consisting of MER front line samples were also analyzed. CNV and SV calling was performed with Sclust and Manta, respectively. Mutational drivers were identified with cDriver. Cosmic mutational signatures were calculated for all samples. Results: For comparison of overall mutation rates between the ndDLBCL and rrDLBCL we analyzed the frequency of mutations in known driver genes [n=211, Chapuy et al, 2018; Schmitz et al, 2018; Reddy et al 2017]. Of these, the top 5 most mutated genes found in the relapse cohort were KMT2D (34.3%), IGLL5 (25.9%), CREBBP (22.4%), TP53 (21.0%), and HIST1H1E (19.6%) all of which were previously reported as mutated genes in rrDLBCL [Morin et al. 2016; Park et al. 2016;Greenawalt et al. 2017]. None of these frequencies changed significantly from the MER ndDLBCL cohort. However, frequencies of 2 genes NCOR1 and ACTB were increased significantly from 3.3% to 13.3% (p= 5.85 x 10-05) and 0.5% to 6.3% (p=2.56 x 10-04), respectively, in the rrDLBCL cohorts. Analysis of the previously reported driver mutations in 18 paired samples revealed a significant increase in the mean variant allele frequency (VAF) from diagnosis to relapse of mutations in BCL2 and BIRC6 from 18.5% to 28.8%(p=0.016) and 4.3% to 13.7%(p=0.011), respectively. In our analysis of novel drivers in rrDLBCL cohorts, 59 significant driver mutations (FDR 0.1) were identified, 23 were novel driver mutations while the remaining 36 have been previously reported as drivers of newly diagnosed DLBCL. These include novel mutations in genes such as IL4R (8.8%), CDK5R1 (6.7%), LAPTM5 (3.7%) and others whose function in lymphoma are unknown. In a secondary analysis, we performed mutational signature analysis for rrDLBCL and identified groups of patients with COSMIC mutational signatures that are associated with aging, defects in homologous end joining resulting in large insertions/deletions, and with defective mismatch repair often found in microsatellite unstable tumors. Together with the high prevalence of TP53 mutations, these results indicate genetic instability as a hallmark of relapse/refractory DLBCL. Analysis of CNV and SV will also be reported. Conclusion: Here we report to our knowledge the largest genomic dataset in relapsed and /or refractory DLBCL using tumor biopsies. Previously identified mutated genes in rrDLBCL were confirmed and novel drivers identified. Paired samples analysis reveals clonal shifts which occur during treatment and at time of relapse. Whether the novel driver mutations appeared as a function of selection of a rare clone over the course of therapy or a therapy induced mutation will be the focus of ongoing work. Novel variants including SNV, CNV and SV can potentially lead to new insights and deepen our understanding of the mechanisms of resistance in DLBCL. * Wenzl, Ortiz and Stong contributed equally **Novak and Gandhi contributed equally Disclosures Ortiz: Celgene Corporation: Employment, Equity Ownership. Stong:Celgene Corporation: Employment, Equity Ownership. Huang:Celgene Corporation: Employment, Equity Ownership. Maurer:Morphosys: Membership on an entity's Board of Directors or advisory committees; Celgene: Research Funding; Nanostring: Research Funding. Towfic:Celgene Corporation: Employment, Equity Ownership. Hagner:Celgene Corporation: Employment, Equity Ownership, Patents & Royalties. Nikbakht:Celgene Corporation: Employment, Equity Ownership. Stokes:Celgene Corporation: Employment, Equity Ownership. Lata:Celgene Corporation: Employment, Equity Ownership. Kamalakaran:Celgene Corporation: Employment, Equity Ownership. Mavrommatis:Celgene Corporation: Employment, Equity Ownership. Trotter:Celgene Corporation: Employment, Equity Ownership. Czuczman:Celgene Corporation: Employment, Equity Ownership. Ansell:Regeneron: Research Funding; LAM Therapeutics: Research Funding; Trillium: Research Funding; Affimed: Research Funding; Bristol-Myers Squibb: Research Funding; Mayo Clinic Rochester: Employment; Affimed: Research Funding; Seattle Genetics: Research Funding; LAM Therapeutics: Research Funding; Mayo Clinic Rochester: Employment; Bristol-Myers Squibb: Research Funding; Trillium: Research Funding; Regeneron: Research Funding; Seattle Genetics: Research Funding. Cerhan:Celgene: Research Funding; NanoString: Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees. Nowakowski:Bayer: Consultancy, Research Funding; Curis: Research Funding; Celgene: Consultancy, Research Funding; MorphoSys: Consultancy, Research Funding; Genentech, Inc.: Research Funding; F. Hoffmann-La Roche Ltd: Research Funding; Selvita: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding. Novak:Celgene Coorperation: Research Funding. Gandhi:Celgene Corporation: Employment, Equity Ownership, Patents & Royalties.

Key Points Using the largest set of patients with newly diagnosed myeloma, we identified 63 mutated driver genes. We identified oncogenic dependencies, particularly relating to primary translocations, indicating a nonrandom accumulation of genetic hits.