Interspecific Competition, Environmental Gradients, Gene Flow, and the Coevolution of Species' Borders

Darwin viewed species range limits as chiefly determined by an interplay between the abiotic environment and interspecific interactions. Haldane argued that species' ranges could be set intraspecifically when gene flow from a species' populous center overwhelms local adaptation at the periphery. Recently, Kirkpatrick and Barton have modeled Haldane's process with a quantitative genetic model that combines density-dependent local population growth with dispersal and gene flow across a linear environmental gradient in optimum phenotype. To address Darwin's ideas, we have extended the Kirkpatrick and Barton model to include interspecific competition and the frequency-dependent selection that it generates, as well as stabilizing selection on a quantitative character. Our model includes local population growth, movements over space, natural selection, and gene flow. It simultaneously addresses the evolution of character displacement and species borders. It reproduces the Kirkpatrick and Barton single-species result that limited ranges can be produced with sufficiently steep environmental gradients and strong dispersal. Further, in the absence of environmental gradients or barriers to dispersal, interspecific competition will not limit species ranges at evolutionary equilibrium. However, interspecific competition can interact with environmental gradients and gene flow to generate limited ranges with much less extreme gradient and dispersal parameters than in the single-species case. Species display character displacement in sympatry, yet the reduction in competition that results from this displacement does not necessarily allow the two species to become sympatric everywhere. When species meet, competition reduces population densities in the region of overlap, which, in turn, intensifies the asymmetry in gene flow from center to margin. This reduces the ability of each species to adapt to local physical conditions at their range limits. If environmental gradients are monotonic but not linear, the transition zone between species at coevolutionary equilibrium occurs where the environmental gradient is steepest. If productivity gradients are also introduced into the model, then patterns similar to Rapoport's rule emerge. Interacting species respond to climate change, as it affects the optimal phenotype over space, by a combination of range shifts and local evolution in mean phenotype, while solitary species respond solely by range shifts. Finally, we compare empirical estimates for intrinsic growth rates and diffusion coefficients for several species to those needed by the single-species model to produce a stable limited range. These empirical values are generally insufficient to produce limited ranges in the model suggesting a role for interspecific interactions.     

On evolution under symmetric and asymmetric competitions

This paper considers the coevolution of phenotypic traits in a community comprising two competitive species subject to strong Allee effects. Firstly, we investigate the ecological and evolutionary conditions that allow for continuously stable strategy under symmetric competition. Secondly, we find that evolutionary suicide is impossible when the two species undergo symmetric competition, however, evolutionary suicide can occur in an asymmetric competition model with strong Allee effects. Thirdly, it is found that evolutionary bistability is a likely outcome of the process under both symmetric and asymmetric competitions, which depends on the properties of symmetric and asymmetric competitions. Fourthly, under asymmetric competition, we find that evolutionary cycle is a likely outcome of the process, which depends on the properties of both intraspecific and interspecific competition. When interspecific and intraspecific asymmetries vary continuously, we also find that the evolutionary dynamics may admit a stable equilibrium and two limit cycles or two stable equilibria separated by an unstable limit cycle or a stable equilibrium and a stable limit cycle.     

A quantitative genetic competition model for sympatric speciation

I use multilocus genetics to describe assortative mating in a competition model. The intensity of competition between individuals is influenced by a quantitative character whose value is determined additively by alleles from many loci. With assortative mating based on this character, frequency- and density-dependent competition can subdivide a population with an initially unimodal character distribution. The character distribution becomes bimodal, and the subpopulations corresponding to the two modes are reproductively separated because mating is assortative. This happens if the resource distribution is unimodal, i.e. even if selection due to phenotypic carrying capacities is not disruptive. The results suggest that sympatric speciation due to frequency-dependent selection can occur in quite general ecological scenarios if mating is assortative. I also discuss the evolution of assortative mating. Since it induces bimodal phenotype distributions, assortative mating leads to a better match of the resources if their distribution is also bimodal. Moreover, in a population with a bimodal phenotype distribution, the average strength of frequency-dependent competition is lower than in a unimodal population. Therefore, assortative mating permits higher equilibrium densities than random mating even if the resource distribution is unimodal. Thus, even though it may lead to a less efficient resource use, assortative mating is favoured over random mating because it reduces frequency-dependent effects of competition.     

