Newcomb, Polly A.

Abstract Background Observational studies have shown higher folate consumption to be associated with lower risk of colorectal cancer (CRC). Understanding whether and how genetic risk factors interact with folate could further elucidate the underlying mechanism. Aggregating functionally relevant genetic variants in set-based variant testing has higher power to detect gene-environment (GxE) interactions and may provide information on the underlying biological pathway. Objective We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. Methods We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based GxE approach. We harmonized total folate intake (mcg/day) based on dietary intake and supplemental use across cohort and case-control studies and calculated sex and study specific quantiles. Analyses were performed using a mixed effects score tests for interactions between folate and genetically predicted expression of 4,839 genes with available genetically predicted expression. We meta-analyzed results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. Results We found suggestive evidence of interaction with folate intake for genes including glutathione S-Transferase Alpha 1 (GSTA1; p=4.3E-4), Tonsuko Like, DNA Repair Protein (TONSL; p=4.3E-4), and Aspartylglucosaminidase (AGA: p=4.5E-4). Glutathione is an antioxidant, preventing cellular damage and is a downstream metabolite of homocysteine and metabolized by GSTA1. TONSL is part of a complex that functions in the recovery of double strand breaks and AGA plays a role in lysosomal breakdown of glycoprotein. Conclusion We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk.

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes’ coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84–1.27 and OR: 1.08, 95% CI: 0.86–1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

Abstract Background: Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis. Methods: We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer–specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer–specific mortality, and somatic mutations. Results: F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer–specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07–3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer–specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations. Conclusions: F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer–specific mortality and specific somatic mutated genes. Impact: Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

Background Tumor-infiltrating lymphocytes (TILs), especially CD8+ TILs, can be used for predicting immunotherapy responsiveness and survival outcome. However, the evaluation of CD8+ TILs currently relies on histopathological methodology with high variability. We therefore aimed to develop a DNA methylation signature for CD8+ TILs (CD8+ MeTIL) that could evaluate immune response and prognosis in colorectal cancer (CRC). Methods A CD8+ MeTIL signature score was constructed by using CD8+ T cell-specific differentially methylated positions (DMPs) that were identified from Illumina EPIC methylation arrays. Immune cells, colon epithelial cells, and two CRC cohorts (n=282 and 335) were used to develop a PCR-based assay for quantitative analysis of DNA methylation at single-base resolution (QASM) to determine CD8 + MeTIL signature score. Results Three CD8+ T cell-specific DMPs were identified to construct the CD8+ MeTIL signature score, which showed a dramatic discriminability between CD8+ T cells and other cells. The QASM assay we developed for CD8+ MeTIL markers could measure CD8+ TILs distributions in a fully quantitative, accurate, and simple manner. The CD8+ MeTIL score determined by QASM assay showed a strong association with histopathology-based CD8+ TIL counts and a gene expression-based immune marker. Furthermore, the low CD8+ MeTIL score (enriched CD8+ TILs) was associated with MSI-H tumors and predicted better survival in CRC cohorts. Conclusions This study developed a quantitative DNA methylation-based signature that was reliable to evaluate CD8+ TILs and prognosis in CRC. This approach has the potential to be a tool for investigations on CD8+ TILs and a biomarker for therapeutic approaches, including immunotherapy.

Abstract Background Individuals who carry pathogenic mutations in DNA mismatch repair (MMR) genes have high risks of cancer, and small studies have suggested that these risks depend on the sex of the parent from whom the mutation was inherited. We have conducted the first large study of such a parent-of-origin effect (POE). Methods Our study was based on all MMR gene mutation carriers and their relatives in the Colon Cancer Family Registry, comprising 18,226 people. The POE was estimated as a hazard ratio (HR) using a segregation analysis approach that adjusted for ascertainment. HR = 1 corresponds to no POE and HR > 1 corresponds to higher risks for maternal mutations. Results For all MMR genes combined, the estimated POE HRs were 1.02 (95% confidence interval (CI) 0.75-1.39, p = 0.9) for male colorectal cancer, 1.12 (95% CI 0.81-1.54, p = 0.5) for female colorectal cancer and 0.84 (95% CI 0.52-1.36, p = 0.5) for endometrial cancer. Separate results for each MMR gene were similar. Conclusions Despite being well-powered, our study did not find any evidence that cancer risks for MMR gene mutation carriers depend on the parent-of-origin of the mutation. Based on current evidence, we don’t recommend that POEs be incorporated into the clinical guidelines or advice for such carriers. Key messages MMR gene mutations inherited from the maternal and paternal side confer similar risks of developing colorectal and endometrial cancer.

