the genetic composition of a natural population is always changing because mutations inevitably arise and rise or fall in frequency due to selection or drift ( i.e. stochasticity in finite populations ) .
genetic changes continue during invasions , when a population is entering and becoming established in new territory .
ecologists recognize that selection can play a crucial role in determining the success , speed and ultimate extent of invasions ( hoffmann and sgr 1995 ; sakai et al .
even neutral genetic change during invasions is of high interest to researchers who use neutral markers to infer the migration history of populations ( semino et al .
2000 ; sunnucks 2000 ; schaal and leverich 2001 ; wang and whitlock 2003 ; kawiecki and ebert 2004 ) .
this is because to test any hypothesis concerning the past spatial distribution of a population , it is necessary to specify a null model that incorporates explicit assumptions about demography and dispersal patterns .
it is further necessary to understand the distribution of neutral allele frequencies under the specified model , in order to compare observed marker data with the distribution expected under the null hypothesis . yet , despite the recognized importance of genetic changes during invasions , theoreticians have paid little attention to this topic until recently ( work on the evolution of dispersal constitutes an important exception ; see , for example , holt 2003 ; simmons and thomas 2004 ) .
recently , several authors have used stochastic simulations to study the fate of a single mutant that arises during an invasion ( edmonds et al .
mutation surfing , in which mutants arise on or advance to the front of a traveling wave of population density and essentially shut out wild types ( i.e. nonmutants ) from a certain point in the invasion onwards ( edmonds et al .
early studies of this topic offered a number of descriptive statements about the fate of mutations in invasions , including a 
midpoint rule for the location of the center of population density of a successful neutral mutation ( edmonds et al .
2004 ) and a proposal that a certain  lumped model parameter has a strong linear correlation with the ultimate fraction of the invading population that carries such a mutation ( klopfstein et al . 2005 ) .
more recent studies have treated surfing with selection , using mathematical models such as annihilating random walks , generalized diffusion equations and stochastic partial differential equations to explain the phenomena observed in surfing simulations and in vitro experiments ( hallatschek et al .
2007 ; hallatschek and nelson 2009 ) ( korolev et al . , 2010 ) .
recently , it has been argued that some genetic differences between human populations that had previously been attributed to selection in fact resulted from surfing by neutral alleles ( hofer et al .
one purpose of the present work is to provide a careful qualitative and quantitative description of neutral mutation surfing as seen in a stochastic model like those studied in edmonds et al .
we focus on a model of a one - dimensional habitat but include some results for two - dimensional habitats .
it should be noted that a one - dimensional habitat is a realistic model for certain types of invasions , such as invasion along a coastline or river ( lubina and levin 1988 ; speirs and gurney 2001 ; pachepsky et al .
2005 ) . as a consequence , study of such models and comparisons between them and
first , neutrality is simpler than selection , and with so little existing theory , it is reasonable to study the simpler case first .
second , neutral genetic markers are of interest because they can provide information about the history of an invasion .
indeed , much of the original interest in mutation surfing was among researchers whose primary goal is to reconstruct range expansion ( e.g. the spread of humans into europe ) with such markers ( edmonds et al .
2004 ) . on the other hand , adaptive change during invasions may be the genetic phenomenon of most practical interest to conservation biologists .
accordingly , we have begun to extend our models to the cases of beneficial and deleterious mutations , and we present some results here .
our main goals are to describe how the probability of surfing depends on model parameters ( with or without selection ) , to explain heuristically the nature of this dependence , and to offer a simple model of the surfing process as a contribution to the development of analytic models that yield quantitative predictions about genetic change during invasions .
our work is based on data obtained from a series of simulations of cellular automata . in
what follows , we state the specifications of the simulations , use statistical methods ( in particular , logistic regression ) to describe the probability of surfing and to assess our analytic model , and offer likely explanations for the quantitative results we obtain .
the model we studied , following edmonds et al . ( 2004 ) , klopfstein et al .
( 2007 ) , is a type of individual - based model known to mathematicians as a contact process ( liggett 1999 ) .
we simulated a contact process in which wild - type ( i.e. nonmutant ) and mutant individuals reproduce asexually and move between adjacent cells in a rectangular grid . for the neutral case ,
grid lengths l used were 100 , 200 and 400 cells ; grid widths w used were 1 , 3 , 7 , 13 and 25 cells . for the case of selection , only 1400 grids were used .
we note that previous studies used only 25  100 grids ( edmonds et al .
we varied grid width in order to study the effect of dimensionality on the probability of surfing .
we varied grid length in order to make sure that numerically ascertained probabilities of surfing on a finite grid came close to asymptotes , which we expect to correspond to probabilities of surfing on an infinitely long grid .
accordingly , all results below pertain to grids of length 400 unless otherwise specified . as in klopfstein
( 2005 ) , each simulation run began with a single wild - type individual placed at the center of the leftmost column of the grid .
for example , the founder could be placed along a side or in the middle of the grid to model colonization beginning other that at the mouth of a river .
first , each individual was replaced in the same cell by a number of offspring chosen from a poisson distribution with mean . for the neutral case , we used   1 = 0.05 , 0.1 , 0.2 , 0.4 , 0.8 , 7 , 15 and 31 ; in each simulation run ,  was the same for mutants and wild types . for the case with selection , we used wt  1 = 0.1 and 0.4 for wild types and mut  1 = (wt  1 ) for mutants , with  = 0.5 , 0.95 , 1.05 and 1.5 .
next , each cell whose population n ( of mutant and wild types combined ) was greater than a carrying capacity k underwent culling , in which nk randomly chosen individuals were discarded .
all individuals within a cell , whether mutant or wild type , had equal probability of being culled . for the neutral case with 
< 2 , we used carrying capacities k = 10 , 50 , 100 and 500 ; we used k = (   1)/2 for all  > 2 . for the case with selection
each individual ( whether mutant or wild type ) had a probability m of migrating in each generation .
for the neutral case with  < 2 , we used mean migration rates m = 0.05 , 0.1 , 0.2 and 0.4 ; with  > 2 , we used m = 0.1 and 0.2 . for the case with selection , we used m = 0.05 and 0.2 .
each individual chosen to migrate then moved to a cell chosen with equal probability from all those adjacent to its original cell .
( cells in the interior of the grid had four neighbors , those on the boundary had three and those in the corners had two ; diagonal moves were not permitted ) .
predictors , estimated coefficients with standard errors and 95% confidence interval lower and upper bounds for odds ratios for one - dimensional habitat with k = 10 when at least one member of the population reached the 10th column from the left , one ( wild - type ) individual chosen randomly with equal probability from all individuals in the 10th column was replaced by a mutant .
thus , each offspring shared the type of its parent , with the sole exception being the original mutant .
a run was stopped after the first generation in which either the entire population went extinct , mutants went extinct after being introduced or at least one individual reached the rightmost column of the grid . after a run was stopped ,
a code for the state of the mutant and wild - type populations was saved .
three states were distinguished for each type : extinction ( no individuals of a given type persisted ) , persistence without advancing ( at least one individual remained , but no individuals of the specified type were found in the rightmost column ) and advancing ( at least one individual of the specified type was found in the rightmost column ) .
experimentation showed that wild types virtually never went extinct after a mutation event had occurred ; so , runs with mutation events were counted as successful , even though they were usually aborted if the mutants subsequently went extinct .
( as the total population size was close to 10k times the habitat width when the mutant was introduced , it is not surprising that extinction almost never occurred after this time . )
surfing was considered to occur in runs where a mutant advanced to the rightmost column .
when this occurred , it was extremely rare for wild types to also be found in the rightmost column .
in addition to outcome codes for each run , we saved the state of the 7h subgrid centered on the 10th column for the three generations immediately following introduction of the mutant .
for the neutral case with a one - dimensional ( 1  400 ) grid , each possible combination of the model parameters  , k and m was used in 600 simulation runs . for the case of selection , each combination of the model parameters was used in 10 000 runs .
for each two - dimensional ( w  400 with w = 3 , 7 , 13 and 25 ) grid , each parameter was used in 50 simulation runs .
we also performed small numbers of runs with parameter combinations other than those listed here ; specifications for these runs are given when they are discussed below . unlike in klopfstein et al .
( 2005 ) , sequences of runs were not stopped when a predetermined number of instances of surfing was reached .
this avoided bias in estimating the probability of surfing , as otherwise a negative binomial distribution would have had to be used to model this probability ( hogg and tanis 1997 ) .
visual inspection of runs in progress gave a strong impression of a dichotomy of outcomes : mutants either vanished or flourished ( see figs 1 and 2 ) . in many runs ,
mutants went extinct . in most runs ( and in nearly all runs on grids of length 400 ) where mutants persisted without advancing , only a few mutants remained on the grid at the end of the simulation , and these few were found close to their point of origin .
longer habitat lengths led to more mutant extinctions and fewer cases of mutant persistence without advancing .
this was apparently because runs on longer grids simply took longer , allowing more time for nonsurfing mutant populations to dwindle and ultimately die out .
wild - type ( circle ) and mutant ( star ) population densities as functions of longitude on a 1  400 grid after 20 , 40 , 80 and 300 generations in a run with mutant surfing .
wild - type ( circle ) and mutant ( star ) population densities as functions of longitude on a 1  400 grid after 20 , 40 , 80 and 300 generations in a run with mutant extinction .
parameters are as in figure 2 . visual inspection also showed that in virtually all runs where mutant surfing occurred , the mutants excluded wild types from the wavefront very soon after the initial mutation arose , and proceeded to nearly monopolize the habitat between the point of mutant introduction and the right - hand boundary .
we thus viewed all successful runs as random trials with two possible outcomes , surfing or mutant extinction . for a given combination of parameters ,
we denote the probability of surfing ( conditioned on nonextinction of the total population ) by psurf .
we note that under neutrality , the left - hand boundary of the habitat occupied by surfing mutants remained close to the point at which the mutation first arose .
surf subsequently lost habitat as wild types advanced to replace the less - fit types .
a major goal of work on mutation surfing is to relate the probability of surfing to model parameters such as cell carrying capacity k , migration rate m and growth rate  ( edmonds et al .
2004 ; klopfstein et al . 2005 ; hallatschek and nelson 2008 ) . for the one - dimensional case ,
plots of mean ( over all runs with given values of  , m and k ) psurf against model parameters showed an interaction between k and  ( fig .
3 ) . for each value of k , psurf increased with . however , the rate of increase was much greater for the smallest carrying capacity , k = 10 , than for the higher values of k. this motivated dividing the runs into two groups , one with k = 10 and one with k = 50 , 100 , 500 , for further analysis .
ln ( psurf ) versus ln ( 1 ) for a one - dimensional habitat , with regression lines for k = 10 ( circles ) and 100 ( squares ) . to quantify the dependence of psurf on k , m and 
, we use logistic regression , a tool which is adapted to processes with binary outcomes for individual trials ( hogg and tanis 1997 ) . in our logistic regression , the response ( or y ) variable is the log - odds of surfing : the predictor ( or x ) variables could in principle be any or all of the model or simulation parameters , or functions of those parameters . in seeking an appropriate statistical model to describe the simulation output , we tried various functions of k , m and  as predictor variables .
we used three different goodness - of - fit tests ( pearson 's chi - squared , deviance and hosmer  lemeshow ) ( hosmer and lemeshow 2000 ) to assess the suitability of each statistical model .
we used three criteria to choose among the models that passed the goodness - of - fit tests at the 1% significance level .
first , we sought a model in which every predictor was associated with an odds ratio that differed significantly from 1 .
third , we sought a model with a low value of the akaike information criterion ( aic ) ( akaike 1974 ) .
we judged the most suitable model to be a function of ln m , ln (   1 ) and ln k : ( standard errors for the coefficients and confidence intervals for odds ratios are given in table 2 ) .
we used mcfadden 's pseudo - r ( r ) to assess the amount of variability explained by equation ( 2 ) ; we found r = 0.31 , indicating reasonably strong explanatory power .
predictors , estimated coefficients with standard errors and 95% confidence interval lower and upper bounds for odds ratios for one - dimensional habitat with k = 10 , 100 , 500 for the same one - dimensional grid with k = 10 , we did not find a satisfactory model .
this may have been because values of psurf with k = 10 were typically very low , so that the relative sampling error in estimating psurf was high enough to obscure associations with the other model parameters ; further simulations are planned to explore this possibility .
it was possible to combine data for all values of k and find a model that satisfied our formal criteria , but this required including so many interaction terms that the model did not seem to be useful for gaining insights .
it is possible to explain heuristically the direction of the dependence of psurf on each predictor (  , m , k ) in equation ( 2 ) .
it is well known for simpler demographic structures that the smaller the total population , the greater is the probability that a mutation initially present in a single copy will rise to high frequencies ( hartl and clark 1997 ) .
correspondingly , in our model the mutant should be more likely to surf if it initially appears in a cell with no other ( wild - type ) occupants .
this is turn should be more likely if the cell carrying capacity k and the migration rate m are low because these conditions should lower the number of individuals at the wavefront , where the mutation occurs .
furthermore , cells at the wavefront should generally contain fewer individuals than cells further back in the wave ; so , fewer offspring will be culled at the front , making the effective growth rate higher there . again , as the mutant must arise at the front , this means that the probability of surfing should increase with  ( as for fixed m and k , a higher value of  should allow more mutants to be born before wild types 
catch up to the mutant population at the leading edge ) . in equation ( 2 ) , it is a straightforward calculus exercise to check that the derivatives   and   are both negative , while   is positive , for the range of parameter values used .
thus , our statistical model predicts that the probability of surfing psurf will decline as m or k increases but rise as  increases , which is consistent with our heuristic argument .
for the ( neutral ) two - dimensional case , we first considered only those runs with a 25  400 grid . using the same criteria as for the one - dimensional case
, we found a simple model to provide the best description of the relationship between psurf and the model parameters : ( standard errors for the three coefficients , in order , were 0.27 , 0.063 and 0.13 ) .
we caution , however , that for this model r = 5% . as in the one - dimensional case , here  
however , including functions of the migration rate m as a predictor did not improve the model in two dimensions .
equations ( not shown ) with coefficients close to those in equation ( 3 ) provided good fits to simulation data with grids of widths 13 or 7 ; the fit was less good for a 3  400 grid .
however , all results for the two - dimensional case should be viewed as preliminary because of the small number of trials performed with each combination of parameters . for the case with selection ,
surfing never occurred when mut was less than 1 ; so , we considered only parameter combinations with mut > 0 ( these did include some deleterious mutations ) .
we did not find a statistical model that simultaneously satisfied goodness - of - fit criteria and provided good explanatory power as measured by r .
the model gave r=0.23 but failed goodness - of fit tests ( p < 0.001 ) ; standard errors for the coefficients , in order , were 0.061 , 0.055 , 0.030 , 0.015 and .0093 .
in equation ( 4 ) , it is noteworthy that psurf depends on mut much more than on wt .
this suggests that for the cases studied , absolute fitness was more important than relative fitness in determining surfing success ( a similar effect was found in an island model analyzed by holt and gomulkiewicz ( 1997 ) ) .
we also note that although the coefficients of ln m and ln k differ from those in the neutral case ( 2 ) , their sign is the same .
thus , for the cases studied , the probability of mutant surfing varies inversely with m and k for both neutral and selected mutations .
beyond merely describing ( even quantitatively ) the dependence of psurf on model parameters , it is clearly important to explain the form of this dependence in a way that sheds light on the process by which a mutation succeeds or fails
. a partial explanation can be obtained from two observations ; we consider the neutral case first . in this case ,
they reproduce and migrate in identical ways , and  count identically toward the carrying capacity of a cell .
the second observation , based on watching generation - by - generation plots of mutant and wild - type population density as a function of longitude , is that in virtually every case when surfing occurs , mutants take the lead and begin to fill the available unoccupied terrain within a very few generations after the original mutation occurs .
as mutant and wild - type individuals behave identically , we may draw a conclusion about the genealogies of individuals in the advancing wave of population density as follows .
imagine drawing a line at the longitude where the mutant first arises , or a short distance past this point .
then when a simulation run terminates , nearly all of the individuals to the right of this line will have descended from a single common ancestor in the generation when the mutant first appeared  regardless of whether mutation surfing occurred .
for when surfing does occur , almost all of the relevant terrain will be occupied by mutants , which must have descended from the single original mutant .
therefore , as mutants and wild types are functionally identical , similar genealogies must arise even when surfing does not occur and the terrain is filled by wild types .
we can thus recharacterize mutation surfing as the scenario in which the common ancestor of nearly all individuals in the population that ultimately invades happens to be the original mutant , rather than any of the wild types present at the time the mutation occurred .
( this line of reasoning informs models of the coalescent with stepping - stone spatial structure ; see , for example , austerlitz et al .
our second observation  that a mutation 's fate appears to be sealed within a few generations of its appearance  suggests a simplified model of the process by which this common ancestor is  
the first step occurs in a early generation ( p , f1 or f2 , where the generation including the original mutant is counted as p ) . in this step ,
one individual is chosen randomly from those in the rightmost cells as a potential common ancestor for the future population .
each individual in the rightmost cell , whether mutant or wild type , has the same probability of being chosen . in the second step ,
the chosen individual succeeds or fails to propagate with probabilities equal to the probabilities of survival or extinction at the very beginning of a run ( when the population consists of exactly one individual , located at the left - hand boundary ) .
motivation for this step comes from the strong association between psurf and prun , the a priori probability that a run will not end in extinction , when  is not very high ( fig .
if the chosen individual fails to propagate ( i.e. its progeny die out ) , the process begins again with a new choice of potential common ancestor from those which were present in the rightmost cells during the p generation ( all of which are assumed to have survived ) . ; ln ( psurf ) versus ln ( prun ) for a one - dimensional habitat and  < 2 .
this model leads to a testable quantitative prediction . to derive an estimate   for psurf under this model , we let n denote the number of individuals present in the rightmost cell in generation p. then the first term on the right - hand side of equation ( 5 ) , ( 1/n)prun , is the probability that the mutant is the first candidate chosen and does succeed in propagating .
the second term is the probability ( n  1)/n(1  prun ) that a wild - type individual is the first chosen but fails to propagate , times the probability 1/(n  1)prun that the mutant is the second chosen and does propagate . the subsequent terms are derived similarly .
algebraic manipulation condenses equation ( 5 ) into the form to test whether equation ( 6 ) was a good fit to the one - dimensional simulation results , we first estimated prun and mean n separately for each combination of the parameters  , m and k that yielded at least one run not ending in extinction ( there were 70 such parameter combinations ) , then ran a linear regression of the observed values of ln ( psurf ) against the predicted values   from equation ( 6 ) . ( taking logarithms rendered the variance of observed psurf approximately equal for all values of , a necessary condition for valid inference in linear regression . )
we caution that in our simulations when  was very high ( i.e. 8 ) , prun was always very close to 1 and so a relationship between prun and psurf was no longer apparent ; therefore , the following discussion considers only our neutral simulations with  < 2 .
the regression equation was ( r=56.1% ) , which is broadly consistent with prediction ( 6 ) , except that the constant was lower than expected ( standard errors for the constant and the coefficient of ln ( prun ) were 0.1766 and 0.1120 respectively ) .
when the regression was run for each value of the migration rate m separately , the discrepancy between predicted and observed psurf was seen to occur only for very high migration rates .
indeed , the constant term and the coefficient of   in the regression equation were not significantly different from the predicted values of 0 and 1 , respectively , except for the highest value of m , namely 0.40 .
agreement between the predicted and observed equations improved as m decreased ( details are shown in table 3 ) .
coefficients for regression of observed on predicted ln ( psurf ) we now briefly consider model ( 6 ) for large . when  and k are both very large , prun approaches 1 , and taking a limit of the right - hand side of equation ( 6 ) with k being a fixed multiple of   1 ( k = a(  1 ) ) and prun = 1  e suggests that psurf should approximately equal 1/. however , our simulations with large  and k = (   1)/2 show far higher values of psurf than this calculation would predict  for example , psurf  0.85 with  = 32 , k = 16 and m = 0.1 .
the explanation appears to be that even though most cells filled to capacity within one or two generations of colonization , the rightmost cell generally contained very few individuals  averaging fewer than two for all conditions simulated  and each of these individuals had a high probability of becoming the common ancestor of the future occupants of the grid from that point onwards .
( it is possible that different behavior might be found if k were held constant as  increased , however . )
analogues of equation ( 6 ) proved to be almost completely unexplanatory in two - dimensional settings under neutrality , yielding r values of less than 5% .
indeed , there appeared to be no meaningful relationship between any powers of psurf and prun in two dimensions . under selection
( i.e. with mut  wt ) , the analogue of equation ( 6 ) is more complicated because the numbers of mutant and wild - type  candidates remaining after each unsuccessful attempt to propagate do not follow a simple linear pattern . rather ,
if the rightmost cell in generation p is occupied by one mutant and n  1 wild types , and if the number of candidates is still assumed to decrease by 1 after each failed attempt , then the expected numbers u(t ) and w(t ) of mutant and wild - type candidates t generations later ( conditioned on mutant survival ) are given by when these expressions are used in an analogue of equation ( 6 ) that also includes different values of prun for mutants and wild types , the result is not readily simplified . however , the first two terms can be computed : where   and   are the values of prun for mutants and wild types , respectively , and ti is the probability of mutant surfing , conditioned on the event that the successful surfing attempt occurs i  1 generations after generation p. for the parameter combinations used in our simulations , mean t1 performed very well as a predictor of mean psurf for beneficial mutations ( r = 96.8% ) , and mean t1 + t2 performed slightly better ( r = 98.5% ) .
specifically , linear regression yielded the equations with standard errors of 0.028 for the constant and 0.056 for the coefficient of , and with standard errors of 0.019 for the constant and 0.037 for the coefficient of t1+t2 .
as equations ( 10 ) and ( 11 ) show that t1 and t2 are linear functions of   and , this indicates that psurf itself is very close to a linear function of these two quantities .
for deleterious mutations , however ,   and t1+t2 drastically overestimated psurf , possibly because deleterious mutations that initially  surfed later lost ground to introgressing wild types .
we have developed descriptive statistical models that quantify the relationship between the probability psurf that a mutation will surf and the parameters  , m and k of a frequently studied cellular automaton model ( edmonds et al .
we have also proposed an analytic model that performs reasonably well at predicting psurf for the case of a neutral or beneficial mutation on a one - dimensional grid .
such simplified models are desirable because a completely rigorous explanation of the dependence of surfing probability on a full range of demographic and genetic parameters currently appears unobtainable .
( this is not to say that the mathematical models that have been developed , such as those in hallatschek and nelson 2008 , 2009 , are uninformative , but simply that they do not yield self - contained formulas for psurf as a function of system parameters in all regimes . ) moreover , a typical rigorous result for an individual - based model of invasion dynamics may simply give conditions under which the probability of persistence is either positive or zero ; see , for example , the recent result on surfing in e. andjel , j. miller and e. pardoux , unpublished manuscript . in the absence of a rigorous theory , descriptive statistical models and simplified analytic models like those developed here can assess the strengths or weaknesses of heuristic explanations and may be of value in making predictions of invasiveness .
thus , such models remain crucial tools for checking and expanding upon nonquantitative verbal models of mutation dynamics during invasions . for the neutral case
, our statistical models differ in several respects from those proposed in earlier work ( klopfstein et al . 2005 ) .
one reason may simply be that previous authors used only grids of length 100 cells for their simulations ( edmonds et al .
2004 ; klopfstein et al . 2005 ; travis et al . 2007 ; hallatschek and nelson 2008 ) .
this may have inflated estimates of the probability of surfing , and especially of mutant survival without surfing ( travis et al .
2007 contains some discussion on this point ) . with the grids of length 400 used here , it is extremely rare for mutants to survive to the end of a simulation unless surfing has occurred . as a result , the bimodality of the final spatial distribution of mutants reported in edmonds et al .
some experimentation ( details not shown ) suggested that extension to still longer grids would have little effect on observed probability of mutant survival and psurf .
we note that with the bimodality removed , it is easy to explain the  midpoint rule for locating the point of origin of a mutation proposed in edmonds et al .
visual inspection of simulations in progress shows that surfing neutral mutants almost invariably occupy nearly all of the terrain from their point of origin to the rightmost boundary of the grid .
it follows that the centroid of the final mutant population is halfway between the point of origin and the boundary , which is a rephrasing of the midpoint rule .
in addition , as we have noted above , when beneficial mutants surf they also introgress backwards into territory previously colonized by wild types . the extent to which weak selective advantage or linkage to selected loci may bias midpoint - rule estimates therefore deserves study . more substantively ,
it is important to address why equations ( 2 ) and ( 3 ) , which we propose as statistical descriptions of the relationship between the log - odds of surfing   and the parameters  , m and k , differ dramatically from the simple linear relationship between  mutant success ( a proxy for psurf ) and the  lumped parameter (   1)/km proposed in klopfstein et al .
( 2005 ) ( which , we note , is also consistent with the heuristics presented above ) . indeed ,
equation ( 2 ) suggests that more appropriate lumped parameters for the one- and two - dimensional cases , respectively , might be ( 1)/km and (   1)/k .
( 2005 ) depends on a linear regression of psurf on (   1)/km .
linear regression is generally not appropriate for response variables , like psurf , that are constrained to lie between 0 and 1 ( hogg and tanis 1997 ) .
( 2005 ) it appears that much of the apparently high r attributed to this regression is due to one or a very few influential data points .
thus , we believe that equations ( 2 ) and ( 3 ) are better justified as statistical descriptions of the dependence of   ( and hence psurf ) on  , m and k. however , perhaps the most noteworthy difference between our statistical model ( 3 ) for the two - dimensional case and that proposed in klopfstein et al .
( 2005 ) is that equation ( 3 ) does not include the migration rate m. in fact , a logistic regression of   on ln ( m ) for a 25  100 grid did show a marginally statistically significant inverse dependence [ , with a standard error of 0.1077 for the coefficient of ln ( m ) ( p = 0.090 ) ] .
( 2005 ) , which was not quantified . however , the important point is the weakness of the dependence .
this is also a major difference between equation ( 3 ) and our proposed model ( 2 ) for the one - dimensional case .
it is reasonable to ask what factors account for the decreased importance of m in the two - dimensional setting .
one such factor may be the irregular profile of the wavefront in two dimensions ( see fig .
first , individuals in cells surrounded by empty or nearly empty cells ( i.e. where the wavefront has positive curvature ) are expected to have more surviving offspring than individuals in cells surrounded by relatively full cells ( i.e. where the wavefront has negative curvature ) .
thus , the probability of surfing will be greater for mutations that happen to arise in reasons of positive curvature .
second , the irregular profile will result in competition between outcroppings of the wavefront at different latitudes in the grid .
unlike in the one - dimensional case , an outcropping of mutants at one latitude may be rivalled by an outcropping of wild types at another latitude .
the mutant and wild types will eventually meet , and only one type will then fill the remaining open terrain ( this has been discussed in detail in hallatschek et al .
thus , the irregular , two - dimensional wavefront introduces two additional random factors affecting the probability of surfing that do not have clear relationships to the migration rate m. contour plot of total population density in part of a 25  100 grid .
black region contains 0 individuals per cell , white region contains eight or more individuals per cell .
parameter values are  = 1.8 , m = 0.4 , k = 10 . a point related to the relative unimportance of m in two dimensions concerns equation ( 6 ) , which predicts the probability of surfing in the one - dimensional case .
as noted in the analytic models section , the model performed better at smaller values of the migration rate m. this may be because the model ignores all individuals present on the grid except for those in the leading cell . in reality , individuals behind the wavefront do have some chance of filling the remaining open terrain with their progeny . to do so , however , their progeny must advance to the wavefront , and the likelihood of this occurring should decrease with m. an analogy may be drawn to solutions of fisher 's equation , a classical model of travelling waves of population density ( fisher 1937 ) .
the diffusion coefficient d in equation ( 14 ) is one - half the mean - squared distance traveled by an individual during one generation ; so , for us d = m/2 .
it is straightforward to check that near the leading edge , traveling wave solutions of equation ( 15 ) have shapes proportional to the exponential function .
thus , the number of individuals originally present one cell behind the leading cell ( say ) whose progeny will travel forward a given distance in a given short span of time is proportional to .
this quantity increases with m. therefore , referring back to our simplified model , when an individual originally in the leading cell fails to propagate , we expect fewer candidates from cells further back to be available at the wavefront as potential propagators . as a result , ignoring all cells other than the leading cell makes equation ( 6 ) a better approximation of psurf when m is small .
the simulations in the present work are necessarily limited in scope ; even for the neutral case , they do not represent all possible scenarios of biological interest . in our simulations ,
a change in the ordering of life cycle events may change the conclusions we draw from the simulations , and it will therefore be necessary to vary this ordering in future work . likewise , our simulations assume nearest neighbor migration in a stepping - stone model with a selection regime that is spatially and temporally invariant . if such a scenario is used as a null model for hypothesis tests regarding a population 's migration history , other spatial structures and migration patterns may need to be considered as alternatives . the present work modeled only mutations arising at the leading edge of an invasion .
some authors have modeled mutations arising further back in an invasion as well ( klopfstein et al . 2005 ; travis et al .
( 2007 ) that the probability of surfing for the model considered here drops off sharply for neutral mutations arising even a few cells behind the front ; in burton and travis ( 2008b ) , it was shown that this holds true for deleterious and , to a lesser extent , beneficial mutations as well ( but different boundary conditions can alter this behavior ) .
as mutants and wild types behave identically in our neutral model , this helps explain why our simplified analytic model achieves good predictive power , even though it assumes that the probability of propagation is zero for individuals behind the front . in all the simulations discussed here ,
a single mutant was introduced at a fixed location once the wild - type invasion had advanced to the 10th grid column from the left .
one could ask how varying mutation times and locations might affect the results . as extinction almost never occurred
after the mutant was introduced , we suspect that increasing the longitude of the mutation would have little effect on psurf ; however , this remains to be checked . on the other hand , changing the model to include multiple mutations , perhaps at random times , would shed light on a number of questions .
banding of the colonized territory by genotype ( r. holt , personal communication ) , allowing for more precise tracking of a population 's migration history and possibly providing the basis for a null distribution in tests for selection at other loci .
simulations with serial selected mutations could be used to study adaptive changes in quantitative traits during invasion . among other questions
, they could be used to study the distribution of effects of quantitative trait loci ( qtl ) in a trait subject to selection during an invasion , extending existing theory which does not take spatial population dynamics into account ( mackay 2001 ; orr 2005 ) . as distance from the front
strongly influences the probability of surfing for neutral mutations , one might ask why we have chosen to study only mutations arising at the front .
the explanation is that we have done so precisely in order to focus on what other factors may affect the probability of surfing .
analogously , one might study a group known to be at high risk of contracting a certain disease , in order to learn what additional factors might predispose an individual to contract the disease or to stay healthy . on the other hand , hallatschek and nelson have shown that as population density is ( almost by definition ) low at the wavefront , if mutation rates are constant over space then the modal point of origin of surfing mutations will be slightly behind the front ( hallatschek and nelson 2008 ) .
thus , studies of adaptive genetic change during invasions via the accumulation of weakly beneficial mutations need to vary their points of origin ; a two - locus model which includes this variation has already been studied in burton and travis ( 2008a , b ) . as have previous studies of mutation surfing ( edmonds et al .
hallatschek and nelson 2008 ) , the present work considers only nearest - neighbor migration , in which an individual can move at most one cell length per generation .
this results in diffusive movement , which can be modeled deterministically by a reaction  diffusion equation ; such equations are the basis for the analysis of surfing in vlad et al .
however , it is well established that models assuming diffusive movement can yield inaccurate predictions when applied to many real invasions .
this may be especially true when the movement rate m is relatively high for a given cell carrying capacity k. not every combination of m and k may be well described by a limiting process on a continuous spatial domain ( see , for example , durrett and neuhauser 1999 ) .
a future task will be to assess , in the first instance through simulation , how patterns that have been observed for discrete mutation surfing models converge or fail to converge to limiting relationships as the scaled island model approaches a continuum .
for accurate modeling ,  heavy - tailed dispersal kernels that model a greater frequency of long - distance dispersal events are often necessary ( see , for example , kot et al . 1996 and references therein ) .
an important area for future work will be to extend studies of mutation surfing to invasion models that incorporate heavy - tailed kernels ( e.g. laplace , or reflected exponential , kernels ) and rare long - range dispersal events ( this was noted in burton and travis 2008b as well ) . yet , as stochastic effects inevitably will be even more important in such models , it may be extremely difficult to develop a rigorous theory that will yield useful quantitative predictions of psurf under these conditions .
simplified models that yield good agreement with simulation and , ultimately , experimental data will be especially necessary in these settings .
the bulk of the present work concerns neutral mutations , but we have obtained some results for beneficial and deleterious mutations ( in one - dimensional habitats ) as well . equation ( 4 ) describes the dependence of psurf on m , k and nonidentical reproduction rates mut and wt .
as noted above , the direction of the dependence on m and k was the same under neutrality and selection .
in addition , the psurf was much more sensitive to the value of mut than to wt , so that absolute fitness was a better predictor of surfing success than relative fitness .
consistent with this finding , equation ( 13 ) , which includes   as a predictor of psurf , performs only slightly better than equation ( 12 ) , which does not . as it is difficult to identify consistent predictors of establishment success from empirical studies ( kolar and lodge 2001 ; hayes and barry 2008 )
, one might ask why analytic models that include prun ( whether for mutants or wild types ) as a predictor of psurf are informative .
first , our results indicate that finding the determinants of surfing and of invasiveness may not be two separate hard problems but only one .
second , they lead to a novel and ( in principle ) empirically testable prediction about adaptive change during invasions , as follows .
our simulations indicate a strong positive association between the probability of establishment ( which in this model leads to invasion )   and the probability psurf of surfing , conditioned on establishment , for neutral and beneficial mutations in one dimension .
( however , this association was not found for parameter regimes in which   was very close to 1 . ) in these simulations , the only difference between mutants and wild types was fecundity , as measured by the growth rate . for most of the parameter combinations studied here , when m and wt are held fixed and cell carrying capacity k and migration rate m are varied , higher   is associated with higher psurf .
this suggests that in the scenarios modeled here , conditions that favor the establishment of a population in a novel habitat might in and of themselves also favor the fixation of beneficial mutations once establishment has occurred and an invasion has begun .
it has been noted previously that high intrinsic growth rates and small deme sizes at the front may promote surfing , and that 
genetic revolutions during invasions can occur because of surfing due to increased drift at the sparsely populated front ( klopfstein et al . 2005 ; travis et al .
2007 ; excoffier and ray 2008 ; hallatschek and nelson 2008 , 2009 ) , but the finding that factors other than  and population density can promote surfing in some invasions is novel , to our knowledge . as k and m do not depend solely on organismal traits ( whether genetic or nongenetic ) but also on features of a given habitat , our results suggest that even abiotic factors that promote establishment might also promote the survival and spread of beneficial mutations ( as well as neutral ones ) , conditioned on establishment .
this finding is complementary to studies that have identified propagule size and number of introductions as the factors with the strongest currently known association with establishment success ( kolar and lodge 2001 ; hayes and barry 2008 ) .
further simulation and empirical studies seem warranted to assess the generality of this effect and its importance in explaining adaptive genetic change during invasions .
 when a neutral mutation arises in an invading population , it quickly either dies out or  surfs 
 , i.e. it comes to occupy almost all the habitat available at its time of origin . 
 beneficial mutations can also surf , as can deleterious mutations over finite time spans . 
 we develop descriptive statistical models that quantify the relationship between the probability that a mutation will surf and demographic parameters for a cellular automaton model of surfing . 
 we also provide a simple analytic model that performs well at predicting the probability of surfing for neutral and beneficial mutations in one dimension . 
 the results suggest that factors  possibly including even abiotic factors  that promote invasion success may also increase the probability of surfing and associated adaptive genetic change , conditioned on such success . 