Abstract

Compared to the vast literature linking competitive interactions and speciation, attempts to understand the role of facilitation for evolutionary diversification remain scarce. Yet, community ecologists now recognize the importance of positive interactions within plant communities. Here, we examine how facilitation may interfere with the mechanisms of ecological speciation. We argue that facilitation is likely to (1) maintain gene flow among incipient species by enabling cooccurrence of adapted and maladapted forms in marginal habitats and (2) increase fitness of introgressed forms and limit reinforcement in secondary contact zones. Alternatively, we present how facilitation may favour colonization of marginal habitats and thus enhance local adaptation and ecological speciation. Therefore, facilitation may impede or pave the way for ecological speciation. Using a simple spatially and genetically explicit modelling framework, we illustrate and propose some first testable ideas about how, when, and where facilitation may act as a cohesive force for ecological speciation. These hypotheses and the modelling framework proposed should stimulate further empirical and theoretical research examining the role of both competitive and positive interactions in the formation of incipient species.

1. Introduction

Over the last 20 years, there has been increasing recognition that besides competition, facilitation, that is, the amelioration of biotic and/or abiotic conditions by neighboring organisms, is a key driver controlling assembly and dynamics of communities ([1], see [2] for recent reviews). The impact of facilitation in a community ecology context has been now widely accepted, but its evolutionary consequences have been rarely explored [2, 3]. There also has been an increasing awareness that positive and negative interactions usually cooccur in nature, and recent works have turned to examine conditions under which the outcome of biotic interactions is positive, negative, or neutral [46]. Yet, while there is a long-standing literature linking competitive interactions and speciation [710], the role of facilitation for evolutionary diversification has been largely overlooked. Very recently, a few empirical studies have attempted to bridge this gap [1113].

The theoretical study of positive interactions in an evolutionary context has an old history (e.g., [14]), but these models only consider mutualism, a two-species interaction that benefits both partners [15]. However, mutualism is characterized by balanced reciprocal benefits for species involved in a one-to-one relationship, while facilitation encompasses a larger spectrum of positive interactions between organisms [2, 16]. Namely, it commonly includes one-to-many interactions where one organism ameliorates habitat conditions for several beneficiary organisms and many-to-one interactions where an organism benefits from diffuse facilitation by neighbors. For example, facilitation between plants occurs when neighboring vegetation or a nurse plant (more generally called a benefactor) ameliorates environmental conditions directly, through the provision of additional resources or shelter from extreme conditions (e.g., high or low temperatures), or indirectly through protection from herbivores or via attraction of pollinators [16]. Since facilitation generates local variation in environmental conditions, it should also affect the manner in which organisms respond and adapt to their environment. However, it is important to realize that facilitation differs intrinsically from microhabitat amelioration due to abiotic components of the system (e.g., rock, ponds, or nutrient patches) because the source of environmental heterogeneity is a living organism. For instance, the distribution of favourable microhabitats may change at the time scale of ecological responses because of the population dynamics of benefactor species [17]. Additionally, the net effect of benefactors can shift from facilitation to competition due to interannual resource dynamics (e.g., change in rainfall pulse [18]).

The growing interest for eco-evolutionary dynamics, that is, the interplay between ecological and evolutionary processes [1922], demonstrates the need to study both the evolutionary consequences of ecological processes and the effect of evolutionary processes on ecological responses. In a recent study, Michalet et al. [17] provided one specific aspect of the problem from the benefactor’s point of view. Namely, they could show that ecotypic variation in benefactor traits may modify the manner in which they affect understory species. Viceversa, the understory also affected fitness of the benefactor. Here, we tackle this problem by looking at the other side of this coin and from a more general perspective, that is, the evolutionary consequences of facilitation on diversification of beneficiaries.

We aim at complementing the vast literature on competitive effects on evolutionary diversification by a conceptual model addressing the role of facilitation on distribution range, local adaptation, and genetic structure of beneficiary species. In particular, we discuss the likely effect of facilitation on microevolutionary processes and ultimately on ecological speciation. Ecological speciation can be defined as a particular type of speciation, where reproductive isolation evolves between populations as a result of ecologically based divergent selection [2325]. In this context, the interplay between biotic and abiotic factors is of paramount importance and, not surprisingly, the attempts to link ecological and evolutionary processes are the most tantalizing [3, 17, 26].

