Citations

Abstract Studies have reported conflicting findings regarding the contribution of germline variants or somatic genomic drivers to racial disparities in multiple myeloma. To comprehensively investigate somatic drivers in relation to inherited genetics in multiple myeloma, we combined newly sequenced whole-genome sequencing data with publicly available datasets (total n = 1,286). Overall, we did not identify germline or somatic genomic differences that explain the different risk of developing multiple myeloma between patients with genetic similarity to African (AFR) or European (EUR) reference populations. A difference in the detectability and timing of APOBEC-associated and germinal center mutational activity was observed. Integrating epidemiologic data and mutational signature-based temporal estimates, we challenge the assumption that individuals in the AFR group develop multiple myeloma at a younger age. Finally, we demonstrate that, with equal access to efficacious therapies, patients in the AFR and EUR groups have equivalent clinical outcomes. Significance: Multiple myeloma is reported to occur at higher rates in individuals who self-identify as non-Hispanic Black. In this large dataset, genomic drivers occur at the same rate among ancestry groups, except for APOBEC mutagenesis. With equivalent therapy, clinical outcomes did not differ for patients grouped by genetic ancestry similarity.

Abstract Background: CLL incidence differs by ancestry, with higher rates in individuals of European ancestry (EA) than African ancestry (AA). The contribution of epigenetic variation to these disparities remains unclear. We investigated chromatin accessibility, histone modifications, and gene expression in CLL patients of EA and AA to identify ancestry-associated epigenetic signatures. Methods: ATAC-seq was performed on blood clonal B cells from 15 AA and 51 EA patients; CUT&Tag profiling of six histone marks (H3K4me1, H3K4me3, H3K27ac, H3K9me3, H3K27me3, H3K36me3) was done in 8 AA and 12 EA cases; RNA-seq was generated for all samples. Ancestry was genetically confirmed. Differential regions were linked to genes using ENCODE and Hi-C datasets. Results: ATAC-seq revealed 4,955 regions with increased and 3,000 with decreased accessibility in AA vs EA (fold change >2, FDR ≤0.05). Regions with increased accessibility mapped to 2,313 genes; 599 were upregulated and enriched in TNF-α/NF-κB (p=1.40E-53) and apoptosis pathways (p=6.90E-11). Regions with decreased accessibility linked to 1,953 genes; 317 were downregulated, enriched in interferon-γ response (p=1.04E-03) and IL-6/JAK/STAT3 signaling (p=6.41E-03). Histone profiling showed ancestry-specific changes across all marks. For H3K27ac (active enhancer), 68 genes with a decreased signal in AA were downregulated and enriched in interferon-γ (p=1.55E-06) and IL-6/JAK/STAT3 (p=9.62E-05) pathways. Increased H3K27me3 (repressive histone modification) in AA affected 399 genes enriched in IL-6/JAK/STAT3 (p=8.13E-09) and interferon-γ (p=2.48E-05) pathways. These findings suggest coordinated regulation by active and repressive chromatin states. Conclusions: Despite limited sample size, this study suggests differences between AA and EA ancestries may exist in epigenetic signatures in CLL, implicating chromatin accessibility and histone modifications in pathways relevant to tumor biology. These signatures may also contribute to the differences in CLL incidence by ancestry. Future studies with more patients and functional validation experiments are needed to clarify the role of these epigenetic differences in CLL pathogenesis. Citation Format: Zhiquan Wang, Huihuang Yan, Mingma Sherpa, Daniel Bihnam, Cristine Allmer, Chantal McCabe, Daniel O’Brien, Dennis P. Robinson, James R. Cerhan, Neil E. Kay, Sameer A. Parikh, Esteban Braggio, Susan Slager. Ancestry-associated epigenetic differences in chronic lymphocytic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 5064.

