Raphaël Forien

I study spatio-temporal stochastic processes arising in population genetics, evolutionary ecology. I try to describe the limiting behaviour and the asymptotic fluctuations of these spatio-temporal processes in different regimes. At the moment, my research focuses on three main aspects.

Spatial models in population genetics

During my PhD, I started to work on the spatial Lambda-Fleming-Viot process, developed by N. Barton, A. Etheridge and A. Véber. This model describes the evolution of the genetic composition of a spatially structured population, taking into account dispersal, mutations, natural selection and the randomness of reproduction events (called genetic drift), for a population evolving in any spatial dimension.

Evolution of neutral allele frequency in the two dimensional SLFV

For example, the figure above shows the evolution of a single (neutral) allele frequency in a two-dimensional space, starting from a uniform 1/2 frequency. After a few reproduction events, the frequency has increased in some regions, and decreased in others, and after many reproduction events, we see some regions where the frequency is close to 1 and other regions where it is close to 0, separated by very wiggly interfaces where the frequency is between 0 and 1.

My work focuses on finding scaling limits of this process (and of the related stepping-stone model in discrete space) in various settings: natural selection, heterogeneous dispersal, with obstacles, etc.

More recently, with Alison Etheridge and Sarah Penington, I have started to study the expansion of an allele under bistable frequency-dependent selection.

Related publications

Adaptation of trait-structured populations under environmental stress

As a part of the ANR project RESISTE, coordinated by Lionel Roques and Guillaume Martin, I study stochastic and deterministic models describing the evolution of the phenotype distribution of adapting populations under environmental stress. This stress can take the form of a moving environment, where the phenotype conferring the maximal fitness gradually moves away from its current position, or of a brutal change in the fitness landscape. In the latter case, some populations are unable to adapt to the new environment, and become extinct, while others are able to evolve in order to survive. This is called evolutionary rescue and is at the centre of the ANR RESISTE project.

Related publications

General epidemic models

The classical SIR epidemic model describes the evolution of the proportions of susceptible (S), infectious (I) and revmoved (R) individuals in a large population during a major epidemic. According to this model, the proportions S(t), I(t) and R(t) solve the following system of ordinary differential equations:

dS(t)/dt = -lambda S(t) I(t) dI(t)/dt = lambda S(t) I(t) - gamma I(t) dR(t)/dt = gamma I(t)

These equation describe the large population limit of a stochastic epidemic model under some crucial assumptions: the population is homogeneous, and infectious contacts take place between uniformly sampled pairs of individuals, and infected individuals remain infectious (with constant infectivity lambda) during i.i.d. exponentially distributed periods (with parameter gamma). In collaboration with Etienne Pardoux and Guodong Pang, I study generalisations of this model which relax one or several of these assumptions.

In particular, we are interested in non-Markovian epidemic models, where infected individuals remain infectious during arbitrarily distributed periods, and where their infectivity is given by a random function of time of the form lambda(t-tau), where tau is the time at which the individual is infected. The large population limit of such a system then takes the form of a system of integral equations, which are non-local in time (also called distributed delay equations).

Related publications