Lin, Yi

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 Protective associations of fruits, vegetables, and fiber intake with colorectal cancer risk have been shown in many, but not all epidemiologic studies. One possible reason for study heterogeneity is that dietary factors may have distinct effects by colorectal cancer molecular subtypes. Here, we investigate the association of fruit, vegetables, and fiber intake with four well-established colorectal cancer molecular subtypes separately and in combination. Nine observational studies including 9,592 cases with molecular subtypes for microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and somatic mutations in BRAF and KRAS genes, and 7,869 controls were analyzed. Both case-only logistic regression analyses and polytomous logistic regression analyses (with one control set and multiple case groups) were used. Higher fruit intake was associated with a trend toward decreased risk of BRAF-mutated tumors [OR 4th vs. 1st quartile = 0.82 (95% confidence interval, 0.65–1.04)] but not BRAF-wildtype tumors [1.09 (0.97–1.22); P difference as shown in case-only analysis = 0.02]. This difference was observed in case–control studies and not in cohort studies. Compared with controls, higher fiber intake showed negative association with colorectal cancer risk for cases with microsatellite stable/MSI-low, CIMP-negative, BRAF-wildtype, and KRAS-wildtype tumors (Ptrend range from 0.03 to 3.4e-03), which is consistent with the traditional adenoma-colorectal cancer pathway. These negative associations were stronger compared with MSI-high, CIMP-positive, BRAF-mutated, or KRAS-mutated tumors, but the differences were not statistically significant. These inverse associations for fruit and fiber intake may explain, in part, inconsistent findings between fruit or fiber intake and colorectal cancer risk that have previously been reported. Significance: These analyses by colorectal cancer molecular subtypes potentially explain the inconsistent findings between dietary fruit or fiber intake and overall colorectal cancer risk that have previously been reported.

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.

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.

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 Genome-wide association studies of colorectal cancer (CRC) have identified over 50 susceptibility loci. These variants represent only a small fraction of total heritability for CRC. Gene-environment (GxE) interaction studies may help identify novel loci and biological interactions that give insight to the pathogenesis of CRC. Previous genome-wide GxE studies with dietary factors have identified interactions between loci and processed meat consumption and alcohol; however, limited statistical power remains a primary concern. Set-based SNP testing has the potential to increase statistical power to detect GxE interactions by aggregating functionally relevant SNPs. In this large pooled analysis using 14 case-control and nested case-control studies, we incorporated functional information from the transcriptome prediction tool, PrediXcan, into a novel set-based approach for testing GxE interactions. We used variant weights from the PrediXcan models of tissue-specific gene expression in the transverse colon as a priori variant information for a set-based GxE approach. We restricted our analysis to variants in the PrediXcan transverse colon gene models (n = 4,842). This discovery phase included 10,360 CRC and advanced adenoma cases, and 11,183 controls of European ancestry from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry. All 14 studies were analyzed together in a pooled data set using the Mixed Effects Score Tests for interactions. We tested for gene interactions with sex- and study- specific quartiles of dietary intake of red meat (servings/day), processed meat (servings/day), vegetables (servings/day), fruits (servings/day), and fiber (g/day). We detected two genes with suggestive interactions (false discovery rate (FDR) &lt; 0.2) with intake of red meat and risk of CRC: Superoxide Dismutase 2 (SOD2) and Ubiquitin Conjugating Enzyme E2 H (UBE2H). No interactions at FDR &lt; 0.2 were observed for processed meat, vegetables, fruits, or fiber. The SOD2 gene, which encodes an enzyme important in apoptotic signaling and clearing of reactive oxygen species, may regulate response to colonic exposure to heme iron or increased bile acid from high-fat content in red meat. UBE2H is part of the ubiquitin-proteasome system that has a role in Wnt signaling, which can be mediated by heme iron in red meat and is commonly found dysregulated in cancer. These findings highlight the efficacy of integrating functional information and set-based testing for novel discovery of genes interacting with known dietary risk factors of CRC. We plan to replicate these findings in additional studies. Citation Format: Paneen S. Petersen, Yu-Ru Su, Sonja I. Berndt, Stephanie A. Bien, Hermann Brenner, Graham Casey, Andrew T. Chan, Jenny Chang-Claude, Jane C. Figueiredo, Steven J. Gallinger, Robert W. Haile, Tabitha A. Harrison, Michael Hoffmeister, Mark A. Jenkins, Amit D. Joshi, Sébastien Küry, Loic Le Marchand, Yi Lin, Noralane M. Lindor, Polly A. Newcomb, John D. Potter, Robert Schoen, Martha L. Slattery, Stephen N. Thibodeau, Emily White, Li Hsu, Ulrike Peters, CCFR, GECCO. Interactions between genetic predictors of gene expression and dietary factors associated with risk of colorectal 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 5268.

