Novak, Anne J.

Introduction: Non-Follicular Low-Grade B-Cell Lymphomas (LGBCL) are a subset of lymphomas that share overlapping clinical and histologic features. While most cases of LGBCL have an indolent course, some experience aggressive disease with no definitive biomarkers to identify aggressive cases at diagnosis. Significant advances have been made in the molecular classification of other lymphomas, yet molecular profiling and identification of mechanisms driving aggressive disease in LGBCL remain an unmet need. We previously used a multi-omic approach to characterize the molecular and immune landscapes of LGBCL tumors and identified a transcriptomic signature associated with cell proliferation and inferior overall survival (OS) in two independent cohorts (HR 7.82; 95% CI 2.4-25.4; p < 0.001 & HR 10.07; 95% CI 2.00-50.61; p = 0.005) and cases that transformed to DLBCL (Hopper et al., 2023). Here, we identify the oncoprotein DEK as a novel transcriptional regulator of this signature and use a functional genomic approach to further interrogate the role of DEK in B-cell lymphoma. Methods: Tumors from 64 newly-diagnosed/untreated LGBCL cases enrolled in the University of Iowa/Mayo Clinic Molecular Epidemiology Resource were included in this study (SMZL = 48, NMZL = 6, LPL = 5, B-NOS = 3, EMZL = 2). Tumor RNA-seq and paired tumor-normal WES were used to assess DEK expression and identify copy number alterations (CNAs), respectively. DEK deficient (KO) cell line models were generated using CRIPSR/Cas9. DEK expression was assessed by western blot (WB) and immunofluorescence (IF) using standard techniques. Proliferation was assessed by MTT assay. Results: To identify transcriptional regulators of our previously identified signature, we used the tool MAGIC to mine publicly available ChIP-seq data and identified DEK, an oncoprotein associated with tumorigenesis and poor prognosis in several cancers, as a key regulator of genes in our signature (p = 0.008). RNA-seq indicated that high levels of DEK expression associated with inferior OS in LGBCL (HR 2.99; 95% CI 1.21-7.36; p = 0.02). Additionally, we identified DEK as a target of copy number gains/amps in LGBCL (7% of cases), which corresponded to significantly higher DEK expression (p = 0.04). The association of DEK with aggressive lymphoma was confirmed when we assessed its expression in public cell line RNA-seq data and observed the highest expression in lines derived from aggressive lymphomas, which we further validated at the protein level in a panel of B-cell lymphoma cell lines by WB and IF confocal microscopy. Additionally, we detected nuclear localization of DEK, consistent with its role as a transcriptional regulator. As a model of DEK overexpression for functional characterization, we used the SUDHL6 lymphoma cell line, which has an amplification of DEK at 6p22.3 and expresses high levels of DEK protein. Using CRISPR technology, we generated DEK deficient SUDHL6 and tested the impact of loss of DEK expression. DEK KO cells proliferated significantly slower than WT SUDHL6 cells in vitro (p < 0.0001)(Panel A), indicating DEK plays a role in regulating proliferative capacity in lymphoma. To identify mechanisms responsible for reduced proliferation in DEK KOs, we assessed baseline protein expression levels of apoptosis proteins in WT versus KO cells. We found near complete abrogation of anti-apoptotic protein Bcl-2 and reduction of Bcl-xL expression in DEK KO cells (Panel B). These results indicate DEK may contribute to tumor proliferative capacity by regulating susceptibility to apoptosis. Conclusion: DEK has been implicated in regulation of gene expression, but its precise role has remained elusive, with little known about its role in lymphoma. Using an integrative molecular and functional genomics approach, we identified DEK as a novel oncoprotein in LGBCL that is upregulated in cases with inferior clinical outcome and is regulated by copy number gains/amps. DEK plays a role in regulating tumor proliferative capacity by altering susceptibility to apoptosis through regulation of Bcl-2 and Bcl-xL. Collectively, this data suggests high-risk LGBCL tumors may enrich for cells with DEK overexpression or CNAs as a mechanism driving oncogene addiction, highlighting DEK, and its downstream targets, as potential therapeutic targets in a subset of high-risk LGBCL.

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.

