No Arabic abstract
We study a mathematical model describing the growth process of a population structured by age and a phenotypical trait, subject to aging, competition between individuals and rare mutations. Our goals are to describe the asymptotic behaviour of the solution to a renewal type equation, and then to derive properties that illustrate the adaptive dynamics of such a population. We begin with a simplified model by discarding the effect of mutations, which allows us to introduce the main ideas and state the full result. Then we discuss the general model and its limitations. Our approach uses the eigenelements of a formal limiting operator, that depend on the structuring variables of the model and define an effective fitness. Then we introduce a new method which reduces the convergence proof to entropy estimates rather than estimates on the constrained Hamilton-Jacobi equation. Numerical tests illustrate the theory and show the selection of a fittest trait according to the effective fitness. For the problem with mutations, an unusual Hamiltonian arises with an exponential growth, for which we show existence of a global viscosity solution, using an uncommon a priori estimate and a new uniqueness result.
We study an equation structured by age and a phenotypic trait describing the growth process of a population subject to aging, competition between individuals, and mutations. This leads to a renewal equation which occurs in many evolutionary biology problems. We aim to describe precisely the asymp-totic behavior of the solution, to infer properties that illustrate the concentration and adaptive dynamics of such a population. This work is a continuation of [38] where the case without mutations is considered. When mutations are taken into account, it is necessary to control the corrector which is the main novelty of the present paper. Our approach consists in defining, by the Hopf transform, a Hamilton-Jacobi equation with an effective Hamiltonian as in homogenization problems. Its solution carries the singular part of the limiting density (typically Dirac masses) and the corrector defines the weights. The main new result of this paper is to prove that the corrector is uniformly bounded, using only the global Lipschitz and semi-convexity estimates for the viscosity solution of the Hamilton-Jacobi equation. We also establish the limiting equation satisfied by the corrector. To the best of our knowledge, this is the first example where such bounds can be proved in such a context.
We introduce and analyze several aspects of a new model for cell differentiation. It assumes that differentiation of progenitor cells is a continuous process. From the mathematical point of view, it is based on partial differential equations of transport type. Specifically, it consists of a structured population equation with a nonlinear feedback loop. This models the signaling process due to cytokines, which regulate the differentiation and proliferation process. We compare the continuous model to its discrete counterpart, a multi-compartmental model of a discrete collection of cell subpopulations recently proposed by Marciniak-Czochra et al. in 2009 to investigate the dynamics of the hematopoietic system. We obtain uniform bounds for the solutions, characterize steady state solutions, and analyze their linearized stability. We show how persistence or extinction might occur according to values of parameters that characterize the stem cells self-renewal. We also perform numerical simulations and discuss the qualitative behavior of the continuous model vis a vis the discrete one.
We study the mathematical properties of a general model of cell division structured with several internal variables. We begin with a simpler and specific model with two variables, we solve the eigenvalue problem with strong or weak assumptions, and deduce from it the long-time convergence. The main difficulty comes from natural degeneracy of birth terms that we overcome with a regularization technique. We then extend the results to the case with several parameters and recall the link between this simplified model and the one presented in cite{CBBP1}; an application to the non-linear problem is also given, leading to robust subpolynomial growth of the total population.
A square lattice is introduced into the Penna model for biological aging in order to study the evolution of diploid sexual populations under certain conditions when one single locus in the individuals genome is considered as identifier of species. The simulation results show, after several generations, the flourishing and coexistence of two separate species in the same environment, i.e., one original species splits up into two on the same territory (sympatric speciation). As well, the mortalities obtained are in a good agreement with the Gompertz law of exponential increase of mortality with age.
We study an S--I type epidemic model in an age-structured population, with mortality due to the disease. A threshold quantity is found that controls the stability of the disease-free equilibrium and guarantees the existence of an endemic equilibrium. We obtain conditions on the age-dependence of the susceptibility to infection that imply the uniqueness of the endemic equilibrium. An example with two endemic equilibria is shown. Finally, we analyse numerically how the stability of the endemic equilibrium is affected by extra-mortality and by the possible periodicities induced by the demographic age-structure.