Abstract
High-grade serous ovarian carcinoma (HG-SOC) is a major cause of cancer-related death. Whether treatment order—primary debulking surgery followed by adjuvant chemotherapy (PDS) or neo-adjuvant chemotherapy with interval surgery (NACT)—affects outcome is controversial. We developed a mathematical framework that holds for hierarchical or stochastic models of tumor initiation and reproduces HG-SOC clinical course. After estimating parameter values, we infer that most patients harbor chemo-resistant HG-SOC cells at diagnosis, PDS is inherently superior to NACT due to better depletion of resistant cells, and earlier diagnosis could improve survival of primary, but not relapsed, patients. Our predictions are supported by primary clinical data from multiple cohorts. Our results have clear implications for these key issues in HG-SOC management.
Significance Statement The optimal order and timing of surgery and chemotherapy, and the potential benefits of earlier diagnosis of HG-SOC, remain controversial. We developed a mathematical framework of tumor dynamics to address such issues, populated the model with primary clinical data and reliably recapitulated clinical observations. Our model prospectively predicts that: (1) PDS is inherently superior to NACT; (2) timely adjuvant chemotherapy is critical for the outcome of PDS with <1mm, but not >1mm, residual tumors; (3) earlier detection of relapse is unlikely to be beneficial with current therapies; (4) earlier detection of primary HG-SOC could have substantial benefit. Our model provides insights into the evolutionary dynamics of HG-SOC, argues for new clinical trials to optimize HG-SOC therapy, and is potentially applicable to other tumor types.
Introduction
Ovarian cancer is the 8th most common cancer in women and the 8th most common cause of cancer death in women worldwide1. High-grade serous ovarian cancer (HG-SOC) constitutes ∼70% of all ovarian malignancies and has the worst prognosis2. Current treatment of HG-SOC consists of cytoreductive surgery and combination chemotherapy with platinum-containing DNA-crosslinking drugs and taxane-based microtubule-stabilizing agents2. Although treatment significantly improves survival, most women with HG-SOC relapse with chemotherapy-refractory disease and eventually succumb3. Multiple mechanisms of chemo-resistance have been documented4,5, including reduced intracellular drug accumulation6, detoxification by increased levels of glutathione7, altered DNA damage repair8,9, dysfunctional apoptotic pathways10,11, and hyper-activation of various cell signaling pathways12–14. These mechanistic studies are consistent with recent genomic analyses that reveal marked clonal evolution of HG-SOC during therapy15. Other evidence, however, supports a hierarchical organization of HG-SOC, featuring intrinsically chemo-resistant “cancer stem cells” (CSC) that can escape initial treatment and seed recurrence16–18.
Although there is uniform agreement that HG-SOC patients should receive surgery and chemotherapy, there is disagreement over the optimal order and timing of these modalities. Two main options exist: primary debulking surgery with adjuvant chemotherapy (PDS), or neo-adjuvant chemotherapy, followed by interval debulking surgery (NACT)19–24. In either case, the surgical standard of care is to seek maximal cytoreduction, with the objective being to leave no visible residual disease. However, the precise definition of such “optimal debulking” varies among different centers, surgeons, and reports19,21,24,25.
Several studies have reported similar outcomes after PDS or NACT, including two highly influential randomized trials (EORTC and CHORUS) carried out across multiple countries22,23,26–28. In both trials, however, there was a bias in patient recruitment, favoring those with more extensive disease, who are less likely benefit from “upfront” surgery23,28. Consistent with this interpretation, overall survival in these trials was significantly shorter than that seen in other HG-SOC cohorts19,24,29,30. Closer examination of these studies revealed additional factors that might have influenced their conclusions. The EORTC study had inconsistencies in optimal debulking rates between participating centers, with the PDS-associated complete debulking data highly influenced by the results from a single institute23. The CHORUS study involved 76 clinical sites, and there were substantial differences in surgery execution and chemotherapy drug selection/dosage between those institutions28.
A review of cases at our center (Princess Margaret Cancer Center) revealed that PDS patients with no visible disease post-resection survived substantially longer (7-year survival >60%) than those receiving NACT (7-year survival ∼10%). Furthermore, although residual tumor post-resection is a critical determinant of survival, its influence on the PDS group was far more dramatic than on NACT group24. Of course, this report suffers from deficiencies common to all retrospective analyses, including lack of randomization to account for tumor burden at diagnosis and other factors; indeed, the NACT group in this study did have more extensive disease.
Another controversy in HG-SOC management focuses on the potential benefit of earlier diagnosis. Earlier diagnosis of primary HG-SOC is generally believed to enhance patient survival and quality of life; certainly, detecting disease confined to the ovary offers the best chance of cure3. Intuitively, one might predict that the same reasoning would apply to recurrent disease; however, survival reportedly is similar in relapsed patients treated earlier, based on increasing serum CA125 levels, versus those treated only when physical symptoms of recurrence appear31. Conceivably, the lead time between CA125 rise and clinical recurrence is too short for earlier chemotherapy to be beneficial; if so, then patient survival might be extended by more sensitive methods, such as testing for circulating tumor DNA (ctDNA)32,33.
To address these issues, we developed a mathematical framework that models the dynamics of HG-SOC progression, response to surgery and chemotherapy, and recurrence. Our results, generated over a wide range of parameters and accounting for hierarchical and stochastic models of tumor initiation, argue that PDS is inherently superior to NACT and suggest that with currently available therapies, the benefits of earlier detection are intrinsically restricted to primary HG-SOC.
RESULTS
Clinical Cohorts
To train and evaluate our model, we used two independent cohorts, comprising a total of 285 FIGO Stage IIIC-IV HG-SOC patients (Table 1 and Figure S1). The first dataset contains information on 148 patients who were treated by PDS or NACT at University Health Network (UHN) between March 2003 and November 201124. In this cohort, “complete (microscopic) debulking” (<1mm residual tumor) was achieved in 97 patients, whereas 51 patients had 1-10mm of residual tumor. Among the 97 patients with complete debulking, 40 received upfront surgery, followed by 6 cycles of platinum and taxane-based chemotherapy (hereafter, “PDS <1mm”), whereas the other 57 patients were treated with upfront chemotherapy, followed by interval surgery after the first 3-4 cycles (hereafter, “NACT <1mm”). For each patient, descriptions of tumor size post-surgery and serum CA125 levels over the course of treatment were recorded (Figure 1a). As reported previously24, overall survival was much better for patients in the PDS <1mm group (Figure 1b). However, patients in the NACT <1mm group tended to have greater tumor burden at diagnosis24, confounding direct comparison of the survival curves. Data from these patients (termed “Training Set” below) were used to train the parameters in our computational model.
The other 51 patients in the UHN dataset underwent debulking, but were left with residual tumor of 1-10mm in diameter. Incomplete debulking usually was due to tumor location. As reported previously24, there was a significant negative association between residual tumor size and patient survival in patients who had undergone PDS (Figure 1c), but no association was found in the NACT group (Figure 1d). Similar trends have been observed in several other studies investigating the association between patient survival and residual tumor size19,20,34. Independent of the <1mm cohort, we used this 1-10mm cohort (“Validation Set 1”, below) in initial tests of the validity of our computational model (Figure S1).
The second dataset, which is independent of the UHN cohort, contains information on 137 patients diagnosed from October 2001 to June 2005 and enrolled in the CAN-NCIC-OV16 trial (hereafter, “NCIC cohort”). These patients all received upfront surgery, followed by 8 cycles of carboplatin/paclitaxel chemotherapy. The median delay for chemotherapy after surgery was 0.8 month (compared with 1 month at UHN). In 20 patients, debulking to <1mm residual tumor was achieved, 36 patients had 1-10mm residual tumor, and 81 patients had residual tumor diameter >10mm. Data from the NCIC cohort was used to further evaluate the validity of our computational model (Figure S1).
Model Overview
We developed a mathematical model of cancer initiation and evolution to investigate the dynamics of HG-SOC growth, the onset of chemo-resistance, and the effects of various treatment strategies on patient survival (Figure 1e). Our initial framework was based on work by the Michor laboratory35,36, and considers exponential expansion of HG-SOC cells starting from a single cancer cell that has all of the genetic alterations needed for proliferation and metastasis, but has not developed chemo-resistance. This framework assumes that during tumor development/progression, any chemo-sensitive HG-SOC cell can acquire mutations and/or epigenetic alterations that enable chemo-resistance; i.e., we assume a stochastic model of tumor initiation. However, the model can be modified easily to accommodate a hierarchical organization of HG-SOC, in which a relatively chemo-resistant CSC gives rise to chemo-sensitive progeny. Importantly, the predictions/implications of the hierarchical and stochastic models are essentially the same (see below and Supplemental Information).
In the absence of chemotherapy, chemo-sensitive cells divide at rate r and die at rate d, genetic or epigenetic changes that cause chemo-resistance occur with probability u per cell division, and chemo-resistant cells divide at rate a and die at rate b (which can be the same or different as the rates of division and death of chemo-sensitive HG-SOC cells). The total number of HG-SOC cells at diagnosis is denoted by M1, with a proportion εM1 metastasized outside the peritoneal cavity and not accessible at debulking surgery. Once cancer has been diagnosed, the (virtual) patient receives PDS or NACT. Post-surgery, the number of HG-SOC cells in the peritoneal cavity is reduced to M; cancer cells at inaccessible locations persist. During chemotherapy, chemo-sensitive cells divide at rate r’ and die at rate d’, with r’ < d’; the probability of conversion from chemo-sensitive to chemo-resistant is u’ per division; and chemo-resistant cells divide at rate a’ and die at rate b’. After first line treatment, persisting sensitive and resistant cells continue to proliferate and reach number M1, at which time relapse is diagnosed, and the patient is again treated with chemotherapy. Note that although tumor recurrence is usually diagnosed before it reaches M1 in the clinic, our analyses (see below and Figure 6), along with clinical evidence31, suggest that the timing of diagnosis of recurrence does not affect treatment outcome. This regression-relapse cycle is repeated until chemotherapy fails to reduce disease burden, at which point tumor growth continues until the patient’s death, with a total number of cancer cells denoted by M2, which is larger than M1.
This framework enables us to estimate the probability that chemo-resistance is present at diagnosis, the dynamics of sensitive and resistant cell numbers along the treatment course, and the length of survival of HG-SOC patients. Analytic approximations for these quantities are shown in the Experimental Procedures section. With these quantities, we estimated the conversion rate (u) and inaccessible proportion (ε) by minimizing the deviations between patient survival data and the corresponding predictions obtained using our model. Using these estimates, we then predicted the survival of patients in test datasets to assess the validity of our model, compared the survival of PDS- or NACT-treated patients when controlled for starting tumor burden, and probed the potential benefits of earlier diagnosis of HG-SOC.
Estimation of Parameter Values
For most model parameters, a clinically relevant range of values could be deduced from clinical data or previous publications (see Extended Experimental Procedures for detailed description of parameter value estimations). We then varied u and ε over a wide range, and computed the expected distribution of survival of patients by Monte-Carlo simulation. We compared the deviation between the Training Set data and the predictions of our model for each combination of u and ε, and from the region of best fit for these parameters, we inferred that 10−9 < u < 10−7 and 10−9 < ε < 10−5 (Figure 2a). Finer investigation of the fit between data and theory in this parameter region showed that u = 10−7.6 and ε = 10−7.4 is the best combination of these parameters, as it minimizes the deviation between data and theory (Figure 2b). Importantly, the trained value for u is consistent with chemo-resistant conversion rates estimated for multiple other cancers37–40. The trained value for ε predicts that inaccessible cancer cells outside the peritoneum are less numerous than cancer cells inside the peritoneum after first-line therapy, which comports with clinical observations that recurrent HG-SOC occurs predominantly inside the abdominal cavity. Using the trained parameter values, we compared the observed (Figure 1b) and predicted (Figure 2c) distributions of patient survival in the Training Set. Reassuringly, the model predictions closely recapitulated the clinical observations (Figures 2d-e).
Using our mathematical framework and these estimated rates, we then calculated the probability (as a function of tumor burden) that chemo-resistant cells are present at diagnosis (Figure S2a). Consistent with clinical and theoretical studies of other malignancies37–39,41–45, our model predicted that at least some chemo-resistant cells are always present at the start of therapy for HG-SOC. We varied the values for each model parameter over a large range to test their influence on the calculated number of resistant cells at diagnosis: within the ranges tested, this number almost always exceeds 103 (Figures S2b-g).
Model Validation
To assess the accuracy of our mathematical framework and the validity of the estimated parameter values, we analyzed the model-predicted distribution of overall survival using data from several test sets (Figure 3). First, we predicted the survival of patients in Validation Set 1 (Figures 3a-b, red curves), by adjusting the value of M. Notably, our model recapitulated the clinical observation that residual tumor size is significantly associated with survival in patients undergoing PDS, but not NACT (Figures 3a-b). Comparison of the observed (Figures 1c-d, red curves) and predicted (Figures 3a-b, red curves) distributions of patient survival revealed no significant differences (Figures 3c-d).
The data for the above validation exercise derived from different patients than those used for the Training Set, but all of these patients were treated over the same time period at the same institution (UHN). To more rigorously test the validity of our model, we analyzed independent patient data, derived from the NCIC cohort. As noted above, the treatment course of these patients also differed somewhat from that of the UHN patients, enabling an even better test of the general applicability of our framework. We incorporated these differences into the model and predicted the survival of patients in the NCIC trial. Again, the model predictions fit very well with clinical observations, and recapitulated the theme that PDS with minimal residual tumor results in the best outcome (Figures 3e-i).
Predicted Outcome of PDS and NACT in Patients with Identical Tumor Burden
Confident in the predictive ability of our mathematical framework, we modeled the expected clinical outcome of PDS and NACT in patients with the same initial tumor burden by imposing the same distribution of M1 on both groups. Our model predicts that PDS patients should survive longer than NACT patients when controlled for residual post-surgery tumor mass. As residual tumor increases, however, the predicted survival advantage of PDS shrinks (Figures 4a-b). To investigate the underlying reason(s) for these predictions, we explored the predicted dynamics of chemo-sensitive (green) and chemo-resistant (blue) HG-SOC cells in patients undergoing treatment by PDS or NACT (Figures 4c-d). Our model predicts that, at diagnosis, a typical HG-SOC patient has low numbers of chemo-resistant cells. In women who undergo PDS, debulking surgery (S in Figure 4c) dramatically reduces the number of chemo-sensitive and chemo-resistant cancer cells, because these cells appear identical to the surgeon and therefore have an equal likelihood of removal. Depending on the (stochastic) distribution of chemo-resistant cells within the abdominal cavity of the HG-SOC patient, all chemo-resistant cells might have been eliminated if, post-surgery, there is <1mm residual tumor. Follow-up chemotherapy (C in Figure 4c) can then reduce the remaining chemo-sensitive cells to very low numbers or even eradicate them.
By contrast, with NACT, neo-adjuvant chemotherapy (C) dramatically enriches for chemo-resistant cells, while killing the sensitive cells; consequently, chemo-resistant cells comprise a large proportion of total tumor cells at surgery (S). Because chemo-sensitive cells are largely depleted by the neo-adjuvant chemotherapy, the amount of residual tumor visible to the surgeon is reduced substantially. Consequently, it is virtually impossible for interval debulking surgery to fully deplete the chemo-resistant cells (Figure 4d). We propose that the relative inability of NACT to deplete chemo-resistant cells underlies the difference in outcome between patients treated with PDS and NACT. Importantly, this conclusion is based on intrinsic properties of the dynamics of cancer proliferation, survival, and death.
We then explored why the predicted superiority of PDS over NACT depends on residual tumor burden post-surgery. By examining the expected distribution of chemo-sensitive and -resistant cell numbers after first-line therapy, we found that PDS with <1mm residual tumor has the potential to deplete all cancer cells in a significant proportion of HG-SOC patients (Figure 4e). This finding can account for the considerable survival difference between PDS and NACT with <1mm residual tumor. By contrast, with >1mm residual tumor, neither PDS nor NACT depletes all malignant cells, even though fewer tumor cells are predicted to remain after PDS (Figure 4f). Consequently, almost all patients with >1mm residual tumor are predicted to relapse and eventually die because of the inability of current agents to kill chemo-resistant cells. This analysis explains why the survival advantage conferred by PDS over NACT shrinks as residual tumor increases and why residual tumor mass is a key determinant of survival after PDS but not NACT (Figures 3a-b).
Predicted Outcome of PDS and NACT in Patients with Identical Tumor Burden
To test the robustness of our predictions and quantitatively examine the influence of each factor on survival, we varied the values of every parameter in our model over a large range. Our conclusions on the relative efficacy of PDS and NACT treatment hold true throughout, although the magnitude of the predicted difference can vary (Figure S3). Systematic modulation of the model parameters also enabled us to assess the contribution of each factor to patient survival, and revealed several interesting features. For example, we found that faster growth of cancer cells does not necessarily lead to worse survival. In particular, although elevated growth rate of chemo-resistant cells in the absence of chemotherapy leads to worse survival, faster-growing, chemo-sensitive cells can sometimes lead to better survival (Figures S3a-d). This result can be attributed to a lower percentage of chemo-resistant cells at diagnosis when chemo-sensitive cells proliferate faster. However, if the proliferation rate of chemo-sensitive cells is too high (e.g., when r = 3), new chemo-resistant cells might be generated between the time of surgery and adjuvant chemotherapy, resulting in worse survival (Figure S3c). We also inferred that the relative importance of chemo-sensitive and chemo-resistant cells in influencing survival might differ at different stages along the clinical course. During treatment-free periods, the growth rate of chemo-resistant cells could play a more dominant role in influencing patient survival (Figures S3a-d), because they underlie ultimate treatment failure. However, during periods of chemotherapy, the growth rate (or depletion rate) of chemo-sensitive cells might play a more dominant role (Figures S3e-h), because the depletion rate of chemo-sensitive cells at this stage determines whether they can be completely eradicated by chemotherapy. By contrast, chemo-resistant cells likely will endure. As a result, drug choice and dose, which likely influences the efficiency of elimination of chemo-sensitive cells, is predicted to be a critical factor in treatment outcome. Consistent with this analysis, residual cancer cell abundance after tumor resection can dramatically influence patient survival, especially in the PDS group (Figure S3i), suggesting that primary cytoreductive surgery should aim for complete removal of cancer cells even though chemotherapy will usually follow. Tumor size at diagnosis can be a critical factor determining whether a patient can be cured (Figure S3j), and will be further addressed below. By contrast, varying the parameter value of tumor size at patient death did not influence the length of survival (Figure S3k).
We also considered two alternative scenarios: (1) that heterogeneous populations of (variably) chemo-resistant cells exist in the same patient (Figures S4a-c); or (2) that HG-SOC initiates from an intrinsically chemo-resistant “cancer stem cell”, which differentiates into chemo-sensitive “tumor progenitor cells” (Figures S4d-f). Either of these scenarios results in the same conclusions about the relative efficacy of PDS and NACT and its dependence on residual tumor post-surgery.
Modeling Alternative Treatment Regimens
We next utilized our model to predict the effects of altering current treatment regimens. For both PDS and NACT, adjuvant chemotherapy typically begins 4-5 weeks post-surgery, an interval chosen to allow patients to recover from their typically aggressive surgical treatment. The length of the post-surgical chemotherapy delay varies between centers and among physicians, but its influence on treatment outcome has not been studied carefully. We varied the length of treatment delay in our model, and tested the predicted effects on patient survival. For PDS with <1mm residual tumor, earlier initiation of chemotherapy should prolong survival, whereas longer treatment delay worsens outcome (Figure 5a, p<0.0001). For PDS with >1mm residual tumor or for NACT with any amount of residual tumor, the length of treatment delay (within the same range) is predicted to have little effect on outcome (Figures 5b-d). These differences arise primarily because PDS that results in <1mm residual tumor potentially can deplete all tumor cells when treatment delay is minimized (Figure 4e), and the probability of depletion decreases with longer delay, largely because chemo-resistant cells can arise during this period (Figure S6a). By contrast, PDS with >1mm residual tumor or NACT is unlikely to deplete all cancer cells (Figures 4e-f and S6b-d), irrespective of treatment delay.
For the NACT regimen, our model predicts that neo-adjuvant chemotherapy enriches for chemo-resistant cells and thus increases the percentage of chemo-resistant tumor cells remaining post-surgery. Conceivably, reducing the number of cycles of neo-adjuvant chemotherapy could attenuate the enrichment for chemo-resistant cells, thereby prolonging patient survival. Indeed, our model predicts that patients undergoing NACT with <1mm or >1mm residual tumor could benefit slightly from reducing the number of pre-surgery chemotherapy cycles (Figures 5e-f). This small predicted improvement arises primarily because of more efficient reduction of the number of chemo-resistant cells at surgery. The benefit is limited, however, because for most patients, chemo-resistant cells are still unlikely to be eradicated after their enrichment during the neo-adjuvant chemotherapy period (Figures S5e-f).
Predicted Effects of Earlier Diagnosis on Survival
We also utilized our model to evaluate the potential benefits of earlier diagnosis. We first modeled the effects of diagnosing relapsed HG-SOC at the earliest possible time enabled by the currently used clinical test, CA-125 detection46. Contrary to the intuitive notion that earlier diagnosis of recurrence should be advantageous, we find that CA-125-based earlier diagnosis is not expected to improve survival (Figures 6a-d, red curves). We also asked if detecting recurrence earlier than is possible with CA-125 monitoring would be advantageous. Yet even with a lead time ∼104 greater than that required for physical symptoms to appear, the current upper limit of sensitivity for ctDNA-based diagnosis of ovarian cancer33, our model predicts no advantage in patient survival with earlier detection of recurrence (Figures 6a-d, green curves).
We explored the reason for this lack of survival advantage by modeling the theoretical numbers of chemo-sensitive and –resistant cells in a virtual patient whose recurrence is diagnosed with increasing levels of sensitivity (Figure 6e). Although earlier diagnosis, followed by prompt re-institution of chemotherapy, can better deplete chemo-sensitive cells, it barely affects chemo-resistant cells, which ultimately expand and cause patient death. Therefore, earlier diagnosis is unlikely to improve survival when applied to relapsed cancers treated with standard cytotoxic chemotherapy. We further interrogated the quantitative relationship between the length of lead time and the potential advantage in survival, and found that to achieve an advantage of 1 month in overall survival, a lead time of over 9 months is probably necessary, and even then, is not sufficient (Methods).
Finally, we used our model to explore the potential benefit of earlier diagnosis of treatment-naïve tumors. Because earlier upfront diagnosis would identify smaller, probably less disseminated tumors, it likely would increase the chance for complete debulking during first-line therapy. We therefore focused our comparison on the predicted effects in patients with <1mm residual tumor. Our analysis argues that earlier diagnosis of treatment-naïve cancer, with concomitant prompt intervention, can improve patient survival whether they undergo PDS or NACT (Figure 7). For PDS, the predicted survival benefit can be dramatic (Figures 7a, 7c, and 7e), primarily because lower volume and less diffuse tumor at presentation can increase the likelihood of disease eradication (Figures S5a, S5c, and S5e). Earlier detection can also provide some benefit for patients undergoing NACT if lead time is sufficient (Figures 7b, 7d, and 7f). In that case, chemo-resistant cells might not yet have arisen, and chemotherapy alone might be sufficient to eradicate disease (Figures S5b, S5d, and S5f).
DISCUSSION
Mathematical modeling has demonstrated potential in the systematic and quantitative assessment of various treatments44,47,48. If the outcomes of different strategies could be modeled accurately a priori, clinical trials could focus on therapeutic combinations that are most likely to succeed, and improvements in patient outcomes could be accelerated. Predictive models also can suggest when specific clinical controversies merit re-examination in the controlled trial setting. Our mathematical model of HG-SOC defines factors that affect the evolution of chemotherapy resistance and can predict patient survival based on the growth dynamics of chemo-sensitive and –resistant cells, whether stochastic or hierarchical models of tumor initiation are assumed. Furthermore, the results of our analyses have important implications for HG-SOC therapy and screening.
We populated our model with clinical data from ∼300 patients receiving PDS or NACT. After estimating the rates of tumor cell proliferation and conversion to chemo-resistance, we concluded that most HG-SOC patients probably harbor chemo-resistant cancer cells at diagnosis. Our outcome predictions closely match patient data from multiple sources and support clinical observations that: (1) PDS that leaves minimal residual tumor is the optimal treatment strategy; (2) residual tumor size is a critical determinant of survival in patients undergoing PDS, but not NACT; (3) earlier diagnosis of relapsed cancer does not – and cannot – lead to better survival with current therapies; and (4) earlier diagnosis of primary (treatment-naïve) HG-SOC could dramatically improve outcomes from this devastating disease, especially if much more sensitive detection methods can be developed.
Our model shows clearly that, when controlled for residual cancer size after surgery, PDS should always lead to a better outcome than NACT. If debulking leaves minimal (<1mm) residual tumor, the elemental properties of tumor dynamics predict that PDS should result in dramatically longer overall survival than NACT. We infer that the reason for the superior predicted outcome of PDS is that upfront surgery can deplete minority chemo-resistant cells more effectively, leaving adjuvant chemotherapy to eradicate residual chemo-sensitive cells. If debulking surgery removes all chemo-resistant cells, the PDS regimen can be curative. By contrast, cure is unlikely after NACT for at least two reasons. First, neo-adjuvant chemotherapy enriches for, and allows the continued expansion of, chemo-resistant cells. Second, by depleting bulk chemo-sensitive cells, chemotherapy removes tumor mass that can mark the location of chemo-resistant cells interspersed in metastatic deposits throughout the peritoneum. We suggest that removing these “sentinel” chemo-sensitive cells renders interval surgery less effective than upfront surgery in depleting chemo-resistant cells, which are the cells that eventually cause death. Put another way, with existing agents, the only way to achieve cure, and the best way to promote long-term survival, is to remove all, or the vast majority of, chemo-resistant cells. At present, the only way to achieve this goal is by surgically removing such cells, which is much more likely with the PDS regime. Of course, the development of drugs that can target these resistant cells might alter these conclusions17,22,49.
Previous clinical studies differed over whether PDS or NACT is superior. A retrospective analysis of patients treated between 1980 and 1997 found no difference in outcome between PDS and NACT26. By contrast, several meta-analyses indicated that NACT is associated with worse prognosis21,50. This controversy occasioned two large, controlled clinical trials to compare PDS and NACT23,28, which found no significant difference in patient survival between PDS and NACT. Three major weaknesses limit the strength of these trials’ conclusions: (1) They only recruited patients with “extensive” HG-SOC, which might explain the notably worse outcome of patients in this study, especially in PDS <1mm group, compared with multiple other reports19,20,24,51. Our computational model predicts that elevated upfront tumor burden can lead to worse overall survival and obscure the difference in outcome between PDS and NACT, even when complete debulking is performed (Figure S3j). This factor alone might contribute to the lack of difference between PDS and NACT in these trials. (2) In the EORTC trial, there was marked inconsistency between the participating institutions in percentage of optimal debulking in the PDS group, with 63% for one center in a single country and <12% for the other 6 countries involved23,52. Our model predicts that residual tumor size is a critical factor for patient survival, especially in the PDS group, and increased abundance of residual cancer cells can decrease the advantage of PDS (Figure S3i). (3) In the CHORUS trial, many patients (30-40%) received carboplatin monotherapy, instead of combination carboplatin and paclitaxel28. Previous studies demonstrated that platinum-based combination therapy in ovarian cancer leads to better survival than platinum therapy alone53, likely by more effectively depleting chemo-sensitive cells or by killing some platinum-resistant cells. Our model predicts that the efficiency of depletion of chemo-sensitive cells during chemotherapy is a critical factor for patient survival (Figure S3e-f). Therefore, the model’s predictions might explain the different findings from existing clinical studies.
A previous mathematical modeling study, based on a Gompertzian growth model, argued that NACT should be superior to PDS54. However, that study did not account for the different dynamics of chemo-sensitive and –resistant cells, and assumed equal efficiency of surgical depletion of large and small tumors, which ignores the intrinsic limitations of surgery. Such assumptions can lead to serious errors in modeling HG-SOC. For example, although aggressive surgery might reduce a tumor containing >1011 cells to a mass of <107 cells, a tumor containing <106 cells is unlikely to be visible during surgery, making it highly unlikely that such a tumor can be surgically reduced to <102 cells. Notably, the conclusions of our study still hold if our model is implemented with Gompertzian growth assumptions (data not shown).
Similarly, a study of the optimal order of surgery and chemotherapy in pancreatic cancer concluded that neo-adjuvant chemotherapy should be superior to upfront surgery36. However, the biology and chemo-responsiveness of pancreatic cancer and HG-SOC differ substantially: whereas HG-SOC is usually quite chemo-responsive, pancreatic cancer is notoriously chemo-resistant. The differential influence on tumor visibility following chemotherapy likely underlies the different predictions in the two diseases. Nevertheless, the difference between our conclusions and those of Haeno et al. argue for caution in extrapolating their conclusions to other types of cancer.
Our model enables quantitative analysis of the dependence of treatment outcome on various factors under different scenarios, providing a powerful tool to assist clinical decision-making. For example, we found that the advantage of PDS over NACT diminishes with larger residual tumor post-surgery, more extensive metastases at unresectable locations, and/or a higher percentage of chemo-resistant cells at diagnosis (Figures 3a-b, and S3a-c). These features might help explain the similar overall survival between PDS and NACT groups in studies involving patients with more extensive upfront disease or less complete surgical removal23,26. Our analysis also argues that for PDS patients with <1mm residual tumor, adjuvant chemotherapy should start as early as possible to provide the best chance of curative outcome; indeed, we predict that differences of even a few weeks might dramatically alter the chance for curative outcome. Conversely, if complete debulking is not achieved at primary surgery, delaying chemotherapy is less likely to affect survival. These predictions are consistent with the results of a meta-analysis of clinical studies55, which found that each extra week of delay was associated with a significantly decreased overall survival in patients with PDS <1mm residual tumor, but not in patients with visible residual disease. We also predict that reducing the number of cycles of neo-adjuvant chemotherapy should result in slightly improved overall survival for NACT patients, a prediction that also is consistent with a meta-analysis of clinical data50. Taken together, these findings argue that even in patients too infirm to undergo immediate surgical debulking, the number of neo-adjuvant cycles before surgery should be minimized.
A final prediction of our model is that the effects of earlier diagnosis of HG-SOC should be quite different in relapsed versus treatment-naïve patients. For the former, earlier diagnosis is unlikely to improve overall survival, primarily because earlier treatment of recurrent tumors with existing agents cannot deplete chemo-resistant cells that have been enriched over the clinical course. This prediction matches very well with earlier clinical observations31. Two scenarios might alter this conclusion: (1) if effective chemotherapy with a resistance profile orthogonal to platinum/taxane-based therapy were to become available at relapse; or (2) if effective debulking could be achieved at relapse. Effective drugs against HG-SOC remain a major clinical limitation56,57. Our findings argue strongly for employing alternative agents as early as possible upon tumor relapse. Surgery is no longer typically performed on recurrent HG-SOC, and its efficacy could be limited by the proportion of cancer cells at unresectable locations at disease relapse. Nevertheless, our model calls for re-evaluation of feasibility and potential efficacy of secondary surgery for the benefit of earlier diagnosis. Consistent with our analysis, several retrospective studies report that secondary surgery achieving complete debulking can be beneficial for HG-SOC patients with platinum-sensitive tumors58,59.
By contrast, our model predicts that earlier diagnosis of treatment-naïve cancer could improve overall survival for at least two reasons: (1) chemo-resistant cells are usually not enriched prior to treatment and earlier intervention can reduce the likelihood that significant numbers will arise; and (2) earlier upfront surgery has a better chance of removing all chemo-resistant cells, assuming that tumor cells have not diffused throughout the peritoneal cavity or seeded unresectable locations. There also are at least two additional benefits of earlier diagnosis and treatment of naïve HG-SOC: (1) patients receiving incomplete debulking with regular diagnosis might achieve complete debulking with earlier diagnosis; and (2) patients considered treatable only by NACT with regular diagnosis might be eligible for PDS with earlier diagnosis.
In summary, our analyses suggest that future randomized clinical trials might consider: (1) the influence of the interval between primary debulking surgery and adjuvant chemotherapy on treatment outcome; (2) the association between the number of neo-adjuvant chemotherapy cycles and treatment outcome; and (3) the effects of alternative chemotherapy and/or complete secondary surgery on relapsed tumor, especially when coupled to earlier diagnosis. Although political and logistical reasons might prevent a new controlled trial of PDS versus NACT, we would argue strongly for such a trial if significant numbers of gynecologic oncologists persist in arguing for the equivalence of these two modalities. Finally, the mathematical abstraction makes our framework potentially applicable to treatment of other malignancies with alternative parameter values.
METHODS
Patient Information
Clinical data from 285 HG-SOC patients were obtained from patient records at University Health Network (148 patients) and from the National Cancer Institute Canada OV16 trial (137 patients). Medical record information, including date of diagnosis, age of patient, timing of treatment, extent of residual disease in diameter (<1mm, 1-10mm, or >10mm), CA125 levels along the treatment course, and survival data (updated in 2014 for UHN data and 2010 for NCIC data) were recorded. Institutional REB approval was obtained through University Health Network and a Data Sharing Agreement was concluded with National Cancer Institute Canada.
Variables Used in Mathematical Model
We denote the following variables for the development of mathematical modeling: r, division rate of chemo-sensitive cells in the absence of chemotherapy; d, death rate of chemo-sensitive cells in the absence of chemotherapy; a, division rate of chemo-resistant cells in the absence of chemotherapy; b, death rate of chemo-resistant cells in the absence of chemotherapy; r’, division rate of chemo-sensitive cells in the presence of chemotherapy; d’, death rate of chemo-sensitive cells in the presence of chemotherapy; a’, division rate of chemo-resistant cells in the presence of chemotherapy; b’, death rate of chemo-resistant cells in the presence of chemotherapy; M, the number of residual cancer cells in the peritoneal cavity after surgery; M1, the total number of cancer cells at diagnosis; M2, the total number of cancer cells at death.
Estimation of Parameter Values
Unless otherwise specified, clinically relevant parameter values were estimated and set as follows:
In the absence of chemotherapy, the proliferation rate of chemo-sensitive cells (r) was obtained from the normal distribution with mean 2 and standard deviation 0.8, based on a clinical study interrogating the doubling time of ovarian cancer cells by BrdU labeling60. Another study compared cancer cell proliferation rates in platinum responders and non-responders by thymidine labeling, and found that responders had significantly higher labeling index than patients with stable or progressive disease61. Accordingly, we set the proliferation rate of chemo-resistant cells (a) to a normal distribution with mean 0.84 and standard deviation 0.42. Death rates (d and b) were set as 10% of proliferation rates (r and a) respectively, which is within physiologically relevant range; in any case, varying this ratio over a large range does not affect the main conclusions of this paper.
During chemotherapy, the proliferation rate of chemo-sensitive cells (r’) was set to 10% of that of randomly growing cells (r); the death rate of chemo-sensitive cells (d’) was set to a normal distribution with mean 4.9 and standard deviation 1. These values were set to match the clinical observation that chemo-responsive relapse occurs on average ∼11 months after 6 cycles of chemotherapy, indicating that cytotoxic reduction of chemo-sensitive cells by 6 cycles of chemotherapy corresponds to ∼11 months of proliferation. Our inference is consistent with published efficacy of chemotherapy in HG-SOC25,62. The proliferation rate of chemo-resistant cells during chemotherapy (a’) was set to the same as that during random growth (a); the death rate of chemo-resistant cells during chemotherapy (b’) was set to twice as that during randomly growing (b), or 20% of a’, to reflect a modest effect of chemotherapy on chemo-resistant cells. Conversion rate during chemotherapy (u’) was set as 10 times of that during randomly growing (u), to reflect the DNA-damaging effect of platinum-based chemotherapy.
The parameters reflecting cancer cell numbers (M, M1, M2) were obtained from normal distributions in base 10 logarithmic scale. The number of residual cancer cells immediately post-surgery (M) was set to mean 6 and standard deviation 0.4 on the log scale for <1mm residual tumor, mean of 8 on the log scale for 1-10mm residual cancer, and mean of 10 on the log scale for >10mm residual cancer. The number of cancer cells at diagnosis (M1) was initially set at different values in patients receiving PDS and NACT, to reflect the clinical observation that patients with NACT tend to have more extensive disease at diagnosis24. M1 for PDS was set at mean 11.5 and standard deviation 0.4 on the log scale, and M1 for NACT was set at mean 12 and standard deviation 0.4; these values result in NACT patients starting with >3x the cancer burden in the model than those receiving PDS. Moreover, we varied the ratio of tumor burden at diagnosis in NACT vs PDS patients from 1 to 101.5, and the main conclusions hold over the entire range. M2 was set to mean 13 and standard deviation 0.4 on the log scale.
We varied the values for all of the above parameters to test the robustness of our conclusions. We found that the main conclusions hold true, though the magnitude of the differences between PDS and NACT treatment outcomes may vary.
Mathematical Deduction of the Expected Number of Chemo-resistant Cells at Diagnosis
We deduced the expected the number of chemo-resistant cells at diagnosis based on previous studies35,36. We first calculated the probability that chemo-resistant cells exist at diagnosis (Pd), and then calculated the expected number of chemo-resistant cells at diagnosis (Yd).
To calculate Pd, we summed the probabilities that the first successful lineage of chemo-resistant cells arises when there are 1, 2, 3, … M1-1 chemo-sensitive cells. P(x) denotes the probability that the first lineage arises when there are x chemo-sensitive cells. P(x) can be expressed as the joint probability that no successful chemo-resistant lineage arises at 1, 2, 3, … x-1 chemo-sensitive cells, and that a successful lineage arises at exactly x chemo-sensitive cells. An expected 1/(1 − d/r) divisions are needed for an effective increase of 1 chemo-sensitive cell, and during these divisions an expected u/(1 − d/r) chemo-resistant cells are generated, among which a proportion (1 − b/a) will successfully persist. Assuming that the number of surviving chemo-resistant cells generated by each division follows a Poisson distribution with mean (1 − b/a)u/(1 − d/r)45, we derive P(x):
Thus the probability that chemo-resistant cells exist at diagnosis, Pd, can be written as the sum of P(x):
If we denote the time between the emergence of a successful chemo-resistant cell and diagnosis as τx, then the expected number of chemo-resistant cells at diagnosis in the patients who have them can be expressed as:
The amount of time τx satisfies:
Computational Estimation of Patient Survival
We applied the following workflow to simulate clinical treatment of HG-SOC. Overall survival is monitored and recorded along the workflow:
Assessment of Prediction Robustness Regarding PDS vs. NACT Regimen
To assess the robustness of the model predictions of outcome following PDS vs. NACT, we varied u, ε, and each of other parameters, and simulated disease progression of mock PDS and NACT patients (N=50 each). Survival curves for PDS and NACT regimens within each combination of parameter values were compared based on the hazard ratio of the two groups, which was calculated by fitting a Cox proportional hazards regression model63 to the simulated survival data using the R survival package.
Alternative Scenario 1 – Heterogeneous Chemo-resistant Cells in a Patient
We initially designed a two-component model where a single type of chemo-resistant cell exists in a patient. We then extended our mathematical framework to include a scenario in which heterogeneous chemo-resistant cells co-exist in a patient. Within this framework (Figure S3d), ovarian cancer starts from a chemo-sensitive cell with division rate r and death rate d. At each division, a chemo-sensitive cell has probability u1 to convert into a Type-1 chemo-resistant cell, which divides at rate a1 and dies at rate b1, and probability u2 to convert into a Type-2 chemo-resistant cell, which divides at rate a2 and dies at rate b2. During chemotherapy, chemo-sensitive cells divide at rate r’ and die at rate d’, Type-1 chemo-resistant cells divide at rate a1’ and die at rate b1’, Type-2 chemo-resistant cells divide at rate a2’ and die at rate b2’, and conversion rates are u1’ and u2’, respectively. We assume that u1 and u2 are small enough so that conversions into Type-1 and Type-2 cells do not take place in the same division, and that the majority of cancer cells at diagnosis are chemo-sensitive.
We could then approximate the expected the numbers of Type-1 (Y1d) and Type-2 (Y2d) chemo-resistant cells at diagnosis:
In the above equations, τ1x and τ2x represent the time from the emergence of a successful Type-1 or Type-2 chemo-resistant cell to diagnosis of disease, respectively, and approximately satisfy:
P1d and P2d represent the probabilities that Type-1 or Type-2 chemo-resistant cells exist at diagnosis, respectively, and can be expressed as:
We then approximated the dynamics of chemo-sensitive (X), Type-1 chemo-resistant (Y1), and Type-2 chemo-resistant (Y2) cells as below:
Alternative Scenario 2 – Cancer Stem Cell Model
We also considered the alternative scenario in which HG-SOC originates from an intrinsically chemo-resistant cancer stem cell, which can differentiate into chemo-sensitive cancer progenitor cells16,17,64. Within this framework (Figure S3g), HG-SOC starts from a chemo-resistant stem cell with division rate a and death rate b per time unit. Each stem cell division has a probability of u to produce a chemo-sensitive differentiated cell, with division rate r and death rate d per time unit. During chemotherapy, chemo-resistant stem cells have division rate a’ and death rate b’, and chemo-sensitive differentiated cells have division rate r’ and death rate d’. Definitions of ε, M, M1, and M2 are the same as in the initial model. The dynamics of chemo-resistant stem cells (Y) and chemo-sensitive differentiated cells can be written as:
Potential Benefit of Earlier Diagnosis of Relapsed Cancer
In order to assess the potential benefit of earlier diagnosis of relapsed cancer more stringently, we mathematically deduced and compared the expected the number of chemo-sensitive and chemo-resistant cells after early or standard diagnosis. For early diagnosis, we assumed that a patient receives chemotherapy from time 0 to TC, and the cancer randomly grows from time TC to TC+TR. For standard diagnosis, we considered that the cancer randomly grows from time 0 to TR, and chemotherapy takes place from TR to TR+TC. TR here represents the “lead time” for early diagnosis. The initial numbers of chemo-sensitive and chemo-resistant cells at the time of early diagnosis are denoted as X0 and Y0, respectively.
We first modeled the effects of early diagnosis. From time 0 to TC, during chemotherapy, the dynamics of chemo-sensitive and –resistant cells can be written as:
Solving these differential equations, we have:
At time TC, the numbers of chemo-sensitive and –resistant cells are:
Then, from time TC to TC+TR, the dynamics of chemo-sensitive and –resistant cells can be written as:
Integrating over time, we derived the numbers of chemo-sensitive and –resistant cells at time TC+TR:
In the equation for the number of chemo-resistant cells, the term represents the number of chemo-resistant cells derived from exponential expansion of those already present at time 0, whereas the other terms represent the number of chemo-resistant cells generated from chemo-sensitive cells during the period of TC+TR.
Using similar methodology, we could deduce the numbers of chemo-sensitive and –resistant cells at TR+TC for standard diagnosis:
We found that Xend(early) = Xend(standard), indicating that one round of early diagnosis is unlikely to affect the number of chemo-sensitive cells. The effect on chemo-sensitive cells in Figure 6e is caused by the shorter time lapse between different lines of therapy, which is facilitated by early diagnosis.
We then calculated the difference between Yend(early) and Yend(standard):
We imputed the clinically relevant values for these parameters (same as the mean values shown in Figure 2), and interrogated the influence of TR on the relative difference in chemo-resistant cell numbers between early and late diagnosis:
Equation (1) demonstrates that the relative difference in chemo-resistant cells is associated with two factors: (1) the initial ratio of chemo-resistant/-sensitive cells at early diagnosis (Y0/ X0), and (2) the length of lead time for early diagnosis (TR). If we impute Y0/X0=0.01, meaning ∼1% cancer cells at early diagnosis of relapsed cancer are chemo-resistant, then TR needs to be longer than 11.6 months in order for the ratio in equation (1) to be less than 0.5, corresponding to an advantage of ∼1 month in overall survival for early diagnosis. If 0.1% cancer cells at early diagnosis of relapsed cancer are chemo-resistant, then TR needs to be longer than 9.4 months for an advantage of ∼1 month in overall survival.
Estimation of Conversion Rate u and Unresectable Proportion ε
For each combination of candidate values for u and ε, we calculated the deviation of model prediction from clinical observations as:
The combination that led to lowest deviation was used for model testing and further predictions.
Statistical Analyses
Comparisons of overall survival between 2 groups were performed by log-rank test. Comparisons of predicted survival between >3 groups were performed by log-rank test for trend when the order of groups is logical. Comparisons of distribution of overall survival between model predictions and clinical observations were performed by Chi-square test.
Data Availability
All data were acquired from published clinical trials or retrospective studies.
AUTHOR CONTRIBUTIONS
S.G. and B.G.N. designed the study; S.G. constructed the computational model and performed the analyses; S.L., L.H.B., I.V., P.C., M.B., B.R., and A.O. provided clinical expertise and helpful discussions; A.S. provided computational and statistical expertise and helpful discussions; L.H.B. and I.V. extracted raw clinical data from UHN dataset; D.T. and W.P. extracted raw clinical data from NCIC dataset; and B.G.N. supervised the research. S.G. and B.G.N. wrote the manuscript.
ACKNOWLEDGMENTS
The authors thank Dr. Hiroshi Haeno from Kyushu University and Drs. Siv Sivaloganthan and Mohammad Kohandel from Waterloo University for helpful discussions. This work was supported by grants from the Terry Fox Foundation (TFPPG 020003), National Cancer Institute (R37 49152), and Department of Defense (60167542 104703 A) (to B.G.N.), and a DCA Award from University of Toronto (to S.G.). B.G.N. was a Canadian Research Chair, Tier 1, and work in his laboratory was supported in part by a grant from the Ontario Ministry of Health and Long Term Care and the Princess Margaret Cancer Foundation.