Abstract Introduction: ASXL1-mutant (mt) Chronic Myelomonocytic Leukemia (CMML) is an aggressive myeloid neoplasm characterized by increased proliferation, resistance to epigenetic therapies, and poor survival. Epigenetic remodeling and up-regulation of leukemogenic driver genes are hallmarks of ASXL1mt CMML. Truncating mutations in ASXL1 activate enhancers that drive the expression of key leukemogenic driver genes such as HOXA9 and its co-factor MEIS1. Here, we investigate the chromatin accessibility of these leukemogenic enhancers by ASXL1 genotype at single-cell resolution (GoT-ChA-seq). Methods: After IRB approval, we performed bulk RNA-seq and ChIP-seq (H3K4me1, H3K4me3, H3K27ac, H3K27me3, H2AK119ub, H2BK20ac) on bone marrow (BM) mononuclear cells from 40 patients with CMML (19 ASXL1mt, 21 -wildtype/wt). After quality control, differential analyses of gene expression and histone modification occupancy between ASXL1mt and wt patients were conducted. A consensus peak calling framework was used to define genomic regions with ASXL1mt-specific histone modifications. Canonical enhancers were defined as genomic regions with a co-occupancy of H3K4me1 and H3K27ac. We then added regions marked by H2BK20ac to capture regulatory regions associated with p300 (Narita et al. Nat Genet 2023). We used an in silico prediction strategy to associate regulatory regions marked by strong histone modification occupancy with highly up-regulated leukemogenic driver genes (HOXA6-10, MEIS1). The ReMap database (Hammal et al. Nucleic Acids Res 2022) was queried for transcription factor binding in the genomic regions of interest. We then performed GoT-ChA-seq (Izzo et al. Nature 2024) to investigate chromatin accessibility for ASXL1mt (n=3) and wildtype BM single cells from CMML patients. In vitro colony forming assays were conducted on primary CMML BM cells (11- ASXL1mt and 9 wt CMML), assessing the efficacy of EP31670, a novel, oral, dual inhibitor of p300/BRD4. Results: This study included 40 patients with CMML, 19 ASXL1mt and 21 wt. The variant allele frequencies for ASXL1 were compatible with heterozygosity (median 0.41) and TET2 mutations were balanced between the two groups. We identified 149 candidate regulatory regions correlating with the expression of up-regulated target genes of interest (HOXA6-10, MEIS1) by bulk ChIP- and RNA-seq. Publicly available ChIP-seq data from myeloid cell lines demonstrated enrichment of p300 (in THP-1 cells, enrichment effect size +2.58, p=0.002) and BRD4 (in K-562 cells, enrichment effect size +1.58, p=0.045) in these candidate regulatory regions. Next, we interrogated the chromatin accessibility of ASXL1mt (n=13620) and ASXL1wt (n=24269) single cells from 3 CMML patients using GoT-ChA-seq. Differential accessibility analysis revealed 144 genomic regions, all of which were preferentially accessible in ASXL1mt cells compared to ASXL1wt cells (p<0.05 for all regions). Four of these preferentially accessible genomic regions (p<0.001 for all regions) overlapped with the previously identified candidate regulatory regions for MEIS1. The H3K27ac and H2BK20ac occupancy levels in these regions was strongly associated with the expression of MEIS1 (median Pearson r=0.73, p=1.16x10-7). These candidate regulatory regions were annotated in ENCODE as distal and proximal enhancers. Targeting this biology with the dual p300/BRD4 inhibitor EP31670 in progenitor colony forming assays led to a preferential therapeutic effect in ASXL1mt BM samples (median IC50 = 25nM) compared to in ASXL1wt samples (median IC50 = 74nM, p=0.03). Treatment with EP31670 was also associated with increased differentiation (increased blast forming units-erythroid; BFU-E) and organization of colonies in ASXL1mt vs wt CMML. Conclusions: ASXL1mt CMML is characterized by the up-regulation of several leukemogenic driver genes including HOXA9 and its co-factor MEIS1, secondary to the activation of genotype-specific enhancers. The preferential accessibility of these enhancers can be defined between ASXL1mt and ASXL1wt cells within patients. The enhancers of interest interact with p300 and BRD4 in myeloid cell lines and targeting them with the novel dual p300/BRD4 inhibitor, EP31670, induces therapeutic responses at low nanomolar drug concentrations in vitro, preferentially in ASXL1-mutant cells. These data solidify the biological rationale for targeting p300/BRD4 in ASXL1mt chronic myeloid neoplasms for therapeutic benefit (NCT05488548).

