Risch, Harvey A.

Abstract Genome wide association studies (GWAS) in 9,013 pancreatic cancer patients and 12,452 controls of European ancestry have discovered over 20 risk loci in the human genome. Here, we fine mapped one such locus on chromosome 16q23.1 in the vicinity of two Chymotrypsinogen precursor genes, CTRB1 and CTRB2. We fine-mapped this locus to rs72802365 (P=2.51x10-17, OR=1.36, 95% CI=1.31-1.40) and identified colocalization (PP=87%) with aberrant exon 5→7 CTRB2 splicing in pancreatic tissues. Imputation of a 584 bp insertion/deletion variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (P=2.83x10-16, OR=1.36, 95% CI=1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsin B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). The intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein may lead to ER stress and pancreatic inflammation, which in turn, may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles. Citation Format: Ashley Jermusyk, Jun Zhong, Naomi Gordon, Katelyn Connelly, Sumeth Perera, Ehssan Abdolalizadeh, Tongwu Zhang, Aidan O'Brien, Jason Hoskins, Irene Collins, Daina Eiser, Chen Yuan, Harvey Risch, Eric J. Jacobs, Donghui Li, Mengmeng Du, Rachael Stolzenberg-Solomon, Alison P. Klein, Jill P. Smith, Brian M. Wolpin, Stephen J. Chanock, Jianxin Shi, Gloria M. Petersen, Christopher Westlake, Laufey T. Amundadottir. An exception to the rule: A coding functional variant at a pancreatic cancer GWAS locus [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2020 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2020;80(22 Suppl):Abstract nr IA-11.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni < 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

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 Introduction: Pancreatic cancer is the seventh leading cause of cancer death worldwide with pancreatic ductal adenocarcinoma (PDAC) being the most common subtype (>90%). Inherited genetic changes and cigarette smoking are established independent risk factors of PDAC. Methods: We conducted a genome-wide gene by smoking interaction analysis of PDAC risk among 7,937 PDAC cases and 12,389 controls of European ancestry from the Pancreatic Cancer Cohort (PanScan) and Pancreatic Cancer Case-Control Consortia (PanC4) with complete genotyping and smoking data. Over 6,769,447 SNPs met our inclusion criteria of imputation score >0.5 and MAF >0.05 after imputation to phase 3 1000G and were included in our analysis. We evaluated the gene-environment interactions using the empirical Bayes method which allows data-adaptive relaxation of the gene-environment independence assumption. The logistic regression model was adjusted for age, sex, top principal components for each phase, study and study region (PanScan only). The analysis was conducted separately for the PanScan and PanC4 GWAS studies, using a score test to efficiently incorporate expected SNP dosages. The estimates were then combined using a fixed-effect meta-analysis. We examined expression quantitative trait locus using data from the NIH Genotype-Tissue Expression (GTEx) database. We considered a P-value < 5.0 × 10−8 statistically significant. Results: We found a genome-wide significant qualitative interaction with the SNP rs1818613 by smoking status (P-value <5 × 10−9) where, compared to the G allele, the T allele was associated with a reduced risk PDAC in never smokers (per T allele OR, 0.87; 95% CI, 0.82-0.93), had no association in former smokers (OR, 1.00; 95% CI, 0.91-1.07) and was associated with an increased risk in current smokers (OR, 1.25; 95% CI, 1.12-1.40). This SNP is located on chromosome 2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) gene and upstream of the cyclin T2 (CCNT2) gene. More than 40 additional SNPs within 100Kbp and in high LD (r2 > 0.8) with rs1818613 provided significant evidence of interaction by smoking status for PDAC risk. rs1818613 was associated with differential expression of TMEM163 in several tissues: for T allele compared with the G allele there was increased expression in the heart, pituitary, whole blood, and esophagus-muscularis and decreased expression in the testis tissues (P-value < 3 × 10−11). It was also associated with decreased CCNT2 expression in the tibial nerve tissue (P-value = 1.9 × 10−9). Conclusions: Our work supports that the genetic variation in the region, in conjunction with cigarette smoking, influences PDAC risk. Future studies are needed to explore this association in other smoking-related cancers and to understand the biological mechanism of such association. Citation Format: Prosenjit Kundu, Evelina Mocci, William Wheeler, Laufey T. Amundadottir, Donghui Li, Eric J. Jacobs, Gloria M. Petersen, Brian M. Wolpin, Harvey A. Risch, Peter Kraft, Nilanjan Chatterjee, Alison P. Klein, Rachael Z. Stolzenberg-Solomon, The PanScan and Panc4 consortia. Genome-wide interaction scan identifies gene by smoking interaction at 2q21.3 for pancreatic cancer 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 1213.

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 Parity is associated with decreased risk of invasive ovarian cancer; however, the relationship between incomplete pregnancies and invasive ovarian cancer risk is unclear. This relationship was examined using 15 case-control studies from the Ovarian Cancer Association Consortium (OCAC). Histotype-specific associations, which have not been examined previously with large sample sizes, were also evaluated. Methods A pooled analysis of 10 470 invasive epithelial ovarian cancer cases and 16 942 controls was conducted. Odds ratios (ORs) and 95% confidence intervals (CIs) for the association between incomplete pregnancies and invasive epithelial ovarian cancer were estimated using logistic regression. All models were conditioned on OCAC study, race and ethnicity, age, and education level and adjusted for number of complete pregnancies, oral contraceptive use, and history of breastfeeding. The same approach was used for histotype-specific analyses. Results Ever having an incomplete pregnancy was associated with a 16% reduction in ovarian cancer risk (OR = 0.84, 95% CI = 0.79 to 0.89). There was a trend of decreasing risk with increasing number of incomplete pregnancies (2-sided Ptrend < .001). An inverse association was observed for all major histotypes; it was strongest for clear cell ovarian cancer. Conclusions Incomplete pregnancies are associated with a reduced risk of invasive epithelial ovarian cancer. Pregnancy, including incomplete pregnancy, was associated with a greater reduction in risk of clear cell ovarian cancer, but the result was broadly consistent across histotypes. Future work should focus on understanding the mechanisms underlying this reduced risk.

Abstract Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two‐sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low‐density lipoprotein [LDL] and high‐density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels ( P < 5 × 10 −8 ) were identified as instrumental variables, and assessed using genome‐wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non‐endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non‐endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non‐endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non‐endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

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 &lt; 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 &lt; 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.