Skip to main content

Beyond KRAS mutation status: influence of KRAScopy number status and microRNAs on clinical outcome to cetuximab in metastatic colorectal cancer patients

Abstract

Background

KRAS mutation is a negative predictive factor for treatment with anti-epidermal growth factor receptor (EGFR) antibodies in metastatic colorectal cancer (mCRC). Novel predictive markers are required to further improve the selection of patients for this treatment. We assessed the influence of modification of KRAS by gene copy number aberration (CNA) and microRNAs (miRNAs) in correlation to clinical outcome in mCRC patients treated with cetuximab in combination with chemotherapy and bevacizumab.

Methods

Formalin-fixed paraffin-embedded primary tumour tissue was used from 34 mCRC patients in a phase III trial, who were selected based upon their good (n = 17) or poor (n = 17) progression-free survival (PFS) upon treatment with cetuximab in combination with capecitabine, oxaliplatin, and bevacizumab. Gene copy number at the KRAS locus was assessed using high resolution genome-wide array CGH and the expression levels of 17 miRNAs targeting KRAS were determined by real-time PCR.

Results

Copy number loss of the KRAS locus was observed in the tumour of 5 patients who were all good responders including patients with a KRAS mutation. Copy number gains in two wild-type KRAS tumours were associated with a poor PFS. In KRAS mutated tumours increased miR-200b and decreased miR-143 expression were associated with a good PFS. In wild-type KRAS patients, miRNA expression did not correlate with PFS in a multivariate model.

Conclusions

Our results indicate that the assessment of KRAS CNA and miRNAs targeting KRAS might further optimize the selection of mCRC eligible for anti-EGFR therapy.

Peer Review reports

Background

Recent advances in our understanding of the specific signalling pathways of cancer cells have introduced targeted therapy into treatment regimes for patients with metastatic colorectal cancer (mCRC). Antibodies against the epidermal growth factor receptor (EGFR), cetuximab and panitumumab, have shown a survival benefit in mCRC patients with KRAS wild-type tumours both as monotherapy [1, 2] and when added to chemotherapy[3, 4]. Patients with a tumour harbouring a KRAS codon 12 or 13 mutation are resistant to anti-EGFR therapy[1, 5]. Therefore the use of these antibodies is restricted to patients with KRAS wild-type tumours. However, within this subset not all patients respond to this treatment, and therefore additional predictive markers are needed. We have previously excluded a discordance in KRAS mutation status between the primary tumour and corresponding metastases as an explanation for the heterogeneous response rate in patients with KRAS wild-type tumours [6]. In routine practice, KRAS mutations in codons 12 and 13 are tested, which comprise approximately 96% of the observed KRAS mutations [7]. Recent data suggest that a codon 13 KRAS mutation has a distinct clinical behaviour and is not associated with cetuximab resistance [8]. Whether other KRAS mutations (like codon 61) result in similar resistance to EGFR monoclonal antibodies remains speculative [9]. A mutation in the BRAF oncogene occurs in approximately 10% of mCRC patients and is restricted to KRAS wild-type tumours, and was first shown to have a negative predictive value for anti-EGFR therapy [10]. Subsequently, we have shown that a BRAF mutation predominantly has a strong negative prognostic value [11]. Other biomarkers in the PI3K and RAS/MAPK pathways [1216], ligands to the EGFR [17, 18], and germline single nucleotide polymorphisms [1921] have not yet shown a predictive value that can be used in clinical practice.

Point mutations in the KRAS oncogene lead to a significantly increased RAS-GTPase activity, ultimately resulting in the stimulation of cell proliferation and the inhibition of apoptosis via the RAS/MAPK pathway [22]. However, in addition to oncogenic mutations, copy number changes of the KRAS gene or posttranslational factors may also be involved in the regulatory mechanism of RAS-GTPase activity. Copy number aberrations (CNA) occur throughout the tumour genome and are an important mechanism in colorectal cancer development[23]. Genome-wide studies in mCRC patients have identified loci that are associated with a poor prognosis [24] and with the prediction of response to chemotherapy[25]. However, little is known about the prevalence and effect of CNA of the KRAS locus on chromosome 12p12.1. By using a TaqMan-based KRAS copy number assay, KRAS amplifications were observed in approximately 2% of the 106 investigated colorectal primary tumours [26]. In CRC cell lines, gains of the KRAS locus were shown to be associated with an eleven-fold increase in RAS-GTPase activity, which is comparable with the twelve-fold increase caused by a codon 12 or 13 mutation [27]. However, no data have been reported on CNA affecting the KRAS locus and their possible association with response to cetuximab.

In recent years, a rapidly expanding interest has manifested on microRNAs (miRNAs). These single stranded RNAs of 19–23 nucleotides regulate gene expression by translational inhibition or mRNA degradation via imperfect base pairing to the 3’-untranslated region (3’UTR) of their target mRNAs [28]. MiRNAs are involved in the development of human cancer, and in case of dysregulation they can act either as oncogenes or tumour suppressors, depending on their target genes [29]. Recently several miRNAs were identified that target KRAS, resulting in the suppression of cancer development [3033]. KRAS contains multiple let-7 complementary sites, allowing the let-7 family of miRNAs to act as a tumour suppressor by regulating the KRAS mRNA [30, 31]. Also mir-18a and miR-143 directly recognize KRAS, and downregulation of these miRNAs accelerates tumourigenesis by reversal of KRAS suppression [32]. The targeting effect of miR-18a on KRAS has been demonstrated in colon cancer cell lines irrespective of KRAS mutation status [33]. MiRNAs interfering with the RAS-signaling pathway may have predictive value or may even serve as targets for treatment. Currently no data are available on the clinical relevance of miRNAs involved in KRAS activity in patients treated with cetuximab.

