Abstract
In metastatic colorectal cancer (mCRC), recent studies have shown the importance to accurately quantify low-abundance mutations of the RAS pathway because anti-EGFR therapy may depend on certain mutation thresholds. We aimed to evaluate the added predictive value of an extended RAS panel testing using two commercial assays and a highly sensitive and quantitative digital PCR (dPCR). Tumor samples from 583 mCRC patients treated with anti–EGFR- (n = 255) or bevacizumab- (n = 328) based therapies from several clinical trials and retrospective series from the TTD/RTICC Spanish network were analyzed by cobas, therascreen, and dPCR. We evaluated concordance between techniques using the Cohen kappa index. Response rate, progression-free survival (PFS), and overall survival (OS) were correlated to the mutational status and the mutant allele fraction (MAF). Concordance between techniques was high when analyzing RAS and BRAF (Cohen kappa index around 0.75). We observed an inverse correlation between MAF and response in the anti-EGFR cohort (P < 0.001). Likelihood ratio analysis showed that a fraction of 1% or higher of any mutated alleles offered the best predictive value. PFS and OS were significantly longer in RAS/BRAF wild-type patients, independently of the technique. However, the predictability of both PFS and OS were higher when we considered a threshold of 1% in the RAS scenario (HR = 1.53; CI 95%, 1.12–2.09 for PFS, and HR = 1.9; CI 95%, 1.33–2.72 for OS). Although the rate of mutations observed among techniques is different, RAS and BRAF mutational analysis improved prediction of response to anti-EGFR therapy. Additionally, dPCR with a threshold of 1% outperformed the other platforms. Mol Cancer Ther; 16(9); 1999–2007. ©2017 AACR.
Introduction
In metastatic colorectal cancer (mCRC), antibodies against the epidermal growth factor receptor (EGFR), cetuximab and panitumumab, vascular endothelial growth factor (VEGF)‐binding antibodies or other VEGF‐trap therapies, including bevacizumab and aflibercept, as well as regorafenib, a multiple tyrosine kinase inhibitor (1–5), have demonstrated to increase survival. However, the efficacy of these targeted agents in the clinic has been variable and limited, adding toxicity and high costs. In this scenario, it is necessary to identify subgroups of patients more likely to benefit from these therapies.
Exon 2 KRAS mutations were firstly identified as mechanisms of resistance to cetuximab and panitumumab (6, 7). However, factors other than exon 2 KRAS mutations (codons 12 and 13) have been further recognized as mechanisms of intrinsic resistance (8). KRAS mutations in exons 3 and 4 (codons 59–61, 117, and 146) as well as mutations in the same regions of the related oncogene NRAS have been recognized (9) and validated as predictive biomarkers to such anti-EGFR therapies (1, 2).
Although Sanger sequencing has been considered the gold standard in clinical genetic testing, no clearly standardized procedures for RAS mutational testing have been established. An increasing number of hotspots targeted qualitative techniques have been developed with different levels of sensitivity and specificity (10–12), although only some of them are commercially available (13). Digital PCR is a highly sensitive and quantitative method that has been recently applied with the purpose of analyzing the role of minor mutated KRAS/NRAS subclones in patients with mCRC treated with anti-EGFR drugs (14, 15), suggesting that a threshold of 1% of mutated subclones is the optimal cutoff to distinguish patients more likely to benefit from cetuximab or panitumumab. However, data reported on RAS mutations analyzed by BEAMing in the KRAS exon 2 wild-type cohort of patients enrolled in the CRYSTAL trial suggested that mCRC patients with tumors bearing mutations between 0.1% and <5% benefited from the addition of cetuximab to FOLFIRI in the first-line setting (1). There is therefore a need for identifying a mutation detection threshold that is clinically relevant to select patients for anti-EGFR therapy.
We previously reported the mutational profile of the RAS pathway genes in a series of 102 patients with mCRC treated with anti-EGFR based therapy in the refractory setting using two platforms with different analytical sensitivity, a quantitative and highly sensitive digital PCR (Fluidigm platform) and a conventional real-time PCR (Lightcycler 480, Roche; ref. 14). Digital PCR identified a higher number of mutations respect to the real-time PCR, but even more interesting, digital PCR helped us to optimize the selection of patients who had better outcomes after anti-EGFR treatment. Patients with tumors harboring RAS and BRAF mutations in a fraction of <1% seemed to benefit from cetuximab or panitumumab. In the present study, we aimed to confirm our previous results in a larger series of mCRC. Our objectives were to evaluate and compare the sensitivity to detect point mutations in KRAS, NRAS, BRAF, and PIK3CA by means of three different mutational tests and additionally, to evaluate whether a cutoff point could be defined for optimal patient selection.
Materials and Methods
Patients and samples
We retrospectively analyzed formalin-fixed paraffin-embedded (FFPE) primary tumors obtained from surgical resections (partial or complete colectomy) or diagnostic biopsies from 583 patients with mCRC treated with anti-EGFR (n = 255) or bevacizumab (n = 328) regimens from retrospective series and several clinical trials from the TTD/RTICC Spanish network (clinical trial identifiers: NCT01071655, NCT00958386, NCT00885885, NCT01126112, and NCT01704703). Patients received treatment between February 2003 and September 2014. Main clinicopathologic characteristics of the patients are described in Table 1. Three patients were excluded because of invalid results in the mutational analyses. Data on KRAS exon 2 mutational status were available before biological treatment in 135 patients. The therascreen KRAS PCR Kit (Qiagen), cobas KRAS Mutation Test kit, and digital PCR (in two cases included in the clinical trial NCT01704703) were used for the analysis. Pathological diagnosis was verified on the basis of H&E-stained sections, and carcinomas were classified according to UICC TNM classification system. Before DNA extraction, stained slides were reviewed by a pathologist to estimate overall neoplastic cell content and a minimum 15% of tumor cell content was considered to be a valid sample for analysis. Tumor response was extracted from medical records evaluated by attending physician, according to Response Evaluation Criteria in Solid Tumors (RECIST) criteria. Patients with stable disease (SD) or disease progression (PD) were defined as nonresponders. Independent ethics committees at each participating center reviewed and approved the protocol and the study was carried out according to the principles of the Declaration of Helsinki. Committees allowed a waiver in individual patient consent for retrospective biomedical research studies which included the permission to use patients' samples without informed consent in case of death.
Clinicopathologic characteristics of the patients included in the efficacy analysis
DNA extraction and primary PCR reactions
DNA from FFPE tissues was extracted using two comparable methods: the cobas DNA Sample Preparation Kit (Roche Molecular Diagnostics) and the QIAamp DNA Mini Kit (Qiagen) according to the manufacturer's instructions. DNA quality and concentration were measured with a NanoPhotometer (Implen GmbH).
Mutations in KRAS (exons 2/3/4), NRAS (exons 2/3/4), BRAF, and PIK3CA were detected by two commercial assays (Roche cobas test and Qiagen therascreen pyrosequencing kit) and by quantitative nanofluidic dPCR (sensitivity of about 0.05%–0.1%) of genomic DNA as previously described (18). For the cobas analysis, DNA extracted by the cobas kit was used, and for therascreen pyro and dPCR analysis DNA extracted by Qiagen kit was used.
The cobas KRAS Mutation Test kit, the cobas 4800 BRAF V600 Mutation Test, and the cobas PIK3CA Mutation Test were used for detection of somatic mutations in KRAS, BRAF, and PIK3CA genes. The cobas KRAS Mutation Test kit is a PCR assay designed to detect 21 somatic mutations in codons 12, 13 (exon 2), and 61 (exon 3). The cobas 4800 BRAF V600 Mutation Test detects the BRAF V600E mutation and the cobas PIK3CA Mutation Test identifies 17 mutations in exons 1, 4, 7, 9, and 20. In brief, 50 ng of DNA (for KRAS and PIK3CA) or 125 ng of DNA (BRAF) was amplified as well as a negative and positive control, according to the manufacturer's instructions. Mutations were detected using the cobas 4800 System (Roche Molecular Systems Inc.) with automated result interpretation software.
For somatic mutations in NRAS codons 12–13, 59–61, 117, and 146 as well as KRAS codons 117 and 146, we used a probe-based assay, LightMix Kit KRAS-NRAS (TIB MOLBIOL). Six different reactions were performed for each sample. Fifteen microliters of every reaction mix was added to 5 μL of sample or control DNA for a final reaction volume of 25 μL according to the manufacturer's instructions. All reactions were performed using the Roche Diagnostics LigthCycler FastStart DNA Master Hybridization Probe reaction mix and run on a LigthCycler 480 instrument (Roche Molecular Systems Inc.). The thermocycling conditions were as follows: Uracil-DNA glycosylase treatment at 40°C for 20 minutes, denaturation at 95°C for 10 minutes, amplification of target DNA; 60 cycles at 95°C for 5 seconds, 58°C for 5-second acquisition, 58°C for 10 seconds and 72°C for 15 seconds; melting at 95°C for 20 seconds, 58°C for 20 seconds, 40°C for 20 seconds and then a gradual increase in temperature to 85°C with continuous acquisition before cooling at 40°C. Genotype calling based on melting curves was performed using the LightCycler 480 software release 1.5.1 (Roche), and mutations are identified by a melting peak only in case of the presence of a mutation.
The PyroMark Q24 system (Qiagen) was used for pyrosequencing analysis of KRAS, NRAS, and BRAF mutations. In brief, PCR templates for pyrosequencing were amplified from 10 ng of DNA using therascreen KRAS, NRAS, RAS, extension and BRAF Pyro Kit according to the manufacturer's instructions. The biotinylated PCR products were placed in 24-well plates and bound to streptavidin-coated sepharose beads (GE Healthcare). The PCR products were denatured, and the nonbiotinylated fragments were washed from the beads using the Pyromark Q24 Vacuum Workstation (Qiagen). The beads were then resuspended in annealing buffer containing 0.4 pmol/L of the sequencing primer. Raw data files were imported into Pyromark Q24 software (version 2.0; Qiagen) for further analysis following pyrosequencing. For mutations in PIK3CA we used the PI3K Mutation Test Kit (Qiagen), a real-time PCR assay for the detection of 4 somatic mutations in exons 9 and 20. The kit uses a real-time PCR assay based on ARMS PCR technology combined with Scorpions detection technology. The PI3K Mutation Test Kit includes a control assay and 3 mutation assays to assess the total DNA in a sample and the presence or absence of mutated DNA. Twenty nanograms of DNA were amplified as well as negative and positive controls, using the LightCycler 480 software release 1.5.1 (Roche) and mutations were identified calculating ΔCt (Delta Ct) according to the manufacturer's instructions.
The dPCR analysis methods are described in detail in our previous study (14), and specific primers and probes are specified in Supplementary Table S1.
Statistical analysis
We performed a survival analysis using the Cox proportional regression model and calculated the hazard ratio (HR) for progression-free survival (PFS) and overall survival (OS), described in detail in our previous study (14).
Independently from the survival analysis, in order to evaluate the best sensitivity threshold for mutations to predict treatment response in the anti-EGFR cohort, we calculated the likelihood ratio of positive test analysis. In this analysis, we only included patients with evaluable radiological response (N = 223). In order to compare the predictive value of different mutational profiles, we defined three different scenarios for analysis: (i) KRAS (exon 2); (ii) RAS (KRAS exons 2/3/4 and NRAS exons 2/3/4); and (iii) RAS and BRAF.
The concordance between techniques was calculated with the Cohen's kappa coefficient. Data were analyzed with statistical software R version 3.1.2.
Results
Mutational profiling detected by cobas, therascreen pyro, and highly sensitivity dPCR
We performed mutational analyses of KRAS, NRAS, BRAF, and PIK3CA through three different platforms. The average time for testing all genes per 24 samples was 14 hours for cobas, 35 hours for therascreen pyro, and 65 hours for digital PCR. Mutations in any of the genes analyzed were detected in 331 tumors (56.77%) by cobas, in 325 tumors (55.74%) by therascreen, and in 361 tumors (61.92%) by dPCR (Table 2). BRAF mutations were most frequently in KRAS exon 2 wild-type tumors, although both mutations were detected in 7.4% tumors by cobas, 8.3% tumors by therascreen, and 13.3% by dPCR. Specific mutations of individual exons in the population of patients are shown in Supplementary Table S2.
Summary of the mutational status of tumors included in the study population detected by the different techniques
In the subset of patients treated with bevacizumab, we did not detect any mutation in 55 of 328 (16.77%) tumors. No valid results were obtained in 3 tumors. As expected, the majority of tumors harbored mutations in KRAS exon 2, independently of the technique. Extended RAS, BRAF, and PIK3CA analysis increased the number of mutated tumors as referred in Supplementary Table S3.
In the cohort of patients treated with cetuximab or panitumumab, no mutations were detected in 102 of 255 (40%) tumors with none of the three techniques used in the study. There was an invalid result in the cobas test and two invalid results in the dPCR test. One hundred and one of 255 (39.6%) tumors harbored one mutation detected by cobas and 18 of them showed multiple mutations. One hundred of 255 (39.2%) tumors harbored one mutation detected by therascreen pyro and 15 of them showed multiple mutations. Mutations analyzed by both techniques were frequently located in KRAS followed by PIK3CA, NRAS, and BRAF as described in Supplementary Table S3. Analysis by dPCR increased the number of patients bearing mutations up to 119 of 255 (46.7%) and identified multiple additional low-frequency mutant alleles in 34 cases. In this setting, KRAS mutations were further overrepresented (81/255; 31.8%) followed by NRAS (33/255; 12.9%), BRAF (17/255; 6.7%), and PIK3CA (8/255, 3.1%; Supplementary Table S3). Sixty-four tumors scored as positive with the three techniques, and the percentage of mutant alleles by quantitative nanofluidic dPCR ranged from 0.2% to 64.5% (median, 46.3%). The percentage of mutant alleles in the 57 tumors detected only by dPCR ranged from 0.04% to 41.5% (median, 2.1%).
The levels of concordance between techniques for mutation detection are shown in Fig. 1. The concordance was good or very good between cobas and therascreen pyro and also with dPCR when mutant allele fraction (MAF) considered is >10%. These levels decreased slightly considering MAF > 1%. The concordance is even lower when we compared cobas or therascreen with dPCR with a threshold of MAF > 0.1%.
Concordance of mutation detection by cobas, therascreen pyro, and highly sensitive dPCR. A, Cohen kappa index measures the agreement between techniques in mutation detection. B, Venn diagram for the results of three tests together. C, Number of mutations detected by cobas and therascreen pyro. D, Number of mutations detected by cobas and digital PCR. E, Number of mutations detected by therascreen pyro and digital PCR. MAF, mutant allele fraction.
Response rate
Radiological tumor response was evaluable in 514 patients (88%). Complete or partial response occurred in 209 patients (35.8%). Tumor stabilization was observed in 216 patients (37%) and disease progression in 89 patients (15.3%). We then assessed the response rate according to biological therapy. In the group of patients treated with bevacizumab, the response rate was not associated with RAS/BRAF mutational status.
In the group of patients treated with anti-EGFR therapy, radiological tumor response was evaluable in 223 of the 255 patients. Regarding the treatment administered in this specific cohort, 30 patients (13.45%) received cetuximab or panitumumab monotherapy (23 of them were chemonaïve for metastatic disease), 75 patients (33.63%) received an irinotecan + anti-EGFR regimen and 118 patients (52.91%) were treated with a fluoropyrimidine-based chemotherapy regimen (78 of them in the first-line or neoadjuvant/post-metastasectomy setting). Overall response rate (ORR) in this subgroup of patients was 40.3%. We further evaluated response rate according to mutational profile (Table 3). KRAS exon 2 wild-type population identified by cobas test presented an ORR of 42.7% [76 responders (89.4%) received anti-EGFR in combination with chemotherapy—19 and 24 patients in the neoadjuvant/post-metastasectomy and first-line setting, respectively—and 9 responders were treated with anti-EGFR monotherapy). Four of the KRAS exon 2 mutated patients responded to anti-EGFR regimens [chemotherapy and anti-EGFR in the neoadjuvant/post-metastasectomy setting (3/4) or in first line]. Extended RAS identified 16 additional RAS mutant mCRC patients with tumor response, 10 of them treated with anti-EGFR in combination with chemotherapy in the neoadjuvant or first-line setting and 2 of them treated only with panitumumab in first line. BRAF genotyping identified 8 additional mutated cases, 2 of them responders (both treated in combination with chemotherapy in first line). RAS/BRAF analysis by cobas increased the response rate up to 45.8% and 26.9% in wild-type and mutated populations, respectively. Analysis of KRAS exon 2 by therascreen pyro kit identified 24 mutated cases, 5 of them responders (one of them treated with anti-EGFR monotherapy in first line, two of them treated in combination with chemotherapy in the neoadjuvant and post-metastasectomy setting, and the other 2 patients treated in second and third line with chemotherapy). This reclassification translated into an ORR of 42.7% in the KRAS exon 2 wild-type subgroup [73 responders (91.2%) received a combination of chemotherapy and anti-EGFR—18 and 24 patients in the neoadjuvant/post-metastasectomy and first-line setting, respectively—and 11 responders received cetuximab or panitumumab alone]. Extended RAS identified 15 additional mutated tumors with radiological response (5 and 6 patients treated in first-line and neoadjuvant/post-metastasectomy setting, respectively). BRAF genotyping identified 11 additional mutated cases (2 of them responders and treated in first line). We have to point out that only one of the RAS/BRAF mutated tumors that responded was treated with anti-EGFR monotherapy. Extended RAS and BRAF panel by therascreen increased response rate up to 48.3% and 25% in wild-type and mutated populations, respectively. Qualitative dPCR genotyping (limit of detection of the technique, 0.05%–0.1% of mutant alleles) of KRAS exon 2 translated into an ORR of 41.7% in the KRAS exon 2 wild-type population. Sixty-seven KRAS exon 2 wild-type responders (91.8%) received a combination of anti-EGFR and chemotherapy (38 of them as neoadjuvant/postmetastasectomy or first-line therapy) and 6 responders received anti-EGFR monotherapy (5 of them in the first-line setting). The extended RAS genotyping reclassified 39 additional mutated cases. Nine of these patients responded to treatment (everyone received anti-EGFR monotherapy), that was administered as first-line or neoadjuvant chemotherapy in 5 of them. Finally, the BRAF analysis identified 11 additional mutated cases, 3 of them responders (all of them treated with anti-EGFR in the first-line setting). The ORR in the RAS/BRAF wild-type subset increased up to 48.8%, whereas ORR in patients reclassified as RAS/BRAF mutated by qualitative dPCR was 29.6%. It is noteworthy that most of the RAS/BRAF-mutated tumors classified as responders were treated with anti-EGFR and chemotherapy combination (26 out of 29 cases), as described with cobas and therascreen pyro. Moreover, approximately 90% of the KRAS exon 2 wild-type patients with a partial or complete response received cetuximab or panitumumab in combination with chemotherapy.
ORR in the anti-EGFR group (N = 223)
An inverse correlation between the proportion of mutant allele and anti-EGFR response is shown in Fig. 2 (P < 0.0001). The median percentage of mutated DNA was 12.5% for responders (range, 0.1%–50%) and 34.2% for nonresponders (range, 0.3%–65%) and of note in 10 of 29 responders (34%) that have been scored as positive considering mutated cases by qualitative dPCR genotyping, the major MAF was below 1%.
Box plot representing the correlation between the proportion of major mutant alleles and the anti-EGFR response In the y axis: the percentage of major mutated allele (RAS + BRAF). In the x axis: patients with an evaluable response (N = 223). CR, complete response; PR, partial response; SD, stable disease; PD, progressive disease. The mean percentage of mutated alleles in each group is shown by a line (the percentage of mutant alleles are shown in Supplementary Table S3).
The likelihood ratio of positive test analysis was performed to estimate the optimal clinically relevant threshold cutoff value of mutant alleles in the prediction of response. For dPCR, the categorization of a fraction of 1% or higher for any mutant allele offered the best likelihood ratio of positive test for all combinations of RAS and BRAF analysis (Supplementary Table S4) delivering the best balance between sensitivity and specificity. The ORR using the 1% cutoff in the population with KRAS exon 2 wild-type tumors was 43.2%, whereas ORR in patients reclassified as KRAS exon 2 mutated tumors was 24.2% (Table 3). The extended RAS genotyping reclassified 37 additional mutated cases (8 responders out of 37) which translate into a higher ORR (48.4%) in the RAS wild-type cases. The ORR among RAS extended mutated reclassified group was 22.9%. Finally, the BRAF analysis identified 12 additional mutated cases. The ORR in this RAS/BRAF wild-type subset was 50.4%, whereas ORR in patients reclassified as RAS/BRAF mutated by qualitative dPCR with 1% cutoff was 23.2%. When compared with qualitative detection of dPCR, the use of 1% cutoff reclassified as wild-type 16 patients: 10 responders and 6 nonresponders. Of the 6 nonresponder cases now reclassified as wild-type, 5 had stable disease while in the remaining 1 disease progressed. We have also calculated ORR with different cutoffs (1% to 5%) but the 1% is the best scenario for the ORR compared with the others (Supplementary Table S5).
Progression-free survival and overall survival
In the overall study population, OS was significantly shorter in the subgroup of patients with RAS or RAS/BRAF mutations assessed by any of the three platforms (Supplementary Table S6). However, a clear and significant benefit to anti-EGFR treatment in terms of both PFS and OS was observed only in patients with wild-type tumors (Table 4). Patients with KRAS exon 2 mutated tumors analyzed by cobas had a significantly increased risk of disease progression and death with an HR of 1.79 (95% CI, 1.12–2.85) and 2.06 (95% CI, 1.27–3.34), respectively. OS was also associated with KRAS exon 2 mutation assessed by dPCR. The risk of disease progression was higher when any mutation assessed by any of the three different techniques was considered (HR of 1.50, P = 0.009, HR of 1.87, P ≤ 0.001 and HR of 1.54, P = 0.003 for cobas, therascreen, and dPCR, respectively). This association was maintained in the OS analysis.
Survival analyses according to mutational status through different techniques in the anti-EGFR cohort
Interestingly, when we analyzed survival according to different cutoffs, patients with KRAS exon 2 mutated tumors with more than 1% of MAF assessed by quantitative nanofluidic dPCR presented an HR of 1.51 (95% CI, 1.02–2.24) for PFS and 1.76 (95% CI, 1.14–2.72) for OS, compared with wild-type patients or with less than 1% of mutant allele fraction. When we extended the analysis to RAS and BRAF and a threshold of 1%, the 97 mutated cases presented an HR of 1.66 (95% CI, 1.23–2.24) and 2.01 (95% CI, 1.41–2.85) for PFS and OS, respectively. In concordance with likelihood ratio of positive test analysis depicting 1% as the best cutoff value, a 5% cutoff resulted in a poorer discrimination of PFS and OS.
Additionally, we differentiated patients treated with anti-EGFR in two groups: patients treated with anti-EGFR in first-line setting and patients treated with anti-EGFR in second-line or beyond, and evaluated the impact of gene mutations on survival (Supplementary Table S7). Eighty-three patients (32.6%) received cetuximab or panitumumab (60 of them in combination with chemotherapy) as first-line treatment and were defined as “chemosensitives,” and 99 patients (38.8%) received anti-EGFR drugs after first-line treatment failure (89 of them in combination with chemotherapy). KRAS exon 2 mutations had no impact on survival in the chemosensitive group. Only RAS mutations assessed by dPCR with a threshold of 1% increased the risk of progression compared with wild-type tumors or with less than 1% of MAF (HR = 1.79, 95% CI, 1.05–3.05). We observed similar results when we extended the mutational analysis to RAS and BRAF. Once again, the cutoff of 1% better discriminated patients likely to benefit from anti-EGFR therapy with an HR of 2.03 (95 % CI, 1.20–3.43). In terms of OS, RAS and BRAF mutations assessed by cobas, therascreen, or dPCR were significantly associated with poor outcome. Similar results were observed in the group of patients treated with anti-EGFR drugs after first-line tumor progression, although in this setting, RAS and BRAF mutations were related to worse PFS, independently of the technique used.
In the subgroup of patients treated with bevacizumab, only BRAF mutations detected by therascreen and dPCR were associated with shorter overall survival.
Discussion
In this retrospective study, we have analyzed the mutational status of RAS pathway genes through different techniques to evaluate and compare the sensitivity of each one. Furthermore, we have demonstrated that anti-EGFR benefit is more evident in patients with RAS and BRAF wild-type tumors.
Nowadays, expanded RAS and BRAF analyses should be performed in mCRC patients to establish the best therapeutic strategy. Several commercial platforms are used in clinical practice for RAS testing, although none of them is considered the standard of care (13). We have used three different techniques, the commercial cobas and therascreen pyro assays and the digital PCR by Fluidigm, for the analysis of mutations in KRAS, NRAS, BRAF, and PIK3CA, and we then compared the sensitivity and predictive capability to anti-EGFR- and bevacizumab-based regimens. As expected, the most common gene mutation was located in KRAS exon 2. We have to point out that the panel of hotspots analyzed by each technique was slightly different. In this context, it was expected that the frequency of gene mutations differed between techniques. Nevertheless, dPCR identified mutations in approximately 62% tumors compared with 55% of mutated tumors identified by cobas or therascreen due to its higher analytical sensitivity. dPCR is a highly sensitive and quantitative method with a limit of detection or sensitivity < 0.1%. Because of this higher technical sensitivity, dPCR increased the mutation rate up to 62% with a range of MAF between 0.04% and 85.4%. When we analyzed the concordance between tests, we observed a high agreement (Cohen's kappa ≈ 0.90) if we considered KRAS exon 2 and a limit of detection of 10%. However, the level of agreement decreased when we compared cobas and therascreen with dPCR with a limit of detection of 1%–10% or <1% and is even lower when we included in the mutational analysis RAS (KRAS exons 3, 4 and NRAS exon, 2, 3, 4) and BRAF. Taking into account that there was no bias in tumor sample selection because we selected one paraffin block per patient with a minimum of tumor cells for DNA extraction and subsequent DNA mutational analysis by the three techniques, our results demonstrated that quantitative dPCR is more sensitive than cobas and therascreen for RAS testing.
From a practical point of view, the average time for DNA extraction was similar for the two methods used and the amounts of DNA required for the three tests were also very similar. However, we found differences between the three techniques regarding the turnaround time: Cobas significantly outperformed the others and the slowest one was the dPCR which needs an improvement in the automation to increase the competitiveness. Finally, costs of the reagents per case analyzed by the three platforms were very similar.
In terms of treatment efficacy, we have observed that the clinical outcomes of patients treated with bevacizumab are independent of RAS mutational status, as previously reported (16). We have also confirmed that RAS mutational status is associated with anti-EGFR treatment outcome in the first-line setting and beyond. This is in line with the results published from different clinical trials in metastatic disease. Pivotal studies with anti-EGFR drugs have analyzed the mutational status of RAS by different techniques (13). The PRIME study used Sanger sequencing and WAVE-based SURVEYOR Scan Kits from Transgenomic for RAS and BRAF testing and detected 52% RAS-mutated tumors (2). RAS mutational analysis in the Crystal study was performed by BEAMing and detected mutations in 43% of tumors (1). Similar RAS mutation rates are reported in the literature, but it is remarkable that the analytical sensitivity of the tests is around 5% to 10%. Highly sensitive methods such as digital PCR or BEAMing are able to identify minor mutated sub-clones present in a low proportion (limit of detection < 0.01%; refs. 17, 18). Both methods have demonstrated to improve the selection of patients more likely to benefit from anti-EGFR therapies; however, a discrepancy in the definition of the optimal cutoff for anti-EGFR treatment was observed (1, 14, 15). To identify the fraction of mutated allele that correlates with resistance to anti-EGFR agents is still a challenge. Our results have confirmed the findings previously reported by Laurent-Puig and colleagues and our own group (14, 15). We have observed that RAS mutational status assessed by dPCR with a threshold of 1% identifies candidates for cetuximab or panitumumab with better clinical outcomes in terms of ORR, PFS, and OS than any of the platforms and sensitivity thresholds, confirming the results we reported from a smaller series of chemorefractory mCRC (14). Overall, the response rate in the RAS/BRAF wild-type group of patients increased around 5% in both series when we compared the dPCR with a cutoff of 1% with a conventional analytical method (either real-time PCR or cobas). However, we need to be cautious because a relevant percentage of patients with mutated tumors (around 20% across different analyses) responded to anti-EGFR regimens. Recent publications of clinical randomized trials have reported results on response rate in the mutated RAS tumors population treated with anti-EGFR and chemotherapy with ORR of 34% to 15% (1, 19). This can be explained by the chemotherapy effect in most of them but not all, and is a fact that must stimulate continuous research for finding better prediction models beyond RAS mutational analyses (20).
A number of limitations of our study need to be highlighted. First, there was no validation cohort to validate the mutation rates for each of the three tests, their concordance, and the prediction of disease outcome. Second, this is a retrospective and pooled analysis of a highly heterogeneous series of patients treated with different chemotherapy regimens and biologicals. Finally, we have analyzed only patients treated with bevacizumab or anti-EGFR containing regimens in a nonrandomized fashion. We have therefore only analyzed clinical outcomes of those patients with tumors with and without mutations assessed by different methods. To validate our results, a similar analysis should be performed in prospective and randomized clinical trials (20). Tumors from patients included in the CAPRI GOIM trial were tested for RAS, BRAF, and PIK3CA by next-generation sequencing and suggested that fraction of neoplastic cells carrying a specific molecular alteration (RAS, BRAF, PIK3CA) could correlate with the level of resistance to anti-EGFR agents (21). Additionally, we have recently completed accrual of the prospective TTD ULTRA clinical trial (clinical trial identifier: NCT01704703) which will provide relevant information about the utility of highly sensitive methods in DNA testing for clinical practice.
Disclosure of Potential Conflicts of Interest
B. Massuti has received speakers bureau honoraria from Roche and Merck and is a consultant/advisory board member for Roche. M. Valladares is a consultant/advisory board member for Merck Serono. J.M. Vieitez reports receiving a commercial research grant from Roche and Amgen and is a consultant/advisory board member for Amgen, Roche, Servier, and Shire. J. Tabernero is a consultant/advisory board member for Amgen, Bayer, Pfizer, Roche, Sanofi, Symphogen, Taiho, Takeda, Boehringer-Ingelheim, Celgene, Chugai, Genentech, Lilly, MSD, Merck Serono, and Novartis. G. Capella reports receiving a commercial research grant member for Roche. No potential conflicts of interest were disclosed by the other authors.
Disclaimer
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Authors' Contributions
Conception and design: C. Santos, D. Azuara, E. Aranda, J. Tabernero, G. Capella, R. Salazar
Development of methodology: C. Santos, D. Azuara, E. Elez, B. Massuti, J. Tabernero, G. Capella, R. Salazar
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): C. Santos, D. Azuara, R. Garcia-Carbonero, P. Garcia Alfonso, A. Carrato, E. Elez, A. Gomez, F. Losa, C. Montagut, B. Massuti, M. Varela, M. Valladares, J.L. Manzano, J.M. Vieitez, X. Sanjuan, J. Tabernero, R. Salazar
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): C. Santos, D. Azuara, R. Garcia-Carbonero, P. Garcia Alfonso, C. Montagut, B. Massuti, V. Navarro, A. Lopez-Doriga, V. Moreno, E. Aranda, J. Tabernero, R. Salazar
Writing, review, and/or revision of the manuscript: C. Santos, D. Azuara, R. Garcia-Carbonero, P. Garcia Alfonso, A. Carrato, E. Elez, F. Losa, B. Massuti, A. Lopez-Doriga, V. Moreno, M. Valladares, E. Aranda, X. Sanjuan, J. Tabernero, R. Salazar
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): C. Santos, D. Azuara, F. Losa, R. Salazar
Study supervision: C. Santos, D. Azuara, P. Garcia Alfonso, A. Carrato, J.M. Vieitez, E. Aranda, J. Tabernero, G. Capella, R. Salazar
Grant Support
Roche Diagnostics funded this study. This study was supported by the Instituto de Salud Carlos III, organismo adscrito al Ministerio de Economía y Competitividad and “Fondo Europeo de Desarrollo Regional (FEDER), una manera de hacer Europa” grant (FIS PI12-01589; to R. Salazar); RTICC: Grupo Cancer Digestivo – Instituto de Salud Carlos III, CIBERONC-CIBERESP; Catalan Health Institute and Autonomous Government of Catalonia DURSI grants 2014SGR1494 (to Ramon Salazar and Josep Tabernero), 2014SGR647 (to Victor Moreno), 2014SGR338 (to Gabriel Capella).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Acknowledgments
The authors would also like to thank Inma Ruiz de Mena and Susana Rodríguez (TTD Data Center) and Dr. Marga Nadal (Research Management Unit, ICO).
The physicians listed below cared for the patients in this study. The authors thank them for their cooperation and support:
Study Chairs: Dr. R Salazar and Dr. J. Tabernero (Barcelona). Central Lab: Dr. G. Capellà and Dr. D. Azuara. Department of Pathology: Dr. X. Sanjuán Garriga and Dr. M. Varela
Dr. Salazar and Dr. Santos-ICO. H. Duran i Reynals (Barcelona); Dr. García Carbonero-H Virgen del Rocío (Sevilla); Dr. García Alfonso-H. Gregorio Marañón (Madrid); Dr. Carrato and Dr. Guillén-H. Ramón y Cajal (Madrid); Dr. Tabernero and Dr. Élez-H. Universitari Vall d'Hebrón (Barcelona); Dr. Aranda -H. Universitario Reina Sofía (Córdoba); Dr. Losa-H. Sant Joan Despí - Moisés Broggi (Barcelona); Dr. Montagut-H. del Mar (Barcelona); Dr. Massutí-H. General Universitario (Alicante); Dr. Valladares and Dr. Reborerdo-C. H. Universitario (La Coruña); Dr. Manzano and Dr. Layos-ICO. H. Germans Trias i Pujol (Barcelona); Dr. Viéitez-H. Universitario Central de Asturias; Dr. Gallego-H. General Universitario de Elche (Alicante); Dr. Robles-H. 12 de Octubre (Madrid); Dra. Safont-H. General Universitario (Valencia); Dr. Falcó-H. Son Llatzer (Palma de Mallorca); Dr. Alonso Orduña-H. Miguel Servet (Zaragoza); Dr. Remón-H. de Mataró (Barcelona); Dr. Escudero-H. C. Universitario Lozano Blesa (Zaragoza); Dr. López -H. Marqués de Valdecilla (Santander); Dr. Pericay-C. S. Parc Taulí (Barcelona); Dr. Conde and Dr. Glez Flores-H. Virgen de las Nieves (Granada); Dr. Páez-H. Santa Creu i Sant Pau (Barcelona); Dr. Provencio-H. Puerta de Hierro (Madrid); Dr. Salud-H. de Lleida Arnau de Vilanova (Lérida); Dr. Llanos, Alonso-H. Universitario de Canarias (Tenerife); Dr. Hernández Yagüe-ICO. H. U. Josep Trueta (Gerona); Dr. Sastre-H. Universitario Clínico San Carlos (Madrid); Dr. Méndez -H. de Móstoles (Madrid); Dr. Vera-C. H. de Navarra.
Footnotes
Note: Supplementary data for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/).
- Received February 16, 2017.
- Revision received May 3, 2017.
- Accepted June 6, 2017.
- ©2017 American Association for Cancer Research.