Braggio, Esteban

Background: Patients with relapsed/refractory multiple myeloma demonstrate an increasing incidence of extramedullary disease (EMD), noted in up to 30% of patients. Presence of EMD portends a markedly inferior prognosis with limited response to salvage therapies, including immune effector therapy. Understanding the molecular alterations in EMD is key to exploring novel therapeutic approaches but remains to be well-characterized. Here we describe the mutational landscape and biologic pathways with alterations in a large cohort of patients with EMD. Methods: We included 18 paraffin-embedded tumor tissues from true EMD biopsies (non-bone contiguous disease) for whole exome sequencing (WES); biopsies from paraskeletal/bone-contiguous soft tissue were not considered. Germline sequencing was performed utilizing cryopreserved peripheral blood mononuclear cells (PBMCs) at initial diagnosis from patients without evidence of circulating tumor cells. Additionally, we studied 12 CD138-sorted cryopreserved bone marrow (BM) aspirates collected at the time of development of EMD. Following standard DNA extraction and library construction practices, we performed WES on the EMD, BM and PBMC samples. Variant calling was performed using consensus calls observed in at least 2 out of the 3 callers (Mutect2, VarScan, Strelka2) and these were included for downstream analysis. Synonymous and low-impact variants were filtered. Pathway analysis was performed using the PanCancer driver gene sets and compared using Fisher's exact test. Results: The median age of patients with EMD was 65.6 years (range 46.2-86 years). Seven (39%) patients had de novo (previously untreated) EMD and 11 (61%) had secondary EMD (developed at relapse). The MAPK pathway was frequently mutated, most commonly involving NRAS (33%), KRAS (22%) and BRAF (22%) genes. These mutations were mutually exclusive. Alterations in TP53 were noted in 11% of EMD samples. Single nucleotide variants (SNVs) in epigenetic modulators were frequently noted in EMD, with the most prevalent alterations seen in ARID1A (33%), KMT2C (28%), EP300 (22%) and KMT2B (22%). Both ARID1A and EP300 are postulated to play a role in IRF4 regulation and plasma cell differentiation. Additionally, inhibition of EP300 can be a viable therapeutic strategy in MM through disruption of IRF4 signaling. We noted frequent mutations in the roundabout (ROBO) genes ROBO1 and ROBO2, observed in 22% and 28% of EMD samples, respectively. ROBO1 has been demonstrated to be necessary for myeloma cell adhesion to BM stromal and endothelial cells, with ROBO1 knockout MM cells infused in mice demonstrating a tendency for extramedullary dissemination. Similarly, variants in LRP1B, a frequently altered large putative tumor suppressor gene in the LDL receptor family, were noted in 33 % of EMD samples. Depletion of LRP1B has been demonstrated to promote anchorage-independent growth, cell migration and invasion in vitro. Notably, in the paired BM samples, no SNVs were noted in ARID1A, LRP1B, ROBO1 and ROBO2 genes (q-value not significant). In paired samples, among the shared mutations of interest, we noted expansion of the BRAF clone in 2 EMD samples [variant allele frequency (VAF) of 3% and 13% in BM and 58% and 29% in the paired samples, respectively]. An expansion of a shared KRAS mutation was noted in another EMD sample (VAF of 15% in BM and 29% in EMD). The frequently altered cellular pathways in EMD samples (n=18) included the MAPK pathway (89%), chromatin histone modifier pathway (83%), transcriptional regulation (83%), immune signaling (72%) and RNA abundance (72%). In the cohort of paired EMD and BM samples (n=24), notable enrichment was evident in EMD for MAPK pathway alterations [91% in EMD vs. 25% in BM q=0.02] and chromatin histone modifier pathway (83% in EMD vs. 42% in BM, q <0.01). Conclusions: Patients withEMD demonstrate a near universal presence of alterations in the MAPK pathway. Variants in chromatin modifiers regulating plasma cell differentiation and genes involved in BM stromal adhesion and homing of plasma cells were frequently noted, potentially facilitating the extramedullary spread of myeloma cells. Functional studies are needed to evaluate these hypotheses. Targeting EP300, ARID1A and downstream MEK/ERK inhibition as part of combination therapies may represent promising therapeutic approaches for further exploration in EMD.

