Expression of pEGFR and pAKT as response-predictive biomarkers for RAS wild-type patients to anti-EGFR monoclonal antibodies in metastatic colorectal cancers.

BACKGROUND: RAS wild-type (RASw/t) tumours have been associated with better outcomes in patients with metastatic colorectal cancer (mCRC) treated with anti-EGFR monoclonal antibodies (mAb). We investigated the expression of EGFR downstream proteins under their active phosphorylated forms as potential markers in response to these patients.
METHODS: One-hundred tumour samples were collected from patients with mCRC refractory to FOLFOX and/or FOLFIRI and treated by a combination of chemotherapy with anti-EGFR mAb. The outcomes were measured on response evaluation criteria in solid tumour (RECIST), progression-free survival (PFS) and overall survival (OS). All samples were assessed for RAS and BRAF mutations, and the key phosphorylated proteins of EGFR downstream signalling were quantitatively analysed using the BioPlex Protein array.
RESULTS: Among the 60 RASw/t patients, 45.0% achieved a complete or partial response when treated with anti-EGFR mAb. Expression of pAKT, pERK1/2 and pMEK1 was significantly lower in RASw/t patients (P=0.0246; P=0.004; P=0.0110, respectively). The response rate was significantly higher for RASw/t patients who express pEGFR and pAKT (P=0.0258; P=0.0277, respectively).
CONCLUSIONS: Overexpression of pEGFR and pAKT may predict the response rate in RASw/t patients treated with anti-EGFR mAb. On the basis of our results, we hypothesise that the association of anti-EGFR mAb and anti-AKT therapies could be of interest.

Colorectal cancer is the third most common cancer worldwide with more than one million patients diagnosed each year, among 50% of whom will develop metastatic disease (Jemal ; Tougeron ). Recent efforts to improve the treatment of metastatic colorectal cancer (mCRC) have led to the development of monoclonal antibodies (mAbs) such as cetuximab and panitumumab that inhibit the activation of the epidermal growth factor receptor (EGFR) and its downstream pathways (namely RAS/RAF/MAPK and PI3K/AKT/mTOR) that promote cell growth, proliferation, inhibition of apoptosis, invasion and metastasis. Evidence has first showed that patients with codon 12 or 13 KRAS-mutant tumours receive little or no benefit from anti-EGFR therapies as single agents or combined with chemotherapy (Di Fiore ; Khambata-Ford ; De Roock ; Lievre ). More recent studies revealed that the presence of mutations on exons 2, 3 and 4 of KRAS or NRAS were also predictive of resistance to anti-EGFR therapies (Douillard ; Sorich ). In other words, cetuximab or panitumumab shows benefits only to ‘RAS wild-type' patients, that is, patients with tumours that are not bearing mutations on exons 2, 3 and 4 of KRAS and NRAS. These findings led the European Medicine Agency to the use of cetuximab and panitumumab for patients with RAS wild-type tumours.

However, from studies including RAS wild-type patients, the response rates to cetuximab or panitumumab therapy ranged from only 40 to 60% (Linardou ), which results in a large fraction of patients being without any known causes for treatment failure. The presence of alterations in other genes such as PIK3CA or BRAF in the EGFR-dependent signalling pathways (Laurent-Puig ; Custodio and Feliu, 2013; Douillard ) is responsible for some of the non-responding cases. Moreover, overexpression or alterations on proteins such as PTEN, PI3K and AKT, involved in the RAS/RAF/MAPK or PI3K/AKT/mTOR signalling pathways, can have a significant impact on cell proliferation or apoptosis. Absence or overexpression of proteins under their phosphorylated forms may be of interest to predict response to anti-EGFR mAb. Ligand fixation of the EGF receptor activates the pathway. The fixation of ligand on the receptor activates the kinase activity that is responsible for the activation of downstream proteins by phosphorylation. The study of the expression of the downstream proteins under their phosphorylated forms is important to understand mechanisms of resistance to anticancer drugs and more precisely in our study of anti-EGFR drugs. In a previous study, we showed that the expression level of EGFR downstream signalling phosphoproteins (pMEK1 and pP70S6K) could predict longer progression-free survival (PFS) and overall survival (OS) in KRAS exon 2 wild-type mCRC using multivariate analysis (Perkins ).

In this paper, we describe the impact of RAS mutations (that is, KRAS and NRAS exon 2, 3 and 4 mutations) and of the expression of proteins involved in RAS/RAF/MAPK and PI3K/AKT/mTOR pathways in their active phosphorylated forms on the response to anti-EGFR mAbs in RAS wild-type patients with an advanced colorectal cancer. The 42 patient population described in our previous study (Perkins ) has been extended to 100 patients, and the research of mutations has been also extended from KRAS exon 2 to RAS (KRAS, NRAS) exons 2, 3 and 4.

