Citations

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.

Purpose The known epithelial ovarian cancer (EOC) susceptibility genes account for less than 50% of the heritable risk of ovarian cancer suggesting that other susceptibility genes exist. The aim of this study was to evaluate the contribution to ovarian cancer susceptibility of rare deleterious germline variants in a set of candidate genes. Methods We sequenced the coding region of 54 candidate genes in 6385 invasive EOC cases and 6115 controls of broad European ancestry. Genes with an increased frequency of putative deleterious variants in cases versus controls were further examined in an independent set of 14 135 EOC cases and 28 655 controls from the Ovarian Cancer Association Consortium and the UK Biobank. For each gene, we estimated the EOC risks and evaluated associations between germline variant status and clinical characteristics. Results The ORs associated for high-grade serous ovarian cancer were 3.01 for PALB2 (95% CI 1.59 to 5.68; p=0.00068), 1.99 for POLK (95% CI 1.15 to 3.43; p=0.014) and 4.07 for SLX4 (95% CI 1.34 to 12.4; p=0.013). Deleterious mutations in FBXO10 were associated with a reduced risk of disease (OR 0.27, 95% CI 0.07 to 1.00, p=0.049). However, based on the Bayes false discovery probability, only the association for PALB2 in high-grade serous ovarian cancer is likely to represent a true positive. Conclusions We have found strong evidence that carriers of PALB2 deleterious mutations are at increased risk of high-grade serous ovarian cancer. Whether the magnitude of risk is sufficiently high to warrant the inclusion of PALB2 in cancer gene panels for ovarian cancer risk testing is unclear; much larger sample sizes will be needed to provide sufficiently precise estimates for clinical counselling.

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.

Abstract Background Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. Methods To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74–421 samples). Results We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. Conclusions By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Abstract Background: Pancreatic cancer is the fourth-leading cause of cancer death in both men and women in the United States. The currently identified common susceptibility loci account for a small fraction of estimated heritability. We sought to estimate overall heritability of pancreatic cancer and partition the heritability by variant frequencies and functional annotations. Methods: Analysis using the genome-based restricted maximum likelihood method (GREML) was conducted on Pancreatic Cancer Case-Control Consortium (PanC4) genome-wide association study (GWAS) data from 3,568 pancreatic cancer cases and 3,363 controls of European Ancestry. Results: Applying linkage disequilibrium- and minor allele frequency-stratified GREML (GREML-LDMS) method to imputed GWAS data, we estimated the overall heritability of pancreatic cancer to be 21.2% (SE = 4.8%). Across the functional groups (intronic, intergenic, coding, and regulatory variants), intronic variants account for most of the estimated heritability (12.4%). Previously identified GWAS loci explained 4.1% of the total phenotypic variation of pancreatic cancer. Mutations in hereditary pancreatic cancer susceptibility genes are present in 4% to 10% of patients with pancreatic cancer, yet our GREML-LDMS results suggested these regions explain only 0.4% of total phenotypic variance for pancreatic cancer. Conclusions: Although higher than previous studies, our estimated 21.2% overall heritability may still be downwardly biased due to the inherent limitation that the contribution of rare variants in genes with a substantive overall impact on disease are not captured when applying these commonly used methods to imputed GWAS data. Impact: Our work demonstrated the importance of rare and common variants in pancreatic cancer risk.

Abstract Objective Genome-wide association studies (GWAS) identify associations of individual SNPs with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants has been shown for many disorders to be a powerful tool in discovering novel networks of susceptibility genes. Design We have conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product (sARTP) method to identify novel gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDA) in 9,040 cases and 12,495 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes identified in PDA-associated gene sets and pathways. Results We identified 14 gene sets and pathways associated with PDA at FDR < 0.05. Five of the strongest signals (P-value ≤ 1.3x10-5 Bonferroni corrected) were in genes from KEGG maturity onset diabetes of the young, Reactome regulation of beta cell development, Biocarta role of epidermal growth factor (EGF) receptor transactivation by G-protein-coupled receptors in cardiac hypertrophy pathway, Nikolsky breast cancer chr17q11-q21 amplicon gene set, and Pujana ATM Pearson correlation coefficient (PCC) network gene set. We identified and validated rs876493 and three correlated SNPs (PGAP3) and rs3124737 (CASP7) as eQTLs in two normal derived pancreas tissue datasets. Conclusion Our agnostic gene set and pathway analysis integrated with functional annotation and eQTL analysis provides novel insight into genes and pathways that may be biologically relevant for risk of PDA. Citation Format: Naomi Walsh, Han Zhang, Paula L. Hyland, Qi Yang, Evelina Mocci, Mingfeng Zhang, Erica J. Childs, Zhaoming Wang, Stephen Chanock, Patricia Hartge, Robert Hoover, Peter Kraft, Donghui Li, Eric J. Jacobs, Gloria M. Petersen, Brian M. Wolpin, Harvey A. Risch, Laufey T. Amundadottir, Kai Yu, Alison P. Klein, Rachael Z. Stolzenberg-Solomon. Large pathway and gene set analysis of GWAS data identifies novel associations for pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 242.

AbstractBackground: Pathogenic germline mutations in the CDKN2A tumor suppressor gene are rare and associated with highly penetrant familial melanoma and pancreatic cancer in non-Hispanic whites (NHW). To date, the prevalence and impact of CDKN2A rare coding variants (RCV) in racial minority groups remain poorly characterized. We examined the role of CDKN2A RCVs on the risk of pancreatic cancer among minority subjects.Methods: We sequenced CDKN2A in 220 African American (AA) pancreatic cancer cases, 900 noncancer AA controls, and 183 Nigerian controls. RCV frequencies were determined for each group and compared with that of 1,537 NHW patients with pancreatic cancer. Odds ratios (OR) and 95% confidence intervals (CI) were calculated for both a case–case comparison of RCV frequencies in AAs versus NHWs, and case–control comparison between AA cases versus noncancer AA controls plus Nigerian controls. Smaller sets of Hispanic and Native American cases and controls also were sequenced.Results: One novel missense RCV and one novel frameshift RCV were found among AA patients: 400G>A and 258_278del. RCV carrier status was associated with increased risk of pancreatic cancer among AA cases (11/220; OR, 3.3; 95% CI, 1.5–7.1; P = 0.004) compared with AA and Nigerian controls (17/1,083). Further, AA cases had higher frequency of RCVs: 5.0% (OR, 13.4; 95% CI, 4.9–36.7; P < 0.001) compared with NHW cases (0.4%).Conclusions: CDKN2A RCVs are more common in AA than in NHW patients with pancreatic cancer and associated with moderately increased pancreatic cancer risk among AAs.Impact: RCVs in CDKN2A are frequent in AAs and are associated with risk for pancreatic cancer. Cancer Epidemiol Biomarkers Prev; 27(11); 1364–70. ©2018 AACR.