Two models for speciation via selection have been proposed. In the well-known model of ‘ecological speciation’, divergent natural selection between environments drives the evolution of reproductive isolation. In a second ‘mutation-order’ model, different, incompatible mutations (alleles) fix in different populations adapting to the same selective pressure. How to demonstrate mutation-order speciation has been unclear, although it has been argued that it can be ruled out when gene flow occurs because the same, most advantageous allele will fix in all populations. However, quantitative examination of the interaction of factors influencing the likelihood of mutation-order speciation is lacking. We used simulation models to study how gene flow, hybrid incompatibility, selective advantage, timing of origination of new mutations and an initial period of allopatric differentiation affect population divergence via the mutation-order process. We find that at least some population divergence can occur under a reasonably wide range of conditions, even with moderate gene flow. However, strong divergence (e.g. fixation of different alleles in different populations) requires very low gene flow, and is promoted when (i) incompatible mutations have similar fitness advantages, (ii) less fit mutations arise slightly earlier in evolutionary time than more fit alternatives, and (iii) allopatric divergence occurs prior to secondary contact.
Natural selection commonly promotes speciation [1–5]. This statement is supported by widespread correlations between levels of ecological divergence and reproductive isolation [3,6], cases of divergence with levels of gene flow too high to allow genetic drift , examples of traits under selection that also cause reproductive isolation , and theoretical studies [1,2,9]. More specifically, reproductive isolation can arise if the different alleles favoured by selection within each population are incompatible with one another when brought together in the genome of a hybrid [1,2]. Indeed, negative interactions between alleles at different loci have been shown to result in hybrid inviability and sterility (Dobzhansky–Muller incompatibilities denoted as DMIs hereafter; [2,10,11]) and molecular studies have demonstrated that selection can drive the evolution of the genes involved [1,12].
However, different types of selection might drive the evolution of reproductive isolation. Consequently, models of speciation via selection have recently been classified into two main categories: ecological speciation and mutation-order speciation [4,13]. In ‘ecological speciation’, divergent natural selection between different environments drives the evolution of reproductive isolation. This model has received substantial support from laboratory experimental evolution studies ( for review), case studies in nature [15–18] and comparative analyses .
In ‘mutation-order speciation’, different and incompatible mutations (alleles) fix in different populations adapting to the same selective pressure (i.e. uniform selection, following Mani & Clarke ; reviewed in Schluter ). The order in which new mutations arise can affect which mutations fix, making ‘mutation-order’ a potentially important stochastic process in evolution. This mechanism has also been referred to within the more inclusive term ‘non-ecological speciation’, but as noted by Sobel et al. , uniform selection can still involve ecology. While laboratory experimental evolution studies have demonstrated that uniform selection on isolated populations can lead to divergence (e.g. [20,21]), mutation-order speciation is less substantiated than ecological speciation, in part because how to demonstrate it in nature is unclear . Thus, it is important to develop predictions for when mutation-order speciation is likely and tests that can distinguish it from ecological speciation. This is the aim of the current article. We focus on how five factors impact mutation-order speciation: (i) rates of gene flow, (ii) the degree of hybrid incompatibility, (iii) the relative fitness of different new mutations, (iv) the timing of origination of different mutations, and (v) an initial period of allopatric divergence preceding secondary contact and gene flow. We consider reproductive isolation (i.e. low hybrid fitness) arising from both within-locus and between-locus incompatibilities. Population differentiation in genes reducing hybrid fitness is our measure of divergence via the mutation-order process, with divergence of more strongly incompatible alleles representing the evolution of stronger reproductive isolation.
It has been argued verbally that mutation-order speciation is difficult, and can perhaps be ruled out, when gene flow occurs [4,13]. The logic is as follows. Suppose that two new adaptive but incompatible alleles arise in two different populations. The alleles will move between populations via gene flow, and because one allele will almost always have at least a slightly higher selective advantage than others, a likely outcome is that the same, most advantageous allele will fix in both populations. Additionally, with gene flow, incompatible genotypic combinations will be created (i.e. ‘hybrids’), and alleles that are part of this combination will be removed by selection . These processes should prevent different incompatible alleles from fixing in different populations which experience gene flow. However, based on verbal arguments alone, it is unclear what level of gene flow would prohibit divergence. It is also unclear how evolutionary factors other than gene flow affect mutation-order speciation, and how different factors interact. In summary, clearer predictions are required about mutation-order speciation.
Previous theoretical work has examined DMIs and mutation-order divergence in a few important ways. Seminal studies demonstrated that pairwise interactions between loci accumulate in a nonlinear fashion, leading to a ‘snowball’ of DMI accumulation [10,11]. However, these studies did not consider gene flow or specify the type of selection driving divergence. Some considered the importance of mutation-order when speciation is driven solely via drift [22–25]. The effect of mutation-order during divergence via selection on new, advantageous alleles in a uniform environment has also been examined in some contexts. Cohan  showed that uniform selection in isolated populations accelerates divergence relative to that observed via drift alone. An insightful paper by Mani & Clarke  used simulations to establish mutation-order as an important stochastic process in evolution, but one distinct from genetic drift, because mutations fix via selection. Similarly, Unckless & Orr  derived analytical expressions for probabilities of the evolution of a DMI when two populations experience identical environments. However, all these studies examined populations that experience no gene flow.
Other models examined the mutation-order process with gene flow [2,28–31]. For example, Gavrilets presents a simple model of mutation-order speciation (, pp. 139–140) and reviews more complex spatial models of the process (, pp. 221–228). The main conclusion to date is that if there is gene flow, then ‘Once established in a local population, an advantageous allele will tend to spread across the whole system’ (, p. 225). Previous models, however, have not ‘competed’ new mutations against one another in the face of gene flow, varying their selective advantage and timing or origination, to test the factors promoting and constraining mutation-order speciation.
The current study extends useful past models in at least three general ways. First, we vary the fitness of new alleles and hybrid genotypes relative to each other and relative to ancestral alleles and genotypes. Second, we vary the timing of origination of selectively advantageous alleles relative to each other. Third, we consider the effects of initial periods of allopatric divergence preceding the onset of gene flow. We use the results to develop predictions about the conditions that affect mutation-order speciation.
2. General methods
We used simulation models to study populations of diploid, sexually reproducing organisms with discrete, non-overlapping generations. We assume that individuals are divided between two identical patches, and we refer to each patch as a population. Within populations, we assume random mating with respect to genotype. Our first model considered a single locus (model 1: ‘single locus model’). The second model considered epistasis between two completely unlinked loci (i.e. a DMI; model 2: ‘epistasis model’). Thus, both models study the origin of an incompatibility, whether it leads to divergence, and if so, what degree of divergence. Although speciation probably often involves divergence at more than two loci , studying a two-locus model is the most logical starting point for a multi-locus model because: (i) the first alleles to diverge are those that initiate the speciation process, and (ii) once one incompatibility is established, the establishment of future incompatibilities via a ‘snowball’ process [10,11] becomes possible. Because qualitatively similar results were observed between models, we present here only the results for the more generally accepted scenario of a DMI (our ‘epistasis model’). Results for the single-locus model are presented in the electronic supplementary material. Unless otherwise noted, descriptions below apply to the epistasis model.
Both models considered the fate of two new mutations that were selectively advantageous, as might be expected for speciation with gene flow [2,32,33]. Thus, we assume no genetic drift. The possible genotypes and their fitnesses are given in table 1. Within a patch, we assume Hardy–Weinberg genotype frequencies after random mating has occurred. Given these assumptions, we can use very straightforward equations (described below) to model allele frequency change from one generation to the next.
The general scheme is as follows (figure 1). Populations are initialized with a single recessive allele initially fixed at each locus of interest. The genotype with only the ancestral, recessive allele(s) is assigned a fitness of 1. We study the spread of two new mutations (alleles), each of which has higher fitness than the ancestral allele but is also incompatible, to some degree, with the other new allele (table 1). When a new allele arises, it is assumed to be dominant to the ancestral allele at its locus, and initially arises with a starting frequency of 0.001 in one of the patches (and 0 in the other patch). Each simulation run then proceeds with the iteration of a two-step cycle for 107 generations or until allele frequencies have stopped changing in both patches (whichever comes first).
In the first step of the cycle, natural selection adjusts allele frequencies according to the selection coefficients for each genotype. In the second step of the cycle, migration between patches occurs. We assume that migration is random with respect to patch and genotype, and we vary the migration rate between simulation runs. Thus, a fraction m of each genotype was moved from each patch to the other patch. Fixation of an allele in a patch was defined to have occurred if the allele reached a frequency greater than 0.9999 in that patch. To determine if ending simulations at 107 generations affected our results, we examined the final 100 timesteps of thousands of cases that reached 107 generations, and in every case the allele-frequency changes were extremely small in magnitude (less than 10−15). That is, the system had essentially reached equilibrium before the simulation was ended, and only infinitesimal changes (near the limit of computational precision for double-precision floating-point decimal numbers) were occurring. All simulations were programmed in C and run on a 64-bit Intel Macintosh computer. Source code and spreadsheets with full results are available upon request.
3. Between-locus incompatibilities: the epistasis model
Here we consider two incompatible alleles that arise at separate loci and interact to reduce hybrid fitness. The ancestral alleles at the two loci are designated a and b, and the two new adaptive but incompatible alleles are A (arises at the locus where a is ancestral) and B (arises at the locus where b is ancestral). Since we are interested in the divergence of populations, we assume that the new alleles always arise in different patches. The patch in which B arises is designated patch 1, and the patch in which A arises is designated patch 2. We let pi represent the frequency of allele A in patch i (i = 1 or 2) and qi represent the frequency of allele B in patch i. Thus, the frequencies of alleles a and b in patch i are, respectively, ri = 1 − pi and si = 1 − qi. We are interested in scenarios in which both new alleles are advantageous when compared with ancestral alleles, but one new allele is more advantageous than the other. Thus, AAbb (or Aabb) was the most favoured genotype and aaBB (or aaBb) was more favoured than aabb, but less favoured than AAbb (table 1 for full fitness scheme).
Allele frequencies changed from the current generation to the next as follows. With the assumptions stated above, we calculate the frequency of A in patch i in the next generation as its proportional weighted fitness (e.g. , p. 42): 3.1 where wAAbb is the fitness of an individual of genotype AAbb, wAabb is the fitness of an individual of genotype Aabb, and so on. The frequency of a in patch i in the next generation was simply 1 − [equation (3.1)]. The frequencies of B and b in the next generation were calculated in an analogous manner.
To explore the effects of the relative timing of origination of the two alleles, we allowed allele B (the less fit of the two new alleles) to arise in patch 1 (with initial frequency q1 = 0.001) and spread to a frequency of qthresh in at least one of the patches, while holding p1 = p2 = 0. In essence, we gave the less advantageous mutation a ‘leg-up’, such that it could spread to various frequencies before the second, more advantageous mutation originated. Following this time-lag period set by qthresh, we then allowed allele A to arise in patch 2 with a frequency p2 = 0.001. In addition to adding biological realism, this feature addresses questions about whether newer alleles that are better adapted can displace incompatible alleles that already have a foothold in a population. We varied qthresh over several values (0, 0.01, 0.1, 0.333333, 0.5 and 0.9), with qthresh = 0 corresponding to simultaneous origination of both mutations.
To explore the effects of an initial period of allopatric divergence preceding gene flow, we ran simulations as above, but set m to zero for the first tcontact number of generations. We varied tcontact over a wide range of values (0, 10, 25, 30, 40, 75, 100, 1000, 10 000). We focused on analysing the role of allopatry using the parameter combinations of qthresh and hybrid incompatibility that yielded the most divergent predictions about evolutionary outcomes (qthresh = 0, wAABB = 0.9 and qthresh = 0.333, wAABB = 0.001, see §4). Exploration of a wider range of parameter values revealed that these two sets of parameter combinations accurately encapsulated the range of effects of allopatric divergence (electronic supplementary material, figures S4 and S5).
Two main sets of simulation runs were conducted that we present here. The first set was used to explore large, biologically realistic areas of parameter space. Simulations in this set were run with all combinations of the following ranges of parameter values (with the constraint that only combinations with wAAbb > waaBB were used): wAAbb = [1.005, 1.01, 1.02, 1.03, 1.04, 1.05, 1.1, 1.15, 1.2, 1.25, 1.3, 1.35, 1.4, 1.45, 1.5]; waaBB = [1.001, 1.005, 1.01, 1.02, 1.03, 1.04, 1.05, 1.1, 1.15, 1.2, 1.25, 1.3, 1.35, 1.4, 1.45, 1.49]; wAABB = [1, 0.999, 0.995, 0.99, 0.98, 0.97, 0.96, 0.95, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1, 0.05, 0.01, 0.005, 0.001, 0]; m = [0, 10−5, 10−4, 10−3, 10−2, 0.0125, 0.025, 0.05, 0.1]; qthresh = [0, 0.01, 0.1, 0.333, 0.5 and 0.9]; and tcontact = 0. For the most part, changing parameter values had moderate, quantitative effects on the results, with the migration rate and ratio of mutant fitnesses having the largest effects. As such, we show a subset of data that illustrate the main trends.
We then performed a second set of simulations that were used to gain higher resolution in regions of parameter space that illustrate the effects of each parameter. For this second set, we fixed waaBB = 1.2 and then varied wAAbb over 100 evenly spaced values between 1.206 and 1.8. This varied the ratio of fitnesses, wAAbb/waaBB, between 1.005 and 1.5 in increments of 0.005. Simultaneously, we varied m over 100 logarithmically spaced values between 10−5 and 10−1. This resulted in a set of 10 000 simulation runs which we display together in a single phase diagram. Different phase diagrams show the results for different combinations of the values of qthresh (0, 0.01 or 0.333), tcontact (0, 10, 25, 40 or 75 generations) and wAABB (0.001 or 0.9).
Using the following classification scheme, we show the locations of boundaries between different evolutionary outcomes. We categorized our results into five outcomes, which correspond to the differently shaded regions shown in phase diagrams (M.O. = mutation-order divergence): (i) ‘extreme M.O.’ = when the two alternative, incompatible alleles each reached and maintained a frequency greater than 0.99 in their respective patches, (ii) ‘moderate M.O.’ = when the alternative, incompatible alleles reached and maintained frequencies greater than 0.9 (but one or both alleles had a frequency less than 0.99) in their respective patches, (iii) ‘weak M.O’ = all cases when A and B were both maintained but divergence was weaker than for ‘moderate M.O.’, (iv) ‘A fixed’ when allele A was fixed in both patches and allele B was eliminated, and (v) ‘B fixed’ when allele B was fixed in both patches and allele A was eliminated. We stress that outcomes (i)–(iii) are not qualitatively different: they all represent population divergence via the mutation-order process. Nonetheless, categorization into three levels is useful for visualizing and discussing the magnitude of population divergence observed, in each case occurring in the absence of ecologically based divergent selection. By contrast, outcomes (iv) and (v) represent cases in which there is no divergence.
4. Epistasis model: results
At least some population divergence via the mutation-order process was observed under a relatively wide range of conditions (figures 2⇓–4; electronic supplementary material, figures S4 and S5). However, strong divergence was unlikely unless gene flow was absent or extremely low. Specifically, moderate M.O. and extreme M.O. were never observed at migration rates above 0.01 and in general consistently occurred only when m < 0.001.
Three factors other than gene flow strongly affected the probability of divergence via the mutation-order process. First, as the difference in selective advantage of the two new mutations decreased, mutation-order divergence became more likely (figures 2⇑–4). This was because small differences in relative fitness caused the most fit allele, A, to have a smaller chance of fixation in both patches.
Second, staggering the timing of origination of mutations affected the outcome, but its effect depended upon other parameters. When the first, less-fit mutation (allele B) originated slightly before the second, mutation-order divergence was promoted; for example, being observed when differences in the selective advantage of new mutations were relatively large. By contrast, when the first mutation, B, originated long before the second mutation (i.e. was given a large ‘leg-up’), the less advantageous allele could actually spread between and even fix in both patches, preventing mutation-order divergence. In part, the latter was owing to the second mutation, A, often finding itself in a low fitness hybrid (AaBb, AaBB, AABb or AABB). The latter effect could only be ameliorated by extremely low migration rates that prevented allele B from rising in frequency in patch 2.
Third, an initial period of allopatric divergence tended to promote mutation-order divergence, which can be seen by noting that all three M.O. regions expand as one moves left to right across the panels in figure 4. Just 25–75 generations of allopatric divergence promoted mutation-order divergence. Periods of allopatry longer than 75 generations has little additional effect (e.g. results for 100–10 000 generations of allopatry are very similar to those for 75 generations; electronic supplementary material, figures S4 and S5).
In contrast to the effects of the other factors, hybrid fitness had relatively weak effects: greater incompatibility caused only slight increases in the conditions favouring mutation-order divergence. Nonetheless, as wAABB decreased from 0.9 (moderate incompatibility; figure 3a,b) to 0.001 (nearly complete incompatibility; figure 3d,e), the M.O. regions all expanded somewhat (figure 3). Overall, these results were qualitatively very similar to those observed for the single locus model (for details see electronic supplementary material, figures S1–S3).
Previous considerations of mutation-order speciation have primarily focused on a single prediction: divergence should be unlikely if there is gene flow . However, there has been a lack of specificity about what level of gene flow would prevent mutation-order divergence, and a lack of predictions about the influence of factors other than gene flow. While our models partly confirm intuitive effects of gene flow (i.e. more gene flow = less divergence), we also show that even with moderate gene flow, the mutation-order process can cause DMIs to become established, thereby beginning the process of population divergence, and potentially sowing the seeds for a ‘snowball’ [10,11] of DMI accumulation. Furthermore, our model offers straightforward predictions about how four other factors affect the likelihood of mutation-order speciation. Here, we discuss the predictions our model generates, how these can be tested empirically, and extensions for future work.
(i) Gene flow
The mutation-order process can lead to at least some divergence of populations even when there is gene flow, but strong population divergence requires little or no gene flow (figures 2⇑–4). There is an upper limit of gene flow beyond which any divergence at all via a mutation-order process is unlikely. This is seen by noting the presence of a boundary on the right-hand sides of the different M.O. regions of figures 3 and 4. Finally, previous verbal arguments about mutation-order speciation (discussed above) suggested that high gene flow would lead to fixation of the most fit allele. Our model shows that this is sometimes, but not always, true. While high gene flow will lead to fixation of a derived allele (rather than divergence of populations), it is not necessarily the most fit allele that fixes in such cases; which allele fixes depends upon other factors (predictions below).
(ii) Relative fitness of incompatible alleles
Not surprisingly, increasing the ratio of fitnesses (wAAbb/waaBB) favoured allele A, causing (i) a transition from mutation-order divergence to complete fixation of A (figure 3a,d; qthresh = 0, i.e. no lag), and (ii) a transition from fixation of allele B to mutation-order divergence (figure 3c,f; qthresh = 0.333). Thus, mutation-order divergence will be more likely when new, incompatible alleles have more similar fitnesses (see also Unckless & Orr ).
(iii) Staggered emergence of new alleles
An intermediate length of time between mutations is the most likely to produce mutation-order divergence (figure 3b,e). When there was no delay, much of the M.O. regions were replaced by the complete fixation of the most favoured allele, A. On the other hand, when the delay was much longer, the M.O. regions were eroded by complete fixation of B. In general, by arising somewhat earlier a less-fit allele could often exclude another, better adapted allele.
(iv) Period of allopatry
Coyne & Orr  and Barton  suggested that allopatric divergence followed by secondary contact could facilitate mutation-order divergence. Our models support this suggestion. In fact, if the allopatric period lasted enough generations, the influences of all parameters except migration rate were nearly eliminated. Intuitively, this is explained by the fact that when tcontact is large, each of the derived alleles has (nearly or completely) gone to fixation in its original patch before any migration occurred; the outcome after secondary contact is then a consequence of migration-selection balance. While this might be expected via intuition alone, a less obvious finding is that relatively short periods of allopatry could promote mutation-order divergence, with longer periods having little additional effect on promoting divergence. However, a long enough period of allopatric divergence eliminated the chance that the less fit of the two new alleles could go to fixation. Periods of allopatric divergence are thus predicted to counter the effects of long time lags between mutations. Some examples of ecological speciation between currently sympatric taxa involved an allopatric element [35,36], but data from cases of mutation-order divergence are required to test the role of allopatry.
(v) Incompatibility between alleles
Somewhat surprisingly, the degree of incompatibility of the two new mutations had effects that were quite small in magnitude, as seen by noting that the top and bottom rows of panels in figure 3 have widely different hybrid fitness values, yet the sizes of the different regions are changed only slightly. Thus, divergence in genes causing weak or strong reproductive isolation occurred, so long as gene flow was low. Intuitively, greater incompatibility did (slightly) increase mutation-order divergence.
(b) Strong inference tests of ecological and mutation-order speciation
A powerful way to test for mutation-order speciation would be to use experimental evolution in controlled laboratory experiments. Some experiments have found divergence or premating isolation to evolve between different replicate populations evolved in the same environment ([14,21]; for review ). However, most such experiments were designed to test ecological speciation, or did not manipulate the variables our model predicts are important for mutation-order divergence (e.g. Cohan & Hoffmann  evolved populations in similar environments without gene flow). Thus, experimental evolution studies that test the mutation-order process and manipulate the factors we have highlighted are warranted.
How might mutation-order speciation versus ecological speciation be supported in nature? Our predictions above generate strong inference (sensu ) tests of these two speciation processes. Namely, one could conduct reciprocal transplant experiments using populations that are known to exhibit some reproductive isolation (e.g. a DMI). Such experiments would rigorously establish (i) whether a given genotype has the same fitness across different environments (i.e. whether selection is nearly uniform or not), and (ii) the relative fitnesses of different genotypes in the same environment. If selection is uniform across environments, as observed in some experiments (Leimu & Fischer  for review), the hypothesis of ecological speciation is eliminated. On the other hand, evidence for divergent selection would falsify a basic assumption of mutation-order speciation.
If transplant experiments eliminated ecological speciation as a viable explanation, further information might add support for mutation-order divergence, by being consistent with the mutation-order process being at work in a specific set of populations. For example, mutation-order speciation would be more strongly supported if the populations under study have low rates of gene flow. Additionally, small fitness differences between genotypes would make mutation-order more likely. While estimating the selective advantage of specific alternative alleles might be difficult, specific mutations underlying adaptation (Hubbard et al.  for review) and speciation (Coyne & Orr  for review) have been detected, and at least one experiment has estimated the strength of selection acting on a specific mutation . A third factor of interest would be to determine the strength of hybrid incompatibility, with a high degree of incompatibility making mutation-order speciation a more likely explanation. Finally, coalescent-based techniques might allow tests of whether a period of allopatric divergence preceded gene flow (Hey  for review), although the efficacy of such techniques requires consideration .
As for the overall prevalence of mutation-order divergence as a driver of speciation, comparative studies across many sets of populations are required. In such studies, the predictions above could be examined by looking at whether the various factors correlate with the incidence of mutation-order divergence.
(c) Future work: extensions to the model
Two obvious extensions are to consider the effects of: (i) many alleles at many loci and (ii) stochastic changes via drift. With regard to the first point: how many loci and alleles tend to cause incompatibility? Empirical evidence is mixed. Strong hybrid incompatibility can be caused by just a few loci [1,42–44], but many different sets of a few loci can potentially generate incompatibility . Rapid evolution (i.e. fixation) of advantageous mutations is possible, once they emerge [1,42,43,46], arguing that numerous advantageous mutations will rarely be simultaneously segregating within a population. At the very least, our scenario of two alleles at two loci is a necessary biological precursor to any greater number of new alleles/loci. Nonetheless, we expect that interesting dynamics could emerge if more alleles/loci are considered, as might occur when gene flow maintains polymorphism. For example, an increased number of alleles could create more possible ways of achieving mutation-order speciation, but could also create more possible incompatibilities, which would inhibit the spread of newer, rare alleles. Analysis of this issue is beyond the scope of this article, but future work could modify past models of the accumulation of DMIs [10,11,27].
Another extension would be stochastic changes in allele frequencies. Our deterministic treatment is justified because we considered the rise and fixation of selectively advantageous mutations [2,32,33]. Adding stochasticity to the dynamics of beneficial allele replacements should not affect the conclusions qualitatively, especially with large population sizes [11,47]. This argument is bolstered by the fact that drift is expected to affect the results in the same direction as the documented trends. That is, when differences in the selective advantage of different mutations are large, the less selectively advantageous mutation should be lost by drift more readily than the more advantageous mutation, preventing mutation-order speciation. This is the exact outcome observed in our deterministic model when differences in selective advantages were great. When such differences are slighter, each mutation would be similarly affected by drift. Nonetheless, stochastic simulations could address how divergent selection, mutation-order and genetic drift might interact. Another avenue would be to examine how divergence of selected loci via the mutation-order process might promote a general barrier to gene flow at neutral loci (cf. [2,48]).
These represent interesting avenues for further research, although they should not alter our main conclusions: some mutation-order divergence can occur under relatively broad conditions, but fixation of alternative alleles via the mutation-order process speciation is extremely difficult unless gene flow is very low, and in such cases where gene flow is low, a number of other factors can promote or inhibit divergence.
We thank R. Safran and D. Schluter for discussion. S. Gavrilets, H. Rundle and two anonymous reviewers provided useful comments on earlier drafts of the manuscript.
- Received June 7, 2010.
- Accepted July 23, 2010.
- © 2010 The Royal Society