Patients and Methods

Patients

We analysed total protein extracts prepared from 100 frozen tumour samples from patients with mCRC treated by a combination of chemotherapy with anti-EGFR therapy, from the collection that has been already described (Lievre ; Laurent-Puig ). Characteristics of the patients are shown in Table 1. All patients were refractory to FOLFOX and/or FOLFIRI regimen and were treated with anti-EGFR antibodies. In summary, 95% of the patients were treated by a combination of anti-EGFR and irinotecan, 2% were treated by a monotherapy using cetuximab or panitumumab, 1% were treated by a combination of anti-EGFR with a chemotherapy associating 5-fluorouracil and oxaliplatin and 1% were treated by 5-fluorouracil. The outcomes were measured on RECIST, PFS and OS. Evaluation was measured every 2 to 3 months according to French standard of care. Patients with complete and partial responses were considered as responders and patients with stable disease or progression as non-responders. Tumours were collected at the time of surgery. All samples were available either as frozen samples or as formalin-fixed paraffin-embedded tissues. This study has been approved by the Ile-de-France ethical committee number 2 (ID-RCB-AO1058-47/2008-135).

Table 1

Characteristics of patients analysed for EGFR pathway protein using the Bioplex protein assay

Median age (years; range)	61 (22–80)
Gender
Male	58
Female	41
NA	1
Number of chemotherapy line used before anti-EGFR
0	0%
1	5%
2	46%
3	36%
4	8%
5	3%
6	2%
Death
Yes	83%
No	17%
Median OS (month–range)	10.0 (0.1–48.1)
Median PFS (month–range)	4.1 (0–20.0)
Response
Complete response	2%
Partial response	26%
Stable disease	34%
Progressive disease	37%
NA	1%
WHO
0	28%
1	48%
2	11%
3	1%
NA	12%
Mutation
KRAS	36%
NRAS	4%
BRAF	2%

Abbreviations: EGFR=epidermal growth factor receptor; PFS=progression-free survival; OS=overall survival; WHO=World Health Organization.

Protein extraction

Whole-cell protein extraction was performed from 5-mm frozen tumour tissues using the Kit RIPA lysis Buffer 1X (Tebubio, Le Perray en Yvelines, France) prepared with inhibitors (PMSF, protease inhibitor cocktail and sodium orthonavate; Dutscher, Brumath, France) according to the manufacturer's recommendations.

Phosphoprotein analysis

The expression of phosphorylated-EGFR (pEGFR(Tyr)) and the key phosphorylated proteins of EGFR downstream signalling pathway such as pAKT(Ser473), pERK1-2(Thr202/Tyr204, Thr185/Tyr187), pGSK3(Ser21/Ser9), pIGF1R(Tyr1131), pMEK1(Ser217/Ser221), pP70S6K(Thr421/Ser424) and pP90RSK(Thr359/Ser363) were quantitatively analysed in freshly frozen samples using BioPlex Protein assay (BPA; Bio-Rad, Marnes-la-coquette, France) as previously described and validated (Chergui ). This technique is on the basis of multiplex sandwich bead immunoassays. Protein extracts were transferred into 96-well dishes and diluted with buffered solution. Fluorescent capture beads coupled to antibodies directed against the phosphoproteins pEGFR(Tyr), pAKT(Ser473), pERK1-2(Thr202/Tyr204, Thr185/Tyr187), pGSK3(Ser21/Ser9), pIGF1R(Tyr1131), pMEK1(Ser217/Ser221), pP70S6K(Thr421/Ser424) and pP90RSK(Thr359/Ser363) were mixed and added into each well and incubated overnight. Following incubation, the plates were washed and incubated with biotinylated antibodies fixing each target protein. Streptavidin–phycoerythrin solution was then added. The analysis consisted of a double laser fluorescence detection allowing simultaneous identification of the target protein through the red fluorescence emission signal of the bead and quantification of the target protein through the fluorescence intensity of phycoerythrin. Results were recorded as the mean fluorescence intensities and compared with negative controls. Positive controls consisting of standard protein extracts from cell lines were added to each series. All results were normalised through the different batches of analyses by the same mutated tumour sample. The expression level of each phosphoprotein was recorded as fluorescence arbitrary unit. BPA assay requires 15–20 mg of tissue containing more than 50% of tumour material or equivalent amount of protein extract (that is, 25 mg total protein per assay in triplicate). All samples were assessed in triplicate to ensure reproducibility of the data.