We do not present an exhaustive list of all potential evolutionary consequences of facilitation but instead propose an alternative to traditional views on the role of biotic interactions for evolutionary processes. Our overarching hypothesis is that, contrary to competition, facilitation may impede ecological speciation by maintaining gene flow between incipient species and preventing niche partitioning. Alternatively, we also discuss under which more specific circumstances facilitation could act as a stepping stone to promote ecological speciation. First, we highlight the role of facilitation in the context of niche theory. Secondly, we present how facilitation is likely to interfere with microevolutionary processes and propose a modeling framework to illustrate our arguments (Box  1). Lastly, we call for a new body of empirical and theoretical works to test the presented hypotheses (Box  2). The ideas and hypotheses developed in this paper are biased towards the plant kingdom (as is the facilitation literature in general, e.g., [2, 16]) but should apply equally to other kingdoms. For example, there is now some evidence showing the importance of facilitation for a wide range of organisms such as terrestrial animals [27], marine invertebrates [28], fungi [29], bacteria [3], and viruses [30].

2. Facilitation and Niche Theory, Moving from Species to Genotypes

Facilitation is likely to be as common as competition along environmental gradients [2]. It is generally hypothesized to increase in importance as environmental severity increases ([1], but see [31]) and is likely to enhance species diversity in harsh habitats [3134]. It is now recognized that facilitation can expand the distribution range of a species at its harsher end, or marginal habitat [35], whereas competition generally reduces it at the favorable end of the environmental gradient [4, 5, 36, 37]. By widening the distribution range of certain species not adapted to stressful conditions, facilitation allows them to persist with subordinate [4, 5, 33, 34, 36] or even dominant status in environments, where without facilitation they would not be able to persist [37]. Therefore, within a particular community, the ability of a species to be facilitated is strongly linked with its own niche [4, 5, 37] in that only species deviated from their niche optima are likely to be facilitated. More specifically, competitive species dominating favorable environments are more affected by facilitation as stress increases than stress-tolerant species adapted to harsher ecological conditions. Conversely, stress-tolerant species are competitively excluded as environmental stress decreases [31] due to a trade-off between competitive ability and stress-tolerance (see [5]).

Here, we apply these concepts developed on a species level to variation among incipient species. Namely, we propose to establish a missing link between facilitation as a force structuring communities and one affecting population structure. Due to the fact that different genotypes within a species may differ in their resource use and stress tolerance [38, 39], they can be characterized by contrasting ecological tolerances and different optima. For example, local adaptation along environmental gradients [40] will lead to ecotype-specific “realized niches.” Note that this concept is commonly incorporated in models of resource competition (e.g., [41, 42]). On the one hand, genotypes exhibiting traits linked to stress tolerance should be selected for, and thus be more common at the margin of the species’ realized niche (harsh conditions). On the other hand, genotypes enabling larger competitive response ability and low stress tolerance should dominate in populations in the center of the species niche, that is, more in favorable environments [43, 44].

Any ecological process that has the potential to significantly impact gene flow in marginal habitats may play a key role in ecological speciation [45, 46]. In that context, dispersal-related processes have been well investigated [43]. Here, we argue that the likelihood that core genotypes can persist in a species’ ecological margin without facilitation is particularly low. In the likely scenario of net dispersal from core (source) to marginal habitats (sinks) [43], facilitation may enhance the establishment of core genotypes in habitats where they otherwise would not be able to persist. This will have crucial consequences for adaptation to marginal habitats, and in turn for species diversification along environmental gradients, which we explain below.

3. Consequence of Facilitation for Gene Flow and Genetic Structure of Beneficiaries

