Evasion of tumours from the control of the immune system: consequences of brief encounters
© Al-Tameemi et al.; licensee BioMed Central Ltd. 2012
Received: 6 February 2012
Accepted: 26 July 2012
Published: 25 September 2012
In this work a mathematical model describing the growth of a solid tumour in the presence of an immune system response is presented. Specifically, attention is focused on the interactions between cytotoxic T-lymphocytes (CTLs) and tumour cells in a small, avascular multicellular tumour. At this stage of the disease the CTLs and the tumour cells are considered to be in a state of dynamic equilibrium or cancer dormancy. The precise biochemical and cellular mechanisms by which CTLs can control a cancer and keep it in a dormant state are still not completely understood from a biological and immunological point of view. The mathematical model focuses on the spatio-temporal dynamics of tumour cells, immune cells, chemokines and “chemorepellents” in an immunogenic tumour. The CTLs and tumour cells are assumed to migrate and interact with each other in such a way that lymphocyte-tumour cell complexes are formed. These complexes result in either the death of the tumour cells (the normal situation) or the inactivation of the lymphocytes and consequently the survival of the tumour cells. In the latter case, we assume that each tumour cell that survives its “brief encounter” with the CTLs undergoes certain beneficial phenotypic changes.
We explore the dynamics of the model under these assumptions and show that the process of immuno-evasion can arise as a consequence of these encounters. We show that the proposed mechanism not only shape the dynamics of the total number of tumor cells and of CTLs, but also the dynamics of their spatial distribution. We also briefly discuss the evolutionary features of our model, by framing them in the recent quasi-Lamarckian theories.
Our findings might have some interesting implication of interest for clinical practice. Indeed, immuno-editing process can be seen as an “involuntary” antagonistic process acting against immunotherapies, which aim at maintaining a tumor in a dormant state, or at suppressing it.
This article was reviewed by G. Bocharov (nominated by V. Kuznetsov, member of the Editorial Board of Biology Direct), M. Kimmel and A. Marciniak-Czochra.
Cancer research, both experimental and theoretical, in last 12 years has been deeply influenced by the paper “The Hallmarks of Cancer” by . In this paper six key aspects (“hallmarks”) of cancer development and growth were identified and examined. Their interplay with other cellular populations was mainly seen as cooperative, and thus positive for the tumour. Recently, the same authors have published an interesting follow-up paper, “Hallmarks of Cancer: The Next Generation” , where their description of cancer development and growth was modified and “up-dated”. In particular, among the various new topics discussed, two important new aspects were added: the role of epigenetic phenomena and the possibility of competitive interplay with the innate and adaptive immune systems. In particular, evasion from immune destruction is explicitly listed as a new “hallmark”.
Tumour cells are characterized by a large number of genetic and epigenetic events leading to the appearance of specific antigens (e.g. mutated proteins, under/over-expressed normal proteins and many others) triggering reactions by the both the innate and the adaptive immune system [3–7]. These observations have provided a theoretical basis to the empirical hypothesis of immune surveillance, i.e. that the immune system may act to eliminate tumours , only recently experimentally and epidemiologically confirmed . Of course, the competitive interaction between tumour cells and the immune system involves a considerable number of events and molecules, and as such is extremely complex.
Moreover, to describe fully these immuno-oncological dynamics, one has to take into account a range of spatial phenomena, which are of outmost relevance in determining the dynamics of both immunogenic and non-immunogenic tumours [10–12]. In particular, the interplay between tumour cells and the immune system is strongly influenced by the spatial mobility of both tumour cells and cells of the immune system i.e. effector cells . Apart from the random motion of both types of cell, a prominent role is played by chemoattraction of effector cells towards the tumour cells. Indeed, chemotactic motion of immune system cells is a hallmark of the defence of the human body against “non-self agents”, including tumours, since cells belonging to both the innate immune system (e.g. Natural Killers, Macrophages, Dendritic Cells, etc. ) and adaptive immune system (e.g. Cytotoxic T Lymphocites, etc. ) are able to reach their targets thanks to the gradients of various kinds of chemicals [3, 14], e.g. inflammation-related substances produced by tumour cells. Thus, chemotaxis is of paramount importance in the interplay between tumours and the immune system, since it influences the control of tumour growth and also the immune surveillance.
However, besides temporal and spatial non-linearities, another important point to stress is that the structure of the above-mentioned interactions is also characterized by a series of evolutionary phenomena. As is self-evident, the immune system is not able to eliminate all neoplasms. In other cases, a dynamic equilibrium may also be established, such that the tumour may survive in a so-called “dormant state” [15–18], which is undetectable (i.e. the tumour persists at a very low, undetectable level of cells but is not completely eliminated by the immune system). Until recently this was largely inferred from clinical data, but  have been able to show experimentally, through an ad hoc mouse model, that adaptive immunity can maintain an occult cancer in an equilibrium state. It is quite intuitive that this equilibrium can be disrupted by sudden events affecting the immune system. If disease-related impairments of the innate and adaptive immune systems, or immuno-suppressive treatments preceding organ transplantations occur, then tumour regrowth occurs [9, 19]. This has been shown both by mouse models and through epidemiological studies [9, 19].
However, there is a major class of causes of disruption of the equilibrium that is not related to immuno-suppression. Over a long period of time , a neoplasm may develop multiple strategies to circumvent the action of the immune system [5, 9], which may allow it to recommence growing [9, 18] into clinically apparent tumours , which theoretically can reach their maximum carrying capacity . From an ecological point of view, we could say that the tumour has adapted to survive in a hostile environment, in which the anti-tumour immune response is activated [9, 18]. For example, the tumour may develop mechanisms to grow and spread by reducing its immunogenicity [5, 9]. In other words, the immunogenic phenotype of the tumour is influenced by the interaction with the immune system of the host. For this reason, the theory of the interactions between a tumour and the immune system has been called immuno-editing theory .
An impressive body of research is accumulating on immuno-evasive strategies, and a recent monograph  has been devoted to some aspects of this fascinating subject and to its close relationship with the effectiveness of immunotherapies. As far as the mathematical modelling of tumour and immune system interactions is concerned, there are many papers in the current literature which use deterministic models [13, 17, 18, 21–30] or stochastic models [31–35], as well as models introduced by Bellomo based on the kinetic theories of nonlinear statistical mechanics [36, 37]. The general approach of Bellomo’s theory is based on the concept of changes of activities of both the tumour cells and the effector cells of the immune system after encounters between them. As far as spatial aspects are concerned, [38, 39] developed a detailed spatio-temporal model focused on the role of macrophages. They showed that the presence of chemoattraction of macrophages towards the tumour cells implies both the onset of traveling waves and a heterogeneous spatial distribution of the tumour cells (see also ). Matzavinos, Chaplain and Kuznetsov proposed a spatiotemporal model of the interactions between tumour cells and cytotoxic T-lymphocytes (CTLs) [13, 16] by including the spatial motility of both tumour cells and CTLs, as well as chemotactic motion of the CTLs. They focused mainly on the role of the immune system in determining dormant states of the tumour, by showing, through a series of simulations, that a dormant state is reached, but the tumour cells are spatially distributed in an irregular pattern, which also temporally oscillates in a non-periodic fashion (see also ). In [18, 28], the immuno-editing phenomenon was empirically modelled by allowing the presence of slowly time-varying generic parameters in deterministic models (with time-scales significantly longer than those typical of the tumour-immune system interaction). Recently, in the framework of the above-mentioned kinetic approach, a generic model has been proposed for the learning ability of effector cells and for the hiding of tumour cells . In this paper, based on the concept introduced by  that the immune system has the ability of “sculpting the phenotype” of tumour cells (i.e. promoting the change towards less imunogenic and more resistant phenotypes), we propose a cell-centered semi-mechanistic approach aimed at describing a possible immunologically realistic kinetic mechanism through which immuno-evasion begins. Since there is strong experimental evidence that type, density and location of CTLs are predictive of the clinical outcome of some tumours, such as colorectal tumours , and since we are interested in the long-time dynamics, here we shall deal with the interplay of a neoplasm with CTLs. In our model, we suppose that the tumour cells that survive an attack by CTLs have a probability of acquiring (through mutations or even by epigenetic changes) a phenotype that is more resistant to future attacks by CTLs. In turn, at each new encounter with a CTL, this resistance can be increased further, and after a finite number of encounters a complete or maximal resistance to specific immunity is acquired. Moreover, specific spatial effects may be linked to the immuno-evasion of neoplasms. Indeed, recently  showed experimentally that tumour cells can produce chemical CXCL12 that, for large concentrations, act as chemorepellent for CTLs, whereas for low concentration it acts as a normal chemoattractor. We integrate - with some simplifications - these experimental findings in our model by permitting in the range of features defining the increasingly resistant tumour cell phenotypes an increasing ability to produce such chemorepulsive substances, whereas in a future model we shall consider the above described “nonmonotone” behavior of the taxis. These two bio-theoretical hypotheses, although new, are in line with the general schema of tumour cell escape from the immune response. Indeed, as stressed by , tumour cells may escape from immune control through two general mechanisms: (a) mechanisms that involve the secretion of soluble factors; (b) mechanisms that are dependent on the contact between the tumour cells and the effectors and that are aimed at reducing antigen recognition/adhesion and apoptotic resistance. Given the current experimental knowledge the above mentioned factors are primarily aimed - apart from, in many cases, their mitogenic action - at inducing the emergence of immunosuppressive networks . In our present model, the factors, in line with the animal model by , are chemicals that repel CTLs. Finally, in the concluding remarks, we shall also discuss, from and evolutionary point of view, the differences of our model with the current evolutionary view of the immuno-editing.
The Mathematical Model
Following the kinetic scheme employed by , in absence of immuno-editing mechanisms the interplay between tumour cells and tumour-infiltrating cytotoxic-T-lymphocytes can be modelled as shown in Figure 1 (see: ), where T denotes a tumour cell, E denotes an effector cell (CTL), C denotes the complex formed, T∗denotes a dead tumour cell and E∗denotes a dead effector cell. The following assumptions are made:
the break-up of complexes can lead to a situation where both the tumour cell and the CTLs are alive with a rate k−
the break-up of complexes can lead to a situation where either the immune cell or the tumour cell survives the encounter with a rate k
the probability that a tumour cell is killed is p, and correspondingly the probability that a CTL is killed (i.e. the tumour cells survives) is (1−p)
The key idea proposed in this paper which develops the work of , is that a proportion of the tumour cells that survive an encounter with a CTL are more resistant to any future attacks by CTLs. Consequently, the phenotypic properties of these new “enhanced” tumour cells will be different from those of the “naive” tumour cells. Specifically, we make the additional assumptions:
their probability of being killed (previously the parameter p) is smaller
their probability of being recognized and also of forming a complex with a CTL (embedded in the parameter k+ ) is smaller
Moreover, we shall also assume that the proliferation rate of CTLs stimulated by the presence of the complexes is also smaller. We denote the naive tumour cells by T0(t) and the non-naive tumour cells by T i , where i stands for the number of previous encounters with the CTLs. We assume that the fitness of tumour cells increases up to a maximum number of encounters N, implying that we consider in total 1 + N “classes” of tumour cells, T0,T1,…,T N .
The new kinetic relationships of our model are illustrated in Figure 2and are characterized by the following new groups of parameters:
the probability that a tumour cell of the i-th class is killed: p i . We assume that p i is decreasing with index i, with p N ≥0;
the probability of transition T i →Ti + 1to the state i: θ i . We assume that θ i is increasing for 0≤i≤N−1. Since we have assumed N classes of tumour cells, θ N =0.
where i=1,…,N, and l=0,…,N.
However, not only the temporal but also the spatio-temporal properties of the “fitter” tumour cells are likely to be different from those of the baseline tumour cells. Namely:
the production rates of chemoattractants stimulated by a complex CTL+’non naive tumour cell’ is assumed to be smaller than that of the naive cells
since recently Vianello et al.  showed in an animal model that tumours produce chemicals that repels the CTLs, here we assume that those chemorepellents are produced by the non-naive cells. For the sake of the precision, Vianello’s findings showed chemoattraction for low concentrations of chemical CXCL12 emitted by tumour cells. For the sake of simplicity here we shall only consider chemorepulsion.
In the following, we provide the full equations for all variables, including the spatial components. As mentioned in the previous section, we have assumed the model of  as our baseline model.
Spatiotemporal Dynamics of the Tumour Cells
where i=1,…,N, and where r1 is the baseline exponential growth rate of the tumour (i.e. its theoretical growth rate when it is ‘small’) and β1 is the inverse of its carrying capacity (in absence of immune reactions).
Spatiotemporal Dynamics of the CTLs
The proliferation of the CTLs, E, stimulated by the complexes C j is embedded in the rate constant q j , which, as a consequence, must be decreasing with the index j, with q N ≥0 (f , g are constant parameters). Note that in absence of immuno-editing this proliferation term reads fC/(g + T), and it has been has been introduced in [22, 45]. It represents the experimentally observed enhanced proliferation of CTLs in response to the tumour. This functional form is consistent with a model in which one assumes that the enhanced proliferation of CTLs is due to signals, such as released interleukins, generated by effector cells in tumour cell-CTL complexes. We note that the growth factors that are secreted by lymphocytes in complexes (e.g IL-2) act mainly in an autocrine fashion. That is to say they act on the cell from which they have been secreted and thus, in our spatial setting, their action can be adequately described by a “local” kinetic term only, without the need to incorporate any additional information concerning diffusivity. The external influx of CTLs is, for the sake of the simplicity, modelled as sh(x), where h(x) is a Heaviside function, taken to be zero over a given subregion of the domain of interest cf. . In other words, we assume that there is a subdomain where lymphocytes are not naturally present and which is penetrated by CTLs only thanks to diffusion and chemotaxis.
where we assume that the production rate constant Π i is decreasing with index i, with Π N ≥0, since we assume that complexes between CTLs and non-naive tumour cells are less and less able to produce such a chemoattractant.
Tumour cell-CTL Complexes
Modeling the transition rates and probabilities
and that their baseline value is sufficiently small: 10−5≤θ0≤10−3. In other words we assume that the probability of acquiring the less immunogenic phenotype is small (in analogy with the smallness of the probability of surviving to an attack by a CTLs).
where we assumed: w MAX ≈Π0(i.e. the maximum production rate of the chemorepellent is equal to the production rate of the chemoattractant in absence of immuno-editing phenomena).
Boundary and initial conditions
Note that E a is the baseline homogenous steady-state value of CTLs in the absence of tumour cells, and that T a is the the baseline homogenous steady-state value of the naive cells in absence of CTLs, whereas C a is the maximum possible density of complexes resulting from initial values of E and T.
Results and Discussion
Moreover  and may be either constant or linearly decreasing, with .
Since the average lifespan of a chimeric mouse is three years, and since we are interested in assessing the possibility (and spatio-temporal modality) of the onset of immunoevasion of a tumour, all simulations (unless stated) represent an interval of length 1100 days ≈3 years.
Spatially Homogenous Case
Moreover, Max(T0) seems to be a decreasing function of p N . These results might be explained as follows: the decrease of the competition between all the tumour cells and the immune system embedded in the decrease of the parameters , might shift the ‘internal’ competition between the naive and the non-naive tumour cells.
Finally, we note that the effect of chemorepulsion, which is well detectable at level of the spatial densities T i , is no more detectable at level of the total density of all tumour cells. However, this effect is again noticeable if we consider the density of T0 versus the density of T1 + ⋯ + T N (see the second plot of Figure 11).
In this paper we have presented a novel mathematical model of the immune response to cancer, focusing on the specific spatio-temporal response of cytotoxic T-lymphocytes to tumour cells. We have developed and extended ideas originally formulated by  by proposing a possible kinetic mechanism leading to tumour evasion from the immune control. Our model is based on the key concept that a tumour cell which survives the formation of a complex with a cytotoxic T-lymphocyte can develop, with a given probability, an increased probability of surviving further attacks by CTLs. We do not specify whether this so-called ‘increased resistance’ is genetic or epigenetic. Indeed, from a kinetic point of view, this is immaterial. However, in order to experimentally validate our hypothesis, this distinction would be of paramount relevance.
Note that the model by  is based on a mass-action law mechanism. The reaction kinetics in a crowded cellular environment could differ from that, as studied (in absence of immunoevasion) in [17, 27, 28]. We shall investigate these more realistic settings in future works, but we think that the basic results showed here should not substantially change.
In this work we have dealt with the spatio-temporal interplay between tumours and a specific immune response from CTLs. We chose this approach because of the experimental evidence on the relevance of CTLs in determining tumour dormancy or the evasion of many important tumours such as melanomas, ovarian carcinomas and colorectal carcinomas, where the presence of infiltrating lymphocytes is a useful prognostic marker [42, 46]. However, tumour immunoevasion from dormancy is a multi-faceted phenomenon. We stress here that by no means do we think that ours is an exhaustive theoretical treatment of a such complex phenomenon.
Concerning spatial issues we also briefly mention that also in case of small dormant tumours the next generation of spatial models of tumour growth should better stress and investigate the interplay of tissue 3D geometry and the tumour vascularization with phenotypic changes in tumour cells. We have built our model based on the tumour dormancy mathematical model of [13, 16], where parameters were fitted to experimental animal (mouse) data. However, embedding the proposed evolutionary mechanism in a more complex setting, where a more detailed description of both adaptive and innate immunity is included, should lead to results qualitatively similar to those here illustrated.
Our simulations suggest that the proposed mechanism is able to mimic various dynamics of immunoevasion during the lifespan of a mouse. We have also highlighted the differential spatiotemporal contributions to evasion due, respectively, to: i) a decrease in the probability p i of being lethally hit; ii) a decrease in the probability, embedded in , that a tumour cell is recognized by a CTL. In particular, our model suggests that a decrease in the parameters p i is needed to produce evasion, which does not occur in the case where p i remains constant at its baseline level inferred from the experimental data. However, the role of the parameters is important since it can greatly accelerate the simulated process. Moreover, our computational simulations also showed that the proposed mechanism can also deeply affect the spatial patterning of the tumour. In particular, our model suggests that to have a uniform invasion profile for the tumour cells necessitates also having a decrease in the recognition rate, embedded in the parameters . These parameters also differentially shape the spatial distribution of the various classes of tumour cells.
Concerning the possible chemorepulsion of CTLs, our computational simulation results showed that, in our biological settings, although it does not affect the spatiotemporal dynamics of the total number of tumour cells, it has a remarkable influence on the spatio-temporal distribution of the different individual classes of tumour cells. Further analysis is needed to ascertain if, with different parameters, the effect of this factor can be different, and in order to understand the behaviour in the current setting.
As far as the key ‘immuno-evasion-related’ parameters such as θ i , p i , and are concerned, we were not able to fit them with experimental data (apart, of course, from the values for p0 and , from [13, 16]) because in the literature, to the best of our knowledge, immuno-evasion of tumours is only illustrated by means of qualitative clinical or molecular experimental findings. In particular, no immuno-evasion-related tumour growth data are available. Indeed, a complete experimental kinetic study of the adaptive evasion from tumour dormancy allowing, for example, the plotting of tumour growth curves would currently be very difficult to undertake. Thus we hope that this theoretical work may contribute to triggering such experimental investigations, which would allow us to validate our model.
From a theoretical point of view, our model, although detailed and focused on a very specific aspect of immuno-oncology, and on some very specific mechanisms, is conceptually in line with the general theories by Bellomo [36, 37, 41], who considers tumour cells and immune system effector cells as “active particles” endowed with activities and properties. Indeed, also in this paper the changes of activities of cells upon encounters between tumour cells and effector cells of the immune system are central in determining the dynamics of the system.
To the best of our knowledge, the evolutionary nature of the immuno-editing process has been studied until now under the framework of the so-called “modern synthesis”, following which the environment (in our case the immune system) is not the “causative agency”  but a mere selective force promoting fixation of adaptive genomic changes . In the case of immunoevasion, a lowly immunogenic clone may appear spontaneously due to the large random mutation rate of tumour cells. This new phenotype is then involuntarily selected by the immune system, which kills the other phenotypes that remain strongly immunogenic. Thus the sculpting of tumour cells phenotypes  mentioned in the introduction, is involuntary, passive.
On the contrary, in our model the more immuno-resistant phenotypes may arise because of genetic or epigenetic causes - due to the interaction between the tumour cells and the immune system. This point of view, which is quasi-Larmackian , is in line with a number of recent discoveries that are leading to a new theory of “extended evolutionary synthesis” , which indeed investigates the impact of both genetic and epigenetic inheritance on evolutionary phenomena in order to decipher the complex interplay between genotypes, epigenotypes, phenotypes and environment. From a biological and biophysical point of view this implies also including two other key components: timescales (see, e.g., the review paper ) and the role of randomness. Following this extended and more modern perspective, we think that two general classes of evolutionary mechanisms might biologically underlie our kinetic model of transitions between phenotypes of tumour cells: stress-induced mutations and epigenetic switches.
In stress-induced mutagenesis [47, 49], some kind of stress induces genomic changes in a cell that, although probabilistic (e.g. in our model θ i ∈(0,1)), are triggered by the cell and its interplay with the “external world”, and that are beneficial for the cell [47, 49]. As a result, the stress-induced mutations are adaptive and beneficial, and thus they are natural candidates to explain biologically the kinetic mechanism of our mathematical model. Moreover, it is important to note that stress-induced genomic changes are a well-known and important mechanism in tumours [47, 50, 51]. Finally, Koonin hypothesized  that quasi-Larmarckian processes are essentially triggered by very strong signals (see e.g. Figure 3 of ), which, we remark, is the case for tumour-CTL interplays.
As far as the epigenetic path is concerned, it is now known that there are inheritable phenotypic changes without underlying genetic variations [48, 49], and that sometimes those changes are unusually rapid . The “epigenome” is dynamic and it reflects an individual’s or a tissue’s environmental exposure during the whole of its lifespan . Among the possible epigenetic changes, we mention specifically the methylation of DNA bases (“epimutations”), and also the switching between two different equilibrium states in the biochemical pathways influencing a set of inter-related phenotypes (e.g. “low immunogenicity” and “large immunogenicity”), due for example, to strong stress signals. Indeed, due to hysteresis phenomena inducing memory, also when an external stress signal is removed the system may not return to its original state. This memory may remain stable for many generations thus providing an example of epigenetic inheritance .
In this work we were interested essentially in the basic facts of the immune response to tumours. However, a number of immunotherapies have been proposed and also theoretically investigated (see [17, 26, 30] and references therein). We believe that both the experimental results concerning immunoevasion of tumours and the theoretical findings we have proposed here might have some implication of interest for clinical practice. Indeed, the immunoevasion mechanisms may be seen as a process counter-acting immunotherapies able to maintain a tumour in a dormant state. In general we share the opinion of , who stressed that recent progress in immuno-oncology has not influenced the way anti-cancer therapies are conceived and applied clinically. Indeed, we think that the immuno-editing process can be seen as an “involuntary” antagonistic process acting against immunotherapies. As far as this point is concerned, in future we will also investigate the possibility that the presence of an immunotherapy is somewhat sensed by a tumour, with a consequential acceleration of immuno-editing.
The model that was proposed here has to be understood as a detailed model at the level of the kinetics of the cellular populations of a possible mechanism that might enable tumour cells to evade from the control of adaptive immunity. The various specific (and tumour-dependent) strategies deployed by those cells in order to reach their aim, are phenomenologically described by means of the model of the dependence of the various parameters on the classes of tumour cells, as well as in the macroscopic modelling of the chemorepulsion. This is a first step in a research effort for a more complete description of tumour cell immunoevasion, which will include the detailed modelling of the biological mechanisms underlying those and other specific evasion strategies. Thus, given the complex network of interplaying between inter-cellular and intra-cellular signalling, and given the various temporal scales (from the rapid dynamics of the intracellular pathways involved, to the relatively slow growth of a tumour, up to the very slow onset of immunoevasion) as well as the spatial ones (from individual cells to visible tumours), a more detailed model will have to be multiscale. This will involve a wide array of computational tools, from those typical of computational biology and bioinformatics, to more classical analytical and numerical methods of statistical mechanics and mathematical physics.
The chosen non-null initial conditions for the naive tumour and immune cells represent a front of tumour cells encountering a front of CTLs, resulting in the formation of CTL-tumour cell complexes. In the absence of a tumour, the homogeneous steady-state density of the CTLs is s/d1 and therefore this is the value we have taken for the initial density E a of CTLs in the initial conditions. Similarly, in the absence of an immune response, the homogeneous steady-state density of the tumour cells is 1/β1 and this is what we take as the initial density of tumour cells T a in the initial conditions. Thus, when the fronts of the two cell populations meet, the maximum density of CTL-tumour cell complexes will be min(E a ,T a ) and hence our choice for C a . The function exp(−1000(x−l)2) was chosen to mimick the onset of sharp but continuous front.
Numerical values of the parameters
As far as the spatiotemporal dynamics of the chemorepellent is concerned, we assume that its diffusion coefficient and decay rate is the same as the chemokine αi.e. D ρ =D α , δ ρ =δ α , μ ρ =μ α , and ξ=γ.
The comments of the referees were reported in italics. All the three referees included minor comments on misprints or undefined parameters or other minor suggestions, which were all implemented. Thus, we only reported the comments of interest to the general readership.
Note for the referees: following your suggestions the revised version now contains an Appendix.
Reviewer #1: Prof. G. Bocharov (nominated by Dr. V. Kuznetsov, member of the Editorial Board of Biology Direct) (Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russian Federation)
The paper presents a theoretical study of the tumour immuno-evasion from dormancy. To this end a mathematical model of reaction-diffusion-chemotaxis type formulated with PDE is proposed. The model considers the spatio-temporal population dynamics of interactions between tumour cells and cytotoxic T cells. The key assumption of the model reflecting recent biological insights into the pathogenesis of the solid tumour growth states that the tumour cell population is heterogeneous with respect to the parameters characterizing the outcome of the interaction with cytotoxic T lymphocytes. The interaction with CTL appears to acts as a selection force shaping the microevolution of the tumour. The heterogeneity assumption enters the model via parameterized dependence of the rate of transition of tumour cells from naive to more mature states, the efficacy of CTL mediated killing and the production of chemicals acting as attractants and repellors for CTLs. Numerical simulations with the model for various parameters combinations show that the immuno-evasion can result from the phenotypic heterogeneity of the tumour cells dynamically adapting to and shaped by the anti-tumour immune response.
1. The issue of spatio-temporal modelling of tumour growth is an area of active research in cell population dynamics (e.g., the studies of Bertuzzi and Gandolfi). The related work on spatio-temporal modelling of tumour growth needs to be refereed to in the Introduction section.
We fully agree, and we added the required references.
Page 6: The authors make use of the Mass Action Law to describe the interaction between CTL and target cells via bilinear terms. Meantime, the reaction kinetics in crowded cell environments can differ from those assumed by the classical chemical kinetics. This issue need to be commented.
This point is very important, and we commented it in the conclusions.
Page 10: The parameter “1000” appearing in the initial functions has to be justified. Please, explain the choice of C a .
We added a full justification of the chosen initial conditions.
Page 13: For the readership, the parameters of the model should be presented in a Table with columns specifying notation, biological meaning, units and uncertainty ranges.
We added a table with average values of the parameters, as obtained in Ref. 13 by means of least squares algorithm, and their meaning.
Page 13, line 20: Can the parameters equality assumption for the repellor and attractant be biologically justified?
Unfortunately there is a lack of experimental results on this important data, so this is a pure assumption.
Page 14, line -5: Please, elaborate more on the specific choice of the parameter values.
We tried your suggestion, but the result was not less clear. Take also into the account that Biology Direct is purely online, thus there are no pages restrictions.
Page 18: Section 2-Dimensional Domain is based upon simulations with the chemotaxis and chemorepulsion not included in the model. The value of the results of the 2D model based simulations is not straightforward. Please, either expand the section by considering the same set of processes as in 1D case or remove it.
We removed it.
Conclusions section: It would be interesting to have the authors opinion on whether the tissue 3D geometry and its vascularisation represent important structural constraints affecting the genomic or phenotypic changes in tumour cells that need to be considered in the next generation of the spatially resolved models of tumour growth?
We fully agree, and we mentioned this important issue in the discussions.
Finally, all suggestions in your minor remarks were implemented.
Reviewer #2: Prof. M. Kimmel (Rice University, Houston, USA)
The paper concerns an extension of Chaplain group’s model (Matzavinos et al. 2004) of response of tumour cells to cytotoxic T-lymphocytes. The extension involves classes from naive through non-naive tumour cells, which have different ability to repel T-lymphocytes and other differing characteristics. The models consists of a complicated system of partial differential equations and is investigated purely by simulation.
In this first investigation we employed numerical analysis for two main reasons: first, we feel the immuno-editing phenomena as mainly a transitory effect, and as such requiring mainly numerical analysis; second, we wanted to write a manuscript aimed to theoretical biologists with biological background, such as the general readership of ’Biology Direct’. However, we also outlined a more compact version of the proposed model where the phenotypes are described by means of a continuum. For this version we soon start a full mathematical analysis.
In my opinion, stratification of tumour cells with respect to their immuno-naivete is a natural concept, which however, leads to many complications and ambiguities in model building. The paper only partially reviews the consequences of varying the hypotheses.
In the above mentioned continuum version we should be able of better exploring these important issues.
In particular, the 2-d model by the end of the paper seems superfluous, as it is not thoroughly investigated and seems to be plagued by boundary artifacts. I suggest removing it for the time being.
Following you suggestions, and the suggestions of another referee, we removed it.
I am not sure if the tumour immuno-editing can be considered a Lamarckian evolutionary process. Even if it is, I am not sure if this is important for the paper.
We respectfully disagree with you on this point. The possibility that immuno-editing might be, at least in some cases, a Lamarckian process, where the interplay between the different kind of cells is the ’driver’, is of some interest, in our opinion.
Page 2. I do not know if “extremely complex” implies “strongly non-linear” as the authors seem to claim.
We removed this observation.
Page 2.(bottom). It would be perhaps helpful to specify which cell types belong to the adaptive immune system and which to the innate one.
Page 5. “Since we are modeling a situation where immuno-evasion of the tumour cells is not considered …” Discussion of what this assumption really implies will be helpful.
We removed the above quite ambiguous phrase. Indeed, we wanted to means that the baseline model, with the considered values of its parameters, referred to a case where the tumour is not able (without immuno-editing) to evade the immune control.
Pages 9-13. Part of this material can probably go into an online Appendix.
Reviewer #3: Prof. A. Marciniak-Czochra (Heidelberg University)
The paper is devoted to modelling of spatio-temporal dynamics of interactions of cytotoxic T-lymphocytes with cells of a homogeneous avascular tumour. It takes into consideration a chemoattraction of the immune cells by the tumour and a chemorepulsion to chemicals produced by non-naive tumour cells. The model is based on the one presented in the paper “Mathematical Modelling of the spatio-temporal Response of Cytotoxic T-Lymphocyte to a Solid Tumour” (Matzavinos A, Chaplain M, Kuznetsov, Mathematical Medicine and Biology, 2004, 21:1-34) The novelty of the current work is in the assumption that the tumour cells, surviving encounter with the immune cells, undergo a functional change, which benefits their immune resistance. This resistance increases with each encounter until the maximal resistance is reached. The resulting system of coupled partial and ordinary di erential equations is solved numerically. The results are discussed in the biological context. I recommend the manuscript for publication in Biology Direct after some minor revisions related to the organisation of the paper. The issues, which should be taken into account before publication are (Note of authors:this referee listed a long series of misprints, which were all fixed. We only report the comments of interest for the reader):
P.8, the second Equation: f and g are not introduced. Although model equations are nearly identical to those in reference (Matzavinos et al, 2004), authors should motivate the term “proliferation and recruitment”.
In the revision, we fully explained the meaning of the term like fC/(g + T).
P.10: It would be helpful to have the initial conditions plotted.
Since the shape of the initial conditions are elementry, and since another referee asked to add some written comemnts on this point, we did not include the requested plots.
P.13: It should be mentioned how parameters are estimated.
The parameters were estimanted in Kuznetsov et al 1991 basically by means of mininmal least squares estimate algorithm.
Please motivate the choice of θ0; p N ; ; how sensitive are the results to the choice of the parameters?
We have choosen a small value for θ0because we think that the probability of acquiring the fitter phenotype is small, in analogy with the small probability of surviving to an attack by a CTL. The value of is taken from (Matzavinos et al 2004).
P. 28, Figure 12: Why (700,1100) and not (800, 1100)?
800 days would have been a too long time with respect to 700 (100 days in mice correspond approximatively to 2300 days in men).
Axis labels of the following figures are not readable: 9, 10, 11, 13, 15, 17, 20, 21.
Concerning the readability of the axis labels, it probably depend on the fact that our eps figures were converted by the publisher in jpeg figures. In our original pdf figures they seems readable. We shall carefully manage this point in case the manuscript should be accepted.
The following were comments on the revised manuscript (other comments on missprints are omitted, and were implemented):
The rescaling of equations have mistakes. For example on p. 30, λ should be,
Indeed, it was a misprint. We apologize.
We had forgot to stress that C a =E a , so what we wrote is correct.
I find the following phrase hardly readable “Finally, we note that the predicted effect on the tumour cells spatial distribution of the induction of chemorepulsion of CTLs is not detectable if we consider the total density of all tumour cells.” Please refomulate it.
We reformulated it.
AD and MC defined the mathematical model; MAT performed the simulations; AD, MC and MAT discussed the results; AD,MC and MAT wrote the manuscript. All authors read and approved the final manuscript.
We thank very much the three referees, whose suggestions helped us to greatly improve this work.
MAJC gratefully acknowledges the support of the ERC Advanced Investigator Grant 227619, ’M5CGS - From Mutations to Metastases: Multiscale Mathematical Modelling of Cancer Growth and Spread’.
The AD’s work has been done in the framework of the Integrated Project “P-medicine - from data sharing and integration via VPH models to personalized medicine” (project ID: 270089), which is partially funded by the European Commission under the 7th framework program.
- Hanahan D, Weinberg RA: The Hallmarks of Cancer. Cell. 2000, 100: 57-70. 10.1016/S0092-8674(00)81683-9.PubMedView ArticleGoogle Scholar
- Hanahan D, Weinberg RA: Hallmarks of Cencer: The Next Generation. Cell. 2011, 144: 646-674. 10.1016/j.cell.2011.02.013.PubMedView ArticleGoogle Scholar
- Delves P, Martin S, Burton D, Roitt I: Roitt’s Essent Immunol. 2011, Chichester: Wiley-BlackwellGoogle Scholar
- Kindt T, Osborne B, Goldsby R: Kuby Immunol. 2006, New York: W. H. FreemanGoogle Scholar
- Pardoll D: Does the Immune System See Tumors as Foreign or Self?. Annu Rev Immununol. 2003, 21: 807-839. 10.1146/annurev.immunol.21.120601.141135.View ArticleGoogle Scholar
- Swann J, Smyth M: Immune Surveillance of Tumors. J Clininical Investigations. 2007, 117: 1137-1146. 10.1172/JCI31405.View ArticleGoogle Scholar
- Yefenov E: Innate and Adaptive Immunity in the Tumor Microenvironment. 2008, New York: SpringerView ArticleGoogle Scholar
- Ehrlich P: Ueber den jetzigen stand der karzinomforschung. Nederlands Tijdschrift voor Geneeskunde. 2009, 53: 273-290.Google Scholar
- Dunn G, Old L, Schreiber R: The Three Es of Cancer Immunoediting. Annu Rev Immununology. 2004, 22: 329-360. 10.1146/annurev.immunol.22.012703.104803.View ArticleGoogle Scholar
- Adams V, Bellomo N: A Survey of Models for Tumor-Immune Syst Dyn. 1996, Boston: BirkhauserGoogle Scholar
- Bertuzzi A, d’ Onofrio A, Fasano A, Gandolfi A: Modelling cell populations with spatial structure: steady state and treatment-induced evolution of tumour cords. Discrete and Continuous Dynamical Syst-Ser B. 2004, 4: 115-134.Google Scholar
- Andasari V, Gerisch A, Lolas G, South A, Chaplain M: Mathematical modeling of cancer cell invasion of tissue: biological insight from mathematical analysis and computational simulation. J Math Biol. 2011, 63: 141-171. 10.1007/s00285-010-0369-1.PubMedView ArticleGoogle Scholar
- Matzavinos A, Chaplain M, Kuznetsov V: Mathematical Modelling of the spatio-temporal Response of Cytotoxic T-Lymphocyte to a Solid Tumour. Math Med Biol: A J of the IMA. 2004, 21: 1-34. 10.1093/imammb/21.1.1.View ArticleGoogle Scholar
- Keener J, Sneyd J: Math Physiol. 2003, New York: Springer-VerlagGoogle Scholar
- Koebel CM, Vermi W, Swann JB, Zerafa N, Rodig SJ, Old LJ, Smyth MJ, Schreiber RD: Adaptive immunity maintains occult cancer in an equilibrium state. Nat. 2007, 450: 903-907. 10.1038/nature06309.View ArticleGoogle Scholar
- Chaplain M, Matzavinos A: Mathematical Modelling of Spatio-temporal Phenomena In Tumour Immunology. Lecture Notes in Math. 2006, 1872: 131-183.View ArticleGoogle Scholar
- d’Onofrio A: A genaral Framework for Modeling Tumor-Immune System Competition and Immunotherapy: Analysis and Medical Inferences. Physica D. 2005, 208: 220-235. 10.1016/j.physd.2005.06.032.View ArticleGoogle Scholar
- d’Onofrio A: Tumor Evasion from Immune Control: Strategies of a MISS to Become a MASS. Chaos, Solitons and Fractal. 2007, 31: 261-268. 10.1016/j.chaos.2005.10.006.View ArticleGoogle Scholar
- Stewart TJ, Abrams SI: How Tumour Escape Mass Destruction. Oncogene. 2008, 27: 5894-5903. 10.1038/onc.2008.268.PubMedView ArticleGoogle Scholar
- Gabrilovich DI, Hurwitz AA: Tumor-Induced Immune Suppression. 2008, New York: SpringerView ArticleGoogle Scholar
- Stepanova N: Course of the Immune Reaction During the Development of a Malignant Tumor. Biophy. 1980, 24: 917-923.Google Scholar
- Kuznetsov VA, Makalkin IA, Taylor MA, Perelson AS: Nonlinear Dynmamics of Immunogenic Tumors: Parameter Estimation and Global Bifercation Analysis. Bull Math Biol. 1994, 56 (2): 295-321.PubMedView ArticleGoogle Scholar
- Kuznetsov V, Knott G: Modeling Tumor Regrowth and Immunotherapy. Math Comput Modell. 2001, 33: 1275-1287. 10.1016/S0895-7177(00)00314-9.View ArticleGoogle Scholar
- de Pillis L, Radunskaya A, Wiseman C: A Validated Mathematical Model of Cell-Mediated Immune Response to Tumor Growth. Cancer Res. 2005, 65: 7950-7958.PubMedGoogle Scholar
- Kogan Y, Forys U, Shukron O, Kronik N, Agur Z: Cellular Immunotherapy for High Grade Gliomas: Mathematical Analysis Deriving Efficacious Infusion Rates Based on Patient Requirement. SIAM J Appl Math. 2010, 70: 1953-1976. 10.1137/08073740X.View ArticleGoogle Scholar
- Arciero JC, Jackson TL, Kirschner DE: A Mathematical Model of Tumour-Immune Evasion and siRNA Treatment. Discrete and Continuous Dynamical Syst Ser B. 2004, 4: 39-58.Google Scholar
- d’Onofrio A: The Role of the Proliferation Rate of Effectors in the Tumor-Immune System Competition. Math Mododels and Methods in Appl Sci. 2006, 16: 1375-1401. 10.1142/S0218202506001571.View ArticleGoogle Scholar
- d’Onofrio A: Metamodeling Tumor-Immune System Interaction, Tumor Evasion and Immunotherapy. Math Comput Modell. 2008, 47: 614-637. 10.1016/j.mcm.2007.02.032.View ArticleGoogle Scholar
- Kronick N, Kogan Y, Elishmereni M, Halevi-Tobias K, Vuk Pavlovic S, Agur Z: Predicting Effect of Prostate Cancer Immunotherapy By Personalized Mathematical Models. PLoS One. 2010, 5: e15482-10.1371/journal.pone.0015482.View ArticleGoogle Scholar
- Kim PS, Lee PP, Levy D: Dynamics and Potential Impact of the Immune Response to Chronic Myelogenous Leukemia. PLoS Computuational Biol. 2008, 4: e1000095-10.1371/journal.pcbi.1000095.View ArticleGoogle Scholar
- Horsthemke W, Lefever R: Phase Transitions Induced By External Noise. Phys Lett A. 1977, 64: 19-21. 10.1016/0375-9601(77)90512-6.View ArticleGoogle Scholar
- Du LC, Mei DC: The Critical Phenomenon and the Re-entrance Phenomenon in the Anti-Tumor Model Induced By the Time Delay. Phys Lett A. 2010, 374: 3275-3279. 10.1016/j.physleta.2010.06.003.View ArticleGoogle Scholar
- Caravagna G, d’Onofrio A, Milazzo P, Barbuti R: Antitumour Immune Surveillance Through Stochastic Oscillations. J Theor Biol. 2010, 265 (3): 336-345. 10.1016/j.jtbi.2010.05.013.PubMedView ArticleGoogle Scholar
- d’Onofrio A: Tumour Evasion from Imune System Control as Bounded-Noise Induced Transition. Phys Rev E. 2010, 81: art. n. 021923-Google Scholar
- Caravagna G, d’Onofrio A, Barbuti R: Fine-tuning anti-tumor immunotherapies via stochastic simulations. BMC Bioinf. 2012, 13 (Suppl 4): art.n. S8-View ArticleGoogle Scholar
- Bellomo N, Bellouquid A, Delitala M: MAthematical Topics on the Modelling Complex Multicellular Systems and Tumor Immune Cell Competition. Math Models and Methods in Appl Sci. 2004, 14: 1683-1733. 10.1142/S0218202504003799.View ArticleGoogle Scholar
- Bellomo N, Delitala M: From the Mathematical Kinetic, and Stochastic Game Theory for Active Particles to Modelling Mutations, Onset, Progression and Immune Competition of Cancer Cells. Phys Life Rev. 2008, 5: 183-206. 10.1016/j.plrev.2008.07.001.View ArticleGoogle Scholar
- Owen M, Sherratt J: Pattern Formation and Spatiotemporal Irregularity in a Model for Macrophage-Tumour Interactions. J Theor Biol. 1997, 189: 63-80. 10.1006/jtbi.1997.0494.PubMedView ArticleGoogle Scholar
- Owen M, Sherratt J: Mathematical Modelling of Macrophage Dynamics in Tumours. Math Methods and Models in Appl Sci. 1999, 9: 513-539. 10.1142/S0218202599000270.View ArticleGoogle Scholar
- d’Onofrio A: Spatiotemporal effects of a possible chemorepulsion of tumor cells by immune system effectors. J Theoreoretical Biol. 2012, 296: 41-48.View ArticleGoogle Scholar
- Bellomo N: Modelling the Hiding–learning Dynamics in Large Living Systems. Appl Math Lett. 2010, 23: 907-911. 10.1016/j.aml.2010.04.008.View ArticleGoogle Scholar
- Galon J, Costes A, Sanchez-Cabo F, Kirilovsky A, Mlecnik B, Lagorce-Pagès C, Tosolini M, Camus M, Berger A, Wind P, Zinzindohoué F, Bruneval P, Cugnenc PH, Trajanoski Z, Fridman WH, Pagès F: Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Sci. 2006, 313: 1960-1964. 10.1126/science.1129139.View ArticleGoogle Scholar
- Vianello F, Papeta N, Chen T, Kraft P, White N, Hart WK, Kircher MF, Swart E, Rhee S, Irimia D, Toner M, Weissleder R, Poznansky MC, Palù G: Murine B16 melanomas expressing high levels of the chemokine stromal-derived factor-1/CXCL12 induce tumor-specific T cell chemorepulsion and escape from immune control. J Immunol. 2006, 176: 2902-2914.PubMedView ArticleGoogle Scholar
- Kim R, Emi M, Tanabe K: Cancer Immunnoediting from Immune Surveillance to Immune Escape. Immunol. 2007, 121: 1-14. 10.1111/j.1365-2567.2007.02587.x.View ArticleGoogle Scholar
- Kuznetsov V: A mathematical model for the interaction between cytotoxic lymphocytes and tumor cells. Analysis of the growth, stabilization and regression of the B cell lymphoma in mice chimeric with respect to the major histocompatibility complex. Biomed Sci. 1991, 2: 465-476.PubMedGoogle Scholar
- Zitvogel L, Tesniere A, Kroemer G: Cancer Despite Immunosurveillance: Immunoselection and Immunosubversion. Nat Rev Immunol. 2006, 6: 715-727. 10.1038/nri1936.PubMedView ArticleGoogle Scholar
- Koonin EV, Wolf YI: Is evolution Darwinian and/or Lamarckian?. Biol Direct. 2009, 4: 42-10.1186/1745-6150-4-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Danchin E, Charmantier A, Champagne F, Mesoudi A, Pujol B, Blanchet S: Beyond DNA: integrating inclusive inheritance into an extended theory of evolution. Nat Rev Genet. 2011, 12: 475-486. 10.1038/nrg3028.PubMedView ArticleGoogle Scholar
- Rando OJ, Verstrepen KJ: Timescales of Genetic and Epigenetic inheritance. Cell. 2007, 128: 655-668. 10.1016/j.cell.2007.01.023.PubMedView ArticleGoogle Scholar
- Bindra RS, Glazer PM: Genetic Instability and the tumour microenviroment towards the concept of microenviroment-induce mutagenesis. Mutation Res. 2005, 569: 75-85. 10.1016/j.mrfmmm.2004.03.013.PubMedView ArticleGoogle Scholar
- Huang LE, Bindra RS, Glazer PM, Harris AL: Hypoxia-induced genetic instability - a calculated mechanism underlying tumour progression. J Mol Med. 2007, 85: 139-148. 10.1007/s00109-006-0133-6.PubMedView ArticleGoogle Scholar
- Handel AE, Ramagopalan SV: Is Larmarckian evolution relevant to medicine?. BMC Med Genet. 2010, 11: 73-[[http://dx.doi.org/10.1186/1471-2350-11-73 ]]PubMedPubMed CentralView ArticleGoogle Scholar
- Uhr J, Marches R: Dormancy in a Model of Murine B Cell Lymphoma. Seminars in Cancer Biol. 2001, 11: 277-283. 10.1006/scbi.2001.0383.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.