Citations

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

Abstract Multiple myeloma (MM) is the most common blood cancer in African American (AA) individuals. AAs have an increased risk of developing Monoclonal Gammopathy of Undetermined Significance (MGUS) and MM compared to non-Hispanic White (NHW) individuals. While several studies have shown a pivotal role for the immune tumor microenvironment (iTME) in MM pathogenesis, the composition of the iTME in AAs with MM is under investigated. Viably frozen whole bone marrow (BM) aspirates from 128 individuals (53 AA and 75 NHW) with either MM (n=65), MGUS (n=28), smoldering MM (SMM, n=13), or healthy controls (n=22) were evaluated by single cell RNA-sequencing (scRNA-seq) or mass cytometry (CyTOF). We characterized 34 (14 AA and 20 NHW) by scRNA-seq and 94 (39 AA and 55 NHW) by CyTOF. CyTOF samples were stained with Maxpar Direct Immune Profiling Assay Kit (Fluidigm), acquired on a Helios mass cytometer and analyzed in Cytobank. Data were analyzed by logistic equivalent regression model with p<0.05 for significance. ScRNA-seq samples were processed as in Yao, 2022. Data were normalized using scvi-tools and clustered by Uniform manifold approximation and projection (UMAP). Differentially expressed markers between clusters were identified using Model-based Analysis of Single-cell Transcriptomics (MAST) in Seurat R package at FDR < 0.05 and log2 fold change > 0.55. ScRNA-seq of the CD45+ non-tumor lymphocytes revealed 17 immune populations (HSC, erythroblast, MDSC, GMP, CD1c+ DC, CD4+ and CD8+ T cell, CD14+ M1 macrophage, CD16+monocyte, monocyte, pDC, naïve B, pre-B, pro-B, memory B, NK/MAIT). CD3E/D+cells were clustered into naïve/central memory (CCR7+, LEF1+, SELL+), memory (GZMK+), effector (GZMB+, GZMH+, FGFBP2+), g/d (TRDC+), MAIT (SLC4A10+, KLRB1+) and IFN CD8+ T cells (IFIT1+, IFIT2+,IFIT3+). Disease progression was associated with increases in median CD4+ and CD8+ GZMK+ T cells, memory B cells and a reduction in all monocyte clusters. Within the diseased cohort in comparison to NHW patients, AA patients had increased CD3+ T cells (CD4 and CD8+) (AA 57.5% vs. NHW 44.3%) and a reduction in overall monocytes (AA 10.5% vs. 26.5%). These results were validated by CyTOF, which confirmed an increase in CD3+ T cells (AA 53.7% vs. NHW 35.7%, p<0.001) including increases in CD8+ (AA 7.0% vs. NHW 3.2%, p<0.05) and CD4+ terminal effector T cells (AA 5.7% vs. NHW 2.7%, p<0.05) and a reduction in monocytes (AA 18.3% vs. NHW 26.6%) among AAs. After adjusting for age and sex, AAs retained an increase in T cells (OR 1.59 (95%CI 1.2-2.01, p=0.00023)) and a reduction in monocytes (OR 0.88 (95%CI 0.78-1, p=0.045)). Despite the increase in CD8+ T cells in AA patients, we identified reduced GZMK+ T cells and reductions in exhaustion markers (TIGIT, LAG3, TOX) in CD8+ T cells in comparison to NHWs. Nevertheless, the distribution of the immune populations was similar between healthy AA and NHW individuals. Our multiomics approach to characterize the iTME identified significant differences among AA compared to NHW patients with MM. Citation Format: Suganti Shivaram, Hongwei Tang, Huihuang Yan, Shulan Tian, Guyu Qin, Emilie Anderson, Neeraj Sharma, Abiola Bolarinwa, Cinthya Zepeda-Mendoza, Stacey Lehman, Felicia Gomez, Amit Mitra, Surendra Dasari, Taxiarchis Kourelis, Shaji Kumar, Linda Baughn. Characterization of the Immune Tumor Microenvironment in a Diverse Cohort of Patients with Multiple Myeloma [abstract]. In: Proceedings of the Blood Cancer Discovery Symposium; 2024 Mar 4-6; Boston, MA. Philadelphia (PA): AACR; Blood Cancer Discov 2024;5(2_Suppl):Abstract nr P06.

Abstract Purpose: Multiple myeloma (MM) is a plasma cell malignancy characterized by a collection of cytogenetic abnormalities that drive clinical presentation, response to treatment, and prognosis. The most common immunoglobulin heavy chain (IGH) translocation found in MM, t(11;14), results in CCND1 overexpression. Although t(11;14) is considered a standard risk cytogenetic abnormality, heterogeneity in patient outcomes and response to therapy within the t(11;14) group has been reported. In contrast to t(11;14) mantle cell lymphoma, where most cases (73%) have a typical, balanced t(11;14) rearrangement, only 47% of t(11;14) MM have a typical, balanced t(11;14) FISH pattern; the remaining 53% of MM are characterized by an atypical t(11;14) FISH result (Dalland, Genes Chromosomes Cancer, 2021). Here, we compared the genetic characteristics of typical vs. atypical t(11;14) MM and the association with other genomic abnormalities, CCND1 expression, IGH breakpoint location and patient outcomes. Methods: We performed a retrospective study of newly diagnosed t(11;14) MM patients seen at Mayo Clinic (N=411) or enrolled in the MMRF CoMMpass trial (N=167). The t(11;14) FISH results of the Mayo Clinic samples were segmented into two groups: typical CCND1/IGH signal pattern (n=179, 3 red, 3 green, 2 fusion signals) and atypical CCND1/IGH signal pattern (n=232, any FISH pattern deviating from typical). The CoMMpass cases had both RNAseq and long-insert WGS data and were analyzed to detect gene expression, structural and copy number variants and IGH breakpoint locations. These data were used to classify the t(11;14) as typical (n=63) or atypical (n=104) and whether there was a single t(11;14) fusion (n=107) or a gain of the t(11;14) fusion (n=60), and whether these cases were associated with changes in gene expression or breakpoint location. Differences between the groups were measured using χ2 and Wilcoxon tests. Survival curves were estimated using Kaplan Meier and compared using the Log-Rank test. Statistical significance was determined at p<0.05. Results: Of the Mayo t(11;14) cohort, 232/411 (56%) had an atypical t(11;14) FISH pattern and 179/411 (44%) had a typical t(11;14) FISH pattern. Although an atypical t(11;14) FISH pattern was not significantly associated with monosomy 17/17p deletion, MYC rearrangements, 1q gain or ISS stage III, it was significantly associated with the presence of any trisomy in comparison to typical t(11;14) cases (40/222=18% vs.12/166=7%, p=0.002). Atypical t(11;14) cases had an increase in CCND1 gene expression compared to typical t(11;14) cases (median 375.0 RPKM vs. 266.4 RPKM, p<0.0001) and cases with a gain of the t(11;14) fusion were associated with even greater CCND1 expression compared to cases with a single t(11;14) fusion (median 459.4 RPKM vs. 263.2 RPKM, p<0.0001). We next evaluated whether the t(11:14) FISH pattern was associated with differences in IGH breakpoint locations. Nearly 80.0% of atypical t(11;14) cases had a breakpoint within the constant region of the IGH locus commonly reported for MM, while only 61.0% of the typical cases had constant region breakpoints (p=0.035); the remaining cases had breakpoints involving the VDJ region. We next evaluated the impact of atypical t(11;14) FISH on OS. In Mayo cases, OS was significantly shorter in patients with atypical FISH patterns compared to typical t(11;14) FISH patterns (6.9 years [95% CI: 5.1-7.8] vs. 9.0 years [95% CI: 7.0-11.6], p=0.019). Although atypical t(11;14) was associated with increased risk of death (RR 1.45, p=0.021), this finding was not retained in a multivariate model including high risk abnormalities such as age, ISS stage, 17p deletion and 1q gain demonstrating atypical t(11;14) is not an independent predictor of poor outcome when treated with standard therapies. Conclusions: Atypical t(11;14) FISH in MM is found in approximately half of t(11;14) cases, is associated with trisomies, higher CCND1 expression and is more likely to have a breakpoint involving the constant region of the IGH locus. The enrichment of VDJ breakpoints within typical t(11;14) MM is consistent with enrichment of VDJ IGH breakpoints in MCL, which are mostly characterized by typical t(11;14) FISH patterns. While not an independent predictor of outcome with standard therapies, it remains to be evaluated whether patients with typical or atypical t(11;14) have a differential response to novel agents, such as venetoclax. Disclosures Fonseca: Janssen: Consultancy; Novartis: Consultancy; Pharmacyclics: Consultancy; Sanofi: Consultancy; Merck: Consultancy; Juno: Consultancy; Caris Life Sciences: Membership on an entity's Board of Directors or advisory committees; Kite: Consultancy; Aduro: Consultancy; OncoTracker: Consultancy, Membership on an entity's Board of Directors or advisory committees; GSK: Consultancy; AbbVie: Consultancy; Patent: Prognosticaton of myeloma via FISH: Patents & Royalties; Mayo Clinic in Arizona: Current Employment; Amgen: Consultancy; BMS: Consultancy; Celgene: Consultancy; Takeda: Consultancy; Bayer: Consultancy; Scientific Advisory Board: Adaptive Biotechnologies: Membership on an entity's Board of Directors or advisory committees. Bergsagel: Pfizer: Consultancy, Honoraria; Celgene: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; GSK: Consultancy, Honoraria; Oncopeptides: Consultancy, Honoraria; Genetech: Consultancy, Honoraria; Novartis: Consultancy, Honoraria, Patents & Royalties: human CRBN mouse. Kumar: Merck: Research Funding; Carsgen: Research Funding; Roche-Genentech: Consultancy, Research Funding; Beigene: Consultancy; KITE: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Bluebird Bio: Consultancy; Antengene: Consultancy, Honoraria; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Consultancy; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Astra-Zeneca: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Research Funding; Adaptive: Membership on an entity's Board of Directors or advisory committees, Research Funding; Tenebio: Research Funding; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Amgen: Consultancy, Research Funding; Abbvie: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; BMS: Consultancy, Research Funding; Sanofi: Research Funding.

Abstract Purpose: Multiple myeloma (MM) is a plasma cell (PC) malignancy with an increasing incidence in the United States. Fluorescence in situ hybridization (FISH) is currently the gold-standard diagnostic assay to detect recurrent genomic abnormalities of prognostic and therapeutic significance in MM. Hyperdiploid MM, a standard risk abnormality, is characterized by gains of odd-numbered chromosomes including chromosome 15. While gains of chromosome 15 have been reported in approximately 53% of standard risk MM patients (Chretien, Blood, 2015), losses of chromosome 15 are rarely reported. Enumeration of chromosome 15 centromere is detected using the alpha satellite FISH probe. We previously discovered ~9% of cases had false-negative chromosome 15 FISH results when compared to Mate Pair Sequencing (MPseq) (Smadbeck, BCJ, 2019). To our knowledge, no previous studies have reported a reduced performance of the alpha satellite D15Z4 FISH probe in PC FISH studies. A false-negative result for D15Z4 could result in an underreporting of trisomy 15 and overreporting of monosomy 15. Here, we investigated the incidence of false-negative D15Z4 FISH results by examining the presence of monosomy 15 from patients with a monoclonal gammopathy. Methods: We identified 1231 samples from patients with a monoclonal gammopathy who had undergone uniform testing to identify MM-specific cytogenetic abnormalities. DNA from bone marrow samples was genotyped on the Precision Medicine Research Array and biogeographical ancestry was assessed using the Geographic Population Structure Origins tool (Elhaik, Nat Commun, 2014; Home DNA 2016). FISH of immunoglobulin (cIg)-stained PCs was performed as described (Baughn, BCJ, 2018). Chromosome 15 centromere probe (Vysis CEP 15 D15Z4, Abbott Molecular, Abbott Park, IL) was used for chromosome 15 enumeration. Samples were divided into two groups: Cases with or without evidence of monosomy 15 by FISH. Mann-Whitney U test was applied to compare difference between groups for data deviating from normal distribution. A Chi-squared test evaluating the differences across these two groups was used. Results: We identified monosomy 15 in 31 cases (2.5% of cohort). Monosomy 15 was observed in nearly 100% of both the PC (cIg+) and non-PC (cIg-) populations, by contrast to the other recurrent PC-specific genomic abnormalities which were restricted to only the PC population. All samples with a banded analysis revealed most metaphases had two normal chromosome 15s demonstrating that the D15Z4 FISH probe failed to correctly enumerate the number of chromosome 15 centromere signals. Of the samples with a monosomy 15 result, the median African ancestry was 76.2% (95% percentiles: 77.5%-89.0%), while the median African ancestry of the non-monosomy 15 cohort was 2.3% (95% percentiles: 0%-87.7%) (). Although not all patients with >50% African ancestry had evidence of monosomy 15, 26/31 (83.9%) of patients with monosomy 15 had >50% African ancestry. Accordingly, patients with African ancestry (similarity ≥ 0.70) had 8.02-fold (95% confidence interval: 3.73-17.25, P = 9.92 X10 -8) increased possibility to have monosomy 15. Of cases with monosomy 15, 22 (71.0%) had an IGH rearrangement with the most prevalent IGH partner being CCND1. The incidence of t(11;14) was greater in the monosomy 15 group (13/31:41.9%), in contrast to the non-monosomy 15 group (330/1200=27.5%) (P = 0.12). Conclusions: We report a polymorphic variant of the chromosome 15 centromere resulting in a false-negative FISH result using the D15Z4 FISH probe. Surprisingly, the false-negative FISH result was not uniformly distributed among patients and was highly enriched in individuals of African ancestry. This may be due to a variation of the chromosome 15 alpha satellite in some individuals that prevents the D15Z4 FISH probe from sufficiently binding (O'Keefe, Genome Research, 2000). Although this study focused on identification of monosomy 15, the incidence of a false-negative result is likely greater than we have reported in this study since we did not evaluate cases that had a true trisomy 15 that was reported as normal chromosome 15 enumeration by FISH. This study highlights a critical need to ensure diagnostic reagents have equal performance for all patients independent of their ancestry. Evaluation of genomic events using unbiased approaches such as whole genome sequencing may circumvent some of these limitations. Disclosures Elhaik: DDC: Consultancy. Baird: DDC: Current Employment. Kumar: Antengene: Consultancy, Honoraria; Bluebird Bio: Consultancy; Beigene: Consultancy; BMS: Consultancy, Research Funding; Tenebio: Research Funding; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Abbvie: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Amgen: Consultancy, Research Funding; Novartis: Research Funding; Astra-Zeneca: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Merck: Research Funding; KITE: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Roche-Genentech: Consultancy, Research Funding; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Oncopeptides: Consultancy; Carsgen: Research Funding; Celgene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Adaptive: Membership on an entity's Board of Directors or advisory committees, Research Funding; Sanofi: Research Funding.

AbstractIt is of great scientific interest to identify interactions between genetic variants and environmental exposures that may modify the risk of complex diseases. However, larger sample sizes are usually required to detect gene‐by‐environment interaction (G × E) than required to detect genetic main association effects. To boost the statistical power and improve the understanding of the underlying molecular mechanisms, we incorporate functional genomics information, specifically, expression quantitative trait loci (eQTLs), into a data‐adaptive G × E test, called aGEw. This test adaptively chooses the best eQTL weights from multiple tissues and provides an extra layer of weighting at the genetic variant level. Extensive simulations show that the aGEw test can control the Type 1 error rate, and the power is resilient to the inclusion of neutral variants and noninformative external weights. We applied the proposed aGEw test to the Pancreatic Cancer Case–Control Consortium (discovery cohort of 3,585 cases and 3,482 controls) and the PanScan II genome‐wide association study data (replication cohort of 2,021 cases and 2,105 controls) with smoking as the exposure of interest. Two novel putative smoking‐related pancreatic cancer susceptibility genes, TRIP10 and KDM3A, were identified. The aGEw test is implemented in an R package aGE.

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: Chronic inflammation is implicated in pancreatic cancer carcinogenesis. Advanced glycation end-products (AGEs), which are rich in cooked red meat and cigarette smoke, can perpetuate inflammation by binding to the receptor for advanced glycation end products (AGER, also known as RAGE). We hypothesized that genetic variation of the AGEs/AGER pathway affects pancreatic cancer risk by modulating chronic inflammation. Methods: We conducted a two-stage case-control study to examine the association between 96 single nucleotide polymorphisms (SNPs) in 21 genes (GLO1, MSR1, PINK1, KIF17, AGPAT1, AKR1B1, DDOST, CD36, SCARB1, AGER, APP, HMGB1, S100A8, S100B, S100P, S100A12, MMP9, ADAM10, PPARG, NOTCH4, and PDGF) of the AGEs/AGER pathway and risk of pancreatic cancer. The discovery study was conducted in 672 pancreatic cancer cases and 1361 controls ascertained from the Women’s Health Initiative (WHI) Study matched on age, race/ethnicity and study arm. The validation study was conducted in a PANC4 pooled hospital-based case-control study of 1,034 women cases and 989 women controls. A pooled analysis of 1706 cases and 2350 controls was also performed. SNP data were obtained or imputed from the previously performed GWASs. Pancreatic cancer risk was calculated as odds ratios (ORs) and 95% confidence intervals (Cis) using logistic regression models adjusting for age, race/ethnicity, body mass index (BMI), type 2 diabetes, and smoking using an additive genetic model. The potential interaction between individual SNP and red meat intake and cigarette smoking was examined. A false discovery rate (FDR) adjusted q value < 0.20 was considered statistically significant in the pooled analyses. We also conducted logistic-Kernel machine (LKM) test to examine the association between sets of SNP and risk of pancreatic cancer in the pooled dataset. Results: We identified 11 SNPs with raw P values < 0.10 in the discovery stage in the WHI study. However, none of the SNPs was validated in the validation and pooled dataset. The LKM test showed that sets of SNPs in CD36 (P = 9.3E-10), GLO1(P =2.0E-7), and DDOST (P = 2.5E-4) were significantly associated with pancreatic cancer risk. We found an interaction effect between the GLO1 SNP rs6932648 and smoking status (never versus ever) in the WHI study (P for interaction = 0.01) and the pooled dataset (P for interaction = 0.049). In the pooled dataset, the variant T allele of rs6932648 (MAF = 0.12) was associated with increased risk of pancreatic cancer among ever smokers (additive OR = 1.26, 95% CI: 1.02-1.55), but not among never smokers (additive OR = 0.94, 95% CI: 0.78-1.14), compared with the C allele. Summary and Conclusion: We did not observe significant association between any of 96 single SNP and pancreatic cancer risk. However, we found combined set of SNPs in the GLO1, DDOST, and CD36 genes were associated with risk of pancreatic cancer. GLO SNP rs6932648 may modify the association between smoking and risk of pancreatic cancer. GLO1, DDOST, and CD36 are involved in detoxification or clearance of AGEs compounds. Further investigation of AGEs and receptors systems in pancreatic cancer is warranted. Citation Format: Guochong Jia, Kathryn Royse, Hongwei Tang, Donghui Li, Liang Chen, Lesley Tinker, Gloria Petersen, Alison Klein, Peter Richardson, Donna White, Haleh Sangi-Haghapeykar, Hashem B. El-Serag, Li Jiao. Genetic variations in the AGEs/AGER pathway and risk of pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-160. doi:10.1158/1538-7445.AM2017-LB-160

Abstract Introduction: Obesity and diabetes are two major modifiable risk factors for pancreatic cancer. Genetic variation affecting the risk of obesity- and diabetes-related pancreatic cancer have yet to be comprehensively investigated at the genome-wide level. Aims: To identify novel gene-environment interactions (GxE) of obesity/diabetes on pancreatic cancer risk. Methods: We conducted a multi-stage genome-wide gene-diabetes and gene-obesity interaction study using GWAS data of 7095 cases and 10710 controls from Pancreatic Cancer Consortium (PanScan I, PanScan II and PanScan III) and Pancreatic Cancer Case Control Consortium (PanC4). About 870,000 SNPs, having MAF ≥ 0.005 and having been genotyped in at least one study, were analyzed. Missing genotypes were imputed using IMPUTE2 and the Haplotype Reference Consortium (HRC) reference panel. Case-control (CC) and case-only (CO) methods were used for SNP-level GxE analysis. Fixed-, random- and joint-effect models built in “rareGE” R package were used for gene-based GxE analysis. About 20,000 genes that have at least one SNP genotyped were included in the analysis. Age, gender and principal components accounting for subpopulation structure were adjusted in both SNP-level and gene-level analyses. Obesity was defined as body mass index ≥30 mg/m2. Diabetics with <3-year duration were excluded for gene-diabetes interaction analysis. The PanScan I, II, III and PanC4 datasets were analyzed separately followed by meta-analysis of the summary statistics. Fixed-effect and Fisher’s p-value based meta-analysis was applied for SNP-level and gene-level GxE analysis, respectively. Results: No genome-wide significant (P<5 x 10-8 ) interaction at the SNP level was detected in any of the PanScan or PanC4 datasets or in the meta-analysis. A significant interaction of diabetes and FAM63A (family with sequence similarity 63 member A) gene (P<2.5 x 10-6) was observed in the meta-analysis (PFixed= 3.4 x 10-5, PRandom= 1.9 x 10-6, PJoint= 6.8 x 10-7) . This interaction showed nominal significance in all individual dataset analyses. The corresponding P values for this interaction were 4.2 x 10-2, 3.8 x 10-2 and 6.8 x 10-2 in PanScan I, 2.8 x 10-2, 3.9 x 10-2 and 6.5 x 10-2 9.1 x 10-3, 2.2 x 10-4 and 8.8 x 10-6 in PanScan III, and 3.1 x 10-3, 3.4 x 10-3 and 8.8 x 10-3 in PanC4, respectively. Conclusions: In the largest genome-wide GxE analysis for pancreatic cancer to date, we identified a significant diabetes interacting gene in pancreatic cancer, which was consistently replicated in 4 individual GWAS studies. This observation identifies FAM63A as a candidate for future functional studies to understand its role in pancreatic cancer etiology. Further increasing the sample size may unveil additional genetic loci that contribute to the susceptibility to pancreatic cancer via G x E. Citation Format: Hongwei Tang, Lai Jiang, Laufey T. Amundadottir, Harvey A. Risch, Rachael S. Stolzenberg-Solomon, Alison P. Klein, Brian M. Wolpin, Gloria Petersen, Nilanjan Chatterjee, Donghui Li, Peter Kraft, Peng Wei. Genomewide gene-diabetes and gene-obesity interaction scan in 7095 cases and 10710 controls from Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case Control Consortium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1295. doi:10.1158/1538-7445.AM2017-1295

Abstract Introduction Pancreatic cancer is a highly lethal malignancy due to late diagnosis and aggressiveness of the tumor. Characterization of high-risk populations for monitoring, intervention and early detection remains challenging. Obesity and diabetes are potentially modifiable contributors to pancreatic cancer. Genetic factors modifying the risk of obesity- and diabetes-induced pancreatic cancer have previously not been fully investigated at the genome-wide level. Hypothesis There are potential genetic factors modifying the associations between obesity/diabetes and pancreatic cancer. Methods Integrating genotype with risk factor data from a GWAS for the Pancreatic Cancer Case Control Consortium, we performed a genome-wide gene-environment interaction (G x E) scan for 457,688 SNPs in a discovery study of 2,028 cases and 2,109 controls using state-of-the-art methods including case-only (CO), case-control (CC), Empirical Hierarchical Bayes (EHB), Empirical-Bayes test (EB) and three 2-step approaches (DG2, EG2, DGEG2). Results We detected a genome-wide significant (P < 1.09 x 10-7 ) interaction of diabetes with rs13061928 (in CNTN4) in CO (P = 8.8 x 10-8 ), which was independent of diabetes in control group (P = 0.08). The minor allele was in a synergism with diabetes on the cancer risk: diabetics carrying GA/AA genotype had a 3.39-fold (95% CI: 2.47-4.65) increased disease risk compared with non-diabetics carrying GG genotype. Consistently, the SNP-diabetes interaction had a high ranking statistic in EHB (No.1) and high p-value ranks for EG2 and DGEG2 methods (No.1, P= 5.44 x 10-5). Conclusions These observations, once validated/confirmed, may help define at-risk subpopulation and initiate targeted intervention and prevention of pancreatic cancer. Genome-wide G x E analysis requires larger sample size than genetic main effects scan; combining the CO method with alternative methods may be the optimal scheme for genome-wide G x E analysis. Citation Format: Hongwei Tang, Eric J. Duell, Harvey A. Risch Risch, Sara H. Olson, H. Bas Bueno-de-Mesquita, Steven Gallinger, Elizabeth A. Holly, Gloria M. Petersen, Paige M. Bracci, Robert R. McWilliams, Mazda Jenab, Elio Riboli, Anne Tjønneland, Marie Christine Boutron-Ruault, Rudolf Kaaks, Dimitrios Trichopoulos, Salvatore Panico, Malin Sund, Petra H. M Peeters, Kay-Tee Khaw, Christopher I Amos, Donghui Li, Peng Wei. Genome-wide gene-diabetes and gene-obesity interaction scan in the pancreatic cancer case control consortium. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 2214. doi:10.1158/1538-7445.AM2014-2214