Citations

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.

AbstractBackgroundPostmenopausal hormone therapy (HT) is associated with a decreased colorectal cancer (CRC) risk. As CRC is a heterogeneous disease, we evaluated whether the association of HT and CRC differs across etiologically relevant, molecularly defined tumor subtypes and tumor location.MethodsWe pooled data on tumor subtypes (microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, pathway: adenoma-carcinoma, alternate, serrated), tumor location (proximal colon, distal colon, rectum), and HT use among 8220 postmenopausal women (3898 CRC cases and 4322 controls) from 8 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CIs) for the association of ever vs never HT use with each tumor subtype compared with controls. Models were adjusted for study, age, body mass index, smoking status, and CRC family history. All statistical tests were 2-sided.ResultsAmong postmenopausal women, ever HT use was associated with a 38% reduction in overall CRC risk (OR =0.62, 95% CI = 0.56 to 0.69). This association was similar according to microsatellite instability, CpG island methylator phenotype and BRAF or KRAS status. However, the association was attenuated for tumors arising through the serrated pathway (OR = 0.81, 95% CI = 0.66 to 1.01) compared with the adenoma-carcinoma pathway (OR = 0.63, 95% CI = 0.55 to 0.73; Phet =.04) and alternate pathway (OR = 0.61, 95% CI = 0.51 to 0.72). Additionally, proximal colon tumors had a weaker association (OR = 0.71, 95% CI = 0.62 to 0.80) compared with rectal (OR = 0.54, 95% CI = 0.46 to 0.63) and distal colon (OR = 0.57, 95% CI = 0.49 to 0.66; Phet =.01) tumors.ConclusionsWe observed a strong inverse association between HT use and overall CRC risk, which may predominantly reflect a benefit of HT use for tumors arising through the adenoma-carcinoma and alternate pathways as well as distal colon and rectal tumors.

AbstractDeficient intake of micronutrients involved in one‐carbon metabolism (eg, choline, methionine, vitamin B12 and folic acid) leads to hepatocellular carcinoma (HCC) development in rodents, but is under‐investigated in humans. We investigated the association between one‐carbon metabolism‐related micronutrient intake and HCC risk in a prospective cohort of 494 860 participants with 16 years of follow‐up in the NIH‐AARP study. Dietary intakes and supplement use were ascertained at baseline using a food‐frequency questionnaire. Total intake (diet plus supplements) of the following one‐carbon metabolism‐related micronutrients were calculated: folate, methionine and vitamins B2 (riboflavin), B3 (niacin), B6 and B12. These micronutrients were examined both individually and simultaneously, with adjustment for covariates. Cox proportional hazard models were used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs). Over the 16‐year follow‐up period, 647 incident HCC cases were diagnosed. When examined individually, higher total vitamin B3 intake was associated with a lower HCC risk (HRQ5 vs Q1 = 0.60; 95% CI = 0.42‐0.85; Ptrend = .008), and the association remained significant when all six micronutrients were examined simultaneously (HRQ5 vs Q1 = 0.32; 95% CI = 0.18‐0.55; Ptrend < .0001). Among participants with >3 years of follow‐up, higher total vitamin B3 intake was again associated with lower risk (HRQ5 vs Q1 = 0.37; 95% CI = 0.20‐0.68; Ptrend = .001), whereas higher total vitamin B6 intake was associated with higher risk (HRQ5 vs Q1 = 2.04; 95% CI = 1.02‐4.07; Ptrend = .04). Restricted cubic spline analyses showed a dose‐response inverse association between total vitamin B3 intake and HCC risk, and dose‐response positive association between total vitamin B6 intake and HCC risk. The study suggests that higher vitamin B3 intake is associated with lower HCC risk, whereas higher vitamin B6 intake is associated with increased risk.

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: 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 Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3′ region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25–2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63–5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71–4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12–11.54; P = 0.0002). No genotype–phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. Significance: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