SPERM COMPETITION PROMOTES ASYMMETRIES IN REPRODUCTIVE BARRIERS BETWEEN CLOSELY RELATED SPECIES

Reproductive barriers between closely related species are often incomplete and asymmetric, but the evolutionary significance of these well-known phenomena remains unsolved. We test the hypothesis that the degree of gametic incompatibility in reciprocal crosses is associated to levels of sperm competition because this selective force favors both increased sperm competitiveness and ovum defensiveness. Using three species of Mus with high, intermediate, and low levels of sperm competition, we examined fertilization rates in competitive and noncompetitive contexts. We found that the influence of sperm competition upon sperm competitiveness is as strong as it is upon ovum defensiveness, revealing an effect upon female gametes so far overlooked. As a result, fertilization success was strongly related to differences in sperm competition levels between species providing sperm and ova, thus generating major asymmetries in reciprocal crosses. When placed in competition, conspecific sperm maintained levels of fertilization success similar to those found in noncompetitive contexts, at the expense of the success of heterospecific sperm. When only heterospecific sperm competed, species with highest levels of sperm competition outcompeted others and asymmetries were exacerbated. We conclude that sperm competition explains both the degree of gametic isolation and the degree of asymmetries between closely related species.     

Parallel evolution and ecological selection: replicated character displacement in spadefoot toads

Ecological character displacement—trait evolution stemming from selection to lessen resource competition between species—is most often inferred from a pattern in which species differ in resource-use traits in sympatry but not in allopatry, and in which sympatric populations within each species differ from conspecific allopatric populations. Yet, without information on population history, the presence of a divergent phenotype in multiple sympatric populations does not necessarily imply that there has been repeated evolution of character displacement. Instead, such a pattern may arise if there has been character displacement in a single ancestral population, followed by gene flow carrying the divergent phenotype into multiple, derived, sympatric populations. Here, we evaluate the likelihood of such historical events versus ongoing ecological selection in generating divergence in trophic morphology between multiple populations of spadefoot toad (Spea multiplicata) tadpoles that are in sympatry with a heterospecific and those that are in allopatry. We present both phylogenetic and population genetic evidence indicating that the same divergent trait, which minimizes resource competition with the heterospecific, has arisen independently in multiple sympatric populations. These data, therefore, provide strong indirect support for competition''s role in divergent trait evolution.     

Ecological character displacement caused by reproductive interference

We carried out a theoretical investigation of whether ecological character displacement can be caused by reproductive interference. Our model assumes that a quantitative character is associated with both resource use and species recognition, and that heterospecific mating incurs costs. The model shows that ecological character displacement can occur as a consequence of evolution of premating isolation; this conclusion is based on the premise that resource competition is less intense between species than within species and that the ecological character also contributes to premating isolation. When resource competition between species is intense, extinction of either species may occur by competitive exclusion before ecological character divergence. Some observational studies have shown that character displacement in body size is associated with not only resources use but also species recognition. We propose that body size displacement can occur as a consequence of evolution of premating isolation. Our results suggest that ecological character displacement results from reproductive character displacement.     

Including competitive asymmetry in measures of local interference in plant populations

