Rothman, Nathaniel

Abstract Pancreatic cancer disproportionately affects Black individuals in the United States, but they have limited representation in genetic studies of pancreatic ductal adenocarcinoma (PDAC). To address this gap, we performed admixture mapping and genome-wide association analysis (GWAS) in genetically inferred African ancestry individuals (1,030 cases and 889 controls). Admixture mapping identified three regions with a significantly higher proportion of African ancestry in cases compared to controls (5q33.3, 10p1, 22q12.3). GWAS identified a genome-wide significant association at 5p15.33 ( CLPTM1L , rs383009:T>C, T Allele Frequency=0.51, OR:1.45, P value=1.24×10 -8 ), a locus previously associated with PDAC. Known loci at 5p15.33, 7q32.3, 8q24.21 and 7q25.1 also replicated (P value <0.01). Multi-ancestral fine-mapping identified two potential causal SNPs (rs3830069 and rs2735940) at 5p15.33. Collectively these findings identified novel PDAC risk loci and expanded our understanding of this deadly cancer in underrepresented populations, emphasizing the multifactorial nature of PDAC risk including inherited genetic and non-genetic factors. Statement of Significance To understand how genetic variation contributes to PDAC risk in Black people in North American, we studied individuals of genetically-inferred African ancestry. We identified novel risk loci and differences in the contribution of known loci. This demonstrates that ancestry-informed genetic analyses improve our understanding of PDAC risk and enhances discovery.

