Brezina, Stefanie

Background: Considering the pathophysiological importance of the intestine-passing microbiota, we hypothesized that the mutational and microbial profiles of colorectal cancer might gradually change along detailed anatomical locations and that these trends might differ by the age of diagnosis and the microsatellite instability (MSI) status. Methods: We assessed somatic mutations and selected bacterial species in colorectal carcinomas along the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum, utilizing targeted sequencing data ( N = 5685) from the Genetics and Epidemiology of Colorectal Cancer Consortium and tissue microbial data ( N = 1389) from the Health Professionals Follow-up Study and Nurses’ Health Study. Results: There were decreasing trends from the cecum to rectum in the prevalence of gene mutations in the IGF/PI3K, receptor tyrosine kinase (RTK)/RAS, TGFB, and WNT signaling pathways (all P trend < 0.0001), as well as high levels of Bacteroides fragilis ( P trend < 0.0001) and high levels of Fusobacterium nucleatum ( P trend = 0.0029) (but not enterotoxigenic Bacteroides fragilis or pks+ Escherichia coli ). Conversely, TP53 pathway mutations showed an increasing cecum-to-rectum trend ( P trend < 0.0001). Similar to later-onset cases diagnosed at ages 50 years and older, in early-onset cases diagnosed under the age of 50 years, we found generally decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, TGFB, and WNT pathway mutations and an increasing cecum-to-rectum trend in TP53 pathway mutations. By contrast, the cecum-to-rectum prevalence-changing trends in IGF/PI3K, RTK/RAS, TGFB, and TP53 pathway mutations were present in non–MSI-high tumors but not in MSI-high tumors (all P interaction < 0.0052). Conclusion: Non–MSI-high colorectal cancers exhibited decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, and TGFB pathway mutations (and an increasing trend of TP53 pathway mutations), in contrast to MSI-high tumors that demonstrated no such trends. Our study highlights the importance of detailed colorectal subsite information combined with tumor molecular and microbial characteristics in colorectal cancer research.

Abstract We used deep learning (DL) to predict multiple genetic biomarkers of colorectal cancer (CRC) from routine histopathology images, aiming to identify morphological changes associated with the mutational profile of tumors. A new Transformer-based prediction model was developed which can simultaneously predict multiple biomarkers from a single histopathology image. We compared the performance of the multi-target approach to conventional single-target models. We analyzed five patient cohorts from the Genetics and Epidemiology of CRC Consortium (GECCO) with colorectal cancer (N=1,385 patients). Our model was trained on an internal data set (N=739) to predict the presence of &gt;200 genetic alterations, focusing on genes with non-silent mutations and mutational signatures, which were assessed using targeted sequencing. We used a 7-fold cross validation and validated the model on an independent, external data set (N=646). Morphological features of microsatellite instability (MSI) are detectable by DL models in histopathology images, making MSI a potential confounding factor. Thus, we also assessed whether our model can predict genetic alterations regardless of MSI status, or whether the model identifies the MSI phenotype. For the majority of biomarkers, the multi-target transformer reached a higher prediction performance as measured by the mean Area Under the Receiver Operating Characteristic (AUROC) Curve and lower standard deviation than single-target DL models. For the external validation set, the mean AUROC of the multi-target versus the single-target model was 0.78 (+/- 0.01) versus 0.72 (+/- 0.06) for BRAF mutations, 0.88 (+/- 0.01) versus 0.86 (+/- 0.03) for hypermutation status, 0.94 (+/- 0.01) versus 0.91 (+/- 0.02) for MSI, 0.86 (+/- 0.01) versus 0.80 (+/- 0.05) for RNF43 mutations and 0.72 (+/- 0.02) versus 0.69 (+/- 0.05) for TP53 mutations, respectively. However, our analysis of the individual prediction scores suggests that the model does not identify the phenotype of the alterations themselves. Instead, the model seems to determine the status of specific genetic alterations based on their correlation with MSI. Consequently, biomarkers which are associated with microsatellite stability have prediction scores negatively correlated with predicted MSI scores. Conversely, biomarkers which are associated with MSI have prediction scores positively correlated with predicted MSI scores. Our results demonstrate potential benefits of multi-target DL models in improving the predictive power and efficacy of histology-based biomarker extraction compared to single-target DL models. However, our results crucially highlight the importance of considering potential confounders, such as MSI status in CRC, and emphasize the need for extensive analysis beyond AUROC in studies focusing on DL-based biomarker detection. Citation Format: Marco Gustav, Marko van Treeck, Zunamys I. Carrero, Chiara M. Loeffler, Nic G. Reitsam, Bruno Märkl, Asier Rabasco Meneghetti, Lisa A. Boardman, Amy J. French, Ellen L. Goode, Andrea Gsur, Stefanie Brezina, Marc J. Gunter, Neil Murphy, Paul Limburg, Stephen Thibodeau, Sebastian Foersch, Robert Steinfelder, Tabitha Harrison, Ulrike Peters, Amanda Phipps, Jakob N. Kather. Assessing microsatellite instability dominance in colorectal cancer phenotype: A multi-study initiative using multi-target transformers for genomic biomarker prediction [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4924.

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.

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