Introduction Lymphoid malignancy (LM) is a heterogenous neoplasm with respect to biology, aggressiveness, and treatment with over 100 subtypes, including Hodgkin lymphoma (HL), non-Hodgkin lymphoma (NHL), and multiple myeloma (MM). Genome-wide association studies have identified common variants associated with risk of specific LM subtypes, but less is known about the contribution of rare inherited pathogenic variants (PV) in the genetic architecture of LM risk. Rare PV have been shown to contribute to the genetic architecture of solid tumor cancer risk, many of which have been found in genes that are important to the DNA-damage repair (DDR) pathway (e.g., BRCA1, BRCA1, ATM). Here, we investigated the prevalence of PV in cancer predisposition genes regularly included on multi-gene panel testing based on NCCN guidelines and correlated with risk of LM overall and by common subtypes. Methods In this association study, cases were newly diagnosed lymphomas (excluding chronic lymphocytic leukemia) from Mayo Clinic who were enrolled between 2002 and 2015 in the Lymphoma SPORE Molecular Epidemiology Resource (MER) and MM from Mayo Clinic enrolled primarily between 1998 and 2022 in a clinical registry. Controls were from the Mayo Clinic Biobank who were enrolled between 2009 and 2016, from which we excluded individuals with a history of hematologic malignancy. All participants provided a blood sample for whole exome sequencing, performed by Regeneron on an Illumina NovaSeq panel (mean coverage of 48X). Variants were called using GATK v4 and variant annotation was performed using BioR. All loss-of-function variants (i.e., nonsense, frameshift, and consensus splice sites) and pathogenic missense in 19 cancer predisposition genes with a minor allele frequency < 0.5% were included in the analyses. Logistic regression was used to estimate odds ratio (OR) and 95% confidence intervals (CI) for the association between mutation status by gene (burden test) and risk of LM overall, and by subtype that had more than 5 mutation carriers. All analyses were adjusted for age (at diagnosis for cases and at enrollment for controls) and sex. Results A total of 6990 cases and 42432 controls were included. Median age was 63 years for both cases and controls (case range 18-94 years, control range 18-99 years). Males accounted for 58.5% (n=4090) of the cases and 41.7% (n=17762) of the controls. MM accounted for 31% (n=2138) of the cases, followed by diffuse large B-cell lymphoma (DLBCL, n=1146) and follicular lymphoma (FL, n=1132) each accounting for 16% of cases. CHEK2 (1.0%), ATM (0.4%), BRCA2 (0.4%), and BRCA1 (0.3%) had the most PV carriers in controls; these frequencies are in the range of published control populations. A total of 1816 (3.7%) individuals had a PV across the 19 genes, with 4.7% of cases and 3.5% of controls carrying any PV (OR=1.34, 95% CI: 1.18-1.51). ATM (OR=1.86, 95% CI: 1.36-2.49), CHEK2 (OR=1.74, 95% CI: 1.42-2.13) and TP53 (OR=9.07, 95% CI: 4.51-18.87) were all significantly associated with increased risk of LM overall (Figure 1). When investigating subtypes (Table 1), CHEK2 (OR=1.93, 95% CI: 1.23-2.89) and TP53(OR=10.97, 95% CI: 3.06-31.43) were significantly associated with DLBCL. ATM (OR=2.22, 95% CI: 1.13-3.90) was significantly associated with FL. ATM (OR=1.97, 95% CI: 1.16-3.12), CHEK2 (OR=1.75, 95% CI: 1.24-2.41) and TP53 (OR=7.67, 95% CI: 2.43-20.72) were significantly associated with MM and ATM (OR=8.62, 95% CI: 4.47-15.12), NBN (OR=6.94, 95% CI: 2.41-15.76) and TP53 (OR=38.30, 95% CI: 10.27-116.69) were significantly associated with mantle cell lymphoma (MCL). When subset to LM cases overall, there was no evidence of enrichment for young onset (age <60 years), male sex, or family history of hematologic malignancy between PV carriers and non-PV carriers. However, FL cases with an ATM PV tended to be less than 60 years (OR=3.21, P=0.09) compared to non- ATM PV carriers. Similarly, MM cases with a TP53 PV tended to be less than 60 years of age (OR=7.39, P=0.07) compared non-carriers. Discussion In this large case-control study, PV in established cancer predisposition genes were associated with increased risk of LM overall, demonstrating that rare inherited variants may play an important etiologic role. This study also suggests that PV in cancer predisposition genes may be subtype specific. These results will help inform LM risk for individuals who undergo genetic testing of cancer predisposition genes.

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 Lymphoma is the sixth most commonly diagnosed cancer in the US. Genome-wide association studies have identified common variants associated with risk of specific lymphoma subtypes, but less is known about the contribution of rare inherited variants in the genetic architecture of lymphoma risk. CHEK2 is important to DNA repair and 2 small studies have found evidence of an association between CHEK2 variants and risk of lymphoma. Here, we investigated loss of function (LoF) variants in CHEK2 with risk of lymphoma (overall and subtypes). The study population included newly diagnosed lymphoma cases from Mayo who were enrolled in the Lymphoma SPORE Molecular Epidemiology Resource (MER). Controls were from the Mayo Clinic Biobank, from which we excluded a prior hematologic malignancy. Whole exome sequencing was performed by Regeneron on an Illumina NovaSeq panel (mean coverage of 48X). Variants were called using GATK v4 and variant annotation was performed using BioR. All LoF variants (nonsense, frameshift and consensus splice sites) in CHEK2 with a minor allele frequency < 0.5% were included in the analyses. Logistic regression was used to estimate odds ratio (OR) and 95% confidence intervals (CI) for the association between CHEK2mutation status and risk of lymphoma overall, and by lymphoma subtypes that had more than 5 mutation carriers, which included diffuse large B-cell (DLBCL), follicular (FL), and T-cell lymphoma (TCL). All analyses were adjusted for age (at diagnosis for cases and at enrollment for controls) and sex. A total of 4,852 lymphoma cases and 49,724 controls were included in this analysis. Median age for both cases and controls was 62 years (range 18-99 years). Males accounted for 57.5% (n=2,789) of the cases and 40.6% (n=20,186) of the controls. The DLBCL (23.6% of cases) and FL (23.3% of cases) subtypes were the most common; 7.6% of cases had TCL. A total of 407 (0.7%) individuals had a LoF variant. The frequency of LoF variants was 1.2% in lymphoma cases and 0.7% in controls. LoF variants were associated with increased risk of lymphoma overall (OR=1.77; 95%CI: 1.32-2.33), DLBCL (OR=2.07; 95%CI: 1.20-3.13), and TCL (OR=2.78; 95%CI: 1.178-5.50), but not FL (OR=1.31; 95%CI: 0.65-2.34). CHEK2 c.1100delC was the most frequently mutated variant in this population, accounting for 76% of all LoF variants. Individuals with a c.1100delC variant had an increased risk of lymphoma overall (OR=1.90; 95%CI: 1.37-2.58), DLBCL (OR=2.06; 95%CI: 1.09-3.54), and TCL (OR=3.17; 95%CI:1.24-6.58), but not FL (OR=1.05; 95%CI: 0.41-2.15). When restricted to cases, there was no significant difference in age at diagnosis (P=0.62) or sex (P=0.53) between LoF variant carriers and non-carriers. In this large lymphoma case-control study, CHEK2 LoF variants were associated with increased risk of lymphoma, demonstrating that rare inherited variants may play an important role in the etiology of lymphoma. Additional work is needed to investigate missense variants in CHEK2 and risk of lymphoma. Citation Format: Nicholas J. Boddicker, Raphael Mwangi, Dennis P. Robinson, Allison C. Rosenthal, Thomas M. Habermann, Andrew L. Feldman, Lisa M. Rimsza, Rebecca L. King, Melissa C. Larson, Brianna J. Gysbers, Stephen M. Ansell, Jithma P. Abeykoon, Grzegorz S. Nowakowski, Thomas E. Witzig, Anne J. Novak, Susan L. Slager, James R. Cerhan. Germline CHEK2 variants and risk of lymphoma. [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 5233.

Abstract Introduction: Cytokines are protein factors that regulate diverse cellular processes. The malignant progression of immunoglobulin M monoclonal gammopathy of undetermined significance (IgM-MGUS) to Waldenstrom macroglobulinemia (WM) is in part cytokine mediated, and we have previously shown IL-6 mediates IgM secretion in WM via the JAK/STAT pathway. Cytokine signaling can be epigenetically modulated by micro-RNAs (miRs). We aimed to identify the role of miRs in the regulation of inflammatory cytokine pathways between these diseases. Methods: Bone marrow samples of patients with IgM-MGUS (n=7) and WM (n=25) were prospectively collected and sorted for CD19+ and/or CD138+ malignant cells. Total RNA extraction and sequencing were performed, miR and mRNA profiling was conducted, and differential expression was calculated between WM and IgM-MGUS using a log2 fold change (FC)>0.5 and ←0.5, false discovery rate<0.05. Experimental differentially expressed miR-mRNA pairing was performed, and pairs were filtered for correlated expression (i.e., upregulated miR/downregulated mRNA, and vice-versa). Filtered pairs were analyzed using the Ingenuity Pathway Analysis to identify pathways between WM and IgM-MGUS with absolute activation Z-score ≥1. Results: We identified 34 differentially expressed miRs experimentally targeting 1132 mRNAs. Pathway analysis showed overall miR-based inactivation of inflammatory cytokine signaling pathways in WM compared to IgM-MGUS, especially of the interleukin (IL) family: IL-2, IL-6, IL-8, IL-9, IL-15, and IL-17. Underlying multiple pathways was upregulation of miRs-146a, 150, 194, all targets of Signal Transducer and Activator of Transcription 1 (STAT1), observed to be downregulated. Of relevance, miR-146a has been previously shown to be an inhibitor of inflammatory cytokines, including IL-6. Next, we observed upregulation of miR-142 which targets phosphoinositide-3-kinase regulatory subunit 6 (PIK3R6), observed to be downregulated. Both the STAT and PIK3 families are important cytokine signal transducers, and we have previously shown the PI3k-AKT-IkB-p65 pathway induces expression of IL-6 via a GLI2-mediated mechanism. We additionally observed granulocyte-macrophage colony-stimulating factor (GM-CSF) signaling to be inactivated in WM. Here, miR-19b and miR-30e were upregulated, both target CSF 2 receptor subunit beta (CSF2RB), found to be downregulated. CSF2RB is the common receptor subunit of GM-CSF, IL-3, and IL-5. Conclusion: Our study identifies multiple cytokines pathways involved in inflammation and the tumor microenvironment to be potentially inactivated by differentially expressed miRs in WM compared to IgM-MGUS. Downregulation of these pathways and targeting of the STAT and PIK3 family indicates a possible miR-based negative feedback loop to down-modulate inflammatory cytokine signaling in WM. Citation Format: Karan Chohan, Jonas Paludo, Surendra Dasari, Jithma P. Abeykoon, Prashant Kapoor, Esteban Braggio, Michelle K. Manske, Aneel Paulus, Craig B. Reeder, Sikander Ailawadhi, Asher A. Chanan-Khan, Robert A. Kyle, Morie A. Gertz, Anne J. Novak, Stephen M. Ansell. Micro-RNA induced inactivation of inflammatory cytokine signaling pathways in Waldenstrom macroglobulinemia compared to IgM-MGUS. [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 3811.

Abstract Non‐follicular low‐grade B‐cell lymphomas (LGBCL) are biologically diverse entities that share clinical and histologic features that make definitive pathologic categorization challenging. While most patients with LGBCL have an indolent course, some experience aggressive disease, highlighting additional heterogeneity across these subtypes. To investigate the potential for shared biology across subtypes, we performed RNA sequencing and applied machine learning approaches that identified five clusters of patients that grouped independently of subtype. One cluster was characterized by inferior outcome, upregulation of cell cycle genes, and increased tumor immune cell content. Integration of whole exome sequencing identified novel LGBCL mutations and enrichment of TNFAIP3 and BCL2 alterations in the poor survival cluster. Building on this, we further refined a transcriptomic signature associated with early clinical failure in two independent cohorts. Taken together, this study identifies unique clusters of LGBCL defined by novel gene expression signatures and immune profiles associated with outcome across diagnostic subtypes.