Abstract Background: Epidemiologic and other studies have identified environmental, occupational, and genetic risk factors for lymphoid and myeloid malignancies, but most studies have been conducted in Western populations. Investigations in populations with differing exposure patterns, distributions of disease subtypes, and genetic architecture are needed to fully understand hematopoietic tumor etiology. Methods: We conducted a large, multicenter, hospital-based case-control study of lymphoid and myeloid neoplasms in East Asia (AsiaLymph), including 5,671 lymphoid cases, 1,879 myeloid cases, and 3,858 controls. Participants completed a computer-assisted personal interview and provided biospecimens. Cases underwent central pathology review and were coded into the WHO Classification. Herein, we describe the study methods in detail and report association results for education, body mass index, and family history. Results: Greater BMI at age 20 but not at age 40 was associated with odds of total lymphoid neoplasm (per 5 kg/m2 increase: OR [95% CI]: 1.18 [1.09-1.27]), total myeloid neoplasm (OR [95% CI]: 1.17 [1.05-1.29]), and several subtypes. Greater educational attainment was associated with increased odds of total lymphoid neoplasm (for college vs. less than primary education: OR [95% CI]: 1.41 [1.21-1.65]) but not myeloid neoplasm (OR [95% CI]: 1.17 [0.95-1.45]; p-heterogeneity=0.02). There were also positive associations between family history of hematologic cancer in first degree relatives and odds of lymphoid neoplasm (OR [95% CI]: 1.53 [1.15-2.03]) and myeloid neoplasm (OR [95% CI]: 1.65 [1.11-2.44]) and specific subtypes. None of the evaluated risk factors were associated with NK/T-cell lymphoma, which supports a distinct etiology for this subtype. Conclusion: The AsiaLymph Study is one of the largest molecular epidemiology studies of both lymphoid and myeloid neoplasms with standardized World Health Organization classification of histopathologic subtypes and will serve as a valuable resource for etiologic investigations into risk factors for these malignancies. Citation Format: Qing Lan, Lauren M. Hurwitz, John K. Chan, Tai Hing Lam, Kexin Chen, Yok Lam Kwong, Xu Caigang, Brian CH Chiu, Raymond Liang, Ip Dennis, Wei Hu, Bryan Bassig, Mark Purdue, Jun Xu, Sarah Locke, Sophia S. Wang, James R. Cerhan, Sonja Berndt, Jonathan N. Hofmann, Jianxin Shi, Kai Yu, Shahinaz Gadalla, Lisa J. McReynolds, Rena Jones, Hongji Dai, Zhangyan Lyu, Lugui Qiu, Wei Liu, Huilai Zhang, Xianhuo Wang, Lindsay M. Morton, Stephen Chanock, Martha Linet, Melissa C. Friesen, Roel Vermeulen, Nathaniel Rothman. A multi-center hospital-based case-control study of lymphoid and myeloid neoplasms (AsiaLymph): Study design and initial findings for education, body mass index, and family history [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 5051.

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&lt;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&lt;3x10-6). Enrichment for DNase I hotspots was observed in B and T lymphocytes (q-value&lt;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&lt;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&lt;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.

Summary Evidence for an association between insulin‐like growth factors (IGF) and multiple myeloma (MM) is inconsistent. We examined total IGF‐I concentrations and risk of MM by combining baseline serological data among UK Biobank participants ( n = 444 187; 732 incident MM) with a two‐sample Mendelian randomisation (MR) analysis using identified genetic variants associated with circulating total IGF‐I and IGF‐binding protein 3 (IGFBP‐3) in the InterLymph consortium (2434 MM and their 2567 controls). Finally, additional lymphoid neoplasm (LN) subtypes were included for comparison with the main hypothesis. Circulating IGF‐I level was positively associated with MM risk Hazard ratio‐HR‐per one standard deviation‐SD‐increase (HR 1‐SD = 1.11, 95% confidence interval [CI]: 1.01–1.22; p ‐value = 0.03), especially closer to diagnosis. Genetically inferred IGF‐I levels were associated with increased MM risk (odds ratio [OR] = 1.27, 95% CI: 1.05–1.54) but not with any other LNs. Genetically inferred IGFBP‐3 levels showed no associations with any LN evaluated. Corroborating previous findings, in a secondary analysis, IGF‐I levels were associated with the risk of chronic lymphocytic leukaemia/small lymphocytic lymphoma (CLL/SLL) in males with higher body mass index (HR 1‐SD in obese male = 1.36, 95% CI: 1.14–1.61). Our serological and MR analyses suggest a contributing role of IGF‐I in the susceptibility of MM; the CLL/SLL findings warrant further investigation considering sex‐specific adiposity.

Abstract Chromosome 5p15.33 harbors several independent association signals which demonstrate antagonistic pleiotropy across cancer types, with causal mechanisms largely unresolved. To identify functional variants and enhancer elements at this locus, we performed statistical fine-mapping followed by massively parallel reporter assays (MPRA) and proliferation based CRISPRi screens. This approach identified eight multi-cancer functional variants (MCFVs) across three GWAS signals. Targeting rs421629 (part of the CLPTM1L signal marked by rs465498) with CRISPRi revealed opposing effects on TERT expression in pancreatic versus lung cancer cells, consistent with the antagonistic pleiotropy observed for this signal. Furthermore, CRISPRi nominated an intronic CLPTM1L variable number tandem repeat (VNTR) as a potent enhancer. Long-read sequencing established VNTR polymorphisms as potential causal variants for the rs465498 signal. We showed that Hippo-pathway transcription factors mediate VNTR enhancer activity in lung and pancreatic cancer cells. Together, these findings indicate that cancer susceptibility at 5p15.33 may be mediated by both SNPs and VNTRs and provide an integrated framework for resolving complex pleiotropic loci.

Abstract Background The associations between different types of diabetes, characterized by distinct pathophysiology and genetic architecture, and pancreatic ductal adenocarcinoma (PDAC) risk are not understood. Methods We investigated associations of genetic susceptibility to type 2 diabetes (T2D), 8 T2D mechanistic clusters, type 1 diabetes (T1D), and maturity-onset diabetes of the young (MODY) with PDAC risk. We used genome-wide association study (GWAS) summary-level statistics for T2D (242 283 cases, 1 569 734 controls), T1D (18 942 cases, 501 638 controls), and PDAC (10 244 cases and 360 535 controls) in individuals of European ancestry. Results Two-sample Mendelian randomization (MR) using the Robust Adjusted Profile Score (MR-RAPS) method indicated that genetically predicted T2D was associated with PDAC risk (OR = 1.10; 95% CI = 1.05 to 1.15), particularly the T2D obesity (OR = 1.28; 95% CI = 1.15 to 1.42) and lipodystrophy (OR = 1.25; 95% CI = 1.03 to 1.51) clusters. No association was observed for T1D with PDAC risk (OR = 1.01; 95% CI = 0.99 to 1.02). Pathway/gene-set analysis using the summary-based Adaptive Rank Truncated Product (sARTP) method revealed a significant association between the MODY gene-sets and PDAC risk (P = 1.5 × 10−8), which remained after excluding 20 known PDAC GWAS loci (P = 7.6 × 10−4). HNF1A, FOXA3, and HNF4A were the top contributing genes after excluding the previously identified GWAS loci regions. Conclusions Our results from this genetic association study support that T2D, particularly the obesity and lipodystrophy mechanistic clusters, and MODY genomic susceptibility regions play a role in the etiology of PDAC.

Abstract Understanding tumor cell dynamics can improve prognosis and treatment but remains limited for lung adenocarcinoma in people who have never smoked (NS-LUAD). With RNA sequencing data from 684 NS-LUAD cases and validation in an independent dataset, we identified three subtypes with distinct phenotypic traits and cell compositions. Additional genomic and histologic data further characterized the subtypes. “Steady,” marked by low proliferation, high alveolar cell fraction, moderate-to-well differentiation, and fewer driver gene alterations, is linked to prolonged survival and low immune evasion. “Proliferative” shows high proliferation markers, TP53 mutations, and gene fusions. “Chaotic,” with high epithelial-to-mesenchymal transition markers, has the worst prognosis, even within stage I tumors. Lacking known molecular or histologic characteristics, this aggressive subtype is solely identified by transcriptomic data. A 60-gene signature recapitulates the classification and predicts survival even within subgroups based on tumor stage or known genomic features, emphasizing its potential for improving early-stage NS-LUAD prognostication in clinical settings. Significance: The transcriptome of 684 NS-LUAD identifies three subtypes with different cellular dynamics and genomic and morphologic features. A 60-gene signature accurately stratifies subjects for mortality risk, even in stage I, offering a potential clinically applicable tool for treatment decision-making in patients with NS-LUAD. See related commentary by Azizi et al., p. 423

Abstract Background Lung adenocarcinoma in never-smokers is a major public health burden, especially among East Asian women. Polygenic risk scores (PRSs) are promising for risk stratification but are primarily developed in European-ancestry populations. We aimed to develop and validate single- and multi-ancestry PRSs for East Asian never-smokers to improve lung adenocarcinoma risk prediction. Methods PRSs were developed using genome-wide association study summary statistics from East Asian (8002 cases; 20 782 controls) and European (2058 cases; 5575 controls) populations. Single-ancestry models included a PRS with 25 single nucleotide variants (PRS-25), PRS clumping and thresholding (CT), and LDpred2; multi-ancestry models included LDpred2 plus European PRS-128, PRS-CSx, and clumping and thresholding, super learning, empirical Bayes (CT-SLEB). Performance was evaluated in independent East Asian data from the Female Lung Cancer Consortium and externally validated in the Nanjing Lung Cancer Cohort. We assessed predictive accuracy via area under the curve operating procedure (AUC), with 10-year and (ages 30-80 years) absolute risks estimates. Results The best multi-ancestry PRS, using East Asian and European data via CT-SLEB, outperformed the best East Asian-only (LDpred2; AUC = 0.629, 95% confidence interval [CI] = 0.618 to 0.641), achieving an AUC of 0.640 (95% CI = 0.629 to 0.653) and odds ratio of 1.71 (95% CI = 1.61 to 1.82) per standard deviation increase. Nanjing Lung Cancer Cohort validation confirmed robust performance (AUC = 0.649, 95% CI = 0.623 to 0.676). The top 20% PRS group had a 3.92-fold higher lung adenocarcinoma risk than the bottom 20%. Further, the top 5% PRS group reached a 6.69% lifetime absolute risk. Notably, this group reached the average population 10-year lung adenocarcinoma risk at age 50 years (0.42%) by age 41 years, 9 years earlier. Conclusions Multi-ancestry PRS approaches enhance lung adenocarcinoma risk stratification in East Asian never-smokers, with consistent external validation, suggesting future clinical utility.

Abstract Obesity has been associated with non‐Hodgkin lymphoma (NHL), but the evidence is inconclusive. We examined the association between genetically determined adiposity and four common NHL subtypes: diffuse large B‐cell lymphoma (DLBCL), follicular lymphoma, chronic lymphocytic leukemia, and marginal zone lymphoma, using eight genome‐wide association studies of European ancestry ( N = 10,629 cases, 9505 controls) and constructing polygenic scores for body mass index (BMI), waist‐to‐hip ratio (WHR), and waist‐to‐hip ratio adjusted for BMI (WHRadjBMI). Higher genetically determined BMI was associated with an increased risk of DLBCL [odds ratio (OR) per standard deviation (SD) = 1.18, 95% confidence interval (95% CI): 1.05–1.33, p = .005]. This finding was consistent with Mendelian randomization analyses, which demonstrated a similar increased risk of DLBCL with higher genetically determined BMI (OR per SD = 1.12, 95% CI: 1.02–1.23, p = .03). No significant associations were observed with other NHL subtypes. Our study demonstrates a positive link between a genetically determined BMI and an increased risk of DLBCL, providing additional support for increased adiposity as a risk factor for DLBCL.

ABSTRACT Understanding lung cancer evolution can identify tools for intercepting its growth. In a landscape analysis of 1024 lung adenocarcinomas (LUAD) with deep whole-genome sequencing integrated with multiomic data, we identified 542 LUAD that displayed diverse clonal architecture. In this group, we observed an interplay between mobile elements, endogenous and exogenous mutational processes, distinct driver genes, and epidemiological features. Our results revealed divergent evolutionary trajectories based on tobacco smoking exposure, ancestry, and sex. LUAD from smokers showed an abundance of tobacco-related C:G>A:T driver mutations in KRAS plus short subclonal diversification. LUAD in never smokers showed early occurrence of copy number alterations and EGFR mutations associated with SBS5 and SBS40a mutational signatures. Tumors harboring EGFR mutations exhibited long latency, particularly in females of European-ancestry (EU_N). In EU_N, EGFR mutations preceded the occurrence of other driver genes, including TP53 and RBM10. Tumors from Asian never smokers showed a short clonal evolution and presented with heterogeneous repetitive patterns for the inferred mutational order. Importantly, we found that the mutational signature ID2 is a marker of a previously unrecognized mechanism for LUAD evolution. Tumors with ID2 showed short latency and high L1 retrotransposon activity linked to L1 promoter demethylation. These tumors exhibited an aggressive phenotype, characterized by increased genomic instability, elevated hypoxia scores, low burden of neoantigens, propensity to develop metastasis, and poor overall survival. Re-activated L1 retrotransposition-induced mutagenesis can contribute to the origin of the mutational signature ID2, including through the regulation of the transcriptional factor ZNF695 , a member of the KZFP family. The complex nature of LUAD evolution creates both challenges and opportunities for screening and treatment plans.