Introduction Monoclonal antibodies (MoAbs) directed against CD38 are a mainstay in multiple myeloma (MM). Downregulation of CD38 expression is a resistance mechanism of antigenic escape, however, cell surface expression may return. While this tumor cell plasticity is likely driven by epigenetic mechanisms, it is unknown if MM cells can develop permanent resistance to anti-CD38 MoAb by acquiring genomic events leading to biallelic disruption of CD38 thus precluding re-treatment. Here, we report biallelic loss of CD38 as a recurrent mechanism of resistance, often via convergent evolution where distinct subclones are advantaged towards the same selective pressure. Methods We interrogated two sets of patients (Mayo Clinic; n = 31, Calgary University; n =19) with relapse following anti-CD38 MoAbs. Two additional cases of interest were selected from Heidelberg University. Whole genome sequencing (WGS, 60-100x) or whole exome sequencing (WXS), flow cytometry (flow), and bulk RNAseq were used to validate the effect of genomic events on CD38. Missense mutations predicted to affect CD38 MoAb binding were validated with site-directed mutagenesis of a CD38 plasmid, transduction into K562 MM cells, and CD38 MoAb binding assays. Results The prevalence of biallelic inactivation of CD38 at anti-CD38 MoAb relapse was 6% (3/50). An additional selected biallelic event was noted in the Heidelberg set. Each case is described in detail: MM-19 (Calgary) received daratumumab (Dara), lenalidomide, and dexamethasone (Dex), and then had a second exposure to Dara, following relapse. At this second relapse, WGS was performed on two independent extramedullary disease localizations in the liver. Both samples had unique, large deletions (>3 Mb) on chromosome 4 coupled with additional unique focal deletions involving CD38. Flow confirmed the presence of CD138+/CD38- plasma cell populations in both sites. A second case of convergent evolution was seen in a patient (Mayo) after seventh-line Dara/Len/Dex. WGS revealed a common deletion encompassing CD38. On the remaining allele, one subclone had a structural variation-mediated loss of CD38(cancer cell fraction; CCF 40%) while another harbored an L153H missense mutation (CCF 60%). Flow confirmed presence of two CD138+ clones; one with diminished CD38 expression. L153H was functionally validated to abrogate Dara binding, while maintaining CD38 expression. The third sample was collected after fifth-line Dara and sixth-line Dara/Bortezomib/Dex (Heidelberg) and represents another case of convergent evolution with a common loss of chromosome 4 and two subclones detected by WGS (P98Lfs*12, CCF 39% resulting in biallelic inactivation; R140G, CCF 61%). Bulk RNAseq revealed, in comparison to NDMM from CoMMpass, CD38 expression below the first decile. Notably, the R140G variant reduced binding affinity to Dara but not to Isatuximab (Isa). A fourth case (Mayo) relapsing after 6 months of Dara-based therapy similarly had a focal deletion of CD38 with an L18Sfs*16 frameshift deletion on the remaining allele associated with a 4.35-fold reduction in CD38 expression by RNAseq compared to pre-treatment. Interestingly, in 3/4 cases the CD38-negative clone emerged after a second exposure to anti-CD38 MoAb, suggesting emergence with the first exposure and then selection during the second. One other monoallelic event (C275Y; Heidelberg) was detected post-Isa. Copy number of the alternate allele is not known given targeted sequencing methodology, but the variant was seen to reduce binding affinity to both Dara and Isa. The prevalence of monoallelic loss of CD38 in newly diagnosed MM from MMRF CoMMpass was 7% (50/701) and two patients had mutations in CD38, but there were no cases of biallelic loss suggesting these events are driven by treatment pressure. Conclusion Distinction of the mode of CD38 loss is valuable for predicting outcome of re-challenge to anti-CD38 MoAbs. We see here an estimated 5-10% of anti-CD38 MoAb relapses are characterized by genomic antigen escape. Results here are in line with antigen escape to BCMA- and GPRC5D-directed immunotherapies (Lee et al, Nat Med 2023), suggesting that acquired antigen escape resulting from prolonged exposure to the selective pressures of immunotherapy is a recurrent mechanism of resistance in MM. Moreover, in select cases, non-overlapping resistance mutations provide rationale for switching agents at second exposure.

