Sylvain Billiard
email: sylvain.billiard@univ-lille1.fr

Regis Ferriere
email: ferriere@biologie.ens.fr

Sylvie Méléard
email: sylvie.meleard@polytechnique.edu

Viet Chi Tran

Régis Ferrière

Stochastic dynamics of adaptive trait and neutral marker driven by eco-evolutionary feedbacks

Keywords: MSC 2000 subject classification: 92D25, 60J80, 92D15, 60J75 Mutation-selection, measure-valued individual-based model, neutral diversity, hitchhiking, selective sweeps, adaptive dynamics, limit theorems for multi-scale processes, Substitution Fleming-Viot Process

de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

The science of biodiversity currently faces the challenge of understanding how ecological processes shape evolutionary change, and reciprocally how evolution affects the structure and function of ecological systems [START_REF] Schoener | The newest synthesis: Understanding the interplay of evolutionary and ecological dynamics[END_REF]. Such eco-evolutionary feedbacks determine the dynamics of so-called adaptive traits -quantitative characters that are heritable yet mutable from parent to offspring [START_REF] Dieckmann | The dynamical theory of coevolution: a derivation from stochastic ecological processes[END_REF][START_REF] Metz | Adaptative dynamics, a geometrical study of the consequences of nearly faithful reproduction[END_REF]. Under the combined assumptions of large population and rare mutation scalings, the time evolution of an adaptive trait can be described as a sequence of mutant invasions, each being driven by positive selection in the ecological context set by the 'resident' value of the adaptive trait [START_REF] Metz | How should we define 'fitness' for general ecological scenarios?[END_REF]. The resulting evolutionary model is a jump process called the Trait Substitution Sequence (TSS): every new mutant trait either goes extinct, or replaces the resident, causing the TSS to jump from the former resident population equilibrium to a new equilibrium [START_REF] Metz | Adaptative dynamics, a geometrical study of the consequences of nearly faithful reproduction[END_REF][START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF][START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF]. In population genetics, these jumps are known as selective sweeps [START_REF] Barton | The effect of hitchhiking on neutral genealogies[END_REF][START_REF] Stephan | The effect of strongly selected substitutions on neutral polymorphism: analytical results based on diffusion theory[END_REF]. Previous works support the view that the TSS as a model of long-term phenotypic evolution is relatively insensitive to the details of the genetic determination of the trait [START_REF] Lloyd | Genetic and phenotypic models of natural selection[END_REF][START_REF] Christiansen | Evolution and intraspecific exploitative competition I. one-locus theory for small additive gene effects[END_REF][START_REF] Hammerstein | Darwinian adaptation, population genetics and the streetcar theory of evolution[END_REF][START_REF] Weissing | Genetic versus phenotypic models of selection: can genetics be neglected in a long-term perspective?[END_REF][START_REF] Matessi | Optimization under frequency-dependent selection[END_REF][START_REF] Eshel | Long-term evolution, short-term evolution, and population genetic theory[END_REF].

