Citations

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

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.

Abstract It was not known whether the polygenic risk scores (PRSs) that predict colorectal cancer could predict colorectal cancer for people with inherited pathogenic variants in DNA mismatch repair genes—people with Lynch syndrome. We tested a PRS comprising 107 established single-nucleotide polymorphisms associated with colorectal cancer in European populations for 826 European-descent carriers of pathogenic variants in DNA mismatch repair genes (293 MLH1, 314 MSH2, 126 MSH6, 71 PMS2, and 22 EPCAM) from the Colon Cancer Family Registry, of whom 504 had colorectal cancer. There was no evidence of an association between the PRS and colorectal cancer risk, irrespective of which DNA mismatch repair gene was mutated, or sex (all 2-sided P > .05). The hazard ratio per standard deviation of the PRS for colorectal cancer was 0.97 (95% confidence interval = 0.88 to 1.06; 2-sided P = .51). Whereas PRSs are predictive of colorectal cancer in the general population, they do not predict Lynch syndrome colorectal cancer.

Objective An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined. Design To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling. Results We identified 13 loci that reached genome-wide significance (p<5×10 −8 ) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer. Conclusion Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.

Abstract Background: Inherited susceptibility is an important contributor to colorectal cancer risk, and rare variants in key genes or pathways could account in part for the missing proportion of colorectal cancer heritability. Methods: We conducted an exome-wide association study including 2,327 cases and 2,966 controls of European ancestry from three large epidemiologic studies. Single variant associations were tested using logistic regression models, adjusting for appropriate study-specific covariates. In addition, we examined the aggregate effects of rare coding variation at the gene and pathway levels using Bayesian model uncertainty techniques. Results: In an exome-wide gene-level analysis, we identified ST6GALNAC2 as the top associated gene based on the Bayesian risk index (BRI) method [summary Bayes factor (BF)BRI = 2604.23]. A rare coding variant in this gene, rs139401613, was the top associated variant (P = 1.01 × 10–6) in an exome-wide single variant analysis. Pathway-level association analyses based on the integrative BRI (iBRI) method found extreme evidence of association with the DNA repair pathway (BFiBRI = 17852.4), specifically with the nonhomologous end joining (BFiBRI = 437.95) and nucleotide excision repair (BFiBRI = 36.96) subpathways. The iBRI method also identified RPA2, PRKDC, ERCC5, and ERCC8 as the top associated DNA repair genes (summary BFiBRI ≥ 10), with rs28988897, rs8178232, rs141369732, and rs201642761 being the most likely associated variants in these genes, respectively. Conclusions: We identified novel variants and genes associated with colorectal cancer risk and provided additional evidence for a role of DNA repair in colorectal cancer tumorigenesis. Impact: This study provides new insights into the genetic predisposition to colorectal cancer, which has potential for translation into improved risk prediction.