Background: Monoclonal B-cell lymphocytosis (MBL) is a common premalignant condition that can be classified into chronic lymphocytic leukemia (CLL)-type (CD5, CD23, and CD20 co-expression), and non-CLL-type (CD5-negative) based on flow cytometry immunophenotypes. CLL-type MBL is relatively well characterized with increased risk of hematologic (especially lymphoid) malignancies, serious infections, and decreased response to vaccinations. However, the prevalence, genetic predispositions, and clinical consequences of non-CLL-type MBL remain to be fully elucidated. Herein, we report findings in non-CLL-type MBL in a large screening cohort. Methods: We screened 10,330 individuals (40 years or older with no prior history of hematologic cancers) for MBL with highly sensitive flow cytometry, in the Mayo Clinic Biobank, with patients recruited from primary care-based clinics. Individuals were genotyped by Regeneron Genetics. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Among those individuals residing in the 27 counties surrounding Rochester, Minnesota, we identified incident hematologic malignancies and melanomas using ICD codes and the Rochester Epidemiology Project, a population based medical records linkage system. All diagnoses were manually reviewed for confirmation. Cox regression was used to calculate hazard ratios (HR) and 95% CIs, adjusted for age and sex. Results: We identified 1,736 individuals with MBL and 8,341 individuals without MBL. Among those with MBL, 246 (14.2%) had non-CLL-type MBL as the major clone. The median age of non-CLL-type MBL was 77 years (range, 47-95). The prevalence of non-CLL-type MBL significantly increased with age, with 0.5% prevalence among those between 40 to 59 years of age, 2.3% prevalence among those between 60 and 79 years, and 7.8% prevalence in those 80 years of age or older. The prevalence of non-CLL-type MBL was significantly higher (P < .001) among males (3.1%) than females (2.0%). Using the genome-wide complex trait analysis tool, the heritability of non-CLL-type MBL was estimated to be ~26%, suggesting inherited genetics plays a role in etiology. To identify common inherited variants, we first evaluated the known CLL susceptibility variants, many of which also associated with CLL-type MBL. Among the 42 variants, none were statistically associated (all P > .05) with non-CLL-type MBL. We next conducted a genome-wide association study to identify novel variants, and the most significant variant resided in 10q26.3 (minor allele frequency [MAF]control = 0.7%, MAFMBL = 3.5%, OR = 6.4, 95% CI 3.6-11.3, P = 5.7 * 10-8). Among the 6,478 individuals (192 with non-CLL-type MBL and 6,286 without MBL) followed for incident hematologic malignancies, 36 developed lymphoid and 23 developed myeloid malignancies after a median of 4.2 years of follow up. Similar to CLL-type MBL, non-CLL-type MBL was associated with a 5.2-fold risk for lymphoid malignancies (95% CI 2.0-13.0, P < .001). Elevated but not significant risks were observed for myeloid malignancies (HR = 2.3, 95% CI 0.5-10.3, P = 0.27). The lymphoid malignancies developed by individuals with non-CLL-type MBL included diffuse large B-cell lymphoma, multiple myeloma, marginal zone lymphoma, and other low-grade lymphomas. We recently reported that individuals with CLL-type MBL are at higher risk for developing melanoma (HR = 1.9, 95% CI 1.3-2.8). Interestingly, non-CLL-type MBL had no evidence of an association with melanoma risk (HR = 0.4, 95% CI 0.1-1.8). We did not observe a significant difference in overall survival between individuals with non-CLL-type MBL vs those without MBL. Conclusions: In the largest screening cohort to date, we evaluated the etiology and clinical outcomes of non-CLL-type MBL, an understudied yet common condition that occurs in ~4 million adults in the United States. Non-CLL-type MBL is more common in men and older individuals, and evidence supports that inherited genetics plays an etiological role with 10q26.3 locus identified as a novel susceptibility locus. Individuals with non-CLL-type MBL had significantly increased risk of lymphoid malignancies, similar to CLL-type MBL, but, distinct from CLL-type MBL, no evidence of an increased risk of melanoma. Our study provides insights into the prevalence, genetic predispositions and clinical consequences of non-CLL-type MBL and how it differs from the classic CLL-type MBL.

Introduction: Multiple myeloma (MM) with t(11;14)(q13;q32), resulting in CCND1;IGH rearrangement, remains the only subset amenable to targeted therapy with the BCL2 inhibitor, Venetoclax. Nevertheless, not all patients with this translocation respond to treatment and some progress early after initial response. The mechanisms behind this clinical heterogeneity are largely unknown. Methods: To investigate mechanisms of heterogeneity in response to Venetoclax, we interrogated 44 whole genome sequencing (WGS) and 24 whole exome sequencing data from a training cohort including 33 patients with t(11;14) MM treated with Venetoclax at the Mayo Clinic. In a validation cohort, we extracted single nucleotide variants (SNV) and indels from 17 MM patients with t(11;14) with available RNAseq data treated with Venetoclax at Emory University. Results: In the training cohort the median number of lines of therapy was 4 and Venetoclax drug regimens included combinations with dexamethasone, proteasome inhibitors, anti-CD38 agents, and immunomodulatory agents. After a median follow up of 30 months, the median progression free survival (PFS) was 9 months, and 13 patients (39%) achieved a complete remission (CR) to therapy, while 9 (27%) were refractory (i.e. no response or progressed within the first 60 days). For our driver discovery analysis, we integrated SNV, indels, copy number variants (CNV) and structural variants (SV). Deletion 17p (del17p) or TP53 mutation, MCL1 gain/amp, and hyper-APOBEC did not correlate with PFS in our cohort (p=0.59, p=0.71, and p=0.05, respectively). In contrast, mutations involving the MAPK pathway (NRAS, KRAS, BRAF, FGFR3) displayed the strongest association with short PFS, which was retained after adjusting for del17p/TP53 mutations, and 1q gain/amp (HR: 4.76, CI: 1.8-12.9, p=0.002). The association between suboptimal response and resistance to Venetoclax in MM patients with MAPK pathway mutations (n=7) was also observed in our validation cohort from Emory University (p=0.046). This is in line with what has been reported in functional studies in acute myeloid leukemia, where RAS pathway promotes MCL1 under the selective pressure of Venetoclax (Zhang et al. Signal Transduct Target Ther 2022). While MAPK mutations explained most cases of refractory patients (86%), the mechanisms through which initially responding patients subsequently progressed remain unclear. To investigate this, we interrogated 10 patients with pre- and post-treatment samples collected at disease progression (PD). We identified genomic events acquired/selected at PD of genes involved in the anti-apoptotic MCL1/BCL2 and MAPK pathways in 6/10 cases. One patient was refractory to Venetoclax and had double minutes involving BCL2, which caused multiple focal amplifications (>30 copies), likely reflecting extrachromosomal DNA. Another patient progressed with an acquired, focal deletion on BCL2L11. In two patients, the progressing clones harbored mutations in MAPK pathways (one BRAF and the other WNK2). Two patients who relapsed after sustained CR had major branching evolution where a previously undetectable focal deletion of PMAIP1 (NOXA) was seen at PD. Importantly, one of these deletions was a biallelic event. Venetoclax sensitivity assays conducted on PMAIP1 knockout in the MM t(11;14) cell line, KMS12PE, demonstrated significantly reduced cell death across all drug concentrations (p<0.05), validating the loss of NOXA as a mechanism of resistance to Venetoclax in MM. Finally, similarly with what was previously reported in chronic lymphocyte leukemia treated with Venetoclax (e.g., Herling et al Nat Comm 2018), CDKN2A/B monoallelic and biallelic focal deletions were selected post-Venetoclax, at PD in 4 patients while TP53 SNV/deletions were selected from subclonal to clonal in 2 patients at PD. Conclusion: In this study, we demonstrate that mutations involving the MAPK pathway identify a population of MM patients with t(11;14)(CCND1;IGH) who are unlikely to respond to Venetoclax. We also show that genomic events involving the MAPK and BCL2/MCL1 pathways are primarily responsible for disease progression after initial response. These findings highlight the critical role of comprehensive genomic profiling through WGS in identifying mechanisms of resistance to targeted therapies in MM, aiding patient selection for Venetoclax, and guiding the development of novel drug combinations.