Determination of KRAS, NRAS and BRAF mutations

KRAS exon 2 and 3 and BRAF exon 15 mutation statuses were assessed by allelic discrimination using TaqMan qPCR probes as previously described (Lievre ; Laurent-Puig ). KRAS exon 4 and NRAS exon 2, 3 and 4 mutation statuses were assessed using Sanger sequencing and next-generation sequencing. Patients with wild-type tumours for exons 2, 3 and 4 of KRAS and NRAS were considered as «RAS wild-type».

Statistical analysis

Statistical analysis was performed using the SAS software (SAS Institute, Cary, NC, USA 25513; version 9.2). A P-value<0.05 was considered as statistically significant.

Quantitative variables were described by the mean, median and interquartile range, and qualitative variables by frequency and percentage. The normality of the distribution of the expression of each phosphoprotein was assessed with the Shaprio–Wilk test.

Overall survival and PFS were described using the Kaplan–Meier method and compared between RAS mutation status with bivariate Cox proportional hazard model. Comparisons of each phosphoprotein according to the RAS mutation status were performed by a Mann–Whitney U-test.

In the subgroup of RAS wild-type patients, the discriminant power of each phosphoprotein to predict the response rate was estimated with the area under the receiver operating characteristic (AUC) and their 95% confidence intervals. An AUC equal to 0.5 corresponds to no discrimination; an AUC between 0.5 and 0.7 means a low discriminant power; an AUC between 0.7 and 0.8 corresponds to an acceptable discrimination and above 0.8 to an excellent discrimination. Each phosphoprotein was then dichotomised according to the best cutoff maximising the sensitivity and the specificity. These thresholds were computed owing to the Youden Index.

The percentage of response rate was calculated for each level of the new binary variables and compared with a χ2-test or a Fisher Exact test. The odds ratio and their 95% confidence interval were computed with a bivariate logistic regression.

The phosphoproteins with a P-value less than 0.2 in bivariate analysis were introduced in a multivariate logistic regression with backward selection (Harrell ). The stability of the selected model was investigated using the bootstrap-resampling method (Sauerbrei and Schumacher, 1992). The discriminant power of the final model was estimated using the AUC with 95% CI. We addressed the optimism bias using a bootstrap-resampling method. The stability of the score was estimated by cross-validation (Altman and Royston, 2000).

Results

Among the 100 tumour samples, 40 were found with a RAS mutation (36% with KRAS mutation and 4% with NRAS mutation) and 60% were RAS wild type. Among the RAS wild-type patients, 45.0% achieved a complete or partial response, and 55.0% had a stable disease or progression (P<0.001) when treated with anti-EGFR. Patients with a RAS mutation had significant lower PFS (HR=3.04 (1.91; 4.83)) and OS (HR=2.49 (1.56; 3.97)). Progression-free survival and OS were significantly higher in RAS wild-type patients (Figure 1). Expression levels of pAKT, pERK1/2 and pMEK1 were significantly lower in RAS wild-type patients than in RAS-mutated patients (P=0.0246; P=0.004; P=0.0110, respectively), and no significant difference was observed between RAS wild-type and RAS-mutated tumours in the expression of pEGFR, pGSK3, pIGFR, pP70S6K and pP90SRK (Table 2).

Figure 1

Kaplan–Meier's curve of

Table 2

Expression of phosphoproteins in all samples, RAS-mutated samples and RAS wild-type samples

	All*	RAS mutated*	RAS wild type*	P
pAKT	104.01; 98.5 (86.04; 124.55)	110.79; 109.46 (93.31; 126.27)	99.49; 92.83 (82.10; 122.07)	0.0246
pEGFR	139.35; 109.39 (86.69; 141.04)	174.69; 105.38 (88.58; 146.54)	115.80; 110.11 (85.16; 139.13)	0.9355
pERK1/2	245.73; 138.98 (101.59; 263.85)	350.22; 178.23 (137.87; 395.01)	176.07; 113.83 (92.04; 169.58)	0.0004
pGSK3	133.03; 101.05 (66.38; 183.82)	135.62; 97.79 (66.38; 168.41)	131.30; 105.88 (65.32; 184.71)	0.7676
pIGFR	83.65; 74.83 (64.29; 96.68)	83.52; 70.81 (63.89; 91.62)	83.74; 77.13 (64.70; 101.09)	0.4024
pMEK1	272.19; 131.82 (80.42; 280.74)	411.25; 209.52 (87.89; 402.14)	179.48; 104.09 (77.70; 201.28)	0.0110
pP70S6K	123.31; 112.49 (89.14; 135.21)	131.39; 124.47 (93.58; 136.83)	117.93; 105.24 (85.00; 135.21)	0.1690
pP90SRK	74.37; 75.61 (49.71; 95.38)	81.03; 80.27 (48.73; 101.78)	69.93; 71.78 (50.11; 88.26)	0.1303