Summary Although considerable evidence exists that plant competition is generally asymmetric or one-sided, with larger plants having a disproportionate competitive effect on smaller plants, currently employed measures of local interference generally assume that competition is two-sided. We describe a simple measure of competitive asymmetry in which the effects of neighbors smaller than a focal individual are discounted by a constant factor, and include this variable in a composite measure of local interference. In this model competition varies between complete asymmetry (the effects of smaller plants are entirely discounted) and complete symmetry (the competitive effect of a neighbor is proportional to its size). The proposed method is applied to two natural populations and one experimental monoculture. In all cases an asymmetric model provides the best fit to the data. Completely two-sided models account for 26–39% of the variance in relative growth rate, while relatively one-sided models account for 44–57%. The increases in r ² values resulting from the inclusion of asymmetry are significant in the two cases in which the data permit randomization tests. Our results suggest that interference is completely asymmetric in a population of Impatiens pallida, a species with very low root allocation and a shallow crown, and somewhat less asymmetric in an experimental monoculture of Ambrosia artemisiifolia and a natural stand of Pinus rigida, cases in which competition for water and nutrient resources is likely to be of greater importance.     

Networks and dominance hierarchies: does interspecific aggression explain flower partitioning among stingless bees?

1. The distribution of consumers among resources (trophic interaction network) may be shaped by asymmetric competition. Dominance hierarchy models predict that asymmetric interference competition leads to a domination of high quality resources by hierarchically superior species. 2. In order to determine the competitive dominance hierarchy and its effect on flower partitioning in a local stingless bee community in Borneo, interspecific aggressions were tested among eight species in arena experiments. 3. All species tested were strongly mutually aggressive in the arena, and the observed interactions were often lethal for one or both opponents. Aggression significantly increased with body size differences between fighting pairs and was asymmetric: larger aggressors were superior over smaller species. Additional aggression tests involved dummies with surface extracts, and results suggest that species‐ and colony‐specific surface profiles are important in triggering the aggressive behaviour. 4. Sixteen stingless bee species were observed foraging on 41 species of flowering plants. The resulting bee–flower interaction network showed a high degree of generalisation (network‐level specialisation H₂’ = 0.11), corresponding to a random, opportunistic distribution of bee species among available flower species. 5. Aggressions on flowers were rare and only occurred at a low level. The dominance hierarchy obtained in the arena experiments did not correlate significantly with plant quality, estimated as the number of flowers per plant or as total bee visitation rate. 6. Our findings suggest that asymmetries in interference competition do not necessarily translate into actual resource partitioning in the context of complex interacting communities.     

Asymmetric interspecific competition drives shifts in signalling traits in fan-throated lizards

Interspecific competition can occur when species are unable to distinguish between conspecific and heterospecific mates or competitors when they occur in sympatry. Selection in response to interspecific competition can lead to shifts in signalling traits—a process called agonistic character displacement. In two fan-throated lizard species—Sitana laticeps and Sarada darwini—females are morphologically indistinguishable and male agonistic signalling behaviour is similar. Consequently, in areas where these species overlap, males engage in interspecific aggressive interactions. To test whether interspecific male aggression between Si. laticeps and Sa. darwini results in agonistic character displacement, we quantified species recognition and signalling behaviour using staged encounter assays with both conspecifics and heterospecifics across sympatric and allopatric populations of both species. We found an asymmetric pattern, wherein males of Si. laticeps but not Sa. darwini showed differences in competitor recognition and agonistic signalling traits (morphology and behaviour) in sympatry compared with allopatry. This asymmetric shift in traits is probably due to differences in competitive abilities between species and can minimize competitive interactions in zones of sympatry. Overall, our results support agonistic character displacement, and highlight the role of asymmetric interspecific competition in driving shifts in social signals.     

Spite and the scale of competition