Background: Immune dysfunction is a hallmark of chronic lymphocytic leukemia (CLL), and several studies suggest that hypogammaglobulinemia at the time of CLL diagnosis may predict shorter time to first therapy (TTFT; Parikh et al, 2015). However, it is unclear if hypogammaglobulinemia can predict TTFT independent of the CLL-International Prognostic Index (CLL-IPI) and the tumor mutational load (TML). We sought to determine if serum immunoglobulin (Ig) levels at diagnosis can independently predict outcomes in patients with CLL and its precursor, high-count monoclonal B-cell lymphocytosis (MBL). Methods: We used the Mayo Clinic CLL Database to identify newly diagnosed CLL/MBL consented within 2 years of diagnosis. Next generation sequencing (NGS) was performed using a 59 gene panel (SureSelect) to detect mutated CLL genes. Serum Ig levels were quantitated by radial immunodiffusion (Immunoplates). Serum Ig levels below the lower limit of normal were defined as hypogammaglobulinemia. We used Chi-squared or Fisher's exact tests to compare discrete variables and the Kruskal Wallis test for continuous variables. TTFT was calculated from date of diagnosis to date of first treatment or date last known to be untreated. We used Cox regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for TTFT and overall survival (OS) associations. Results: We identified 895 newly diagnosed patients (588 CLL and 307 MBL) who had undergone NGS for recurrent somatic mutations and had serum Ig levels available at the time of diagnosis. The median absolute B-cell count was 6.7 x109/L; 613 (69%) were male, 375 (42%) had unmutated IGHV genes, and 81 (9%) had TP53 disruption. Based on CLL-IPI, 412 (46%) patients had low risk, 285 (32%) had intermediate risk, and 198 (22%) had high/very high risk. The most commonly mutated genes were NOTCH1 (14%), SF3B1 (11%), TP53 (10%), NFKBIE (8%), ATM (6%), and BIRC3 (6%). Based on the number of genes with high impact or hotspot mutations, we classified 40%, 29%, and 31% patients with TML scores of 0, 1, or 2+, respectively. At CLL/MBL diagnosis, the median serum IgG, IgA, and IgM was 891 mg/dL (111-2620), 121 mg/dL (1-1150), and 44 mg/dL (5-1320), respectively; 284 (32%) patients had low serum IgG, 128 (17%) had low serum IgA, and 322 (42%) had low serum IgM. Low serum IgA was associated with mutations in NOTCH1, SF3B1, ATM, BIRC3 and NFKBIE, low IgM was associated with SF3B1 and BIRC3 mutations (all comparisons, P < .05), and low IgG was associated with BIRC3 mutation (P = .05). The median follow-up for the cohort was 8.4 years. A total of 359 patients progressed requiring therapy (299 CLL and 60 MBL), and 265 patients died (177 CLL and 88 MBL). The TTFT was significantly shorter for patients with low serum IgG compared to normal serum IgG at diagnosis (median 4.4 vs 10.8 years) and similar findings were observed in individuals with low serum IgA (1.7 vs 10.8 years) and those with low serum IgM (4.6 vs 11.1 years). On univariable analysis, the following were associated with a shorter TTFT: low serum IgG (HR: 1.8, 95% CI 1.4-2.2), low serum IgA (HR: 3.5, 95% CI 2.7-4.4), low serum IgM (HR: 1.9, 95% CI 1.5-2.3), TML score of 1 (HR: 1.7, 95% CI 1.3-2.2), TML of 2+ (HR: 4.4, 95% CI 3.4-5.7), CLL-IPI intermediate risk (HR: 3.7, 95% CI 2.8-4.8), and CLL-IPI high/very high risk (HR: 7.6, 95% CI 5.8-10.1) (all comparisons, P < .0001). In multivariable analyses after adjusting for CLL-IPI and TML, low serum levels of IgG (HR: 1.8, 95% CI: 1.5-2.3), IgA (HR: 2.4, 95% CI: 1.8-3.0), and IgM (HR: 1.8, 95% CI: 1.5-2.3) were independently associated with a shorter TTFT (all comparisons, P < .0001). No significant differences in OS were detected between patients with low and normal serum IgG, IgM, and IgA on univariate analysis. Similar results for TTFT and OS were observed after adding clonal B-cell count to the Cox model. Conclusions: Our study shows that low serum Ig levels at the time of initial diagnosis is predictive of TTFT independently of the CLL-IPI and TML, which can further risk stratify individuals with newly diagnosed CLL/MBL. Additionally, low serum IgA, IgM, and IgG levels at diagnosis are significantly associated with specific somatic mutations in individuals with newly diagnosed CLL/MBL.

