Next Article in Journal
Association of a Palliative Surgical Approach to Stage IV Pancreatic Neuroendocrine Neoplasms with Survival: A Systematic Review and Meta-Analysis
Next Article in Special Issue
The FMS like Tyrosine Kinase 3 (FLT3) Is Overexpressed in a Subgroup of Multiple Myeloma Patients with Inferior Prognosis
Previous Article in Journal
Genomic Mapping Identifies Mutations in RYR2 and AHNAK as Associated with Favorable Outcome in Basal-Like Breast Tumors Expressing PD1/PD-L1
Previous Article in Special Issue
Utility of Circulating Tumor DNA for Detection and Monitoring of Endometrial Cancer Recurrence and Progression
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 12
Issue 8
10.3390/cancers12082245
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Somatic Genetic Variants as Predictors of Resistance to EGFR-Targeted Therapies in Metastatic Colorectal Cancer Patients

by
Pau Riera
1,2,3,4,5,*,
Benjamín Rodríguez-Santiago
1,4,5,
Adriana Lasa
1,4,5,
Lidia Gonzalez-Quereda
1,4,5,
Berta Martín
6,
Juliana Salazar
7,
Ana Sebio
6,
Anna C. Virgili
6,
Jordi Minguillón
1,4,5,8,
Cristina Camps
1,3,4,
Jordi Surrallés
1,4,5,8 and
David Páez
5,6,*
1
Genetics Department, Hospital de la Santa Creu i Sant Pau, 08041 Barcelona, Spain
2
Pharmacy Department, Hospital de la Santa Creu i Sant Pau, 08041 Barcelona, Spain
3
Faculty of Pharmacy and Food Sciences, Universitat de Barcelona (UB), 08028 Barcelona, Spain
4
Join Research Unit on Genomic Medicine UAB-IR Sant Pau, Biomedical Research Institute, Hospital de la Santa Creu i Sant Pau, 08041 Barcelona, Spain
5
U705 and U745, ISCIII Center for Biomedical Research on Rare Diseases (CIBERER), 08041 Barcelona, Spain
6
Medical Oncology Department, Hospital de la Santa Creu i Sant Pau, 08041 Barcelona, Spain
7
Translational Medical Oncology Laboratory, Institut de Recerca Biomèdica Sant Pau (IIB-Sant Pau), 08041 Barcelona, Spain
8
Department of Genetics and Microbiology, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
*
Authors to whom correspondence should be addressed.
Cancers 2020, 12(8), 2245; https://doi.org/10.3390/cancers12082245
Submission received: 18 June 2020 / Revised: 17 July 2020 / Accepted: 6 August 2020 / Published: 11 August 2020
(This article belongs to the Collection Cancer Biomarkers)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: About 40% of RAS/BRAF wild-type metastatic colorectal cancer (mCRC) patients undergoing anti-EGFR-based therapy have poor outcomes. Treatment failure is not only associated with poorer prognosis but higher healthcare costs. Our aim was to identify novel somatic genetic variants in the primary tumor and assess their effect on anti-EGFR response. Patients and Methods: Tumor (somatic) and blood (germline) DNA samples were obtained from two well-defined cohorts of mCRC patients, those sensitive and those resistant to EGFR blockade. Genetic variant screening of 43 EGFR-related genes was performed using targeted next-generation sequencing (NGS). Relevant clinical data were collected through chart review to assess genetic results. Results: Among 61 patients, 38 were sensitive and 23 were resistant to treatment. We identified eight somatic variants that predicted non-response. Three were located in insulin-related genes (I668N and E1218K in IGF1R, T1156M in IRS2) and three in genes belonging to the LRIG family (T152T in LRIG1, S697L in LRIG2 and V812M in LRIG3). The remaining two variants were found in NRAS (G115Efs*46) and PDGFRA (T301T). We did not identify any somatic variants related to good response. Conclusions: This study provides evidence that novel somatic genetic variants along the EGFR-triggered pathway could modulate the response to anti-EGFR drugs in mCRC patients. It also highlights the influence of insulin-related genes and LRIG genes on anti-EGFR efficacy. Our findings could help characterize patients who are resistant to anti-EGFR blockade despite harboring RAS/BRAF wild-type tumors.

1. Introduction

The epidermal growth factor receptor (EGFR), usually overexpressed in colorectal cancer (CRC), plays a pivotal role in tumor growth and progression [1,2]. Multiple proteins are involved in the EGFR signaling pathway, including other receptors, ligands, intracellular downstream effectors, and regulators [2,3,4]. In the metastatic setting, anti-EGFR targeted antibodies, cetuximab and panitumumab, are commonly used, but response rates are variable [5,6]. Several somatic mutations along the EGFR-triggered pathway, such as activating RAS mutations and the BRAF V600E mutation, are validated predictors of primary resistance to anti-EGFR-based therapies [7,8,9,10,11,12]. Other promising biomarkers of non-response are PIK3CA or PTEN mutations, although the level of evidence is lower [13,14]. In addition, some studies have observed that right-sided and mesenchymal tumors show worse outcomes to EGFR-targeted therapies regardless of RAS mutation status [15,16,17,18]. As about 40% of RAS-wild-type patients undergoing anti-EGFR therapy do not benefit from this treatment [19,20], we hypothesized that other mutations in the EGFR pathway could act as additional mechanisms of resistance to EGFR blockade.
Next-generation sequencing (NGS) technologies have revolutionized research in cancer genomics. NGS allows for the simultaneous analysis from several genes to complete genomes with higher sensitivity and cost-effectiveness than the previously used Sanger sequencing methods. In addition, its high sensitivity has detected somatic variants in the tumor at a low allelic fraction.
In this study we used NGS technology to analyze the exons and intron boundaries of 43 EGFR pathway-related genes. We genotyped both germline and tumor DNA samples to optimize the identification of somatic genetic variants. The main objective was to identify novel genetic variants in two cohorts of extreme responder patients with RAS wild-type metastatic CRC (mCRC). Extreme responders were either primary resistant or highly sensitive to anti-EGFR therapy. We also aimed to increase knowledge of EGFR-related genes as this could lead to the identification of new biomarkers of anti-EGFR response and promising therapeutic targets for mCRC.

2. Patients and Methods

2.1. Patient Population

We conducted this case-control study with mCRC patients from the Hospital de la Santa Creu i Sant Pau (HSCSP, Barcelona, Spain). We retrospectively analyzed patients who underwent any anti-EGFR-containing regimen between 2012 and 2017. All tumor samples had been previously classified as RAS wild-type using the therascreen KRAS test (Qiagen, Hilden, Germany) or the TruSight Tumor 15 panel (Illumina, San Diego, CA, USA).
Two well-defined cohorts of patients were studied, those sensitive (control group) and those resistant (case group) to EGFR blockade. Response to the anti-EGFR-based treatment was determined by total body CT scan. The first CT scan reassessment was performed 2–3 months after treatment. Patients showing disease progression in the first CT scan were classified as resistant. Conversely, those achieving a complete or partial response at this time point, or disease stabilization lasting at least 6 months were considered sensitive. Response to anti-EGFR treatment was assessed according to RECIST (Response Evaluation Criteria in Solid Tumors) v1.1 [21]. An Eastern cooperative oncology group (ECOG) performance status ≤ 2 and age ≥ 18 was required for inclusion in the study. We excluded patients for whom tumor DNA was not available.
Clinical data collected from hospital records included gender, age, performance status (PS) according to the ECOG scale, smoking habit, primary tumor location, number of metastatic sites, time to metastases, resection of the primary tumor, previous lines of chemotherapy, type of anti-EGFR administered, and concomitant chemotherapy. Formalin-fixed paraffin-embedded (FFPE) primary tumor tissues were available from all patients. Germline DNA was also available for 92% of the patients. The study was approved by the Institutional Ethics Committee at HSCSP (ethical code: 22/2012) and all study participants gave written informed consent.

2.2. Gene Selection and Primer Design