Abstract Colorectal cancer (CRC) is a heterogenous disease that develops through somatic mutations in driver genes, leading to activation of diverse neoplastic pathways. To systematically examine if somatically mutated genes and pathways impact survival, we sequenced tumor and normal DNA samples for 4,512 CRC cases using a targeted panel. We performed Cox regression to estimate hazard ratios (HR) and 95% confidence intervals (CIs) for associations of disease-specific (DS)-survival with somatically mutated genes and pathways, adjusting for age and sex, and stratifying baseline hazards by study population. We assessed statistical significance using Bonferroni p-value thresholds to account for multiple testing of 214 genes (2.34x10-4) and 6 key CRC pathways (8.3x10-3). We limited analyses to non-silent mutations. We observed that DS-survival was significantly more favorable among individuals with hypermutated (HM) tumors, primarily consisting of microsatellite unstable and POLE-mutated tumors, as compared to those with non-hypermutated (NHM) tumors (HR=0.4, 95% CI: 0.3-0.5, p=1.2x10-17). BRAF V600E mutations were associated with poorer survival (HR=2.0, 95% CI: 1.6-2.5, p=1.6x10-10). This association was more pronounced among NHM tumors (HR=2.3, 95% CI: 1.8-2.9, p=4.3x10-12). We identified suggestive associations (p-values < 5.0x10-3 in overall or stratified analyses) between DS-survival and mutations in B2M, TP53, and SMAD4. Mutations in B2M may provide more favorable prognosis for survival (HR=0.5, 95% CI: 0.3-0.8, p=4.4x10-3), with similar effect sizes in HM and NHM tumors. Poorer survival may be associated with mutations in TP53 (HR=1.2, 95% CI:1.2-1.4, p=9.3x10-4) and SMAD4 (HR=1.3, 95% CI:1.1-1.6, p=3.2x10-3) in NHM tumors. We further observed statistically significant associations between survival and three mutated pathways: TP53/ATM (HR=1.2, 95% CI:1.1-1.4, p=9.0x10-4), receptor tyrosine kinases (RTK) and RAS (HR=1.3, 95% CI:1.1-1.5, p=4.5x10-5), and TGF-beta (HR=1.3, 95% CI:1.1-1.5, p=2.9x10-4). Findings for TP53/ATM and RTK/RAS were primarily due to mutations in one gene within each pathway (TP53 and BRAF, respectively). However, the TGF-beta pathway finding, which was more pronounced in NHM tumors, included two genes with p-values less than 0.05 (SMAD4 p=3.2x10-3 and TGFBR2 p=0.01). Despite our large sample size, relatively few somatically mutated genes were significantly associated with DS-survival. It may be that pathway-level testing affords more power for analyses and that larger sample sizes are needed. It is of interest that non-silent mutations in B2M, whose product is a component of the class I major histocompatibility complex, was associated with improved DS-survival, though our finding was not statistically significant after multiple testing correction so further exploration and replication is needed. Citation Format: Tabitha A. Harrison, Syed H. Zaidi, Conghui Qu, Amanda I. Phipps, Robert S. Steinfelder, Quang M. Trinh, Sonja I. Berndt, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Kim F. Doheny, David A. Drew, Jane C. Figueiredo, Steven J. Gallinger, Marios Giannakis, Andrea Gsur, Marc J. Gunter, Michael Hoffmeister, Wen-Yi Huang, Paul J. Limburg, Victor Moreno, Polly A. Newcomb, Shuji Ogino, Ross L. Prentice, Tameka Shelford, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB090.

Background. Lynch syndrome is the most common genetic predisposition for hereditary cancer. Carriers of pathogenic changes in mismatch repair (MMR) genes have an increased risk of developing colorectal (CRC), endometrial, ovarian, urinary tract, prostate, and other cancers, depending on which gene is malfunctioning. In Lynch syndrome, differences in cancer incidence (penetrance) according to the gene involved have led to the stratification of cancer surveillance. By contrast, any differences in penetrance determined by the type of pathogenic variant remain unknown. Objective. To determine cumulative incidences of cancer in carriers of truncating and missense or aberrant splicing pathogenic variants of the MLH1 and MSH2 genes. Methods. Carriers of pathogenic variants of MLH1 (path_MLH1) and MSH2 (path_MSH2) genes filed in the Prospective Lynch Syndrome Database (PLSD) were categorized as truncating or missense/aberrant splicing according to the InSiGHT criteria for pathogenicity. Results. Among 5199 carriers, 1045 had missense or aberrant splicing variants, and 3930 had truncating variants. Prospective observation years for the two groups were 8205 and 34,141 years, respectively, after which there were no significant differences in incidences for cancer overall or for colorectal cancer or endometrial cancers separately. Conclusion. Truncating and missense or aberrant splicing pathogenic variants were associated with similar average cumulative incidences of cancer in carriers of path MLH1 and path_MSH2.