Background: Multiple myeloma (MM) is an incurable plasma cell malignancy with several recurrent primary cytogenetic abnormalities including t(11,14)(q13;q32), resulting in increased expression of CCND1. Although 40% of patients with t(11;14) MM respond favorably to the BCL-2 inhibitor venetoclax (ven), 60% fail to achieve an objective response. Previous studies have suggested upregulation of PMAIP1, encoding the pro-apoptotic NOXA protein, in association with ven sensitivity (Leblay, Blood, 2024). Here we sought to evaluate the role for NOXA in ven response in t(11;14) MM. Methods: We analyzed whole exome, whole genome and RNA sequencing data of CD138+ tumor cells from 1,359 patients with MM in the CoMMpass and Mayo Clinic cohorts. The histone mark H3K27ac was assessed in t(11;14)+ KMS12PE and U266B1 human myeloma cell lines using ChIP-seq. Gene expression and protein levels of PMAIP1 (NOXA) were also detected by qRT-PCR and Western blot, respectively. Live cell detection was evaluated by Annexin V and PI using flow cytometry and by YoYo3 dye using IncuCyte following treatment with ven. Results: The t(11;14) MM subtype was associated with a 1.6-5.5-fold increased expression of PMAIP1 compared to non-t(11,14) MM in the CoMMpass and Mayo Clinic cohorts (p<0.0001). These findings are consistent with a significant correlation between CCND1 and PMAIP1 expression in both cohorts (p<0.001). To evaluate whether cyclin D1 is directly associated with increasing expression of PMAIP1 in t(11;14) MM, we generated a cyclin D1 knock out (KO) version of the U266B1 line using CRISPR/cas9. Two independent U266B1 KO cell lines showed reduced expression of PMAIP1 and NOXA protein compared to the WT U266B1 consistent with a role for cyclin D1 in regulating PMAIP1 expression. In addition, the KO lines had a 4.8-fold reduced H3K27ac signal at the PMAIP1 locus in comparison to the WT control suggesting a role for cyclin D1 in promoting H3K27ac. We previously observed increased expression of CCND1 in association with the G allele compared to the A allele of rs9344 in both the CoMMpass and Mayo Clinic cohorts (Leukemia, in press). We next asked whether the G risk allele of rs9344 within CCND1 impacted gene expression and H3k27ac levels at the PMAIP1 locus. To study the allele specific effect of CCND1, we used the CRISPR/cas9 edited KMS12PE WT (rs9344 AA) and Knock In (KI) (rs9344 GG) cell lines. In the KI cells, the G allele was accompanied by an ~2-fold increased expression of PMAIP1 and NOXA protein detected by western blot (<0.05). Increased PMAIP1 expression in KI lines was also accompanied by a 2.2-fold increased H3K27ac further suggesting a role for the G allele of rs9344 in promoting H3K27ac and gene expression at the PMAIP1 locus. Finally, we evaluated the consequence of PMAIP1 expression levels in relation to ven response. Increased expression of PMAIP1 and elevated levels of NOXA protein in the KMS12PE KI lines was associated with increased sensitivity to ven resulting in ~50% increase in cell death of KI vs. WT lines (p<0.02). Conclusion: We demonstrate a role for cyclin D1 in regulating PMAIP1 expression through H3K27ac resulting in increased NOXA protein. Levels of NOXA are associated with the G allele of rs9344 and increased venetoclax response.