In recent years there has been a large body of theoretical work examining how local competition can reduce and even remove selection for altruism between relatives. However, it is less well appreciated that local competition favours selection for spite, the relatively neglected ugly sister of altruism. Here, we use extensions of social evolution theory that were formulated to deal with the consequences for altruism of competition between social partners, to illustrate several points on the evolution of spite. Specifically, we show that: (i) the conditions for the evolution of spite are less restrictive than previously assumed; (ii) previous models which have demonstrated selection for spite often implicitly assumed local competition; (iii) the scale of competition must be allowed for when distinguishing different forms of spite (Hamiltonian vs. Wilsonian); (iv) local competition can enhance the spread of spiteful greenbeards; and (v) the theory makes testable predictions for how the extent of spite should vary dependent upon population structure and average relatedness.     

Does competition drive character differences between species on a macroevolutionary scale?

Sympatric sister species generally have a degree of phenotypic differentiation that allows them to coexist. It has been well documented that phenotypic similarity results, through resource competition, in one of two major outcomes: local extinction of either competitor or character displacement. Limiting similarity suggests that there is a maximum degree of phenotypic niche overlap with which similar species may coexist. Breaching that maximum would result in exclusion. Character displacement, on the other hand, implies that the species differentiate phenotypically so that resource competition is reduced to the point where coexistence is possible. While it has been suggested that these theories have the potential to accelerate (character displacement) or limit phenotypic evolution (competitive exclusion) on microevolutionary time scales, their effects on macroevolution remain under‐studied. If competition accelerates evolution on a macroevolutionary scale, one would expect that phenotypic diversity increases as novel species ‘push aside’ existing species. On the other hand, one might also expect that phenotypic evolution comes to a halt as novel species are trapped in the (ever decreasing) phenotypic space not yet occupied by existing species, except at the extremes of the phenotypic spectrum. Studying the current geographical ranges of more than 3000 extant species representing 29 mammalian families and their respective body masses, I found little evidence of competition accelerating body size differentiation between species.     

Evolution of Dispersal in a Structured Metapopulation Model in Discrete Time

In this article, a structured metapopulation model in discrete time with catastrophes and density-dependent local growth is introduced. The fitness of a rare mutant in an environment set by the resident is defined, and an efficient method to calculate fitness is presented. With this fitness measure evolutionary analysis of this model becomes feasible. This article concentrates on the evolution of dispersal. The effect of catastrophes, dispersal cost, and local dynamics on the evolution of dispersal is investigated. It is proved that without catastrophes, if all population–dynamical attractors are fixed points, there will be selection for no dispersal. A new mechanism for evolutionary branching is also found: Even though local population sizes approach fixed points, catastrophes can cause enough temporal variability, so that evolutionary branching becomes possible.     

Year-class coexistence in biennial plants

I extend the well known and biologically well motivated Skellam model of plant population dynamics to biennial plants. The model has two attractors: either one year class competitively excludes the other, resulting in 2-cycles with only vegetative vs only flowering plants in alternating years, or the two year classes coexist at an interior equilibrium. Contrary to earlier models, these two attractors can exist also simultaneously. I investigate the robustness of the model by including delayed flowering, a common phenomenon in plants, and provide a full numerical bifurcation analysis of the generalized model. High fecundity implies strong competition within year classes and promotes coexistence, whereas high survival results in strong competition between year classes and promotes competitive exclusion. Delayed flowering tends to stabilize the interior equilibrium, but (unlike in density-independent matrix models) the population cycles are robust with respect to some delay in flowering.     

Monarch butterfly spatially discrete advection model

We study the population cycles of the Monarch butterfly using one of the simplest systems incorporating both migration and local dynamics. The annual migration of the Monarch involves four generations. Members of Generations 1-3 (occasionally 4) migrate from the over-wintering site in Central Mexico to breeding grounds that extend as far north as the Northern United States and Southern Canada. A portion of the Generation 3 and all members of the Generation 4 butterflies begin their return to the over-wintering grounds in August through October where they enter reproductive diapause for several months. We developed a simple discrete-time island chain model in which different fecundity functions are used to model the reproductive strategies of each generation. The fecundity functions are selected from broad classes of functions that capture the effects of either contest or scramble intraspecific competition in the Monarch population. The objectives of our research are multiple and include the study of the generationally dependent intraspecific competition and its effect on the pool size of migrants as well as the persistence of the overall butterfly populations. The stage structure used in modeling the Monarch butterfly dynamics and their generationally dependent reproductive strategies naturally support fluctuating patterns and multiple attractors. The implications of these fluctuations and attractors on the long-term survival of the Monarch butterfly population are explored.     