In this study we analyzed the KRAS copy number status and the expression of miRNAs targeting KRAS in relation to clinical outcome in mCRC patients treated with first-line cetuximab-containing therapy.

Methods

Patients

The patients included in this study participated in the CAIRO2 trial (CKTO 2005–02; ClinTrials.gov NCT00208546) of the Dutch Colorectal Cancer Group (DCCG) [34]. In this multicenter phase III trial, 755 mCRC patients were randomized between first-line treatment with capecitabine (1000 mg/m2 bid.), oxaliplatin (130 mg/m2), and bevacizumab (7.5 mg/kg), or the same schedule with the addition of weekly cetuximab (250 mg/m2, after initial 400 mg/m2). Translational research on tumour tissue was part of the informed consent procedure. The primary end point of the study was progression free survival (PFS), and secondary end points were overall survival, response rate, and toxicity. The median PFS in patients treated with cetuximab was 9.4 months (95% CI 8.4-10.5 months), which was significantly shorter than the PFS of patients treated in the group without cetuximab (median PFS 10.7 months, 95% CI 9.7-12.3 months, p = 0.01). Patients in the cetuximab-group with a KRAS mutated tumour had a significantly decreased median PFS compared to patients with a KRAS wild-type tumour (8.1 versus 10.5 months, respectively, p = 0.04).

For the current analysis we selected patients who had been randomized to the cetuximab treatment arm, received at least three treatment cycles, did not discontinue treatment for other causes than disease progression, had a normal serum lactate dehydrogenase at randomisation, and of whom formalin-fixed paraffin-embedded (FFPE) material of the primary tumour as well as normal tissue was available. Patients with a rectal carcinoma and patients who had received preoperative radiotherapy on the pelvis have been excluded for these analyses. From this group the 17 best and 17 worst responding patients were selected based on both extremes of the PFS time. Throughout the article the terms good and poor responders are used, which does not apply to response according to RECIST, but to the patients with the longest and shortest PFS on cetuximab-based treatment. This outcome parameter was chosen for the current study because it is the best reflection of the clinical trial upon which this analysis is based. Next, especially with respect to targeted agents PFS appears to be superior of response rate in terms of clinical outcome.

DNA extraction and mutation analysis

Genomic DNA was extracted from 4–8 manually micro dissected 50 μm sections of FFPE tissue as previously described [35]. DNA concentration was determined using the Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies Inc., Wilmington, USA). DNA quality was assessed by performing a multiplex PCR using 4 primer sets, resulting in fragments of 100, 200, 300 and 400 base pairs [36]. The KRAS mutation status [35] and BRAF mutation status [11, 12] were assessed by sequencing analysis as previously described.

Assessment of the KRASgene copy number and data analysis