Two custom panels were created, one for blood samples and the other for tumor samples. Custom amplicon oligonucleotides were designed for each region of interest following the manufacturer’s instructions (Illumina, San Diego, CA, USA). As FFPE tumor DNA is more degraded and fragmented than germline DNA, amplicons of less length are needed to achieve good quality sequencing reads. The panel for germline DNA therefore contained 771 amplicons with an average size of 250–300 base pairs (bp) whereas the panel for tumor DNA contained 1124 amplicons of ~175 bp. Both panels included the 43 candidate genes related to the EGFR pathway and had an expected coverage >99% for all the coding regions (~200,000 bp). These genes mainly encoded receptors, ligands, intracellular downstream effectors or proteins involved in EGFR turnover. We included the most relevant genes of the pathway and also those related to anti-EGFR response in previous studies [3,22,23]. Table 1 provides information about the selected genes and the function of their corresponding encoded proteins.

2.3. Isolation and Quantification of DNA

Germline DNA was automatically extracted from peripheral whole-blood samples (Autopure, Qiagen, Hilden, Germany). Tumor DNA was extracted from primary tumor tissue samples using the GeneRead DNA FFPE Kit (Qiagen, Hilden, Germany). This kit purifies tumor DNA by removing artificial C > T mutations. DNA concentrations were measured using Qubit™ dsDNA HS Assay Kits with the Qubit 3.0 Fluorometer (Invitrogen, Carlsbad, CA, USA).

2.4. Library Preparation, Sequencing Runs, and NGS Analysis

Library preparation was carried out following the manufacturer’s instructions (TruSeq Custom Amplicon Low Input Library Prep, Illumina, San Diego, CA, USA). We used ~15 ng of germline DNA or ~100 ng of tumor DNA, due to the low quality of DNA from FFPE samples. Library sizes were determined using QIAxcel DNA Screening Gel Cartridge on QIAxcel capillary electrophoresis system (Qiagen GmbH, Hilden, Germany). Library concentrations were measured with Qubit 3.0 Fluorometer (Invitrogen, Carlsbad, CA, USA). All libraries were subsequently diluted and pooled equimolarly.
Germline DNA sequencing was performed on a MiSeq platform (Illumina) to obtain 150-bp paired-end reads. Samples were sequenced using the MiSeq v2-300 Reagent Kit. To achieve higher coverage, tumor DNA sequencing was performed on a NextSeq 500 platform (Illumina), also obtaining 150-bp paired-end reads. Samples were sequenced using the NextSeq Mid v2-300 Reagent Kit.
Sequence reads were aligned against the human reference genome (version GRCh37) using the Burrows-Wheeler Aligner (BWA, version 0.7.12) [24]. Single nucleotide and indel variants were called by means of the Mutect2 tool (version 4.0.12.0) that can use both input tumor and germline data. It can also manage high coverage data from tumor sequences [25]. Alignment and calling were performed following software development best practices (Broad Institute, Cambridge, MA, USA). After alignment, a panel of normal variation (PoN) obtained from germline DNA sequencing data was used to exclude rare germline variants and individual-specific artifacts. The number of somatic variants per patient was calculated by filtering tumor DNA sequencing data with the PoN. We excluded variants located in intronic, intergenic or UTR regions, and also polymorphisms (GnomAD allele frequency >0.001). Results were inspected using the Integrative Genomics Viewer (IGV). Resistant patients whose tumors harbored a known non-response mutation (in KRAS, NRAS or BRAF V600E) of over 5% of mutated clones were excluded. The presence of these non-response mutations was confirmed by a second technique (Sanger sequencing or the TruSight Tumor 15 panel) in some patients. Finally, only those variants that were over 1% of mutant clones in all patients harboring them were kept for further examination. COSMIC cancer database v91 was used to assess whether the candidate somatic variants identified had been previously described [26]. Additional functional annotation of variants was performed using ANNOVAR [27]. Variant analyses and interpretation were performed using Alamut® Visual v2.15 software (SOPHiA GENETICS, Boston, MA, USA) and the Cancer Genome Interpreter platform (Institute for Research in Biomedicine, Barcelona, Spain), which is publicly available at http://www.cancergenomeinterpreter.org [28].

2.5. Statistical Analyses

