Asymptotic behavior of an epidemic model with infinitely many variants

We investigate the long-time dynamics of a SIR epidemic model with infinitely many pathogen variants infecting a homogeneous host population. We show that the basic reproduction number $\mathcal{R}_0$ of the pathogen can be defined in that case and corresponds to a threshold between the persistence ($\mathcal{R}_0>1$) and the extinction ($\mathcal{R}_0\leq 1$) of the pathogen. When $\mathcal{R}_0>1$ and the maximal fitness is attained by at least one variant, we show that the systems reaches an equilibrium state that can be explicitly determined from the initial data. When $\mathcal{R}_0>1$ but none of the variants attain the maximal fitness, the situation is more intricate. We show that, in general, the pathogen is uniformly persistent and any family of variants that have a fitness which is uniformly lower than the optimal fitness, eventually gets extinct. We derive a condition under which the total pathogen population converges to a limit which can be computed explicitly. We also find counterexamples that show that, when our condition is not met, the total pathogen population may converge to an unexpected value, or the system can even reach an eternally transient behavior where the total pathogen population between several values. We illustrate our results with numerical simulations that emphasize the wide variety of possible dynamics.

[1]  J. Demongeot,et al.  A Kermack–McKendrick model with age of infection starting from a single or multiple cohorts of infected patients , 2022, Proceedings of the Royal Society A.

[2]  H. Brüssow COVID‐19: Omicron – the latest, the least virulent, but probably not the last variant of concern of SARS‐CoV‐2 , 2022, Microbial biotechnology.

[3]  Zhihua Liu,et al.  Differential Equations and Population Dynamics I , 2022, Lecture Notes on Mathematical Modelling in the Life Sciences.

[4]  J. Burie,et al.  Asymptotic and transient behaviour for a nonlocal problem arising in population genetics , 2018, European Journal of Applied Mathematics.

[5]  Q. Griette Singular measure traveling waves in an epidemiological model with continuous phenotypes , 2017, Transactions of the American Mathematical Society.

[6]  Avraham Adler,et al.  Lambert-W Function , 2015 .

[7]  B. Perthame,et al.  Long-term behaviour of phenotypically structured models , 2014, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[8]  B. Perthame,et al.  Invasion fronts with variable motility: Phenotype selection, spatial sorting and wave acceleration , 2012, 1207.2355.

[9]  P. Jabin,et al.  On selection dynamics for competitive interactions , 2011, Journal of mathematical biology.

[10]  Horst R. Thieme,et al.  Global stability of the endemic equilibrium in infinite dimension: Lyapunov functions and positive operators , 2011 .

[11]  Gaël Raoul,et al.  Long Time Evolution of Populations under Selection and Vanishing Mutations , 2011 .

[12]  G. Barles,et al.  Concentration in Lotka-Volterra Parabolic or Integral Equations: A General Convergence Result , 2009, 0903.4952.

[13]  N. Mideo,et al.  Virulence evolution and the trade‐off hypothesis: history, current state of affairs and the future , 2009, Journal of evolutionary biology.

[14]  S. Mischler,et al.  On selection dynamics for continuous structured populations , 2008 .

[15]  S. Gandon,et al.  Applying population-genetic models in theoretical evolutionary epidemiology. , 2007, Ecology letters.

[16]  J. K. Hunter,et al.  Measure Theory , 2007 .

[17]  Xiao-Qiang Zhao,et al.  Global Attractors and Steady States for Uniformly Persistent Dynamical Systems , 2005, SIAM J. Math. Anal..

[18]  M. J.,et al.  Two Species Competition in a Periodic Environment * , 2004 .

[19]  Azmy S Ackleh,et al.  Competitive exclusion and coexistence for pathogens in an epidemic model with variable population size , 2003, Journal of mathematical biology.

[20]  Pierre Magal Mutation and recombination in a model of phenotype evolution , 2002 .

[21]  G. Webb,et al.  Mutation, selection, and recombination in a model of phenotype evolution , 1999 .

[22]  Bingtuan Li,et al.  Global Asymptotic Behavior of the Chemostat: General Response Functions and Different Removal Rates , 1998, SIAM J. Appl. Math..

[23]  Gail S. K. Wolkowicz,et al.  Global Asymptotic Behavior of a Chemostat Model with Discrete Delays , 1997, SIAM J. Appl. Math..

[24]  Sze-Bi Hsu,et al.  Competitive exclusion and coexistence for competitive systems on ordered Banach spaces , 1996 .

[25]  Gail S. K. Wolkowicz,et al.  Global dynamics of a mathematical model of competition in the chemostat: general response functions and differential death rates , 1992 .

[26]  Gail S. K. Wolkowicz,et al.  A MATHEMATICAL MODEL OF THE CHEMOSTAT WITH A GENERAL CLASS OF FUNCTIONS DESCRIBING NUTRIENT UPTAKE , 1985 .

[27]  P. Ewald,et al.  Host-Parasite Relations, Vectors, and the Evolution of Disease Severity , 1983 .

[28]  Hal L. Smith Competitive Coexistence in an Oscillating Chemostat , 1981 .

[29]  Sze-Bi Hsu,et al.  Limiting Behavior for Competing Species , 1978 .

[30]  Sze-Bi Hsu,et al.  A Mathematical Theory for Single-Nutrient Competition in Continuous Cultures of Micro-Organisms , 1977 .

[31]  J. P. Lasalle Some Extensions of Liapunov's Second Method , 1960 .

[32]  G. Tullock,et al.  Competitive Exclusion. , 1960, Science.

[33]  G. Hardin The competitive exclusion principle. , 1960, Science.

[34]  W. O. Kermack,et al.  A contribution to the mathematical theory of epidemics , 1927 .