Dynamic maximum entropy provides accurate approximation of structured population dynamics

Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.

[1]  T. Schoener The Newest Synthesis: Understanding the Interplay of Evolutionary and Ecological Dynamics , 2011, Science.

[2]  H. P. de Vladar,et al.  The statistical mechanics of a polygenic character under stabilizing selection, mutation and drift , 2010, Journal of The Royal Society Interface.

[3]  Lisa Borland,et al.  Ito-Langevin equations within generalized thermostatistics , 1998 .

[4]  L. A. Segel,et al.  The Quasi-Steady-State Assumption: A Case Study in Perturbation , 1989, SIAM Rev..

[5]  W. Bialek,et al.  Maximum entropy models for antibody diversity , 2009, Proceedings of the National Academy of Sciences.

[6]  C. Gardiner Handbook of Stochastic Methods , 1983 .

[7]  E. Jaynes The Minimum Entropy Production Principle , 1980 .

[8]  J. N. Thompson,et al.  Rapid evolution as an ecological process. , 1998, Trends in ecology & evolution.

[9]  Kingshuk Ghosh,et al.  The Maximum Caliber Variational Principle for Nonequilibria. , 2020, Annual review of physical chemistry.

[10]  J. Polechová Is the sky the limit? On the expansion threshold of a species’ range , 2017, bioRxiv.

[11]  Prügel-Bennett,et al.  Analysis of genetic algorithms using statistical mechanics. , 1994, Physical review letters.

[12]  C. DeWitt-Morette,et al.  Techniques and Applications of Path Integration , 1981 .

[13]  K. Šišková,et al.  Extension and Justification of Quasi-Steady-State Approximation for Reversible Bimolecular Binding , 2015, Bulletin of mathematical biology.

[14]  M. Rattray,et al.  Cumulant dynamics of a population under multiplicative selection, mutation, and drift. , 1999, Theoretical population biology.

[15]  J. Kingsolver,et al.  Evolution of plasticity and adaptive responses to climate change along climate gradients , 2017, Proceedings of the Royal Society B: Biological Sciences.

[16]  K. Dill,et al.  Principles of maximum entropy and maximum caliber in statistical physics , 2013 .

[17]  A. R. Plastino,et al.  Statistical treatment of autonomous systems with divergencelless flows , 1996 .

[18]  W. Ewens Mathematical Population Genetics : I. Theoretical Introduction , 2004 .

[19]  Polygenic local adaptation in metapopulations: A stochastic eco‐evolutionary model , 2021, Evolution; international journal of organic evolution.

[20]  Kingshuk Ghosh,et al.  Perspective: Maximum caliber is a general variational principle for dynamical systems. , 2017, The Journal of chemical physics.

[21]  Michael J. Berry,et al.  Weak pairwise correlations imply strongly correlated network states in a neural population , 2005, Nature.

[22]  Sebastian Walcher,et al.  Quasi-Steady State: Searching for and Utilizing Small Parameters , 2013 .

[23]  G. Tkačik,et al.  A General Approximation for the Dynamics of Quantitative Traits , 2015, Genetics.

[24]  A. Plastino,et al.  The nonlinear Fokker-Planck equation with state-dependent diffusion - a nonextensive maximum entropy approach , 1999 .

[25]  A. R. Plastino,et al.  MINIMUM KULLBACK ENTROPY APPROACH TO THE FOKKER-PLANCK EQUATION , 1997 .

[26]  H. P. de Vladar,et al.  The contribution of statistical physics to evolutionary biology. , 2011, Trends in ecology & evolution.

[27]  P. Markowich,et al.  Well posedness and maximum entropy approximation for the dynamics of quantitative traits , 2017, Physica D: Nonlinear Phenomena.

[28]  E. Jaynes Information Theory and Statistical Mechanics , 1957 .

[29]  S. Gueron,et al.  The Equilibrium Behavior of Reversible Coagulation-Fragmentation Processes , 1999 .

[30]  J. Pandolfi,et al.  Predicting evolutionary responses to climate change in the sea. , 2013, Ecology letters.

[31]  J. Hermisson,et al.  Evolutionary Rescue in Structured Populations , 2013, The American Naturalist.

[32]  W. Bialek,et al.  Statistical mechanics for natural flocks of birds , 2011, Proceedings of the National Academy of Sciences.

[33]  Tsallis,et al.  Anomalous diffusion in the presence of external forces: Exact time-dependent solutions and their thermostatistical basis. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[34]  Himani Sachdeva Effect of partial selfing and polygenic selection on establishment in a new habitat , 2019, bioRxiv.

[35]  H. Haken,et al.  Synergetics , 1988, IEEE Circuits and Devices Magazine.

[36]  H. P. de Vladar,et al.  Statistical Mechanics and the Evolution of Polygenic Quantitative Traits , 2009, Genetics.