Background: Chronic lymphocytic leukemia (CLL) is one of the most common leukemias in the United States (US), with an estimated incidence of 20,700 cases in 2024. Population-based studies conducted in Olmsted County (OC), Minnesota, reported an increase in CLL incidence over a 55-year period (1935 through 1989), which was attributed to increased detection of early-stage cases in individuals over age 50. In an updated study from 2000 to 2010 in OC, we reported no further increase in the incidence of CLL using the 2008 International Workshop on CLL (IWCLL 2008) guidelines. Additionally, when we compared the 1996 National Cancer Institute Working Group (NCI-WG 96) criteria with the IWCLL 2008 guidelines incorporating the monoclonal B-cell lymphocytosis designation, we found that the use of the IWCLL 2008 guidelines reduced the overt CLL incidence, modified the distribution of Rai Stages at diagnosis, and shortened the median time to first treatment (TTFT). Over the last two decades, 5-year survival for CLL has increased due to new therapies and improved supportive care. Herein we provide updated data on the incidence, prevalence, and natural progression of CLL in OC and studied similar features in an expanded 6-county region around Mayo Clinic, Rochester. Methods: We utilized the Rochester Epidemiology Project (REP), a population-based medical records linkage system, to identify CLL cases from 2000 to 2020 using International Classification of Diseases codes. We included individuals residing in OC and the expanded 7-county region (REP-7, including OC), for which data became available starting in 2010. CLL cases were included if they met the IWCLL 2008 criteria. All CLL diagnoses were confirmed by manual review of medical records and data on CLL-related factors were collected, including treatment status, prior history of other cancers, and overall survival (OS). Age- and sex-specific incidence rates were derived from REP census data and adjusted to the 2020 US white population. OS analyses were conducted using the Kaplan-Meier method. TTFT was calculated with death as a competing risk. Cox regression was used to evaluate the effect of decade on TTFT and OS among patients in OC. Results: From 2010-2020, 242 patients were diagnosed with CLL in REP-7, and between 2000-2020, 174 patients were diagnosed in OC. Among the 242 CLL patients in REP-7, the median age at diagnosis was 71 years (range 30-96),158 (65%) were male, 236 (99%) were white, and 74 (31%) had a history of other cancer, with non-melanoma skin cancer being the most common. After age and sex adjustment to the 2020 US white population, the CLL incidence in REP-7 and OC from 2010 to 2020 was 7.8 and 8.1 per 100,000 individuals, respectively. In OC, the age- and sex-adjusted CLL incidence from 2000 to 2020 was 7.5 per 100,000 individuals, showing no significant trend changes over time. Among the 242 patients diagnosed with CLL in REP-7 between 2010 and 2020, 102 received treatment with a median TTFT of 9.8 years. In OC from 2000 to 2020, 79 out of 174 patients received treatment with a median TTFT of 9.4 years. When comparing across two decades in OC, the median TTFT was 12.8 years for the period from 2000 to 2009, decreasing to 5.9 years between 2010 and 2020, but was not significantly different after accounting for age and sex. The median OS for patients diagnosed with CLL in REP-7 from 2010 to 2020 was 11 years. The median OS for patients in OC between 2000 to 2020 was 12.6 years. When comparing across decades in OC, the median OS was 8.9 years for the period from 2000 to 2009, increasing to 13.6 years for those diagnosed between 2010 and 2020, but this trend was not significantly different after accounting for age and sex. Prevalence data for REP-7 and OC were analyzed as of 01/01/2020. Additionally, prevalence data for OC were available for 01/01/2010. In REP-7, the prevalence of CLL was 97 per 100,000 individuals on 01/01/2020. In OC, the CLL prevalence was 86 per 100,000 individuals on 01/01/2010, and 100 per 100,000 individuals on 01/01/2020. Conclusions: We observed no substantial increase in the incidence rates of CLL in REP-7 from 2010 to 2020, nor in the 20 years of available data from OC. However, the prevalence of CLL in OC increased by 16% from 2010 to 2020. Additionally, evidence indicated a longer OS among newly diagnosed CLL cases in OC compared to previous decades, suggesting a potential impact of improved therapeutic approaches on OS.