Multiple attractors in a discrete competition model

The existence of multiple attractors in a competition model implies that the question of coexistence vs. extinction can depend on initial conditions. A discrete stage-structured model of two competing species is derived from a well-tested single-species model of insect populations, and is shown to exhibit multiple attractors for parameter values similar to those used in laboratory experiments which demonstrated chaos in population dynamics. The corresponding basins of attraction are investigated and shown to have very complex structures, and the initial stage structure of the populations is shown to have a significant impact on final outcomes.     

Two-patch dispersal-linked compensatory-overcompensatory spatially discrete population models

We study the role of asynchronous and synchronous dispersals on discrete-time two-patch dispersal-linked population models, where the pre-dispersal local patch dynamics are of mixed compensatory and overcompensatory types. Single-species dispersal-linked models behave as single-species single-patch models whenever all pre-dispersal local patch dynamics are compensatory and dispersal is synchronous. However, the dynamics of the corresponding two-patch population model connected by asynchronous dispersal depends on the dispersal rates. The species goes extinct on at least one patch when the asynchronous dispersal rates are high, while it persists when the rates are low. We use numerical simulations to show that in both synchronous and asynchronous mixed compensatory and overcompensatory systems, symmetric and asymmetric dispersals can control and impede the onset of cyclic population oscillations via period-doubling reversal bifurcations. Also, we show that in mixed systems both asynchronous and synchronous dispersals are capable of altering the pre-dispersal local patch dynamics from overcompensatory to compensatory dynamics. Dispersal-linked population models with ‘unstructured’ overcompensatory pre-dispersal local dynamics connected by synchronous dispersal can generate multiple attractors with fractal basin boundaries. However, mixed compensatory and overcompensatory systems appear to exhibit single attractors and not coexisting (multiple) attractors.     

Metapopulation dynamics with quasi-local competition

Stepping-stone models for the ecological dynamics of metapopulations are often used to address general questions about the effects of spatial structure on the nature and complexity of population fluctuations. Such models describe an ensemble of local and spatially isolated habitat patches that are connected through dispersal. Reproduction and hence the dynamics in a given local population depend on the density of that local population, and a fraction of every local population disperses to neighboring patches. In such models, interesting dynamic phenomena, e.g. the persistence of locally unstable predator-prey interactions, are only observed if the local dynamics in an isolated patch exhibit non-equilibrium behavior. Therefore, the scope of these models is limited. Here we extend these models by making the biologically plausible assumption that reproductive success in a given local habitat not only depends on the density of the local population living in that habitat, but also on the densities of neighboring local populations. This would occur if competition for resources occurs between neighboring populations, e.g. due to foraging in neighboring habitats. With this assumption of quasi-local competition the dynamics of the model change completely. The main difference is that even if the dynamics of the local populations have a stable equilibrium in isolation, the spatially uniform equilibrium in which all local populations are at their carrying capacity becomes unstable if the strength of quasi-local competition reaches a critical level, which can be calculated analytically. In this case the metapopulation reaches a new stable state, which is, however, not spatially uniform anymore and instead results in an irregular spatial pattern of local population abundance. For large metapopulations, a huge number of different, spatially non-uniform equilibrium states coexist as attractors of the metapopulation dynamics, so that the final state of the system depends critically on the initial conditions. The existence of a large number of attractors has important consequences when environmental noise is introduced into the model. Then the metapopulation performs a random walk in the space of all attractors. This leads to large and complicated population fluctuations whose power spectrum obeys a red-shifted power law. Our theory reiterates the potential importance of spatial structure for ecological processes and proposes new mechanisms for the emergence of non-uniform spatial patterns of abundance and for the persistence of complicated temporal population fluctuations.     

