BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sylvie Méléard

Sylvie Méléard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.PRPopulations and Evolutionmath.APApplicationsmath.DS

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Filling the gap between individual-based evolutionary models and Hamilton-Jacobi equations

We consider a stochastic model for the evolution of a discrete population structured by a trait with values on a finite grid of the torus, and with mutation and selection. Traits are vertically inherited unless a mutation occurs, and influence the birth and death rates. We focus on a parameter scaling where population is large, individual mutations are small but not rare, and the grid mesh for the trait values is much smaller than the size of mutation steps. When considering the evolution of the population in a long time scale, the contribution of small sub-populations may strongly influence the dynamics. Our main result quantifies the asymptotic dynamics of sub-population sizes on a logarithmic scale. We establish that under the parameter scaling the logarithm of the stochastic population size process, conveniently normalized, converges to the unique viscosity solution of a Hamilton-Jacobi equation. Such Hamilton-Jacobi equations have already been derived from parabolic integro-differential equations and have been widely developed in the study of adaptation of quantitative traits. Our work provides a justification of this framework directly from a stochastic individual based model, leading to a better understanding of the results obtained within this approach. The proof makes use of almost sure maximum principles and careful controls of the martingale parts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Time reversal of spinal processes for linear and non-linear branching processes near stationarity

We consider a stochastic individual-based population model with competition, trait-structure affecting reproduction and survival, and changing environment. The changes of traits are described by jump processes, and the dynamics can be approximated in large population by a non-linear PDE with a non-local mutation operator. Using the fact that this PDE admits a non-trivial stationary solution, we can approximate the non-linear stochastic population process by a linear birth-death process where the interactions are frozen, as long as the population remains close to this equilibrium. This allows us to derive, when the population is large, the equation satisfied by the ancestral lineage of an individual uniformly sampled at a fixed time $T$, which is the path constituted of the traits of the ancestors of this individual in past times $t\leq T$. This process is a time inhomogeneous Markov process, but we show that the time reversal of this process possesses a very simple structure (e.g. time-homogeneous and independent of $T$). This extends recent results where the authors studied a similar model with a Laplacian operator but where the methods essentially relied on the Gaussian nature of the mutations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bacterial Metabolic Heterogeneity: from Stochastic to Deterministic Models

We revisit the modeling of the diauxic growth of a pure microorganism on two distinct sugars which was first described by Monod. Most available models are deterministic and make the assumption that all cells of the microbial ecosystem behave homogeneously with respect to both sugars, all consuming the first one and then switching to the second when the first is exhausted. We propose here a stochastic model which describes what is called "metabolic heterogeneity". It allows to consider small populations as in microfluidics as well as large populations where billions of individuals coexist in the medium in a batch or chemostat. We highlight the link between the stochastic model and the deterministic behavior in real large cultures using a large population approximation. Then the influence of model parameter values on model dynamics is studied, notably with respect to the lag-phase observed in real systems depending on the sugars on which the microorganism grows. It is shown that both metabolic parameters as well as initial conditions play a crucial role on system dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inference with selection, varying population size and evolving population structure: Application of ABC to a forward-backward coalescent process with interactions

Genetic data are often used to infer demographic history and changes or detect genes under selection. Inferential methods are commonly based on models making various strong assumptions: demography and population structures are supposed \textit{a priori} known, the evolution of the genetic composition of a population does not affect demography nor population structure, and there is no selection nor interaction between and within genetic strains. In this paper, we present a stochastic birth-death model with competitive interactions and asexual reproduction. We develop an inferential procedure for ecological, demographic and genetic parameters. We first show how genetic diversity and genealogies are related to birth and death rates, and to how individuals compete within and between strains. {This leads us to propose an original model of phylogenies, with trait structure and interactions, that allows multiple merging}. Second, we develop an Approximate Bayesian Computation framework to use our model for analyzing genetic data. We apply our procedure to simulated data from a toy model, and to real data by analyzing the genetic diversity of microsatellites on Y-chromosomes sampled from Central Asia human populations in order to test whether different social organizations show significantly different fertility.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaling limits of bisexual Galton-Watson processes