We used the chi-square test to compare the baseline clinical characteristics between the two cohorts of patients and Fisher’s exact test to compare the prevalence of alterations between sensitive and resistant patients. The Mann-Whitney test was carried out to compare the number of somatic genetic variants between sensitive and resistant patients. All statistical tests were performed using R software (version 3.3.2., https://cran.r-project.org/bin/windows/base/old/3.3.2/) and IBM SPSS® statistics software (version 25, https://www.ibm.com/support/pages/downloading-ibm-spss-statistics-25). A 95% confidence level was set for all tests of significance.

3. Results

3.1. Patient Population

A total of 168 mCRC patients were treated with anti-EGFR-containing regimens between 2012 and 2017. Sixty-one of the patients (38 sensitive and 23 resistant to EGFR blockade) were included in the study, as they met the inclusion criteria and quality of tumor DNA was good (Figure 1). They were all diagnosed with stage IV CRC (37.7% metachronous, 50.8% ≥ 2 metastatic locations). Regarding treatment, 33 patients received cetuximab and 28 received a panitumumab-containing regimen, with no significant differences between the two cohorts (p = 0.814). There were more females and more patients with worse ECOG PS in the non-responder cohort (p < 0.001 and p = 0.013, respectively). No differences were observed between the two cohorts in respect to primary tumor location (p = 0.634). A higher percentage of patients in the sensitive group received the anti-EGFR-containing regimen as first-line treatment (47.4% in the sensitive group vs. 26.1% in the resistant group) and presented synchronous metastases (71.1% in the sensitive group vs. 47.8% in the resistant group). Baseline clinical features are described in Table 2.

3.2. Genetic Analyses

The mean target coverage for tumor samples was 3600, achieving 100× or greater coverage for 88% of the bases. For germline samples, the mean target coverage was 640, achieving 30× or greater coverage for 87% of the bases. As shown in Figure 2, we identified 46,512 somatic genetic variants. We found a mean of 707 variants per sample in sensitive patients and 854 variants per sample in resistant patients (p = 0.19).

3.3. KRAS, NRAS, BRAF, and PIK3CA Assessment and Patient Selection

Five anti-EGFR resistant patients were finally excluded from the analyses because a validated mutation of non-response (in RAS or BRAFV600E) was detected at an allelic fraction over 5%. One of the five patients (P11) had the KRAS G12C mutation. The remaining four patients (P3, P39, P55, and P64), all with right-sided tumors, presented the BRAF V600E mutation. In one patient (P64) this mutation had not been previously detected by Sanger. Patient P55 presented the mutations BRAF V600E and NRAS G13D at a frequency over 5% for both mutations. As for patient P39, BRAF V600E and KRAS Q61L mutations coexisted, although the allelic fraction for KRAS Q61L was only 0.8%. No patients harbored somatic mutations in both KRAS and NRAS genes. One sensitive patient (P28) who underwent FOLFIRI plus panitumumab as a second-line treatment presented the KRAS A146V mutation at an allelic fraction of 5.0%. This mutation had not been identified previously. Table 3 and Figure S1 show all the somatic mutations found in these genes and their allelic fractions.
We also analyzed PIK3CA mutations and BRAF mutations other than V600E. Two patients, one resistant (P63) and one sensitive (P66), presented the PIK3CA E545K mutation, with an allelic fraction of 10.1% and 6.1%, respectively. In addition, two BRAF mutations previously related to favorable prognosis (D594N and G466A) [29,30,31,32] were identified at around 20% in two sensitive patients (P45 and P59). Table 3 shows all the findings concerning BRAF and PIK3CA assessment.

3.4. Identification of Novel Genetic Variants Related to Anti-EGFR Response

We identified eight potential somatic variants of non-response at a frequency of over 1% in 12 out of 18 (66.7%) non-responder patients (Table 4). In 8 cases, two or more of these variants coexisted. Mutant allele fractions differed substantially among the patients (Figure S2). No potential resistance variants were found in 6 non-responders. Additionally, no variants of good response were identified.
Three of the eight variants detected were missense variants located in insulin-related genes, such as IGF1R I668N and E1218K or IRS2 T1156M. Three others were found in genes belonging to the LRIG family: LRIG1 (T152T), LRIG2 (S697L), and LRIG3 (V812M). The remaining two variants were found in NRAS (G115Efs*46) and PDGFRA (T301T). All variants were non-synonymous, except for LRIG1 and PDGFRA variants. According to Alamut software, these two variants may create a novel cryptic acceptor site identifiable by the splicing complex.

4. Discussion

We sought to investigate the existence of novel somatic variants in EGFR-related genes as predictive markers of response to anti-EGFR antibodies. We found eight potential somatic variants that could explain the lack of response to these agents, highlighting the variants in the insulin-related and LRIG family genes. In contrast, we did not find somatic variants related to good response.
Accurate identification of somatic variants is challenging. In the past, Sanger sequencing was the only technique available to detect somatic mutations in tumor samples. Currently, NGS technology is replacing Sanger method as it allows the sequencing of hundreds of genes simultaneously and shows high sensitivity [35]. These advantages have enabled the detection of concomitant mutations in several genes of interest, including low-allele-fraction variants not previously found by Sanger sequencing. Mutant allele fractions may influence the response to targeted therapies. In this sense, the mutant allele fraction that determines anti-EGFR response continues to be debated. A major point of discussion is whether the optimal threshold of the RAS mutant allele fraction to identify patients likely to benefit from anti-EGFR drugs should be 1% or 5%. We used a threshold of 5% as it has been reported that reducing the threshold to 1% does not improve outcomes [36,37]. We found one sensitive patient (P28) who harbored the A146V KRAS mutation at 5%. The good response in this patient could be related to the chemotherapy scheme concomitantly given with the anti-EGFR drug. We also found a novel somatic variant of non-response in the NRAS gene. NRAS variants routinely tested prior to anti-EGFR initiation are normally missense mutations (≈95%) [7]. Conversely, the new NRAS variant we found, G115Efs*46, is a truncating mutation. It consists of a deletion located in exon 4 that leads to a premature stop codon and a truncated protein [26]. The Cancer Genome Interpreter predicts that it is a passenger mutation with a highly deleterious effect, but its role in anti-EGFR resistance should be further explored before a solid conclusion can be reached.
Similarly to RAS genes, BRAF and PIK3CA are driver oncogenes involved in colorectal carcinogenesis. It has been reported that BRAF mutations implying a high kinase activity (such as V600E) confer a poor prognosis, whereas those implying a low kinase activity (such as those located in codons 594 and 596) confer a favorable prognosis [29,30,31]. Our results reinforce the differential prognostic role of BRAF mutations as we only detected mutations implying a high kinase activity (V600E) in patients who were resistant to treatment, and we only found mutations causing low kinase activity in sensitive patients (D594N and G466A). As for location, several studies show that V600E mutations are more common in right colon cancers [38,39]. Accordingly, we only found V600E mutations in right-sided tumors. In contrast, we found D594N and G466A mutations in left-sided tumors. This differential distribution of BRAF mutations may affect the prognosis of left-sided tumors vs. right-sided tumors and define a clinically distinct subtype of CRC with an excellent prognosis. In addition, one resistant patient harbored NRAS and BRAF mutations over 5%, indicating that they are not always mutually exclusive. Regarding PIK3CA, the role of its activating mutations on anti-EGFR response remains under discussion [14,40,41]. In the present study, two patients, one responder and one non-responder, harbored the PIK3CA E545K mutation with an allelic fraction over 5%. Our results therefore strengthen the notion that this mutation is not critical for anti-EGFR response.
Our findings suggest some insulin-related genes (IGF1R and IRS2) have a substantial influence on anti-EGFR outcomes. A growing amount of evidence indicates that the insulin-like growth factor-1 receptor (IGF1R) is frequently overexpressed in CRC and that its activation is related to poorer outcomes [42]. IGF1R is a tyrosine kinase receptor that can also activate the RAS pathway and promote proliferation of cancer cells, resistance to apoptosis, and epithelial-mesenchymal transition [43]. Furthermore, it has been reported that a functional IGF1R receptor is required for EGFR-mediated growth and transformation [44] and that IGF1R expression modulates anti-EGFR efficacy in mCRC patients [45]. Consequently, in the same way that activating mutations in the EGFR receptor potentiate the RAS pathway, IGF1R activating mutations could play a similar role, with worse responses to EGFR-targeted therapies. To our knowledge, this is the first time that the missense IGF1R variants identified in our study (I668N and E1218K) have been related to anti-EGFR resistance. The variant E1218K has been previously described in a patient with leiomyosarcoma [33], but the variant I668N has not been reported previously [26]. We also identified a missense variant in the IRS2 gene (T1156M) associated with a lack of response to anti-EGFR drugs. Interestingly, Bertotti et al. demonstrated that IRS2 knockdown reduced sensitivity to cetuximab [23]. Truncating variants in this gene could therefore imply a lack of response to this drug. In our study, we found a missense variant in this gene harbored by three non-responders (Table 4). This variant had already been reported in a patient with colon cancer [33]. Like EGFR and IGF1R signaling pathways, the PDGFRA pathway is also involved in cell proliferation and migration [46]. The platelet-derived growth factor receptor A (PDGFRA) is a tyrosine kinase receptor that is frequently mutated in gastrointestinal stromal tumors [47]. In the present study, we found a splicing variant located in PDGFRA. Mutations in this gene have been correlated with a lack of response to anti-EGFR agents, although evidence is still scarce. In this line, Bertotti et al. described novel missense mutations located in/near the catalytic domain as mechanisms of anti-EGFR resistance [23].
We also found that LRIG1-3 genes could play a role in anti-EGFR response. Little is known about the role of somatic variants in these genes in responsiveness to anti-EGFR agents. Several studies have demonstrated that LRIG1 acts as a tumor suppressor by down-regulating ErbB and Met receptors, including EGFR [48,49,50,51,52]. In contrast, no definitive conclusions concerning the contribution of LRIG2 and LRIG3 to EGFR levels can be drawn as results reported to date are contradictory [48,53,54,55,56]. Gelfo et al. found that a lower LRIG1 expression predicted resistance to cetuximab therapy in CRC xenopatients, an effect that was not observed with LRIG3 [57]. In our study, one somatic variant in each LRIG gene was significantly associated with resistance to anti-EGFR blockade. Interestingly, the alteration located in LRIG1 had already been described in two cancer patients, one with a colon adenocarcinoma and the other with bladder cancer [34]. Conversely, the variants found in LRIG2 and LRIG3 have not been previously reported [26].
We wish to emphasize that no somatic variants related to good response were found in our study. This result is not striking as resistant mCRC phenotypes tend to be more heterogeneous than sensitive phenotypes [58]. Our results are in keeping with most papers to date describing mutations of non-response to anti-EGFR antibodies [59,60,61]. Only three of the eight somatic variants identified had been previously described [26], strengthening the extreme phenotype approach as a useful strategy to identify rare alterations. However, our study has some limitations. First, the two cohorts of patients are relatively small. The low number of patients showing an extremely poor response could be the result of good selection by clinicians after assessing RAS and BRAF mutational status and considering patient’s characteristics. The lack of statistically significant differences in the number of somatic mutations between the two cohorts of patients and between sidedness and drug response could also be attributable to the small sample size. Likewise, the association observed between sex and responsiveness could be a false positive predictive marker. Second, patients were included over a five-year period (2012–2017), during which time scientific knowledge regarding RAS and BRAF mutations increased significantly. This could explain why some patients with a BRAF V600E mutation already detected in the first assessment received an anti-EGFR agent. Larger, prospective and functional studies are needed to confirm the validity of our findings.
Therapies targeting different members and downstream effectors of the EGFR signaling pathway, such as KRAS or BRAF, have previously been explored [62,63]. However, most patients develop resistance to these agents. Our findings reveal novel mechanisms of resistance to anti-EGFR targeted therapies in mCRC patients that could provide new avenues for therapeutic intervention involving insulin-related and/or LRIG proteins. The utility of IGF1R as a possible therapeutic target has already been evaluated. In non-small cell lung cancer, IGF1R hyperactivity is related to acquired resistance to erlotinib, and simultaneous inhibition of EGFR and IGF1R is effective to prevent and also to overcome erlotinib resistance [64]. In CRC, the simultaneous inhibition of EGFR together with the IGF1R antagonist dalotuzumab has been tested in KRAS wild-type mCRC patients showing negative results [65]. Our results, together with the evidences thus far reported, could promote studies assessing insulin-related and also LRIG1-3 proteins as possible therapeutic targets in CRC. The impact of somatic variants in the genes encoding these proteins on anti-EGFR efficacy should also be confirmed.

5. Conclusions

In conclusion, this study shows that novel genetic variants along the EGFR-triggered pathway could affect the response to anti-EGFR drugs in mCRC patients. Our findings may help to better identify patients who are resistant to anti-EGFR drugs despite harboring RAS/BRAF-wild-type tumors. The eight genetic variants predictive of non-response could help guide clinical decision-making and improve outcomes by tailoring anti-EGFR drugs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/8/2245/s1, Figure S1: KRAS, NRAS, BRAF, and PIK3CA mutations, visualized with the Integrative Genomics Viewer (IGV), Figure S2: Genetic variations significantly associated with lack of response to anti-EGFR blockade, visualized with the Integrative Genomics Viewer (IGV).

Author Contributions

Conceptualization, P.R., J.S. (Jordi Surrallés), and D.P.; Formal analysis, P.R., B.R.-S., A.L., L.G.-Q., B.M. and D.P.; Funding acquisition, J.S. (Jordi Surrallés) and D.P.; Investigation, P.R., L.G.-Q. and C.C.; Methodology, P.R., B.R.-S., J.S. (Jordi Surrallés) and D.P.; Project administration, J.M., J.S. (Jordi Surrallés) and D.P.; Resources, P.R., B.R.-S., A.L., L.G.-Q., B.M., J.S. (Juliana Salazar), A.S., A.C.V., C.C. and D.P.; Software, B.R.-S.; Supervision, J.M., J.S. (Jordi Surrallés) and D.P.; Validation, P.R.; Writing–original draft, P.R.; Writing–review & editing, P.R., B.R.-S., A.L., L.G.-Q., B.M., J.S. (Juliana Salazar), A.S., A.C.V., J.M., C.C., J.S. (Jordi Surrallés) and D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Asociación Española Contra el Cáncer (Beca Campus CIMA de Sanitas Hospitales) and the Instituto de Salud Carlos III (grant CM18/00207 to Pau Riera]. Jordi Surrallés is supported by ICREA Academia Programme.

Acknowledgments

The authors thank Carolyn Newey for English language editing.

Conflicts of Interest

Ana Sebio reports receiving travel grants from Merck, Amgen, Sanofi, and Roche. Anna C. Virgili has honoraria from Speakers Bureau of Amgen, Sanofi and Kyowa Kirin, declares a scientific advisory role for Roche and Amgen, and reports receiving travel grants from Merck, Roche, Amgen, Sanofi, MSD and Servier. David Páez has honoraria from Speakers Bureau of Merck Serono and F. Hoffmann-La Roche Ltd., and declares a scientific advisory role for Amgen and Sanofi. All remaining authors have declared no conflicts of interest.

References

  1. Saif, M.W. Colorectal cancer in review: The role of the EGFR pathway. Expert Opin. Investig. Drugs 2010, 19, 357–369. [Google Scholar] [CrossRef] [PubMed]
  2. Normanno, N.; De Luca, A.; Bianco, C.; Strizzi, L.; Mancino, M.; Maiello, M.R.; Carotenuto, A.; De Feo, G.; Caponigro, F.; Salomon, D.S. Epidermal growth factor receptor (EGFR) signaling in cancer. Gene 2006, 366, 2–16. [Google Scholar] [CrossRef] [PubMed]
  3. Citri, A.; Yarden, Y. EGF-ERBB signalling: Towards the systems level. Nat. Rev. Mol. Cell Biol 2006, 7, 505–516. [Google Scholar] [CrossRef]
  4. Avraham, R.; Yarden, Y. Feedback regulation of EGFR signalling: Decision making by early and delayed loops. Nat. Rev. Mol. Cell Biol. 2011, 12, 104–117. [Google Scholar] [CrossRef] [PubMed]
  5. Douillard, J.Y.; Oliner, K.S.; Siena, S.; Tabernero, J.; Burkes, R.; Barugel, M.; Humblet, Y.; Bodoky, G.; Cunningham, D.; Jassem, J.; et al. Panitumumab-FOLFOX4 treatment and RAS mutations in colorectal cancer. N. Engl. J. Med. 2013, 369, 1023–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Van Cutsem, E.; Köhne, C.H.; Hitre, E.; Zaluski, J.; Chang Chien, C.R.; Makhson, A.; D’Haens, G.; Pintér, T.; Lim, R.; Bodoky, G.; et al. Cetuximab and chemotherapy as initial treatment for metastatic colorectal cancer. N. Engl. J. Med. 2009, 360, 1408–1417. [Google Scholar] [CrossRef] [Green Version]
  7. Allegra, C.J.; Rumble, R.B.; Hamilton, S.R.; Mangu, P.B.; Roach, N.; Hantel, A.; Schilsky, R.L. Extended RAS Gene Mutation Testing in Metastatic Colorectal Carcinoma to Predict Response to Anti-Epidermal Growth Factor Receptor Monoclonal Antibody Therapy: American Society of Clinical Oncology Provisional Clinical Opinion Update 2015. J. Clin. Oncol. 2015, 34, 179–185. [Google Scholar] [CrossRef] [PubMed]
  8. Sorich, M.J.; Wiese, M.D.; Rowland, A.; Kichenadasse, G.; McKinnon, R.A.; Karapetis, C.S. Extended RAS mutations and anti-EGFR monoclonal antibody survival benefit in metastatic colorectal cancer: A meta-analysis of randomized, controlled trials. Ann. Oncol. 2015, 26, 13–21. [Google Scholar] [CrossRef] [PubMed]
  9. Ali, M.; Kaltenbrun, E.; Anderson, G.R.; Stephens, S.J.; Arena, S.; Bardelli, A.; Counter, C.M.; Wood, K.C. Codon bias imposes a targetable limitation on KRAS-driven therapeutic resistance. Nat. Commun. 2017, 8, 15617. [Google Scholar] [CrossRef] [Green Version]
  10. Pietrantonio, F.; Petrelli, F.; Coinu, A.; Di Bartolomeo, M.; Borgonovo, K.; Maggi, C.; Cabiddu, M.; Iacovelli, R.; Bossi, I.; Lonati, V.; et al. Predictive role of BRAF mutations in patients with advanced colorectal cancer receiving cetuximab and panitumumab: A meta-analysis. Eur. J. Cancer 2015, 51, 587–594. [Google Scholar] [CrossRef]
  11. Rowland, A.; Dias, M.M.; Wiese, M.D.; Kichenadasse, G.; McKinnon, R.A.; Karapetis, C.S.; Sorich, M.J. Meta-analysis of BRAF mutation as a predictive biomarker of benefit from anti-EGFR monoclonal antibody therapy for RAS wild-type metastatic colorectal cancer. Br. J. Cancer 2015, 112, 1888–1894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. van Brummelen, E.M.J.; de Boer, A.; Beijnen, J.H.; Schellens, J.H.M. BRAF Mutations as Predictive Biomarker for Response to Anti-EGFR Monoclonal Antibodies. Oncologist 2017, 22, 864–872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Therkildsen, C.; Bergmann, T.K.; Henrichsen-Schnack, T.; Ladelund, S.; Nilbert, M. The predictive value of KRAS, NRAS, BRAF, PIK3CA and PTEN for anti-EGFR treatment in metastatic colorectal cancer: A systematic review and meta-analysis. Acta Oncol. 2014, 53, 852–864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Karapetis, C.S.; Jonker, D.; Daneshmand, M.; Hanson, J.E.; O’Callaghan, C.J.; Marginean, C.; Zalcberg, J.R.; Simes, J.; Moore, M.J.; Tebbutt, N.C.; et al. PIK3CA, BRAF, and PTEN status and benefit from cetuximab in the treatment of advanced colorectal cancer—results from NCIC CTG/AGITG CO.17. Clin. Cancer Res. 2014, 20, 744–753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Arnold, D.; Lueza, B.; Douillard, J.Y.; Peeters, M.; Lenz, H.J.; Venook, A.; Heinemann, V.; Van Cutsem, E.; Pignon, J.P.; Tabernero, J.; et al. Prognostic and predictive value of primary tumour side in patients with RAS wild-type metastatic colorectal cancer treated with chemotherapy and EGFR directed antibodies in six randomized trials. Ann. Oncol. 2017, 28, 1713–1729. [Google Scholar] [CrossRef] [PubMed]
  16. Moretto, R.; Cremolini, C.; Rossini, D.; Pietrantonio, F.; Battaglin, F.; Mennitto, A.; Bergamo, F.; Loupakis, F.; Marmorino, F.; Berenato, R.; et al. Location of Primary Tumor and Benefit From Anti-Epidermal Growth Factor Receptor Monoclonal Antibodies in Patients With RAS and BRAF Wild-Type Metastatic Colorectal Cancer. Oncologist 2016, 21, 988–994. [Google Scholar] [CrossRef] [Green Version]
  17. De Sousa EMelo, F.; Wang, X.; Jansen, M.; Fessler, E.; Trinh, A.; de Rooij, L.P.; de Jong, J.H.; de Boer, O.J.; van Leersum, R.; Bijlsma, M.F.; et al. Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions. Nat. Med. 2013, 19, 614–618. [Google Scholar] [CrossRef] [PubMed]
  18. Linnekamp, J.F.; Wang, X.; Medema, J.P.; Vermeulen, L. Colorectal cancer heterogeneity and targeted therapy: A case for molecular disease subtypes. Cancer Res. 2015, 75, 245–249. [Google Scholar] [CrossRef] [Green Version]
  19. Khattak, M.A.; Martin, H.; Davidson, A.; Phillips, M. Role of first-line anti-epidermal growth factor receptor therapy compared with anti-vascular endothelial growth factor therapy in advanced colorectal cancer: A meta-analysis of randomized clinical trials. Clin. Colorectal Cancer 2015, 14, 81–90. [Google Scholar] [CrossRef] [PubMed]
  20. Lee, M.S.; Kopetz, S. Current and Future Approaches to Target the Epidermal Growth Factor Receptor and Its Downstream Signaling in Metastatic Colorectal Cancer. Clin. Colorectal Cancer 2015, 14, 203–218. [Google Scholar] [CrossRef]
  21. Eisenhauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.H.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M.; et al. New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1). Eur. J. Cancer 2009, 45, 228–247. [Google Scholar] [CrossRef] [PubMed]
  22. Bertotti, A.; Sassi, F. Molecular Pathways: Sensitivity and Resistance to Anti-EGFR Antibodies. Clin. Cancer Res. 2015, 21, 3377–3383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bertotti, A.; Papp, E.; Jones, S.; Adleff, V.; Anagnostou, V.; Lupo, B.; Sausen, M.; Phallen, J.; Hruban, C.A.; Tokheim, C.; et al. The genomic landscape of response to EGFR blockade in colorectal cancer. Nature 2015, 526, 263–267. [Google Scholar] [CrossRef] [PubMed]
  24. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  25. Cibulskis, K.; Lawrence, M.S.; Carter, S.L.; Sivachenko, A.; Jaffe, D.; Sougnez, C.; Gabriel, S.; Meyerson, M.; Lander, E.S.; Getz, G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 2013, 31, 213–219. [Google Scholar] [CrossRef]
  26. Tate, J.G.; Bamford, S.; Jubb, H.C.; Sondka, Z.; Beare, D.M.; Bindal, N.; Boutselakis, H.; Cole, C.G.; Creatore, C.; Dawson, E.; et al. COSMIC: The Catalogue Of Somatic Mutations In Cancer. Nucleic Acids Res. 2019, 47, D941–D947. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  28. Tamborero, D.; Rubio-Perez, C.; Deu-Pons, J.; Schroeder, M.P.; Vivancos, A.; Rovira, A.; Tusquets, I.; Albanell, J.; Rodon, J.; Tabernero, J.; et al. Cancer Genome Interpreter annotates the biological and clinical relevance of tumor alterations. Genome Med. 2018, 10, 25. [Google Scholar] [CrossRef]
  29. Jones, J.C.; Renfro, L.A.; Al-Shamsi, H.O.; Schrock, A.B.; Rankin, A.; Zhang, B.Y.; Kasi, P.M.; Voss, J.S.; Leal, A.D.; Sun, J.; et al. Non-V600BRAF Mutations Define a Clinically Distinct Molecular Subtype of Metastatic Colorectal Cancer. J. Clin. Oncol. 2017, 35, 2624–2630. [Google Scholar] [CrossRef]
  30. Cremolini, C.; Di Bartolomeo, M.; Amatu, A.; Antoniotti, C.; Moretto, R.; Berenato, R.; Perrone, F.; Tamborini, E.; Aprile, G.; Lonardi, S.; et al. BRAF codons 594 and 596 mutations identify a new molecular subtype of metastatic colorectal cancer at favorable prognosis. Ann. Oncol. 2015, 26, 2092–2097. [Google Scholar] [CrossRef]
  31. Yao, Z.; Yaeger, R.; Rodrik-Outmezguine, V.S.; Tao, A.; Torres, N.M.; Chang, M.T.; Drosten, M.; Zhao, H.; Cecchi, F.; Hembrough, T.; et al. Tumours with class 3 BRAF mutants are sensitive to the inhibition of activated RAS. Nature 2017, 548, 234–238. [Google Scholar] [CrossRef] [PubMed]
  32. Dankner, M.; Rose, A.A.N.; Rajkumar, S.; Siegel, P.M.; Watson, I.R. Classifying BRAF alterations in cancer: New rational therapeutic strategies for actionable mutations. Oncogene 2018, 37, 3183–3199. [Google Scholar] [CrossRef] [PubMed]
  33. Zehir, A.; Benayed, R.; Shah, R.H.; Syed, A.; Middha, S.; Kim, H.R.; Srinivasan, P.; Gao, J.; Chakravarty, D.; Devlin, S.M.; et al. Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat. Med. 2017, 23, 703–713. [Google Scholar] [CrossRef] [PubMed]
  34. Giannakis, M.; Mu, X.J.; Shukla, S.A.; Qian, Z.R.; Cohen, O.; Nishihara, R.; Bahl, S.; Cao, Y.; Amin-Mansour, A.; Yamauchi, M.; et al. Genomic Correlates of Immune-Cell Infiltrates in Colorectal Carcinoma. Cell Rep. 2016, 15, 857–865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sabour, L.; Sabour, M.; Ghorbian, S. Clinical Applications of Next-Generation Sequencing in Cancer Diagnosis. Pathol. Oncol. Res. 2017, 23, 225–234. [Google Scholar] [CrossRef] [PubMed]
  36. Santos, C.; Azuara, D.; Viéitez, J.M.; Páez, D.; Falcó, E.; Élez, E.; López-López, C.; Valladares, M.; Robles-Díaz, L.; García-Alfonso, P.; et al. Phase II study of high-sensitivity genotyping of KRAS, NRAS, BRAF and PIK3CA to ultra-select metastatic colorectal cancer patients for panitumumab plus FOLFIRI: The ULTRA trial. Ann. Oncol. 2019, 30, 796–803. [Google Scholar] [CrossRef] [Green Version]
  37. Vidal, J.; Bellosillo, B.; Santos Vivas, C.; García-Alfonso, P.; Carrato, A.; Cano, M.T.; García-Carbonero, R.; Élez, E.; Losa, F.; Massutí, B.; et al. Ultra-selection of metastatic colorectal cancer patients using next-generation sequencing to improve clinical efficacy of anti-EGFR therapy. Ann. Oncol. 2019, 30, 439–446. [Google Scholar] [CrossRef]
  38. Isnaldi, E.; Garuti, A.; Cirmena, G.; Scabini, S.; Rimini, E.; Ferrando, L.; Lia, M.; Murialdo, R.; Tixi, L.; Carminati, E.; et al. Clinico-pathological associations and concomitant mutations of the RAS/RAF pathway in metastatic colorectal cancer. J. Transl. Med. 2019, 17, 137. [Google Scholar] [CrossRef]
  39. Salem, M.E.; Weinberg, B.A.; Xiu, J.; El-Deiry, W.S.; Hwang, J.J.; Gatalica, Z.; Philip, P.A.; Shields, A.F.; Lenz, H.J.; Marshall, J.L. Comparative molecular analyses of left-sided colon, right-sided colon, and rectal cancers. Oncotarget 2017, 8, 86356–86368. [Google Scholar] [CrossRef] [Green Version]
  40. Prenen, H.; De Schutter, J.; Jacobs, B.; De Roock, W.; Biesmans, B.; Claes, B.; Lambrechts, D.; Van Cutsem, E.; Tejpar, S. PIK3CA mutations are not a major determinant of resistance to the epidermal growth factor receptor inhibitor cetuximab in metastatic colorectal cancer. Clin. Cancer Res. 2009, 15, 3184–3188. [Google Scholar] [CrossRef] [Green Version]
  41. De Roock, W.; Claes, B.; Bernasconi, D.; De Schutter, J.; Biesmans, B.; Fountzilas, G.; Kalogeras, K.T.; Kotoula, V.; Papamichael, D.; Laurent-Puig, P.; et al. Effects of KRAS, BRAF, NRAS, and PIK3CA mutations on the efficacy of cetuximab plus chemotherapy in chemotherapy-refractory metastatic colorectal cancer: A retrospective consortium analysis. Lancet Oncol. 2010, 11, 753–762. [Google Scholar] [CrossRef]
  42. Lee, J.; Jain, A.; Kim, P.; Lee, T.; Kuller, A.; Princen, F.; Gu, I.; Kim, S.H.; Park, J.O.; Park, Y.S.; et al. Activated cMET and IGF1R-driven PI3K signaling predicts poor survival in colorectal cancers independent of KRAS mutational status. PLoS ONE 2014, 9, e103551. [Google Scholar] [CrossRef] [PubMed]
  43. Dallas, N.A.; Xia, L.; Fan, F.; Gray, M.J.; Gaur, P.; van Buren, G., 2nd; Samuel, S.; Kim, M.P.; Lim, S.J.; Ellis, L.M.; et al. Chemoresistant colorectal cancer cells, the cancer stem cell phenotype, and increased sensitivity to insulin-like growth factor-I receptor inhibition. Cancer Res. 2009, 69, 1951–1957. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Coppola, D.; Ferber, A.; Miura, M.; Sell, C.; D’Ambrosio, C.; Rubin, R.; Baserga, R. A functional insulin-like growth factor I receptor is required for the mitogenic and transforming activities of the epidermal growth factor receptor. Mol. Cell. Biol. 1994, 14, 4588–4595. [Google Scholar] [CrossRef] [Green Version]
  45. Alonso, V.; Escudero, P.; Fernández-Martos, C.; Salud, A.; Méndez, M.; Gallego, J.; Rodriguez, J.R.; Martín-Richard, M.; Fernández-Plana, J.; Manzano, H.; et al. Coexpression of p-IGF-1R and MMP-7 Modulates Panitumumab and Cetuximab Efficacy in RAS Wild-Type Metastatic Colorectal Cancer Patients. Neoplasia 2018, 20, 678–686. [Google Scholar] [CrossRef] [PubMed]
  46. Hayashi, Y.; Bardsley, M.R.; Toyomasu, Y.; Milosavljevic, S.; Gajdos, G.B.; Choi, K.M.; Reid-Lombardo, K.M.; Kendrick, M.L.; Bingener-Casey, J.; Tang, C.M.; et al. Platelet-Derived Growth Factor Receptor-α Regulates Proliferation of Gastrointestinal Stromal Tumor Cells With Mutations in KIT by Stabilizing ETV1. Gastroenterology 2015, 149, 420–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Heinrich, M.C.; Corless, C.L.; Duensing, A.; McGreevey, L.; Chen, C.J.; Joseph, N.; Singer, S.; Griffith, D.J.; Haley, A.; Town, A.; et al. PDGFRA activating mutations in gastrointestinal stromal tumors. Science 2003, 299, 708–810. [Google Scholar] [CrossRef]
  48. Rafidi, H.; Mercado, F., 3rd; Astudillo, M.; Fry, W.H.; Saldana, M.; Carraway, K.L., 3rd; Sweeney, C. Leucine-rich repeat and immunoglobulin domain-containing protein-1 (Lrig1) negative regulatory action toward ErbB receptor tyrosine kinases is opposed by leucine-rich repeat and immunoglobulin domain-containing protein 3 (Lrig3). J. Biol. Chem. 2013, 288, 21593–21605. [Google Scholar] [CrossRef] [Green Version]
  49. Faraz, M.; Herdenberg, C.; Holmlund, C.; Henriksson, R.; Hedman, H. A protein interaction network centered on leucine-rich repeats and immunoglobulin-like domains 1 (LRIG1) regulates growth factor receptors. J. Biol. Chem. 2018, 293, 3421–3435. [Google Scholar] [CrossRef] [Green Version]
  50. Wang, Y.; Poulin, E.J.; Coffey, R.J. LRIG1 is a triple threat: ERBB negative regulator, intestinal stem cell marker and tumour suppressor. Br. J. Cancer 2013, 108, 1765–1770. [Google Scholar] [CrossRef]
  51. Wong, V.W.; Stange, D.E.; Page, M.E.; Buczacki, S.; Wabik, A.; Itami, S.; van de Wetering, M.; Poulsom, R.; Wright, N.A.; Trotter, M.W.; et al. Lrig1 controls intestinal stem-cell homeostasis by negative regulation of ErbB signalling. Nat. Cell. Biol. 2012, 14, 401–408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Powell, A.E.; Wang, Y.; Li, Y.; Poulin, E.J.; Means, A.L.; Washington, M.K.; Higginbotham, J.N.; Juchheim, A.; Prasad, N.; Levy, S.E.; et al. The pan-ErbB negative regulator Lrig1 is an intestinal stem cell marker that functions as a tumor suppressor. Cell 2012, 149, 146–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wang, B.; Han, L.; Chen, R.; Cai, M.; Han, F.; Lei, T.; Guo, D. Downregulation of LRIG2 expression by RNA interference inhibits glioblastoma cell (GL15) growth, causes cell cycle redistribution, increases cell apoptosis and enhances cell adhesion and invasion in vitro. Cancer Biol. Ther. 2009, 8, 1018–1023. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yang, H.-K.; Chen, H.; Mao, F.; Xiao, Q.-G.; Xie, R.-F.; Lei, T. Downregulation of LRIG2 expression inhibits angiogenesis of glioma via EGFR/VEGF-A pathway. Oncol. Lett. 2017, 14, 4021–4028. [Google Scholar] [CrossRef] [PubMed]
  55. Wu, X.; Hedman, H.; Bergqvist, M.; Bergström, S.; Henriksson, R.; Gullbo, J.; Lennartsson, J.; Hesselius, P.; Ekman, S. Expression of EGFR and LRIG proteins in oesophageal carcinoma with emphasis on patient survival and cellular chemosensitivity. Acta Oncol. 2012, 51, 69–76. [Google Scholar] [CrossRef] [Green Version]
  56. Guo, D.; Yang, H.; Guo, Y.; Xiao, Q.; Mao, F.; Tan, Y.; Wan, X.; Wang, B.; Lei, T. LRIG3 modulates proliferation, apoptosis and invasion of glioblastoma cells as a potent tumor suppressor. J. Neurol. Sci. 2015, 350, 61–68. [Google Scholar] [CrossRef]
  57. Gelfo, V.; Pontis, F.; Mazzeschi, M.; Sgarzi, M.; Mazzarini, M.; Solmi, R.; D’Uva, G.; Lauriola, M. Glucocorticoid Receptor Modulates EGFR Feedback upon Acquisition of Resistance to Monoclonal Antibodies. J. Clin. Med. 2019, 8, 600. [Google Scholar] [CrossRef] [Green Version]
  58. Dagogo-Jack, I.; Shaw, A.T. Tumour heterogeneity and resistance to cancer therapies. Nat. Rev. Clin. Oncol. 2018, 15, 81–94. [Google Scholar] [CrossRef]
  59. Hammond, W.A.; Swaika, A.; Mody, K. Pharmacologic resistance in colorectal cancer: A review. Ther. Adv. Med. Oncol. 2016, 8, 57–84. [Google Scholar] [CrossRef] [Green Version]
  60. Bignucolo, A.; De Mattia, E.; Cecchin, E.; Roncato, R.; Toffoli, G. Pharmacogenomics of Targeted Agents for Personalization of Colorectal Cancer Treatment. Int. J. Mol. Sci. 2017, 18, 1522. [Google Scholar] [CrossRef] [Green Version]
  61. Cremolini, C.; Morano, F.; Moretto, R.; Berenato, R.; Tamborini, E.; Perrone, F.; Rossini, D.; Gloghini, A.; Busico, A.; Zucchelli, G.; et al. Negative hyper-selection of metastatic colorectal cancer patients for anti-EGFR monoclonal antibodies: The PRESSING case–control study. Ann. Oncol. 2017, 28, 3009–3014. [Google Scholar] [CrossRef]
  62. Canon, J.; Rex, K.; Saiki, A.Y.; Mohr, C.; Cooke, K.; Bagal, D.; Gaida, K.; Holt, T.; Knutson, C.G.; Koppada, N.; et al. The clinical KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. Nature 2019, 575, 217–223. [Google Scholar] [CrossRef] [PubMed]
  63. Kopetz, S.; Grothey, A.; Yaeger, R.; Van Cutsem, E.; Desai, J.; Yoshino, T.; Wasan, H.; Ciardiello, F.; Loupakis, F.; Hong, Y.S.; et al. Encorafenib, Binimetinib, and Cetuximab in BRAF V600E-Mutated Colorectal Cancer. N. Engl. J. Med. 2019, 381, 1632–1643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Vazquez-Martin, A.; Cufí, S.; Oliveras-Ferraros, C.; Torres-Garcia, V.Z.; Corominas-Faja, B.; Cuyàs, E.; Bonavia, R.; Visa, J.; Martin-Castillo, B.; Barrajón-Catalán, E.; et al. IGF-1R/epithelial-to-mesenchymal transition (EMT) crosstalk suppresses the erlotinib-sensitizing effect of EGFR exon 19 deletion mutations. Sci. Rep. 2013, 3, 2560. [Google Scholar] [CrossRef] [PubMed]
  65. Sclafani, F.; Kim, T.Y.; Cunningham, D.; Kim, T.W.; Tabernero, J.; Schmoll, H.J.; Roh, J.K.; Kim, S.Y.; Park, Y.S.; Guren, T.K.; et al. A Randomized Phase II/III Study of Dalotuzumab in Combination With Cetuximab and Irinotecan in Chemorefractory, KRAS Wild-Type, Metastatic Colorectal Cancer. J. Natl. Cancer Inst. 2015, 107, djv258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 12 02245 g001 550
