Abstract
When examining some factors that contribute to the growth or decline of a population or tumor, it is essential to consider a random hypothesis. By analyzing the effects of stress on a population (or volume of tumor growth) in a random environment, we develop stochastic models describing the dynamics of the population (or tumor growth) based on random adjustments to the population’s intrinsic growth rate, carrying capacity, and harvesting efforts (or tumor treatments). Apart from the models’ ability to capture fluctuations, the availability of a shape parameter in the models gives it the flexibility to describe a variety of population/tumor data with different shapes. The distribution of the stressed population size with or without harvesting (or treatments) is derived and used to calculate the maximum expected amount of harvests that can be taken from the population without depleting resources in the long run (or the minimum amount of chemotherapy needed to cause shrinkage or eradication of a tumor). The work done is applied to analyze tumor growth using published data comprising of the volume of breast tumor obtained by orthotopically implanting LM2-\(4^{LUC+}\) cells into the right inguinal mammary fat pads of 6- to 8-week-old female Severe Combined Immuno-Deficient mice.
Similar content being viewed by others
Data Availibility
The data used in this study are published data and can be found in the work of Vaghi et al. (2020).
References
Aggrey SE (2002) Comparison of three nonlinear and spline regression models for describing chicken growth curves. Poult Sci 81:1782–8 (PMID: 12512566)
Angeli S, Stylianopoulos T (2016) Biphasic modeling of brain tumor biomechanics and response to radiation treatment. J Biomech 49:1524–1531
Asadi M, Crescenzo AD, Sajadi FA, Spina S (2020) A generalized Gompertz growth model with applications and related birth-death processes. Ricerche mat. https://doi.org/10.1007/s11587-020-00548-y
Bateman H (1953) Higher transcendental functions, vol 1. McGraw-Hill, New York
Beddington JR, May RM (1977) Harvesting natural populations in a randomly fluctuating environment. Science 197:463–465
Bhattacharya R, Waymire EC (2021) Random walk, brownian motion, and martingales, graduate texts in mathematics book series, Springer Nature Switzerland, vol. 292, pp 191–198 (2021)
Chapman DG. Stochastic models in animal population ecology. Animal Population Ecology; Fifth Berkeley Symposium, pp 147–162
Chau NP (2000) Destabilizing effect of periodic harvest on population dynamics. Ecol Modell 127:1–9
Coren RL (2001) The importance of randomness in population studies. BioScience 51(8):604–605
da SilveiraCosta MI (2007) Harvesting induced fluctuations: insights from a threshold management policy. Mat. Biosci 205:77–82
Demicheli R (1980) Growth of testicular neoplasm lung metastases: tumor-specific relation between two Gompertzian parameters. Eur J Cancer 16(12):1603–1608. https://doi.org/10.1016/0014-2964(80)90034-1
Dempster JP, Pollard E (1981) Fluctuations in resource availability and insect populations. Oecologia 50:412–416
Dorini Fabio A, Cecconello MS, Dorini Leyza B (2016) On the logistic equation subject to uncertainties in the environmental carrying capacity and initial population density. Commun Nonlinear Sci Numer Simulat 33:160–173
Dorini Fabio A, Bobko Nara, Dorini Leyza B (2018) A note on the logistic equation subject to uncertainties in parameters. Comput Appl Math 37:1496–1506
Eric-Jan W, Simon F (2004) AIC model selection using Akaike weights. Psychon Bul Rev 11(1):192–196
Fornalski KW, Reszczynska J, Dobrzynski L, Wysocki P, Janiak MK (2020) Possible source of the Gompertz law of proliferating cancer cells: mechanistic modeling of tumor growth. Acta Phys Pol A. https://doi.org/10.12693/APhysPolA.138.854
Friberg S, Mattson S (1997) On the growth rates of human malignant tumors: implications for medical decision making. J Surg Oncol 65:284–297
Gard TC (1992) Stochastic models for toxicant-stressed populations. Bull Math Biol 54(5):827–837
Gatenby RA, Artzy-Randrup Y, Epstein T, Reed DR, Brown JS (2020) AEradicating metastatic cancer and the eco-evolutionary dynamics of anthropocene extinctions. Cancer Res 80(3):613–623
Goel NS, Maitra SC, Montroll EW (1971) On the volterra and other nonlinear models of interacting populations. Rev Mod Phys 43(2):231
Goh BS (1977) Stability in a stock-recruitment model of an exploited fishery. Math Biosci 33:359
Hallam TG, Zhien M (1986) Persistence in population models with demographic fluctuations. J Math Biol 24:327–339
Halmi MIE, Shukor MS, Johari WLW, Shukor MY (2014) Evaluation of several mathematical models for fitting the growth of the algae Dunaliella tertiolecta. Asian J Plant Biol 2(1):1–6
Harkos C, Svensson SF, Emblem KE, Stylianopoulos T (2022) Inducing biomechanical heterogeneity in brain tumor modeling by MR elastography: Effects on tumor growth, vascular density and delivery of therapeutics. Cancers 14:884. https://doi.org/10.3390/cancers14040884
Hartvigsen G. Carrying capacity, concept of. In: Levin, SA (ed) Encyclopedia of biodiversity, Academic Press, pp 641–649
Horsthemke W, Lefever R (1984) Noise-induced transitions. Theory and applications in physics, chemistry, and biology. Springer-Verlag, Berlin, Heidelberg
Jazwinski AH (1970) Stochastic processes and filtering theory. Academic, New York
Jonzén N, Ranta E, Lundberg P, Kaitala V, Lindén H (2003) Harvesting induced fluctuations? Wildl Biol 9:59
Kloeden PE, Platen E (1995) Numerical solution of stochastic differential equations. Springer-Verlag, New York
Lagarias JC, Reeds JA, Wright MH, Wright PE (1998) Convergence properties of the Nelder–Mead simplex method in low dimensions. SIAM J Optim 9(1):112–147
Laird AK (1964) Dynamics of tumor growth. Br J Cancer 13:490–502. https://doi.org/10.1038/bjc.1964
Levenberg K (1944) A method for the solution of certain non-linear problems in least squares. Q Appl Math 2(2):164–168
Levins R (1969) The effect of random variations of different types of population growth. Zoology 62:1061–1065
Lewontin RC, Cohen D (1969) On population growth in a randomly varying environment. Zoology 62:1056–1060
Li Wenxue, Wang Ke (2010) Optimal harvesting policy for general stochastic logistic population model. J Math Anal Appl 368:420–428
Liu M, Bai C (2015) Optimal harvesting of a stochastic logistic model with time delay. J Nonlinear Sci 25:277–289
Lungu EM, Oksendal B (1997) Optimal harvesting from a population model in a stochastic crowded environment. Math Biosc 145:47–75
Marquardt DW (1963) An algorithm for least-squares estimation of nonlinear parameters. J Soc Indust Appl Math 11(2):431–441
Moummou EK, Gutierrez-Sanchez R, Melchor MC, Ramos-Ábalos E (2014) A stochastic Gompertz model highlighting internal and external therapy function for tumour growth. Appl Math Comput 246:1–11
Mummert A, Otunuga OM (2019) Parameter identification for a stochastic SEIRS epidemic model: case study influenza. J Math Biol 79(2):705–729. https://doi.org/10.1007/s00285-019-01374-z
Norton L (1988) A Gompertzian model of human breast cancer growth. Cancer Res 48(24):7067–7071
Norton L, Simon R, Brereton HD, Bogden AE (1976) Predicting the course of Gompertzian growth. Nature 264(5586):542–545
Øksendal B (1987) Stochastic differential equations: an introduction with applications, 5th edn. Springer-Verlag, Heidelberg, New York
Olver FW, Lozier DW, Boisvert RF, Clark CW (2010)NIST Handbook of mathematical functions, National Institute of Sciences and Technology, U.S. Department of Commerce and Combridge University Press, (Dubai, Tokyo and New York U.S.A.)
Otunuga OM (2021) Time-dependent probability distribution for number of infection in a stochastic SIS model: case study COVID-19. Chaos Solit Fract 147:110983
Otunuga OM (2021) Time-dependent probability density function for general stochastic logistic population model with harvesting effort. Physica A 573:125931
Otunuga OM (2022) Stochastic modeling and forecasting of Covid-19 deaths: analysis for the fifty states in the United States. Acta Biotheoretica 70:25. https://doi.org/10.1007/s10441-022-09449-z
Parfitt AM, Fyhrie DP (1997) Gompertzian growth curves in parathyroid tumours: further evidence for the setpoint hypothesis. Cell Prolif 30(8–9):341–349
Prajneshu (1980) Time dependent solution of the logistic model for population growth in random environment. J Appl Prob 17
Pray L, Hartwell’s yeast LH (2008) A model organism for studying somatic mutations and cancer. Nat Educ 1(1):183
Protter P (2005) Stochastic integration and differential equations. Springer-Verlag, vol. 21
Renshaw E (1991) Modeling biological populations in space and time, Cambridge University Press, pp 6–9
Richards FJ (1959) A flexible growth function for empirical use. J Exp Bot 10(2):290–301. https://doi.org/10.1093/jxb/10.2.290
Ricker WE (1979) Growth rates and models. Fish physiology. Academic Press, London, pp 677–743
Schaefer MB (1954) Some aspects of the dynamics of populations important to the management of commercial marine fisheries. Bull Inter-Am Trop Tuna Comm 1(2):27–56
Silvert W (1978) Anomalous enhancement of mean population levels by harvesting. Math Biosci 42:253
Steel GG (2002) Growth kinetics of tumors. Oxford: Clarendon Press. ISBN 0-19-857388-X. Math. Biosci. 180, 29–48
Tang CY, Chen SX (2009) Parameter estimation and bias correction for diffusion processes. J Econom 149:65–81
Tjørve E, Tjørve KMC (2010) A unified approach to the Richards-model family for use in growth analyses: Why we need only two model forms. J Theor Biol 267:417–25. https://doi.org/10.1016/j.jtbi.2010.09.008. (PMID: 20831877)
Tjørve KMC, Tjørve E (2017) The use of Gompertz models in growth analyses, and new Gompertz-model approach: An addition to the Unified-Richards family. PLOS ONE 12(6):e0178691. https://doi.org/10.1371/journal.pone.0178691
Tvedt J (1995) Market structure, freight rates and assets in bulk shipping. Dr. Ocean dissertation, Norwegian school of economics and business administration, Bergen, Norway
Uyenoyama M (2004) The Evolution of population biology, Cambridge University Press. pp 1–19
Vaghi C, Rodallec A, Fanciullino R, Ciccolini J, Mochel J, Mastri M, Poignard C, Ebos JM, Benzekry S (2020) Population modeling of tumor growth curves and the reduced Gompertz model improve prediction of the age of experimental tumors. PLoS Comput Biol 16(2):e1007178
van Kampen NG (1981) Itô versus Stratonovich. J Stat Phys 24:175
West BJ, Bulsara AR, Lindenberg K, Seshadri V, Shuler KE (1979) Stochastic processes with non-additive fluctuations: I. Itô and Stratonovich calculus and the effects of correlations. Physica A 97(2):211–233
Wong E, Zakai M (1965) On the convergence of ordinary integrals to stochastic integrals. Ann Math Stat 36(5):1560–1564
Wu K, Wang DQ, Sornette D (2020) Generalized logistic growth modeling of the COVID-19 outbreak: comparing the dynamics in the 29 provinces in China and in the rest of the world. Nonlinear Dyn 101:1561–1581
Funding
No funding was received for conducting this study. The research does not include any human or animal participants.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
There is no conflict of interest concerning the work done here. All authors certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A: Simulation results for the volume of tumor for remaining specimens using models A, B, C, D
Figures 11, 12 and 13 show the simulation results for the remaining mice specimens using models A, B, C, and D.
Appendix B: Table containing goodness of fit estimates for models A, B, C, and D for all 66 specimens
Here, we use \(L=100\) sample paths in the LS-LM estimation process to estimate the parameters in models A, B, C, D used in describing the volume of tumor for the remaining specimens not included in Fig. 10. The tumor growth for these specimens are reported in Figs. 11, 12 and 13, and their corresponding goodness-of-fit estimates reported on Tables 5 and 6.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Otunuga, O.M. Tumor growth and population modeling in a toxicant-stressed random environment. J. Math. Biol. 88, 18 (2024). https://doi.org/10.1007/s00285-023-02035-y
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00285-023-02035-y