A first mechanism by which facilitation may impact adaptive diversification is in enabling maladapted genotypes/forms to persist and reproduce outside of their niche optima. Figure 1 depicts two habitat types with contrasting environmental conditions (core and marginal), where different genotypes are distributed. Genotypes are supposed to differ by a key trait under selective pressure (e.g., a stress-tolerance-related trait). Local adaptation to marginal habitats is possible if gene flow from core habitats is limited due to pre- or postzygotic barriers ([47], but see [48]). For example, ecological filtering is a prezygotic barrier that will select genotypes best adapted to the local conditions and will reduce viability of the immigrants [49]. Postzygotic barriers may include reduced fitness of hybrids related to stress intolerance or to competitive displacement by more adapted forms [50, 51]. If these mechanisms are sufficiently strong to reduce gene flow between adapted and maladapted forms, trait divergence and ultimately reproductive isolation between populations may be observed (Figure 1(a)). An alternative scenario is depicted in Figure 1(b), where maladapted forms are facilitated by neighboring vegetation, that is, benefactors allow establishment, growth, and reproductive success of nonoptimal genotypes. In this case, facilitation may cause allele frequencies in the sink (marginal) to be similar to those of the source (core) habitat. The result would be stable coexistence between two genotypes, a core and a marginal genotype, and thus an overall increased genetic diversity and decreased inbreeding but also greater genetic swamping, that is, a decreased potential to adapt to the marginal conditions. Facilitation may thus be a cohesive force that tends to maintain gene flow between these populations and counteracts diversifying selection. To our knowledge, only one experimental field study has tested this hypothesis and supported our conjecture. This study was conducted on the annual grass Brachypodium distachyon in the Middle East. B. distachyon has a broad distribution range and occurs from Mediterranean (core habitat) to semiarid habitats (marginal habitat). Ecotypic differentiation has been observed in this species, with the semiarid ecotypes being more stress-tolerant [44, 52, 53] and the Mediterranean ecotype having a better ability to cope with competition [44]. When both ecotypes where transplanted to the semiarid environment, both were facilitated by the presence of shrubs but the Mediterranean ecotype could only survive to reproduction underneath the nurse [13]. Facilitation, therefore, enabled a maladapted morph to persist and reproduce in the marginal habitat as described above.

Another scenario would consider the secondary contact of two populations that exhibit trait differentiation. Reinforcement, that is, selection against hybridization, and introgression are two opposing forces that operate during this secondary contact phase, the former enhancing genetic divergence and the latter having the potential to erase it [23]. Secondary contact zones often represent marginal habitats for incipient species (e.g., [54]). Facilitation could favor the cooccurrence of ecotypes according to the mechanism described above and ultimately lead to increased gene flow between ecotypes. Facilitation may also circumvent postzygotic barriers by increasing the fitness of introgressed forms. By increasing the chance of introgression and by preventing reinforcement, we suggest that facilitation could strongly impede further genetic differentiation between ecotypes.

Persistence of core genotypes in marginal habitats increases genetic diversity and thus increases the potential for adapting to changes [46], at least if these occur towards conditions more similar to the core habitat. Higher genetic diversity in marginal habitats compared to core habitats, that is, arid versus Mediterranean environments [55], could support this prediction if this diversity is maintained through time. The trend toward higher genetic variability may be enhanced when the marginal environment varies in time. In that case, the different genotypes inhabiting patches with and without facilitation could have an advantage in different years. Facilitation would then increase the role of temporal variation in supporting stable coexistence of genotypes (e.g., [56]), similar to the storage effect [57] that increases diversity at a species level.

Recent numerical models have been explored under which conditions ecological speciation is likely to occur [51, 5861]. The hypothesis that facilitation acts as a cohesive force and increases genetic diversity in marginal habitats can be tested with a spatially and genetically explicit modeling framework derived from Kirkpatrick and Barton [58] and Bridle et al. [62] (Box  1). Our simulations confirm that this modeling framework represents a flexible tool to explore the important set of parameters and scenarios that we believe should be investigated (Box  2). We chose here to explore a reduced set of parameters, that is, size of the patch created by the benefactors and the environmental difference between the core and the marginal habitat (Box  1). Our simulations show that, by producing a mosaic of mild conditions in a harsh environment, benefactors increase gene flow between populations from core and marginal habitats [63], that is, lower FST between populations for both neutral [64] and selected loci (Figures 2(a) and 2(b), resp.). They thus ultimately prevent local adaptation in the harsh environment, as expressed in lower fitness in the marginal environment with facilitation compared to a scenario without facilitation (see Figure 2(c) comparing without benefactors versus smaller patch size benefactors).

4. Benefactors and Disruptive Selection