Bisexual Galton-Watson processes are discrete Markov chains where reproduction events are due to mating of males and females. Owing to this interaction, the standard branching property of Galton-Watson processes is lost. We prove tightness for conveniently rescaled bisexual Galton-Watson processes, based on recent techniques developed by Bansaye, Caballero and M{é}l{é}ard. We also identify the possible limits of these rescaled processes as solutions of a stochastic system, coupling two equations through singular coefficients in Poisson terms added to square roots as coefficients of Brownian motions. Under some additional integrability assumptions, pathwise uniqueness of this limiting system of stochastic differential equations and convergence of the rescaled processes are obtained. Two examples corresponding to mutual fidelity are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Horizontal gene transfer: numerical comparison between stochastic and deterministic approaches

Horizontal gene Transfer (HT) denotes the transmission of genetic material between two living organisms, while the vertical transmission refers to a DNA transfer from parents to their offspring. Consistent experimental evidence report that this phenomenon plays an essential role in the evolution of certain bacterias. In particular, HT is believed to be the main instrument of developing the antibiotic resistance. In this work, we consider several models which describe this phenomenon: a stochastic jump process (individual-based) and the deterministic nonlinear integrod-ifferential equation obtained as a limit for large populations. We also consider a Hamilton-Jacobi equation, obtained as a limit of the deterministic model under the assumption of small mutations. The goal of this paper is to compare these models with the help of numerical simulations. More specifically, our goal is to understand to which extent the Hamilton-Jacobi model reproduces the qualitative behavior of the stochastic model and the phenomenon of evolutionary rescue in particular.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non local Lotka-Volterra system with cross-diffusion in an heterogeneous medium

We introduce a stochastic individual model for the spatial behavior of an animal population of dispersive and competitive species, considering various kinds of biological effects, such as heterogeneity of environmental conditions, mutual attractive or repulsive interactions between individuals or competition between them for resources. As a consequence of the study of the large population limit, global existence of a nonnegative weak solution to a multidimensional parabolic strongly coupled model of competing species is proved. The main new feature of the corresponding integro-differential equation is the nonlocal nonlinearity appearing in the diffusion terms, which may depend on the spatial densities of all population types. Moreover, the diffusion matrix is generally not strictly positive definite and the cross-diffusion effect allows for linearly growing influences of the subpopulations' sizes. We prove uniqueness of the finite measure-valued solution and give conditions under which the solution takes values in a functional space. We then make the competition kernels converge to a Dirac measure and obtain the existence of a solution to a locally competitive version of the previous equation. The large population limit is obtained by means of weak convergence tools for measure valued processes. The techniques employed in the study of the liming equation are novel, and are essentially based on the underlying stochastic flow related to the dispersive part of the dynamics, together with the use of suitable dual distances in the space of finite measures. This is an extended version of the paper "Non local Lotka-Volterra system with cross-diffusion in a heterogeneous medium" appeared in J. Math. Biol (2015) 70:829-854.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Speed of coming down from infinity for birth and death processes

We finely describe the speed of "coming down from infinity" for birth and death processes which eventually become extinct. Under general assumptions on the birth and death rates, we firstly determine the behavior of the successive hitting times of large integers. We put in light two different regimes depending on whether the mean time for the process to go from $n+1$ to $n$ is negligible or not compared to the mean time to reach $n$ from infinity. In the first regime, the coming down from infinity is very fast and the convergence is weak. In the second regime, the coming down from infinity is gradual and a law of large numbers and a central limit theorem for the hitting times sequence hold. By an inversion procedure, we deduce that the process is a.s. equivalent to a non-increasing function when the time goes to zero. Our results are illustrated by several examples including applications to population dynamics and population genetics. The particular case where the death rate varies regularly is studied in details.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of competition and horizontal trait inheritance on invasion, fixation and polymorphism

