The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers.

Colorectal tumours that are wild type for KRAS are often sensitive to EGFR blockade, but almost always develop resistance within several months of initiating therapy. The mechanisms underlying this acquired resistance to anti-EGFR antibodies are largely unknown. This situation is in marked contrast to that of small-molecule targeted agents, such as inhibitors of ABL, EGFR, BRAF and MEK, in which mutations in the genes encoding the protein targets render the tumours resistant to the effects of the drugs. The simplest hypothesis to account for the development of resistance to EGFR blockade is that rare cells with KRAS mutations pre-exist at low levels in tumours with ostensibly wild-type KRAS genes. Although this hypothesis would seem readily testable, there is no evidence in pre-clinical models to support it, nor is there data from patients. To test this hypothesis, we determined whether mutant KRAS DNA could be detected in the circulation of 28 patients receiving monotherapy with panitumumab, a therapeutic anti-EGFR antibody. We found that 9 out of 24 (38%) patients whose tumours were initially KRAS wild type developed detectable mutations in KRAS in their sera, three of which developed multiple different KRAS mutations. The appearance of these mutations was very consistent, generally occurring between 5 and 6 months following treatment. Mathematical modelling indicated that the mutations were present in expanded subclones before the initiation of panitumumab treatment. These results suggest that the emergence of KRAS mutations is a mediator of acquired resistance to EGFR blockade and that these mutations can be detected in a non-invasive manner. They explain why solid tumours develop resistance to targeted therapies in a highly reproducible fashion.

One major barrier to testing any hypothesis about the nature of acquired resistance to anti-EGFR antibodies is limited access to post-treatment tumor tissue. Even when post-treatment tumor tissue is available, sampling bias confounds interpretation because only a small portion of one tumor is usually biopsied, precluding assessment of genetic heterogeneity within or among lesions. To circumvent the tissue access problem, we have examined circulating, cell-free DNA - a form of “liquid biopsy”. It has been previously shown that circulating tumor DNA (ctDNA) can be found in the majority of patients with metastatic colorectal cancers[7-9]. Analysis of ctDNA is informative because it not only can identify a specific mutant genotype but can also provide a measurement of the total tumor burden[7]. If tumors became resistant to anti-EGFR antibodies as a result of the emergence of KRAS mutations in their tumors, we expected that mutant KRAS genes would be released into the circulation in a time frame consistent with the emergence of resistance.

We retrospectively analyzed longitudinal serum samples from 28 patients with chemorefractory metastatic colorectal cancer (CRC) receiving single-agent therapy with panitumumab[10]. Four patients with KRAS mutant tumors, who never achieved disease control, were selected as controls. As expected, these four patients were found to have progressive disease at the time of first tumor assessment, 7 ± 2 weeks (mean ± 1 standard deviation) after initiating treatment with panitumumab (Supplementary Table 1)[1,2]. The other 24 patients with WT KRAS tumors achieved a partial response (n=8), had prolonged stable disease (n=14), or had retrospectively-determined progressive disease but remained on study for an extended period (n=2). These 24 patients developed clinically evident progressive disease 23 ± 10 weeks (mean ± 1 standard deviation) following initiation of treatment (Supplementary Table 1) as determined by radiographic imaging.

Serum samples obtained from patients prior to the initiation of therapy were evaluated for all common mutations at codons 12 and 13 of KRAS using a digital ligation assay with a detection limit of one mutant molecule per ml of serum (examples in Supplementary Fig. 1)[11]. Mutations were independently confirmed in a second aliquot of the same serum and the results quantified via a PCR assay that can digitally enumerate the fraction of rare variants in a complex mixture of DNA template molecules (examples in Supplementary Fig. 1 and Supplementary Table 2)[12].