Though facilitation is likely to be a cohesive force, we would also like to touch upon possible conditions under which positive interactions may contribute to, rather than impede, ecological speciation. This is most easily illustrated with particularly stressful environments, where benefactors such as cushion plants in alpine ecosystems or shrubs in water-limited ecosystems [16], represent islands of milder conditions or islands of fertility. At the community scale, this results in a mosaic of microhabitats determined by the distribution and abundance of benefactors. Such systems have provided the most conclusive evidence for facilitation [65], even if the outcome of biotic interactions may strongly depend upon the particular location of the focal species within this mosaic. Arguably, the patchy distribution of biotically engineered favorable microhabitats within a harsher matrix (e.g., shrubs creating islands of fertility in deserts, see Figure 3) may represent a case at hand for disruptive selection between adjacent subpopulations (see [66] for the example of Bromus erectus adapting to different Thymus vulgaris chemotypes in calcareous grasslands). The benefactors could promote local adaptation in the understory habitats as observed in our simulation when they were forming large patches (Figure 2). Therefore, facilitation could also pave the way to ecological speciation, particularly if directional selection leads to differences in traits that affect the reproductive system, such as flowering time [67]. While this scenario is certainly a possibility, we suggest that it is less likely than the opposite case due to the small spatial scale of environmental heterogeneity usually created by benefactors. Namely, the size of these islands of mild conditions roughly corresponds to the size of the nurse’s crown or to a clump of shrubs and is usually small compared to the dispersal kernel and pollination distance [68]. Furthermore, in extreme habitats the nurses also function as a sink for seeds [69]. As a result, it is unlikely that gene flow will be severely restricted between populations of the understory habitats and those of the harsher matrix. Nonetheless, further empirical and theoretical studies (Box  2) are needed to test for the overwhelming importance of facilitation as a cohesive, rather than a disruptive, force under these particular environmental conditions.

In the context of colonization of a new habitat, where individuals from the population of the core habitat could not persist away from a benefactor in the first place (Figure 4), benefactors could first promote local adaptation at the niche margin and then in the harsh matrix. Recent studies have pointed out the role of intermediate quality habitats in ecotypic differentiation and in the colonisation of a new niche [70]. In the Littorina saxatilis case study, it may be suggested that facilitation played a role in this process given the patchy distribution of mussels and barnacles (in intermediate habitats) associated with the two ecotypes observed in Spain for these species (SU and RB, resp.). However, this hypothesis needs to be properly tested. Thus, on the one hand, positive interactions in intermediate habitats may favor niche broadening and local adaptation. On the other hand, facilitation may contribute to gene flow maintenance between the core and the newly colonized habitats as described in Figure 2. To our opinion, this model system provides excellent opportunities to examine how the spatial and temporal dynamics of positive interactions could interfere with ecological speciation.

5. Facilitation and Niche Conservatism

Linkage between facilitation and macroevolutionary processes has recently arisen from experimental work on semiarid Mexican plant communities [11, 71]. Here, benefactor species were phylogenetically unrelated to beneficiary species; the former belonged to recent Quaternary lineages and were well adapted to the present dry climate, the latter showed closer affinities with Tertiary lineages and were drought intolerant [11]. It was hypothesized that less adapted “Tertiary” species survive because of facilitation [11]. In this scenario, our model of facilitation acting as a cohesive force is likely to explain niche conservatism of Tertiary species and the fact that they did not become more stress-tolerant or go extinct through evolutionary time. Transposed to the microevolutionary scale, these results suggest that facilitation may cause a genotype’s niche to remain unchanged over time despite environmental changes (Figure 1(c)). Namely, one may imagine that core and marginal habitats are temporally segregated, that is, a population experiences increasingly stressful conditions with time and adaptation to the new conditions is hampered by the persistence of facilitated and maladapted genotypes from more favorable times. Recently, there has been a renewed interest in evaluating niche conservatism across lineages [72], in particular by using species distribution models [73]. To our knowledge, none of these models has properly included the balance between positive and negative interactions along gradients to address the underlying mechanisms and the resulting patterns of niche conservatism.

6. Conclusion

We have argued that facilitation can be viewed as a cohesive force limiting genetic and phenotypic differentiation in marginal habitats. Facilitation may also explain the lack of adaptive diversification over time in certain lineages. On the other hand, under very specific conditions, facilitation could also operate as a stepping stone and pave the way to ecological speciation. This calls for a new body of empirical or theoretical studies addressing the strength of these two possible effects of facilitations on ecological speciation (Box  2). Our focus was on the population scale to provide the linkage with mechanisms of speciation. We thus complement recent opinions about the role of various types of facilitation on macro-evolutionary processes [74] and about the effects of evolutionary processes on facilitative interactions [17]. So far, there is no empirical study that has simultaneously examined genetic population structure of beneficiaries and facilitation for contrasting genotypes/phenotypes in natural conditions. Therefore, there is an urgent need to design novel studies that more closely associate community ecologists and evolutionary biologists to tackle these questions (Box  2).

