Citations

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: Evidence for aspirin's chemopreventative properties on colorectal cancer (CRC) is substantial, but its mechanism of action is not well-understood. We combined a proteomic approach with Mendelian randomization (MR) to identify possible new aspirin targets that decrease CRC risk. Methods: Human colorectal adenoma cells (RG/C2) were treated with aspirin (24 hours) and a stable isotope labeling with amino acids in cell culture (SILAC) based proteomics approach identified altered protein expression. Protein quantitative trait loci (pQTLs) from INTERVAL (N = 3,301) and expression QTLs (eQTLs) from the eQTLGen Consortium (N = 31,684) were used as genetic proxies for protein and mRNA expression levels. Two-sample MR of mRNA/protein expression on CRC risk was performed using eQTL/pQTL data combined with CRC genetic summary data from the Colon Cancer Family Registry (CCFR), Colorectal Transdisciplinary (CORECT), Genetics and Epidemiology of Colorectal Cancer (GECCO) consortia and UK Biobank (55,168 cases and 65,160 controls). Results: Altered expression was detected for 125/5886 proteins. Of these, aspirin decreased MCM6, RRM2, and ARFIP2 expression, and MR analysis showed that a standard deviation increase in mRNA/protein expression was associated with increased CRC risk (OR: 1.08, 95% CI, 1.03–1.13; OR: 3.33, 95% CI, 2.46–4.50; and OR: 1.15, 95% CI, 1.02–1.29, respectively). Conclusions: MCM6 and RRM2 are involved in DNA repair whereby reduced expression may lead to increased DNA aberrations and ultimately cancer cell death, whereas ARFIP2 is involved in actin cytoskeletal regulation, indicating a possible role in aspirin's reduction of metastasis. Impact: Our approach has shown how laboratory experiments and population-based approaches can combine to identify aspirin-targeted proteins possibly affecting CRC risk.

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

Abstract Background: Regular use of nonsteroidal anti-inflammatory drugs (NSAID) is associated with lower risk of colorectal cancer. Genome-wide interaction analysis on single variants (G × E) has identified several SNPs that may interact with NSAIDs to confer colorectal cancer risk, but variations in gene expression levels may also modify the effect of NSAID use. Therefore, we tested interactions between NSAID use and predicted gene expression levels in relation to colorectal cancer risk. Methods: Genetically predicted gene expressions were tested for interaction with NSAID use on colorectal cancer risk among 19,258 colorectal cancer cases and 18,597 controls from 21 observational studies. A Mixed Score Test for Interactions (MiSTi) approach was used to jointly assess G × E effects which are modeled via fixed interaction effects of the weighted burden within each gene set (burden) and residual G × E effects (variance). A false discovery rate (FDR) at 0.2 was applied to correct for multiple testing. Results: Among the 4,840 genes tested, genetically predicted expression levels of four genes modified the effect of any NSAID use on colorectal cancer risk, including DPP10 (PG×E = 1.96 × 10−4), KRT16 (PG×E = 2.3 × 10−4), CD14 (PG×E = 9.38 × 10−4), and CYP27A1 (PG×E = 1.44 × 10−3). There was a significant interaction between expression level of RP11-89N17 and regular use of aspirin only on colorectal cancer risk (PG×E = 3.23 × 10−5). No interactions were observed between predicted gene expression and nonaspirin NSAID use at FDR < 0.2. Conclusions: By incorporating functional information, we discovered several novel genes that interacted with NSAID use. Impact: These findings provide preliminary support that could help understand the chemopreventive mechanisms of NSAIDs on colorectal cancer.

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.

AbstractBackgroundBody mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology.MethodsTo improve statistical power and interpretation for gene‐environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype‐Tissue Expression Project for all genes with heritability ≥1%. We used a mixed‐effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI‐CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2.ResultsAmong 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10−5), PSMC5 (P = 4.51 × 10−4) and CD33 (P = 2.71 × 10−4) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10−5) and SCN1B (P = 2.76 × 10−4) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10−5).ConclusionsAggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk.

Abstract Background: Lipid traits have been inconsistently linked to risk of non-Hodgkin lymphoma (NHL). We examined the association of genetically predicted lipid traits with risk of diffuse large B-cell lymphoma (DLBCL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and marginal zone lymphoma (MZL) using Mendelian randomization (MR) analysis. Methods: Genome-wide association study data from the InterLymph Consortium were available for 2,661 DLBCLs, 2,179 CLLs, 2,142 FLs, 824 MZLs, and 6,221 controls. SNPs associated (P &lt; 5 × 10−8) with high-density lipoprotein (HDL, n = 164), low-density lipoprotein (LDL, n = 137), total cholesterol (TC, n = 161), and triglycerides (TG, n = 123) were used as instrumental variables (IV), explaining 14.6%, 27.7%, 16.8%, and 12.8% of phenotypic variation, respectively. Associations between each lipid trait and NHL subtype were calculated using the MR inverse variance–weighted method, estimating odds ratios (OR) per standard deviation and 95% confidence intervals (CI). Results: HDL was positively associated with DLBCL (OR = 1.14; 95% CI, 1.00–1.30) and MZL (OR = 1.09; 95% CI, 1.01–1.18), while TG was inversely associated with MZL risk (OR = 0.90; 95% CI, 0.83–0.99), all at nominal significance (P &lt; 0.05). A positive trend was observed for HDL with FL risk (OR = 1.08; 95% CI, 0.99–1.19; P = 0.087). No associations were noteworthy after adjusting for multiple testing. Conclusions: We did not find evidence of a clear or strong association of these lipid traits with the most common NHL subtypes. While these IVs have been previously linked to other cancers, our findings do not support any causal associations with these NHL subtypes. Impact: Our results suggest that prior reported inverse associations of lipid traits are not likely to be causal and could represent reverse causality or confounding.

Abstract Background: Results from epidemiologic studies examining polyunsaturated fatty acids (PUFA) and colorectal cancer risk are inconsistent. Mendelian randomization may strengthen causal inference from observational studies. Given their shared metabolic pathway, examining the combined effects of aspirin/NSAID use with PUFAs could help elucidate an association between PUFAs and colorectal cancer risk. Methods: Information was leveraged from genome-wide association studies (GWAS) regarding PUFA-associated SNPs to create weighted genetic scores (wGS) representing genetically predicted circulating blood PUFAs for 11,016 non-Hispanic white colorectal cancer cases and 13,732 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Associations per SD increase in the wGS were estimated using unconditional logistic regression. Interactions between PUFA wGSs and aspirin/NSAID use on colorectal cancer risk were also examined. Results: Modest colorectal cancer risk reductions were observed per SD increase in circulating linoleic acid [ORLA = 0.96; 95% confidence interval (CI) = 0.93–0.98; P = 5.2 × 10−4] and α-linolenic acid (ORALA = 0.95; 95% CI = 0.92–0.97; P = 5.4 × 10−5), whereas modest increased risks were observed for arachidonic (ORAA = 1.06; 95% CI = 1.03–1.08; P = 3.3 × 10−5), eicosapentaenoic (OREPA = 1.04; 95% CI = 1.01–1.07; P = 2.5 × 10−3), and docosapentaenoic acids (ORDPA = 1.03; 95% CI = 1.01–1.06; P = 1.2 × 10−2). Each of these effects was stronger among aspirin/NSAID nonusers in the stratified analyses. Conclusions: Our study suggests that higher circulating shorter-chain PUFAs (i.e., LA and ALA) were associated with reduced colorectal cancer risk, whereas longer-chain PUFAs (i.e., AA, EPA, and DPA) were associated with an increased colorectal cancer risk. Impact: The interaction of PUFAs with aspirin/NSAID use indicates a shared colorectal cancer inflammatory pathway. Future research should continue to improve PUFA genetic instruments to elucidate the independent effects of PUFAs on colorectal cancer.

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