Horizontal transfer (HT) of heritable information or `traits' (carried by genetic elements, endosymbionts, or culture) is widespread among living organisms. Yet current ecological and evolutionary theory addressing HT is limited. We present a modeling framework for the dynamics of two populations that compete for resources and exchange horizontally (transfer) an otherwise vertically inherited trait. Competition influences individual demographics, affecting population size, which feeds back on the dynamics of transfer. We capture this feedback with a stochastic individual-based model, from which we derive a deterministic approximation for large populations. The interaction between horizontal transfer and competition makes possible the stable (or bi-stable) polymorphic maintenance of deleterious traits (including costly plasmids). When transfer rates are of a general density-dependent form, transfer stochasticity contributes strongly to population fluctuations. For an initially rare trait, we describe the probabilistic dynamics of invasion and fixation. Acceleration of fixation by HT is faster when competition is weak in the resident population. Thus, HT can have a major impact on the distribution of mutational effects that are fixed, and our model provides a basis for a general theory of the influence of HT on eco-evolutionary dynamics and adaptation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Singular limits for reaction-diffusion equations with fractional Laplacian and local or nonlocal nonlinearity

We perform an asymptotic analysis of models of population dynamics with a fractional Laplacian and local or nonlocal reaction terms. The first part of the paper is devoted to the long time/long range rescaling of the fractional Fisher-KPP equation. This rescaling is based on the exponential speed of propagation of the population. In particular we show that the only role of the fractional Laplacian in determining this speed is at the initial layer where it determines the thickness of the tails of the solutions. Next, we show that such rescaling is also possible for models with non-local reaction terms, as selection-mutation models. However, to obtain a more relevant qualitative behavior for this second case, we introduce, in the second part of the paper, a second rescaling where we assume that the diffusion steps are small. In this way, using a WKB ansatz, we obtain a Hamilton-Jacobi equation in the limit which describes the asymptotic dynamics of the solutions, similarly to the case of selection-mutation models with a classical Laplace term or an integral kernel with thin tails. However, the rescaling introduced here is very different from the latter cases. We extend these results to the multidimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptation in a stochastic multi-resources chemostat model

We are interested in modeling the Darwinian evolution resulting from the interplay of phenotypic variation and natural selection through ecological interactions, in the specific scales of the biological framework of adaptive dynamics. Adaptive dynamics so far has been put on a rigorous footing only for direct competition models (Lotka-Volterra models) involving a competition kernel which describes the competition pressure from one individual to another one. We extend this to a multi-resources chemostat model, where the competition between individuals results from the sharing of several resources which have their own dynamics. Starting from a stochastic birth and death process model, we prove that, when advantageous mutations are rare, the population behaves on the mutational time scale as a jump process moving between equilibrium states (the polymorphic evolution sequence of the adaptive dynamics literature). An essential technical ingredient is the study of the long time behavior of a chemostat multi-resources dynamical system. In the small mutational steps limit this process in turn gives rise to a differential equation in phenotype space called canonical equation of adaptive dynamics. From this canonical equation and still assuming small mutation steps, we prove a rigorous characterization of the evolutionary branching points.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How do birth and death processes come down from infinity?

We finely describe the "coming down from infinity" for birth and death processes which eventually become extinct. Our biological motivation is to study the decrease of regulated populations which are initially large. Under general assumptions on the birth and death rates, we describe the behavior of the hitting time of large integers. We let two regimes appear and derive an expression of the speed of coming down from infinity. In the case of death rates with regular variations, we also get a central limit theorem and the asymptotic probability of extinction in small times. Finally, we apply our results to birth and death processes in varying environment in whose the environment influences the competition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear historical superprocess approximations for population models with past dependence

We are interested in the evolving genealogy of a birth and death process with trait structure and ecological interactions. Traits are hereditarily transmitted from a parent to its offspring unless a mutation occurs. The dynamics may depend on the trait of the ancestors and on its past and allows interactions between individuals through their lineages. We define an interacting historical particle process describing the genealogies of the living individuals; it takes values in the space of point measures on an infinite dimensional càdlàg path space. This individual-based process can be approximated by a nonlinear historical superprocess, under the assumptions of large populations, small individuals and allometric demographies. Because of the interactions, the branching property fails and we use martingale problems and fine couplings between our population and independent branching particles. Our convergence theorem is illustrated by two examples of current interest in biology. The first one relates the biodiversity history of a population and its phylogeny, while the second treats a spatial model with competition between individuals through their past trajectories.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A host-parasite multilevel interacting process and continuous approximations

We are interested in modeling some two-level population dynamics, resulting from the interplay of ecological interactions and phenotypic variation of individuals (or hosts) and the evolution of cells (or parasites) of two types living in these individuals. The ecological parameters of the individual dynamics depend on the number of cells of each type contained by the individual and the cell dynamics depends on the trait of the invaded individual. Our models are rooted in the microscopic description of a random (discrete) population of individuals characterized by one or several adaptive traits and cells characterized by their type. The population is modeled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation and death for individuals and birth and death for cells. The interaction between individuals (resp. between cells) is described by a competition between individual traits (resp. between cell types). We look for tractable large population approximations. By combining various scalings on population size, birth and death rates and mutation step, the single microscopic model is shown to lead to contrasting nonlinear macroscopic limits of different nature: deterministic approximations, in the form of ordinary, integro- or partial differential equations, or probabilistic ones, like stochastic partial differential equations or superprocesses. The study of the long time behavior of these processes seems very hard and we only develop some simple cases enlightening the difficulties involved.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A rigorous model study of the adaptative dynamics of Mendelian diploids

Adaptive dynamics so far has been put on a rigorous footing only for clonal inheritance. We extend this to sexually reproducing diploids, although admittedly still under the restriction of an unstructured population with Lotka-Volterra-like dynamics and single locus genetics (as in Kimura's 1965 infinite allele model). We prove under the usual smoothness assumptions, starting from a stochastic birth and death process model, that, when advantageous mutations are rare and mutational steps are not too large, the population behaves on the mutational time scale (the 'long' time scale of the literature on the genetical foundations of ESS theory) as a jump process moving between homozygous states (the trait substitution sequence of the adaptive dynamics literature). Essential technical ingredients are a rigorous estimate for the probability of invasion in a dynamic diploid population, a rigorous, geometric singular perturbation theory based, invasion implies substitution theorem, and the use of the Skorohod $M_1$ topology to arrive at a functional convergence result. In the small mutational steps limit this process in turn gives rise to a differential equation in allele or in phenotype space of a type referred to in the adaptive dynamics literature as 'canonical equation'.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lévy flights in ecology

We are interested in modeling Darwinian evolution resulting from the interplay of phenotypic variation and natural selection through ecological interactions. The population is modeled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation, and death, as influenced by each individual's trait values, and interactions between individuals. An offspring usually inherits the trait values of her progenitor, except when a random mutation causes the offspring to take an instantaneous mutation step at birth to new trait values. In the case we are interested in, the probability distribution of mutations has a heavy tail and belongs to the domain of attraction of a stable law. We investigate the large-population limit with allometric demographies: larger populations made up of smaller individuals which reproduce and die faster, as is typical for micro-organisms. We show that depending on the allometry coefficient the limit behavior of the population process can be approximated by nonlinear Lévy flights of different nature: either deterministic, in the form of nonlocal fractional reaction-diffusion equations, or stochastic, as nonlinear super-processes with the underlying reaction and a fractional diffusion operator. These approximation results demonstrate the existence of such nontrivial fractional objects; their uniqueness is also proved.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quasi-stationary distributions and population processes

This survey concerns the study of quasi-stationary distributions with a specific focus on models derived from ecology and population dynamics. We are concerned with the long time behavior of different stochastic population size processes when 0 is an absorbing point almost surely attained by the process. The hitting time of this point, namely the extinction time, can be large compared to the physical time and the population size can fluctuate for large amount of time before extinction actually occurs. This phenomenon can be understood by the study of quasi-limiting distributions. In this paper, general results on quasi-stationarity are given and examples developed in detail. One shows in particular how this notion is related to the spectral properties of the semi-group of the process killed at 0. Then we study different stochastic population models including nonlinear terms modeling the regulation of the population. These models will take values in countable sets (as birth and death processes) or in continuous spaces (as logistic Feller diffusion processes or stochastic Lotka-Volterra processes). In all these situations we study in detail the quasi-stationarity properties. We also develop an algorithm based on Fleming-Viot particle systems and show a lot of numerical pictures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Slow and fast scales for superprocess limits of age-structured populations

A superprocess limit for an interacting birth-death particle system modelling a population with trait and physical age-structures is established. Traits of newborn offspring are inherited from the parents except when mutations occur, while ages are set to zero. Because of interactions between individuals, standard approaches based on the Laplace transform do not hold. We use a martingale problem approach and a separation of the slow (trait) and fast (age) scales. While the trait marginals converge in a pathwise sense to a superprocess, the age distributions, on another time scale, average to equilibria that depend on traits. The convergence of the whole process depending on trait and age, only holds for finite-dimensional time-marginals. We apply our results to the study of examples illustrating different cases of trade-off between competition and senescence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Uniform estimates for metastable transition times in a coupled bistable system

We consider a coupled bistable N-particle system driven by a Brownian noise, with a strong coupling corresponding to the synchronised regime. Our aim is to obtain sharp estimates on the metastable transition times between the two stable states, both for fixed N and in the limit when N tends to infinity, with error estimates uniform in N. These estimates are a main step towards a rigorous understanding of the metastable behavior of infinite dimensional systems, such as the stochastically perturbed Ginzburg-Landau equation. Our results are based on the potential theoretic approach to metastability.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Polymorphic evolution sequence and evolutionary branching

We are interested in the study of models describing the evolution of a polymorphic population with mutation and selection in the specific scales of the biological framework of adaptive dynamics. The population size is assumed to be large and the mutation rate small. We prove that under a good combination of these two scales, the population process is approximated in the long time scale of mutations by a Markov pure jump process describing the successive trait equilibria of the population. This process, which generalizes the so-called trait substitution sequence, is called polymorphic evolution sequence. Then we introduce a scaling of the size of mutations and we study the polymorphic evolution sequence in the limit of small mutations. From this study in the neighborhood of evolutionary singularities, we obtain a full mathematical justification of a heuristic criterion for the phenomenon of evolutionary branching. To this end we finely analyze the asymptotic behavior of 3-dimensional competitive Lotka-Volterra systems.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Invasion and adaptive evolution for individual-based spatially structured populations

The interplay between space and evolution is an important issue in population dynamics, that is in particular crucial in the emergence of polymorphism and spatial patterns. Recently, biological studies suggest that invasion and evolution are closely related. Here we model the interplay between space and evolution starting with an individual-based approach and show the important role of parameter scalings on clustering and invasion. We consider a stochastic discrete model with birth, death, competition, mutation and spatial diffusion, where all the parameters may depend both on the position and on the trait of individuals. The spatial motion is driven by a reflected diffusion in a bounded domain. The interaction is modelled as a trait competition between individuals within a given spatial interaction range. First, we give an algorithmic construction of the process. Next, we obtain large population approximations, as weak solutions of nonlinear reaction-diffusion equations with Neumann's boundary conditions. As the spatial interaction range is fixed, the nonlinearity is nonlocal. Then, we make the interaction range decrease to zero and prove the convergence to spatially localized nonlinear reaction-diffusion equations, with Neumann's boundary conditions. Finally, simulations based on the microscopic individual-based model are given, illustrating the strong effects of the spatial interaction range on the emergence of spatial and phenotypic diversity (clustering and polymorphism) and on the interplay between invasion and evolution. The simulations focus on the qualitative differences between local and nonlocal interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Estimates for the density of a nonlinear Landau process

The aim of this paper is to obtain estimates for the density of the law of a specific nonlinear diffusion process at any positive bounded time. This process is issued from kinetic theory and is called Landau process, by analogy with the associated deterministic Fokker-Planck-Landau equation. It is not Markovian, its coefficients are not bounded and the diffusion matrix is degenerate. Nevertheless, the specific form of the diffusion matrix and the nonlinearity imply the non-degeneracy of the Malliavin matrix and then the existence and smoothness of the density. In order to obtain a lower bound for the density, the known results do not apply. However, our approach follows the main idea consisting in discretizing the interval time and developing a recursive method. To this aim, we prove and use refined results on conditional Malliavin calculus. The lower bound implies the positivity of the solution of the Landau equation, and partially answers to an analytical conjecture. We also obtain an upper bound for the density, which again leads to an unusual estimate due to the bad behavior of the coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Individual-based probabilistic models of adaptive evolution and various scaling approximations

We are interested in modelling Darwinian evolution, resulting from the interplay of phenotypic variation and natural selection through ecological interactions. Our models are rooted in the microscopic, stochastic description of a population of discrete individuals characterized by one or several adaptive traits. The population is modelled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation, and death, as influenced by each individual's trait values, and interactions between individuals. An offspring usually inherits the trait values of her progenitor, except when a mutation causes the offspring to take an instantaneous mutation step at birth to new trait values. We look for tractable large population approximations. By combining various scalings on population size, birth and death rates, mutation rate, mutation step, or time, a single microscopic model is shown to lead to contrasting macroscopic limits, of different nature: deterministic, in the form of ordinary, integro-, or partial differential equations, or probabilistic, like stochastic partial differential equations or superprocesses. In the limit of rare mutations, we show that a possible approximation is a jump process, justifying rigorously the so-called trait substitution sequence. We thus unify different points of view concerning mutation-selection evolutionary models.

0 citations0 reviewsNo code linkedOpen access