Introduction: Cyclin D1, encoded by the CCND1 gene, is primarily known as a major regulator of cell cycle progression. Emerging studies have demonstrated a role for cyclin D1 in regulating gene transcription. In mantle cell lymphoma with t(11;14), cyclin D1 strongly binds to the promoters of highly transcribed genes and its overexpression led to genome-wide transcriptional downregulation. However, the transcriptional role for cyclin D1 in t(11;14) MM remains undefined. This study aims to elucidate the global chromatin and transcriptional changes induced by cyclin D1 in t(11;14) MM. Methods: Knockout (KO) of CCND1 in U266B1, a MM cell line carrying t(11;14), was performed by the Synthego (Redwood City, CA) using CRISPR/Cas9, which generated a frameshift insertion of base pair (bp) “A” after 69,641,478 bp within exon 1. Two independent KO clones were used for downstream analysis. ChIP-seq and ATAC-seq were performed for KO (n=2) and wild-type clones (WT, n=2). ChIP-seq used antibodies against H3K4me3 (promoter) and H3K4me1 and H3K27ac (enhancer). Libraries were sequenced to 2x51 bp on the Illumina NextSeq 2000 platform. Total RNA libraries (in triplicate) were sequenced to 2x150 bp. ChIP-seq and ATAC-seq reads were mapped to hg38 with BWA, and peaks were identified by MACS2 at FDR 0.01. RNA-seq reads were mapped to hg38 with STAR and gene raw counts were estimated using featureCounts. Differential peaks and differentially expressed genes between KO and WT were identified using edgeR at FDR 0.05 and fold-change 2. Only peaks present in at least 2 samples and genes having reads-per-million value 0.1 in at least 1 sample were included, excluding those from chrY. Pathway analysis was performed with the Enrichr package. Results: Differential expression analysis identified 2-fold more genes that were down-regulated than up-regulated in KO compared to WT (467 vs. 233). Pathway analysis revealed an enrichment of TNF-alpha signaling via NF-kappaB, inflammatory response, and apoptosis in association with down-regulated genes, while up-regulated genes are enriched in myogenesis, KRAS signaling, and epithelial mesenchymal transition pathway. In parallel, differential chromatin analysis identified 2.5% (1,161) of the peaks from H3K27ac and 0.5-1% (359-610) of the peaks for the other 3 marks that showed significant changes of signal levels between KO and WT. For the promoter mark H3K4me3, there are largely a comparable number of peaks showing increased (191 peaks) and decreased (168 peaks) signals, respectively, in KO compared to WT. However, for ATAC-seq and H3K4me1, about 70% of the differential peaks showed an increased signal in KO, in contrast with H3K27ac for which 60% of differential peaks showed decreased signal in KO. The preferential loss of H3K27ac occupancy is consistent with the observation of a higher proportion of genes being down-regulated in KO. These results indicated that knockout of CCND1 led to a global increase of chromatin accessibility, but a reduction of H3K27ac and gene expression. Further annotation revealed that, for H3K27ac and H3K4me1/ATAC-seq, 92% and 97% of the differential peaks were located >2kb away from the transcription start sites (TSSs) of expressed genes. Even for H3K4me3, 88% of the differential peaks are >2kb away from TSSs. The results suggest that the changes of chromatin upon CCND1 knockout predominantly occurred at distal regulatory regions. For up- and down-regulated genes, we estimated ChIP-seq and ATAC-seq signals in the TSS 2kb regions, revealing that up-regulated genes were overall associated with increased signals. Further analysis found that, for ATAC-seq and H3K4me3, 44% and 59% of the differential peaks located in the promoter regions were associated with differentially expressed genes, which dropped to ~20% for the two enhancer marks. These data indicate that changes in chromatin state impact the expression of the associated genes. Conclusion: Knockout of CCND1 impacts local chromatin state, particularly at enhancer regions, and the transcriptional program. In promoter-proximal regulatory regions, the changes of chromatin accessibility and H3K4me3 are strongly correlated with the expression changes of nearby genes. Further studies will identify the genes targeted by the differential enhancers and their possible roles in MM pathogenesis mediated by cyclin D1 following the t(11;14) event.