Expression of pAKT, pERK1/2 and pMEK1 is significantly lower in RAS wild-type patients.

*Mean; median (lower quartile; upper quartile).

In RAS wild-type patients, the AUC of the expression of each phosphoprotein according to the response rate were under 0.80. After dichotomising the expression of each phosphoprotein with the best threshold to maximise the sensitivity and the specificity, the response rate was significantly higher for tumours that express pEGFR and pAKT above the calculated threshold (P=0.0258 and P=0.0277, respectively). No significant relation was found between the response rate and the level of expression of the other phosphoproteins (Table 3a).

Table 3

Expression of phoshoproteins and response rate

(a)
pProt	AUC and 95% CI	Thresholda	CR or PR	P-value	OR and 95% CI
pAKT	0.588 (0.436; 0.74)	<120	36.36% (16)	0.0258	1
		⩾120	68.75% (11)		3.85 (1.13; 13.07)
pEGFR	0.634 (0.491; 0.776)	<100	25.00% (5)	0.0277	1
		⩾100	55.00% (22)		3.67 (1.12; 12.03)
pERK1/2	0.497 (0.348; 0.646)	<110	38.45% (10)	0.373	1
		⩾110	50.00% (17)		1.60 (0.57; 4.51)
pGSK3	0.504 (0.355; 0.653)	<130	38.24% (13)	0.228	1
		⩾130	53.85% (14)		1.88 (0.67; 5.31)
pIGFR	0.513 (0.365; 0.662)	<60	33.33% (3)	0.445	1
		⩾60	47.06% (24)		1.78 (0.40; 7.89)
pMEK1	0.497 (0.348; 0.646)	<100	38.46% (10)	0.373	1
		⩾100	50.00% (17)		1.60 (0.57; 4.51)
pP70S6K	0.585 (0.439; 0.731)	<105	36.67% (11)	0.194	1
		⩾105	53.33% (16)		1.97 (0.70; 5.54)
pP90SRK	0.577 (0.43; 0.723)	<70	35.71% (10)	0.176	1
		⩾70	53.13% (17)		2.04 (0.72; 5.77)
(b)
Multivariate
pProt	OR and 95% CI
pAKT	1
	3.70 (1.04;13.18)
pEGFR	1
	3.53 (1.03;12.11)

Abbreviations: AUC=area under the receiver operating characteristic; CI=confidence interval; CR=complete response; OR=odds ratio; PR=partial response.

(a) Patients with expression of pAKT over 120 and pEGFR over 100 achieved significantly better response rate. (b) These two phosphoproteins remained statistically significant in multivariate analysis to predict response after bootstrap validation. Statistically significant p-values are presented in bold.

Optimal threshold according to the Youden index.

These two phosphoproteins remained statistically significant in multivariate analysis to predict response after bootstrap validation (Table 3b). The AUC of the final model was 0.673 (95% CI (0.540; 0.805)). After correction for the optimism bias, the AUC was 0.654. The cross-validation covariate was associated with a (beta) value of 1.08 close to 1.

No statistically significant difference was found between expression of phosphoproteins and PFS or OS in RAS wild-type tumour-bearing patients (data not shown).

Discussion

The development of molecular biology and knowledge about signalling pathways in cancer have been a real opportunity to improve treatment for patients with mCRC. The history of RAS and the use of anti-EGFR therapies are characteristic of the importance of identifying response-predictive biomarkers and molecular testing. Unluckily, a large proportion of RAS wild-type patients do not reach benefit from anti-EGFR-targeted therapies and it is important to explore new targets and identify patients who will benefit from targeted therapies. The selection of these patients is crucial to avoid unworthy toxicities or inefficient treatment, but also for economic issues.

In our study, the repartition of RAS mutation in the studied population is consistent with previously published data (Sorich ). The observed response rate is also consistent with known data with few or no response in RAS-mutant tumour (Bokemeyer ). Complete or partial responses were achieved in ∼45.0% of RAS wild-type tumour-bearing patients, which is higher than that described in previously published data for patients treated for second or further line with an association for anti-EGFR and FOLFIRI (Iwamoto ). This difference can probably be related to the lower median age of our population and the better WHO status.

