Abstract
Background:
In colorectal cancer the most frequently analyzed markers are RAS/RAF mutations which have negative predictive role for anti-EGFR antibody therapy. Although the hotspot mutations of the K- and NRAS exon 2-4 are well defined, the predictive role of rare mutations of other codons are questionable and the current guidelines are less precise in this respect.
Methods:
We have retrospectively re-analyzed the full spectrum of RAS/RAF mutations by next-generation sequencing (NGS) in a small cohort of nineteen cetuximab treated metastatic colorectal cancer patients whose enrollment was based on KRAS exon 2 mutation analysis.
Results:
We have detected 13 rare RAS/RAF mutations where 5 occurred together with hotspot (HSP) mutations. In the five non-HSP mutant cases we have detected 8 rare/non-HSP RAS/RAF mutations: KRAS-E37.,C51Y, NRAS-D57N,E62K,G60E/V14I, HRAS-P34L and BRAF-G464R. Kaplan-Meier analysis indicated that the progression free survival (PFS) of non-HSP cases ran parallel with wild type (WT) cases separately from the HSP mutant cases without statistical significance (p = 0.156). However, the combination of WT group with non-HSP cases resulted in a near significant survival difference from HSP cases (p = 0.057, HR = 3.501). Kaplan-Meier analysis of the overall survival (OS) again demonstrated that the survival curve of WT cases ran parallel with non-HSP cases, clearly separated from the HSP cases (p = 0.024). Overall survival of WT cases was significantly different as compared to HSP cases (p = 0.005) with HR of 6.362 (CI 1.496–27.055) as expected, while it was not statistically different from non-HSP cases (p = 0.092) with a HR of 0.891 (CI 0.248–3.205).
Conclusion:
In this small retrospective study we showed that the analyzed non-HSP mutations of RAS/RAF in colorectal cancer treated with cetuximab do not have negative effect on the efficacy. Since our cohort is very small and do not have accordingly the necessary statistical power of the survival analyses, therefore can only be considered as hypothesis generating one of the issue of the predictive role of non-hotspot RAS mutations on cetuximab efficacy.
Introduction
RAS oncogenes are the most significant drivers of colorectal cancer with a ∼40% incidence of KRAS, <10% incidence of NRAS and a rare occurrence of HRAS [, ]. Exon 2 mutations at codon 12 and 13 are the most prevalent followed by codon 59/61 of exon 3 and codons 117/146 in codon 4. Concerning allelic variants, G12D, G12V and G13D are the most frequent „hotspot” mutations of RAS that drastically affect its functionality. Notably, switching of RAS proteins between ON and OFF states requires facilitation by regulatory proteins (guanin nucleotide exchange factor,GEF and GTPase activating protein, GAPs), like SOS1 and NF1, respectively [, ]. Among these, GAPs are huge, multidomain proteins that trigger hydrolysis of GTP to GDP upon binding to RAS proteins. Of note, RAS proteins possess intrinsic hydrolytic activity, though it is a very slow process []. Thus, switching OFF RAS signaling depends on the activity of GAP proteins. As a proof of principle, hotspot mutations affect primarily codons that are implicated in GAP binding (G12) and GTP hydrolysis (Q61) and render these proteins insensitive to GAP-mediated inactivation [, ]. Mutation at these sites also drastically diminishes the intrinsic GTPase activity of RAS proteins.5 Interestingly, the KRAS G12C mutation was found to be an exception, as this mutational variant retains ∼75% of the intrinsic GTPase activity of wild-type proteins []. Furthermore, G13D, A59T and K117 show significantly higher nucleotide exchange rates compared to other hotspot mutations, while diminishing GAP sensitivity and/or intrinsic hydrolytic activity [, ]. Mutations in A146 are also implicated in enhanced nucleotide exchange, without affecting GAP sensitivity [].
Another member of the RAS-RAF-MAPK signaling pathway, BRAF, is also relatively frequently mutated in colorectal cancer (predominantly in the right-sided variants) with an incidence of ∼10%, predominantly at the codon V600 of the exon 15, which results in a constitutive elevation of its kinase activity [, ].
Anti-EGFR antibody therapies were introduced to colorectal cancer management almost two decades ago, with a controversial predictive marker of EGFR protein expression [, ]. Initial observational studies have shown that metastatic colorectal cancer (mCRC) patients with KRAS mutations in exon 2 (codons 12 and 13) are less likely to respond to anti-EGFR therapies. Further research revealed several additional mutations in exons 2, 3, and 4 of KRAS and NRAS that influence responsiveness to anti-EGFR therapies. Consequently, in July 2013, the European Medicines Agency (EMA) limited the use of panitumumab to patients with metastatic colorectal cancer (mCRC) whose tumors are RAS wild type (i.e., no mutations in exons 2, 3, or 4 of KRAS and NRAS, including codons 12, 13, 59, 61, 117, and 146). Subsequently, in November 2013, the EMA ruled that cetuximab should only be administered to metastatic CRC patients who do not carry RAS mutations. Mutational status should be determined by an experienced laboratory using validated test methods for the detection of KRAS and NRAS (exons 2, 3, and 4) mutations [–17]. It was another important step in this field when clinical studies (mostly retrospective ones) discovered that the BRAF V600 mutation is also a negative predictive marker of the anti-EGFR antibody therapies [, ].
With the widespread use of modern NGS technologies, molecular diagnostics become more complex and reveal interesting novel mutations in tumors including colorectal cancer. A retrospective analysis of 4122 anti-EGFR antibody-treated colorectal cancer patients was performed, where various more sensitive NGS-based tests revealed 54 “rare” mutant cases, which demonstrated a lower response rate (RR) and disease control rate (DCR) [18]. However, this cohort contained nine KRAS G12-, five A59 and one K117 hot spot mutants and no PFS and OS analysis was performed.
We examined a retrospective cohort of KRAS exon 2 wild type Cetuximab-treated metastatic colorectal cancer patients, where previous analysis found, on the contrary, what would be expected, that low EGFR expression of tumors was associated with better efficacy of anti-EGFR antibody therapy [19]. A portion of these cases was suitable for NGS sequencing to reveal the complete RAS/RAF mutation profiles and to test the cetuximab efficacy.
Materials and methods
Patients
In a previous study, we collected clinical data of colorectal cancer patients (all adenocarcinomas) treated with cetuximab at the Hungarian Defense Forces Medical Center [19]. Initial molecular diagnostics included KRAS exon 2 mutation testing according to the then-existing guidelines [20]. In this observational clinical study a retrospective next-generation sequencing (NGS) analysis was performed, where we have collected good-quality formalin-fixed paraffin embedded (FFPE) tissue samples from 19 anonymized metastatic patients. All samples were taken before any anti-EGFR therapies. The tumors were all KRAS exon 2 wild-type. The analysis was approved by the Institutional Review Board (IRB) of Semmelweis University (19/1043) and the full study protocol can be obtained from the IRB Office of Semmelweis University. Clinical data of the patients are shown in Table 1.
TABLE 1
| Case no. | Gender | Age | Location/primary | T/N stage/primary | Grade | Metastasis | Cetuximab treatment |
|---|---|---|---|---|---|---|---|
| 1 | Female | 40 | Sigmoidal | T3N1 | 2 | Ovarian | 2L |
| 2 | Female | 42 | Sigmoidal | T4N3 | 2 | Peritoneal | 2L |
| 3 | Male | 49 | Descending colon | T3N1 | 3 | Liver | 1L |
| 4 | Male | 54 | Rectal | TxNx | 2 | brain,ag,bone | 2L |
| 5 | Male | 56 | Colon transvers | T3N2 | 3 | Lung | 2L |
| 6 | Male | 57 | Colon transvers | T4N0 | 2 | Liver | 2L |
| 7 | Male | 60 | Coecal | T3N1 | 2 | Lung | 2L |
| 8 | Female | 61 | Rectal | T3N2 | 2 | Liver | 2L |
| 9 | Female | 63 | Coecal | T4N0 | 3 | Soft tissue | 2L |
| 10 | Male | 66 | Descending colon | T3N1 | 1 | Liver | 2L |
| 11 | Female | 67 | Flexura hepatica | T4N2 | 2 | Cutaneous | 2L |
| 12 | Male | 69 | Flexura hepatica | T3N2 | 3 | Liver | 2L |
| 13 | Male | 70 | Sigmoidal | T3N0 | 2 | Lung | 2L |
| 14 | Male | 72 | Rectosigmoidal | T3N2 | 2 | liver(ag) | 2L |
| 15 | Female | 73 | Rectal | T3N2 | 2 | Liver | 2L |
| 16 | Male | 74 | Sigmoidal | T3N0 | 2 | Liver | 2L |
| 17 | Male | 75 | Rectosigmoidal | T3N1 | 2 | Liver | 1L |
| 18 | Male | 78 | Coecal | T4N0 | 2 | Liver | 2L |
| 19 | Male | 78 | Sigmoidal | T3N1 | 2 | Liver | 1L |
Cetuximab treated colorectal cancer patient cohort.
1L = first line therapy, 2L = second line therapy, ag = adrenal gland.
NGS sequencing
NGS testing was performed in a blind manner concerning the clinical information of the involved patients. Samples were assessed by light microscopy and predominantly tumor cell-containing areas were macrodissected. The HighPur PCR Template Kit (cat:11796828001, Hoffman LaRoche, Basel, Switzerland) was used to isolate genomic DNA. Quantity and purity of DNA were assessed in a Qubit 4 Fluorometer using the Qubit dsDNA HS Assay Kit (cat,Q32850, Thermo Fisher Scientific, Kalden, Germany). Library preparation was performed with the Ion AmpliSeq Cancer Hotspot Panel.v2 (cat:4475348,Thermo Fisher Scientific). Ten ng of genomic DNA was subjected to multiplex PCR amplification using gene-specific primer pools, followed by enzymatic digestion and ligation of Ion Xpress Barcode Adapters. Barcoded libraries were quantified with the Ion Library TaqMan Quantitation Kit (cat:4468802, Thermo Fisher Scientific) and processed on the Ion Chef System (cat:4484177, Thermo Fisher Scientific) for automated template preparation and chip loading. Sequencing was carried out on the Ion PGM Sequencing System with Ion 318 chips (cat:4484354, Thermo Fisher Scientific) following the manufacturer’s protocol with 5% variant allel frequency (VAF) sensitivity. Coverage uniformity: the percentage of bases covered at ≥20% of the mean coverage were >95%. On-target reads (percentage of reads that mapped to target regions, out of total mapped reads per run) were >90%. Average depth of coverage was >300x. SNP detection sensitivity: 98% detection rate for 5% variant frequency at positions with average amplicon length was 115 bp. Primary data processing, including base calling, alignment, and variant calling, was performed using Ion Torrent Suite Software. Variant annotation and filtering were conducted in Ion Reporter Software (v5.20) using the NCBI clin Var, and the Franklin (Qiagen, Hilden, germany) database. Variant alleles have been categorized for pathogenicity [21].
Statistical analysis
Overall survival (OS) and progression-free survival (PFS) were defined as the time from start of cetuximab therapy to death and progression in days, respectively. OS and PFS were evaluated using the Kaplan–Meier method and log-rank statistics. Univariate risk factor analysis was performed by using Cox’s regression model. Differences were considered statistically significant when the p-value proved to be under 0.05. All statistical calculations were performed by Statistica 14.0.1.25 (TIBCO Software Inc, Palo Alto, CA, USA).
Results
This small cohort is composed of nineteen KRAS exon 2 wild type metastatic colorectal cancer patients (mCRC) where in the majority of cases cetuximab + FOLFIRI treatment was administered in second line after bevacizumab + FOLFIRI regimen (Table 1), while in three cases cetuximab + FOLFIRI was applied as first line therapy: wild type cases 17. and 19. and one hot spot (HSP) case 3. Therefore the non-HSP cases were all second line treatments. In most cases, primary tumors were available for NGS analysis, but in several cases, only metastatic samples could be tested due to various quality issues of the primary tumor samples. We have confidently used metastatic samples since previous large-scale genomic studies confirmed that the major driver gene profiles are maintained in metastases of colorectal cancer [22, 23].
We have used a hotspot panel containing K-/N-/H-RAS and BRAF for NGS sequencing; raw data are available in Supplementary Table 1. Mutational profiles of the cases are presented in Table 2. The classification of the mutants was done according to various databases and classified as pathogenic or likely pathogenic. Accordingly, this cohort contained six confirmed RAS/RAF wild type (WT) cases, The rest of the cohort contained four KRAS-, two NRAS-, and three BRAF mutant cases. Furthermore, there was one KRAS/NRAS-, one KRAS/NRAS/BRAF-, one KRAS/HRAS- and one NRAS/HRAS co-mutant case. For the classification of mutant types, we have considered hotspot (HSP) mutations in RAS at the codons 12/13/59/61/117/146 while in BRAF the codon V600 [, 24]. In case of co-occurring mutations, when any of them was in the hotspot codon, the case was considered HSP. We have detected 13 rare RAS/RAF mutations, where 5 occurred together with HSP mutations. The other 8 RAS/RAF rare mutations occurred in five cases. With such stratification the cohort contained eight HSP cases, five non-HSP cases, and six WT cases (Table 2). This cohort stratified by HSP versus non-HSP mutations and the wild type cases is clinically relatively heterogenous. At first, two wild-type cases (17,19) and one HSP case (3) were treated by cetuximab in first line but all the non-HSP cases obtained second line treatment. Secondly, the staging of the cases were also relatively heterogenous (Table 3.). Furthermore, the metastatic site pattern was also relatively heterogenous, since WT cases mostly had liver metastases similar to non-HSP cases but the HSP cases were heterogenous in this respect (Table 3.).
TABLE 2
| Case no. | Tumor type | KRAS | NRAS | HRAS | BRAF | Mutation category |
|---|---|---|---|---|---|---|
| 1 | Primary | - | D57N/E62K | P34L | - | Non-HSP |
| 2 | Primary | Q61H | - | G77N | - | HSP |
| 3 | Primary | A146T | - | - | | HSP |
| 4 | Adrenal gland | - | - | - | - | wt |
| 5 | Primary | - | G12N | - | - | HSP |
| 6 | Liver | - | - | - | - | wt |
| 7 | Primary | Q61R/E63K | G12N/S145L | - | - | HSP |
| 8 | Primary | E37 | - | - | - | Non-HSP |
| 9 | Primary | - | - | - | V600E | HSP |
| 10 | Liver | T35A | G12N | - | L584E | HSP |
| 11 | Primary | - | - | - | V600E | HSP |
| 12 | Primary | Q61L | - | - | - | HSP |
| 13 | Primary | - | G60E/V14I | | - | Non-HSP |
| 14 | Liver | - | - | - | - | wt |
| 15 | Primary | - | - | - | G464R | Non-HSP |
| 16 | Primary | - | - | - | | wt |
| 17 | Primary | - | - | - | - | wt |
| 18 | Liver | C51Y | - | - | - | Non-HSP |
| 19 | Primary | - | - | - | - | wt |
RAS/RAF mutation status of the cohort.
HSP, hotspot mutation (bold), non-HSP, non-hotspot mutation, wt, wild type.
TABLE 3
| | WT | | | HSP | | | Non-HSP | |
|---|---|---|---|---|---|---|---|---|
| Case | T/N stage | met | Case | T/N stage | met | Case | Stage | met |
| 4 | TxNx | B, Br | 2 | T4N3 | P | 1 | T3N1 | O |
| 6 | T4N0 | Li | 3 | T3N1 | Li | 8 | T3N2 | Li |
| 14 | T3N2 | Li | 5 | T3N2 | L | 13 | T3N0 | L |
| 16 | T3N0 | Li | 7 | T3N1 | L | 15 | T3N2 | Li |
| 17 | T3N1 | Li | 9 | T4N0 | S | 18 | T4N0 | Li |
| 19 | T3N1 | Li | 10 | T3N1 | Li | | | |
| | | | 11 | T4N2 | C | | | |
| | | | 12 | T3N2 | Li | | | |
Clinical parameters of the cohort by the mutational statuses of RAS/RAF.
B = bone, Br = brain, C = cutaneous, HSP = hotspot, L = lung, Li = liver, O = ovarian, P = peritoneal, S = soft tissue, WT = wild-type.
Kaplan-Meier analysis of the PFS indicated that the WT and non-HSP cases show similar survival curves, while HSP cases separated, however, without statistical significance (p = 0.156) (Figure 1). The direct comparison of the WT cases to HSP cases indicated a statistically significant difference (p = 0.033) with a hazard ratio (HR) of 4.161, confidence intervalls (CI) 1.496–27.055 as expected but the comparison of wild type cases to non-HSP cases revealed no statistically significant difference (p = 0.544) with a HR of 1.453 (CI 0.414–5.096). When we have combined the WT cases with the non-HSP cases and compared them to HSP cases, the Kaplan-Meier analysis indicated a near-significant difference in PFS (p = 0.057) with an HR of 3.501 (CI 1.12–10.942) (Figure 2). The combination of HSP cases with those with non-HSP mutations, the comparison did not result in a significant survival difference from the WT cases (p = 0.137) (Supplementary Figure 1).
FIGURE 1
FIGURE 2
Kaplan-Meier analysis of the OS again demonstrated that the survival curve of WT overlaps on non-HSP cases but clearly separated from the HSP cases (p = 0.024) (Figure 3). Overall survival of WT cases was significantly different as compared to HSP cases (p = 0.005) with HR of 6.362 (CI 1.496–27.055) as expected, while it was not statistically different from non-HSP cases (p = 0.092) with a HR of 0.891 (CI 0.248–3.205) (Figure 3). This analysis was repeated with a combination of WT cases with the non-HSP cases into one group, which indicated a significant difference from HSP cases (p = 0.005) with a HR of 6.693 (CI 1.749–25.61) (Figure 4). Again, the combination of HSP cases with non-HSP cases did not result in a significant survival difference from the WT cases (p = 0.262) (Supplementary Figure 2).
FIGURE 3
FIGURE 4
Discussion
When BRAF inhibitors entered the clinic for the treatment of BRAF mutant malignancies it became clear that not all mutants are equal, as some of the mutations confer sensitivity to these treatments, but others do not [25]. Later on, the classification of BRAF mutations was established based on genomic and functional analyses into class I-III variants based on the kinase activities and sensitivity to BRAF inhibitors. Class-I BRAF mutations are codon V600 of the exon 15 with very high kinase activity, class-II mutations cause moderate elevation of kinase activity, while class-III mutations result in decreased kinase activity [26]. This issue became critical in colorectal cancer when it was discovered that class-I BRAF mutations confer resistance to anti-EGFR antibody therapies [27, 28], later on leading to the development of a combination strategy using BRAF inhibitors [29]. In our small retrospective cohort, there were two Class-I and one Class-II BRAF mutant case, as well as one case of unclassified non-HSP mutation, co-occurring with HSP mutation of NRAS.
It is interesting that in case of RAS mutations, the most frequent mutation type in colorectal cancer with strong clinical relevance in relation to anti-EGFR antibody therapies, such a classification has not been developed yet. Although it is established that the K- and N-RAS hotspot mutations are at the codons 12/13/59/61/117/146 of the exons 2–4, other mutations are not tested rigorously for their significance in relation to anti-EGFR antibody therapies of colorectal cancer. Such rare mutations became more and more frequent using sensitive molecular techniques involving NGS [18, 24, 30]. In our cetuximab-treated cohort, we have detected 11 K-, N-, and HRAS rare mutations, out of which 7 mutations were found in four so-called non-HSP cases, depicted on Figure 5, while one non-HSP case was a rare BRAF mutant one. Notably, although no oncogenic role was attributed to the mutations found in our non-HSP group, their localization shows a clear pattern as they are present in key structural positions, such as P-loop (GTP/GDP binding), switch I and II (protein-protein interactions), as shown in Figure 5 [].
FIGURE 5
The rare mutations of NRAS, D57N and G60E have been observed before [18, 24]. Functional assay indicated that the D57N variant does not change NRAS activity while the G60E results in a moderate increase [24]. It is of note that this variant was previously detected in Noonan syndrome, a special form of RASopathy [31, 32]. The possible functional activity of the five other rare variants: KRAS-E37/C51Y, NRAS-V14I/E62K, and HRAS-P34L are not known. There is a recent biochemical and bioinformatic study which is based on 935 HSP and non-HSP RAS variants, but further functional analysis was conducted only on a small portion of these mutations, which did not contain the rare variants found in our cohort [33]. Treatment of KRAS-mutant cancer has changed recently with the introduction of G12C allosteric irreversible inhibitors [34]. Although this was a real breakthrough, in the case of both clinically used inhibitors, sotorasib and adagrasib, developing resistance poses a significant challenge, involving several novel KRAS mutations not only in the hotspot regions but in other locations [35, 36]. However, none of those „rare” mutations have been observed in our presented colorectal cancer cases.
In a previous retrospective analysis of 54 „rare” KRAS mutant colorectal cancer cases, anti-EGFR antibody efficacy (RR and DCR) was significantly affected. Meanwhile, this cohort with “rare” mutations also involved 15 well-established hotspot mutations, which have possibly affected the analysis [18]. The study only described one variant, an NRAS G60 mutation, that occurred in our cohort. Another preclinical study analyzing rare KRAS mutations for the possible effect on cetuximab treatment demonstrated that these mutations have a wide ranges of effects on the cetuximab efficacy from resistance through neutrality or preserved sensitivity [30]. Unfortunately, none of those mutations that were included in the cetuximab testing have been present in our cohort.
Molecular diagnostics must evaluate the found genetic alteration(s) according to the therapeutic level of evidence, either following the guidance of OncoKB [37] or ESCAT [38]. Since the labels of both anti-EGFR antibodies suggest excluding RAS exon 2–4 mutant patients, the issue of rare RAS mutations is clinically relevant. Our analysis of the clinical effects of the analyzed „rare” mutations of KRAS on the cetuximab efficacy of colorectal cancer patients demonstrated that the PFS and OS of these patients were not different from the WT cases. The limitation of our retrospective analysis is evident: the cohort is small, especially the „rare” mutational group, which also included a class-II BRAF mutant case. Meanwhile, this analysis suggests that similar retrospective analysis must be performed on large cancer series since the previous notion, that all exon 2–4 mutations of K- or N-RAS cases are resistant to cetuximab therapy is challenged by our data and is also supported by recent experimental analysis of the effect of rare mutations on the KRAS function and/or cetuximab sensitivity [22, 28, 31]. In summary, without careful evaluation of the functional consequences of the individual mutations, the clinical outcome cannot be properly predicted.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by IRB of Semmelweids University 19/1043. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because the study is retrospective, patient personal data are not used and the findings do not affect treatment.
Author contributions
JT was responsible for the design execution of the project as well as for the drafting of the manuscript. NGS and raw data analysis were performed by ER and TB. RAS mutation analysis was done by MB. Statistical analysis was performed by IK. Patient cohort was selected by AU and clinical data also have been collected and evaluated by ZP and AU. All authors contributed to the article and approved the submitted version.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National RD Innovation Office, Hungary (K-135540, A-152516, 2020.1.1.6.-Jovo-2021-0004, TKP/EGA-2025) to provide consumables and to cover the salary of molecular biologists.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.por-journal.com/articles/10.3389/pore.2026.1612560/full#supplementary-material
References
1.
TímárJKashoferK. Molecular epidemiology and diagnostics of KRAS mutations in human cancer. Cancer Metastasis Rev (2020) 39:1029–38. 10.1007/s10555-020-09915-5
2.
ScharpRBBalanARicciutiBFikselJCherryCWangCet alGenomic landscapes and hallmarks of mutant RAS in human cancer. Cancer Res (2022) 82:4058–78. 10.1158/0008-5472.CAN-22-1731
3.
BosJLRehmannHWittinghoferA. GEFs and GAPs: critical elements in the control of small G proteins. Cell (2007) 129(5):865–77. 10.1016/j.cell.2007.05.018
4.
SimanshuDKNissleyDVMcCormickF. RAS proteins and their regulators in human disease. Cell (2017) 170(1):17–33. 10.1016/j.cell.2017.06.009
5.
HunterJCManandharACarrascoMAGurbaniDGondiSWestoverKD. Biochemical and structural analysis of common cancer-associated KRAS mutations. Mol Cancer Res (2015) 13(9):1325–35. 10.1158/1541-7786.MCR-15-0203
6.
PriorIALewisPDMattosC. A comprehensive survey of ras mutations in cancer. Cancer Res (2012) 72(10):2457–67. 10.1158/0008-5472.CAN-11-2612
7.
HaigisKM. KRAS alleles: the devil is in the detail. Trends Cancer (2017) 3:686–97. 10.1016/j.trecan.2017.08.006
8.
OrlandiACalegariAInnoABerenatoRCaporaleMNigerMet alBRAF in metastatic colorectal cancer: the future starts now. Pharmacogenomics (2015) 16:2069–81. 10.2217/pgs.15.140
9.
ChenYSunHDengYMaYHuangHLiuYet alThe clinical and genomic distinctions of class1/2/3 BRAF-Mutant colorectal cancer and differential prognoses. Biomarker Res (2023) 11:11. 10.1186/s40364-022-00443-8
10.
HechtJRmitchellEneubauerMABurrisHASwansenPLopezTet alLack of correlation between epidermal growth factor receptor status and response to panitumumab monotherapy in metastatic colorectal cancer. Clin Cancer Res (2010) 16:2205–13. 10.1158/1078-0432.CCR-09-2017
11.
LicitraLStörkelSKerrKMVan CutsemEPirkerRHirschFRet alPredictive value of epidermal groth factor receptor expression for first line chemotherapy plus cetuximab in patients with head and neck and colorectal cancer. Analysis of data from the EXTREME and crystal studies. Eur J Cancer (2013) 49:1161–8. 10.1016/j.ejca.2012.11.018
12.
DufraingKKeppensCTackVSiebersAGKafatosGDubeSet alEvolution of RAS testing over time: factors influencing mutation rates in metastatic colorectal cancer patients. Colorectal Cancer (2020) 9:CRC14. 10.2217/crc-2019-0013
13.
DeanLKaneM. Cetuximab therapy and RAS and BRAF genotype. In: Medical Genetics Summaries-NCBI Bookshelf (2020). p. p1–15.
14.
BensonABVenookAPAdamMChangGChenY-JCiomborKKet alColon Cancer, Version 3.2024. J Natl Compr Cancer Netw (2024). p. 22. 10.6004/jnccn.2024.0029
15.
Erbitux. EPAR – product information-tracked changes- EMA (2025). Available online at: www.ema.europa.eu\documents\erbitux (Accessed August 07, 2026).
16.
Vectibix,INN-panitumumab-EMA (2025). Available online at: www.ema.europa.eu\health\anx_110919_hu (Accessed August 07, 2026).
17.
CremoliniCChalabiMFassanMGelliMGoereDElezEet alMetastatic colorectal cancer: clinical practice guideline for diagnosis, treatment and follow up. Ann Oncol (2026) 37:759–76. 10.1016/j.annonc.2026.03.005
18.
SuzukiSSaitoYSaitoKYamadaYTakahashiKKumanishiRet alLimited efficacy of anti-EGFR monoclonal antibodies in colorectal cancer patients with rare RAS variants: analysis of the C-CAT database. Curr Issues Mol Biol (2024) 46:144476–14486. 10.3390/cimb46120869
19.
UhlyarikAPiurkoVPápaiZRasoELahnEKissEet alEGFR protein expression in KRAS wild-type metastatic colorectal cancer is another negative predictive factor of the cetuximab therapy. Cancers (2020) 12:614. 10.3390/cancers12030614
20.
CserepesMOstorosGLohinaiZRasoEBarbaiTTimarJet alSubtype specific KRAS mutations in advanced lung adenocarcinoma. a retrospective study of patients treated with platinum-based chemotherapy. Eur J Cancer (2014) 50:1819–28. 10.1016/j.ejca.2014.04.001
21.
StrandeNTBrnichSERomanTSBergJS. Navigating the nuances of clinical sequence variant interpretation in Mendelian disease. Genet Med (2018) 20:918–26. 10.1038/s41436-018-0100-y
22.
IshaqueNAbbaMLHauserCPatilNParamasivamNHuebschmannDet alWhole genome sequencing puts forward hypothesis on metastasis evolution and therapy in colorectal cancer. Nat Com (2018) 19:4782. 10.1038/s41467-018-07041-z
23.
HuZDingJMaZSunRSeoaneJAScott ShafferJet alQuantitative evidence for early metastatic seeding in colorectal cancer. Nat Genet (2019) 51:1113–22. 10.1038/s41588-019-0423-x
24.
LoreeJMWangYSyedMASorokinAVCokerOXiuJet alClinical and functional characterization of atypical KRAS/NRAS mutations in metastatic colorectal cancer. Clin Cancer Res (2021) 27:4587–98. 10.1158/1078-0432.CCR-21-0180
25.
GoudaMASubbianV. Precision oncology for BRAF-mutant cancers with BRAF and MEK inhibitors: from melanoma to tissue-agnostic therapy. ESMO Open (2023) 8:1–14. 10.1016/j.esmoop.2023.100788
26.
OwsleyJSteinMKPorterJInGKSalemMO'DaySet alPrevalence of classI-III BRAF mutations among 114,662 cancer patients in a large genomic database. Exp Biol Med (2021) 246:31–9. 10.1177/1535370220959657
27.
SchirripaMBiasonPLonardiSPellaNPinoMSUrbanoFet alClass1,2 and 3 BRAF mutated metastatic colorectal cancer: a detailed clinical, pathologic, and molecular characterization. Clin Cancer Research (2019)25:3954–60. 10.1158/1078-0432.CCR-19-0311
28.
SahinIH. Klostergaard. BRAF mutations as actionable targets: a paradigm shift in the management of colorectal cancer and novel avenues. JCO Pres Oncol Parct (2021) 17:723–31. 10.1200/OP.21.00160
29.
ElezEYoshinoTShenLLonardiSVan CutsemEEngCet alEncorafenib, cetuximab and mFOLFOX6 in BRAF-mutated colorectal cancer. N Engl J Med (2025) 392:2425–37. 10.1056/NEJMoa2501912
30.
ParkSIkegamiMKitadaRTakamochiKSuzukiKManoHet alFunctional analysis of rare RAS variants of unknown significance. Cancer Res Commun (2025) 5:1747–57. 10.1158/2767-9764.CRC-25-0188
31.
RauenKA. Defining RASopathy. Dis Models Mech (2022) 15:1–5. 10.1242/dmm.049344
32.
RazzaqueAKomolkeYNishizawaTInaiKFurutaniMHigashinagawaTet alCharacterization of a novel KRAS mutation identified in Noonan syndrome. Am J Med Genet (2021) 158A:524–32. 10.1002/ajmg.a.34419
33.
TripathiSDsouzaNRMathisonAJLeverenceEUrrutiaRZimmermannMT. Enhanced interpretation if 935 hotspot and non-hotspot RAS variants using evidence-based structural bioinformatics. Comp Struct Biotech J (2022) 20:117–37. 10.1016/j.csbj.2021.12.007
34.
RiedlJMMatsubaraHMcNeilRPatelPde la CruzFFGulhanDCet alEmerging landscape of KRAS inhibitors in cancer treatment. Cancer Cell (2026) 44:471–97. 10.1016/j.ccell.2026.01.001
35.
AwadMMLiuSRybkinIIArbourKCDillyJZhuVWet alAcquired resistance to KRASG12C inhibition in cancer. N Engl J Med (2021) 384:2382–93. 10.1056/NEJMoa2105281
36.
RiedlJMFece de la CruzFLinJJParseghianCKimJEMatsubaraHet alGenomic landscape of clinically acquired resistance alterations in patients treated with KRASG12C inhibitors. Ann Oncol (2025) 6:682–92. 10.1016/j.annonc.2025.01.020
37.
ChakravartyDGaoJPhillipsSMKundraRZhangHWangJet alOncoKB: a precision oncology knowledge base. JCO Precis Oncol (2017) 1:1–16. 10.1200/PO.17.00011
38.
MateoJChakravartyDDienstmannRJezdicSGonzalet-PerezALopez-BigasNet alA framework to rank genomic alterations as targets for cancer precision medicine: the ESMO Scale for clinical actionability of molecular targets (ESCAT). Ann Oncol (2018) 29:1895–902. 10.1093/annonc/mdy263
Summary
Keywords
BRAF mutation, cetuximab therapy, colorectal cancer, rare mutation, RAS mutation
Citation
Tímár J, Rásó E, Barbai T, Baranyi M, Kenessey I, Pápai Z and Uhlyarik A (2026) A retrospective analysis of cetuximab efficacy based on the comprehensive RAS/RAF mutation profiles of colorectal cancer. Pathol. Oncol. Res. 32:1612560. doi: 10.3389/pore.2026.1612560
Received
23 July 2026
Revised
09 October 2026
Accepted
23 October 2026
Published
03 November 2026
Volume
32 - 2026
Edited by
Andrea Ladányi, National Institute of Oncology, Hungary
Updates
Copyright
© 2026 Tímár, Rásó, Barbai, Baranyi, Kenessey, Pápai and Uhlyarik.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: József Tímár, jtimar@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.