×

Reaction-difusion model of early carcinogenesis: the effects of influx of mutated cells. (English) Zbl 1337.92047

Summary: In this paper we explore a new model of field carcinogenesis, inspired by lung cancer precursor lesions, which includes dynamics of a spatially distributed population of pre-cancerous cells \(c(t,x)\), constantly supplied by an influx \(\mu\) of mutated normal cells. Cell proliferation is controlled by growth factor molecules bound to cells, \(b(t,x)\). Free growth factor molecules \(g(t,x)\) are produced by precancerous cells and may diffuse before they become bound to other cells. The purpose of modelling is to investigate the existence of solutions, which correspond to formation of multiple spatially isolated lesions of pre-cancerous cells or, mathematically, to stable spike solutions. These multiple lesions are consistent with the field theory of carcinogenesis. In a previous model published by these authors, the influx of mutated cells was equal to zero, \(\mu=0\), which corresponded to a single pre-malignant colony of cells. In that model, stable patterns appeared only if some of the growth factor was supplied from outside, arguably, a biologically tenuous hypothesis. In the present model, when \(\mu>0\), that hypothesis is no more required, which makes this model more realistic. We present a range of results, both mathematical and computational, which taken together allow understanding the dynamics of this model. The equilibrium solutions in the current model result from the balance between new premalignant colonies being initiated and the old ones dying out.

MSC:

92C37 Cell biology
35K57 Reaction-diffusion equations
35B35 Stability in context of PDEs
35B40 Asymptotic behavior of solutions to PDEs
92C50 Medical applications (general)
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] H. Awaya, Y. Takeshima, O. Furonaka, N. Kohno, K. Inai. Gene amplification and protein expression of EGFR and HER2 by chromogenic in situ hybridisation and immunohistochemistry in atypical adenomatous hyperplasia and adenocarcinoma of the lung. J. Clin. Pathol., 58 (2005), 1076-1080.
[2] R. Axelrod, D. E. Axelrod, K. J. Pienta. Evolution of cooperation among tumour cells. Proc. Nat. Acad. Sci., 103 (2006), 13474-13479.
[3] K. Chueh, C. Conley, J. Smoller. Positively invariant regions for systems of nonlinear diffusion equations. Ind. Univ. Math. J., 26 (1977), 373-392. · Zbl 0368.35040
[4] D. Henry. Geometric theory of semilinear parabolic equations. Springer-Verlag, 1981. · Zbl 0456.35001
[5] T. Hillen. A classification of Spikes and Plateaus. SIAM Review, 2007, in press. · Zbl 1160.35033
[6] E. Kamke. Differentialgleichungen: Losungsmethoden und Losungen. In: Gewohnliche Differentialgleichungen B. G. Teubner, 1977.
[7] T. Kawai, S. Hiroi, K. Nakanishi, A. K. Meeker. Telomere length and telomerase expression in atypical adenomatous hyperplasia and small bronchioloalveolar carcinoma of the lung. Am. J. Clin. Pathol., 127 (2007), 254-262.
[8] K. M. Kerr. Pulmonary preinvasive neoplasia. J. Clin. Pathol., 54 (2001), 257-271.
[9] M. Kimmel, D. E. Axelrod. Fluctuation test for two-stage mutations: application to gene amplification. Mutat. Research, 306 (1994), 45-60.
[10] K. Kishi, S. Homma, A. Kurosaki, S. Tanaka, H. Matsushita, K. Nakata. Multiple atypical adenomatous hyperplasia with synchronous multiple primary bronchioloalveolar carcinomas. Intern. Med., 41 (2002), 474-477.
[11] H. Kitagawa, A. Goto, T. Niki, M. Hironaka, J. Nakajima, M. Fukayama. Lung adenocarcinoma associated with atypical adenomatous hyperplasia. A clinicopathological study with special reference to smoking and cancer multiplicity. Pathol. Int., 53 (2003), 823827.
[12] A. Marciniak-Czochra, M. Kimmel. Modelling of early lung cancer progression: Influence of growth factor production and cooperation between partially transformed cells. Math. Mod. Meth. Appl. Sci., 17 (2007), 1693-1719. · Zbl 1135.92019
[13] A. Marciniak-Czochra, M. Kimmel. Dynamics of growth and signalling along linear and surface structures in very early tumours. Comp. Math. Meth. Med., 7 (2006), 189-213. 113 A. Marciniak-Czochra and M. KimmelModel of early carcinogenesis · Zbl 1111.92031
[14] A. Marciniak-Czochra, M. Kimmel. Reaction-diffusion approach to modelling of the spread of early tumours along linear or tubular structures. J. Theor. Biol., 244 (2006), 375-387.
[15] A. Marciniak-Czochra, M. Kimmel. Mathematical model of tumour invasion along linear or tubular structures. Math. Comp. Modelling, 41 (2005), 1097-1108. · Zbl 1080.92046
[16] A. Marciniak-Czochra, M. Ptashnyk. Derivation of a macroscopic receptor- based model using homogenisation techniques. SIAM J. Math. Anal., 2008, to appear. · Zbl 1178.35050
[17] A. Marciniak-Czochra. Developmental models with cell surface receptor densities defining morphological position. PhD thesis, University of Heidelberg, 2004. · Zbl 1104.92304
[18] L. Morandi, S. Asioli, A. Cavazza, A. Pession, S. Damiani. Genetic relationship among atypical adenomatous hyperplasia, bronchioloalveolar carcinoma and adenocarcinoma of the lung. Lung Cancer, 56 (2007), 35-42.
[19] J. D. Murray. Mathematical biology. Springer-Verlag, Berlin, 2003. · Zbl 1006.92002
[20] F. Rothe. Global solutions of reaction-diffusion systems. Springer-Verlag, Berlin, 1984. · Zbl 0546.35003
[21] T. Sano, Y. Kitayama, H. Igarashi, M. Suzuki, F. Tanioka, K. Chida, K. Okudela, H. Sugimura. Chromosomal numerical abnormalities in early stage lung adenocarcinoma. Pathol. Int., 56 (2006), 117-125.
[22] J. Smoller. Shock-waves and reaction-diffusion equations. Springer-Verlag, Berlin, 1994. · Zbl 0807.35002
[23] K. Takamochi, T. Ogura, K. Suzuki, H. Kawasaki, Y. Kurashima, T. Yokose, A. Ochiai, K. Nagai, Y. Nishiwaki, H. Esumi. Loss of heterozygosity on chromosomes 9q and 16p in atypical adenomatous hyperplasia concomitant with adenocarcinoma of the lung. Am. J. Pathol., 159 (2001), 1941-1948. 114
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.