Of the four cases whose archival tumors harbored KRAS mutations, three had detectable levels of mutant KRAS in the serum prior to treatment with panitumumab (Supplementary Table 2). In these three patients, the KRAS mutations found in the circulation were identical to those found in the patients’ tumor tissues even though the time of serum assessment was, on average, 88 weeks after the diagnosis of metastatic disease and even longer after the initial tumor excision (Supplementary Tables 1 and 2). No mutations in KRAS were detected in the pre-treatment serum DNA from patients whose archival tumor tissue was WT for KRAS (Supplementary Table 2).

Next, we examined 169 serially acquired serum samples from the 28 patients for the presence of mutant KRAS fragments (Supplementary Fig. 2). These samples were collected at approximately four week intervals until disease progression (Supplementary Table 2). Serum was assessed for mutations at KRAS codons 12 and 13 as described above. When sufficient serum was available (23 of 28 patients), it was assessed for BRAF mutations at codons 600 and 601 using the identical assay. Of the 24 patients who did not have KRAS mutations at baseline, nine (38%) were found to develop KRAS mutations during the course of therapy (Supplementary Table 2) while none developed BRAF mutations. In three cases (Patients #1, 12, and 22), multiple KRAS mutations appeared in the circulation - two different mutants in one case and four different mutants in the two other cases (examples in Fig. 1). In each of these cases, the time to appearance of all mutations in the circulation was very similar (Fig. 2). Circulating mutant KRAS templates were identified prior to radiographic evidence of disease progression in three of the nine cases (Patients #1, 10 and 24). The lead time (i.e., the interval between detectable ctDNA and radiographic evidence of disease progression) averaged 21 weeks (Supplementary Table 2). The level of ctDNA generally paralleled that of CEA, the standard biomarker used for following disease progression in metastatic CRC (Supplementary Table 3 and Supplementary Fig. 2).

Fig. 1

Emergence of circulating mutant KRAS.

The time course of circulating mutant KRAS alleles, CEA and tumor burden are depicted in two patients where fragments of circulating DNA containing mutant KRAS were detected. (A) Demonstrates the emergence of four different mutant KRAS alleles in codon 12 (cDNA nt 35T, 34T, 35C and 34C) in the serum of Patient #1 and (B) demonstrates the emergence of two different mutant KRAS alleles in codon 12 (cDNA nt 34T and 35C) in the serum of Patient #12. Tumor burden refers to the aggregate cross-sectional diameter of the index lesions.

Fig. 2

Predicted probability distribution of times from when treatment starts until resistance mutations become observable in circulating DNA. Overlaid is the observed time to detection of mutant KRAS fragments (23 ± 5 weeks, mean ± 1 standard deviation) and time to clinical progression (25 ± 10 weeks, mean ± 1 standard deviation) in the studied patients. Predictions were based on the Lea-Coulson model with death, as introduced by Dewanji et al.[14], or equivalently, the branching process model of Iwasa et al[14,16].

In the three patients (Patients #8, 20 and 28) with detectable mutant KRAS in their tumors as well as in their circulation, no new mutations in KRAS emerged (Supplementary Table 2).

The progression-free survival (Hazard ratio 0. 9; 95% CI 0.3366 to 2.453; p=0.85; Log-Rank Test) and overall survival (Hazard ratio 0.42; 95% CI 0.1599 to 1.144; p=0.09; Log-Rank Test) of patients was not significantly different whether they developed secondary KRAS mutations or not (Supplementary Fig. 3)

The availability of serially collected serum samples provided a unique opportunity to model the tumor evolutionary process in those patients that responded to panitumumab. The first question we addressed was whether KRAS mutations were likely to be present prior to the initiation of therapy with panitumumab. To estimate the average tumor growth rate, we used the ctDNA data (Supplementary Table 2) from the nine patients that developed KRAS mutations during therapy. The average tumor growth rate was found to be 0.069 - i.e., the number of tumor cells resistant to panitumumab doubled approximately every ten days (=ln2/0.069 days; see Supplementary Appendix). The growth rate represents the difference between the cell birth rate b and the cell death rate d, i.e., b – d. Previous studies[13] have shown that b = ~0.25 for colorectal cancer cells, corresponding to one cell division every four days, yielding a value of 0.18 (=0.25-0.069) for d.

Using these data-derived estimates of b and d, a branching process model was used to test the null hypothesis that there were no cells with KRAS mutations in the tumors prior to the initiation of therapy. We calculated the probability that the number of mutant KRAS cells could have grown to at least the observed levels at the times they were measured, assuming this null hypothesis (see Supplementary Appendix). Using the known tumor burdens and pre-treatment ctDNA levels measured in Patients 8 and 20 who harbored KRAS mutations in their tumors prior to therapy, as well as data obtained on previously studied patients with metastatic disease[7], we calculated that one mutant KRAS template per ml of serum corresponds to a tumor containing 44 million KRAS-mutant cells. We performed the statistical valuation separately for each patient. In the three patients that developed more than one circulating KRAS mutation, we assumed that each of the detected KRAS mutations originated from a different lesion (see below). Thus a total of 16 lesions from nine patients were analyzed. For each lesion, the possibility that the observed mutations were absent at the start of treatment could be rejected (thus confirming the presence of pre-existing KRAS mutations) with >95% confidence (in most cases, >99.9% confidence; Supplementary Table 4). [14,15] Furthermore, we varied the birth rate b around the previously calculated value of 0.25, allowing it to be as small as 0.15 or as large as 0.35, and found that this did not affect our conclusion that the mutations were present prior to therapy, though for one lesion the confidence interval fell below the 95% threshold as the birth rate was increased (Supplementary Table 4).

Next, we estimated how many mutations gave rise to resistance to panitumumab in patients in vivo. The total number of lesions large enough to detect by imaging averaged 7.0 lesions/patient in the 24 patients studied with WT KRAS tumors (Supplementary Table 1). Assuming that each individual KRAS mutation detected in the serum emanated from a single metastatic lesion, the 16 lesions responsible for contributing ctDNA (Supplementary Table 2) accounted for 9.6% of the total 167 lesions. The other 90% of lesions presumably had developed genetic alterations that were not assessed in our study. As we assessed four different mutated nucleotides (nt) of KRAS (at codons 12 and 13, Supplementary Table 5), this analysis suggests that a total of ~42 (= 4 × 167/16) nt gave rise to resistance in our 24 patients. These remaining 38 nt (42 minus 4) would represent a maximum of 38 genes (one potential mutant nt per gene), or ~ 10 genes if, like KRAS, there was an average of four nt per gene that could give rise to resistance when mutated.

Would at least one of these 42 mutant nt be expected to be present in a metastatic lesion initiated by a single cell that was WT at all 42 nt? We analyzed this question using a Luria-Delbruck model generalized to incorporate cell death[14-16]. If we conservatively assume that cancer cells have the same mutation rate as normal cells and the same birth rate described above (0.25)[15,17], then a tumor detectable by CT-scanning (one billion cells) is almost certain (probability > 1 × 10−32) to contain at least one cell with a mutation at one of these 42 nt. The expected number of such cells is ~3,200 (Supplementary Appendix). [16,17] These 3,200 cells are distributed among various clones that arose during the growth of the metastatic lesion prior to therapy with panitumumab. The largest clonal subpopulation is in most cases the progeny of the subpopulation containing the first resistance mutation to arise and survive stochastic drift. Using the equations described in the Supplementary Appendix, we calculated the expected size of this clonal subpopulation as 2,200 cells; thus, 69% (2,200/3,200) of the resistant cells should be derived from a single clone. As described in the Supplementary Appendix, varying the number of possible mutations conferring resistance from 4 to 100 did not change this expectation: there are hundreds to thousands of cells with resistance mutations in each metastasis prior to the initiation of panitumumab therapy and more than half of these resistant cells are expected to be derived from a single clone.

Our mathematical model also accounts for the very similar times at which circulating KRAS fragments were observed across patients and lesions. We calculated the probability distribution for the times at which circulating KRAS fragments should become detectable (Supplementary Appendix). As shown in Fig. 2, the mean of this distribution was 22 weeks with a 95% confidence interval ranging from 18 to 25 weeks. These predicted times are strikingly similar to those actually observed in the patients analyzed in this study (23 ± 5 weeks, mean ± one standard deviation) and when clinical progression was observed (25 ± 10 weeks, mean ± one standard deviation). Finally, we performed simulations to further validate these analytical findings. The results were similar: mutant KRAS fragments were predicted to become evident at 22 ± 1.5 weeks (mean ± one standard deviation) following treatment (Supplementary Appendix).

Our study indicates that the resistance mutations in KRAS and other genes were highly likely to be present in a clonal subpopulation within the tumors prior to the initiation of panitumumab therapy. Though we base this conclusion on new data and mathematical analyses, it is consistent with prior experimental data on other targeted agents [18-21] as well as theoretical predictions[22-24]. As noted in several previous studies of targeted therapeutic agents [18-21], most lesions recur at approximately the same time following therapy. This was also true in our selected group of patients where disease progression occurred at the same time in patients who did or did not develop circulating KRAS mutations (Fig. 2).

Our study suggests that only a limited number of nucleotide mutations and genes are likely to be able to exert a resistance phenotype. If there were many more genes with this capacity, then it would be unlikely that such a high fraction of patients (38%) and lesions (9.6%) would develop mutations in a single gene (KRAS). This finding is consistent with prior studies in that only a small number of genes has been found to be expressed differentially in resistant tumors and an even smaller number to be mutated - none as frequently as KRAS[18,25,26]. That only a small number of genes can confer resistance is encouraging; if there were a very large number of genes (and nucleotides) that had the capacity to produce resistance to panitumumab, there would be little hope of combining this drug with others to circumvent resistance.

In sum, our results suggest the following scenario for the development of resistance to panitumumab. Each relatively large metastatic lesion is expected to contain a subclone containing hundreds or thousands of cells with one of ~42 mutations conferring resistance to the antibody. Resistance is therefore a fait accompli - the time to recurrence is simply the interval required for the subclone to re-populate the lesion. This generally takes five to six months (Fig. 2) and is due to the rapid expansion of the resistant subclone immediately following treatment initiation. To make these remissions last longer than five to six months, combination therapies targeting at least two different pathways will be required.

METHODS

Patients and specimens

Patients with chemorefractory metastatic colorectal carcinoma were enrolled into one of two panitumumab monotherapy studies (NCT00089635 and NCT00083616)[10]. The study protocols were approved by the institutional review boards and all patients signed a written consent form. A subset of patients was selected for this analysis from a total of 388 patients. Patients received panitumumab 6mg/kg every two weeks until disease progression. Tumor scans were read centrally by a panel of at least two blinded independent radiologists using a modification of the WHO criteria. Assessments were performed at four week intervals through week 28 and every three months thereafter until progression of disease. Responses were confirmed at least four weeks after response criteria were first met. KRAS mutational status in the tissue was predetermined using the DxS assay (Qiagen).

DNA purification

DNA was purified from of 0.2–1mL of banked serum using the QIAamp circulating nucleic acid kit (Qiagen, cat. # 55114) following the manufacturer’s recommendations. Amplifiable DNA was quantified with qPCR, employing primers and conditions as previously described[27].

Assessments of circulating KRAS-mutant DNA

The ligation assays were performed as previously described[11] using the primers and probes indicated in Supplementary Table 5. In brief, KRAS fragments containing codons 12 and 13 were amplified with primers designed to yield a small PCR product (106 bp) to accommodate the degraded DNA found in serum[28]. Note that some of the probes contained locked-nucleic acid (LNA) linkages (Exiqon). To confirm and further quantify samples containing KRAS mutations, BEAMing assays were employed[12] using the primers and probes listed in Supplementary Table 5. The number of mutant fragments per ml serum was determined from the fraction of alleles containing the mutant allele (determined by BEAMing[7]) and the number of alleles assessed per ml of plasma (determined by qPCR[27] and reported in Supplementary Table 6).