Abstract Background: Smoking has been associated with colorectal cancer (CRC) risk; but limited evidence has shown the association between smoking and molecular subtypes of CRC. Methods: We analyzed 9,422 CRC cases and 9,950 controls from 10 observational studies. Multinomial logistic regression analysis was performed to assess the association between sex-study-specific quartiles of pack years of smoking and risk of CRC molecular subtypes, using non-smokers as reference group, adjusting for age, sex, and study. Known oncogenic mutations in four somatic colorectal tumor markers were assessed individually and in combination, including BRAF mutations, KRAS mutations, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. Case-only analysis was also performed to estimate heterogeneity in risk of molecular subtypes of CRC. Results: Compared with controls, higher pack years of smoking were statistically significantly associated with increased CRC risk when stratified, individually, by all four markers, and the associations got stronger with higher quartiles pack years (p-trend&lt;0.001). Associations between smoking and CRC risk also differed significantly among molecular subtypes. Compared to nonsmokers, the risk of BRAFmut CRC was 83% higher for smokers within the highest quartile of pack-years (OR=1.83; 95% CI: 1.50, 2.25), and 29% higher for BRAFwt CRC (OR=1.29; 95% CI: 1.17, 1.43; Ratio of ORs (ROR)=1.45; 95% CI: 1.18, 1.17; p=3.45x10-4). Similarly, heavy pack-years of smoking was associated with almost two times higher risk of CIMP-high CRC (OR=1.93; 95% CI: 1.60, 2.31), but only 33% higher risk of CIMP-low/negative CRC (OR=1.33; 95% CI: 1.19, 1.48; ROR=1.49; 95% CI: 1.24, 1.79; p=1.72x10-5). The association between smoking and CRC was also stronger in MSI-high tumors (ORMSI-H=1.65; 95% CI: 1.36, 2.00; ORMSI-L/MSS=1.37; 95% CI: 1.23, 1.52; ROR=1.22; 95% CI: 1.00, 1.48; p=0.046). In contrast, the association between smoking and CRC risk was stronger for KRASwt (OR=1.43; 95% CI: 1.27, 1.60), than KRASmut tumors (OR=1.18; 95% CI: 1.02, 1.37; ROR=0.83; 95% CI: 0.71, 0.97; p =0.016). When combining tumor markers, smoking was found to be significantly associated with higher risk of colorectal tumors from the serrated pathway. Conclusion: In this largest study with a total of 19,372 subjects, we found that heavier pack years of smoking was associated with increased risk of all CRC molecular subtypes. Smokers with heavier pack-years of smoking had particularly higher risk of CRC subtypes with BRAF mutation and CIMP-high, suggesting smoking may be particularly involved in the development of these subtypes of colorectal tumor. Citation Format: Xiaoliang Wang, Efrat Amitay, Barbara L. Banbury, Hermann Brenner, Daniel D. Buchanan, Peter T. Campbell, Jenny Chang-Claude, Steven J. Gallinger, Graham G. Giles, Tabitha A. Harrison, John L. Hopper, Mark A. Jenkins, Yi Lin, Reiko Nishihara, Polly A. Newcomb, Shuji Ogino, Lori C. Sakoda, Robert E. Schoen, Martha L. Slattery, Steven N. Thibodeau, Bethany Van Guelpen, Michael O. Woods, Li Hsu, Michael Hoffmeister, Ulrike Peters. Smoking is associated with risks of molecular subtypes of colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 630.

Interindividual differences in DNA repair systems may play a role in modulating the individual risk of developing colorectal cancer. To better ascertain the role of DNA repair gene polymorphisms on colon and rectal cancer risk individually, we evaluated 15,419 single nucleotide polymorphisms (SNPs) within 185 DNA repair genes using GWAS data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), which included 8,178 colon cancer, 2,936 rectum cancer cases and 14,659 controls. Rs1800734 (in MLH1 gene) was associated with colon cancer risk (p‐value = 3.5 × 10−6) and rs2189517 (in RAD51B) with rectal cancer risk (p‐value = 5.7 × 10−6). The results had statistical significance close to the Bonferroni corrected p‐value of 5.8 × 10−6. Ninety‐four SNPs were significantly associated with colorectal cancer risk after Binomial Sequential Goodness of Fit (BSGoF) procedure and confirmed the relevance of DNA mismatch repair (MMR) and homologous recombination pathways for colon and rectum cancer, respectively. Defects in MMR genes are known to be crucial for familial form of colorectal cancer but our findings suggest that specific genetic variations in MLH1 are important also in the individual predisposition to sporadic colon cancer. Other SNPs associated with the risk of colon cancer (e.g., rs16906252 in MGMT) were found to affect mRNA expression levels in colon transverse and therefore working as possible cis‐eQTL suggesting possible mechanisms of carcinogenesis.

Abstract Over 50 genetic variants have been associated with colorectal cancer (CRC) risk through genome-wide association studies (GWAS), yet these variants represent only a fraction of the total estimated heritability. CRC is a heterogenous disease with diverse tumor etiology. Assessing genetic risk in molecular subtypes may help to identify novel loci and characterize genetic risk among tumor subtypes. We used microsatellite instability (MSI), an established CRC classifier with etiological and therapeutic relevance, to define CRC subtypes for GWAS analyses. We conducted a case-case analysis to estimate odds ratios (OR) and 95% confidence intervals (CI) for association of genome-wide variants with microsatellite stable (MSS) versus unstable (MSI) carcinomas. We ran an inverse-variance weighted fixed-effects meta-analysis across GWAS in a discovery set of 4,163 population-based CRC cases with harmonized microsatellite instability (MSI) marker and imputed genotype data. For each analysis, we used log-additive logistic regression, adjusting for age, sex, and principal components to account for population substructure. We then followed up with replication of 102 SNPs that reached p-values less than 5x10-6 in 1,698 cases. A total of 845 (20.3%) cancer cases were microsatellite unstable in the discovery population and 174 (10.2%) were unstable in the replication population. No variants reached the genome-wide significance level of 5x10-8 in the discovery set. However, we identified two variants that reached a Bonferroni corrected p-value of 4.0x10-4 in the replication set. This included one variant in MLH1 (Replication: OR=1.74, 95% CI=1.53-1.98, p=1.63x10-5; Discovery+Replication: OR=1.45, 95% CI=1.37-1.54, p=9.76x10-11) and one variant in LOC105377645 (Replication: OR=1.70, 95% CI=1.49-1.94, p=5.13x10-5; Discovery+Replication: OR=1.45, 95% CI=1.37-1.54, p=9.76 x 10-11). The MLH1 gene is a DNA mismatch repair gene implicated in Lynch Syndrome, the hallmark of which is microsatellite instability. This is the first genome-wide scan to identify a common variant in MLH1 that is associated with CRC. This variant (minor allele frequency, MAF = 23% in this all European ancestry population) is located in the 5'-untranslated region of MLH1 and is thought to act as a long-range regulator of DCLK3, a potential tumor driver gene. The second variant, located in LOC105377645 with an MAF of 22%, is in an uncharacterized region of the genome and has not previously been implicated in cancer development. These findings suggest that accounting for molecular heterogeneity is important for discovery and characterization of genetic variants associated with CRC risk. We plan to run polytomous regression analyses, increase our sample size, and further investigate CRC subtypes by CIMP, BRAF mutation, KRAS mutation status. Citation Format: Tabitha A. Harrison, Yiwen Lu, Chenjie Zeng, Flora Qu, Kristin Anderson, Hermann Brenner, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Jenny Chang-Claude, Graham G. Giles, Bethany Van Guelpen, Michael Hoffmeister, Mark A. Jenkins, Noralane M. Lindor, Roger L. Milne, Polly A. Newcomb, Reiko Nishihara, Michael O. Woods, Shuji Ogino, John D. Potter, Martha L. Slattery, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Genome-wide association study by colorectal carcinoma subtype [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 229.