Figure 1. Flow chart of patient selection. * Inclusion criteria: patients ≥ 18 years, with an Eastern cooperative oncology group (ECOG) performance status ≤ 2 and RAS wild-type tumors. Patients had to be sensitive (patients achieving complete or partial response at the first CT-scan or disease stabilization lasting at least 6 months) or resistant to anti-EGFR blockade (patients showing disease progression at the first CT-scan).
Figure 1. Flow chart of patient selection. * Inclusion criteria: patients ≥ 18 years, with an Eastern cooperative oncology group (ECOG) performance status ≤ 2 and RAS wild-type tumors. Patients had to be sensitive (patients achieving complete or partial response at the first CT-scan or disease stabilization lasting at least 6 months) or resistant to anti-EGFR blockade (patients showing disease progression at the first CT-scan).
Cancers 12 02245 g001
Cancers 12 02245 g002 550
Figure 2. Number of somatic genetic variants and response to anti-EGFR agents.
Figure 2. Number of somatic genetic variants and response to anti-EGFR agents.
Cancers 12 02245 g002
Table
Table 1. Selected genes classified according to the function of their encoding proteins.
Table 1. Selected genes classified according to the function of their encoding proteins.
LigandsReceptorsIntracellular Downstream EffectorsProteins Involved in EGFR TurnoverOthers
AREGEGFR (ERBB1 or HER1)AKT1AGR2TP53
BTCERBB2 (HER2)BRAFCBLYAP1
EGFERBB3 (HER3)HRASLRIG1
EPGNERBB4 (HER4)IRS2LRIG2
EREGFGFR1KRASLRIG3
HBEGFIGF1RMAP2K1NEDD8
HGFMETNRASERRFI1 (RALT or MIG6)
IGF1-2PDGFRAPIK3CASOCS4
NRG1-4 PTENSOCS5
TGFα SPRY2
Abbreviations: AGR2, Anterior Gradient 2, Protein Disulphide Isomerase Family Member; AKT1, AKT Serine/Threonine Kinase 1; AREG, Amphiregulin; BRAF, B-Raf Proto-Oncogene, Serine/Threonine Kinase; BTC, Betacellulin; CBL, Cbl Proto-Oncogene; EGF, Epidermal Growth Factor; EGFR, Epidermal Growth Factor Receptor; EPGN, Epithelial Mitogen; ERBB2, Erb-B2 Receptor Tyrosine Kinase 2; ERBB3, Erb-B3 Receptor Tyrosine Kinase 3; ERBB4, Erb-B4 Receptor Tyrosine Kinase 4; EREG, Epiregulin; ERRFI1, ERBB Receptor Feedback Inhibitor 1; FGFR1, Fibroblast Growth Factor Receptor 1; HBEGF, Heparin Binding EGF Like Growth Factor; HGF, Hepatocyte Growth Factor; HRAS, HRas Proto-Oncogene, GTPase; IGF1, Insulin Like Growth Factor 1; IGF2, Insulin Like Growth Factor 2; IGF1R, Insulin Like Growth Factor 1 Receptor; IRS2, Insulin Receptor Substrate 2; KRAS, KRAS Proto-Oncogene, GTPase; MAP2K1, Mitogen-Activated Protein Kinase Kinase 1; LRIG1, Leucine Rich Repeats And Immunoglobulin Like Domains 1; LRIG2, Leucine Rich Repeats And Immunoglobulin Like Domains 2; LRIG3, Leucine Rich Repeats And Immunoglobulin Like Domains 3; MET, MET Proto-Oncogene, Receptor Tyrosine Kinase; NEDD8, Neural Precursor Cell Expressed, Developmentally Down-Regulated 8; NRAS, NRAS Proto-Oncogene, GTPase; NRG1, Neuregulin 1; NRG2, Neuregulin 2; NRG3, Neuregulin 3; NRG4, Neuregulin 4; PDGFRA, Platelet Derived Growth Factor Receptor Alpha; PIK3CA, Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha; PTEN, Phosphatase And Tensin Homolog; SOCS4, Suppressor Of Cytokine Signaling 4; SOCS5, Suppressor Of Cytokine Signaling 5; SPRY2, Sprouty RTK Signaling Antagonist 2; TGFα, Transforming Growth Factor Alpha; TP53, Tumor Protein P53; YAP1, Yes Associated Protein 1.
Table
Table 2. Baseline patient characteristics.
Table 2. Baseline patient characteristics.
CharacteristicStudy Population
(n = 61)		Sensitive Patients
(n = 38)		Resistant Patients
(n = 23)		p-Value
N (%)N (%)N (%)
Sex
Male40 (65.6%)32 (84%)8 (34.8%)<0.001
Female21 (34.4%)6 (16%)15 (65.2%)
Age
<7543 (70.5%)26 (68.4%)17 (73.9%)0.649
≥7518 (29.5%)12 (31.6%)6 (26.1%)
Mean Age 67.766.9
Performance status (ECOG)
031 (50.8%)24 (63.2%)7 (30.4%)0.013
1–230 (49.2%)14 (36.8%)16 (69.6%)
Smoking habit
Never smokers24 (39.3%)13 (34.2%)11 (47.8%)0.291
Current or former smokers37 (60.7%)25 (65.8%)12 (52.2%)
Tumor side
Right24 (39.4%)15 (39.5%)9 (39.2%)
Left20 (32.7%)13 (34.2%)7 (30.4 %)0.634
Rectal16 (26.3%)10 (26.3%)6 (26.1%)
Jejunum1 (1.6%)01 (4.3%)
Number of metastatic sites
130 (49.2%)18 (47.4%)12 (52.3%)0.716
≥231 (50.8%)20 (52.6%)11 (47.7%)
Time to metastases
Synchronous38 (62.3%)27 (71.1%)11 (47.8%)0.070
Metachronous23 (37.7%)11 (28.9%)12 (52.2%)
Primary resected
Yes49 (80.3%)32 (84.2%)17 (73.9%)0.327
No12 (19.7%)6 (15.8%)6 (26.1%)
Previous lines of treatment
024 (39.3%)18 (47.4%)6 (26.1%)
132 (52.5%)17 (44.7%)15 (65.2%)0.246
≥25 (8.2%)3 (7.9%)2 (8.7%)
Type of anti-EGFR
Cetuximab33 (54.1%)21 (55.3%)12 (52.2%)0.814
Panitumumab28 (45.9%)17 (44.7%)11 (47.8%)
Combination QT
FOLFOX18 (29.5%)14 (36.8 %)4 (17.4%)
Irinotecan scheme40 (65.6%)23 (60.6%)17 (73.9%)0.192
Monotherapy3 (4.9%)1 (2.6 %)2 (8.7%)
PFS (months) 18.84.7
OS (months) 41.217.2
Abbreviations: ECOG, Eastern Cooperative Oncology Group; OS, overall survival; PFS, progression-free survival. p-values below 0.05 are highlighted in bold.
Table
Table 3. Assessment of KRAS, NRAS, BRAF, and PIK3CA mutational status by next-generation sequencing.
Table 3. Assessment of KRAS, NRAS, BRAF, and PIK3CA mutational status by next-generation sequencing.
Treatment OutcomeGeneMutationPatient% of MutationCoverageMutational Status Prior to Anti-EGFR Prescription
ResistantKRASG12CP1145.7%2316KRAS wild-type
KRASQ61LP390.8%14,637KRAS exon 3 not tested
NRASG12SP512.6%5093NRAS wild-type
NRASG13DP558.1%12,265NRAS not tested
NRASG13DP570.5%41,758NRAS wild-type
BRAFV600EP324.4%10,995BRAF V600E mutated
BRAFV600EP3910.6%7906BRAF V600E mutated
BRAFV600EP5520.8%22,185BRAF V600E mutated
BRAFV600EP6412.1%1316BRAF V600E wild-type
PIK3CAE545KP6310.1%17,051PIK3CA not tested
SensitiveKRASA146VP285.0%2743KRAS wild-type
BRAFD594NP4519.7%23,584BRAF codon 594 not tested
BRAFG466AP5921.0%10,234BRAF codon 466 not tested
PIK3CAE545KP666.1%12,888PIK3CA not tested
Abbreviations: BRAF, B-Raf Proto-Oncogene, Serine/Threonine Kinase; KRAS, KRAS Proto-Oncogene, GTPase; NRAS, NRAS Proto-Oncogene, GTPase; PIK3CA, Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha.
Table
Table 4. Genetic variations significantly associated with lack of response to anti-EGFR blockade.
Table 4. Genetic variations significantly associated with lack of response to anti-EGFR blockade.
GeneGenetic VariantPatients with the Variant% of the Somatic VariantCoveragep- Value *Presence in COSMIC Cancer Database v91
IGF1RNM_001291858.1: c.2003T > A; p.(I668N)P43.5%80080.029Not described
P94.2%1745
P1010.6%3093
NM_001291858.1: c.3652G > A; p.(E1218K)P16.3%27400.008Mutation Id: 6919417
Patient with a leiomyosarcoma (n = 1)
[33]
P211.8%1866
P92.4%1238
P611.5%10,017
IRS2NM_003749.2: c.3467C > T; p.(T1156M)P127.0%7430.029Mutation Id: 6974893
Patient with colon cancer (n = 1)
[33]
P145.1%1093
P573.4%11,118
LRIG1NM_015541: c.456G > A; p.(T152T)P1289.3%1960.008Mutation Id: 4005617
Patient with colon cancer (n = 1)
Patient with bladder cancer (n = 1)
[34]
P572.7%12,073
P637.9%6106
P674.0%3218
LRIG2NM_014813.2: c.2090C > T; p.(S697L)P12.2%88020.029Not described
P104.2%1464
P127.9%1794
LRIG3NM_001136051.2: c.2434G > A; p.(V812M)P11.7%12,3780.029Not described
P24.8%4477
P94.6%2266
NRASNM_002524.3: c.344del; p.(G115Efs*46)P91.6%22730.029Not described
P103.1%3140
P632.2%20,507
PDGFRANM_001347828: c.903G > A; p.(T301T)P5222.0%19,0780.029Not described
P633.8%26,774
P671.2%55,158
* p-value was obtained by Fisher’s exact test. Abbreviations: IGF1R, Insulin-like Growth Factor 1 Receptor; IRS2, Insulin Receptor Substrate 2; LRIG1, Leucine-rich Repeats and Immunoglobulin-like Domains 1; LRIG2, Leucine-rich Repeats and Immunoglobulin-like Domains 2; LRIG3, Leucine-rich Repeats and Immunoglobulin-like Domains 3; NRAS, NRAS Proto-Oncogene, GTPase; PDGFRA, Platelet-derived Growth Factor Receptor Alpha.

