Citations

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.

AbstractIntroduction/AimsALS is a heterogeneous disease that may be complicated or in part driven by inflammation. NP001, a regulator of macrophage activation, was associated with slowing disease progression in those with higher levels of the plasma inflammatory marker C‐reactive protein (CRP) in phase 2A studies in ALS. Here, we evaluate the effects of NP001 in a phase 2B trial, and perform a post hoc analysis with combined data from the preceding phase 2A trial.MethodsThe phase 2B trial enrolled 138 participants within 3 y of symptom onset and with plasma hs‐CRP values >1.13 mg/L. They were randomized 1:1 to receive either placebo or NP001 for 6 mo. Change from baseline ALSFRS‐R scores was the primary efficacy endpoint. Secondary endpoints included vital capacity (VC) change from baseline and percentage of participants showing no decline of ALSFRS‐R score over 6 mo (non‐progressor).ResultsThe phase 2B study did not show significant differences between placebo and active treatment with respect to change in ALSFRS‐R scores, or VC. The drug was safe and well tolerated. A post hoc analysis identified a 40‐ to 65‐y‐old subset in which NP001‐treated patients demonstrated slower declines in ALSFRS‐R score by 36% and VC loss by 51% compared with placebo. A greater number of non‐progressors were NP001‐treated compared with placebo (p = .004).DiscussionAlthough the phase 2B trial failed to meet its primary endpoints, post hoc analyses identified a subgroup whose decline in ALSFRS‐R and VC scores were significantly slower than placebo. Further studies will be required to validate these findings.

