Gene surfing in expanding populations

Large scale genomic surveys are partly motivated by the idea that the neutral genetic variation of a population may be used to reconstruct its migration history. However, our ability to trace back the colonization pathways of a species from their genetic footprints is limited by our understanding of the genetic consequences of a range expansion. Here, we study, by means of simulations and analytical methods, the neutral dynamics of gene frequencies in an asexual population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a "surfing" mechanism [Edmonds, C.A., Lillie, A.S., Cavalli-Sforza, L.L., 2004. Mutations arising in the wave front of an expanding population. Proc. Natl. Acad. Sci. 101, 975-979]. We quantify this phenomenon in terms of (i) the spatial distribution of lineages that reach fixation and, closely related, (ii) the continual loss of genetic diversity (heterozygosity) at the wave front, characterizing the approach to fixation. Our stochastic simulations show that an effective population size can be assigned to the wave that controls the (observable) gradient in heterozygosity left behind the colonization process. This effective population size is markedly higher in the presence of cooperation between individuals ("pushed waves") than when individuals proliferate independently ("pulled waves"), and increases only sub-linearly with deme size. To explain these and other findings, we develop a versatile analytical approach, based on the physics of reaction-diffusion systems, that yields simple predictions for any deterministic population dynamics. Our analytical theory compares well with the simulation results for pushed waves, but is less accurate in the case of pulled waves when stochastic fluctuations in the tip of the wave are important.