Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study

preprint OA: gold CC-BY-4.0
📄 Open PDF Full text JSON View at publisher

Abstract

Background Deep Learning (DL) has emerged as a powerful tool to predict genetic biomarkers directly from digitized Hematoxylin and Eosin (H&E) slides in colorectal cancer (CRC). However, few studies have systematically investigated the predictability of biomarkers beyond routinely available alterations such as microsatellite instability (MSI), and BRAF and KRAS mutations. Methods Our primary dataset comprised H&E slides of CRC tumors across five cohorts totaling 1,376 patients who underwent comprehensive panel sequencing, with an additional 536 patients from two public datasets for validation. We developed a DL model using a single transformer model to predict multiple genetic alterations directly from the slides. The model’s performance was compared against conventional single-target models, and potential confounders were analyzed. Findings The multi-target model was able to predict numerous biomarkers from pathology slides, matching and partly exceeding single-target transformers. The Area Under the Receiver Operating Characteristic curve (AUROC, mean ± std) on the primary external validation cohorts was: BRAF (0·78 ± 0·01), hypermutation (0·88 ± 0·01), MSI (0·93 ± 0·01), RNF43 (0·86 ± 0·01); this biomarker predictability was mirrored across metrics and co-occurrence analyses. However, biomarkers with high AUROCs largely correlated with MSI, with model predictions depending considerably on MSI-associated morphology upon pathological examination. Interpretation Our study demonstrates that multi-target transformers can predict the biomarker status for numerous genetic alterations in CRC directly from H&E slides. However, their pre-dictability is mainly associated with MSI phenotype, despite indications of slight biomarker-inherent contributions to a phenotype. Our findings underscore the need to analyze confounders in AI-based oncology biomarkers. To enable this, we developed a validated model applicable to other cancers and larger, diverse datasets. Funding The German Federal Ministry of Health, the Max-Eder-Programme of German Cancer Aid, the German Federal Ministry of Education and Research, the German Academic Exchange Service, and the EU.
Full text 81,436 characters · extracted from preprint-html · click to expand
Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study View ORCID Profile Marco Gustav , Marko van Treeck , Nic G. Reitsam , Zunamys I. Carrero , Chiara M. Loeffler , Asier Rabasco Meneghetti , Bruno Märkl , Lisa A. Boardman , Amy J. French , Ellen L. Goode , Andrea Gsur , Stefanie Brezina , Marc J. Gunter , Neil Murphy , Pia Hönscheid , Christian Sperling , Sebastian Foersch , Robert Steinfelder , Tabitha Harrison , Ulrike Peters , Amanda Phipps , View ORCID Profile Jakob Nikolas Kather doi: https://doi.org/10.1101/2025.02.04.25321660 Marco Gustav 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Marco Gustav Marko van Treeck 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Nic G. Reitsam 2 Pathology, Faculty of Medicine, University of Augsburg , Augsburg, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Zunamys I. Carrero 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Chiara M. Loeffler 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany 3 Department of Medicine I, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Asier Rabasco Meneghetti 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Bruno Märkl 2 Pathology, Faculty of Medicine, University of Augsburg , Augsburg, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Lisa A. Boardman 4 Division of Gastroenterology and Hepatology, Mayo Clinic , Rochester, Minnesota, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Amy J. French 5 Division of Laboratory Genetics, Department of Laboratory Medicine and Pathology, Mayo Clinic , Rochester, Minnesota, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Ellen L. Goode 6 Department of Quantitative Health Sciences, Division of Epidemiology, Mayo Clinic , Rochester, Minnesota, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Andrea Gsur 7 Center for Cancer Research, Medical University of Vienna , Vienna, Austria Find this author on Google Scholar Find this author on PubMed Search for this author on this site Stefanie Brezina 7 Center for Cancer Research, Medical University of Vienna , Vienna, Austria Find this author on Google Scholar Find this author on PubMed Search for this author on this site Marc J. Gunter 8 Nutrition and Metabolism Branch, International Agency for Research on Cancer, World Health Organization , Lyon, France 9 Cancer Epidemiology and Prevention Research Unit, School of Public Health, Imperial College London , London, United Kingdom Find this author on Google Scholar Find this author on PubMed Search for this author on this site Neil Murphy 8 Nutrition and Metabolism Branch, International Agency for Research on Cancer, World Health Organization , Lyon, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site Pia Hönscheid 10 Institute of Pathology, University Hospital Carl Gustav Carus (UKD), Technical University Dresden (TUD) , Dresden, Germany 11 National Center for Tumor Diseases (NCT), Partner Site Dresden, German Cancer Research Center Heidelberg , Dresden, Germany 12 German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ) , Heidelberg, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Christian Sperling 10 Institute of Pathology, University Hospital Carl Gustav Carus (UKD), Technical University Dresden (TUD) , Dresden, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sebastian Foersch 13 Institute of Pathology, University Medical Center Mainz , Mainz, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site Robert Steinfelder 14 Division of Public Health Sciences, Fred Hutchinson Cancer Center , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Tabitha Harrison 14 Division of Public Health Sciences, Fred Hutchinson Cancer Center , Seattle, WA, USA 15 Department of Epidemiology, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Ulrike Peters 14 Division of Public Health Sciences, Fred Hutchinson Cancer Center , Seattle, WA, USA 15 Department of Epidemiology, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Amanda Phipps 14 Division of Public Health Sciences, Fred Hutchinson Cancer Center , Seattle, WA, USA 15 Department of Epidemiology, University of Washington , Seattle, WA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jakob Nikolas Kather 1 Else Kroener Fresenius Center for Digital Health, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany 3 Department of Medicine I, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology , 01307 Dresden, Germany 16 Medical Oncology, National Center for Tumor Diseases (NCT), University Hospital Heidelberg , Heidelberg, Germany 17 Pathology & Data Analytics, Leeds Institute of Medical Research at St James’s, University of Leeds , Leeds, United Kingdom Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jakob Nikolas Kather For correspondence: jakob-nikolas.kather{at}alumni.dkfz.de Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Background Deep Learning (DL) has emerged as a powerful tool to predict genetic biomarkers directly from digitized Hematoxylin and Eosin (H&E) slides in colorectal cancer (CRC). However, few studies have systematically investigated the predictability of biomarkers beyond routinely available alterations such as microsatellite instability (MSI), and BRAF and KRAS mutations. Methods Our primary dataset comprised H&E slides of CRC tumors across five cohorts totaling 1,376 patients who underwent comprehensive panel sequencing, with an additional 536 patients from two public datasets for validation. We developed a DL model using a single transformer model to predict multiple genetic alterations directly from the slides. The model’s performance was compared against conventional single-target models, and potential confounders were analyzed. Findings The multi-target model was able to predict numerous biomarkers from pathology slides, matching and partly exceeding single-target transformers. The Area Under the Receiver Operating Characteristic curve (AUROC, mean ± std) on the primary external validation cohorts was: BRAF (0·78 ± 0·01), hypermutation (0·88 ± 0·01), MSI (0·93 ± 0·01), RNF43 (0·86 ± 0·01); this biomarker predictability was mirrored across metrics and co-occurrence analyses. However, biomarkers with high AUROCs largely correlated with MSI, with model predictions depending considerably on MSI-associated morphology upon pathological examination. Interpretation Our study demonstrates that multi-target transformers can predict the biomarker status for numerous genetic alterations in CRC directly from H&E slides. However, their pre-dictability is mainly associated with MSI phenotype, despite indications of slight biomarker-inherent contributions to a phenotype. Our findings underscore the need to analyze confounders in AI-based oncology biomarkers. To enable this, we developed a validated model applicable to other cancers and larger, diverse datasets. Funding The German Federal Ministry of Health, the Max-Eder-Programme of German Cancer Aid, the German Federal Ministry of Education and Research, the German Academic Exchange Service, and the EU. Introduction Exome sequencing, including targeted panel sequencing, plays a key role in precision oncology of colorectal cancer (CRC) 1 , 2 but remains inaccessible to many CRC patients worldwide due to the need for cost-intensive laboratory equipment and complex data analysis. 3 In contrast, Hematoxylin and Eosin (H&E) stained histopathology slides are a standard diagnostic tool available for almost every cancer patient globally. The advent of deep learning (DL) has unlocked these slides as a quantifiable data resource. Recent studies have extensively suggested that DL can predict molecular biomarkers directly from digitized H&E slides, including microsatellite instability (MSI) 4 – 10 , hypermutation status, and gene mutations such as TP53 , BRAF , and KRAS. 7 , 8 , 10 – 13 When applied as pre-screening tools, such DL systems can streamline the diagnostic workflow 9 , 14 , identifying those cases that need further testing and ruling out others. 15 Previous DL studies in CRC have predominantly focused on specific genetic alterations as potential biomarkers, referred to as ‘prediction targets’ in computational pathology, with only few studies adopting a pan-molecular alteration approach. 11 , 16 Such efforts have often been constrained by the challenge of obtaining diverse datasets with extensive sequencing, and even then, traditional methods require training separate models for individual prediction targets 5 – 8 , 12 , 13 , making them labor- and resource-intensive. Addressing both challenges, this study introduces the first pan-biomarker DL approach for CRC, utilizing a single model to predict multiple molecular targets. The transformer-based model was trained and validated on a unified dataset from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), which consolidates sequencing data from diverse cohorts. 17 Furthermore, its generalizability was evaluated using two extensively studied CRC cohorts. As DL has demonstrated robust capability in linking phenotype to genotype, particularly for pre-dicting MSI 5 – 9 , we leverage the strengths of our dataset and model to extend this analysis to multiple prediction targets. Specifically, we analyzed the co-occurrence of genetic alterations with MSI, evaluated the ability of DL to predict these alterations, and examined the corresponding slide morphology, relating it to features known to be relevant for MSI detection. Our study included genetic alterations investigated in prior DL-based studies and clinically relevant mutations, such as BMPR2 , RNF43 and BRAF as well as MSI and hypermutation status. Materials and Methods Patient Samples GECCO is an international collaboration utilizing genotype and sequencing data from a growing resource of over 60 studies. 17 Our experiments were conducted on five cohorts within GECCO, comprising N=1,376 patients with complete data and digitized whole-slide images (WSIs): European Prospective Investigation into Cancer (EPIC, N=183) 18 , Colorectal Cancer Study of Austria (CORSA, N=158) 19 , Iowa Women’s Health Study (IWHS, N=390) 20 , Cancer Risk Assessment study (CRA, N=321) and Women’s Health Initiative (WHI, N=324) 21 ( Fig. 1A-C ). These studies provided digitized WSIs with unified demographic, clinical, and lifestyle data (relevant data for this study in Tab. S1), establishing them as the primary dataset for this study. Centralized targeted tumor sequencing of up to 356 genes was conducted on all samples, with MSI and hypermutation status also assessed. The analysis targeted non-silent mutations and mutational signatures using panel sequencing of 1·8 Megabases (Mb), with tumor and normal sequencing coverage of 975× and 273×, respectively. To ensure generalizability of our model and comparability with existing literature, we utilized a secondary external dataset with two well-characterized publicly available CRC cohorts from The Cancer Genome Atlas (TCGA, N=426) 22 and the Clinical Proteomic Tumor Analysis Consortium (CPTAC, N=110, fresh-frozen tissue samples) 23 ( Fig. 1C ). We included patients with WSIs and ground truth data on the primary genetic alterations analyzed in this study. The MSI status for TCGA was defined per Liu et al. 24 , with high grade MSI categorized as MSI and low grade MSI and microsatellite stable (MSS) grouped as MSS 25 . Ground truth labels were derived from the clinical data accessible at https://portal.gdc.cancer.gov/ (accessed on Nov 20, 2024). Download figure Open in new tab Fig. 1: Experimental design, cohort characterization, and schematic for predictive analysis. A . Tissue samples from colorectal cancer (CRC) patients across five independent cohorts were obtained via surgical resection, with associated demographic, clinical, and sequencing data collected. Upon Hematoxylin and Eosin (H&E) staining, tumor tissues are digitized into Whole Slide Images (WSIs) for profiling genetic alterations. The WSIs are then used to train and test a deep learning (DL) algorithm for biomarker detection, to simultaneously predict multiple mutational statuses and provide heatmap explanations. B . The DL pipeline tessellates the WSIs into smaller tiles while rejecting background and blurry areas, extracting n feature vectors from n tiles. Feature vectors are compressed and processed in a multi-target transformer, employing an attention mechanism in an encoder-decoder structure for class token learning. The transformer generates individual scores for the respective amount of classes per target. The code is able to comprise positional tile embedding (dashed lines), which did not result in improved performance and were therefore excluded from our study. C . Overview of the five GECCO and two public cohorts, including patient numbers, slides, extracted features, and MSI case proportions. The cohorts are divided into train datasets and test datasets. D . Schematic for interpreting result plots and statistics, delineating dataset partitioning based on microsatellite (MSS: microsatellite stability, MSI: microsatellite instability) and gene mutational status (MUT: mutated, WT: wild type). The diagram illustrates distinct groups by color, with the left side representing MSI prediction scores and the right side for prediction target scores. True ground truth labels of samples guide the group organization, with model-generated scores depicted in corresponding colors. Experimental Design Our study included genetic alterations investigated in prior DL-based studies, clinically relevant mutations, and genetic alterations strongly associated with MSI in the TCGA CRC cohort 26 , including APC , BMPR2 , BRAF , KRAS , RNF43 , TP53 , ZNRF3 . 7 , 8 , 11 – 13 , 26 , 27 Together with MSI status and hypermutation, these exemplify the main prediction targets among the broad target set utilized from the GECCO cohorts (Tab. S2). Continuous prediction targets were discretized into classes using predefined thresholds to achieve balanced case distribution. Only prediction targets with at least 20 samples per class were included in model training to ensure robust analysis. Model training involved 731 patients from three GECCO cohorts (EPIC, CORSA, IWHS; Fig. 1C ) using seven-fold cross-validation to balance training and validation sets, ensure class distribution, and identify a median-performing model. The resulting seven models were deployed on external test sets. The primary test set consisted of 645 patients from two GECCO cohorts (CRA, WHI; Fig. 1C ), ensuring a diverse and representative split while addressing missing cases (Tab. S1). Exclusively female cohorts were included in both datasets (IWHS for training, WHI for testing). To assess generalizability and enable comparison to the literature, selected models were validated on a secondary test set comprising TCGA and CPTAC cohorts. ‘External validation’ and ‘test set’ typically refer to the primary test set unless stated otherwise; cohort-wise analyses are explicitly noted. In total, we trained eleven models using the same training set: a primary model incorporating the broad target set (Tab. S2) and a secondary model excluding MSI as a target to evaluate its impact on other predictions—both tested on the primary and secondary test sets. Additionally, we trained nine single-target models, each dedicated to one of the main prediction targets, and evaluated them on the primary test set. All models were evaluated in direct comparison to the primary multi-target model and compared to the literature. Data Analysis of Genetic Alterations In CRC, there is a pronounced co-occurrence of genetic alterations with MSI. 26 , 27 We employed two distinct methods to investigate the interplay and dependencies between various genetic alterations within the five GECCO cohorts: First, we used hierarchical clustering to group molecular alterations in our dataset using the ‘Euclidean’ metric and ‘Ward’ procedure 28 , using only alterations with complete data (Tab. S1). Subsequently, we applied association rule mining 29 to reveal the co-occurrence of gene mutations and conditions such as MSI, considering an initiating genetic alteration (’antecedent’) and a potentially resulting alteration (’consequent’). The presence of an initiating alteration statistically increases the probability of observing a resulting alteration, without implying a cause-effect biological sequence. We used six complementary metrics to assess co-occurrence: ‘Support’, ‘Confidence’, ‘Lift’, ‘Leverage’, ‘Conviction’, ‘Zhang’s Metric’. Image Processing and Deep Learning Techniques According to the framework of Wagner et al. 25 , the digitized WSIs were tessellated into tiles of 224×224 pixels, corresponding to 256×256 micrometers ( Fig. 1B ). Tiles predominantly containing background (brightness value ≥ 224) and blurred regions, identified using Canny edge detection 30 (thresholds: 40, 100) with fewer than two percent of the tile’s pixels identified as edges were excluded. We used the pre-trained CTransPath feature extractor 25 , 31 to extract a 768-dimensional feature representation for each tile, which we used as an input for a subsequent transformer model ( Fig. 1B ). Advancing upon previous work 25 we developed a model with an encoder-decoder architecture to simultaneously predict multiple targets from tile embeddings ( Fig. 1B ). Specifically, the 768-dimensional tile features are projected into a 512-dimensional space using a fully connected layer to reduce complexity and enhance computational efficiency for the subsequent model. The encoded tokens are decoded in class tokens 32 , each with a dimension of 1×512, with one token dedicated for each prediction target. The decoded tokens are fed through a fully connected layer to generate target-specific predictions with scores ranging from zero (negative prediction) to one (positive prediction) for each class and patient. To address class imbalance during training, cross-entropy loss is calculated per target, weighted by the inverse mutation 5 frequency, and summed to ensure proportional importance of rare mutations. The training and deployment procedures were performed on a NVIDIA RTX A6000 featuring 48 GB of GPU memory. All source codes for preprocessing and DL are available in our open source repositories (Tab. S3). TRIPOD reporting guidelines were used. Explainability To investigate the model’s detection of relevant regions for predicting alterations in CRC, we generated heatmaps using Grad-CAM 33 , focusing on the fold with the median AUROC for MSI detection to ensure consistent comparisons. These heatmaps highlight the contribution of each tile to patient-level predictions, enabling analysis of whether the model relies on distinct patterns and areas in WSIs for different targets or similar features across targets. Notably, a high number of positively contributing tiles does not necessarily correspond to a high final score due to the nonlinear aggregation process, and the model may also incorporate global cues not captured by the heatmap. To examine morphology in greater detail, we identified highly predictive tiles (‘top tiles’) for representative cases, selected based on their scores and attention values assigned by the model. For direct comparability, we included top tiles for each prediction target alongside those for MSI from the same slide. As an additional mode of explainability, we analyzed class token interactions in the decoder during deployment of the primary model to assess overlap for the main prediction targets. Grad-CAM was used to compute activations, capturing class token-score interactions, which were aggregated into a cross-correlation matrix, averaged across patients, and visualized as a heatmap. Statistical Analysis We used descriptive statistics to summarize sociodemographic and clinicopathologic characteristics across cohorts (Tab. S1). Model performance was evaluated using AUROC during training and testing, with additional metrics and target-specific prediction scores (ranging from 0·00 to 1·00) analyzed. For the main prediction targets we applied a two-sided DeLong test to compare AUROCs of single-target models and the primary multi-target model on the primary external dataset, utilizing mean prediction scores across seven folds per target. Mean and median AUROC values for each target were calculated over the seven folds. DeLong tests were also applied to compare the primary model (including MSI as a target) with the secondary model (excluding MSI) within each of the four external cohorts, based on mean prediction scores from the seven models. Given the rarity of mutations and the limitations of metrics such as AUROC and the Area Under the Precision-Recall Curve (AUPRC) in capturing imbalanced distributions 34 , we employed diverse metrics and analyzed individual prediction scores for the main prediction targets, stratified by MSI status. Average prediction scores across seven folds were calculated per patient and categorized into four subgroups per target based on microsatellite and mutational status: MSS/WT and MSS/MUT represent MSS tumors with wild-type (WT) and mutated (MUT) targets, respectively, while MSI/WT and MSI/MUT represent MSI tumors with wild-type and mutated targets ( Fig. 1D ). Model mutation detection performance was evaluated through comparisons on training and testing sets. After testing for normality with the Shapiro-Wilk test, the Mann-Whitney U test was used to assess statistical significance between MSI and prediction target scores within each subgroup. Additionally, four comparisons using the Wilcoxon test evaluated the model’s ability to identify prediction target mutations by comparing score distributions for MSS/WT vs MSS/MUT and MSI/WT vs MSI/MUT. Finally, subgroup comparisons based on prediction target status were performed: wild type (MSS/WT vs MSI/WT) or mutated (MSS/MUT vs MSI/MUT). Role of the funding source The funders of the study had no role in data collection, analysis, interpretation, writing of the manuscript and the decision to submit. Results Multi-Target Transformers Match and Partly Exceed Single-Target Models While Enabling Simultaneous Prediction of Genetic Alterations Our approach aimed to reproduce prior DL-studies predicting genetic alterations in CRC from H&E slides 5 – 9 , 11 – 13 , 35 , 36 , expanding the scope to a broader set of alterations using a multi-target transformer architecture ( Fig. 1B ). We compared the performance of our transformer for selected targets including the main prediction targets ( Fig. 2 ) with the external validation AUROCs reported in the literature for single-target models, which used varying model architectures. Download figure Open in new tab Fig. 2: Analysis of genetic alterations co-occurrence in CRC for GECCO cohorts. Hierarchical clustering analysis was conducted on the ground truth of genetic alterations with fully available mutational information. Each row corresponds to a genetic alteration, and each column represents a patient from the dataset. The top row indicates the distribution of patients from various cohorts within genetic clusters. The distance calculation was performed using the ‘Euclidean’ metric, and the ‘Ward’ method was applied to clustering. Three unique genetic clusters were created and marked. The patient clustering shows a diverse distribution of samples across all five cohorts and genetic clusters (top row). For MSI detection, single- and multi-target transformers achieved AUROCs of 0·91 (±0·02) and 0·93 (±0·01), respectively, on the primary test set (p<0·05; refer to Tab. S4, Fig. 3A here and in the following). The primary model achieved AUROCs from 0·87 (±0·01, TCGA) to 0·94 (±0·01, WHI) (Tab. S5-S8, Fig. S1; referred to here and in the following), consistent with the reported literature range of 0·77–0·96. 4 , 7 , 8 , 10 , 12 , 13 , 25 , 36 The following summarizes model performance for external validation for selected targets, with AUROCs for single- and multi-target transformers (Tab. S4, Fig. 3A ) and comparisons between primary (multi-target with MSI) and secondary (multi-target without MSI) models (Tab. S5-S12, Fig. S2). In BRAF MUT detection, the the AUROC was 0·72 (±0·06, single-target) and 0·78 (±0·01, multi-target, p0·05)) to highest on CRA (0·83 (±0·02) vs. 0·82 (±0·03, p>0·05)). These results align with the literature range of 0·66 to 0·88. 7 , 8 , 10 – 13 , 16 , 25 , 36 In RNF43 MUT detection, the AUROC was 0·80 (±0·05, single-target) and 0·86 (±0·01, multi-target, p0·05) on TCGA to 0·87 (±0·02) vs. 0·85 (±0·01, p>0·05) on CRA. These results exceed the literature range of 0·63 to 0·72. 11 , 16 For KRAS MUT detection, the AUROC was 0·65 (±0·02) for both single- and multi-target models (p>0·05). The primary vs. secondary models ranged from 0·56 (±0·02) vs. 0·55 (±0·02, TCGA, p>0·05) to 0·69 (±0·06) vs. 0·72 (±0·03, CPTAC, p>0·05). These results fall below and within the literature range of 0·60 to 0·80. 8 , 11 , 12 , 16 , 25 , 36 Interestingly, CPTAC had no KRAS MUT cases with MSI, while TCGA showed the highest number of such cases among all cohorts, with up to 19 instances. For hypermutation and TP53 MUT detection, the results fell within the literature ranges of 0·81–0·87 8 , 10 , 13 and 0·60–0·75 8 , 10 , 11 , 16 , with no significant differences between models. For APC mutational status prediction, the results ranged both above and below the single available literature AUROC of 0·67 11 , 16 , with no statistical differences between models, highlighting the limited number of comparable values. For BMPR2 and ZNRF43 there was no comparable data in the literature (AUROCs shown in Tab. S4, Fig. 3A , Fig. S2). Additional prediction targets with AUROCs ≥0·75 are detailed in Tab. S4. Download figure Open in new tab Fig. 3: Evaluation of the performance of the Multi-Target Transformer on selected prediction targets for the external cohorts from GECCO. A. The comparison of Single-Target Transformer versus Multi-Target Transformers shows the Area Under the Receiver Operating Characteristic curve (AUROC) from each of the 7 folds of external cross-validation, with the median value highlighted with a horizontal line in each box. The figure includes selected representative potential biomarkers of genetic alterations associated with MSS ( Fig. 2 , genetic Cluster 1) and MSI ( Fig. 2 , genetic Cluster 2). The test set cohorts consist of CRA and WHI ( Fig. 1C ). The horizontal line positioned at an AUROC of 0·50 represents a random guess of the model. Significance was determined through a two-sided DeLong test with a p-value threshold of less than 0·05. B. Performance metrics of Multi-Target Transformers for external validation. The mean (center of dot) and standard deviation (diameter of dot) for relevant selected prediction targets for the whole external set, as well as the MSI and MSS subgroups, are displayed based on the 7 folds of cross-validation. The threshold for binary classification is pre-defined as 0·50· The evaluation metrics include the Area Under the Receiver Operating Characteristic Curve (AUROC), and the Area Under the Precision-Recall Curve (AUPRC), along with the corresponding mutation rates in external cohorts. The Mutation Rate refers to the fraction of instances with a specific mutation in the subgroup. The MSI & Genetic Alteration Co-Occurrence Ratio is the fraction of cases harboring MSI among all cases with a particular genetic mutation. The data is sorted for AUROC and shown in Tab. S14 and Tab. S17-S18. An extended version of this panel with more metrics is shown in Fig. S3B. C. Distribution of Areas under the Receiver Operating Characteristic Curve (AUROCs, mean ± standard deviation) for selected prediction targets and their co-occurrence with MSI with corresponding values and further metrics shown in Tab. S14. An extended version of this panel with MSS/MSI-subgroup specific AUROCs is shown in Fig. S3C. Enhanced Predictive Performance for Genetic Alterations Associated with MSI Hierarchical clustering identified two primary genetic clusters: Cluster 2, comprising MSI-associated genes (e.g., BRAF , BMPR2 , ZNRF3 , RNF43 ) with strong co-occurrence with MSI and hypermutation, and Cluster 1, including MSS-associated genes such as TP53 , KRAS , and APC ( Fig. 2 ). Association rule mining reinforced robust MSI correlations in Cluster 2 (e.g., BMPR2 , RNF43 ), contrasting with inverse relationships in Cluster 1 MSS-linked genes (e.g., TP53 , KRAS ) (Fig. S3A, Tab. S13). MSI-associated alterations showed higher AUROC values for predicting genetic mutations compared to MSS-associated alterations (Cluster 1) ( Fig. 3B-C , Tab. S14, Fig. S1). Additional metrics ( Fig. 3B , Fig. S3B, Tab. S14) and patient-level prediction score distributions were analyzed for MSI and genetic alterations within Clusters 1 and 2 ( Fig. 1D , Fig. 4 , Tab. S15-S16). Prediction scores quantified model certainty, with high scores indicating mutation presence, low scores indicating absence, and scores near 0·5 reflecting indecision. Comprehensive summaries supported internal (Fig. S4) and external validation (Fig. S5). For MSS-associated Cluster 1 genes ( TP53 , APC , KRAS ) the external validation AUROCs ranged from 0·65 to 0·72, with MSI scores accurately reflecting ground truth: low for MSS and high for MSI cases ( Fig. 2 , Fig. 4A , Tab. S4). High MSI scores tended to align with WT target predictions, while low scores aligned with MUT predictions. As a result, discrepancies were observed in MSS/WT and MSI/MUT subgroups, where target predictions deviated from ground truth. Genetic Cluster 2, including MSI, hypermutation, and the genes BMPR2 , ZNRF3 , RNF43 , and BRAF , demonstrated higher AUROCs (0·75–0·88) during external validation ( Fig. 2 , Fig. 4B , Tab. S4). In MSS patients, mutations in biomarkers like BMPR2 (n=3) and ZNRF3 (n=9) were rare, and MSS/MUT subgroups showed less distinct prediction scores for RNF43 (n=20) and BRAF (n=39), reflecting higher prediction uncertainty. MSI scores aligned with the ground truth (MSS/MSI), but unlike Cluster 1, alteration scores in all subgroups correlated with MSI scores. This produced accurate prediction trends for MSS/WT and MSI/MUT subgroups (low scores for WT, high for MUT), while MSS/MUT and MSI/WT subgroups showed deviations from the target ground truth. These findings, partially reflected in AUROC results ( Fig. 3B , Fig. S3B), suggest less effective differentiation between MUT and WT in MSS and MSI subgroups compared to the combined group. Despite the strong influence of MSI-associated morphology on predictions, aligning the prediction target scores with MSS (Cluster 1) or MSI-high (Cluster 2), intrinsic phenotypes for some alterations cannot be fully confirmed or excluded. Analyses of MSI-high and MSS subgroups (Tab. S17-S18, Fig. 3B , Fig. S3B-S3C) revealed AUROCs of 0·60–0·70 for genes like BRAF , RNF43 , and TP53 in MSS patients, indicating modest discrimination between MUT and WT. Statistical differences in score distributions for MSI and these genetic alterations further supported these findings ( Fig. 4A-B ). Download figure Open in new tab Fig. 4: Evaluation of prediction scores based on the multi-target transformer in external validation on the GECCO test set subgrouped by the co-occurrence of the prediction targets with MSI. A.-B. Violin plots representing individual patient scores from the test set cohorts for MSI and representative genetic alterations in four subgroups based on microsatellite and alteration mutational status. The left y-axis represents the MSI score scale (left violin halfs) and the right y-axis corresponds to the prediction target scores (right violin halfs). The legend displays gray horizontal lines in the concept violins that represent the optimal position of the prediction scores based on ground truth. The selection of prediction targets includes TP53 , APC , and KRAS from genetic Cluster 1 ( A .), and BMPR2 , ZNRF3 , Hypermutation (HM), RNF43 , and BRAF from genetic Cluster 2 ( B .) ( Fig. 2 ). The data encompasses both external cohorts CRA and WHI ( Fig. 1C ). Each dot represents the mean value of individual patient prediction scores calculated from 7 folds, with the horizontal line on each side of the violin indicating the median of all individual mean patient scores. A horizontal line at 0·50 denotes the line of model uncertainty. The sample count for each subgroup is indicated below the violins. Statistical significance is denoted in the figures as follows: * for p < 0·05, ** for p < 0·01, *** for p < 0·001, with more details provided in Fig. 1D . After testing for normal distribution (Tab. S20), the Mann-Whitney U test was used for within-group comparisons, and the Wilcoxon test was used for between-group comparisons. Abbreviations: HM: Hypermutation; MSI: Microsatellite instability; MSS: Microsatellite stability; MUT: Mutated; WT: Wild type. MSI-Associated Morphological Patterns Drive Predictions Across Prediction Targets To investigate the morphological patterns underlying the prediction of genetic alterations in CRC from H&E slides using DL, we analyzed WSI heatmaps ( Fig. 5 , Fig. S6-S10) and highly predictive tiles (‘top tiles’, Fig. 6 , Fig. S11-S23) for key prediction targets (Cluster 1: TP53 , APC , KRAS ; Cluster 2: RNF43 , BRAF , hypermutated) in 25 cases from CRA and WHI test cohorts. This was assessed via manual histopathological review of WSI heatmaps and 20 top tiles per target (Tab. S19). We found that our model predominantly focused on tumor regions, with minimal relevance attributed to pen marks or non-tumor areas (e.g., Fig. 5A-D , Fig. S6B-D, Fig. S7A-D). Pen marks, present on most slides, were rarely highlighted and appeared faintly in top tiles (e.g., Fig. S13, Fig. S22), though slight highlighting in heatmaps occurred in rare cases (e.g., Fig. 5C , Fig. S9D). Regions of high model attention primarily consisted of tumor tissue, not normal intestinal mucosa or uninvolved connective/adipose tissue, despite the model not being explicitly trained for tumor detection. Download figure Open in new tab Fig. 5: Heatmaps of representative samples for prediction of MSI, KRAS, BRAF and Hypermutation (HM) from the external GECCO validation dataset. The heatmaps are derived from the model with the median AUROC for MSI detection and the majority of prediction targets evaluated by sevenfold cross-validation. The cohort, Sample-ID, ground truth and prediction scores for MSI, along with the individual mutational status of the target, a brief pathological evaluation and magnified views of specific areas are provided for in-depth analysis. The heatmaps indicate relevant areas for the various predictions. The red areas are of high importance and indicate a mutant type (MUT), while the blue areas are of low importance and indicate a wild type (WT). The color intensity showcases the model’s attention to that distinct area. A. The tumor exhibits both gland-forming and more solid components and extremely high numbers of tumor-infiltrating lymphocytes (TILs) with dense lymphoid aggregates. The pathological examination confirms the plausibility of a high MSI score indicating MSI which is also the ground truth. A low KRAS score indicates KRAS WT but the ground truth is KRAS MUT. The heatmap highlights similar tumor areas but with diverging scores: where MSI map is red indicating high score, KRAS map is blue indicating low score. B. The presence of mucinous differentiation in a MSI, BRAF WT case results in high MSI and BRAF scores. The MSI score is pathologically plausible whereas the BRAF score indicates a contrary prediction tendency than the ground truth holds. For both predictions, the model focuses on similar tumor areas with similar scores indicating MSI/ BRAF MUT. C. Partly mucinous morphology indicates the possibility of MSI, with a high score predicting MSI. HM is also predicted MUT with a high score, even though HM is WT for this sample. Both heatmaps primarily label the tumor and the same region with comparable significance. D. Villous adenoma with high grade dysplasia is a common precursor lesion associated with high frequency of KRAS mutations 41 , 42 The heatmaps highlight similar large scale tumor areas but with converging scores: where the MSI map is red indicating a high score, the KRAS map is blue indicating low score. E . The tumor area appears to be mainly MSS, and the heatmap predicts a low score, indicating this. Although it is being mutated in the ground truth, it is still predicted as non-hypermutated. This is a rare MSS case with HM 54 Both heatmaps predominantly mark the tumor area and the same region with comparable relevance. Abbreviations: HM: Hypermutation; MSI: Microsatellite instability; MSS: Microsatellite stability; MUT: Mutated; WT: Wild type; w/: with Download figure Open in new tab Fig. 6: Top tiles for prediction of genetic alterations (left column) and MSI (right column) for two selected slides from the GECCO test set. A.-B. WHI, 1031792: Medullary carcinoma with sheets of tumor cells, low stroma content, high number of tumor-infiltrating lymphocytes in a KRAS WT and MSI case, leading to high MSI prediction scores ( B. ) and low KRAS MUT prediction scores ( A. ). KRAS mutations show a lower frequency in MSI CRCs. Medullary carcinoma is a key morphological feature of MSI CRCs. C.-D. WHI 1031557: Top tiles for BRAF MUT as well as MSI prediction, both predictions with high prediction scores, both displaying a mixed morphology with partly medullary, partly mucinous, partly gland-forming histology, and high number of tumorinfiltrating/associated lymphocytes. Medullary growth pattern with lymphocytic infiltration and mucinous differentiation are typical features of a MSI-like morphology. Accordingly, the case was correctly predicted as MSI with a really high prediction score ( D. ). As BRAF MUT and MSI often co-occur and share morphologic overlap, the case was misclassified with regards to BRAF -status, resulting in high prediction scores for BRAF MUT ( C. ), even though the ground truth was BRAF WT. Abbreviations: CRC: Colorectal cancer; MSI: Microsatellite instability; MSS: Microsatellite stability; MUT: Mutated; WT: Wild type. Frequent MSI-associated patterns in most MSI cases and Cluster 2 gene mutations included medullary growth, high tumor-infiltrating lymphocytes (TILs), and mucinous differentiation, consistent with prior studies. 37 – 39 However, some MSI cases in our cohort and the literature display atypical morphologies not commonly linked to MSI (e.g., Fig. S8D, Fig. S9D). 40 In MSS-related Cluster 1, KRAS mutation predictions highlighted luminal tumor regions, especially villous adenomas with high-grade dysplasia, as equally important as invasive adenocarcinoma areas (Fig. S16A-B, Fig. S17A-B). This aligns with the known prevalence of KRAS mutations in villous adenomas 41 and carcinomas with close proximity to polyps 42 , reflecting the algorithm’s ability to capture intratumoral heterogeneity. Tiles with high MSI relevance displayed medullary carcinoma patterns, including tumor cell sheets and high TILs, which corresponded to low KRAS MUT prediction scores ( Fig. 6A-B ). Alterations such as TP53 and APC , less common in MSI-high cases, relied on MSS-like morphology (e.g., gland-forming conventional adenocarcinoma with dirty necrosis 40 ; Fig. S18-S20), while MSI-associated alterations (e.g., BRAF , RNF43 , hypermutation) depended on MSI-related features, including medullary patterns, mucinous differentiation (with signet-ring cells; Fig. S23B), and high TILs 40 ( Fig. 6C-D , Fig. S12-S15, Fig. S21, Tab. S19). To explore target-specific patterns beyond MSI, a pathological review of top tiles from selected slides (Fig. S22-S23) identified tumor budding as a morphological feature associated with BRAF mutations in MSS cases (Fig. S23A-C). 43 Discussion Previous research in CRC has focused on high AUROCs for biomarker detection, showing MSI as the most detectable, while often neglecting co-occurrence with other mutations and lacking indepth analyses of multiple genetic alterations and their links to MSI-associated morphologies. To address these limitations, we developed a DL model to predict multiple genetic alterations in CRC from H&E slides, using five GECCO cohorts with unified sequencing data with additional validation on two public CRC-datasets, TCGA and CPTAC. The multi-target model achieved performance within the literature range, outperforming single-target models for some alterations while enabling the efficient and sustainable simultaneous prediction of numerous prediction targets. This in-depth analysis was made possible by the comprehensive dataset and novel model architecture, both of which are shared open-source with the scientific community. Our cohort, comprising multiple sub-cohorts, is broadly representative of CRC patients 44 , particularly in Western populations, with data primarily from the US and Europe. It includes a higher proportion of female patients (IWHS, WHI) and disproportionately includes White individuals, while TCGA had the highest proportion of Black/African American patients, and CPTAC was limited to fresh-frozen colon samples with the smallest cohort. MSI frequency ranged from 7–35%, aligning with reported rates of 15–20% 45 (Tab. S1). Variability in target mutation frequencies across cohorts, remained within plausible ranges, as did the experimental results. However, small sample sizes for some prediction targets likely contributed to higher variability (Fig. S2), underscoring the need for larger, more diverse datasets. The genetic clusters analyzed included key CRC biomarkers (e.g., TP53 , APC , KRAS , BRAF ) 7 , 8 , 11 – 13 , with mutational profiles and MSI co-occurrences aligning with prior studies 26 , 27 , 46 . In line with the literature MSI was the most predictable alteration 7 , 8 , 12 , 13 , 25 , 36 , though some MSI cases scored low due to a lack of associated morphological features, which is plausible and pathological assessment has shown. For the detection of multiple genetic alterations apart from MSI, our multi-target model outperformed several of our own single-target models and models from the literature, such as for detection of RNF43 , using the AUROC as a common performance measure. 4 – 9 , 35 , 36 Notably, hypermutation prediction demonstrates a strong interdependence with 10 MSI due to their frequent co-occurrence and close biological linkage, with shared morphological features leading to lower scores for hypermutated MSS cases. BRAF mutations, critical in CRC for their prognostic and predictive significance 47 , showed detectable phenotypic changes, including mucinous differentiation and poorly differentiated clusters. The histopathological evaluation revealed features such as mucinous differentiation and stroma-rich patterns dominating BRAF MUT and BRAF WT tiles, respectively 12 , 47 Interestingly, we could identify mucinous/signet-ring differentiation relevant for BRAF predictions. 47 However, these distinct features overlapped with MSI-associated morphology, confounding specificity and occasionally leading to misclassification in MSS cases. In contrast, some top tiles for BRAF MUT prediction showed poorly differentiated clusters and tumor budding, which has already been linked to BRAF mutations 48 , and is associated with MSS. 43 While mutations like BRAF induce distinct phenotypic changes detectable by DL, MSI morphology often overshadows fine-grained, target-specific patterns. Excluding MSI as a prediction target did not affect the predictability of other prediction targets, and independent class token behavior with no cross-token information flow in the decoder, a key feature of the model architecture, was demonstrated (Fig. S24). In summary, we have shown that CRC mutation predictability is primarily driven by MSI-associated morphology, the most discriminative feature for mutation detection, though subgroup analyses of MSI and MSS suggest subtle mutation-specific patterns warrant further study. This study has several limitations. Despite its extent, the dataset was limiting the detection of rare alterations and subtle mutation-specific morphologies, highlighting the need for larger and more diverse datasets. Underrepresentation of Black/African-American individuals in the GECCO dataset and missing or incomplete annotations further constrain generalizability and accuracy. While monochrome tiles were excluded during preprocessing, some pen markings remained. These markings, though occasionally attracting slight model attention and potentially introducing bias, were largely transparent, limited to a few tiles, and may hold biologically relevant information, particularly near tumor borders, warranting further investigation. Additionally, our study emphasizes the need for advanced explainability methods to validate DL-detected pathological patterns and link them to specific prediction targets, fostering a bidirectional exchange between pathology and DL. This could clarify whether MSI-associated morphology dominates model pre-dictions due to the absence of other target-specific features or its overwhelming influence. Approaches such as modified loss functions could help mitigate MSI dominance and enhance detection of subtle features 49 and integrating multimodal data could improve generalizability, accuracy, and understanding of mutation-specific morphologies. Our model expands on single-target transformers and raises important questions regarding potential bias toward MSI-associated targets due to their larger subgroup size; however, as our model performs within the range reported in the literature, further investigation is required to elucidate the influence of MSI- and non-MSI-associated features, providing insights into the mechanisms of single- and multi-target transformers. Similar to prior research 34 , 50 , we show that AUROC can be a misleading metric for evaluating biomarker detection; thus, we report it for comparability but include additional metrics. For future studies, we recommend using extensive, diverse datasets, analyzing multiple prediction targets and their interactions with various metrics (e.g. Accuracy, Precision, Sensitivity, Specificity, F1 Score), and addressing potential confounding factors by leveraging approaches such as our proposed model architecture. From a practical point of view, our multitarget DL-based biomarker prediction provides significant value for research by identifying which prediction targets can be reliably predicted and the morphological basis for these predictions, guiding future studies and highlighting prediction targets that merit further investigation. For some alterations, such as BRAF , existing surrogate immunohistochemical assays could be applied 51 to predicted regions in follow-up studies to deepen understanding of the biological mechanisms behind predictions. Clinically, this approach streamlines colorectal cancer testing by eliminating the need for separate workflows and serving as a costeffective prescreening tool, particularly in resource-limited settings or earlier CRC stages. Biologically, it sheds light on how genetic alterations influence tumor morphology and underscores the dominant role of surrogate markers like MSI in CRC histology. As the comprehensive analysis of co-occurrences in CRC is scarcely explored in relation to DL and is rare for other cancer types 52 , this efficient, scalable method can be extended to any cancer type, enabling simultaneous analysis of multiple prediction targets and advancing both research and precision oncology. In conclusion, our study highlights the utility of multi-target transformers for detecting biomarker-specific patterns in CRC H&E images, with MSI phenotype as the dominant factor influencing predictability. The model’s reliance on MSI- and MSS-associated morphology underscores the importance of morphology and co-occurrence patterns over AUROC metrics. By enabling the efficient investigation of numerous prediction targets simultaneously, this model facilitates a more comprehensive understanding of biomarker interactions. These findings emphasize the need to consider individual prediction scores and tumor-specific confounders in biomarker studies, laying the foundation for future research to address these effects across other cancer types. Declarations Ethics approval and consent to participate This study was performed in accordance with the Declaration of Helsinki. This study is a retrospective analysis of scanned images of anonymized tissue samples of various cohorts of cancer patients. Data were collected and anonymized and ethical approval was obtained. The overall analysis was approved by the Ethics board of the Medical Faculty of Technical University Dresden under the ID BO-EK-444102022· Contributions MG and JNK conceptualized the study. CS, PH, LAB, AJF, ELG, AG, MJG, PL, NM, ST and UP provided clinical and digital histopathological data. MG curated the source data and conducted the literature research. MvT implemented the source code for the deep learning pipeline and the heatmaps. MG developed the code for data analysis and visualization. MG planned and conducted the experiments, visualizations and analysis of data and results. ARM assisted with the statistical evaluation. MG, CMLL, ZIC and NGR interpreted the data. NGR and BM did the pathological examination of the samples and analyzed the heatmaps and top/bottom tiles. MG wrote the first draft of the manuscript. All authors revised the manuscript draft, jointly interpreted the data and agreed to the submission of this article. All authors had access to all the data, and MG and JNK have verified the data. Declaration of interest JNK declares consulting services for Bioptimus, France; Owkin, France; DoMore Diagnostics, Norway; Panakeia, UK; AstraZeneca, UK; Mindpeak, Germany; and MultiplexDx, Slovakia. Furthermore, he holds shares in StratifAI GmbH, Germany, Synagen GmbH, Germany; has received a research grant by GSK; and has received honoraria by AstraZeneca, Bayer, Daiichi Sankyo, Eisai, Janssen, Merck, MSD, BMS, Roche, Pfizer, and Fresenius. MG has received honoraria for lectures sponsored by Techniker Krankenkasse (TK) and AstraZeneca. SF has received honoraria for lectures by BMS and MSD. UP declares consulting services for AbbVie and her husband is holding individual stocks for the following companies: BioNTech SE – ADR, Amazon, CureVac BV, NanoString Technologies, Google/Alphabet Inc Class C, NVIDIA Corp, Microsoft Corp.. No other potential conflicts of interest are reported by any of the authors. Data sharing All source code we used to conduct this study is publicly available with publication. The code for image preprocessing is publicly available at https://github.com/KatherLab . The tesselation script is available at https://github.com/KatherLab/preprocessing-ng . Extraction of CTransPath features was conducted with scripts from https://github.com/KatherLab/marugoto . The scripts for the multi-target transformer and the heatmaps used for explainability can be accessed at https://github-com/LocalToasty/barspoon-transformer/ . The saved checkpoints for the trained models with the lowest validation loss for all seven folds can be obtained from https://github.com/gustavmarco/barspoon-transformer/releases/tag/gustav2024 · Links to the exact repository versions used can be found in Tab. S3. The digitized whole slide images (WSIs) for the TCGA 22 cohort are publicly accessible at https://portal.gdc.cancer.gov/ and https://www.cbioportal.org/ . The digitized WSIs for the CPTAC cohort are publicly accessible at The Cancer Imaging Archive (TCIA) 53 . The molecular and clinical data for TCGA and CPTAC is publicly accessible at https://portal.gdc.cancer.gov/ and https://www.cbioportal.org/ . The datasets from the GECCO consortium are available from the corresponding author on reasonable request. All data generated or analyzed during this study are included in this published article and its supplementary information files. Acknowledgements JNK is supported by the German Federal Ministry of Health (DEEP LIVER, ZMVI1-2520DAT111), the Max-Eder-Programme of German Cancer Aid (grant #70113864), the German Federal Ministry of Education and Research (PEARL, 01KD2104C; CAMINO, 01EO2101; SWAG, 01KD2215A; TRANSFORM LIVER, 031L0312A; TANGERINE, 01KT2302 through ERA-NET Transcan), the German Academic Exchange Service (SECAI, 57616814), the German Federal Joint Committee (Transplant.KI, 01VSF21048) the European Union’s Horizon Europe and innovation programme (ODELIA, 101057091; GENIAL, 101096312) and the National Institute for Health and Care Research (NIHR, NIHR213331) Leeds Biomedical Research Centre. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. SF is supported by the German Federal Ministry of Education and Research (SWAG, 01KD2215C), the German Cancer Aid (DECADE, 70115166 and TargHet, 70115995) and the German Research Foundation (504101714). The Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) is funded by: National Cancer Institute, National Institutes of Health, U.S. Department of Health and Human Services (U01 CA137088, R01 CA488857, P20 CA252733). Genotyping/Sequencing services were provided by the Center for Inherited Disease Research (CIDR) contract number HHSN268201700006I. This research was funded in part through the NIH/NCI Cancer Center Support Grant P30 CA015704· Scientific Computing Infrastructure at Fred Hutch funded by ORIP grant S10OD028685· The CORSA study was funded by Austrian Research Funding Agency (FFG) BRIDGE (grant 829675, to Andrea Gsur), the “Herzfelder’sche Familienstiftung” (grant to Andrea Gsur) and was supported by COST Action BM1206· CRA was supported by the National Institutes of Health grant R01 CA068535· The coordination of EPIC is financially supported by the International Agency for Research on Cancer (IARC) and also by the Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London which has additional infrastructure support provided by the NIHR Imperial Biomedical Research Centre (BRC). The national cohorts are supported by: Danish Cancer Society (Denmark); Ligue Contre le Cancer, Institut Gustave Roussy, Mutuelle Générale de l’Education Nationale, Institut National de la Santé et de la Recherche Médicale (INSERM) (France); German Cancer Aid, German Cancer Research Center (DKFZ), German Institute of Human Nutrition Potsdam-Rehbruecke (DIfE), Federal Ministry of Education and Research (BMBF) (Germany); Associazione Italiana per la Ricerca sul Cancro-AIRC-Italy, Compagnia di SanPaolo and National Research Council (Italy); Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); Health Research Fund (FIS) - Instituto de Salud Carlos III (ISCIII), Regional Governments of Andalucía, Asturias, Basque Country, Murcia and Navarra, and the Catalan Institute of Oncology - ICO (Spain); Swedish Cancer Society, Swedish Research Council and and Region Skåne and Region Västerbotten (Sweden); Cancer Research UK (14136 to EPIC-Norfolk; C8221/A29017 to EPIC-Oxford), Medical Research Council (1000143 to EPIC-Norfolk; MR/M012190/1 to EPIC-Oxford) (United Kingdom). The IWHS study was supported by NIH grants CA107333 (R01 grant awarded to P.J. Limburg) and HHSN261201000032C (N01 contract awarded to the University of Iowa). The WHI program is funded by the National Heart, Lung, and Blood Institute, National Institutes of Health, U.S. Department of Health and Human Services through contracts 75N92021D00001, 75N92021D00002, 75N92021D00003, 75N92021D00004, 75N92021D00005. We kindly thank all individuals who agreed to participate in the CORSA study. Furthermore, we thank all cooperating physicians and students and the Biobank Graz of the Medical University of Graz. We also acknowledge the TCGA Research Network and the Clinical Proteomic Tumor Analysis Consortium (CPTAC), which generated the data on which some of the results shown in this study are based. Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy or views of the International Agency for Research on Cancer/World Health Organization. The authors thank the WHI investigators and staff for their dedication, and the study participants for making the program possible. A full listing of WHI investigators can be found at: http://www.whi.org/researchers/Documents%20%20Write%20a%20Paper/WHI%20Investigator%20Short%20List.pdf . The readability and language of the work were improved using ChatGPT-4o and DeepL. List of abbreviations AUROC Area under the Receiver Operating Characteristic Curve AUPRC Area Under the Precision-Recall Curve CIN Chromosomal instability CPTAC Clinical Proteomic Tumor Analysis Consortium CRC Colorectal cancer DL Deep Learning H&E Hematoxylin and eosin Mb Megabases MSI Microsatellite instability MSS Microsatellite stable MUT Mutated NOS Not otherwise specified px Pixel ROC Receiver Operating Characteristic Curve TCGA The Cancer Genome Atlas TILs Tumor infiltrating lymphocytes ViT Vision Transformer WSI Whole Slide Image WT Wild type References 1. ↵ Tsimberidou AM , Fountzilas E , Nikanjam M , Kurzrock R . Review of precision cancer medicine: Evolution of the treatment paradigm . Cancer Treat Rev . 2020 Jun ; 86 : 102019 . 2. ↵ Xiao W , Ren L , Chen Z , Fang LT , Zhao Y , Lack J , et al. Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing . Nat Biotechnol . 2021 Sep ; 39 ( 9 ): 1141 – 50 . OpenUrl CrossRef PubMed 3. ↵ Phillips KA , Douglas MP , Wordsworth S , Buchanan J , Marshall DA . Availability and funding of clinical genomic sequencing globally . BMJ Glob Health [Internet ]. 2021 Feb ; 6 ( 2 ). Available from : doi: 10.1136/bmjgh-2020-004415 OpenUrl Abstract / FREE Full Text 4. ↵ Kather JN , Pearson AT , Halama N , Jäger D , Krause J , Loosen SH , et al. Deep learning can predict microsatellite instability directly from histology in gastrointestinal cancer . Nat Med . 2019 Jul ; 25 ( 7 ): 1054 – 6 . OpenUrl CrossRef PubMed 5. ↵ Echle A , Grabsch HI , Quirke P , van den Brandt PA , West NP , Hutchins GGA , et al. Clinical-Grade Detection of Microsatellite Instability in Colorectal Tumors by Deep Learning . Gastroenterology . 2020 Oct ; 159 ( 4 ): 1406 – 16.e11 . OpenUrl PubMed 6. Lee SH , Song IH , Jang HJ . Feasibility of deep learning-based fully automated classification of microsatellite instability in tissue slides of colorectal cancer . Int J Cancer . 2021 Aug 1 ; 149 ( 3 ): 728 – 40 . OpenUrl CrossRef PubMed 7. ↵ Saldanha OL , Quirke P , West NP , James JA , Loughrey MB , Grabsch HI , et al. Swarm learning for decentralized artificial intelligence in cancer histopathology . Nat Med . 2022 Jun ; 28 ( 6 ): 1232 – 9 . OpenUrl PubMed 8. ↵ Bilal M , Raza SEA , Azam A , Graham S , Ilyas M , Cree IA , et al. Development and validation of a weakly supervised deep learning framework to predict the status of molecular pathways and key mutations in colorectal cancer from routine histology images: a retrospective study . Lancet Digit Health . 2021 Dec ; 3 ( 12 ): e763 – 72 . OpenUrl 9. ↵ Saillard C , Dubois R , Tchita O , Loiseau N , Garcia T , Adriansen A , et al. Validation of MSIntuit as an AI-based pre-screening tool for MSI detection from colorectal cancer histology slides . Nat Commun . 2023 Nov 6 ; 14 ( 1 ): 6695 . OpenUrl PubMed 10. ↵ Zamanitajeddin N , Jahanifar M , Bilal M , Eastwood M , Rajpoot N . Social network analysis of cell networks improves deep learning for prediction of molecular pathways and key mutations in colorectal cancer . Med Image Anal . 2024 Apr ; 93 ( 103071 ): 103071 . OpenUrl PubMed 11. ↵ Kather JN , Heij LR , Grabsch HI , Loeffler C , Echle A , Muti HS , et al. Pan-cancer imagebased detection of clinically actionable genetic alterations . Nat Cancer . 2020 Aug ; 1 ( 8 ): 789 – 99 . OpenUrl PubMed 12. ↵ Niehues JM , Quirke P , West NP , Grabsch HI , van Treeck M , Schirris Y , et al. Generalizable biomarker prediction from cancer pathology slides with self-supervised deep learning: A retrospective multi-centric study . Cell Rep Med . 2023 Apr 18 ; 4 ( 4 ): 100980 . OpenUrl PubMed 13. ↵ Guo B , Li X , Yang M , Jonnagaddala J , Zhang H , Xu XS . Predicting microsatellite instability and key biomarkers in colorectal cancer from H&E-stained images: achieving state-of-the-art predictive performance with fewer data using Swin Transformer . Hip Int . 2023 May ; 9 ( 3 ): 223 – 35 . OpenUrl 14. ↵ Bilal M , Tsang YW , Ali M , Graham S , Hero E , Wahab N , et al. Development and validation of artificial intelligence-based prescreening of large-bowel biopsies taken in the UK and Portugal: a retrospective cohort study . Lancet Digit Health . 2023 Nov 1 ; 5 ( 11 ): e786 – 97 . OpenUrl PubMed 15. ↵ Kacew AJ , Strohbehn GW , Saulsberry L , Laiteerapong N , Cipriani NA , Kather JN , et al. Artificial Intelligence Can Cut Costs While Maintaining Accuracy in Colorectal Cancer Genotyping [Internet] . Vol. 11 , Frontiers in Oncology. 2021 . Available from : doi: 10.3389/fonc.2021.630953 OpenUrl CrossRef 16. ↵ Saldanha OL , Loeffler CML , Niehues JM , van Treeck M , Seraphin TP , Hewitt KJ , et al. Self-supervised attention-based deep learning for pan-cancer mutation prediction from histopathology . NPJ Precis Oncol . 2023 Mar 28 ; 7 ( 1 ): 35 . OpenUrl PubMed 17. ↵ Common Variant GWAS , Genetics & Epidemiology of Colorectal Cancer Consortium (GECCO) [Internet] . [cited 2023 Sep 28 ]. Available from: https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs001078.v1.p1 18. ↵ Riboli E , Kaaks R . The EPIC Project: rationale and study design. European Prospective Investigation into Cancer and Nutrition . Int J Epidemiol . 1997 ; 26 Suppl 1 : S6 – 14 . OpenUrl CrossRef PubMed Web of Science 19. ↵ Gsur A , Baierl A , Brezina S . Colorectal Cancer Study of Austria (CORSA): A Population-Based Multicenter Study . Biology [Internet ]. 2021 Jul 28 ; 10 ( 8 ). Available from : doi: 10.3390/biology10080722 OpenUrl CrossRef 20. ↵ Folsom AR , Kushi LH , Anderson KE , Mink PJ , Olson JE , Hong CP , et al. Associations of general and abdominal obesity with multiple health outcomes in older women: the Iowa Women’s Health Study . Arch Intern Med . 2000 Jul 24 ; 160 ( 14 ): 2117 – 28 . OpenUrl CrossRef PubMed Web of Science 21. ↵ Hays J , Hunt JR , Hubbell FA , Anderson GL , Limacher M , Allen C , et al. The Women’s Health Initiative recruitment methods and results . Ann Epidemiol . 2003 Oct ; 13 ( 9 Suppl): S18 – 77 . OpenUrl CrossRef PubMed Web of Science 22. ↵ Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer . Nature . 2012 Jul 18 ; 487 ( 7407 ): 330 – 7 . OpenUrl CrossRef PubMed Web of Science 23. ↵ Edwards NJ , Oberti M , Thangudu RR , Cai S , McGarvey PB , Jacob S , et al. The CPTAC Data Portal: A resource for cancer proteomics research . J Proteome Res . 2015 Jun 5 ; 14 ( 6 ): 2707 – 13 . OpenUrl CrossRef PubMed 24. ↵ Liu Y , Sethi NS , Hinoue T , Schneider BG , Cherniack AD , Sanchez-Vega F , et al. Comparative Molecular Analysis of Gastrointestinal Adenocarcinomas . Cancer Cell . 2018 Apr 9 ; 33 ( 4 ): 721 – 35.e8 . OpenUrl CrossRef PubMed 25. ↵ Wagner SJ , Reisenbüchler D , West NP , Niehues JM , Zhu J , Foersch S , et al. Transformer-based biomarker prediction from colorectal cancer histology: A large-scale multicentric study . Cancer Cell . 2023 Sep 11 ; 41 ( 9 ): 1650 – 61.e4 . OpenUrl CrossRef PubMed 26. ↵ Zhuang Y , Wang H , Jiang D , Li Y , Feng L , Tian C , et al. Multi gene mutation signatures in colorectal cancer patients: predict for the diagnosis, pathological classification, staging and prognosis . BMC Cancer . 2021 Apr 9 ; 21 ( 1 ): 380 . OpenUrl PubMed 27. ↵ Bao X , Zhang H , Wu W , Cheng S , Dai X , Zhu X , et al. Analysis of the molecular nature associated with microsatellite status in colon cancer identifies clinical implications for immunotherapy . J Immunother Cancer [Internet ]. 2020 Oct ; 8 ( 2 ). Available from : doi: 10.1136/jitc-2020-001437 OpenUrl Abstract / FREE Full Text 28. ↵ Murtagh F , Legendre P. Ward’s Hierarchical Clustering Method: Clustering Criterion and Agglomerative Algorithm [Internet] . arXiv [stat.ML]. 2011 . Available from: http://arxiv.org/abs/1111.6285 29. ↵ Han J , Kamber M , Pei J. Data Mining: Concepts and Techniques . 3rd ed . Oxford, England : Morgan Kaufmann ; 2011 . 744 p. (The Morgan Kaufmann Series in Data Management Systems). 30. ↵ Canny J . A computational approach to edge detection . IEEE Trans Pattern Anal Mach Intell . 1986 Jun ; 8 ( 6 ): 679 – 98 . OpenUrl CrossRef PubMed Web of Science 31. ↵ Wang X , Yang S , Zhang J , Wang M , Zhang J , Yang W , et al. Transformer-based unsupervised contrastive learning for histopathological image classification . Med Image Anal . 2022 Oct ; 81 : 102559 . 32. ↵ Dosovitskiy A , Beyer L , Kolesnikov A , Weissenborn D , Zhai X , Unterthiner T , et al. An Image is Worth 16×16 Words: Transformers for Image Recognition at Scale [Internet] . arXiv . 2020 . Available from: http://arxiv.org/abs/2010.11929 33. ↵ Selvaraju RR , Cogswell M , Das A , Vedantam R , Parikh D , Batra D. Grad-CAM: Visual Explanations from Deep Networks via Gradient-based Localization [Internet] . arXiv [cs.CV]. 2016 . Available from: http://arxiv.org/abs/1610.02391 34. ↵ McDermott MBA , Hansen LH , Zhang H , Angelotti G , Gallifant J. A Closer Look at AUROC and AUPRC under Class Imbalance [Internet] . arXiv [cs.LG]. 2024 . Available from: http://arxiv.org/abs/2401.06091 35. ↵ Yamashita R , Long J , Longacre T , Peng L , Berry G , Martin B , et al. Deep learning model for the prediction of microsatellite instability in colorectal cancer: a diagnostic study . Lancet Oncol . 2021 Jan ; 22 ( 1 ): 132 – 41 . OpenUrl CrossRef PubMed 36. ↵ Schrammen PL , Ghaffari Laleh N , Echle A , Truhn D , Schulz V , Brinker TJ , et al. Weakly supervised annotation-free cancer detection and prediction of genotype in routine histopathology . J Pathol . 2022 Jan ; 256 ( 1 ): 50 – 60 . OpenUrl PubMed 37. ↵ Tuppurainen K , Mäkinen JM , Junttila O , Liakka A , Kyllönen AP , Tuominen H , et al. Morphology and microsatellite instability in sporadic serrated and non-serrated colorectal cancer . J Pathol . 2005 Nov ; 207 ( 3 ): 285 – 94 . OpenUrl CrossRef PubMed Web of Science 38. Greenson JK , Huang SC , Herron C , Moreno V , Bonner JD , Tomsho LP , et al. Pathologic predictors of microsatellite instability in colorectal cancer . Am J Surg Pathol . 2009 Jan ; 33 ( 1 ): 126 – 33 . OpenUrl CrossRef PubMed 39. ↵ Malik A , Bhatia JK , Sahai K , Boruah D , Sharma A . Evaluating morphological features for predicting microsatellite instability status in colorectal cancer . Armed Forces Med J India . 2022 Sep ; 78 ( Suppl 1 ): S96 – 104 . OpenUrl 40. ↵ Shia J , Schultz N , Kuk D , Vakiani E , Middha S , Segal NH , et al. Morphological characterization of colorectal cancers in The Cancer Genome Atlas reveals distinct morphology-molecular associations: clinical and biological implications . Mod Pathol . 2017 Apr ; 30 ( 4 ): 599 – 609 . OpenUrl CrossRef PubMed 41. ↵ Zauber P , Marotta S , Sabbath-Solitare M . KRAS gene mutations are more common in colorectal villous adenomas and in situ carcinomas than in carcinomas . Int J Mol Epidemiol Genet . 2013 Mar 18 ; 4 ( 1 ): 1 – 10 . OpenUrl PubMed 42. ↵ Rosty C , Young JP , Walsh MD , Clendenning M , Walters RJ , Pearson S , et al. Colorectal carcinomas with KRAS mutation are associated with distinctive morphological and molecular features . Mod Pathol . 2013 Jun ; 26 ( 6 ): 825 – 34 . OpenUrl CrossRef PubMed Web of Science 43. ↵ Hatthakarnkul P , Quinn JA , Matly AAM , Ammar A , van Wyk HC , McMillan DC , et al. Systematic review of tumour budding and association with common mutations in patients with colorectal cancer . Crit Rev Oncol Hematol . 2021 Nov ; 167 : 103490 . 44. ↵ Siegel RL , Wagle NS , Cercek A , Smith RA , Jemal A . Colorectal cancer statistics, 2023 . CA Cancer J Clin . 2023 May 1 ; 73 ( 3 ): 233 – 54 . OpenUrl CrossRef PubMed 45. ↵ Nojadeh JN , Behrouz Sharif S , Sakhinia E . Microsatellite instability in colorectal cancer . EXCLI J . 2018 Jan 22 ; 17 : 159 – 68 . OpenUrl PubMed 46. ↵ Giannakis M , Hodis E , Jasmine Mu X , Yamauchi M , Rosenbluh J , Cibulskis K , et al. RNF43 is frequently mutated in colorectal and endometrial cancers . Nat Genet . 2014 Dec ; 46 ( 12 ): 1264 – 6 . OpenUrl CrossRef PubMed 47. ↵ Pai RK , Jayachandran P , Koong AC , Chang DT , Kwok S , Ma L , et al. BRAF-mutated, microsatellite-stable adenocarcinoma of the proximal colon: an aggressive adenocarcinoma with poor survival, mucinous differentiation, and adverse morphologic features . Am J Surg Pathol . 2012 May ; 36 ( 5 ): 744 – 52 . OpenUrl CrossRef PubMed 48. ↵ Trinh A , Lädrach C , Dawson HE , Ten Hoorn S , Kuppen PJK , Reimers MS , et al. Tumour budding is associated with the mesenchymal colon cancer subtype and RAS/RAF mutations: a study of 1320 colorectal cancers with Consensus Molecular Subgroup (CMS) data . Br J Cancer . 2018 Nov ; 119 ( 10 ): 1244 – 51 . OpenUrl CrossRef PubMed 49. ↵ Nguyen TT , Huynh TT , Le Nguyen P , Liew AWC , Yin H , Nguyen QVH . A Survey of Machine Unlearning [Internet] . arXiv [cs.LG] . 2022 . Available from: http://arxiv.org/abs/2209.02299 50. ↵ Kleppe A . Area under the curve may hide poor generalisation to external datasets . ESMO Open . 2022 Apr ; 7 ( 2 ): 100429 . OpenUrl PubMed 51. ↵ Bledsoe JR , Kamionek M , Mino-Kenudson M . BRAF V600E immunohistochemistry is reliable in primary and metastatic colorectal carcinoma regardless of treatment status and shows high intratumoral homogeneity . Am J Surg Pathol . 2014 Oct ; 38 ( 10 ): 1418 – 28 . OpenUrl CrossRef PubMed 52. ↵ Dawood M , Eastwood M , Jahanifar M , Young L , Ben-Hur A , Branson K , et al. Cross-linking breast tumor transcriptomic states and tissue histology . Cell Rep Med . 2023 Dec 19 ; 4 ( 12 ): 101313 . OpenUrl PubMed 53. ↵ Clark K , Vendt B , Smith K , Freymann J , Kirby J , Koppel P , et al. The Cancer Imaging Archive (TCIA): maintaining and operating a public information repository . J Digit Imaging . 2013 Dec 25 ; 26 ( 6 ): 1045 – 57 . OpenUrl CrossRef PubMed Web of Science 54. ↵ Zaidi SH , Harrison TA , Phipps AI , Steinfelder R , Trinh QM , Qu C , et al. Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival . Nat Commun . 2020 Jul 20 ; 11 ( 1 ): 3644 . OpenUrl PubMed View the discussion thread. Back to top Previous Next Posted February 08, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study Marco Gustav , Marko van Treeck , Nic G. Reitsam , Zunamys I. Carrero , Chiara M. Loeffler , Asier Rabasco Meneghetti , Bruno Märkl , Lisa A. Boardman , Amy J. French , Ellen L. Goode , Andrea Gsur , Stefanie Brezina , Marc J. Gunter , Neil Murphy , Pia Hönscheid , Christian Sperling , Sebastian Foersch , Robert Steinfelder , Tabitha Harrison , Ulrike Peters , Amanda Phipps , Jakob Nikolas Kather medRxiv 2025.02.04.25321660; doi: https://doi.org/10.1101/2025.02.04.25321660 Share This Article: Copy Citation Tools Assessing Genotype-Phenotype Correlations with Deep Learning in Colorectal Cancer: A Multi-Centric Study Marco Gustav , Marko van Treeck , Nic G. Reitsam , Zunamys I. Carrero , Chiara M. Loeffler , Asier Rabasco Meneghetti , Bruno Märkl , Lisa A. Boardman , Amy J. French , Ellen L. Goode , Andrea Gsur , Stefanie Brezina , Marc J. Gunter , Neil Murphy , Pia Hönscheid , Christian Sperling , Sebastian Foersch , Robert Steinfelder , Tabitha Harrison , Ulrike Peters , Amanda Phipps , Jakob Nikolas Kather medRxiv 2025.02.04.25321660; doi: https://doi.org/10.1101/2025.02.04.25321660 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Pathology Subject Areas All Articles Addiction Medicine (568) Allergy and Immunology (863) Anesthesia (299) Cardiovascular Medicine (4422) Dentistry and Oral Medicine (443) Dermatology (382) Emergency Medicine (607) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1507) Epidemiology (15219) Forensic Medicine (30) Gastroenterology (1123) Genetic and Genomic Medicine (6584) Geriatric Medicine (667) Health Economics (997) Health Informatics (4524) Health Policy (1367) Health Systems and Quality Improvement (1611) Hematology (540) HIV/AIDS (1264) Infectious Diseases (except HIV/AIDS) (15910) Intensive Care and Critical Care Medicine (1103) Medical Education (623) Medical Ethics (145) Nephrology (667) Neurology (6587) Nursing (345) Nutrition (998) Obstetrics and Gynecology (1143) Occupational and Environmental Health (956) Oncology (3331) Ophthalmology (970) Orthopedics (369) Otolaryngology (420) Pain Medicine (435) Palliative Medicine (129) Pathology (663) Pediatrics (1690) Pharmacology and Therapeutics (691) Primary Care Research (710) Psychiatry and Clinical Psychology (5437) Public and Global Health (9218) Radiology and Imaging (2195) Rehabilitation Medicine and Physical Therapy (1369) Respiratory Medicine (1195) Rheumatology (593) Sexual and Reproductive Health (709) Sports Medicine (529) Surgery (709) Toxicology (99) Transplantation (289) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'9ff7b3617ef141e2',t:'MTc3OTQwOTcwNQ=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

Ask this paper AI returns verbatim quotes from the full text · source: preprint-html

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2025) — citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00
unpaywall
last seen: 2026-05-21T05:10:58.409756+00:00
License: CC-BY-4.0