Abstract Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10–8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82–0.93; former smokers 1.00, 95% CI, 0.91–1.07; current smokers 1.25, 95% CI 1.12–1.40, Pinteraction = 3.08 × 10–9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. Significance: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

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 Pancreatic ductal adenocarcinoma (PDAC) is one of most lethal malignancies with few known risk factors and biomarkers. Identification of disease biomarkers is critical for understanding the pathogenesis of this cancer and identifying high risk individuals for close surveillance. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers usually in small samples. To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, by using genetic instruments. Protein quantitative trait loci (pQTLs) for 1,226 plasma proteins identified in a large INTERVAL study of 3,301 healthy European descendants were used as instruments to evaluate associations between genetically predicted protein levels and PDAC. For proteins showing a significant association, we further conducted conditional analysis with adjustments for previously identified risk variants to assess whether the observed associations between genetically predicted protein concentrations and PDAC risk were independent of the risk variants identified in genome-wide association studies (GWAS). Furthermore, for the proteins that were associated with PDAC risk, we performed an enrichment analysis of the genes encoding these proteins to examine whether they are enriched in specific pathways, functions or networks. We observed associations between predicted concentrations of 38 proteins and PDAC risk at a false discovery rate of < 0.05, including those of 23 proteins that showed a significant association even after Bonferroni correction (4.08 × 10−5). These include Histo-blood group ABO system transferase encoded by ABO, which has been previously implicated as a potential target gene of PDAC risk variant identified in GWAS. Eight of the identified proteins (Beta-crystallin B2, Dedicator of cytokinesis protein 9, VIP36-like protein, Erythrocyte band 7 integral membrane protein, Tensin-2, Transmembrane protease serine 11D, Alcohol dehydrogenase 1B, and C-X-C motif chemokine 10) were associated with PDAC risk after conditioning on previously reported pancreatic cancer risk variants (odds ratios ranged from 0.79 to 1.52, P-values from 1.28 × 10−3 to 6.47 × 10−4). Pathway enrichment analysis showed that the encoding genes for the implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL-15 production. In conclusion, we identified 38 protein biomarker candidates for PDAC risk, which if validated by additional studies, may contribute to the etiological understanding of PDAC tumor development. Citation Format: Jingjing Zhu, Xiang Shu, Xingyi Guo, Duo Liu, Jiandong bao, Roger Milne, Graham G Giles, Chong Wu, Mengmeng Du, Emily White, Harvey A Risch, Nuria Malats, Eric J. Duell, Phyllis J. Goodman, Donghui Li, Paige Bracci, Verena Katzke, Rachel E Neale, Steven Gallinger, Stephen Van Den Eeden, Alan Arslan, Federico Canzian, Charles Kooperberg, Brian Wolpin, Laura Beane-Freeman, Ghislaine Scelo, Kala Visvanatha, Christopher A. Haiman, Loïc Le Marchand, Herbert Yu, Gloria M Petersen, Rachael Stolzenberg-Solomon, Alison P Klein, Laufey T Amundadottir, Qiuyin Cai, Jirong Long, Xiao-Ou Shu, Wei Zheng, Lang Wu. Associations between genetically predicted blood protein biomarkers and pancreatic ductal adenocarcinoma risk [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1200.

Abstract Background: Genome-wide association studies (GWAS) revealed that 25 regions in 24 genes are associated with adult diffuse glioma development. These regions were identified by performing GWAS of glioma overall and by pathology (GBM and nonGBM). The 2016 WHO Classification of Tumors of the Central Nervous System utilizes two somatic alterations to molecularly-classify adult diffuse glioma: IDH mutation and 1p/19q codeletion. We hypothesized that germline variants may increase susceptibility to, or interact with, these somatic alterations to accelerate the development of specific molecular subtypes of glioma. We further hypothesize that germline variants associated with IDH-mutated glioma might be associated with other IDH-mutated tumors, namely, cholangiocarcinoma, acute myeloid leukemia (AML) and melanoma. Methods: We performed a GWAS by glioma molecular subtype - as defined by presence or absence of IDH somatic mutation and 1p/19q codeletion. A total of 1320 glioma cases and 1889 controls were used in the discovery set, and 799 glioma cases and 808 controls in the validation set. A meta-analysis was performed with a genome-wide p-value threshold of 5 × 10−8. GTEx data were used to perform an expression quantitative trait loci (eQTL) analysis. For germline variants that were significantly associated with IDH-mutated glioma, we evaluated pleiotropy with cholangiocarcinoma, AML and melanoma using TCGA and Mayo Biobank controls. Results: Variants in or near D2HGDH on chromosome 2 were genome-wide significant in IDH-mutated glioma (meta p-value = 2.82 × 10−10). TCGA reported that the D2HGDH region was commonly deleted in IDH-mutated gliomas that do not have 1p/19q codeletion. In TCGA data for IDH-mutated, non-codeleted glioma, we observed that the D2HGDH variant was inversely associated with tumor deletions of D2HGDH (odds ratio=0.57, p-value=0.015). The eQTL analyses demonstrated significant associations between D2HGDH germline variant and expression of D2HGDH (p=2.2 × 10−11). Further stratifying IDH-mutated glioma by 1p/19q codeletion status, one variant near FAM20C on chromosome 7 was genome-wide significant in gliomas that have IDH mutation and 1p/19q codeletion (meta p-value=9.56 × 10−9). Analyses are currently underway to evaluate pleiotropy of these IDH-mutated glioma germline variants with other IDH-mutated tumors including cholangiocarcinoma, AML and melanoma. Variants in or near GMEB2 on chromosome 20 were genome-wide significant in IDH wild-type glioma (meta p-value=2.60 × 10−10). The most significant variant in the GMEB2 region remained significant after adjustment for the known RTEL1 glioma risk variant nearby on chromosome 20 (p=0.029). Conclusions: We identified and validated novel germline variants in two genes that are associated with etiology of IDH-mutated and one gene that is associated with IDH wild-type adult diffuse glioma. Citation Format: Jeanette E. Eckel-Passow, Kristen L. Drucker, Thomas M. Kollmeyer, Matthew L. Kosel, Paul A. Decker, Annette M. Molinaro, Terri Rice, Corrine E. Praska, Lauren E. Clark, Alissa A. Caron, Alexej Abyzov, Anthony Batzler, Jun S. Song, Melike Pekmezci, Helen M. Hansen, Lucie S. McCoy, Paige M. Bracci, Joseph Wiemels, John K. Wiencke, Stephen Francis, Terence C. Burns, Caterina Giannini, Daniel H. Lachance, Margaret Wrensch, Robert B. Jenkins. Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH, FAM20C and GMEB2 [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1193.

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 Registry-based epidemiologic studies suggest associations between chronic inflammatory intestinal diseases and pancreatic ductal adenocarcinoma (PDAC). As genetic susceptibility contributes to a large proportion of chronic inflammatory intestinal diseases, we hypothesize that the genomic regions surrounding established genome-wide associated variants for these chronic inflammatory diseases are associated with PDAC. We examined the association between PDAC and genomic regions (±500 kb) surrounding established common susceptibility variants for ulcerative colitis, Crohn's disease, inflammatory bowel disease, celiac disease, chronic pancreatitis, and primary sclerosing cholangitis. We analyzed summary statistics from genome-wide association studies data for 8,384 cases and 11,955 controls of European descent from two large consortium studies using the summary data-based adaptive rank truncated product method to examine the overall association of combined genomic regions for each inflammatory disease group. Combined genomic susceptibility regions for ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis were associated with PDAC at P values < 0.05 (0.0040, 0.0057, 0.011, and 3.4 × 10−6, respectively). After excluding the 20 PDAC susceptibility regions (±500 kb) previously identified by GWAS, the genomic regions for ulcerative colitis, Crohn disease, and inflammatory bowel disease remained associated with PDAC (P = 0.0029, 0.0057, and 0.0098, respectively). Genomic regions for celiac disease (P = 0.22) and primary sclerosing cholangitis (P = 0.078) were not associated with PDAC. Our results support the hypothesis that genomic regions surrounding variants associated with inflammatory intestinal diseases, particularly, ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis are associated with PDAC. Significance: The joint effects of common variants in genomic regions containing susceptibility loci for inflammatory bowel disease and chronic pancreatitis are associated with PDAC and may provide insights to understanding pancreatic cancer etiology.

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: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, with few known risk factors and biomarkers. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers, usually in small samples. In this study, we evaluated associations of circulating protein levels and PDAC risk using genetic instruments. Methods: To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, using genetic instruments of protein quantitative trait loci. Results: We observed associations between predicted concentrations of 38 proteins and PDAC risk at an FDR of < 0.05, including 23 of those proteins that showed an association even after Bonferroni correction. These include the protein encoded by ABO, which has been implicated as a potential target gene of PDAC risk variant. Eight of the identified proteins (LMA2L, TM11D, IP-10, ADH1B, STOM, TENC1, DOCK9, and CRBB2) were associated with PDAC risk after adjusting for previously reported PDAC risk variants (OR ranged from 0.79 to 1.52). Pathway enrichment analysis showed that the encoding genes for implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL15 production. Conclusions: We identified 38 candidates of protein biomarkers for PDAC risk. Impact: This study identifies novel protein biomarker candidates for PDAC, which if validated by additional studies, may contribute to the etiologic understanding of PDAC development.