Abstract Background Smoking is associated with colorectal cancer (CRC) risk. Previous studies suggested this association may be restricted to certain molecular subtypes of CRC, but large-scale comprehensive analysis is lacking. Methods A total of 9789 CRC cases and 11 231 controls of European ancestry from 11 observational studies were included. We harmonized smoking variables across studies and derived sex study–specific quartiles of pack-years of smoking for analysis. Four somatic colorectal tumor markers were assessed individually and in combination, including BRAF mutation, KRAS mutation, CpG island methylator phenotype (CIMP), and microsatellite instability (MSI) status. A multinomial logistic regression analysis was used to assess the association between smoking and risk of CRC subtypes by molecular characteristics, adjusting for age, sex, and study. All statistical tests were 2-sided and adjusted for Bonferroni correction. Results Heavier smoking was associated with higher risk of CRC overall and stratified by individual markers (Ptrend < .001). The associations differed statistically significantly between all molecular subtypes, which was the most statistically significant for CIMP and BRAF. Compared with never-smokers, smokers in the fourth quartile of pack-years had a 90% higher risk of CIMP-positive CRC (odds ratio = 1.90, 95% confidence interval = 1.60 to 2.26) but only 35% higher risk for CIMP-negative CRC (odds ratio = 1.35, 95% confidence interval = 1.22 to 1.49; Pdifference = 2.1 x 10-6). The association was also stronger in tumors that were CIMP positive, MSI high, or KRAS wild type when combined (Pdifference < .001). Conclusion Smoking was associated with differential risk of CRC subtypes defined by molecular characteristics. Heavier smokers had particularly higher risk of CRC subtypes that were CIMP positive and MSI high in combination, suggesting that smoking may be involved in the development of colorectal tumors via the serrated pathway.

PURPOSE This study assessed the joint association of pathogenic variants (PVs) in breast cancer (BC) predisposition genes and polygenic risk scores (PRS) with BC in the general population. METHODS A total of 26,798 non-Hispanic white BC cases and 26,127 controls from predominately population-based studies in the Cancer Risk Estimates Related to Susceptibility consortium were evaluated for PVs in BRCA1, BRCA2, ATM, CHEK2, PALB2, BARD1, BRIP1, CDH1, and NF1. PRS based on 105 common variants were created using effect estimates from BC genome-wide association studies; the performance of an overall BC PRS and estrogen receptor–specific PRS were evaluated. The odds of BC based on the PVs and PRS were estimated using penalized logistic regression. The results were combined with age-specific incidence rates to estimate 5-year and lifetime absolute risks of BC across percentiles of PRS by PV status and first-degree family history of BC. RESULTS The estimated lifetime risks of BC among general-population noncarriers, based on 10th and 90th percentiles of PRS, were 9.1%-23.9% and 6.7%-18.2% for women with or without first-degree relatives with BC, respectively. Taking PRS into account, more than 95% of BRCA1, BRCA2, and PALB2 carriers had > 20% lifetime risks of BC, whereas, respectively, 52.5% and 69.7% of ATM and CHEK2 carriers without first-degree relatives with BC, and 78.8% and 89.9% of those with a first-degree relative with BC had > 20% risk. CONCLUSION PRS facilitates personalization of BC risk among carriers of PVs in predisposition genes. Incorporating PRS into BC risk estimation may help identify > 30% of CHEK2 and nearly half of ATM carriers below the 20% lifetime risk threshold, suggesting the addition of PRS may prevent overscreening and enable more personalized risk management approaches.

Abstract Background: A positive association between circulating C-reactive protein (CRP) and colorectal cancer survival was reported in observational studies, which are susceptible to unmeasured confounding and reverse causality. We used a Mendelian randomization approach to evaluate the association between genetically predicted CRP concentrations and colorectal cancer–specific survival. Methods: We used individual-level data for 16,918 eligible colorectal cancer cases of European ancestry from 15 studies within the International Survival Analysis of Colorectal Cancer Consortium. We calculated a genetic-risk score based on 52 CRP-associated genetic variants identified from genome-wide association studies. Because of the non-collapsibility of hazard ratios from Cox proportional hazards models, we used the additive hazards model to calculate hazard differences (HD) and 95% confidence intervals (CI) for the association between genetically predicted CRP concentrations and colorectal cancer–specific survival, overall and by stage at diagnosis and tumor location. Analyses were adjusted for age at diagnosis, sex, body mass index, genotyping platform, study, and principal components. Results: Of the 5,395 (32%) deaths accrued over up to 10 years of follow-up, 3,808 (23%) were due to colorectal cancer. Genetically predicted CRP concentration was not associated with colorectal cancer–specific survival (HD, −1.15; 95% CI, −2.76 to 0.47 per 100,000 person-years; P = 0.16). Similarly, no associations were observed in subgroup analyses by stage at diagnosis or tumor location. Conclusions: Despite adequate power to detect moderate associations, our results did not support a causal effect of circulating CRP concentrations on colorectal cancer–specific survival. Impact: Future research evaluating genetically determined levels of other circulating inflammatory biomarkers (i.e., IL6) with colorectal cancer survival outcomes is needed.