To conclude, we suggest that plant ecologists should invest more into studying biotic interactions at the population level. First, shifting from species to populations and from population to genotypes in the design of ecological experiments is essential to examine the role of facilitation in ecological speciation [12, 13]. Second, population models should gain more realism by acknowledging the balance between positive and negative interactions along environmental gradients and their underlying ecological mechanisms. This research agenda should help to bridge the gap between community ecology and evolutionary biology. Bruno et al. [36] claimed that “including facilitation in niche theory will challenge some of our most cherished paradigms.” Ecological speciation is another paradigm that is worth revisiting by examining the impact of both negative and positive interactions on the evolutionary history of incipient species. Facilitation is certainly challenging our understanding of the origin of new species.

Boxes
Box  1: Description of the Modeling Framework
The model of Bridle et al. [62] (initially developed to study adaptation on an environmental gradient) was modified as follows to explore the effect of facilitation on ecological speciation.
The Environment
We used a grid of 8192 × 4096 cells wrapped on a torus. Without facilitation, the left half is the core habitat with an optimum ( in Bridle et al. [62]) of and the right half is the marginal habitat with optimum of and reduced carrying capacity (, the carrying capacity if defined on a circle with radius of 50 cells). In the presence of benefactors, the half right-hand side of the grid is fragmented in a checkerboard-like pattern with alternating benefactor and marginal regions. The environments produced by the presence of benefactors have an optimum intermediate between the core and the marginal habitat () and the same carrying capacity () as the core environment.
The Individuals
Individuals are diploid with 128 bi-allelic unlinked loci (64 under selection and 64 neutral). Allele values are either zero or one with a symmetric mutation probability between them of 10−4. The phenotype () is the sum of those loci and thus range from zero to 128. The fitness of males and females is determined by a logistic growth and local adaptation where if the fitness, is the maximum rate of increase (set to 0.8), is the local density in a radius of 50 cells, s is the strength of stabilizing selection (set at 0.125, for the two scenarios presented in Figures 2 and 4), and and are the carrying capacity and local optimum at the position of the individual. The markers start with no polymorphism (all alleles are zero) and have a mutation rate of 10−3 that change their value by plus or minus one, so each allele is represented by an integer. We used FST [75] to infer gene flow among populations using neutral markers and selected loci separately.
Reproduction
To find a mate, females scan in a radius of mating distance (MD) set at 150 cells and choose randomly a male with probability proportional to its fitness (). Each couple leaves an average of offsprings, drawn from a Poisson distribution. Each offspring is placed at a distance from the mother drawn from a normal distribution of mean zero and standard deviation D (set at 100 cells) in a direction randomly selected from a uniform distribution.
Initial Conditions
We started each simulation with 500 individuals randomly located in the centre 700 × 700 of the core environment, and we follow their evolution over 2000 generations. Three patch sizes were used to simulate the presence of benefactors in the harsh environment, smaller patches (512 cells, isolated benefactors in the landscape), intermediate sizes (1024 cells), and large patch sizes (2048 cells). Each parameter set was replicated ten times (the boxplots in Figures 2 and 4 reflect the ten replicates).
Example of Two Contrasting Scenarios
The two sets of simulations that we performed, corresponding to two different scenarios, are presented in Figures 2 and 4. The first scenario presents the case where facilitation inhibits ecological speciation (, , , Figure 2), and the second scenario presents the case where facilitation promotes ecological speciation (, , , Figure 4). Following (1), for the first scenario at carrying capacity in both environments, , the fitness for an individual locally adapted to the core environment () is equal to 2, but the same individual would have a fitness equal to 1 in the marginal environment, representing the equivalent of a selection coefficient of 50%. This value is consistent with the average measure of local adaptation recently reported for reciprocal transplantation experiments [76]. In the second scenario, we set the environmental difference between the core () and the marginal environment () such that individuals locally adapted to the core environment () could not maintain a viable population in the marginal environment.
Box  2: Facilitation and Ecological Speciation: Ways Forward
Stress Gradients as a Model System
For empirically studying the role of facilitation in ecological speciation, one should focus on target species with populations widely distributed along stress gradients, and for which facilitation and trait divergence in marginal habitats have been documented. Common examples are gradients of water availability in semiarid landscapes or gradients of temperature along altitudinal gradients. Population genetic structure and inferences about gene flow could be evaluated first descriptively by using neutral or selected markers and by QTL analyses [77]. The next step would be to investigate adaptive genetic variation between microhabitats with and without benefactor. These should be accompanied by experiments testing whether genotypes from “favorable” habitats rely on facilitation to persist and reproduce in marginal habitats [13]. These questions can be addressed with the same experimental designs that are carried out at the species level, that is, reciprocal transplants or sowing experiments coupled with neighbor manipulations.
Which Traits to Measure?
A particular challenge to exploring the effects of facilitation is to identify adaptive genetic variation in traits (1) that are under disruptive selection pressure along the investigated gradient and (2) for which trait expression strongly depends on the balance between stress tolerance and success in the core habitat. Seed size, for example, is a good candidate because it is related to seedling size and tolerance of hazards for young seedlings [78]. As competitive ability and stress tolerances ability are related to the probability of being facilitated [5], traits linked with the resource capture and utilization, such as leaf traits, plant stature traits, and root traits, need to be considered as well (see [79]).
A key trait would be both under strong selective pressure and closely related to mate choice, the so-called “magic traits” [45, 51, 80]. Because the benefactor may modify the length of the favorable season, for instance, phenological traits such as time to flowering or length of the life cycle [44, 81] may fall in this category and are worth examining. Depending on how benefactors affect these key traits, facilitation could potentially either hinder, accelerate or even lead to ecological speciation (see “Section 4”).
Theory and Simulations
There are large open avenues for tackling these questions with spatially explicit models of population dynamics, building on a rich set of various modeling frameworks [26, 51, 62, 70, 8284]. We present here an example (Figures 2 and 4) but such models should be refined to explore the specific feature of facilitation and its implication on adaptive diversification along gradients and other spatiotemporal configurations. Namely, the crucial role of dispersal for adaptation to marginal habitats should be investigated in the presence of facilitation. For example, an important consequence of facilitation may be that it enhances “effective” dispersal from the source to sinks, because establishment probabilities of maladapted genotypes are enhanced. Furthermore, the sensitivity of model outputs to parameters controlling (i) the shape of the competition versus facilitation balance along ecological gradients, that is still widely debated [6, 31], (ii) the spatial distribution and dynamics of nurse plants in marginal habitats, (iii) the complex interactions among benefactors and beneficiaries [17], and (iv) the dispersal ability and life history traits of focal species, has to be investigated. Ultimately, these models should help us sharpen and phrase testable hypotheses about the relative importance of facilitation on adaptation and ecological speciation. Finally, a spatially and genetically explicit model [51, 5961] could couple both benefactor effect on beneficiaries’ microevolutionary responses (this study) and the feedback effect of beneficiaries on population and evolutionary dynamics of benefactor species (cost of facilitation, [17]). Such an approach would shed light on how the evolutionary consequence of facilitation affects community assembly and dynamics.