Abstract Background Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10−8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk. Methods We conducted a GWAS in European descent CRC cases and control subjects using a discovery–replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10−8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided. Results The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10−8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0. Conclusions This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

Abstract Purpose: Clustering of breast and colorectal cancer has been observed within families in population-based studies giving rise to speculation that there are “breast-colon” cancer susceptibility genes. We performed a genome-wide association study (GWAS) to identify genetic markers associated with a potential breast-colon cancer phenotype. Methods: Cases and controls were ascertained from the Breast and Colon Cancer Family Registries’ (CFR) GWAS subjects. “Breast-colon phenotype” definition was based on 2 or more breast and colon cancer cases diagnosed in first- or second-degree relatives within a family. Cases (n = 985) were women with a history of breast cancer and a family history of colorectal cancer, or persons with colorectal cancer with a family history of breast cancer. Unrelated controls (n = 1769) were frequency matched to cases for age and gender. Following standard quality control measures, 6,220,060 directly measured and imputed single nucleotide polymorphisms (SNPs) were included in the discovery set. SNPs associated at p&lt;1×10−5 were analyzed in a replication dataset of cases (n = 293) and controls (n = 2103) from the Genetics and Epidemiology of Colorectal Cancer Consortium. Results: In a regression model that included gender, age, CFR center and 6 principal components, the top association in the discovery set was in the chromosome 8q22.3 region overlying the BAALC gene (top SNP rs12548629, P = 3.63×10−7). In the replication dataset, of 341 SNPs tested, the top signal was found in multiple correlated SNPs overlying the ROBO1 gene on chromosome 3p12 (rs7429100, P = 2.8×10−3), however, the signal on chromosome 8, the BAALC gene did not replicate. In the discovery set, P values for the SNPs overlaying ROBO1 gene ranged from 2.2×10−5 to 9.7×10−5 (rs7429100, P = 3.9×10−5). The combined meta-analysis showed strongest association at the ROBO1 gene (rs7429100, Pfixed effects 1.84×10−6, Pheterogeneity &gt;0.05 for testing for heterogeneity between the overall discovery set and the replication set). ROBO1 (roundabout, axon guidance receptor, homolog 1) is a transmembrane receptor of the immunoglobulin family and is differentially expressed in human cancers, with a possible role as a tumor suppressor gene. Low ROBO1 expression has been shown to be an adverse prognostic factor for invasive ductal breast cancer and may also play a role in pathogenesis of colorectal cancer. BAALC (brain and acute leukemia, cytoplasmic) is predominantly expressed by neural and blood cells and is largely implicated in acute leukemia. Conclusion: In this exploratory analysis, to elucidate genes/regions associated with pleiotropic effect for breast and colorectal cancer risk we identified germline variation in the region of ROBO1 and BAALC. Validation in a larger dataset and functional characterization of the loci is warranted to elucidate mechanisms by which these genes/SNPs may contribute to the development of breast and colorectal cancer. Citation Format: Mala Pande, Aron Joon, Sanjay Shete, Abenaa M. Brewster, Cathy Eng, Wei V. Chen, Habibul Ahsan, Irene L. Andrulis, Esther M. John, Yi Lin, Polly A. Newcomb, Noralane M. Lindor, Christopher I. Amos, John Hopper, Patrick M. Lynch. Risk loci for a breast-colon cancer phenotype: results from a genome-wide association study. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2782. doi:10.1158/1538-7445.AM2015-2782