Whereas eco-evolutionary feedbacks can result in variation of adaptive traits among populations (and even within populations when evolutionary branching occurs, [START_REF] Geritz | Evolutionary singular strategies and the adaptive growth and branching of the evolutionary tree[END_REF], much of the molecular diversity measured by population geneticists involve DNA sequences of no known adaptive value, i.e. selectively neutral. A neutral sequence that is physically linked in the genome to the sequence that codes for the adaptive trait is called a marker of that trait. A longstanding question in evolutionary theory is understanding how variation in such molecular markers evolves, and how patterns of neutral molecular evolution can be used to infer the history of trait mutation that have driven past adaptation.

When adaptive mutations are rare, adaptation proceeds as a series of selective sweeps: a trait mutation occurs while the population is monomorphic for the trait, and increases rapidly in frequency toward fixation. Following on from [START_REF] Kojima | Survival process of linked mutant genes[END_REF], [START_REF] Smith | The hitchhiking effect of a favourable gene[END_REF] pointed out that selective sweeps purge genetic variation at linked sites: a particular marker allele goes to fixation as a consequence of linkage with the selected allele, a phenomenon they dubbed the 'hitchhiking effect'. Maynard Smith and Haigh's deterministic model was revisited in a stochastic approach by [START_REF] Ohta | Effects of selected linked locus on heterozygosity of neutral alleles (hitchhiking effect)[END_REF]. These seminal studies of hitchhiking focused on the short-term dynamics of an interaction between two alleles at the locus under selection and two alleles at the neutral locus. Longterm dynamics were considered first by [START_REF] Kaplan | The hitchhiking effect revisited[END_REF] who developed a stochastic model for finite populations to describe the effect of recurrent hitchhiking. In order to describe stationary levels of nucleotide diversity at the marker locus, they used the infinite site model and a coalescent approach under the assumption of constant population size and constant selection coefficients. This has generated an abundant theoretical literature on modeling the impact of selection on neutral polymorphism [START_REF] Barton | Genetic hitchhiking[END_REF][START_REF] Etheridge | An approximate sampling formula under genetic hitchhiking[END_REF]Durrett and Schweinsberg, 2004, and references therein). Recent deterministic models have relaxed the assumption of constant selection either because of the presence of genetic backgrounds (e.g. assuming a quantitative trait, [START_REF] Chevin | Selective sweep at a quantitative trait locus in the presence of background genetic variation[END_REF] or in the case of a parasite, because of the complexity of the demographic events involved in the life cycle [START_REF] Schneider | An analytical model for genetic hitchhiking in the evolution of antimalarial drug resistance[END_REF]. All previous models assume constant population size and constant selection, or that the population size is independent of the selective value of the individuals.

In this article, our goal is to relax these key assumptions. Under general ecological scenarios, eco-evolutionary feedbacks operate: as the adaptive trait evolves, population size and selection co-vary. The eco-evolutionary process of adaptive trait and neutral marker dynamics requires a rigorous mathematical framework, the foundation of which we establish here. We start with a 'microscopic', individual-based model where individuals have two heritable characteristics: (i) an adaptive trait that influences their intrinsic demographic rates and ecological interactions, and (ii) a genetic marker that has no demographic or ecological effects, hence, is selectively neutral. This work focuses entirely on asexual populations and short genomic regions that remain perfectly linked to the loci under selection, neglecting recombination. The population is described by a measure according to which each individual is represented by a Dirac mass that weights its characters. This leads to study the population eco-evolutionary dynamics as a measure-valued stochastic process.

The dynamics are driven by competition between individuals, asexual reproduction without or with mutation, and death. Variation in population size and selection as the trait evolves are mediated by the demographic effects of change in the trait. These effects are expected to influence the generation and maintenance of neutral variation.

The effect of mutation on the marker can be continuous or discrete. Our framework thus encompasses a variety of conventional mutation models such as the two-allele model, the stepwise mutation model, and the continuous state mutation model. Our distinctive assumption here is that the marker mutation process is much faster than the trait mutation process but much slower than the ecological time-scale of birth and death events. This is supported by the fact that most mutations are neutral or nearly neutral (such as mutations involved in microsatellite variation). Therefore, there are three time scales in the model: the fast ecological time scale of birth and death events, the slow time scale of trait mutation, and an intermediate time scale of marker mutation. We study the joint process of trait and marker dynamics on the trait mutation time scale.

We are interested in limit theorems when the population carrying capacity goes to infinity. Then, the population size stabilizes in a neighborhood of the ecological equilibrium and jumps to another equilibrium when a successful trait mutant goes to fixation in the population. This is the TSS dynamics of the adaptive trait. It does not depend on the marker and has been mathematically proved by [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF]. The novelty in the model and in the proofs comes from the time-scales difference for the marker and trait mutations. The study of the marker distribution during the invasion period requires careful consideration of the individual process and of the different scales involved. In a first period, starting with the single invading mutant, we prove that the marker distribution remains close to a Dirac mass at the value of the initial mutant. Until the next jump of the TSS, the marker evolves as a stochastic distribution-valued process. In the case where the marker mutation effects are continuous and small, this is a Fleming-Viot process whose drift and covariance depend on the resident adaptive trait. In every cases, for any marker mutation model, the collated dynamics define a measure-valued diffusive process with jumps that we call Substitution Fleming-Viot Process (SFVP). The convergence of the microscopic process to the SFVP is shown both in the sense of finite dimensional distributions and in the sense of convergence of trait-marker-time measures, thus improving previous results of [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF].

From a biological standpoint, we recover the conventional hitchhiking phenomenon: when a new mutant trait appears and sweeps through the population to fixation, the marker carried by the mutant individual is hitchhiked, and the marker distribution undergoes a genetical bottleneck. The mathematical construction of the SFVP process has new implications of biological relevance. Neutral diversity is restored after each adaptive jump, but as the adaptive trait evolves, population size, the mutation rate, genetic drift and demographic fluctuations change, which causes the rate of neutral polymorphism build-up and the moments of the marker distribution to change too. This suggests that the nature and structure of the whole eco-evolutionary feedback loop (i.e. how adaptive traits influence demographic rates and ecological interactions, and how ecological processes shape selection pressures on adaptive traits) may be important to explain the extreme disparities in genetic neutral diversity observed among species, even closely related ones and in the absence of differences in recombination profiles [START_REF] Cutter | Genomic signatures of selection at linked sites: unifying the disparity among species[END_REF]. In fact, it is well-known that demographic differences due to external causes (demographic bottleneck or population expansion due to environmental changes) can affect neutral diversity of a population and that closely related species can show very different neutral diversity patterns. Here, we show that internal causes of demographic variation involved in adaptation can also affect species differently.

The article is organized as follows. In Section 2, we start with the model description. The stochastic individual-based process and its key assumptions are carefully described and examples are provided. A key parameter is K, an integer that gives the order of the population size and is used to rescale the mutation rates and kernels. By letting K go to infinity we study the large population limit of the stochastic process. The main theorem is enounced and discussed in Section 3, where biological implications are also highlighted. Time scale separations implied by the dependence in K of the trait and marker mutations lead to homogenization phenomena and then to the SFVP. Our mathematical analysis provides a precise description of the genetical bottleneck that occurs at each trait substitution. We show that the marker of the initial mutant individual dominates in the marker distribution of the mutant population until this population reaches a neighborhood of the new ecological equilibrium. Then, we present two numerical examples based on an ecological model adapted from [START_REF] Dieckmann | On the origin of species by sympatric speciation[END_REF]. In the first example, marker mutation is described by a continuous state model that leads to a piecewise Fleming-Viot process (section 3.3) for the marker. In the second example, marker mutation follows a discrete two-allele model and the classical Wright-Fisher diffusion (3.5) is recovered. Further generalizations are discussed. The proof of the main theorem in the adaptive dynamics scaling is in Section 4. After having introduced a semi-martingale decomposition of our stochastic measure-valued process, we start with recalling and refining the result of [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] for the convergence of trait-marginals. For this purpose, we introduce the M1-topology on the Skorokhod space where the TSS lives, using some ideas of [START_REF] Collet | A rigorous model study of the adaptative dynamics of Mendelian diploids[END_REF]. The second part of the proof focuses on the marker distribution in an invading mutant population. This gives the result on the genetical bottleneck. Then, between two trait substitutions, the dynamics of the marker converges to a diffusive measure-valued process. As a conclusion of the proof, we show the convergence to the SFVP in the space of trait-marker-time measures.

The stochastic model

We consider an asexual population driven by births and deaths where each individual is characterized by hereditary types: a phenotypic trait under selection and a neutral marker. The trait and marker spaces X and U are assumed to be compact subsets of R. The type of individual i is thus a pair (x i , u i ), x i ∈ X being the trait value and u i ∈ U its neutral marker. The individual-based microscopic model from which we start is a stochastic birth and death process with density-dependence whose demographic parameters are functions of the trait under selection and are independent of the marker. We assume that the population size scales with an integer parameter K tending to infinity while individuals are weighted with 1 K . At any time t ≥ 0, we have a finite number N K t of individuals, each of them holding trait and marker values in X × U. Let us denote by ((x 1 , u 1 ), . . . , (x N K t , u N K t )) the trait and marker values of these individuals. The state of the population at time t ≥ 0, rescaled by K, is described by the point measure

ν K t = 1 K N K t i=1 δ (x i ,u i ) , (2.1)
where δ (x,u) is the Dirac measure at (x, u). This measure belongs to the set of finite point measures on X × U with mass 1/K. This set is a subset of the set M F (X × U) of finite measures on X × U, which is embedded with the weak convergence topology. We denote by ν, f the integral of the measurable function f with respect to the measure ν and by Supp(ν) the support of ν. Then ν K t , 1 = N K t K . For any t ≥ 0, we also introduce the trait marginal of the measure ν K t on X , denoted by X K t and defined by

X K t = 1 K N K t i=1 δ x i .
Therefore, the population measure ν K t writes

ν K t (dx, du) = X K t (dx) π K t (x, du) (2.2)
where π K t (x, du) is the marker distribution for a given trait value x defined by

π K t (x, du) = N K t i=1 1l x i =x δ u i N K t i=1 1l x i =x . (2.3)
Our purpose is to study the asymptotic behavior of the measure-valued process ν K at large times, when the trait and marker are inherited but mutations occur. The main interest of our model is that these mutations happen at different time scales for trait and marker, but both longer than the individuals lifetime scale. The trait mutates much slower than the marker and drives the evolution time scale. Thus, the limiting behavior results from the interplay of three time scales: births and deaths, trait mutations and marker mutations.

We describe the individuals' life history. The trait has an influence on the ability of individuals to survive (including competition with other ones) and to reproduce but the marker is neutral. The demographic parameters are thus functions of the trait only and are defined on X .

Assumption 2.1

• An individual with trait x and marker u reproduces with birth rate given by 0 ≤ b(x) ≤ b, the function b being continuous.

• Reproduction produces a single offspring which usually inherits the trait and marker of its ancestor except when a mutation occurs. Mutations on trait and marker occur independently with probabilities p K and q K respectively. Mutations are rare and the marker mutates much more often than the trait. We assume that

q K = p K r K , with p K = 1 K 2 , q K (log K) 2 → K→∞ 0 , r K → K→∞ +∞. (2.4)
• When a trait mutation occurs, the new trait of the descendant is x + k ∈ X with k chosen according to the probability measure m(x, k)dk.

• When a marker mutation occurs, the new marker of the descendant is u + h ∈ U with h chosen according to the probability measure G K (u, dh).

For any u ∈ U, G K (u, .) is approximated as follows when K tends to infinity:

lim K→+∞ sup u∈U r K K U (φ(u + h) -φ(u))G K (u, dh) -Aφ(u) = 0, (2.5)
where (A, D(A)) is the generator of a Feller semigroup and φ ∈ D(A) ⊆ C b (U , R), the set of continuous bounded real functions on U .

• An individual with trait x and marker u dies with intrinsic death rate 0 ≤ d(x) ≤ d, the function d being continuous. Moreover the individual experiences competition the effect of which is an additional death rate η

(x) C * ν K t (x) = η(x) K N K t i=1 C(x -x i ).
The quantity C(x -x i ) describes the competition pressure exerted by an individual with trait x i on an individual with trait x. We assume that the functions C and η are continuous and that there exists η > 0 such that ∀x, y ∈ X , η(x) C(x -y) ≥ η > 0.

(2.6)

A classical choice of competition function is C ≡ 1 which is called "mean field case" or "logistic case". In that case the competition death rate is η(x)N K t /K.

Remark 2.2 Let us insist on the generality of Assumption (2.5) which allows a larger set of possible dynamics.

• Equation (2.5) is for example true for

U = [u 1 , u 2 ], G K a centered Gaussian law (conditioned to U ) with variance σ K → 0 such that lim K σ 2 K r K K = σ 2 and Aφ = σ 2 2 φ ′′ for φ ∈ C 2 with φ ′ (u 1 ) = φ ′ (u 2 ) = 0. Choosing for example r K = K 3/2 , q K = 1/ √ K and σ 2 K = 1/ √ K works.
This choice can be seen as a continuous state generalization of the stepwise mutation model [START_REF] Ohta | Model of mutation appropriate to estimate number of electrophoretically detectable alleles in a finite population[END_REF].

• If in addition the distribution G K has a non zero mean µ K such that r K µ K K → µ > 0 corresponding to a mutational directional drift, then the operator A will be defined by Aφ = σ 2 2 φ ′′ + µφ ′ .

• If we relax the compactness of U and assume that U = R, a third example consists in taking for G K the law of a Pareto variable with index α ∈ (1, 2) divided by

K η/α , for η ∈ (0, 1].
Then it has been proved by [START_REF] Jourdain | Lévy flights in evolutionary ecology[END_REF] that

lim K sup u K η R (φ(u + h) -φ(u))G K (u, dh) - α 2 D α φ(u) = 0, where D α φ(u) = R (φ(u + h) -φ(u) -hφ ′ (u)1l |h|≤1 ) dh |h| 1+α
is the fractional Laplacian with index α. Thus if we take r K such that r K K 1+η converges as K tends to infinity, and choose A = D α in (2.5), Assumptions (2.4)-(2.5) will be satisfied as soon as η < 1.

• Another very interesting case is the discrete case when U = {a, A} is a set of two alleles. The mutation kernel is given by

G K (u, dv) = 1l u=a q a δ A (dv) + 1l u=A q A δ a (dv).
(2.7)

In this case, (2.5) implies that r K /K has a limit when K → +∞. Let r be this limit, then

Aφ(u) = r 1l u=a q a φ(A) -φ(a) + 1l u=A q A φ(a) -φ(A) .
(2.8)

We see that the ratio between the two mutation probabilities r K = q K /p K that allows convergence is highly dependent on the mutation distribution.

Note that since the demographic rates do not depend on the marker, the dynamics of the population distribution of the trait is independent of the marker distribution. But the dynamics of the marker distribution cannot be separated from the trait distribution as we shall see.

The process (ν K t , t ≥ 0) is a càdlàg M F (X × U)-valued Markov process. Existence and uniqueness in law of the process can be adapted from [START_REF] Fournier | A microscopic probabilistic description of a locally regulated population and macroscopic approximations[END_REF] and [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF] under the assumption that E( ν K 0 , 1 ) < +∞. Moreover, Assumption (2.6) allows to prove as by Champagnat citet

[Lemma 1]champag- nat06 that if for p ≥ 1, sup K∈N * E( ν K 0 , 1 p ) < +∞, then sup t∈R + ,K∈N * E ν K t , 1 p < +∞ (2.9)
which will be useful to study the tightness and convergence of the sequence.

Convergence to the Substitution Fleming-Viot Process

The adaptive trait mutation time scale is the slowest, equal to 1 Kp K = K by Assumptions (2.4). It scales the evolutionary time. So we shall consider the limiting behavior of (ν K Kt , t ≥ 0). We will see in section 4.2.1 that p K of order 1/K 2 is the only choice which leads to a non-trivial or non-degenerate marker dynamics.

Before stating our main result, we introduce several important ingredients which are used to describe the limit of (ν K Kt , t ≥ 0) when K → +∞. We conclude the section with extensions and simulations.

Invasion fitness function

The large population behavior of the process (ν K t , t ≥ 0) as K tends to infinity, can be studied by classical arguments and is given in the appendix. At the ecological time scale (of order 1), no mutation occurs in the asymptotic K → +∞. If the initial population has a single adaptive trait x, then, in the limit K → +∞, the trait distribution remains δ x since p K and q K vanish in the limit. The rescaled population size process (N K t /K, t ≥ 0) converges to the solution (n t , t ≥ 0) of the ordinary differential equation

dn t dt = b(x) -d(x) -η(x)C(0)n t n t (3.1)
which converges when t tends to infinity to the equilibrium

n x = b(x) -d(x) η(x)C(0) . (3.2)
Conversely, at the adaptive trait-mutation time scale Kt, new mutant traits can invade.

If they replace the previous traits, then the corresponding event is called "fixation". The probability of fixation of a mutant trait y in a trait resident population x at equilibrium depends on the invasion fitness function f (y; x):

f (y; x) = b(y) -d(y) -η(y) C(y -x) n x . (3.3)
This fitness function describes the initial growth of the mutant population. It does not depend on the neutral marker.

By simplicity we work under the assumption of 'invasion implies fixation', but this assumption will be relaxed in Section 3.5. When a mutant trait appears, either its line of descent replaces the resident population or it disappears. As a consequence, two traits can not coexist in the long term.

Assumption 3.1 ("Invasion implies fixation") For all x ∈ X and for almost every

y ∈ X , either b(y) -d(y) η(y)C(y -x) < b(x) -d(x) η(x)C(0) , or b(y) -d(y) η(y)C(y -x) > b(x) -d(x) η(x)C(0) and b(x) -d(x) η(x)C(x -y) < b(y) -d(y) η(y)C(0) .
Remark 3.2 In the case of logistic populations with C ≡ 1, this assumption is satisfied as soon as x → n x is strictly monotonous.

Main theorem

Let us first give the definition of the Fleming-Viot process which will appear in our setting (see e.g. [START_REF] Dawson | Wandering random measures in the Fleming-Viot model[END_REF][START_REF] Dawson | Mesure-valued markov processes[END_REF][START_REF] Donnelly | A countable representation of the Fleming-Viot measure-valued diffusion[END_REF][START_REF] Etheridge | An introduction to superprocesses[END_REF]. We recall that the operator A has been introduced in (2.5).

In the sequel, we denote by P(U ) and P(X ×U ) the probability measure spaces respectively on U and on X × U.

Definition 3.3 Let us fix x ∈ X and u ∈ U. The Fleming-Viot process (F u t (x, .), t ≥ 0) indexed by x, started at time 0 with initial condition δ u and associated with the mutation operator A is the P(U )-valued process whose law is characterized as the unique solution of the following martingale problem. For any φ ∈ D(A),

M x t (φ) = F u t (x, .), φ -φ(u) -b(x) t 0 F u s (x, .), Aφ ds (3.4)
is a continuous square integrable martingale with quadratic variation process

M x (φ) t = b(x) + d(x) + η(x)C(0) n x n x t 0 F u s (x, .), φ 2 -F u s (x, .), φ 2 ds = 2b(x) n x t 0 F u s (x, .), φ 2 -F u s (x, .), φ 2 ds. (3.5)
Let us now state our main theorem that describes the slow-fast dynamics of adaptive traits and neutral markers at the (trait) evolutionary time scale.

Theorem 3.4 We work under Assumptions 2.1 and 3.1. The initial conditions are ν K 0 (dy, dv) = n K 0 δ (x 0 ,u 0 ) (dy, dv) with lim K→∞ n K 0 = n x 0 and for any ǫ > 0, sup

K∈N * E((n K 0 ) 2+ǫ ) < +∞.
Then, the population process (ν K Kt , t ≥ 0) converges in law to the M F (X × U ))-valued process (V t (dy, dv), t ≥ 0) defined by

V t (dy, dv) = n Yt δ Yt (dy) F Ut t (Y t , dv), (3.6)
where the process

((Y t , U t ), t ≥ 0) on X × U, started at (x 0 , u 0 ), jumps at time t from (x, u) to (x + k, v) with the jump measure b(x) n x [f (x + k; x)] + b(x + k) F u t (x, dv) m(x, k)dk. (3.7)
The convergence holds in the sense of finite dimensional distributions on M F (X × U).

In addition, the convergence also holds in the space of trait-marker-time measures, i.e. the measure ν K Kt (dy, dv)dt on X ×U ×[0, T ] converges weakly to the measure n Yt δ Yt (dy) F Ut t (Y t , dv)dt for any T > 0.

Definition 3.5 The limiting measure-valued process (V t (dy, dv), t ≥ 0) is called Substitution Fleming-Viot Process. It generalizes the Trait Substitution Sequence (TSS) introduced by [START_REF] Metz | Adaptative dynamics, a geometrical study of the consequences of nearly faithful reproduction[END_REF].

We observe that the Substitution Fleming-Viot Process includes the three qualitative behaviors due to the three different time scales: deterministic equilibrium for the size of the population (driven by the ecological birth and death events), transitory diffusive behavior for the marker distribution (driven by marker mutation), jump process for the trait distribution (driven by adaptive trait mutation).

Remark 3.6 Equations (3.4)-(3.5) have important biological implications regarding neutral genetic diversity. Once the fixation of a favorable mutation has occurred and the population is monomorphic for the selected trait, the evolution of the neutral marker distribution is described by a Fleming-Viot process whose law is given by the martingale in (3.4). The bracket of the martingale in (3.5) shows that the stochastic fluctuations with time of the marker distribution are due to randomness in births and deaths and mutations. The multiplicative factor 2b(x)/ n x in (3.5) depends on the trait value x, and on the assumed ecological model which determines the relationships between x, the death and birth rates and the competition kernel. Notice that 2b(x)/ n x corresponds to the quotient of variance (here 2b(x)) divided by effective size N e (here n x ) that appears in the usual Wright Fisher equation. The quantity n x corresponds to the mass of the population when there is an infinite number of small individuals; if the size of the population is of order K, it means that there is approximately n x K individuals of weights 1/K. The right term in (3.4), (i.e. the drift term in a mathematical sense) involves the generator A and is associated with the mutation model as seen in Assumption (2.5). The generator A describes the speed at which the neutral diversity is restored. For instance in a continuous state model, if Aφ = σ 2 2 φ ′′ , we recover the heat equation whose solutions have a variance in t. In a discrete state model similar to (2.8), this equation gives the growth of the support. In short, (3.4)-(3.5) shows that the distribution of the neutral marker depends on ecological processes and their parameters: every changes in x will result in changes in the distribution of the neutral marker, through changes in birth, death and mutation rates, in competition and equilibrium population size. This result is biologically relevant and important since it differs from the assumptions of classical genetic hitchhiking models, in which selection and population size remain constant, leading to the fact that the neutral diversity restoration will not depend on the trait substitution and its history. In examples below, we will give more detailed results regarding the distribution of the neutral marker changes.

The proof of Theorem 3.4 is the subject of Section 4. The trait dynamics in the limit of Theorem 3.4 is the Trait Substitution Sequence obtained by the Theorem 1 in [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] whose assumptions are satisfied. Our main contribution in Theorem 3.4 is to prove that at the adaptive trait mutation time scale, a homogeneization phenomenon takes place. There is a deterministic limit for the fastest process (the births and deaths leading to n x ), and stochastic limits for the two slower processes. The limiting process (V t , t ≥ 0) is a measure-valued process with jumps (corresponding to trait mutations) and diffusion (corresponding to marker dynamics). If the population is trait-monomorphic with trait x, the jump measure is

b(x) n x X -{x} [f (x + k; x)] + b(x + k) m(x, k)dk,
where X -{x} = {y -x, y ∈ X }. When a jump occurs at t, the process jumps from (x, u) to

(x + k, v)
where k is chosen in m(x, k)dk and v is chosen at time t in the marker distribution F u t (x, dv). The marker distribution is the second fastest-evolving component, but marker mutations are assumed small (2.5), allowing to recover a non-degenerate Fleming-Viot superprocess parameterized by the trait of the population but with jumps. Between the jumps, this superprocess is the pathwise limit of the marker dynamics where traits are fixed. The jumps are hitchhiking events due to the trait mutations (see in another context [START_REF] Etheridge | An approximate sampling formula under genetic hitchhiking[END_REF]. There is a bottleneck at each successful invasion-fixation of mutant traits. Indeed, the individuals present at the fixation time are all descendants of the successful initial mutant. The trait and marker of the latter alone determine the state of the new mutant population, hence creating the bottleneck for the whole population genealogy. This result is biologically intuitive since we assume that the neutral marker and the trait are completely linked, but the mathematical proof of these phenomena is the hardest part of the proof of Theorem 3.4, and we will show that our results still have biological interest. Extending this model to the case of recombination is a challenging problem for future work (see Smadi 2014 in this direction). It is also worth to notice that contrarily to other extensions of the TSS (e.g. the TSS with age-structure of [START_REF] Méléard | Trait substitution sequence process and canonical equation for age-structured populations[END_REF] or the Polymorphic Evolution Sequence for a multi-resource chemostat in Champagnat et al 2014) that usually jump from an equilibrium to another, the marker distribution is here described by a stochastic process and not an equilibrium measure. This is due to the fact that the time scales of the trait and marker mutations are assumed different: in the time scale of marker mutations, the trait mutations are too rare and not seen.

An illustration of the invasion and fixation phenomena is summed up in Fig. 3.1.

✲ ✻ 0 ε n x+k n x n x δ x ⊗ F u t 2 (x, dv) 1 K δ (x+k,v) population size t 2 ✲ log K σ 1 t τ 1 n x+k δ (x+k,v)
Figure 3.1: Invasion and fixation of a successful trait mutant. In the population of resident trait x and marker distribution F u t (x, dv), a mutant trait x + k appears at time τ 1 . Let v be the marker of the mutant individual. As in [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF], the fluctuations of the resident population can be neglected in first approximation and the mutant population evolves as a birth and death process with rates b(x + k) and d(x + k) + η(x + k)C(k) n x , independent of the marker distribution. When the mutant population reaches a sufficient size ε at time t 2 , with probability [f (x + k, x)] + /b(x + k), the 'invasion implies fixation' assumption leads to the replacement of the former population in a time t K such that t K / log(K) → ∞. This time interval is too short to allow other marker mutant to appear in non-negligible proportion, with large probability. Thus, when the mutant population has fixed, at time σ 1 , it is close to n x+k δ (x+k,v) . Before the next adaptive trait mutation occurs, the marker mutates a lot, since marker mutations happen on a faster scale. The dynamics of the marker distribution is then the one of a Fleming-Viot superprocess started at δ v and with statistics depending on x + k.

An example from Dieckmann and Doebeli (1999)

Let us first illustrate our model by simulations based on an example inspired from Roughgarden (1979) and [START_REF] Dieckmann | On the origin of species by sympatric speciation[END_REF]. Here X = [-1, 1], U = [-2, 2] and K = 1000 in all the simulations. The individual dynamics is characterized by

• the birth rate b(x) = exp(-x 2 /2σ 2 b )
with σ b = 0.9. The probability of mutation of the trait and marker are respectively p K = 1/K 2 and q K = 1/ √ K. The adaptive trait mutation kernel m(x, k)dk is a Gaussian law with mean 0 and variance 0.1, conditioned to [-1, 1]. The marker mutation kernel G K (u, dh) is a Gaussian law with mean 0 and variance

σ 2 K = 1/ √ K, conditioned to [-2, 2].
• symmetric competition for resources, with η(x) = 1 and C(x -y) = exp(-(xy) 2 /2σ 2 C ), σ C = 0.8.

Here, the 'optimal trait' is x = 0 where the birth rate has its maximum and the population is governed by local competition. We start with the initial condition: x 0 = -1, u 0 = 0. The marker value of the initial mutant is indicated by the red line. When the mutant trait appears, the resident population is quickly invaded by the mutant population during a transition period. In (b), we can see that if the support of the marker distribution for the resident population remains wide (see also (a)), the size of the resident population decreases quickly. In the second column of (b), we see that the marker distribution in the mutant population remains spiked at the marker value of the first mutant individual during the whole transition period. After invasion (see (a)), the spread of the marker distribution follows the Fleming-Viot process (3.4)-(3.5). On (a), we see that for the Fleming-Viot process, the support of the marker distribution spreads slowly.

The simulations (see Fig. 3.2) illustrate Theorem 3.4. They show the replacement of a resident population by a mutant population. In Fig. 3.2 (a), the dynamics of the support of the marker distribution is represented. The mutant and resident populations are pictured together and separately to better observe the extinction of the resident population (black) and the expansion of the mutant population from one individual (orange). The invasion started around time 3175, and after time 3250, the mutant population has totally replaced the resident one.

In Fig. 3.2 (b), the histograms of traits and markers at three times during the invasion are represented simultaneously, to underline the hitchhiking effect of the marker during the 'invasion implies fixation' phase. We can see that the distribution of the marker during the fixation remains close to a Dirac mass at the marker value of the initial mutant (red line). This illustrates the bottleneck phenomenon, the existence of which we prove rigorously (Equation (4.8) of Proposition 4.5).

Let now focus on the Dieckmann-Doebeli's example and highlight the biological implications regarding the eco-evolutionary feedback on the distribution of the neutral marker.

Here, n x = b(x)-d(x) and therefore the Fleming-Viot process F u t (x, .) is the solution of the martingale problem (3.4) with Aφ = σ 2 2 φ ′′ and with bracket (3.5) given for all φ ∈ C(U , R) by

2 b(x) b(x)-d(x) t 0 F u s (x, .), φ 2 -F u s (x, .
), φ 2 ds. If the death rate is a constant d(x) = d, then the multiplicative factor in the bracket (3.5), b(x)/(b(x) -d), decreases when b(x) increases. Heuristically we expect that the stochastic fluctuations in time of the distribution of the neutral marker decrease when the trait x approaches the evolutionary stable strategy (ESS, see [START_REF] Smith | Evolution and the Theory of Games[END_REF]) and b(x) increases, since the equilibrium size is greater and the diffusion coefficient is lower. The drift term is b(x) σ 2 2 t 0 F u s (x, .), φ ′′ ds and thus the multiplicative factor b(x) increases when approaching the ESS, contrarily to the multiplicative factor of the bracket (3.5). In the case d ≡ 0, the Fleming-Viot process has a constant diffusion coefficient and the bracket (3.5) does not depend on x. The Fleming-Viot process depends only on the trait x through the drift term. Notice that this is true for any mutation model satisfying (2.5). This simple result illustrates how the ecological processes can shape the neutral diversity.

Corollary: the Wright-Fisher Evolutionary Process

There exists a version of the SFVP in the case when the marker space U is discrete. Assume for instance that there exist only two alleles of the marker trait, denoted by a and A, so that U = {a, A}. In this case, we apply Theorem 3.4 with the mutation kernel G K defined in (2.7) and r K /K → r > 0 when K → +∞.

Proposition 3.7 We work under Assumptions 2.1 and 3.1 with probabilities q A and q a to mutate from marker A to marker a and from marker a to marker A. Moreover, we consider similar initial conditions ν K 0 as in Theorem 3.4. Then, the population process

(ν K Kt , t ≥ 0) converges in law to the M F (X × {a, A})-valued process ( n Yt W a t δ (Yt,a) (dy, du) + (1 -W a t ) δ (Yt,A) (dy, du) , t ≥ 0), where (Y t , t ≥ 0) is the TSS process that jumps from x to x + k in X with the jump measure b(x) n x [f (x+k;x)] + b(x+k)
m(x, k)dk and where (W a t , t ≥ 0) is the following Wright-Fisher jump process that represents the proportion of alleles a in the population of trait Y t at time t. Between jumps, it satisfies the usual Wright-Fisher equation with mutations

dW a t = r b(Y t ) q A (1 -W a t ) -q a W a t dt + 2b(Y t ) n Yt W a t 1 -W a t dB t (3.8) (B t , t ≥
0) being a standard Brownian motion. It jumps with the TSS and at jump time t, the process (W a t , 1 -W a t ) goes to (1, 0) with probability W a t and to (0, 1) with probability 1 -W a t .

An illustration of this proposition is given in Fig. 3.3.

This result can be generalized to discrete marker spaces U = {a 1 , . . . a m }, by introducing the transition probabilities q ij to mutate from a i to a j , i, j ∈ {1, . . . , m}. An application is when the marker corresponds to the genetical sequence of n nucleotides (A, T , G or C for each position). In this case, m = Card U = 4 n .

Traditionally in a population genetics framework, the evolution in finite populations of the diversity at a neutral marker is described as a diffusion process with two fixed parameters: the population size and the mutation "rate" (e.g. [START_REF] Crow | An Introduction to Population Genetics Theory[END_REF]. The population size is related to what is called the "genetic drift" and generally refers to the random sampling of gametes performed for reproduction at the beginning of each generation, and the higher the population size, the lower the genetic drift. Under this framework, genetic drift induces stochastic fluctuations in the frequencies of the alleles A and a and can cause the decrease of the neutral genetic diversity when an allele is randomly lost. On the other hand, mutation introduces continuously alleles A and a in the population and thus allows the restoration and the maintenance of the neutral genetic diversity. It is important to note that under the population genetics framework, mutation rates and population size are fixed and do not depend on the ecological processes and their parameters, neither on the trait value when the population is monomorphic for the trait under selection. As a consequence, those parameters do not change as successive selective sweeps occur especially during the adaptation process. Here we can use Equation (3.8) and try to compare the classical population genetics results about the distribution of neutral diversity and the one in our model. In an eco-evolutionary framework, (3.8) first shows that mutation rates and population size, i.e. the genetic drift, are not fixed and depend on the ecological processes and on the trait value x. The mutation rates are r b(Y t )q A and r b(Y t )q a in our framework while it is only q A and q a under a population genetics framework (e.g. [START_REF] Crow | An Introduction to Population Genetics Theory[END_REF]. The genetic drift, i.e. the equilibrium population size, is given by 1/ n Yt while it is a constant 1/n in population genetics framework. Second,(3.8) shows that extra ecological processes affect the distribution of the neutral marker since in the left-hand side there is the term 2b(Y t ). This term can be interpreted as the effect of demographic stochasticity, which is not taken into account in population genetics.

Extensions to co-existing traits

The work of [START_REF] Champagnat | Polymorphic evolution sequence and evolutionary branching[END_REF] generalizes the TSS to the case of coexisting trait values, when Assumption 3.1 is relaxed. They define a polymorphic TSS called polymorphic evolutionary sequence (PES) and denoted by (X t ) t≥0 ∈ D(R + , M F (X )). When a mutant trait y appears in a resident population of trait x 0 at time t 1 , either its descendent line is killed with probability 1 -[f (y; x 0 )/b(y)] + , or it survives. In that case, we can have coexistence of y and x 0 when there is a positive globally stable non-trivial equilibrium (n * x 0 ,y , n * y,x 0 ) to the Lotka-Volterra system defined in (A.4). Therefore the population jumps from X t 1 -= n x 0 δ x 0 to

X t 1 = n *
x 0 ,y δ x 0 (dx) + n * y,x 0 δ y (dx). Fisher diffusions while the size of the population stabilizes around the equilibrium given by the trait value. A trait mutant appears around time 18290, invades and fixes into the population. Before the appearance of this mutant trait, fluctuations in the marker distribution are due to (fast) marker mutation, birth and death stochastic events. At the time when the mutant trait appears, the A-allele frequency is 85%, giving a high probability for an A-allele hitchhike. This is the case in the simulation. After fixation time (around time 18490), the a-allele population is extinct. It is regenerated by mutations of the marker but get extinct three times before taking up around time 19600.

For a probability π, a trait measure X and x ∈ X , let us denote by F t (π, x, X, du) the Fleming-Viot process started at π, evolving in the trait distribution X and parameterized by x.

Let π 0 be the initial marker distribution of the monomorphic population of trait x 0 . Before the time t 1 of appearance of the first mutant, the marker distribution evolves as (F t (π 0 , x 0 , n x 0 δ x 0 , du)) t≥0 . Let π t 1 = F t 1 (π 0 , x 0 , n x 0 δ x 0 , du) be the marker distribution at t 1 and let V 1 be a random variable drawn in the distribution π t 1 . After t 1 and before the occurence of the second trait-mutation at t 2 , the population evolves as

n * x 0 ,y δ x 0 (dx)F t-t 1 (π t 1 , x 0 , X t 1 , du) + n * y,x 0 δ y (dx)F t-t 1 (δ V 1 , y, X t 1 , du).
The processes F t (π t 1 , x 0 , X t 1 , du) and F t (δ V 1 , y, X t 1 , du) are independent generalizations of the Fleming-Viot process defined in Definition 3.3, conditionally on π t 1 , X t 1 and V 1 . Indeed their semimartingale decompositions are respectively:

F t (π t 1 , x 0 , X t 1 , .), φ = π t 1 , φ + b(x 0 ) t 0 F s (π t 1 , x 0 , X t 1 , .), Aφ ds + M 1 t (φ) ; F t (δ V 1 , y, X t 1 , .), φ = φ(V 1 ) + b(y) t 0 F s (δ V 1 , y, X t 1 , .), Aφ ds + M 2 t (φ), (3.9)
where M 1 (φ) and M 2 (φ) are independent square integrable martingales such that

M 1 (φ) t = b(x 0 ) + d(x 0 ) + η(x 0 )C(0)n * x 0 ,y + η(x 0 )C(x 0 -y)n * y,x 0 n * x 0 ,y + n * y,x 0 t 0 F s (π t 1 , x 0 , X t 1 , .), φ 2 -F s (π t 1 , x 0 , X t 1 , .), φ 2 ds = 2b(x 0 ) n * x 0 ,y + n * y,x 0 t 0 F s (π t 1 , x 0 , X t 1 , .), φ 2 -F s (π t 1 , x 0 , X t 1 , .), φ 2 ds, M 2 (φ) t = b(y) + d(y) + η(y)C(y -x 0 )n * x 0 ,y + η(y)C(0)n * y,x 0 n * x 0 ,y + n * y,x 0 t 0 F s (δ V 1 , y, X t 1 , .), φ 2 -F s (δ V 1 , y, X t 1 , .), φ 2 ds = 2b(y) n * x 0 ,y + n * y,x 0 t 0 F s (δ V 1 , y, X t 1 , .), φ 2 -F s (δ V 1 , y, X t 1 , .), φ 2 ds. (3.10)
At time t 2 , when a third trait appears in the population, the system can evolve to three two or just one coexisting traits, depending on the new trait equilibrium of the Lotka equations that is reached. For each of the traits, the marker distribution evolves as a generalization of the Fleming-Viot processes above.

Remark 3.8 The above equations show that, when there is coexistence of two traits in the population, the markers in the subpopulations defined by the two traits evolve independently but with parameters depending on the two co-existing traits. Thus, when there is a diversification event in the population, the distribution of the neutral diversity in one of the two subpopulations does not evolve as completely forgetting the other one. The parameters of the underlying Fleming-Viot process depend on the complete trait distribution.

We present in Figure 3.4 simulations in the case of coexistence, with the same model and parameters as in Section 3.3, except σ C = 0.7 and the initial condition: x 0 = -0.1. The simulations (see Fig. 3.4) show the appearance of a new mutant trait (yellow) in a population of two coexisting traits (black and blue).

4 Proof of Theorem 3.4

Let us sketch the proof. In this section, we will suppose that Assumptions 2.1, 3.1 are satisfied and the initial conditions are ν K 0 (dy, dv) = n K 0 δ (x 0 ,u 0 ) (dy, dv) with lim K→∞ n K 0 = n x 0 and sup K∈N * E((n K 0 ) 3 ) < +∞. First, we recall results due to [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF] that provide the finite marginal convergence of the trait process (X K Kt ; t ≥ 0). We extend these results to obtain the weak convergence of the measures (X K Kt (dx)dt; K ≥ 0) in M F (X × [0, T ]) embedded with the weak convergence topology. This corresponds to the convergence of (X K Kt ; t ≥ 0) as a trait-marker-time measure, as developed by [START_REF] Kurtz | Averaging for martingale problems and stochastic approximation In Applied stochastic analysis[END_REF]. Secondly, we include the fast component (the marker) and prove the tightness of the sequence (ν K Kt (dx, du)dt; K ≥ 0) in M F (X ×U ×[0, T ]). We then consider a subsequence, again denoted by (ν K Kt (dx, du)dt, K ≥ 0) with an abuse of notation, that converges to a limit Γ(dt, dx, du) ∈ M F ([0, T ] × X × U) that we have to identify. This derivation is done in several steps. When a successful mutant appears in the monomorphic population with trait x, the transition period to fixation is to be considered carefully. It has been proved by [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] that these transitions are of order log(K). We prove that during this time interval, the marker distribution in the mutant subpopulation remains a Dirac mass at the value of the initial mutant. This results from the combined effects of small or rare marker mutations, large population and slow take-off of the new mutant population. Then, we show that in a trait monomorphic population with value x, the marker distribution converges to a Fleming-Viot superprocess parameterized by x.

Semimartingale decomposition of ν K

Let us introduce some notation to keep forthcoming formula simple. For ν ∈ M F (X × U ) and φ(x, u) ∈ C(X × U, R), we define the (nonlinear) generators B K and D K (ν) such that

B K φ(x, u) =(1 -p K )(1 -q K )b(x)φ(x, u) +p K (1 -q K )b(x) X φ(x + k, u)m(x, k)dk +q K (1 -p K )b(x) U φ(x, u + h)G K (u, dh) +p K q K b(x) X ×U φ(x + k, u + h)m(x, k)dk G K (u, dh) (4.1) D K (ν)φ(x, u) = d(x) + η(x)C * ν(x) φ(x, u). (4.
2)

The process ν K . , φ is a square integrable semi-martingale and we give its characteristics.

Proposition 4.1 For a continuous bounded function φ(x, u) on X × U, the process

M K,φ t = ν K t , φ -ν K 0 , φ - t 0 ds X ×U ν K s (dx, du) B K -D K (X K s ) φ(x, u) (4.3)
is a square integrable martingale with previsible quadratic variation

M K,φ t = 1 K t 0 ds X ×U ν K s (dx, du) B K + D K (X K s ) φ 2 (x, u). (4.4)
Proof The dynamics being given in Section 2, the proof can be adapted from (Fournier and Méléard, 2004, Lemma 5.2). One main step consists in showing that there exists a Poisson point measure driving the measure-valued processes ν K for all K ∈ N * .

Convergence of the trait-marginal in the trait mutation time scale

As previously emphasized, the trait dynamics is described by the measure-valued process X K which does not depend on the markers. This process has been fully studied in [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] and [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF]. In this section, we recall the finite marginal convergence result obtained in these papers. We give some additional properties concerning the topology involved. This result shows a time scale separation with successive fixations of successful mutants, under Assumptions 2.1 and 3.1. Notice that the time scale assumption is

∀V > 0, log K ≪ 1 Kp K ≪ exp(V K), as K → ∞, (4.5)
which is realized in our case for p K = 1/K 2 .

Theorem 4.2 Under Assumptions 2.1 and 3.1, let us also assume that the initial population is trait-monomorphic: X K 0 = n K 0 δ x for x ∈ X and n K 0 → n x in probability and sup K∈N * E((n K 0 ) 3 ) < +∞. Then, the sequence (X K Kt ; t ≥ 0) converges to the pure jumps singleton measure-valued Markov process ( n Yt δ Yt ; t ≥ 0) defined as follows: Y 0 = x, and the process Y jumps from

x to x + k with jump measure b(x) n x [f (x+k;x)] + b(x+k) m(x, k)dk.
The convergence holds in the sense of finite dimensional distributions on M F (X ) equipped with the topology of total variation.

This theorem has been proved by [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] for the logistic case and generalized by [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF].

The trait-marginal process (X K Kt ; t ∈ [0, T ]) does not converge in D([0, T ], M F (X )) embedded with the Skorokhod topology. Indeed, the size of jumps is upperbounded by 1 K and nevertheless the limiting total mass process has jumps, preventing trajectorial tightness (at least in the J1-topology). Following the idea of [START_REF] Kurtz | Averaging for martingale problems and stochastic approximation In Applied stochastic analysis[END_REF] and as developed in [START_REF] Méléard | Slow and fast scales for superprocess limits of age-structured populations[END_REF] and [START_REF] Gupta | A new proof for the convergence of an individual based model to the trait substitution sequence[END_REF], a weaker topology consists in forgetting the process point of view and considering the measure X K Kt (dx)dt in M F ([0, T ] × X ) embedded with the topology of weak convergence. This convergence strengthens the result of Theorem 4.2 but in a topology weaker than the Skorohod topology. To achieve this, as in [START_REF] Collet | A rigorous model study of the adaptative dynamics of Mendelian diploids[END_REF], we first introduce the M 1 -topology on D([0, T ], R + ). It is weaker than the usual J 1 -topology and allows monotonous processes with jumps tending to 0 to converge to processes with jumps (see [START_REF] Skorohod | Limit theorems for stochastic processes[END_REF]. For a càdlàg function h on [0, T ], the continuity modulus for the M 1 -topology is given by

w δ (h) = sup 0≤t 1 ≤t≤t 2 ≤T ; 0≤t 2 -t 1 ≤δ d(h(t), [h(t 1 ), h(t 2 )]). (4.6)
Note that if the function h is monotone, then w δ (h) = 0.

Proposition 4.3 Let us consider a continuous function g. Then, under Assumptions 2.1 and 3.1, the process (R K t , t ∈ [0, T ]) defined by

R K t = g(x)X K Kt (dx)
converges in law in the sense of the Skorohod M 1 -topology to the process (R t , t ∈ [0, T ]) where R t = n Yt g(Y t ).

Proof From Theorem 4.2, finite dimensional distributions of (R K t , t ∈ [0, T ]) converge to those of ( n Yt g(Y t ), t ∈ [0, T ]). By Skorohod (1956, Theorem 3.2.1), it remains to prove that for all η > 0, lim

δ→0 lim sup K→∞ P(w δ (R K . ) > η) = 0,
where w δ has been defined in (4.6).

The mutation rate in (R K t , t ∈ [0, T ]) being bounded, the probability that two mutations occur within a time less than δ is o(δ). It is therefore enough to study the case where there is at most one mutation in the time interval [0, δ]. Following [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF], the path of R K can be decomposed into several subpaths, each of them being closed to a large population deterministic measure-valued function ξ (See Proposition A.1 in the appendix) with a probability tending to 1. Away from invading mutations and for a traitmonomorphic population with trait x, ξ Kt , g = g(x)n Kt (x) where n . (x) is the solution of the logistic equation (3.1). We can easily check that t → n t (x) converges monotonously to its stable equilibrium n x and then ξ Kt , g is monotonous and the modulus of continuity tends to 0. Around an invading mutant y , ξ Kt , g is close to n Kt (x)g(x) + n Kt (y)g(y) where (n t (x), n t (y)) is solution of the Lotka-Volterra system (A.4) with an initial condition close to ( n x , 0). The mutant y invades if the fitness function f (y; x) is positive (and f (x; y) is negative). From Assumption 3.1, an easy study of the Lotka-Volterra system (see for example the appendix in Champagnat 2004, Figure (b) p.187), shows that either n t (x) and n t (y) are increasing or ṅt (x) < 0; ṅt (y) > 0. In that case, n t (x)g(x) + n t (y)g(y) is the sum of two monotonous functions and the modulus of continuity tends also to 0.

Corollary 4.4 The sequence of random measures X K Kt (dx)dt converges in law to the random measure n Yt δ Yt (dx)dt in M F ([0, T ]×X ) embedded with the weak convergence topology.

Proof It is enough to prove the convergence in law of h(t)e -qx X K Kt (dx)dt to h(t)e -qx nYt δ Yt (dx)dt for a measurable bounded function h and q ∈ Q. In [START_REF] Skorohod | Limit theorems for stochastic processes[END_REF], it is proved that if x K converges to x in D([0, T ], R) embedded with the M 1 -topology, then for t outside a denumerable set, x K (t) converges to x(t). Then it follows by Lebesgue's Theorem that T 0 H(t, x K (t))dt converges to T 0 H(t, x(t))dt, as soon as H is bounded and continuous. We apply this result to the process ( X e -qx X K Kt (dx), t ≥ 0) and the function

H M (t, y) = h(t)(y ∧ M ),
for any M > 0. Estimate (2.9) (with p = 1) allows to conclude.

Marker distribution in a new adaptive trait mutant population

In this section, we study the transition of the marker distribution when a new mutant adaptive trait appears in a monomorphic population with trait x 0 . We consider this phenomenon at the ecological time scale and we prove that the fixation of the mutant trait creates a genetical bottleneck.

Let K be fixed. Initially we have a trait monomorphic population with trait x 0 and a marker distribution π K (x 0 , du). Then an individual (x 0 , v) from this population gives birth to an individual with mutant trait y and marker v (v has been chosen according to π K (x 0 , du)). We consider the process (ν K t ; t ≥ 0) started at

ν K 0 (dx, du) =X K 0 (dx)π K 0 (x, du) = 1 K δ (y,v) (dx, du) + N K 0 -1 K δ x 0 (dx)π K 0 (x 0 , du).
Proposition 4.5 Under Assumptions 2.1 and 3.1, let us consider a mutant (y, v) appearing in a monomorphic population with trait x 0 and marker distribution π K 0 (x 0 , du). Let us assume that f (y; x 0 ) > 0, where the fitness function has been defined in (3.3). There exists ε > 0 such that for any sequence (t K ; K ∈ N * ) with lim K→+∞ t K / log K = +∞ and lim K→+∞ t K /K = 0 (for example t K = (log K) 2 ), we have

lim K→+∞ P ν K t K , 1l y > ε = f (y; x 0 ) b(y) and lim K→+∞ P ν K t K , 1l y = 0 = 1 - f (y; x 0 ) b(y) . (4.7)
Further, for the marker distribution, we can prove that

lim K→+∞ P π K t K (y, du) = δ v (du) = f (y; x 0 ) b(y) . (4.8)
The equation (4.8) tells us that when the mutant trait survives in the resident population of trait x 0 , then by the time t K it needs to reach a non-negligible size, its marker distribution is still a Dirac mass at y. The assumption q K (log K) 2 in Assumption 2.1 ensures this. This assumption is not very restrictive as (log K) 2 is a very slow growth rate. Additional comments are given after the proof.

Proof Properties (4.7) have been proved in [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] and [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF] and depend only on the trait distribution. We consider test functions φ(x, u) of the form 1l y (x)g(u) with g ∈ C 2 (U , R) such that g ∞ + g ′′ ∞ ≤ 1. Starting from Proposition 4.1 and using Itô's formula with jumps, we obtain as soon as the population with trait y survives,

U g(u)π K t K (y, du) = ν K t K , 1l y g ν K t K , 1l y =g(v) + M K,g t K + q K (1 -p K ) b(y) t K 0 1 - 1 K ν K s , 1l y + 1 (4.9) × U π K s (y, du) U g(u + h) -g(u) G K (u, dh) ds
where M K,g is a square integrable martingale with previsible quadratic variation:

M K,g t K = 1 K t K 0 ds b(y)(1 -q K )(1 -p K ) ν K s , 1l y ν K s , 1l y + 1 K 2 U g(u) -π K s , g 2 π K s (y, du) + d(y) + η(y)C * ν K s (y) ν K s , 1l y ν K s , 1l y -1 K 2 U g(u) -π K s , g 2 π K s (y, du) +b(y)q K (1 -p K ) ν K s , 1l y ν K s , 1l y + 1 K 2 U π K s (y, du) U G K (u, dh) g(u + h) -π K s , g 2 .
(4.10)

The third term in the right hand side of (4.9) is of order t K /K. Indeed thanks to (2.4) and (2.5), it is upper bounded by t K K b Ag ∞ .

Equation (4.10) needs more attention. As soon as the mass ν K s , 1l y of the mutant population is of order 1, the variance of M K,g t K is in t K /K which tends to zero when K → +∞. However, since we start from 1 individual, we have to separate the time interval [0, t K ] into 2 parts. Let us introduce a sequence (s K ) such that s K ≤ t K for any K and

log K ≪ s K ≪ 1 √ q K .
This is possible thanks to the assumption that q K (log K) 2 → 0. Notice that s K can be equal to t K . Using Assumption 3.1, we can prove as in Champagnat (2006, Lemma 3) that there exists ε 0 > 0 such that

lim K→∞ P ∀s ∈ [s K , t K ], ν K s , 1l y ≥ ε 0 = f (y; x 0 ) b(y) .
It turns immediately out that

E 1l {∀s∈[s K ,t K ], ν K s ,1ly >0} K t K s K b(y) + d(y) + η(y)C * ν K s (y) ν K s , 1l y U g(u) -π K s , g 2 π K s (y, du)ds ≤ C t K K . (4.11)
Before time s K , the population size with trait y is not large enough and 1 K ν K s ,1ly can only be upper bounded by 1. Therefore we have to control the expectation of the variance of g under π K s . The expected number of marker mutations at time s along a lineage is sq K and the variance of such mutation is bounded by g 2 ∞ = sup{g(h) 2 , h ∈ U}. Then E U g(u) -π K s , g 2 π K s (y, du) ≤ s q K g 2 ∞ , (4.12)

and

E 1 K s K 0 b(y) + d(y) + η(y)C * ν K s (y) ν K s , 1l y U g(u) -π K s , g 2 π K s (y, du)ds ≤ C (s K ) 2 r K K 2 .
The upper bound converges to 0 by the assumption on q K . The third term of (4.10) is treated similarly. This concludes the proof.

Remark 4.6 For q K = 1/ √ K, let us notice that the rate of appearance of mutant markers in a population of size K is of order Kq K = √ K which does not tend to zero. This means that many mutant markers appear in the population of trait y during the t K time interval following the first mutant (y, v). However, heuristically, since in a tree the mass is concentrated around the leaves, the mutants do not appear with the same probability along the time interval and mutations are mostly observed after the time s K when the mutant population (y, v) is already large. Moreover, using that the marker mutation step and/or marker mutation frequency is small we obtain that the mutant markers remain in negligible proportion between s K and t K . tends to 0 in L 1 when K → +∞. Thus, lim K→+∞ E Ψ s,t ( π K,x 0 K. ) = 0. (4.19)

On the other hand, using (2.9) and the convergence of ( π K,x 0 K(.∧τ K 1 ) (du); K ∈ N * ) to π ∈ C([0, T ], M F (U )), we get E Ψ s,t (π) = lim K→+∞ E Ψ s,t ( π K,x 0 K. (du)) .

(4.20)

This shows that E Ψ s,t (π) = 0 and hence the process M x 0 (φ) defined in (3.4) is a martingale obtained as the uniform limit in time of H K,φ Kt , when K → +∞. Moreover, the bracket (4.17 Indeed, the integral in (4.17) can be separated into two integrals, one between 0 and log K K ∧ t ∧ τ K 1 and the other between log K K ∧ t ∧ τ K 1 and t ∧ τ K 1 . The second integral converges to (4.21), but some caution is needed for the first integral since the ratios ν K Ks , 1l x 0 / ν K Ks , 1 ± 1 K 2 are of order K. Using the same arguments as for (4.12), we can upper bound the integral between 0 and log K K ∧ t ∧ τ K 1 by

CK log K/K 0 sq K ds = C (log K) 2 K q K → K→+∞ 0.
Using Theorem 3.12 p. 432 of [START_REF] Jacod | Limit Theorems for Stochastic Processes[END_REF], that provides the convergence of H K,φ K. to the solution of the martingale problem (3.4)-(3.5) with x = x 0 . By the independence of π K,x 0 K. (du) and τ K 1 , τ 1 is independent of π. (du) and π.∧τ 1 = π u 0

.∧τ 1 (x 0 , du). This concludes the proof.

Conclusion

Using Theorem 4.2, Proposition 4.5 and Proposition 4.7, we prove the first part of Theorem 3.4, for the convergence in finite distribution. Let us now consider the convergence in the space of trait-marker-time measures.

Corollary 4.8 The family (ν K Kt (dx, du) dt, K ∈ N * ) is uniformly tight in M F (X × U × [0, T ]) embedded with the weak convergence topology and converges in distribution to the measure V t (dx, du)dt, where V is defined in Theorem 3.4.

Proof Since the space X × U × [0, T ] is compact, it is sufficient to prove that sup

K∈N * E T 0 ν K Kt , 1 dt < +∞

Figure 3 . 2 :

 32 Figure 3.2: We consider a resident trait-monomorphic population (black) in which a mutant trait (orange) appears and goes to fixation. Here the intrinsic death rate d(x) = 0. (a) Evolution of the support of the marker distribution with respect to time; the support of the resident trait-monomorphic population is in black while the support of the mutant population is in orange. The mutant and resident populations are shown separately and together. (b) Distributions of the traits (left) and markers (right) in the population at three times during the invasion. The marker and trait values are in abscissa and frequencies are in ordinate.The marker value of the initial mutant is indicated by the red line. When the mutant trait appears, the resident population is quickly invaded by the mutant population during a transition period. In (b), we can see that if the support of the marker distribution for the resident population remains wide (see also (a)), the size of the resident population decreases quickly. In the second column of (b), we see that the marker distribution in the mutant population remains spiked at the marker value of the first mutant individual during the whole transition period. After invasion (see (a)), the spread of the marker distribution follows the Fleming-Viot process (3.4)-(3.5). On (a), we see that for the Fleming-Viot process, the support of the marker distribution spreads slowly.

Figure 3 . 3 :

 33 Figure3.3: Evolution of sizes of the subpopulations with markers a and A. The simulation uses individual-based algorithms. The proportions of marker alleles a and A follow Wright-Fisher diffusions while the size of the population stabilizes around the equilibrium given by the trait value. A trait mutant appears around time 18290, invades and fixes into the population. Before the appearance of this mutant trait, fluctuations in the marker distribution are due to (fast) marker mutation, birth and death stochastic events. At the time when the mutant trait appears, the A-allele frequency is 85%, giving a high probability for an A-allele hitchhike. This is the case in the simulation. After fixation time (around time 18490), the a-allele population is extinct. It is regenerated by mutations of the marker but get extinct three times before taking up around time 19600.

  Figure 3.4: Neutral marker dynamics in a trait-dimorphic population. Evolution of markers with respect to time. (a) The mutant population (yellow) and resident populations (black and blue) are shown separately and together. (b) Distributions of the traits and markers in the population at three times during the invasion. The marker value of the initial mutant is indicated by the red line.

  0 ) + d(x 0 ) + η(x 0 ) n x 0 n x 0 U π(du) φ(u) -πs , φ 2 ds = t 0 2b(x 0 ) n x 0 πs , φ 2 -π, φ 2 ds. (4.21)

Proof The convergence in large population of (ν K ; K ∈ N * ) is a consequence of Proposition A.1.

S.B., S.M. and V.C.T. have been supported by the ANR MANEGE (ANR-09-BLAN-0215), the Chair "Modélisation Mathématique et Biodiversité" of Veolia Environnement-Ecole Polytechnique-Museum National d'Histoire Naturelle-Fondation X. V.C.T. also acknowledges support from Labex CEMPI (ANR-11-LABX-0007-01). R.F. acknowledges support from National Science Foundation Award EF-0623632, the Institut Universitaire de France, and the Agence Nationale de la Recherche ("EVORANGE" grant).

Convergence of the marker distribution process in a trait-monomorphic population

For K ∈ N * , we introduce, as in [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF], the following sequence of stopping times τ K k and θ K k :

The times τ K k 's are the times of appearance of the successive mutant traits in the population and the θ K k 's are the times at which the population returns to monomorphic state. These times are possibly infinite, if the corresponding sets are empty. It has been proved in [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] that for t K be such that lim K→+∞ t K / log(K) = +∞ and lim K→+∞ t K /K = 0,

In the trait-mutation time scale, the time of first mutation converges in distribution as follows:

where τ 1 is an exponential time with parameter b(x 0 ) n x 0 .

(ii) Let us consider the processes (π K K(t∧τ K 1 ) ; t ∈ [0, T ]) stopped at the time of first mutation. When K → +∞, this sequence converges in distribution in D([0, T ], P(U )) to the Fleming-Viot process F u 0 (x 0 , du) (see Definition 3.3) and stopped at the independent exponential time τ 1 .

Proof First of all, the trait and marker mutations are independent. Thus, the stopping time τ K 1 is independent of the marker distribution π K (x 0 , du). The results of Champagnat and coauthors [START_REF] Champagnat | A microscopic interpretation for adaptative dynamics trait substitution sequence model[END_REF] and [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF] are unchanged and give (4.14). Moreover, by Champagnat et al (2008, Lemma 5.4)

Since the population is trait-monomorphic with trait x 0 , then

Kt ,1 . Thus, from Proposition 4.1 and Itô's formula, we get that in the time scale Kt

where H K,φ is a square integrable martingale with quadratic variation is

The computation shows that the order of the quadratic variation of π K t is 1 K . Thus at time scale Kt, this order will be 1. That justifies Assumption (2.4) for p K which is the only choice to get a non degenerate diffusive limit.

Let us introduce a process ( π K,x 0 t (du), t ≥ 0) coupled with (π K t (x 0 , du), t ≥ 0), on the same probability space and driven by the same Poisson point measures, that satisfies the following properties. The dynamics of π K,x 0 .

(du) is given by (4.16)-(4.17) but without the stopping times τ K 1 and we have that ∀t ≥ 0,

In a nutshell, ( π K,x 0 t , t ≥ 0) corresponds to the process (π K t (x 0 , .), t ≥ 0) that is obtained by setting the trait mutation kernel to the Dirac mass at 0. Thanks to (2.5), (4.15) and using that π K,x 0 is a probability-valued process, (4.16) and (4.17) imply that for any φ ∈ C(U , R), the distribution sequence of ( π K,x 0 K. , φ ; K ∈ N * ) is uniformly tight in D([0, T ], R). By Roelly's criterion (Roelly, 1986, Theorem 2.1), this implies the tightness of the sequence of the laws of ( π K,x 0 K. ; K ∈ N * ) in D([0, T ], P(U )). Let us consider a limiting value (π t (du); t ∈ [0, T ]) of the tight sequence and a subsequence, again denoted by π K,x 0 K. (du), that converges to π. (du). By Assumption (4.15) and since individuals have weight 1/K, the limiting laws only charge C([0, T ], P(U )).

It remains to identify π. (du). Let

On the one hand, using (4.16), we obtain that

which is a consequence of Fubini's theorem since

Estimate (2.9) concludes the proof of tightness.

Let us now consider continuous functions φ ∈ C(X × [0, T ], R) and g ∈ C(U , R), where the stopping times τ K k and θ K k have been introduced in Section 4.3. Then

The limit (4.13) implies that the second term of the right hand side of (4.22) converges to 0. Given X K , the processes

, for k ≥ 0, in the first term of the r.h.s. of (4.22) are independent and, by Proposition 4.7, they converge in distribution in D([0, T ], R) to the Fleming-Viot processes (3.4)-(3.5) with the initial conditions described by the jumps of the extended TSS (Y, U ). Corollary 4.4 and dominated convergence theorem allows us to conclude the proof.

A Limit theorems for the trait and marker distributions in the ecological time scale

If we let K → +∞ without changing the time scale, we obtain:

Proposition A.1 Assume that the sequence (ν K 0 (dx, du); K ∈ N * ) converges in probability to the measure ξ 0 (dx, du) when K → +∞. Then the sequence of processes (ν K t (dx, du); t ≥ 0) K∈N * converges in D(R + , M F (X × U)) to the deterministic process ξ ∈ C(R + , M F (X × U )) defined for every φ(x, u) ∈ C(X × U, R) by:

Proof The proofs proceed in a classical tightness-uniqueness way (cf. [START_REF] Fournier | A microscopic probabilistic description of a locally regulated population and macroscopic approximations[END_REF].

Notice that no mutation can be seen at this scale. To see the trait mutations, we have to consider the process at the mutation scale Kt (cf. [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF].

Trait-monomorphic case From Proposition A.1:

Corollary A.2 Assume that the initial population is trait-monomorphic ν K 0 (dx, du) = n K 0 δ x 0 (dx)π K 0 (x 0 , du) where lim K→+∞ n K 0 = n 0 and lim K→+∞ π K (x 0 , du) = π 0 (x 0 , du) in probability. Then the sequence (ν K ; K ∈ N * ) converges, in probability and uniformly on every compact time interval [0, T ] with T > 0, to (ν t (dx, du) = n t (x)δ x 0 (dx)π 0 (x 0 , du); t ≥ 0) where n t (x) is the deterministic solution of the logistic equation

which converges when t tends to infinity to

Proof For the part of the proof dealing with ξ, we refer to [START_REF] Champagnat | From individual stochastic processes to macroscopic models in apative evolution[END_REF].

Trait-dimorphic case From Proposition A.1:

where n K 0 (x) is the number of individuals with trait x renormalized by K. We assume that for x ∈ {x 1 , x 2 }, lim K→+∞ n K 0 (x) = n 0 (x) > 0 and lim K→+∞ π K (x, du) = π 0 (x, du) in probability. We also assume that (x 1 , x 2 ) satisfies the Assumption 3.1. Then the sequence (ν K ; K ∈ N * ) converges, in probability and uniformly on every compact time interval [0, T ] with T > 0, to (ν t (dx, du) = n t (x 1 )δ x 1 (dx)π 0 (x 1 , du) + n t (x 2 )δ x 2 (dx)π 0 (x 2 , du); t ≥ 0) where (n t (x 1 ), n t (x 2 )) solves the system dn t (x 1 ) dt = b(x 1 ) -d(x 1 ) -η(x 1 )(C(0)n t (x 1 ) + C(x 1 -x 2 )n t (x 2 )) n t (x 1 ) dn t (x 2 ) dt = b(x 2 ) -d(x 2 ) -η(x 2 )(C(x 2 -x 1 )n t (x 1 ) + C(0)n t (x 2 )) n t (x 2 ). (A.4)

whose only stable equilibrium is (0, n x 2 ), with n x 2 defined in (3.2).