Abstract Introduction: Protein truncating somatic ASXL1 mutations in CMML are associated with high-risk features and poor outcomes. We have previously demonstrated that ASXL1-mutant (mt) CMML is characterized by permissive chromatin states and overexpression of leukemia stemness genes (Binder et al. Nat Commun; 2022). DLK1 (delta-like homolog 1) is a pseudo-NOTCH like ligand that is maternally imprinted and paternally expressed and is physiologically not expressed beyond embryogenesis. Aberrant DLK1 expression has been associated with poor outcomes in solid tumors and we have shown differential overexpression in ASXL1mt versus wild-type (wt) CMML. Both the promoter and the last exon of DLK1 contain CpG islands that play a role in regulating DLK1. The CpG island within the fifth exon of DLK1, known as the DLK1-differentially methylated region (DMR) [chr14:100,726,892-100,738,224] is completely unmethylated on the maternal allele and partially methylated on the paternal allele. Previous work (Developmental Cell; 2021) has also identified a DLK1 imprinting control region (ICR) known as the IG-DMR [chr14:100,811,001–100,811,037]. DNA methylation state of the IG-DMR has been shown to impact the expression of MEG3 mRNA (lncRNA) which functions as a cofactor to direct PRC2 to the DLK1 promoter, suppressing DLK1 expression. Since the regulation of DLK1 in ASXL1mt CMML is not well understood, we conducted this multiomic analyses to define epigenetic mechanisms of DLK1 overexpression in ASXL1mt CMML. Methods: After IRB approval, we performed bulk RNA-seq, ChIP-seq (H3K4me1/me3, H3K27ac, H3K27me3) and Illumina MethylationEPIC array on bone marrow (BM) mononuclear cells of 40 CMML patients (19 ASXL1mt, 21 wt). scRNA-seq was conducted on 5 ASXL1mt and 3 wt cases. All libraries were sequenced on an Illumina NovaSeq instrument. We assessed these samples for differential gene expression (RNA-seq) and changes in histone modifications (ChIP-seq) between ASXL1mt and wt CMML, specifically at the DLK1-DMR and IG-DMRregulatory loci. We conducted ddPCR (digital droplet PCR) to determine allelic status of DLK1 in ASXL1mtand wt cases. Global and allele/sequence specific DNA methylation changes were determined using Biomodal's duet EvoC methodology (6 base sequencing, including 5mC and 5hmC) in 5 ASXL1mt and 5 wt CMML cases. Results: We validated the overexpression of DLK1 in 19 ASXL1mt versus 21 wt (log2 fold change +4.85, FDR 4.67e-13) CMML BM samples. scRNA-seq on 5 ASXL1mt versus 3 wt samples demonstrated that DLK1 is only expressed in mutant cases, with the expression being limited to HSPCs. ChIP-seq analysis at the DLK1 regulatory loci validated (a) reduced occupancy of H3K27me3, correlating with ASXL1mt associated disruption of the PRC2-complex, and (b) increased co-occupancy of H3K4me1 and H3K27ac in mt cases, validating a permissive chromatin state in ASXL1mt cases at the DLK1 regulatory regions. Through correlation analysis of RNA-seq and ChIP-seq, we defined the promoter for DLK1, characterized by high Pearson correlation coefficients (R > 0.5) and enrichment of H3K27ac and H3K4me3. Differential methylation analysis identified global hypomethylation around the DLK1-DMR and promoter regions. We observed a mean reduction of ~40% in methylation status at the DLK1 promoter in ASXL1mt versus wt cases, suggesting a loss of epigenetic silencing at this locus. This was corroborated by significant enrichment of 5hmC at the fifth exon of DLK1 indicating active DNA demethylation. Conversely, we observed a gain in methylation at the IG-DMR, suggesting loss of imprinting control through MEG3 silencing and DLK1 overexpression. To further determine the impact of imprinting deregulation on overexpression of DLK1 in ASXL1mt CMML, allele-specific methylation was profiled using the Biomodal whole genome methylation data. This demonstrated that the methylation patterns in ASXL1mt CMML cases were significantly deregulated at the DLK1-DMR, with >30% loss of methylation compared to GRCh38 reference alleles. Furthermore, ddPCR was able to confirm bi-allelic expression of DLK1 in ASXL1mt CMML, suggesting that the observed methylation dysregulation leads to loss of imprinting of the maternal allele. Conclusion: In summary, we demonstrate that ASXL1 mutations in CMML disrupt epigenetic pathways, promoting DNA demethylation, resulting in permissive chromatin states and an aberrant and oncogenic re-expression of maternal DLK1, due to loss of imprinting control.