Introduction MYD88 L265P mutations (MUT) are found in 90% of patients with Waldenstrom macroglobulinemia (WM) and lead to the activation of pro-survival signaling. Limited literature has demonstrated that patients with MYD88 wild-type (WT) WM have poorer responses to conventional therapies and a higher likelihood of histological transformation. Increasingly the role of non-coding genomics in lymphomagenesis has demonstrated importance and currently, the epigenetic landscape of MYD88 WT WM remains underexplored. This study utilized a multi-omics approach to assess DNA methylation in MYD88 WT WM compared to MYD88 MUT patients and controls. Methods Bone marrow samples were prospectively collected from 20 subjects [3 MYD88 WT/CXCR4 WT WM, 12 MYD88 MUT WM and 5 normal controls] and positive selection for CD19+ and/or CD138+ B cells was performed. The analysis focused on the comparison between MYD88 WT WM to normal controls and MYD88-WT to MUT WM. Genome-wide DNA methylation analysis was performed, and CpG methylation ratios were segmented into 200-bp regions. Differentially methylated regions (DMRs) were determined using a Q-value <0.05 and absolute methylation difference ≥10%. Next, RNA-seq data was interrogated and differentially expressed genes (log2fold change>0.5 or <-0.5, false discovery rate<0.05) with corresponding promoter hypo- or hypermethylation were assessed. Ingenuity Pathway Analysis was used to determine potentially epigenetically regulated pathways using a corrected p<0.05. Results Comparing MYD88 WT WM samples to controls, a statistically greater number of hypermethylated DMRs were observed across genomic regions (hyper- vs hypomethylated): promoter (29129 vs 6300), 3'untranslated region (UTR) (2068 vs 1551), 5'UTR (10169 vs 735), intron (30344 vs 24511), exon (2224 vs 6855), CpG island (30372 vs 1240) and CpG shore (14477 vs 5872) (all p<0.001). Next, mRNA analysis demonstrated that 214 genes were differentially up-/downregulated with concordant promoter hypo-/hypermethylation. On the assessment of promoter-methylation-based pathways, multiple pathways involved in cell signaling, proliferation, gene regulation, cytokine and immune signaling and metabolism were found to be deregulated. Most pathways were found to be negatively enriched in MYD88 WT patients compared to controls including several members of interleukin (IL) signaling (IL-2, IL-3, IL,7) and regulators of cell proliferation and signaling including mTOR, PI3K/AKT, and p53. Positive enrichment of the PTEN, FGFR3 and ERK/MAPK pathways were observed. Underlying multiple pathways was promoter hypermethylation with corresponding mRNA downregulation of AKT3 and PIK3CA. Furthermore, promoter hypomethylation with transcript upregulation of PTEN was observed. Next, in comparison of MYD88 WT to MUT samples, a statistically greater number of hypermethylated DMRs were observed across genomic regions (hyper- vs hypomethylated): promoter (24829 vs 2544), 3'UTR (3534 vs 378), 5'UTR (6111 vs 605), intron (52292 vs 5360), exon (21066 vs 2187), CpG island (18334 vs 1551) and CpG shore (16351 vs 1957) (all p<0.001). Assessing mRNA expression, 354 genes were up-/downregulated with concordant promoter hypo-/hypermethylation. Pathway analysis revealed positive enrichment of the mitotic G1 phase to G1/S transition pathway with upregulation of the CDK4, LYN, and PSMA6 transcripts in MYD88 WT compared to MUT samples. Further, negative enrichment of cell signaling regulators including the RHO GTPase family, IL pathways (IL-4, IL-8, IL-13, IL-17) and integrin signaling were found in MYD88 WT compared to MUT samples. Underlying multiple pathways were downregulation with promoter hypermethylation of BCL2, an important regulator of apoptosis. Conclusion Genome-wide methylation analysis of MYD88 WT WM revealed a unique DNA methylation profile, with potential epigenetic alteration of key oncogenic transcripts and pathways. The negative enrichment of cell signaling and cytokine/inflammatory pathways in MYD88 WT samples compared to controls and MYD88 MUT WM suggests an epigenetic mechanism for immune evasion and potentially underlies the reduced response to conventional WM-directed therapies. Overall, this study reveals that MYD88 WT WM is epigenetically distinct from MYD88 MUT WM and regulators such as methylation may underlie the unique clinical and phenotypic properties observed.

PURPOSE Chronic lymphocytic leukemia (CLL)–phenotype monoclonal B-cell lymphocytosis (MBL) is a premalignant condition that is roughly 500-fold more common than CLL. It is unknown whether the two-fold increased risk of developing melanoma associated with CLL extends to individuals with MBL. METHODS Using the Mayo Clinic Biobank, we identified participants who were 40 years or older with no previous hematological malignancies, who resided in the 27 counties around Mayo Clinic, and who had available biospecimens for screening. Eight-color flow cytometry was used to screen for MBL. Individuals with MBL were classified as low-count MBL (LC-MBL) or high-count MBL on the basis of clonal B-cell percent. Incident melanomas were identified using International Classification of Diseases codes and confirmed via medical records review. Cox regression models were used to estimate hazard ratios (HRs) and 95% CI. RESULTS Of the 7,334 participants screened, 1,151 were identified with a CD5-positive MBL, of whom 1,098 had LC-MBL. After a median follow-up of 3.2 years (range, 0-13.5), 131 participants developed melanoma, of whom 36 individuals were positive for MBL. The estimated 5-year cumulative incidence of melanoma was 3.4% and 2.0% among those with and without MBL, respectively. After adjusting for age, sex, and history of previous melanoma, individuals with MBL exhibited a 1.86-fold (95% CI, 1.25 to 2.78) risk of melanoma. This elevated risk persisted when analysis was restricted to those without a history of melanoma (HR, 2.05 [95% CI, 1.30 to 3.23]). Individuals with LC-MBL had a 1.92-fold (95% CI, 1.29 to 2.87) increased risk of developing melanoma overall and a 2.74-fold increased risk (95% CI, 1.50 to 5.03) of melanoma in situ compared with those without MBL. CONCLUSION LC-MBL is associated with an approximately two-fold increased risk of melanoma overall and a 2.74-fold increased risk of melanoma in situ.