Park, Hanla A.

Abstract Background: Mantle cell lymphoma (MCL) is an aggressive and rare B-cell non-Hodgkin lymphoma with poor prognosis, even with treatment. Despite its clinical importance, the etiology of MCL remains largely unknown. Some studies have reported a positive association with family history, suggesting that genetic factors could contribute to risk. Methods: To identify germline genetic variants associated with risk, we conducted the first genome-wide association studies (GWAS) and GWAS meta-analysis for MCL, comprising 1,163 cases and 61,271 controls. The most significant loci were taken forward for replication in three independent studies, including 576 cases and 771,773 controls. To identify potential target genes and pathways of discovered loci, we conducted F-MAGMA, colocalization, and transcriptome-wide association study (TWAS) analyses. To gain insight into underlying regulatory mechanisms, we performed integrative epigenomic analyses using DNase I hotspots and histone mark ChIP-seq data across diverse cell types applying the FORGE2 framework. Results: In the joint analyses, we identified eight genome-wide significant loci (P<5x10-8) associated with MCL risk, several of which were located near genes with known DNA repair (ATM), telomere (TERT), or RNA-binding functions (RBM20). Further analyses using F-MAGMA, which uses DNase-seq data to prioritize genes from GWAS, identified genes that were significantly associated with risk, including SP140, a chromatin reader and immune regulator of pathogen response. This gene also displayed strong colocalization with our lead variant in whole blood (COLOC posterior probability PP4=0.97). TWAS demonstrated significant associations with SP140 and ACTA2, which has been linked to cell proliferation (P<3x10-6). Enrichment for DNase I hotspots was observed in B and T lymphocytes (q-value<0.05), providing evidence for a role in immune cell-specific regulatory regions. Enrichment for enhancer-associated histone mark H3K4me1 was observed in B cells (q-value<0.01), highlighting a role for immune cell-specific enhancers. Transcription factor (TF) motif enrichment analysis pointed to TFs critical for B-cell development and differentiation, notably TCF3 and EP300 (q<0.01 for both), suggesting an association between B-cell lineage regulation and MCL risk. Finally, genome-wide SNP-based heritability was estimated to be 17%, underscoring the importance of common variants in MCL risk. Conclusion: Our study provides novel insight into the inherited susceptibility of MCL by identifying eight significant loci and estimating substantial common variant heritability. These findings highlight a role for several distinct biological pathways -specifically DNA repair, telomerase function, and B-cell differentiation- in the etiology of MCL, offering new targets for mechanistic investigation. Citation Format: Charles E. Breeze, Alyssa Clay-Gilmour, Hanla A. Park, Ólafur B. Davíðsson, Angelica Macauda, Brenda M. Birmann, Elizabeth E. Brown, Murat Güler, Michelle A. T. Hildebrandt, Alexandra Nieters, Krystle Ong, Jojo Biel-Nielsen Dietz, Karin E. Smedby, Karl Smith-Byrne, Rosalie Griffin, Sophia S. Wang, Susan Slager, James R. Cerhan, Jonathan N. Hofmann, Qing Lan, Nathaniel Rothman, Ingrid Glimelius, Henrik Hjalgrim, James Mckay, Sonja I. Berndt. Genome-wide association study of mantle cell lymphoma identifies novel loci suggesting a critical role for B-cell chromatin readers and DNA repair mechanisms [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 3602.

Abstract Introduction: With over 80, 000 new cases of non-Hodgkin's lymphoma (NHL) diagnosed annually, modifiable risk factors remain unclear. Circadian disruption, linked to cancers like breast and prostate, affects immune cells (e.g., natural killer cells and T-helper cells) by altering their trafficking and proliferation. Some studies (PMID: 30566672, 27611440) suggest long sleep may increase overall NHL risk, including a Mendelian randomization (MR) study (PMID: 32895918). However, no research has examined the causal relationship between sleep traits, including long sleep, and NHL risk, by subtypes. We hypothesize that sleep traits are causally associated with NHL subtypes in European populations, testing this using MR with data from the International Lymphoma Epidemiology Consortium (InterLymph). Methods: We conducted two-sample MR using genome-wide association studies (GWAS) summary data for sleep traits and NHL subtypes. Only independent, genome-wide significant (p < 5×10-8) Single Nucleotide Polymorphisms (SNPs) were selected as valid instruments, sourced from the UK Biobank. Sleep traits included chronotype (153 SNPs), insomnia (48 SNPs), sleep duration (78 SNPs), excessive daytime sleepiness (37 SNPs), short sleep (27 SNPs), and long sleep (8 SNPs). NHL subtype data were from InterLymph: follicular lymphoma (6, 508 cases / 64, 183 controls), mantle cell lymphoma (1, 169 cases / 61, 603 controls), Waldenstrom macroglobulinemia/lymphoplasmacytic lymphoma (WM/LPL: 1, 697 cases / 59, 333 controls), and chronic lymphocytic leukemia (8, 522 cases / 67, 653 controls). The primary analysis used inverse-variance weighted (IVW) random-effects, with MR-Egger for pleiotropy adjustment. Results: Using the IVW method, we did not identify any statistically significant causal association between sleep traits and the risk of NHL subtypes, but did yield several intriguing trends. Long sleep (sleep duration >= 9 h per night) showed a marginally increased NHL risk [Odds ratio (OR) range: 1.51 - 247.9] across each subtype. Short sleep (sleep duration < 7 h per night) showed a similar trend [OR range: 1.25 - 1.81], except for WM/LPL, where we noticed a negative trend. The confidence limit for long sleep was imprecise due to the small number of variants. We did not notice any significant directional horizontal pleiotropy or weak instrument bias. Conclusions: We found no conclusive evidence of a relationship between sleep traits and four NHL subtypes. However, we see suggestive trends of long and short sleep associated with higher risk of some NHL subtypes. These findings align with the previous studies, which showed long sleep was associated with higher risk of NHL overall. Future directions include expanding our analysis to additional NHL subtypes, generating polygenic risk scores, and conducting stratified analyses. Citation Format: Pankhil Shah, Brittany Crawford, Anwar Merchant, Brenda Birmann, Angelica Macauda, Michelle A. Hildebrandt, Aaron Norman, Neil E. Caporaso, Meredith Yeager, Michael Dean, Immaculata De Vivo, Lynn Goldin, Nicola J. Camp, Rosalie Griffin, Delphine Casabonne, Federico Canzian, Pelin Unal, Elad Ziv, Catherine R. Marinac, Alexandra Nieters, Stephen J. Chanock, Mitchell J. Machiela, Michael Conry, Charlie Zhong, Hanla A. Park, Simon Cheah, Jonathan N. Hofmann, Elizabeth E. Brown, Celine Vachon, Susan Slager, Sonja Berndt, Sophia S. Wang, Vijai Joseph, James McKay, Henrik Hjalgrim, Lara Sucheston-Campbell, Karl Smith-Byrne, Alyssa Clay-Gilmour. Association between sleep traits and risk of non-Hodgkin's lymphoma subtypes: a mendelian randomization study in the International Lymphoma Epidemiology (InterLymph) Consortium [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 7415.

Introduction: WM is a rare subtype of lymphoma, identified by the presence of IgM paraprotein, infiltration of clonal, small lymphoplasmacytic B cells in the bone marrow, and the somatic mutation, MYD88 L265P, which is observed in over 90% of WM patients. We previously identified two rare variants at 6p25.3 (EXOC2) and 14q32.13 (TCL1) that were significantly associated with an increased risk of WM/LPL (PMID: 30305637). However, the genetic basis for WM/LPL predisposition remains largely unknown. We conducted the largest WM/LPL GWAS to discover further germline susceptibility loci for WM/LPL. Methods: Our meta-analysis discovery set included 1,707 cases and 58,328 controls of European ancestry from the InterLymph Consortium, genotyped on a high-density SNP array and imputed to TOPMed, providing over 10 million common variants for association testing. Germline genetic loci with suggestive significance (P<5x10-6) were selected for replication in four independent sets (570 cases and 771,770 controls; UK Biobank, Copenhagen Hospital Biobank, FinnGen, Mayo Clinic), giving us a total of 2,277 WM/LPL cases and 831,098 controls. Conditional and sensitivity analyses were performed considering family history and sex. A transcriptome-wide association study (TWAS) was performed using FUSION with expression quantitative trait loci (eQTL) data from GTEx and Young Finns Study in whole blood, lymphocytes, and other tissues. Functional annotations were done using FUMA, F-MAGMA, FORGE, and OpenTarget. Colocalization was measured using CoLoC with whole blood gene expression data. Results & Discussion: We identified 11 germline genetic loci associated with WM/LPL risk that replicated and reached genome-wide significance (P< 5x10-8), including the two known loci on 6p25.3 and 14q32.13 and nine novel loci in the following regions: 6p25.3 (n=2 independent loci), 2q37.1, 9p13.2, 6p21.32, 10q23.31, 16q21, 1p31.1, and 15q15.1. Four of the discovered variants showed strong associations with WM/LPL [odds ratios (ORs) >10]. Enrichment was observed for enhancer and promoter elements in B-cells (FDR q < 0.01), and gene-based and colocalization analyses identified several potential biological pathways. Colocalization analyses show the intronic PAX5 risk variant on 9p (OR=1.5, P= 2.3x10-25) associated with whole blood eQTL expression of PAX5, which plays an essential role in regulating lymphoid progenitors to the B-lymphocyte lineage, as well as RP11-297B17.3 (RNA gene, antisense to PAX5). In addition, GWAS have shown this variant is in linkage disequilibrium (LD) with GWAS loci for CD20 on IgD+ CD38- B cell that represents a specific immunophenotype that helps identify and characterize B cells involved with WM and other similar conditions. The 2q37.1 (SP140) locus (OR=1.6, P=1.2x10-33) contains variants linked to chronic lymphocytic leukemia (CLL) and autoimmune diseases, including multiple sclerosis and Crohn's disease. This locus is important in chromatin-mediated gene regulation, suggesting a significant role in the immune system. The specific variant influences exon-skipping and potentially protein expression, altering the gene's function. The variants on 10q23.31 in FAS (OR=1.3, P=4.4x10-13) and 15q15.1 (OR=1.2, P=8.5x10-11) in BCL2 modifying factor (BMF), are in strong LD with known CLL susceptibility variants. The 15q variant is also linked with multiple lymphoid leukemia GW significant SNPs. Sex-specific analyses identified one novel locus at 20q13.33 associated with the expression of ZGPAT and UCKL1 in males (OR=1.9, P=9.5x10-10). TWAS identified associations with CRHR1 and CDKN2A; overexpression of both are shown to be associated with an increased risk of WM/LPL. Expression of CDKN2A has been linked to childhood acute lymphoblastic leukemia, while CRHR1 has been associated with breast cancer risk. Conclusions: Our study, the largest GWAS of WM to date, has identified 11 susceptibility loci, including two previously known and nine novel regions. The findings underscore the genetic complexity of WM, with the PAX5 and SP140 loci implicated in both chromatin-mediated gene regulation and immune system functions. These insights not only enhance our understanding of WM predisposition but also pave the way for future research focused on functional understanding of these loci, potentially leading to new therapeutic strategies.

Introduction: Follicular lymphoma (FL) is the most common indolent subtype of non-Hodgkin lymphoma. Family history is a risk factor. Genome-wide association studies (GWAS) have identified several germline variants, especially within the HLA region, to be associated with an increased risk of FL, but these variants only explain a fraction (16%) of heritability. Here, we present novel findings from GWAS of FL uncovering multiple germline genetic variants contributing to disease susceptibility. Methods: A GWAS was performed comprised of a discovery analysis in 6,508 FL cases and 64,183 controls of European ancestry followed by a replication of the most significant variants (P<5x10-7) in 2,745 FL cases and 782,226 controls from the biobanks of UK, Copenhagen Hospital, FinnGen and Mayo Clinic . At 9,253 FL cases and 846,409 FL free controls imputed genome-wide, this is the largest FL GWAS of European genetic ancestry to date. HLA alleles were imputed, and Class I and Class II signals were analysed at the amino acid level. A transcriptome wide association study (TWAS) was performed using FUSION investigate relationship to expression in GTEx and Young Finns Study in whole blood and cross-tissue . Functional annotations were done using FUMA, F-MAGMA, Open Targets and FORGE. Colocalization was conducted using CoLoC with whole blood gene expression data to understand causal variants. Mendelian randomization (MR) was performed on blood cell traits and leukocyte telomere lengths (LTL) to infer their relationships with FL genetic risk. Results: We replicated the previously published susceptibility loci within and outside the HLA region, finding two additional independent SNPs at two known loci (P < 5 x10-8) and discovered thirteen novel, genome-wide significant loci located at chromosome 2q13, 3p13, 3p24.1,3q13, 3q24.1, 4q31.21, 6p21, 6p22, 7q31, 10p14, 10q22, 10q23, 17q12, 19p13 and 20p11 with P < 5 x10-8. These SNPs were located near genes that show relevance for B-cell biology, such as BCL2L11, FOXP1, CD86, EOMES, IZKF3, PRF1, KLF2 NKX2-2 and HLA Class I and Class II genes. TWAS uncovered a role for H2AX, IGKV2-29, IMP, UBE4A and IKZF3 and provided directionality of effect for associated transcripts. Colocalization analysis corroborated a positive association with IKZF3 expression showing a robust correlation with FL risk variants (posterior probability PP4>0.7). In addition, several genome-wide significant variants in IKZF3 colocalize with eQTLs across tissues. IKZF3 (Aiolos) is part of the Ikaros family of transcription factors, critically important to the regulation of B-cell development and function. Therapeutic strategies leveraging Immunomodulatory imide drugs (IMiDs) to degrade IKZF3 have shown efficacy in multiple myeloma, and we speculate may also have potential benefits in FL by enhancing anti-tumor immunity and disrupting malignant B-cell function. Gene-set analysis showed blood and spleen to be the most relevant tissues for FL risk variants with lymphocyte activation and immune cell type differentiation genes demonstrating enrichment based on gene-ontology terms. Associated traits included allergies, respiratory diseases (e.g., asthma phenotypes), inflammatory bowel diseases, type 1 diabetes (T1D), lipid levels and commonly measured blood cell traits, with some SNPs showing evidence of colocalization with select blood traits. Interestingly, no causal association to blood traits was observed using MR; however, MR analyses did reveal a strong association with LTL. For the HLA loci, FL shows a diverse set of amino acid associations in HLA Class I and Class II genes that extends beyond what is known from autoimmune disease associations, such as rheumatoid arthritis or T1D. Our study also reveals interesting biological insights about FL susceptibility, such as a novel locus intragenic to FOXP1, an activated oncogene preferentially expressed in FL and DLBCL tumors, but rarely observed in Hodgkin's lymphoma. High levels of FOXP1 are known to result in poor survival in FL. Conclusion: This large GWAS of FL has identified multiple susceptibility variants, uncovering interesting gene candidates for further pathophysiological characterization of FL risk and B-cell development. It suggests new biological pathways and genes that could be potential leads for drug development and targeting treatment, in addition to identifying genetically high-risk populations for FL susceptibility.

AbstractWe conducted the first large genome‐wide association study to identify novel genetic variants that predict better (or poorer) prognosis in colorectal cancer patients receiving standard first‐line oxaliplatin‐based chemotherapy vs chemotherapy without oxaliplatin. We used data from two phase III trials, NCCTG N0147 and NCCTG N9741 and a population‐based patient cohort, DACHS. Multivariable Cox proportional hazards models were employed, including an interaction term between each SNP and type of treatment for overall survival (OS) and progression‐free survival. The analysis was performed for studies individually, and the results were combined using fixed‐effect meta‐analyses separately for resected stage III colon cancer (3098 patients from NCCTG N0147 and 549 patients from DACHS) and mCRC (505 patients from NCCTG N9741 and 437 patients from DACHS). We further performed gene‐based analysis as well as in silico bioinformatics analysis for CRC‐relevant functional genomic annotation of identified loci. In stage III colon cancer patients, a locus on chr22 (rs11912167) was associated with significantly poorer OS after oxaliplatin‐based chemotherapy vs chemotherapy without oxaliplatin (Pinteraction < 5 × 10−8). For mCRC patients, three loci on chr1 (rs1234556), chr12 (rs11052270) and chr15 (rs11858406) were found to be associated with differential OS (P < 5 × 10−7). The locus on chr1 located in the intronic region of RCSD1 was replicated in an independent cohort of 586 mCRC patients from ALGB/SWOG 80405 (Pinteraction = .04). The GWA gene‐based analysis yielded for RCSD1 the most significant association with differential OS in mCRC (P = 6.6 × 10−6). With further investigation into its biological mechanisms, this finding could potentially be used to individualize first‐line treatment and improve clinical outcomes.

Abstract Background: Associations between candidate germline genetic variants and treatment outcome of oxaliplatin, a drug commonly used for patients with colorectal cancer, have been reported but not robustly established. This study aimed to construct polygenic hazard scores (PHSs) as predictive markers for oxaliplatin treatment outcome by using a supervised principal component approach (PCA). Methods: Genome-wide association analysis for overall survival, including interaction terms (SNP*treatment type) was carried out using two phase III trials, 3,098 resected stage III colon cancer (rCC) patients of NCCTG N0147 and 506 metastatic colorectal cancer (mCRC) patients of NCCTG N9741, separately. SNPs showing interaction with genome-wide significance (P &lt; 5 × 10–8) were selected for PCA to derive a PHS. PHS interaction with treatment was included in Cox regression models to predict outcome. Replication of prediction models was performed in an independent cohort, DACHS. Results: The two PHSs based on the first two principal components of selected SNPs (15SNPs for rCC and 13SNPs for mCRC) were used to construct interaction terms with treatment type and included in models adjusted for clinical covariables. However, in the DACHS study, the addition of the two PHS terms to clinical models did not improve the prediction error in either patients with rCC or mCRC. PHS interaction was also not replicated. Conclusions: The PHSs derived using principal components efficiently combined multiple predictive SNPs for estimating likelihood of benefit from oxaliplatin versus other treatment but could not be replicated. Impact: These results highlight the potential but also challenges in generating evidence for a predictive polygenic score for oxaliplatin efficacy.