Abstract Cyclin D1, encoded by CCND1, is a cell cycle regulator, overexpressed in ~15% of multiple myeloma (MM) cases via CCND1::IGH super-enhancer juxtaposition through t(11;14). Cyclin D1 also plays a non-canonical role in transcriptional regulation, but in t(11;14) MM its function in gene regulation remains undefined. We investigated the global transcriptional and epigenomic consequences of CCND1overexpression in t(11;14) MM. To define a CCND1 signature, RNA-seq from two independent primary MM patient cohorts, CoMMpass and Mayo, was utilized to compare high CCND1 t(11;14) cases (≥300 transcripts-per-million [TPM], n=271) vs. low CCND1 hyperdiploid cases (<30 TPM, n=233). Two CCND1 knockout (KO) clones in t(11;14) MM U266B1 were also generated by CRISPR/Cas9, introducing a frameshift insertion (A) after chr11:69641478 (hg38) in exon 1. ATAC-seq and ChIP-seq (H3K4me3, H3K4me1, and H3K27ac) were performed in duplicates, and RNA-seq in triplicates in U266B1 wild-type (WT) and KO. Libraries were sequenced on Illumina NextSeq 2000 (2x51 bp for ChIP/ATAC, 2x150 bp for RNA-seq). Reads were mapped to hg38 by BWA (ChIP/ATAC) or STAR (RNA-seq), peaks called by MACS2 (FDR< 0.01), peak and gene raw counts estimated by featureCounts. Differential peaks (fold-change ≥1.5) and differentially expressed genes (fold-change ≥2) were identified by edgeR (FDR<0.05), including peaks detected in ≥2 samples and genes with TPM ≥0.1 in ≥1 sample. Interactions between promoters and distal regulatory regions were identified using publicly available chromatin interaction data in blood cell types. Pathway analysis was performed using Enrichr. Co-immunoprecipitation was performed in t(11;14) KMS12BM cells. Differential analysis comparing CCND1-high t(11;14) and CCND1-low MM cases in the two patient cohorts identified 21 upregulated and 21 downregulated genes associated with CCND1 expression which were also captured in the wild-typevs. CCND1 KO U266 cells. Enriched pathways (p<0.05) included Notch signaling, early/late estrogen response, epithelial-mesenchymal transition and reduced mitotic spindle and myogenesis. Dysregulated genes included those related to apoptosis, ubiquitin-proteosome, extracellular matrix, adhesion, and immune response, demonstrating a broader role for cyclin D1 in MM beyond cell cycle control. To explore the epigenomic basis of this transcriptional regulation, for each of the 42 differentially expressed genes we examined the changes of histone marks and chromatin accessibility in promoter and the interacting distal regulatory regions in the U266 KO model. Ten of the 21 genes upregulated in CCND1-high (downloaded in KO; e.g., argininosuccinate synthase 1 [ASS1], PMAIP1, MYEOV, CCR7, and DTX1) and 12 of the 21 genes downregulated in CCND1-high (upregulated in KO; e.g., UACA, RASAL2, PTPRZ1, PPT2, and IL23R) showed concordant changes in histone marks and/or chromatin accessibility, predominantly in distal regulatory regions, suggesting a chromatin-mediated mechanism of gene regulation by cyclin D1. The most prominent changes involved H3K27ac, a mark defining active enhancers. As H3K27ac is deposited by histone acetyltransferases CBP, we tested for the possible cyclin D1-CBP interaction. Co-IP confirmed this interaction, suggesting a role for this complex in chromatin remodeling and transcriptional control in t(11;14) MM. One gene notably upregulated in association with CCND1 is ASS1, which encodes a urea cycle enzyme responsible for endogenous arginine synthesis. Metabolomic studies demonstrated KO cells had elevated aspartate and citrulline and lowered arginine and ornithine compared to WT cells (p<0.01) demonstrating impaired urea cycle function. Further, KO cells displayed significantly reduced ASS1 protein expression, cell growth and S-phase entry compared to WT cells in arginine-replete conditions, with both WT and KO cell types further inhibited by arginine depletion (p<0.05). Together, a cyclin D1-associated gene signature was identified supporting its role in t(11;14) MM beyond the CDK4/6 mediated cell cycle control, including a role in metabolic plasticity and survival under arginine depleted conditions. These findings highlight cyclin D1 as a therapeutic vulnerability and open new avenues for personalized treatment strategies.

Abstract Motivation Efforts to address health disparities are often limited by the lack of robust computational tools for inferring genetic ancestry by calculating an individual’s genetic similarity to continental groups. We have already shown that a preferred alternative to self-described race is using ancestry-informative markers (AIMs) that can be classified into ancestral components and used to estimate their similarity to those of known populations to identify continental groups. However, real-world genomic data can present challenges, including limited availability of germline DNA, a small number of AIMs for each sample, and the use of different variant calling software, limiting the application of existing solutions. Results Here, we describe a novel supervised machine-learning tool AncestryGeni, which infers genetic ancestry for samples with even a hundred markers and is applicable to any genomic data, including whole exome sequencing (WES) and RNA sequencing (RNA-Seq) data. Applying AncestryGeni to a real-world genomic dataset obtained from the Multiple Myeloma Research Foundation (MMRF) CoMMpass study, we show that it is more accurate than the commonly used FastNGSadmix when using nonstandard genomic material. We also demonstrate that when using AncestryGeni, the tumor-derived sequence obtained from WES and RNA-Seq can be a robust data source to accurately estimate an individual’s genetic similarity to a continental group. Availability and implementation AncestryGeni pipeline is available at https://github.com/eelhaik/AncestryGeni/tree/main.

