Population modeling of tumor growth curves and the reduced Gompertz model improve prediction of the age of experimental tumors

Tumor growth curves are classically modeled by means of ordinary differential equations. In analyzing the Gompertz model several studies have reported a striking correlation between the two parameters of the model, which could be used to reduce the dimensionality and improve predictive power. We analyzed tumor growth kinetics within the statistical framework of nonlinear mixed-effects (population approach). This allowed the simultaneous modeling of tumor dynamics and inter-animal variability. Experimental data comprised three animal models of breast and lung cancers, with 833 measurements in 94 animals. Candidate models of tumor growth included the exponential, logistic and Gompertz models. The exponential and—more notably—logistic models failed to describe the experimental data whereas the Gompertz model generated very good fits. The previously reported population-level correlation between the Gompertz parameters was further confirmed in our analysis (R2> 0.92 in all groups). Combining this structural correlation with rigorous population parameter estimation, we propose a reduced Gompertz function consisting of a single individual parameter (and one population parameter). Leveraging the population approach using Bayesian inference, we estimated times of tumor initiation using three late measurement timepoints. The reduced Gompertz model was found to exhibit the best results, with drastic improvements when using Bayesian inference as compared to likelihood maximization alone, for both accuracy and precision. Specifically, mean accuracy (prediction error) was 12.2% versus 78% and mean precision (width of the 95% prediction interval) was 15.6 days versus 210 days, for the breast cancer cell line. These results demonstrate the superior predictive power of the reduced Gompertz model, especially when combined with Bayesian estimation. They offer possible clinical perspectives for personalized prediction of the age of a tumor from limited data at diagnosis. The code and data used in our analysis are publicly available at https://github.com/cristinavaghi/plumky.Author summary: Mathematical models for tumor growth kinetics have been widely used since several decades but mostly fitted to individual or average growth curves. Here we compared three classical models (exponential, logistic and Gompertz) using a population approach, which accounts for inter-animal variability. The exponential and the logistic models failed to fit the experimental data while the Gompertz model showed excellent descriptive power. Moreover, the strong correlation between the two parameters of the Gompertz equation motivated a simplification of the model, the reduced Gompertz model, with a single individual parameter and equal descriptive power. Combining the mixed-effects approach with Bayesian inference, we predicted the age of individual tumors with only few late measurements. Thanks to its simplicity, the reduced Gompertz model showed superior predictive power. Although our method remains to be extended to clinical data, these results are promising for the personalized estimation of the age of a tumor from limited measurements at diagnosis.