High-resolution genome-wide DNA copy number profiles were generated by array-based comparative genomic hybridization (array CGH) using 720 k Whole-Genome Tiling CGH arrays (Roche NimbleGen Inc., Madison, USA). Optimal signal-to-noise ratios were obtained by hybridizing test (tumour) and reference (normal colon) DNA of similar quality, which was determined by giving similar yield in a Bioscore Screening and Amplification kit (ENZO diagnostics Inc., Farmingdale, USA). For hybridization, 500 ng of amplified DNA from test and reference samples were labelled with Cy3 and Cy5, respectively, using random-primed labelling (Bioprime genomic DNA labelling kit, Invitrogen, Breda, the Netherlands), and hybridized for 48 hours at 42°C using a MAUI hybridization system (Biomicro Systems, Salt Lake City, USA). After washing, arrays were scanned in an Axon Genepix 4200AL microarray scanner. The NimbleScan 2.4 software package (NimbleGen Systems Inc., Madison, USA) was used to calculate log2 ratios after performing spatial correction, normalization and a 25 kb average smoothing window on the data. Further data interpretation and CNA calling was done with Nexus Copy Number 5.0 software (Biodiscovery, El Segundo, USA) using the Rank Segmentation Algorithm. In 26 patients, hybridizations were performed against normal DNA from the same patient to normalize for germline copy number changes. In the other 8 cases, germline copy number changes were excluded using both public (http://projects.tcag.ca/variation/) and private CNA databases. The cut-off value for gene copy number gain and loss were manually set for each sample to adjust for differences in signal strength and incorrectly centered baselines.

Prevalence of KRASlocus gene copy number changes

In order to assess the clinical relevance in terms of prevalence of KRAS gene copy number changes in mCRC patients, we assessed the KRAS gene copy number status in FFPE primary tumour tissue of 225 unselected mCRC patients who participated to our previous phase III study which did not involve the use of targeted agents [37]. In these patients a 250 k oligonucleotide array CGH was performed as previously described[25].

Multiplex Ligation-dependent Probe Amplification (MLPA)

MLPA was performed according to the manufacturer’s instruction using the SALSA P145-A2 kit (MRC Holland, Amsterdam, the Netherlands), containing 40 probes including the 12p12.1 KRAS probe. Briefly, 200 ng DNA was denatured and allowed to hybridize for 16 h at 60°C in a thermocycler. Then SALSA Ligase-65 enzyme was added and ligation was allowed at 54°C. After heat inactivation of the ligase enzyme at 98°C, primers, dNTPs and polymerase were added and PCR amplification was performed for 35 cycles (60s at 95°C, 30s at 60°C, and 90s at 72°C). Reactions were performed on a PTC 200 thermal cycler (MJ Research Inc., Waltham, Massachusetts, USA). One microliter of PCR product was analysed by capillary electrophoresis on an ABI 3730 Analyzer (Applied Biosystems), and quantitative data were obtained by Genemapper analysis (Applied Biosystems).

MLPA data analysis

For each tumour sample, the peak area of the 12p12.1 and reference probes were determined in duplicate for further analysis. The reference peak area was obtained from blood samples from three different individuals, each of which were analysed at least two times independently. In every sample, for every probe, a tumour to normal DNA copy number ratio was calculated by dividing the median area under the peak for the 12p12.1 probe by the value for the reference DNA. Subsequently, all ratios were normalized by setting the median tumour to normal DNA copy number ratio of the reference genes in de probe mixture to 1.0. A ratio lower than 0.8 was considered a loss and a ratio higher than 1.2 a gain.

MiRNA selection

Selection of miRNAs regulating KRAS was performed using PicTar (http://pictar.mdc-berlin.de/), TargetScan (http://www.targetscan.org/) and miRNA targets (http://cbio.mskcc.org/mirnaviewer/). Venn diagram analysis was used to select 14 miRNAs who were identified by at least two algorithms (Additional file 1: Figure S1). In addition to the prediction programs we also selected six extra miRNAs (Let-7, miR-18a, miR-21, miR-133a, miR-133b, miR-205) which have been shown to target KRAS in previous studies [30, 33, 38], resulting in a total test series of 20 miRNAs. Two Taqman microRNA assay were not available (mir-18a, mir-200c), resulting in 18 miRNAs to analyze.

Total RNA extraction, miRNA reverse transcription and real-time PCR

Total RNA was isolated from FFPE tissue of 34 primary tumours and matched normal tissue using the RecoverAll™ Total Nucleic Acid Isolation Kit (Applied Biosystems, Foster city, USA). In brief, four tissue slices of 20 μm were micro dissected and incubated with 100% xylene at 50°C to remove paraffin excess, followed by ethanol washes. Proteins were degraded by protease at 50° and 80°C. The RNA was extracted followed by nuclease digestion. Total RNA quantity and quality were determined using the Nanodrop 26 ND-1000 spectrophotometer (Nanodrop Technologies Inc., Wilmington, USA).

To determine the expression levels of miRNAs, Taqman microRNA assays directed to seventeen miRNAs and the endogenous reference gene (RNU 6B) were used following the manufacture’s protocol (Applied Biosystems, Foster City, USA). Firstly, cDNA was synthesized in duplicate from total RNA using miRNA specific stem loop primers. Reverse transcriptase reactions were conducted using 10 ng total RNA, 1 mM dNTPs, 50 U MultiScribeTM Reverse Transcriptase, 1 x RT buffer, 3.8 U RNase inhibitor and 1 x Taqman® MicroRNA RT Primer (Applied Biosystems, Foster city, USA). The 15 μl reactions were incubated at 16°C for 30 minutes at 42°C for 30 minutes and at 85°C for 5 minutes.

Secondly, the quantitative PCR was performed in which the total mixture of 20 μl included 1.33 μl RT product (1:5 diluted from RT reaction), 1 x Taqman® Universal PCR Master Mix (No AmpErase® UNG, Applied Biosystems, Foster City, USA) and 1 x the dedicated primer and probe mix. The reactions were incubated in a 96-well optical plate at 95°C for 10 minutes, followed by 40 cycles at 95°C for 15 seconds and at 60°C for 1 minute. All reactions were carried out in duplicate in a 7500 Real Time PCR System (Applied Biosystems, Foster City, USA). The threshold cycle (Ct) was defined as the fractional cycle number at which the fluorescence passes the fixed threshold. Relative quantification of miRNA expression was calculated using the ΔΔCt method as described previously [39].

Statistical analysis

PFS was defined as the interval from the date of randomization to the date of first documented disease progression or death, whichever occurred first. Statistical differences of clinical and pathological parameters between good and poor responders were evaluated using the Student’s t-test, Pearson’s χ2 test or Fisher’s exact test where appropriate. The miRNA expression in colorectal tumours was described by the relative quantity (RQ) of the target miRNA, normalized in respect to RNU6B and relative to matched normal tissue. Box plots were used to appreciate the descriptive statistics of the data. Differences in expression of the target miRNA between good and poor responders were evaluated on the log scale (ΔΔCt scale) to obtain normally distributed data. The Student’s t-tests was used in exploratory analyses on the miRNA expression in relation to response and to KRAS mutation status. When focusing on the actually observed PFS, we investigated by Cox regression analysis the influence of each miRNA on PFS, using KRAS mutation status, the interaction term between miRNA and KRAS mutation, and differentiation grade as covariates. Due to the limited number of patients and the ensuing risk for overfitting, it was not possible to assess the influence of all miRNAs together (i.e. correct the influence of miRNA for each other), nor to correct for other baseline characteristics.

Results

Patients

Of the 34 patients selected for this analysis, the median PFS was 22.5 months (range 14.8-39.8 months) in the 17 good responders, and 6.0 months (range 2.3-7.2 months) in the 17 poor responders. Clinical and pathological characteristics of the primary tumour were well balanced between the 17 good and 17 poor responders. Only poor differentiation grade of the primary tumour was more frequently observed in the poor responders.

A KRAS mutation was demonstrated in the primary tumour of 15 patients (6 good responders and 9 poor responders), and KRAS wild-type in the primary tumour of 19 patients (11 good responders and 8 poor responders). KRAS codon 12 mutation was observed in 14 patients, and one poor responder had a codon 13 mutation. Of the KRAS wild-type patients, 4 had a BRAF mutated tumour (1 good responder and 3 poor responders) (Table 1).

Table 1 Clinical and histopathological characteristics of patients and their respective tumours

12p12.1 copy number changes in good and poor responders

By using high resolution array CGH, two copy number gains (of which one amplification) and 5 losses were detected at the 12p12.1 locus where KRAS is localized (Additional file 2: Figure S2). Both copy number gains, which were confirmed by MLPA, were observed in poor responders with a KRAS wild-type tumour. Of these tumours one sample contained a gain of the complete p-arm of chromosome 12 and the other sample contained a high copy number gain of a region including the KRAS locus.

A 12p12.1 copy number loss, detected by array CGH, was observed in the tumour of 5 patients with a good response. One tumour contained a loss of the whole chromosome, three tumours included a loss of the short arm of the chromosome and one tumour contained a loss of a 27.5 Mb region of the short arm of chromosome 12 including the KRAS locus (Additional file 2: Figure S2). Of these 5 tumours with loss of the 12p12.1 locus, 2 tumours harboured a KRAS mutation, and one tumour had a BRAF mutation, suggesting that the mechanism of gene copy number loss is independent of the KRAS and BRAF mutation status (Figure 1).

Figure 1
figure 1

12p12.1 copy number changes in good and poor responders according to the KRAS mutation status. Abbreviations: CNA = copy number aberration. MT = mutant, WT = wild-type.

12p12.1 gene copy number changes in a control group of mCRC patients

In an unselected group of 222 mCRC patients from our previous trial with comparable baseline characteristics[37], the prevalence of 12p12.1 copy number changes was assessed. In this group three amplifications (1.4%), 32 copy number gains (14.4%), and 12 losses (5.4%) of the 12p12.1 locus were observed. There was no effect of KRAS copy number gain or loss on prognosis in these patients treated with first-line chemotherapy without cetuximab (p = 0.97 and p = 0.75, respectively, data not shown).

MiRNA expression in good and poor responders

To assess the role of miRNA expression in relation to clinical outcome, the expression levels of 18 miRNAs targeting KRAS were determined by real-time RT-PCR in 32 primary colorectal tumours relative to their matched normal tissue. Two patients (1 good and 1 poor responder) were not accessible for miRNA expression due to an insufficient RNA amount in normal mucosa. MiR-205 expression was undetectable in both tumour and normal mucosa, therefore 17 miRNAs were included in our final analysis. By using NormFinder [40] and GeNorm [41], the use of RNU6B as a reference gene was justified.

The expression level of 14 miRNAs showed a trend towards a higher expression in patients with a good response compared to patients with a poor response, however this trend was not statistically significant. MiR-143, miR-133a and miR-133b expression was decreased in patients with a good response, of which miR-143 showed a relative expression in good versus poor responders of 0.49 (p = 0.07) (Table 2).

Table 2 MiRNA expression in good versus poor responders

MiRNA expression in good and poor responders according to KRASmutation status

In patients with a wild-type KRAS tumour, the expression level of miR-181a showed a 1.87-fold increase in good responders compared to poor responders (p = 0.04), which was not observed in patients with mutated KRAS tumours (0.91-fold increase, p = 0.69). A higher expression of miRNAs in wild-type KRAS good responders compared with wild-type KRAS poor responders was also observed for MiR-200b (2.48-fold increase, p = 0.01) and miR-21 (1.66-fold increase, p = 0.06).

A difference between the expression of miR-143 in good versus poor responders was more obvious in mutated KRAS tumours. The relative expression level of miR-143 showed a 0.30 fold increase in mutated KRAS good responders versus mutated KRAS poor responders (p = 0.11) (Figure 2).

Figure 2
figure 2

Box plots of the expression levels of miR-181a, miR-200b and miR-143 in mCRC patients according to clinical outcome and KRAS mutation status. Abbreviations: G = good responders, P = poor responders, MT = mutant, WT = wild-type.

Multivariate model of PFS in relation to miRNA expression and KRASmutation status

Each miRNA was analyzed individually together with differentiation grade as a covariate for PFS in wild-type KRAS and mutated KRAS patients treated with first-line cetuximab-containing therapy (Table 3). Differentiation grade was used as a covariate in the Cox regression model because this pathological feature is a well known prognostic factor and differentially distributed between good and poor responders.

Table 3 A multivariate model in which each miRNA was analyzed individually together with differentiation grade as a predictor for PFS in wild-type- KRAS and mutated- KRAS patients treated with chemotherapy, bevacizumab and cetuximab

Elevated expression of mir-200b was associated with a better PFS in patients with a mutated KRAS tumour (HR 0.56 (0.28-1.15); p = 0.10). This trend was not present in patients with a wild-type KRAS tumour. Surprisingly, increased expression of miR-143 resulted in a shorter PFS in patients with a mutated KRAS tumour (HR 1.59 (1.01-2.50); p = 0.04). The hazard ratio for PFS was not influenced by miR-143 expression in wild-type KRAS tumours.

Discussion

We demonstrated that regulation of the KRAS oncogene at several levels might affect clinical outcome in a selected group of cetuximab-treated mCRC patients treated in a phase III trial [34]. Copy number loss of the KRAS locus was restricted to good responders, whereas a copy number gain was associated with a poor PFS in patients with wild-type KRAS tumours. Increased expression of miR-200b that targets KRAS was associated with improved PFS in patients with a mutated KRAS tumour. Surprisingly, decreased miR-143 expression was correlated with improved PFS in these patients.

The predictive strength of KRAS mutation status stresses the importance of RAS-GTPase activity for the response to cetuximab. Therefore, other regulatory mechanisms of RAS-GTPase activity are obvious novel candidate markers. CNA of the KRAS locus occur independently of the KRAS mutation status in a considerable percentage of colorectal tumours (21.2%) as assessed in a large and unselected mCRC population. Previously, it has been shown that KRAS copy number gains are correlated with increased RAS-GTPase activity in colorectal cell lines and with worse clinical outcome in lung adenocarcinomas [27]. Our results suggest that KRAS copy number gains are associated with worse clinical outcome in wild-type KRAS mCRC patients who are treated with a cetuximab-containing first-line regimen. This influence of KRAS copy number gain on prognosis was absent in mCRC patients treated without cetuximab, suggesting a predictive effect on cetuximab response. The correlation between miRNAs targeting KRAS and PFS was absent in wild-type KRAS patients. Inhibition of KRAS translation by miRNAs is probably only relevant when the KRAS expression levels are high, which is not the case in absence of an activating KRAS mutation.

KRAS mutations occur in approximately 38% of mCRC patients [7], and these patients are currently excluded from treatment with anti-EGFR antibodies. However, in our selected good responders, 6 patients (35%) had a KRAS mutated tumour. A recent publication showed that patients with codon 13-mutated tumours might benefit from cetuximab treatment [8]. In the current series none of the good responders had a tumour with a KRAS codon 13 mutation. Our data show that the presence of KRAS copy number loss in two of the mutated KRAS mCRC patients might justify treatment with cetuximab. Decreased expression of KRAS caused by loss of gene copies in correlation with response to cetuximab has not been described earlier. Despite the limitations in sample size and concomitant treatment our results indicate that patients with KRAS copy number loss might benefit from treatment with an anti-EGFR antibody although their tumour is KRAS mutated.

Next, we demonstrated that increased expression of miR-200b was associated with improved PFS in mutated KRAS patients. We hypothesize that reducing KRAS protein levels in the presence of a mutation might improve clinical outcome in patients treated with cetuximab. KRAS is not the only target of the miR-200 family, miR-200b is also capable of reducing ERRFI-1 mRNA and subsequent activation of EGFR [42]. Adam et al. showed that increased expression of miR-200b facilitates optimal EGFR functionality, resulting in an efficient response of bladder cancer cells to cetuximab. To our knowledge, our results are the first data in vivo suggesting that in the presence of a KRAS mutation, an increased miR-200b expression is associated with an improved PFS in cetuximab-treated mCRC patients. Surprisingly a decreased expression of miR-143 was associated with improved PFS in patients with mutated KRAS tumours. MiR-143 is thought to inhibit KRAS translation and thereby to suppress tumour cell growth during tumourigenesis [32]. In an established tumour, modulation of KRAS by miR-143 may be differentially regulated which could possibly explain our findings. However, since miRNAs are capable of repressing over a hundred different mRNAs [28], miR-143 could also target mRNAs which may be relevant in response to capecitabine, oxaliplatin and bevacizumab. Previous studies on biomarkers have shown divergent results which stresses the importance of the results of our current hypothesis-generating study being confirmed in a larger, independent series of mCRC patients preferably treated with cetuximab monotherapy.

The patients used in this study were derived from a clinical trial, and the observed outcome is also influenced by the effect of the other agents used. The relative contribution of cetuximab to this outcome is therefore unclear. The phase III CAIRO2 trial showed that cetuximab plus chemotherapy and bevacizumab resulted in a significantly decreased median PFS compared to treatment with chemotherapy and bevacizumab alone. The explanation of this detrimental outcome is unclear [43], and complicates the interpretation of the current analysis. Excessive toxicity in the cetuximab group does not appear to be cause of these results. Negative interaction between the antibodies or between antibodies and chemotherapy might have influenced the outcome although preclinical observations supporting this hypothesis are not yet available. The interpretation of the current analysis is complicated by the detrimental outcome of the trial. Whether this outcome also affects the PFS in the good responders remains unclear.

In conclusion, the analysis of KRAS CNA and miRNAs targeting KRAS may optimize the selection of mCRC patients eligible for anti-EGFR therapy. Elevated expression of miR-200b, decreased miR-143 level and copy number losses may identify patients with mutated KRAS tumours who benefit from anti-EGFR therapy, whereas copy number gains in wild-type KRAS patients could predict resistance to cetuximab. Our results are relevant for the development of predictive biomarkers for anti-EGFR therapy, and suggest that the clinical effects of KRAS are the result of a complex interaction of several regulatory mechanisms beyond the KRAS point mutation status.

Conclusions

KRAS activity, an important regulator of response to anti-EGFR therapy, can be influenced by genetic and epigenetic regulation. CNA and specific miRNAs may provide important additional information to KRAS mutation status and their use could further improve the selection of mCRC patients for anti-EGFR therapy. The hypothesis-generating nature of our study urges for our results to be confirmed in larger series.

References

  1. Karapetis CS, Khambata-Ford S, Jonker DJ, O'Callaghan CJ, Tu D, Tebbutt NC, Simes RJ, Chalchal H, Shapiro JD, Robitaille S, et al: K-ras mutations and benefit from cetuximab in advanced colorectal cancer. N Engl J Med. 2008, 359: 1757-1765. 10.1056/NEJMoa0804385.

    Article  CAS  PubMed  Google Scholar 

  2. Amado RG, Wolf M, Peeters M, Van Cutsem E, Siena S, Freeman DJ, Juan T, Sikorski R, Suggs S, Radinsky R, et al: Wild-type KRAS is required for panitumumab efficacy in patients with metastatic colorectal cancer. J Clin Oncol. 2008, 26: 1626-1634. 10.1200/JCO.2007.14.7116.

    Article  CAS  PubMed  Google Scholar 

  3. Van Cutsem E, Kohne CH, Hitre E, Zaluski J, Chang Chien CR, Makhson A, D'Haens G, Pinter 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. 10.1056/NEJMoa0805019.

    Article  CAS  PubMed  Google Scholar 

  4. Douillard JY, Siena S, Cassidy J, Tabernero J, Burkes R, Barugel M, Humblet Y, Bodoky G, Cunningham D, Jassem J, et al: Randomized, phase III trial of panitumumab with infusional fluorouracil, leucovorin, and oxaliplatin (FOLFOX4) versus FOLFOX4 alone as first-line treatment in patients with previously untreated metastatic colorectal cancer: the PRIME study. J Clin Oncol. 2010, 28: 4697-4705. 10.1200/JCO.2009.27.4860.

    Article  CAS  PubMed  Google Scholar 

  5. Lievre A, Bachet JB, Boige V, Cayre A, Le Corre D, Buc E, Ychou M, Bouche O, Landi B, Louvet C, et al: KRAS mutations as an independent prognostic factor in patients with advanced colorectal cancer treated with cetuximab. J Clin Oncol. 2008, 26: 374-379. 10.1200/JCO.2007.12.5906.

    Article  CAS  PubMed  Google Scholar 

  6. Knijn N, Mekenkamp LJ, Klomp M, Vink-Borger ME, Tol J, Teerenstra S, Meijer JW, Tebar M, Riemersma S, van Krieken JH, et al: KRAS mutation analysis: a comparison between primary tumours and matched liver metastases in 305 colorectal cancer patients. Br J Cancer. 2011, 104: 1020-1026. 10.1038/bjc.2011.26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Oliveira C, Westra JL, Arango D, Ollikainen M, Domingo E, Ferreira A, Velho S, Niessen R, Lagerstedt K, Alhopuro P, et al: Distinct patterns of KRAS mutations in colorectal carcinomas according to germline mismatch repair defects and hMLH1 methylation status. Hum Mol Genet. 2004, 13: 2303-2311. 10.1093/hmg/ddh238.

    Article  CAS  PubMed  Google Scholar 

  8. De Roock W, Jonker DJ, Di Nicolantonio F, Sartore-Bianchi A, Tu D, Siena S, Lamba S, Arena S, Frattini M, Piessevaux H, et al: Association of KRAS p.G13D mutation with outcome in patients with chemotherapy-refractory metastatic colorectal cancer treated with cetuximab. JAMA. 2010, 304: 1812-1820. 10.1001/jama.2010.1535.

    Article  CAS  PubMed  Google Scholar 

  9. Loupakis F, Ruzzo A, Cremolini C, Vincenzi B, Salvatore L, Santini D, Masi G, Stasi I, Canestrari E, Rulli E, et al: KRAS codon 61, 146 and BRAF mutations predict resistance to cetuximab plus irinotecan in KRAS codon 12 and 13 wild-type metastatic colorectal cancer. Br J Cancer. 2009, 101: 715-721. 10.1038/sj.bjc.6605177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Di Nicolantonio F, Martini M, Molinari F, Sartore-Bianchi A, Arena S, Saletti P, De Dosso S, Mazzucchelli L, Frattini M, Siena S, et al: Wild-type BRAF is required for response to panitumumab or cetuximab in metastatic colorectal cancer. J Clin Oncol. 2008, 26: 5705-5712. 10.1200/JCO.2008.18.0786.

    Article  CAS  PubMed  Google Scholar 

  11. Tol J, Nagtegaal ID, Punt CJ: BRAF mutation in metastatic colorectal cancer. N Engl J Med. 2009, 361: 98-99. 10.1056/NEJMc0904160.

    Article  CAS  PubMed  Google Scholar 

  12. Tol J, Dijkstra JR, Klomp M, Teerenstra S, Dommerholt M, Vink-Borger ME, van Cleef PH, van Krieken JH, Punt CJ, Nagtegaal ID: Markers for EGFR pathway activation as predictor of outcome in metastatic colorectal cancer patients treated with or without cetuximab. Eur J Cancer. 2010, 46: 1997-2009. 10.1016/j.ejca.2010.03.036.

    Article  CAS  PubMed  Google Scholar 

  13. Sartore-Bianchi A, Martini M, Molinari F, Veronese S, Nichelatti M, Artale S, Di Nicolantonio F, Saletti P, De Dosso S, Mazzucchelli L, et al: PIK3CA mutations in colorectal cancer are associated with clinical resistance to EGFR-targeted monoclonal antibodies. Cancer Res. 2009, 69: 1851-1857. 10.1158/0008-5472.CAN-08-2466.

    Article  CAS  PubMed  Google Scholar 

  14. Sartore-Bianchi A, Di Nicolantonio F, Nichelatti M, Molinari F, De Dosso S, Saletti P, Martini M, Cipani T, Marrapese G, Mazzucchelli L, et al: Multi-determinants analysis of molecular alterations for predicting clinical benefit to EGFR-targeted monoclonal antibodies in colorectal cancer. PLoS One. 2009, 4: e7287-10.1371/journal.pone.0007287.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Frattini M, Saletti P, Romagnani E, Martin V, Molinari F, Ghisletta M, Camponovo A, Etienne LL, Cavalli F, Mazzucchelli L: PTEN loss of expression predicts cetuximab efficacy in metastatic colorectal cancer patients. Br J Cancer. 2007, 97: 1139-1145. 10.1038/sj.bjc.6604009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. De Roock W, Claes B, Bernasconi D, De Schutter J, Biesmans B, Fountzilas G, Kalogeras KT, 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. 10.1016/S1470-2045(10)70130-3.

    Article  CAS  PubMed  Google Scholar 

  17. Khambata-Ford S, Garrett CR, Meropol NJ, Basik M, Harbison CT, Wu S, Wong TW, Huang X, Takimoto CH, Godwin AK, et al: Expression of epiregulin and amphiregulin and K-ras mutation status predict disease control in metastatic colorectal cancer patients treated with cetuximab. J Clin Oncol. 2007, 25: 3230-3237. 10.1200/JCO.2006.10.5437.

    Article  CAS  PubMed  Google Scholar 

  18. Jacobs B, De Roock W, Piessevaux H, Van Oirbeek R, Biesmans B, De Schutter J, Fieuws S, Vandesompele J, Peeters M, Van Laethem JL, et al: Amphiregulin and epiregulin mRNA expression in primary tumours predicts outcome in metastatic colorectal cancer treated with cetuximab. J Clin Oncol. 2009, 27: 5068-5074. 10.1200/JCO.2008.21.3744.

    Article  CAS  PubMed  Google Scholar 

  19. Pander J, Gelderblom H, Antonini NF, Tol J, van Krieken JH, van der Straaten T, Punt CJ, Guchelaar HJ: Correlation of FCGR3A and EGFR germline polymorphisms with the efficacy of cetuximab in KRAS wild-type metastatic colorectal cancer. Eur J Cancer. 2010, 46: 1829-1834. 10.1016/j.ejca.2010.03.017.

    Article  CAS  PubMed  Google Scholar 

  20. Bibeau F, Lopez-Crapez E, Di Fiore F, Thezenas S, Ychou M, Blanchard F, Lamy A, Penault-Llorca F, Frebourg T, Michel P, et al: Impact of Fc{gamma}RIIa-Fc{gamma}RIIIa polymorphisms and KRAS mutations on the clinical outcome of patients with metastatic colorectal cancer treated with cetuximab plus irinotecan. J Clin Oncol. 2009, 27: 1122-1129. 10.1200/JCO.2008.18.0463.

    Article  CAS  PubMed  Google Scholar 

  21. Zhang W, Gordon M, Schultheis AM, Yang DY, Nagashima F, Azuma M, Chang HM, Borucka E, Lurje G, Sherrod AE, et al: FCGR2A and FCGR3A polymorphisms associated with clinical outcome of epidermal growth factor receptor expressing metastatic colorectal cancer patients treated with single-agent cetuximab. J Clin Oncol. 2007, 25: 3712-3718. 10.1200/JCO.2006.08.8021.

    Article  CAS  PubMed  Google Scholar 

  22. Scaltriti M, Baselga J: The epidermal growth factor receptor pathway: a model for targeted therapy. Clin Cancer Res. 2006, 12: 5268-5272. 10.1158/1078-0432.CCR-05-1554.

    Article  CAS  PubMed  Google Scholar 

  23. Shih IM, Zhou W, Goodman SN, Lengauer C, Kinzler KW, Vogelstein B: Evidence that genetic instability occurs at an early stage of colorectal tumourigenesis. Cancer Res. 2001, 61: 818-822.

    CAS  PubMed  Google Scholar 

  24. Poulogiannis G, Ichimura K, Hamoudi RA, Luo F, Leung SY, Yuen ST, Harrison DJ, Wyllie AH, Arends MJ: Prognostic relevance of DNA copy number changes in colorectal cancer. J Pathol. 2010, 220: 338-347. 10.1002/path.2640.

    Article  CAS  PubMed  Google Scholar 

  25. Postma C, Koopman M, Buffart TE, Eijk PP, Carvalho B, Peters GJ, Ylstra B, van Krieken JH, Punt CJ, Meijer GA: DNA copy number profiles of primary tumours as predictors of response to chemotherapy in advanced colorectal cancer. Ann Oncol. 2009, 20: 1048-1056. 10.1093/annonc/mdn738.

    Article  CAS  PubMed  Google Scholar 

  26. Smith G, Bounds R, Wolf H, Steele RJ, Carey FA, Wolf CR: Activating K-Ras mutations outwith 'hotspot' codons in sporadic colorectal tumours - implications for personalised cancer medicine. Br J Cancer. 2010, 102: 693-703. 10.1038/sj.bjc.6605534.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Soh J, Okumura N, Lockwood WW, Yamamoto H, Shigematsu H, Zhang W, Chari R, Shames DS, Tang X, MacAulay C, et al: Oncogene mutations, copy number gains and mutant allele specific imbalance (MASI) frequently occur together in tumour cells. PLoS One. 2009, 4: e7464-10.1371/journal.pone.0007464.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Slaby O, Svoboda M, Michalek J, Vyzula R: MicroRNAs in colorectal cancer: translation of molecular biology into clinical application. Mol Cancer. 2009, 8: 102-10.1186/1476-4598-8-102.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, et al: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995.

    Article  CAS  PubMed  Google Scholar 

  30. Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014.

    Article  CAS  PubMed  Google Scholar 

  31. Akao Y, Nakagawa Y, Naoe T: let-7 microRNA functions as a potential growth suppressor in human colon cancer cells. Biol Pharm Bull. 2006, 29: 903-906. 10.1248/bpb.29.903.

    Article  CAS  PubMed  Google Scholar 

  32. Chen X, Guo X, Zhang H, Xiang Y, Chen J, Yin Y, Cai X, Wang K, Wang G, Ba Y, et al: Role of miR-143 targeting KRAS in colorectal tumourigenesis. Oncogene. 2009, 28: 1385-1392. 10.1038/onc.2008.474.

    Article  CAS  PubMed  Google Scholar 

  33. Tsang WP, Kwok TT: The miR-18a* microRNA functions as a potential tumour suppressor by targeting on K-Ras. Carcinogenesis. 2009, 30: 953-959. 10.1093/carcin/bgp094.

    Article  CAS  PubMed  Google Scholar 

  34. Tol J, Koopman M, Cats A, Rodenburg CJ, Creemers GJ, Schrama JG, Erdkamp FL, Vos AH, van Groeningen CJ, Sinnige HA, et al: Chemotherapy, bevacizumab, and cetuximab in metastatic colorectal cancer. N Engl J Med. 2009, 360: 563-572. 10.1056/NEJMoa0808268.

    Article  CAS  PubMed  Google Scholar 

  35. Tol J, Dijkstra JR, Vink-Borger ME, Nagtegaal ID, Punt CJ, van Krieken JH, Ligtenberg MJ: High sensitivity of both sequencing and real-time PCR analysis of KRAS mutations in colorectal cancer tissue. J Cell Mol Med. 2010, 14: 2122-2131.

    Article  CAS  PubMed  Google Scholar 

  36. van Dongen JJ, Langerak AW, Bruggemann M, Evans PA, Hummel M, Lavender FL, Delabesse E, Davi F, Schuuring E, Garcia-Sanz R, et al: Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations: report of the BIOMED-2 Concerted Action BMH4-CT98-3936. Leukemia. 2003, 17: 2257-2317. 10.1038/sj.leu.2403202.

    Article  CAS  PubMed  Google Scholar 

  37. Koopman M, Antonini NF, Douma J, Wals J, Honkoop AH, Erdkamp FL, de Jong RS, Rodenburg CJ, Vreugdenhil G, Loosveld OJ, et al: Sequential versus combination chemotherapy with capecitabine, irinotecan, and oxaliplatin in advanced colorectal cancer (CAIRO): a phase III randomised controlled trial. Lancet. 2007, 370: 135-142. 10.1016/S0140-6736(07)61086-1.

    Article  CAS  PubMed  Google Scholar 

  38. Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M, et al: Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoural tissues. Mol Cancer. 2006, 5: 29-

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.

    Article  CAS  PubMed  Google Scholar 

  40. Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.

    Article  CAS  PubMed  Google Scholar 

  41. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034 Epub 2002 June 18

    Google Scholar 

  42. Adam L, Zhong M, Choi W, Qi W, Nicoloso M, Arora A, Calin G, Wang H, Siefker-Radtke A, McConkey D, et al: miR-200 expression regulates epithelial-to-mesenchymal transition in bladder cancer cells and reverses resistance to epidermal growth factor receptor therapy. Clin Cancer Res. 2009, 15: 5060-5072. 10.1158/1078-0432.CCR-08-2245.

    Article  CAS  PubMed  Google Scholar 

  43. Punt CJ, Tol J: More is less – combining targeted therapies in metastatic colorectal cancer. Nat Rev Clin Oncol. 2009, 6: 731-733. 10.1038/nrclinonc.2009.168.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We thank Begona Diosdado for her practical support and expertise in analyzing the miRNA experiments. This study was supported by the Dutch Colorectal Cancer Group and the Netherlands Organization for Health Research and Development (ZonMw). No potential conflicts of interest were disclosed.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Iris D Nagtegaal.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LJMM, JT, JRD, IdK, MEVB, SvV, EK, and EV carried out the experiments. LJMM, JT, JRD, EV, ST, RK and IDN analysed and interpreted data. LJMM and JT were responsible for the study design and wrote the paper. MK, GAM, JHJMvK, RK, CJAP and IDN revised the manuscript critically. All authors were involved in editing the paper and had final approval of the submitted and published versions.

Leonie JM Mekenkamp, Jolien Tol contributed equally to this work.

Electronic supplementary material

Additional file 1: Figure S1. Selection of MiRNAs that were identified by multiple algorithms. (TIFF 2 MB)

12885_2011_3420_MOESM2_ESM.tiff

Additional file 2: Figure S2. Heat map representation and individual array CGH plots of patients with KRAS copy number aberrations. A: Heat map representation of the 7 patients with CNA of the KRAS locus. Each row represents a patient with a CNA of the KRAS locus (loss, gain, amplification). Whole chromosome 12, containing the KRAS locus, is depicted on the horizontal axis. B: Amplification of a genomic region in 12p12.1 detected by array CGH in two patients, as confirmed by MLPA. The DNA log2 ratios and whole chromosome 12 are represented on the vertical and horizontal axis, respectively. C: Array CGH plot of chromosome 12 of a patient with a deletion of KRAS, which could not be validated using MLPA, but was detected by the Nexus copy number algorithm. Most of the genomic deletions detected by array CGH appeared to be present subclonal, below the detection threshold of MLPA. The DNA log2 ratios and whole chromosome 12 are represented on the vertical and horizontal axis, respectively. Abbreviations: ampl = amplification. (TIFF 12705 kb) (TIFF 12 MB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mekenkamp, L.J., Tol, J., Dijkstra, J.R. et al. Beyond KRAS mutation status: influence of KRAScopy number status and microRNAs on clinical outcome to cetuximab in metastatic colorectal cancer patients. BMC Cancer 12, 292 (2012). https://doi.org/10.1186/1471-2407-12-292

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-12-292

Keywords