Authors’ Contribution

P. Liancourt wrote the first draft of the manuscript, X. Thibert-Plante performed the modelling work. All authors made a substantial intellectual contribution to the paper and its revision. P. Choler, N. Gross and X. Thibert-Plante had equal contribution to the paper.

Acknowledgment

P. Liancourt and K. Tielbörger were supported by the GLOWA Jordan River Project funded by the German Ministry of Science and Education (BMBF), and further support was granted to K. Tielbörger by the German Research Foundation (TI338/10-1). P. Choler was partly funded by a Marie-Curie fellowship (contract CASOAR MOIF-CT-2006-039688). X. Thibert-Plante conducted this work while a Postdoctoral Fellow at the National Institute for Mathematical and Biological Synthesis, an institute sponsored by the National Science Foundation, the U.S. Department of Homeland Security, and the U.S. Department of Agriculture through NSF Award no. EF-0832858, with additional support from The University of Tennessee, Knoxville, and Le Fonds Québécois de la Recherche sur la Nature et les Technologies (FQRNT) directly to X. Thibert-Plante. This research was also partially supported by NIH Grant GM-56693 to S. Gavrilets. The authors thank R. Glade for editing the paper to improve the English, J. Bridle and J. Polechová for sharing with them details of their model, S. Gavrilets, A. P. Henry, R. M. Callaway, and R. Michalet for feedback on the paper, and R. Faria and the anonymous reviewers for helpful comments on the paper and suggesting the analysis of the fixation index on the selected loci.