Share and Cite

MDPI and ACS Style

Riera, P.; Rodríguez-Santiago, B.; Lasa, A.; Gonzalez-Quereda, L.; Martín, B.; Salazar, J.; Sebio, A.; Virgili, A.C.; Minguillón, J.; Camps, C.; et al. Novel Somatic Genetic Variants as Predictors of Resistance to EGFR-Targeted Therapies in Metastatic Colorectal Cancer Patients. Cancers 2020, 12, 2245. https://doi.org/10.3390/cancers12082245

AMA Style

Riera P, Rodríguez-Santiago B, Lasa A, Gonzalez-Quereda L, Martín B, Salazar J, Sebio A, Virgili AC, Minguillón J, Camps C, et al. Novel Somatic Genetic Variants as Predictors of Resistance to EGFR-Targeted Therapies in Metastatic Colorectal Cancer Patients. Cancers. 2020; 12(8):2245. https://doi.org/10.3390/cancers12082245

Chicago/Turabian Style

Riera, Pau, Benjamín Rodríguez-Santiago, Adriana Lasa, Lidia Gonzalez-Quereda, Berta Martín, Juliana Salazar, Ana Sebio, Anna C. Virgili, Jordi Minguillón, Cristina Camps, and et al. 2020. "Novel Somatic Genetic Variants as Predictors of Resistance to EGFR-Targeted Therapies in Metastatic Colorectal Cancer Patients" Cancers 12, no. 8: 2245. https://doi.org/10.3390/cancers12082245

APA Style

Riera, P., Rodríguez-Santiago, B., Lasa, A., Gonzalez-Quereda, L., Martín, B., Salazar, J., Sebio, A., Virgili, A. C., Minguillón, J., Camps, C., Surrallés, J., & Páez, D. (2020). Novel Somatic Genetic Variants as Predictors of Resistance to EGFR-Targeted Therapies in Metastatic Colorectal Cancer Patients. Cancers, 12(8), 2245. https://doi.org/10.3390/cancers12082245

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop