Stong, Nicholas

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.

ABSTRACTDiffuse large B‐cell lymphoma (DLBCL) patients that fail to achieve a complete metabolic response with frontline immunochemotherapy have a poor prognosis. Genomic profiling has led to a broader understanding of the molecular drivers in DLBCL, but it is unknown how well current classifiers identify patients that will experience primary treatment resistance (PTR). Using whole exome and RNA sequencing data from newly diagnosed DLBCL patients, we evaluated the genomic landscape of PTR and compared it to that of non‐PTR DLBCL. We found a significant increase in the frequency of TP53 (34% vs. 15%, p = 0.005) and ARID1A mutations (21% vs. 7%, p = 0.007) in PTR cases, with pathway analysis further demonstrating a downregulation of TP53 and an increase in chromatin modifying pathways. These results suggest that TP53 and ARID1A may be key mediators of PTR and important pathways contributing to the poor outcomes. We found that the current molecular classifiers were unable to identify PTR cases at diagnosis. However, our newly identified high‐risk signature identified 46% of PTR cases at diagnosis. Overall, these results contribute to our understanding of the genomic landscape of patients with primary treatment resistance.

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.

Presence of deletion17p (del17p) in frontline (1L) and relapsed/refractory (R/R) CLL has been associated with shorter time to treatment initiation, early relapse after therapy and inferior overall survival (OS) (Mato et.al., Haemotologica 2022). However, there exists notable clinical heterogeneity in terms of outcome within this group (Tam et.al., Blood 2009, Stephens et.al; AJH 2015). Previous reports in 1L and R/R CLL have identified percentage of del17p nuclei by fluorescence in situ hybridization (FISH), having TP53 mutations on the second allele, copy number alterations and complex karyotype (CK) as additional risk factors within the del17p population (Yu et al., CCR 2017; Baliakas et al., Blood 2019). Here we profile a cohort of wait and watch (W&W, n=14), 1L treated (n=36) and R/R CLL (n=3) patients by whole genome sequencing (WGS) and RNA sequencing (RNAseq) to gain insights into genomic factors and gene expression pathways contributing to the high-risk nature of del17p across the spectrum of CLL patients. Using the Mayo Clinic CLL Database, we identified 53 patients with del17p (identified by FISH testing in routine clinical practice at diagnosis) who had stored peripheral blood mononuclear cells (PBMCs) available. RNA and DNA from matched peripheral blood (tumor) and DNA from CD19- (germline) samples at diagnosis (n=50) and relapse (n=3) was extracted using Qiagen RNeasy (RNA) and Purene (DNA) kits and subjected to RNAseq and WGS respectively. Additional FISH assays (del11q, trisomy 12, del13q) and CK by chromosomal banding analysis (CBA) were performed. Baseline clinical characteristics, time to first treatment, types of treatment, and survival were abstracted from the CLL Database. The Mayo Clinic IRB approved this study. Copy number analysis identified deletions in 3p, 4p, 8p, 9p, 9q and gains in 2p, 8q, 15q and 17q chromosomal lesions occurring in patients. We explored the association of clinical and genetic factors, and genomic instability signatures with survival to identify high-risk features in del17p patients. Both total number of deletions by WGS and number of chromosomal aberrations by CBA as a continuous variable were significantly associated with unfavorable OS in the entire cohort (n=53) in a univariate analysis (Hazard ratio, HR= 2.29, p= 0.00054 and HR= 2.33, p=0.0054 respectively). However, there was no significant association of number of gained genes by WGS or aneuploidy scores with survival (HR=0.74, p=0.37). In patients with CK data (n=41), deletion burden, but not gains, were highly correlated to chromosomal aberrations by CBA. Using cutpointr package for identifying deletion threshold corresponding to CK≥ 3, we defined patients with deletion burden ≥ 514 as enhanced-high risk (EHR; median OS = 4.83 years) compared to the group with high-risk (HR; median OS=9.94 years) (p=0.025). There was no significant difference in number of deleted genes between WW and treated patients. The EHR group is enriched in patients with high clonal cell fraction of del17p by FISH, unmutated IGHV and CK≥ 3. The EHR group remains significant in a multi-variable analysis with IGHV-mutation status, ZAP70, and age (HR: 5.16, 95% CI: 1.04 - 25.6, p=0.04). Mutational analysis identified 35 known CLL driver mutations with TP53 (86%), ATM (15%), NOTCH1(13%), MGA (13%), CHD2 and MED12(11%) as the top mutated genes. Total number of driver genes mutated was not significantly different between EHR and HR subgroups, but the EHR subset was significantly enriched with mutations in NOTCH1 (21%), and MGA (21%) as compared to no mutations in the HR group in these genes. Deleted genes (≥ 514) was prognostic in bi-allelic TP53 aberration patients as well (p=0.021) Differential expression gene-set enrichment analysis indicated a significant enrichment of DNA repair, cellular response to DNA damage and downregulation of inflammatory responses in the EHR group. CK is a known prognostic marker which includes both chromosomal losses and gains. In this small cohort of del17p we were able to show the prognostic value of deletions only, but not gains, as captured by WGS. Using genome wide sequencing, we identify increasing genomic deletions as a feature of enhanced-high risk del17p. While deletion burden cut-offs identified here are specific to our research method and require further validation in additional independent cohorts, EHR subgroup remained significant after adjusting for other known prognostic variables .