Abstract Clonal hematopoiesis (CH) of indeterminate potential (CHIP), driven by somatic mutations in leukemia-associated genes, confers increased risk of hematologic malignancies, cardiovascular disease, and all-cause mortality. In blood of healthy individuals, small CH clones can expand over time to reach 2% variant allele frequency (VAF), the current threshold for CHIP. Nevertheless, reliable detection of low-VAF CHIP mutations is challenging, often relying on deep targeted sequencing. Here, we present UNISOM, a streamlined workflow for enhancing CHIP detection from whole-genome and whole-exome sequencing data that are underpowered, especially for low VAFs. UNISOM utilizes a meta-caller for variant detection, in couple with machine learning models which classify variants into CHIP, germline, and artifact. In whole-exome sequencing data, UNISOM recovered nearly 80% of the CHIP mutations identified via deep targeted sequencing in the same cohort. Applied to whole-genome sequencing data from Mayo Clinic Biobank, it recapitulated the patterns previously established in much larger cohorts, including the most frequently mutated CHIP genes and predominant mutation types and signatures, as well as strong associations of CHIP with age and smoking status. Notably, 30% of the identified CHIP mutations had < 5% VAFs, demonstrating its high sensitivity toward small mutant clones. This workflow is applicable to CHIP screening in population genomic studies. The UNISOM pipeline is freely available at https://github.com/shulanmayo/UNISOM and https://ngdc.cncb.ac.cn/biocode/tool/7816.

Abstract Most disease-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWASs) reside in noncoding regions. There is increasing evidence that they play regulatory roles through epigenetic mechanisms by altering local chromatin state, 3D chromatin structure, and looping interactions. We and others have identified 42 CLL risk loci from GWASs. At 12 risk loci located outside of gene promoters, we previously demonstrated that functional SNPs are highly enriched in enhancers and super-enhancers that mostly target BCL2 family genes and transcriptional regulators. However, the study was limited by using only public epigenomic data generated in lymphoblastoid cell lines (LCL) and normal immune related cell types. To gain a better understanding of how the functional SNPs alter local chromatin state and modulate target gene expression in CLL, we generated RNA-seq data (n=588) and transposase-accessible chromatin with sequencing (ATAC-seq) data (n=69) in CLL tumor B cells. We performed molecular quantitative trait loci (QTL) analysis with index and correlated SNPs (R2 0.5) at 41 risk loci (excluding the HLA locus). The analyses identified cis expression QTLs (eQTLs), within 1Mb of genes’ transcription start site, at 25 risk loci and chromatin accessibility QTLs (caQTLs, within peak) at 11 risk loci, with both eQTLs and caQTLs at 8 risk loci. We identified a total of 54 eGenes, for which the expression is affected by nearby SNPs within 1Mb, including 40 novel eGenes at 16 CLL risk loci not previously reported. Notably, based on publicly available chromosome interaction data from CLL, LCL, and normal B cells, for over 70% of the novel eGenes, their promoters showed looping interactions with the regions carrying the eQTLs. Further, 13 SNPs were identified as both caQTL and eQTL, with the associated novel eGenes including CD81 at 11p15.5, IPCEF1 at 6q25.2, as well as CXXC1 and MBD1 at 18q21.1. The result suggests that these SNPs play a role in altering chromatin accessibility contributed to gene expression variation. At 11p15.5 risk locus, SNP rs2651823 (A/G) represents a novel caQTL and eQTL. The alternative G allele is associated with increased chromatin accessibility and CD81 expression. Hi-C data from GM12878 and CD20+ B cells supported interaction between the open chromatin region carrying rs2651823 and CD81 promoter. CD81 facilitates the trafficking of CD19 to the surface of activated B cells, and CD81-null mice showed impaired immune response. In a mouse model, CD81 contributes to tumor progression and metastasis in Lewis lung carcinoma. In summary, our molecular QTL analysis revealed that functional SNPs, located in CLL risk loci, modulate chromatin accessibility and expression of downstream target genes that could impact CLL risk. Citation Format: Huihuang Yan, Zhiquan Wang, Mingma Hoel, Cristine Allmer, Chantal McCabe, Daniel O’Brien, Dennis Robinson, James Cerhan, Sameer Parikh, Esteban Braggio, Susan Slager. Molecular QTL analyses identify functional variants at inherited risk loci in chronic lymphocytic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2763.

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.