The level of expression of phosphoproteins in the RAS-mutated patients was found to be significantly higher than in RAS wild-type patients as we expected. The accumulation of RAS protein under its activated form causes the activation of RAF/RAS/MAPK and PI3K/AKT/mTOR signalling pathways, resulting in an overexpression of the downstream proteins involved in these pathways (Newbold, 1984).

In our study, the level of expression of phosphoproteins was measured with the BPA and not with immunohistochemistry. BPA is a validated method and is comparable to western blot analysis (Chergui ). Most of the studies published in the literature are using western blot or immunohistochemistry, which are qualitative and semiquantitative assays, respectively. BPA is a quantitative assay and allows more precision in the measurement of the level of phosphorylated proteins. However, immunohistochemistry allows the determination of the expression according to the localisation of the target—that is, nucleus, cytoplasm, stroma and so on—whereas BPA is using the whole tumour and is unable to determine the localisation of the overexpression in the cell. This difference of technical approach may be responsible to the discrepancies we will discuss further in this paragraph.

Overexpression of pEGFR and pAKT was significantly associated with higher response rate in RAS wild-type patients; however, no statistically significant relation was found between the expression of phosphoproteins and PFS and OS. These date do not confirm our previously published data (Perkins ). In the study published by Perkins ), PFS was related to a lower expression of pMEK and pP70S6K in wild-type KRAS exon 2 tumour-bearing patients. Among the 42 patients included in this study, only 23 were wild-type KRAS tumour-bearing. The lack of subjects in the wild-type KRAS class and the impact of KRAS exon 3, 4 and NRAS exon 2, 3 and 4 mutations may explain this difference with our present data. Moreover, although being in the same order as sample sizes in other papers investigating phosphorylated proteins as putative response-predictive biomarkers in mCRC (Scartozzi ; Cejas. ), the statistical power would have been better with a larger sample size and these results have to be confirmed in the future on a largest population.

In our paper, phosphoproteins have been extracted and analysed from primary tumours; however, some studies suggest that the correlation of the expression of EGFR and AKT between the primary and the metastasis does not appear in 20–40% of cases (Scartozzi ; Cejas. ). However, in these two papers pEGFR and pMAPK have been investigated using IHC considering weakly stained samples as negative. In our paper, phosphorylated protein expression was quantified in frozen-specimen protein extracts and considered as positive when exceeding 3 × the background noise value, as validated in Chergui . In addition, in these two papers, the analyses were performed in FFPE specimens, which does not warrant optimal preservation of phosphoproteins as the freezing procedure we used. According to these published data, it would be of interest to compare the expression of pEGFR and pAKT in primary and metastasis in a larger series of samples and to compare these data to response rate, PFS and OS.

The relation between the overexpression of pEGFR and response rate has already been described (Razis ) and can be explained by the avidity of the tumour cell for its ligand: the more the cell is dependent on its ligand, the more EGFR will be expressed and the more anti-EGFR will be effective against tumour cells. Considering this ligand avidity, it is not surprising that a better response rate was achieved with anti-EGFR with RAS wild-type tumours overexpressing pEGFR. More surprisingly, a significant relation has been found between higher response rate and overexpression of pAKT. It is known that activation of AKT is dependent on the activation of EGFR (Mirzoeva ). It has also been described that expression of pAKT was associated with lower response rate and lower PFS (Scartozzi ); however, the studied population was selected between 2007 and 2011, and no indication mentions whether the patients treated by anti-EGFR were only RAS wild-type or not. In our study, the expression of pAKT was associated with a higher response rate, but no statistically significant difference in PFS and OS was found (HR=1.004; 95% CI (0.995; 1.013) and HR=1.006; 95% CI (0.996; 1.016), respectively). In the study of Saglam , an increasing in the pAKT expression has been described during colorectal cancer progression. We hypothesise that the overexpression of pEGFR induces better response rate for RAS wild-type patients treated with anti-EGFR; however, this overexpression of pEGFR induces an overexpression of pAKT that leads to progression and annihilates anti-EGFR benefits and explains that no statistically significant difference in PFS was found. As already suggested, AKT may be an interesting target in colorectal cancers (De Roock ; Jeong ). Ongoing trials are currently testing an allosteric inhibitor of AKT, the MK-2206 (Agarwal ; Jeong ) in patients with mCRC. On the basis of our results, we hypothesise that associating anti-EGFR and anti-AKT therapies may be of interest in RAS wild-type patients.

In conclusion, these data suggest that the overexpression of pEGFR and pAKT may be response-predictive biomarkers for RAS wild-type patients treated with anti-EGFR, and, on the basis of our results, we hypothesise that the association of anti-EGFR and anti-AKT therapies could be of interest.