Brenner, Hermann

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.

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.

Abstract We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast 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.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

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: 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.

1562 Background: In Caucasian-predominant populations, overweight or obese NSCLC pts (BMI≥25-kg/m2) have better prognosis while underweight (BMI≤18.5-kg/m2) pts have a worse prognosis. A large pooled sample allowed us to evaluate the role of ethnicity in this BMI-NSCLC OS relationship. Methods: Using individual data, survival analysis was performed on 15 ILCCO studies, assessing the interactions between ethnicity and BMI on overall survival (OS). Adjusted Hazard Ratios (aHR) from Cox models and adjusted penalized smoothing spline plots were generated. Results: Among 13416 (77%) Caucasian, 2975 (17%) Asian, and 935 (5%) Black NSCLC pts analyzed, we confirmed that for all pts, being underweight at NSCLC diagnosis was associated with worse OS (aHR=1.68, CI 1.5-1.8, p<0.001), while overweight/obese pts had improved survival (aHR 0.89, CI 0.8-0.9, p<0.001), when compared to pts with normal BMI. In general, Black pts had poorer OS than Caucasian pts (aHR 1.26 CI 1.1-1.4, p<0.001). However, the BMI-OS relationship differed according to ethnicity (BMI-ethnicity interaction, p=0.009): Caucasian underweight pts had poorer OS (aHR 1.67, CI 1.5-1.8, P<0.001) while overweight/obese pts had improved OS (aHR 0.89, CI 0.8-0.9, P<0.001). In Asian pts, these two associations were aHR 1.15, CI 0.8-1.5, P=0.33, and aHR 0.91, CI 0.7-1.1, P=0.32, respectively. In Black pts, these two associations were aHR 1.06, CI 0.7-1.5, p=0.73 and 0.75, CI 0.6-0.8 p<0.001, respectively. For overweight/obese patients, as BMI rises from 25-kg/m2 through 45-kg/m2, the prognosis worsens for Caucasian pts, remains stable for Asian pts, but improves for Black pts. Conclusions: InCaucasian pts, being underweight had a greater negative impact on OS than for Asian or Black pts, while being obese had a greater beneficial impact in Blacks than in other ethnic groups. These ethnic differences likely reflect genetically-informed muscle/adipose tissue distributions, where Black pts may have less sarcopenic obesity than other ethnicities. In future prognostic studies, BMI relationships must account for ethnic differences.