Rothman, Nathaniel

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 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.

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.

AbstractNon‐Hodgkin lymphoma (NHL) is composed of a heterogeneous collection of subtypes with considerable differences in genetics, biology and aetiology. Studies to date on physical activity and NHL risk have not had sufficient sample size to evaluate whether associations differ by subtype. We pooled data from nine case‐control studies to examine the association between moderate‐to‐vigorous intensity physical activity (MVPA) and risk of NHL overall and by subtype (diffuse large B‐cell lymphoma, follicular lymphoma, chronic lymphocytic leukaemia/small lymphocytic lymphoma, marginal zone lymphoma and mature T‐cell lymphoma). A total of 5653 cases and 9115 controls were included in the pooled analysis. Physical activity was harmonised across nine studies and modelled as study‐specific tertiles. Multinomial logistic regression was used to estimate the association between physical activity and NHL, adjusting for confounders. The overall odds of NHL was 13% lower among participants in the most active tertile of MVPA compared to the least active tertile (adjusted odds ratio = 0.87, 95% CI = 0.80, 0.95). Similar decreases were observed across NHL subtypes. In summary, in this pooled analysis of case‐control studies, physical activity was associated with a modest risk reduction for each NHL subtype examined and with overall NHL.

Abstract Background: A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non–Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes. Methods: In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell–mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression. Results: We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08–1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59–2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell–mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction. Conclusions: Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions. Impact: Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.

Abstract Background: Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome. Methods: We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case–Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data. Results: Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98–1.02; arachidonic acid OR = 1.00, 95% CI = 0.99–1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87–1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex. Conclusions: Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk. Impact: These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer.

Abstract Background: Obesity and diabetes are major modifiable risk factors for pancreatic cancer. Interactions between genetic variants and diabetes/obesity have not previously been comprehensively investigated in pancreatic cancer at the genome-wide level. Methods: We conducted a gene–environment interaction (GxE) analysis including 8,255 cases and 11,900 controls from four pancreatic cancer genome-wide association study (GWAS) datasets (Pancreatic Cancer Cohort Consortium I–III and Pancreatic Cancer Case Control Consortium). Obesity (body mass index ≥30 kg/m2) and diabetes (duration ≥3 years) were the environmental variables of interest. Approximately 870,000 SNPs (minor allele frequency ≥0.005, genotyped in at least one dataset) were analyzed. Case–control (CC), case-only (CO), and joint-effect test methods were used for SNP-level GxE analysis. As a complementary approach, gene-based GxE analysis was also performed. Age, sex, study site, and principal components accounting for population substructure were included as covariates. Meta-analysis was applied to combine individual GWAS summary statistics. Results: No genome-wide significant interactions (departures from a log-additive odds model) with diabetes or obesity were detected at the SNP level by the CC or CO approaches. The joint-effect test detected numerous genome-wide significant GxE signals in the GWAS main effects top hit regions, but the significance diminished after adjusting for the GWAS top hits. In the gene-based analysis, a significant interaction of diabetes with variants in the FAM63A (family with sequence similarity 63 member A) gene (significance threshold P < 1.25 × 10−6) was observed in the meta-analysis (PGxE = 1.2 ×10−6, PJoint = 4.2 ×10−7). Conclusions: This analysis did not find significant GxE interactions at the SNP level but found one significant interaction with diabetes at the gene level. A larger sample size might unveil additional genetic factors via GxE scans. Impact: This study may contribute to discovering the mechanism of diabetes-associated pancreatic cancer.

Abstract Background. The strongest established prognostic factors for follicular lymphoma are age and clinical factors. Biologic studies suggest a major role for the host immunologic environment in follicular lymphomagenesis, which is in part determined by host genetic background. Cytokines are key regulators of immune function and regulation, are highly polymorphic, and have been implicated in lymphoma etiology and prognosis. We evaluated the hypothesis that inherited variability in cytokine and related immune genes impact overall survival in follicular lymphoma. Methods. We genotyped 73 single nucleotide polymorphisms (SNPs) from 44 candidate cytokine and immune genes in 278 follicular lymphoma patients who participated in a population-based case-control study conducted from 1998–2000 in the Detroit, Iowa, Los Angeles and Seattle SEER (Surveillance, Epidemiology and End Results) cancer registries. Baseline clinical data and survival through early 2005 were obtained from cancer registry files. Cox proportional hazards analysis was used to estimate hazard ratios (HR) and corresponding 95% confidence intervals for the association between individual SNPs and overall survival, adjusting for age, demographic and clinical factors. We also used parallel modeling strategies to identify the best summary multi-SNP risk score to predict survival, and applied these results in a time-dependent receiver-operator characteristics (ROC) analysis. Results. The median age at diagnosis was 57 years (range, 25–74), and 59 (21%) of the patients died during follow-up, with a median follow-up of 59 months (range, 27 – 78 months) for surviving patients. In multivariate modeling, SNPs in IL8 (rs4073; HRTT=2.14, 1.26–3.63), IL2 (rs2069762; HRGT/TT=1.80, 1.06–3.05), IL12B (rs3212227; HRAA/CC=1.83, 1.06–3.06), and IL1RN (rs454078; HRAA=1.93, 1.11–3.34) were the strongest and most robust predictors of survival. A summary score of the number of deleterious genotypes from these four genes was strongly associated with survival (p=1.4 × 10−5) after accounting for demographic and clinical variables (HR=2.06 per deleterious genotype, 1.52–2.79). A combination of the four SNPs and the demographic and clinical risk scores was strongly associated with survival (p=1.8 × 10−11); the 5-year Kaplan Meier survival estimates were 96% (93%–100%), 72% (62%–83%) and 58% (48%–72%) for low, intermediate, and high risk groups respectively. In a time-dependent ROC analysis the three risk groups had an area under the curve (AUC) of 0.83 at 72 months of follow-up (95% CI 0.72–0.93). Conclusion: Host genetic variability in immune genes, particularly IL8, IL2, IL12B, and IL1RN, individually and particularly in combination, was associated with overall survival in follicular lymphoma after accounting for demographic and clinical factors. Host immunogenetics is a promising class of prognostic factors in follicular lymphoma that warrants further evaluation.