Background: Tumor mutational burden (TMB), defined as the number of nonsynonymous mutations per megabase (Mb) in the coding regions of a tumor genome, has been associated with improved response to immune checkpoint blockade treatments and more broadly with overall survival (OS) in many cancers. Limited prior studies have found that TMB is higher in diffuse large B-cell lymphoma (DLBCL) compared to many other cancers. However, only two studies have evaluated TMB and prognosis in newly diagnosed DLBCL patients. The prognostic role of TMB in diffuse large B-cell lymphoma (DLBCL) remains unclear as prior studies have been limited to panel sequencing of cancer genes and/or lacked matched germline-tumor samples for filtering. Methods: We characterized TMB from whole exome sequencing (WES) in newly diagnosed DLBCL patients with matched diagnostic tumor-normal samples prior to treatment with standard immunochemotherapy (IC). Tumor DNA was extracted from formalin-fixed, paraffin embedded tissue sections and germline DNA from peripheral blood. We evaluated the association between TMB and event-free (EFS) and OS using Cox regression adjusted for the international prognostic index (IPI) and other factors including cell of origin (COO), MYC double hit (DH) with BLC2 or BCL6 by FISH, and TP53 mutation status. EFS and OS were defined as time from diagnosis to earliest of treatment failure (EFS only), death, or last follow-up; EFS at 24 months (EFS24) was also modeled as a marker of early aggressive DLBCL. An optimal high/low TMB cutpoint was determined based on EFS using maximally selected rank statistics; sensitivity analyses used an 20/80 percentile to define the threshold. We also evaluated the association of driver tumor mutational load (driverTML) with EFS and OS, defined as the total number of mutations in 292 DLBCL driver genes. Results: In the analysis cohort of 394 DLBCL patients, the median age at diagnosis was 65 years (interquartile range: 56-73), 43% were female, 62.2% had IPI 0-2, and 57.4% germinal center B-cell COO. MYC was DH in 7.5% and TP53 was mutated in 16.2% of patients. Through a median follow-up of 7.5 years, there were 141 events and 122 deaths. EFS24 was achieved in 75.1% of the cohort; 5-year EFS was 69% (95% CI 64-74%) and 5-year OS was 75% (95% CI 71-80%). Median TMB in our cohort of DLBCL patients was 3.6 mutations/Mb, which was the fifth highest compared with 32 cancers in The Cancer Genome Atlas. As a continuous variable, increasing TMB (mutation/Mb) was associated with superior EFS (HR=0.89, 95% CI 0.81-0.98) and OS (HR=0.89, 95% CI 0.81-0.99) after adjustment for IPI. Spline analyses confirmed linear associations of TMB with EFS and OS. The optimal modeled cutpoint based on EFS was 4.5 mutations/Mb, with 31% classified as high TMB (>4.5). No significant ( P >.05) differences were observed between high and low TMB patients with clinical (i.e., age, sex, IPI category) or molecular (cell of origin, MYC DH, TP53 mutated) factors. In Cox regression, high-TMB was significantly protective for EFS and for OS compared to low-TMB after adjustment for each clinical factor ( Figure 1), with simultaneous adjustment for all these factors associated with superior EFS (HR=0.52, 95% CI 0.33-0.84) and OS (HR=0.59, 95% CI 0.36-1.00). TMB (high/low) combined with IPI (0-2, 3-5) was significantly ( P =.001) associated with EFS in Kaplan-Meier analysis ( Figure 2). Among patients with IPI 3-5, high TMB was significantly associated with superior EFS (HR=0.47, 95% CI 0.27-0.81) and OS (HR=0.57, 95% CI 0.33-0.99). Among patients with IPI 0-2, we also observed longer EFS (HR=0.62, 95% CI 0.34-1.11) and OS (HR=0.62, 95% CI 0.32-1.20) in high TMB patients, although the associations were not statistically significant ( P >.05). In sensitivity analyses, a 20/80 percentile TMB threshold found similar results for EFS (HR=0.64, 95% CI 0.41-1.03) and OS (HR=0.69, 95% CI 0.42-1.12) as in our main analyses. TMB and driverTML were significantly correlated ( R=0.65, P <10 -16) and we found similar effect sizes for EFS (HR=0.45, 95% CI 0.22-0.91) and OS (HR=0.55, 95% 0.27-1.13) as in our main analyses. Conclusions: In a large and well-defined cohort of IC treated DLBCL patients, we found DLBCL patients with higher TMB had superior EFS and OS after adjustment for other clinical and molecular prognostic factors, and the prognostic association for TMB was robust to threshold and gene set.

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

Sonic hedgehog signaling regulates processes of embryonic development across multiple tissues, yet factors regulating context-specific Shh signaling remain poorly understood. Exome sequencing of families with polymicrogyria (disordered cortical folding) revealed multiple individuals with biallelic deleterious variants in TMEM161B , which encodes a multi-pass transmembrane protein of unknown function. Tmem161b null mice demonstrated holoprosencephaly, craniofacial midline defects, eye defects, and spinal cord patterning changes consistent with impaired Shh signaling, but were without limb defects, suggesting a CNS-specific role of Tmem161b. Tmem161b depletion impaired the response to Smoothened activation in vitro and disrupted cortical histogenesis in vivo in both mouse and ferret models, including leading to abnormal gyration in the ferret model. Tmem161b localizes non-exclusively to the primary cilium, and scanning electron microscopy revealed shortened, dysmorphic, and ballooned ventricular zone cilia in the Tmem161b null mouse, suggesting that the Shh-related phenotypes may reflect ciliary dysfunction. Our data identify TMEM161B as a regulator of cerebral cortical gyration, as involved in primary ciliary structure, as a regulator of Shh signaling, and further implicate Shh signaling in human gyral development.