UNEXPLAINED SPLIT SEX RATIOS IN THE NEOTROPICAL PLANT-ANT, ALLOMERUS OCTOARTICULATUS VAR. DEMERARAE (MYRMICINAE): A TEST OF HYPOTHESES

We investigated sex allocation in the Neotropical ant Allomerus octoarticulatus var. demerarae . Because Allomerus is a plant symbiont, we could make geographically extensive collections of complete colonies and of foundresses in saplings, allowing us to estimate not only population- and colony-level sex allocation but also colony resource levels and the relatednesses of competing ant foundresses. This species exhibits a strongly split sex ratio, with 80% of mature colonies producing ≥90% of one sex or the other. Our genetic analyses (DNA microsatellites) reveal that Allomerus has a breeding system characterized by almost complete monogyny and a low frequency of polyandry. Contrary to theoretical explanations, we find no difference in worker relatedness asymmetries between female- and male-specialist colonies. Furthermore, no clear link was found between colony sex allocation and life history traits such as the number of mates per queen, or colony size, resource level, or fecundity. We also failed to find significant support for male production by workers, infection by Wolbachia , local resource competition, or local mate competition. We are left with the possibility that Allomerus exhibits split sex ratios because of the evolution of alternative biasing strategies in queens or workers, as recently proposed in the literature.     

Caught in the middle: Asymmetric competition causes high variance in intermediate trait abundances

In asymmetric competition between two individuals of the same or different species, one individual has a distinct advantage over the other due to a particular beneficial trait. An important trait that induces asymmetric competition is size (body size in animals, height in plants). There is usually a trade-off between fecundity and the trait that leads to competitive superiority (e.g. seed number vs seed size), enabling coexistence of populations with different trait values. These predictions on coexistence are based on classic deterministic models. Here, we explore the behaviour of a stochastic model of asymmetric competition where stochasticity is assumed to be demographic. We derive approximations for the temporal variance and covariance of the population sizes of the coexisting species. The derivations highlight that the variability of the population size of a species strongly depends on the stochastic fluctuations of species with higher trait values, while they are less influenced by species with lower trait values. Particularly, species with intermediate trait values are strongly affected resulting in relatively high variability. As a result these species have a relative high probability of extinction even though they have a larger population size than species with high trait values. We confirm these approximations with individual-based simulations. Thus, our analysis can explain gaps in size distributions as an emergent property of systems with a fecundity–competition trade-off.     

How does the magnitude of genetic variation affect ecological and reproductive character displacement?

While previous studies on character displacement tended to focus on trait divergence and convergence as a result of long-term evolution, recent studies suggest that character displacement can be a special case of evolutionary rescue, where rapid evolution prevents species extinction by weakening interspecific competition. Here we analyzed a simple model to examine how the magnitude of genetic variation affects evolutionary rescue via ecological and reproductive character displacement that weakens interspecific competition in exploitation of shared resources (i.e., resource competition) and in the mating process caused by incomplete species recognition (i.e., reproductive interference), respectively. We found that slow trait divergence due to small genetic variance results in species extinction in reproductive character displacement but not in ecological character displacement. This is because one species becomes rare in slow character displacement, and this causes deterministic extinction due to minority disadvantage of reproductive interference. On the other hand, there is no deterministic extinction in the process of ecological character displacement. Furthermore, species extinction becomes less likely in the case of positive covariance between ecological and reproductive traits as divergence of the ecological trait (e.g., root depths) increases the divergence speed of the reproductive trait (e.g., flower colors) and vice versa. It will be interesting to compare intraspecific genetic (co)variance of ecological and reproductive traits in future studies for understanding how ecological and reproductive character displacement occur without extinction.