Huang, C. Chris

Introduction: Clonal hematopoiesis (CH) is a common age-related premalignant condition defined by the presence of an acquired somatic mutation in a gene frequently mutated in myeloid malignancies characterized by a variant allele frequency (VAF) of ≥2%, without fulfillment of other diagnostic criteria for myeloid neoplasm. Incidence of CH increases with age and found in 10%-30% of people > 70 years of age. CH-related mutations confer an increased risk of myeloid neoplasms, cardiovascular diseases, and are linked to inferior survival in solid tumor patients treated with chemo/radiation therapy. Such mutations can act as a driver/amplifier of inflammation, especially through myeloid cells. To distinguish true disease drivers from background genetic events related to aging or environmental factors, we need a deeper understanding of the impact of CH in different tissue microenvironments and under therapeutic pressure for better risk stratification of CH carrier patients. Recent reports suggested high prevalence of CH-related mutations in lymphoma. Current analyses aimed to investigate: i) the prevalence of CH-related mutations in two real-world lymphoma cohorts from Molecular Epidemiology Resource (MER) MER-FL and MER-DLBCL, ii) CH-related mutational landscape in DLBCL and FL, iii) the impact of CH-related mutations on clinical outcomes. Methods: Whole exome sequencing (WES) derived from blood/PBMC from lymphoma patients with paired tumor samples was analyzed to identify CH-related myeloid drivers in 42 previously reported genes. We used Sention-Tnhaplotype2 somatic pipeline with a custom panel of normals from whole blood samples. To avoid tumor contamination, we filtered variants for which VAF in the tumor was higher than that in the germline sample. Other filtering criteria included HIGH/MODERATE impact; read count > 1 in both forward and reverse reads; gnomAD_AF < 0.01. We used the R “survival” package Cox regression implementation for outcome and other statistical analysis. Results: Genomic and clinical data of 507 patients from MER-FL and 383 patients from MER-DLBCL cohorts, with median age of 60 and 64 respectively, were available for analysis. CH-related mutations were present in 17.16% (87/507) in FL and 18.28% (70/383) in DLBCL cohorts with a VAF cut off >2%. The most common CH-related mutations were DNMT3A (5%), CREBBP (2%), TET2 (1%), ASXL1 (1%) and PPM1D (1%) in MER-FL while DNMT3A (5%), TET2 (3%), TP53 (3%), ASXL1 (2%) and PPM1D (1%) in MER-DLBCL cohort. Both the prevalence of CH-related mutations and their average VAF increased with age, consistent with literature. Finally, cox regression analysis suggested that CH-related mutations were associated with significantly inferior overall survival (OS) in both FL and DLBCL cohorts (HR 1.6 p=0.028 and HR 2.3 p=<0.001 respectively); however only in DLBCL this association was significant after age correction (HR 1.8 p=0.008). When separated by individual gene mutations, genes of interest with inferior OS included TP53 (n=4, HR 3.1, p=0.057) in FL and PPM1D (n=3 HR 2.6 p=0.1), TP53 (n=8 HR 2.1 p=0.1) and DNMT3A (n=14, HR=2.0, p=0.048) in DLBCL; however, these estimates did not quite reach statistical significance. Conclusions: Taken together, our exploratory analysis indicates that CH-related mutations associate with worse OS in lymphoma patients, supporting the rationale for evaluation of CH for the risk assessment in lymphoma. Furthermore, DLBCL displays of larger size effect, which may be related to different treatment regimes. Additional studies will be needed to elucidate the interaction between high dose chemotherapy regimes and CH-associated adverse events.

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.

Abstract Introduction: Deletion in the short arm of chromosome 17 (del17p) is a strong, negative prognostic biomarker found in about eight percent of patients with untreated chronic lymphocytic leukemia (CLL). CLL patients with del17p exhibit resistance to therapy, accelerated disease progression, and global immunodeficiency. It is unclear the extent to which global immunodeficiency impacts disease course and therapy outcomes. The goal of this study was to document the blood immune repertoire of patients with untreated CLL who have del17p and associate with clinical features, time to first treatment (TFT), and overall survival (OS). Methods: Biobanked blood samples from treatment-naïve del17p CLL patients (n = 58) and age-matched controls (n = 5) were obtained. B cells and T cells were analyzed with comprehensive, cell type-specific flow cytometry panels and visualized with unbiased dimensionality reduction. Lab values were log-transformed and compared using t-tests and one-way ANOVAs. TFT was calculated from date of sample collection to date of first CLL-directed therapy or last known untreated date; OS was calculated from date of sample collection to death or last known alive status. Cox proportional hazards regression analyses (results reported as hazard ratio [HR] and 95% confidence interval [CI]) were employed to identify associations of lab values with OS and TFT. Results: Patients with del17p CLL had increased proportions of T regulatory cells (Tregs) among total CD4+ T cells compared to controls (CLL 11.1%, control 5.7%, p = 0.05). A higher proportion of CD8+ T cells was associated with longer TFT (HR 0.57, 95%CI 0.36-0.89); in contrast, higher proportions of TIGIT+ (HR 1.83, 95%CI 1.03-3.28) and TIM3+ conventional (non-Treg) CD4+ T cells were associated with shorter TFT (HR 1.70, 95%CI 1.02-2.84). A higher proportion of PD1+ CD8+ T cells was associated with shorter OS (HR 1.79, 95%CI 1.05-3.06), in addition to PD1+ conventional CD4+ T cell proportion (HR 2.21, 95%CI 1.06-4.60). CLL B cells expressed PD-1 and inducible T cell costimulator (ICOS) at significantly lower median fluorescence intensities (MFIs) than control B cells (p < 0.0001 for both). The CXCR5+ MFI (HR 1.0002, 95%CI 1.0001-1.003) and CD23+ MFI (HR 1.88, 95%CI 1.01-3.50) of CLL B cells were associated with shorter TFT. Discussion: In this cohort of treatment-naïve del17p CLL patients, we determined that patients with higher proportions of exhausted CD4+ and CD8+ T cell phenotypes had significantly shorter OS. Cytotoxic T cell proportion positively associated with TFT, whereas proportions of CD4+ T cell exhaustion phenotypes negatively associated with TFT. Higher CXCR5+ and CD23+ MFIs of CLL B cells also associated with shorter TFT. The preliminary findings from our study suggest that the blood immune repertoire has strong associations with clinically relevant outcomes such as time to first treatment and overall survival in patients with untreated del17p CLL. Citation Format: Ethan J. Katcher, Kim A. Gwin, Sameer A. Parikh, Kari G. Rabe, Petra K. Lothert, Preeti Trisal, Esteban Braggio, Susan L. Slager, Nicholas Strong, Daniel L. Van Dyke, Zhiquan Wang, C. Chris Huang, Anita K. Gandhi, Kay L. Medina, Neil E. Kay. Correlating the blood immune repertoire and clinical features of del17p CLL patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB345.

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.