Evolutionary dynamics through multispecies competition

Disruptive selection, emerging from frequency-dependent intraspecific competition can have very exciting evolutionary outcomes. One such outcome is the origin of new species through an evolutionary branching event. Literature on theoretical models investigating the emergence of disruptive selection is vast, with some investigating the sensitivity of the models on assumptions of the competition and carrying capacity functions’ shapes. What is seldom modeled is what happens once the population escapes its effect via increase phenotypic or genotypic variance. The expectation is mixed: disruptive selection could diminish and ultimately disappear or it could still exist leading to further speciation events through multiple evolutionary branching events. Here, we derive the conditions under which disruptive selection drives two subpopulations that originated at a branching point to other points in trait space where each subpopulation again experiences disruptive selection. We show that the general pattern for further branchings require that the competition function to be even narrower than what is required for the first evolutionary branching. However, we also show that the existence of disruptive selection in higher dimensional systems is also sensitive to the shapes of the functions used.     

Determinants of the strength of disruptive and/or divergent selection arising from resource competition

We investigate how the intensity of competition for resources affects the strength of disruptive selection on a resource acquisition trait. This is done by analyzing several consumer–resource models in which consumers use a linear array of resources. We show that disruptive selection can be diminished under both strong and weak competition, making disruptive selection a unimodal function of the strength of competition. Weak selection under strong competition arises when competition causes the extinction (for self-reproducing resources) or depletion (for abiotic resources) of the most rapidly caught resources. Weak selection under weak competition is a consequence of minimal effects of consumers on resources. The precise relationship between intensity of competition and strength of disruptive selection is sensitive to the shape of the consumer's resource utilization curve and the nature of resource growth. The most strongly unimodal competition–selection relationships result from utilization curves with long tails. Our results show that a simple comparison of the width of the resource abundance distribution and the consumer's utilization function is not sufficient to determine whether selection is disruptive. The results may explain some contradictory experimental findings regarding the effect of consumer mortality on the strength of disruptive selection.     

Additive genetic variation under intraspecific competition and stabilizing selection: a two-locus study

A diallelic two-locus model is investigated in which the loci determine the genotypic value of a quantitative trait additively. Fitness has two components: stabilizing selection on the trait and a frequency-dependent component, as induced, for instance, if the ability to utilize different food resources depends on this trait. Since intraspecific competition induces disruptive selection, this model leads to a conflict of selective forces. We study how the underlying genetics (recombination rate and allelic effects) interacts with the selective forces, and explore the resulting equilibrium structure. For the special case of equal effects, global stability results are proved. Unless the locus effects are sufficiently different, the genetic variance maintained at equilibrium displays a threshold-like dependence on the strength of competition. For loci with equal effects, the equilibrium fitnesses of genotypic values exhibit disruptive selection if and only if competition is strong enough to maintain a stable two-locus polymorphism. For unequal effects, disruptive selection can be observed for weaker competition and in the absence of a stable polymorphism.     

The conditions for speciation through intraspecific competition

Abstract It has been shown theoretically that sympatric speciation can occur if intraspecific competition is strong enough to induce disruptive selection. However, the plausibility of the involved processes is under debate, and many questions on the conditions for speciation remain unresolved. For instance, is strong disruptive selection sufficient for speciation? Which roles do genetic architecture and initial composition of the population play? How strong must assortative mating be before a population can split in two? These are some of the issues we address here. We investigate a diploid multilocus model of a quantitative trait that is under frequency‐dependent selection caused by a balance of intraspecific competition and frequency‐independent stabilizing selection. This trait also acts as mating character for assortment. It has been established previously that speciation can occur only if competition is strong enough to induce disruptive selection. We find that speciation becomes more difficult for very strong competition, because then extremely strong assortment is required. Thus, speciation is most likely for intermediate strengths of competition, where it requires strong, but not extremely strong, assortment. For this range of parameters, however, it is not obvious how assortment can evolve from low to high levels, because with moderately strong assortment less genetic variation is maintained than under weak or strong assortment sometimes none at all. In addition to the strength of frequency‐dependent competition and assortative mating, the roles of the number of loci, the distribution of allelic effects, the initial conditions, costs to being choosy, the strength of stabilizing selection, and the particular choice of the fitness function are explored. A multitude of possible evolutionary outcomes is observed, including loss of all genetic variation, splitting in two to five species, as well as very short and extremely long stable limit cycles. On the methodological side, we propose quantitative measures for deciding whether a given distribution reflects two (or more) reproductively isolated clusters.     

Competition and the Effect of Spatial Resource Heterogeneity on Evolutionary Diversification

A model is presented to explore how the form of selection arising from competition for resources is affected by spatial resource heterogeneity. The model consists of a single species occupying two patches connected by migration, where the two patches can differ in the type of resources that they contain. The main goal is to determine the conditions under which competition for resources results in disruptive selection (i.e., selection favoring a polymorphism) since it is this form of selection that will give rise to the evolutionary diversification of resource exploitation strategies. In particular, comparing the conditions giving rise to disruptive selection when the two patches are identical to the conditions when they contain different resources reveals the effect of spatial resource heterogeneity. Results show that when the patches are identical, the conditions giving rise to disruptive selection are identical to those that give rise to character displacement in previous models. When the patches are different, the conditions giving rise to disruptive selection can be either more or less stringent depending upon demographic parameters such as the intrinsic rate of increase and the migration rate. Surprisingly, spatial resource heterogeneity can actually make forms of evolutionary diversification such as character displacement less likely. It is also found that results are dependent on how the resource exploitation strategies and the spatial resource heterogeneity affect the population dynamics. One robust conclusion however, is that spatial resource heterogeneity always has a disruptive effect when the migration rate between patches is low.     

Sexual dimorphism and adaptive speciation: two sides of the same ecological coin

Models of adaptive speciation are typically concerned with demonstrating that it is possible for ecologically driven disruptive selection to lead to the evolution of assortative mating and hence speciation. However, disruptive selection could also lead to other forms of evolutionary diversification, including ecological sexual dimorphisms. Using a model of frequency-dependent intraspecific competition, we show analytically that adaptive speciation and dimorphism require identical ecological conditions. Numerical simulations of individual-based models show that a single ecological model can produce either evolutionary outcome, depending on the genetic independence of male and female traits and the potential strength of assortative mating. Speciation is inhibited when the genetic basis of male and female ecological traits allows the sexes to diverge substantially. This is because sexual dimorphism, which can evolve quickly, can eliminate the frequency-dependent disruptive selection that would have provided the impetus for speciation. Conversely, populations with strong assortative mating based on ecological traits are less likely to evolve a sexual dimorphism because females cannot simultaneously prefer males more similar to themselves while still allowing the males to diverge. This conflict between speciation and dimorphism can be circumvented in two ways. First, we find a novel form of speciation via negative assortative mating, leading to two dimorphic daughter species. Second, if assortative mating is based on a neutral marker trait, trophic dimorphism and speciation by positive assortative mating can occur simultaneously. We conclude that while adaptive speciation and ecological sexual dimorphism may occur simultaneously, allowing for sexual dimorphism restricts the likelihood of adaptive speciation. Thus, it is important to recognize that disruptive selection due to frequency-dependent interactions can lead to more than one form of adaptive splitting.     

The evolution of genetic architecture under frequency-dependent disruptive selection

We propose a model to analyze a quantitative trait under frequency-dependent disruptive selection. Selection on the trait is a combination of stabilizing selection and intraspecific competition, where competition is maximal between individuals with equal phenotypes. In addition, there is a density-dependent component induced by population regulation. The trait is determined additively by a number of biallelic loci, which can have different effects on the trait value. In contrast to most previous models, we assume that the allelic effects at the loci can evolve due to epistatic interactions with the genetic background. Using a modifier approach, we derive analytical results under the assumption of weak selection and constant population size, and we investigate the full model by numerical simulations. We find that frequency-dependent disruptive selection favors the evolution of a highly asymmetric genetic architecture, where most of the genetic variation is concentrated on a small number of loci. We show that the evolution of genetic architecture can be understood in terms of the ecological niches created by competition. The phenotypic distribution of a population with an adapted genetic architecture closely matches this niche structure. Thus, evolution of the genetic architecture seems to be a plausible way for populations to adapt to regimes of frequency-dependent disruptive selection. As such, it should be seen as a potential evolutionary pathway to discrete polymorphisms and as a potential alternative to other evolutionary responses, such as the evolution of sexual dimorphism or assortative mating.     

Oligomorphic dynamics for analyzing the quantitative genetics of adaptive speciation

Ecological interaction, including competition for resources, often causes frequency-dependent disruptive selection, which, when accompanied by reproductive isolation, may act as driving forces of adaptive speciation. While adaptive dynamics models have added new perspectives to our understanding of the ecological dimensions of speciation processes, it remains an open question how best to incorporate and analyze genetic detail in such models. Conventional approaches, based on quantitative genetics theory, typically assume a unimodal character distribution and examine how its moments change over time. Such approximations inevitably fail when a character distribution becomes multimodal. Here, we propose a new approximation, oligomorphic dynamics, to the quantitative genetics of populations that include several morphs and that therefore exhibit multiple peaks in their character distribution. To this end, we first decompose the character distribution into a sum of unimodal distributions corresponding to individual morphs. Characterizing these morphs by their frequency (fraction of individuals belonging to each morph), position (mean character of each morph), and width (standard deviation of each morph), we then derive the coupled eco-evolutionary dynamics of morphs through a double Taylor expansion. We show that the demographic, convergence, and evolutionary stability of a population’s character distribution correspond, respectively, to the asymptotic stability of frequencies, positions, and widths under the oligomorphic dynamics introduced here. As first applications of oligomorphic dynamics theory, we analytically derive the effects (a) of the strength of disruptive selection on waiting times until speciation, (b) of mutation on conditions for speciation, and (c) of the fourth moments of competition kernels on patterns of speciation.     

On a genetic model of intraspecific competition and stabilizing selection

A genetic model is investigated in which two recombining loci determine the genotypic value of a quantitative trait additively. Two opposing evolutionary forces are assumed to act: stabilizing selection on the trait, which favors genotypes with an intermediate phenotype, and intraspecific competition mediated by that trait, which favors genotypes whose effect on the trait deviates most from that of the prevailing genotypes. Accordingly, fitnesses of genotypes have a frequency-independent component describing stabilizing selection and a frequency- and density-dependent component modeling competition. We study how the underlying genetics, in particular recombination rate and relative magnitude of allelic effects, interact with the conflicting selective forces and derive the resulting, surprisingly complex equilibrium patterns. We also investigate the conditions under which disruptive selection on the phenotypes can be observed and examine how much genetic variation can be maintained in such a model. We discovered a number of unexpected phenomena. For instance, we found that with little recombination the degree of stably maintained polymorphism and the equilibrium genetic variance can decrease as the strength of competition increases relative to the strength of stabilizing selection. In addition, we found that mean fitness at the stable equilibria is usually much lower than the maximum possible mean fitness and often even lower than the fitness at other, unstable equilibria. Thus, the evolutionary dynamics in this system are almost always nonadaptive.     

Predictable patterns of disruptive selection in stickleback in postglacial lakes

Disruptive selection is often assumed to be relatively rare, because it is dynamically unstable and hence should be transient. However, frequency-dependent interactions such as intraspecific competition may stabilize fitness minima and make disruptive selection more common. Such selection helps explain the maintenance of genetic variation and may even contribute to sympatric speciation. There is thus great interest in determining when and where disruptive selection is most likely. Here, we show that there is a general trend toward weak disruptive selection on trophic morphology in three-spine stickleback (Gasterosteus aculeatus) in 14 lakes on Vancouver Island. Selection is inferred from the observation that, within a lake, fish with intermediate gill raker morphology exhibited slower growth than phenotypically extreme individuals. Such selection has previously been shown to arise from intraspecific competition for alternate resources. However, not all environments are equally conducive to disruptive selection, which was strongest in intermediate-sized lakes where both littoral and pelagic prey are roughly balanced. Also, consistent with theory, we find that sexual dimorphism in trophic traits tends to mitigate disruptive selection. These results suggest that it may be possible to anticipate the kinds of environments and populations most likely to experience disruptive selection.     

Speciation through competition: a critical review

We examined causes of speciation in asexual populations in both sympatry and parapatry, providing an alternative explanation for the speciation patterns reported by Dieckmann and Doebeli (1999) and Doebeli and Dieckmann (2003). Both in sympatry and parapatry, they find that speciation occurs relatively easily. We reveal that in the sympatric clonal model, the equilibrium distribution is continuous and the disruptive selection driving evolution of discrete clusters is only transient. Hence, if discrete phenotypes are to remain stable in the sympatric sexual model, there should be some source of nontransient disruptive selection that will drive evolution of assortment. We analyze sexually reproducing populations using the Bulmer's infinitesimal model and show that cost-free assortment alone leads to speciation and disruptive selection only arises when the optimal distribution cannot be matched--in this example, because the phenotypic range is limited. In addition, Doebeli and Dieckmann's analyses assumed a high genetic variance and a high mutation rate. Thus, these theoretical models do not support the conclusion that sympatric speciation is a likely outcome of competition for resources. In their parapatric model (Doebeli and Dieckmann 2003), clustering into distinct phenotypes is driven by edge effects, rather than by frequency-dependent competition.     

Competition and evolution along environmental gradients: patterns, boundaries and sympatric divergence

The present study examined how competitive interactions and environmental conditions generate species boundaries and determine species distributions. A spatially explicit, quantitative genetic, two-species competition model was used to manipulate the strengths of competition, gene flow and local adaptation along environmental gradients. This allowed us to assess the long-term persistence of each species and whether the ranges they inhabited had boundaries in space or were unlimited. We found that a species boundary arises along less steep environmental gradients when the strength of stabilizing selection and diversifying selection are similar. We also found that a species boundary may arise along shallow environmental gradients if interspecific competition is more intense than intraspecific, which relaxes previous requirements for steep gradients for generating range limits. We determined an analytical form for the critical environmental gradient as a function of ecological and genetic parameters at which a species boundary is expected to arise by competition. Results suggest an alternative to resource competition as an explanation for phenotypic divergence between sympatric competitors. Competitors sharing a trait that is under stabilizing selection along an environmental gradient may segregate spatially and evolve in different regions, with phenotypic sympatric divergence reflecting the resulting clines. Along various types of environmental gradients, variation in stabilizing selection intensities could lead to contrasting patterns in the distribution of species. For stabilizing selection strengths in accord with field data estimates, this study predicts that the level of sympatric character divergence would be limited along environmental gradients.     

Variation,selection and evolution of function-valued traits

We describe an emerging framework for understanding variation, selection and evolution of phenotypic traits that are mathematical functions. We use one specific empirical example – thermal performance curves (TPCs) for growth rates of caterpillars – to demonstrate how models for function-valued traits are natural extensions of more familiar, multivariate models for correlated, quantitative traits. We emphasize three main points. First, because function-valued traits are continuous functions, there are important constraints on their patterns of variation that are not captured by multivariate models. Phenotypic and genetic variation in function-valued traits can be quantified in terms of variance-covariance functions and their associated eigenfunctions: we illustrate how these are estimated as well as their biological interpretations for TPCs. Second, selection on a function-valued trait is itself a function, defined in terms of selection gradient functions. For TPCs, the selection gradient describes how the relationship between an organism's performance and its fitness varies as a function of its temperature. We show how the form of the selection gradient function for TPCs relates to the frequency distribution of environmental states (caterpillar temperatures) during selection. Third, we can predict evolutionary responses of function-valued traits in terms of the genetic variance-covariance and the selection gradient functions. We illustrate how non-linear evolutionary responses of TPCs may occur even when the mean phenotype and the selection gradient are themselves linear functions of temperature. Finally, we discuss some of the methodological and empirical challenges for future studies of the evolution of function-valued traits.     

A multilocus analysis of intraspecific competition and stabilizing selection on a quantitative trait

The equilibrium structure of an additive, diallelic multilocus model of a quantitative trait under frequency- and density-dependent selection is derived. The trait is under stabilizing selection and mediates intraspecific competition as induced, for instance, by differential resource utilization. It is assumed that stabilizing selection is weak, but the strength of competition may be arbitrary relative to it. Density dependence is caused by population regulation, which may be of a very general kind. The number and effects of loci are arbitrary, and stabilizing selection is not necessarily symmetric with respect to the range of phenotypic values. All previously studied models of intraspecific competition for a continuum of resources known to the author reduce to a special case of the present model if overall selection is weak. Therefore, in this case our results are applicable as approximations to all these models. Our central result is the (nearly) complete characterization of the equilibrium and stability structure in terms of all parameters. It is derived under the sole assumption that selection is weak enough relative to recombination to ignore linkage disequilibrium. In particular, necessary and sufficient conditions on the strength of competition relative to stabilizing selection are found that ensure the maintenance of multilocus polymorphism and the occurrence of disruptive selection. In this case, explicit formulas for the number of polymorphic loci at equilibrium, the allele frequencies, the genetic variance, and the strength of disruptive selection are obtained. For two loci, the effects of linkage are investigated analytically; for several loci, they are studied numerically.     

Genetics of male nuptial colour divergence between sympatric sister species of a Lake Victoria cichlid fish

The hypothesis of sympatric speciation by sexual selection has been contentious. Several recent theoretical models of sympatric speciation by disruptive sexual selection were tailored to apply to African cichlids. Most of this work concludes that the genetic architecture of female preference and male trait is a key determinant of the likelihood of disruptive sexual selection to result in speciation. We investigated the genetic architecture controlling male nuptial colouration in a sympatric sibling species pair of cichlid fish from Lake Victoria, which differ conspicuously in male colouration and female mating preferences for these. We estimated that the difference between the species in male nuptial red colouration is controlled by a minimum number of two to four genes with significant epistasis and dominance effects. Yellow colouration appears to be controlled by one gene with complete dominance. The two colours appear to be epistatically linked. Knowledge on how male colouration segregates in hybrid generations and on the number of genes controlling differences between species can help us assess whether assumptions made in simulation models of sympatric speciation by sexual selection are realistic. In the particular case of the two sister species that we studied a small number of genes causing major differences in male colouration may have facilitated the divergence in male colouration associated with speciation.     

Density-Dependent Selection Incorporating Intraspecific Competition. II. a Diploid Model

A diploid model is introduced and analyzed in which intraspecific competition is incorporated within the context of density-regulated selection. It is assumed that each genotype has a unique carrying capacity corresponding to the equilibrium population size when only that type is present. Each genotypic fitness at a single diallelic autosomal locus is a decreasing function of a distinctive effective population size perceived as a result of intraspecific competition. The resulting fitnesses are both density and frequency dependent with selective advantage determined by a balance between genotypic carrying capacity and sensitivity to intraspecific competition. A major finding is that intergenotypic interactions may allow genetic variation to be more easily maintained than in the corresponding model of purely density-dependent selection. In addition, numerical study confirms the possible existence of multiple interior equilibria and that neither overdominance in fitness nor carrying capacity is necessary for stability. The magnitude of the equilibrium population size and optimization principles are also discussed.     

Estimating the distribution of selection coefficients from phylogenetic data with applications to mitochondrial and viral DNA

The distribution of selection coefficients of new mutations is of key interest in population genetics. In this paper we explore how codon-based likelihood models can be used to estimate the distribution of selection coefficients of new amino acid replacement mutations from phylogenetic data. To obtain such estimates we assume that all mutations at the same site have the same selection coefficient. We first estimate the distribution of selection coefficients from two large viral data sets under the assumption that the viral population size is the same along all lineages of the phylogeny and that the selection coefficients vary among sites. We then implement several new models in which the lineages of the phylogeny may have different population sizes. We apply the new models to a data set consisting of the coding regions from eight primate mitochondrial genomes. The results suggest that there might be little power to determine the exact shape of the distribution of selection coefficient but that the normal and gamma distributions fit the data significantly better than the exponential distribution.     

Under which conditions is character displacement a likely outcome of secondary contact?

Sympatric character displacement is one possible mechanism that prevents competitive exclusion. This mechanism is thought to be behind the radiation of Darwin's finches, where character displacement is assumed to have followed secondary contact of ecologically similar species. We use a model to evaluate under which ecological and environmental conditions this mechanism is likely. Using the adaptive dynamics theory, we analyse different ecological models embedded in the secondary contact scenario. We highlight two necessary conditions for character displacement in sympatry: (i) very strong premating isolation between the two populations, and (ii) secondary contact to occur at an evolutionary branching point. Character displacement is then driven by adaptation to interspecific competition. We determine how ecological and environmental parameters influence the probability of ecological divergence. Finally, we discuss the likelihood of sympatric character displacement under disruptive selection in natural populations.     

Inference in MCMC step selection models

Habitat selection models are used in ecology to link the spatial distribution of animals to environmental covariates and identify preferred habitats. The most widely used models of this type, resource selection functions, aim to capture the steady-state distribution of space use of the animal, but they assume independence between the observed locations of an animal. This is unrealistic when location data display temporal autocorrelation. The alternative approach of step selection functions embed habitat selection in a model of animal movement, to account for the autocorrelation. However, inferences from step selection functions depend on the underlying movement model, and they do not readily predict steady-state space use. We suggest an analogy between parameter updates and target distributions in Markov chain Monte Carlo (MCMC) algorithms, and step selection and steady-state distributions in movement ecology, leading to a step selection model with an explicit steady-state distribution. In this framework, we explain how maximum likelihood estimation can be used for simultaneous inference about movement and habitat selection. We describe the local Gibbs sampler, a novel rejection-free MCMC scheme, use it as the basis of a flexible class of animal movement models, and derive its likelihood function for several important special cases. In a simulation study, we verify that maximum likelihood estimation can recover all model parameters. We illustrate the application of the method with data from a zebra.     

Density compensation,isocline shape and single-level competition models

Criticism of the Lotka-Volterra competition model implies that the theory of competition should be based upon more general concepts. It is suggested that the shape of the competitor isoclines can provide this basis.The relationship between the total density of a competitive community and species number depends crucially upon isocline shape. This has immediate relevance to the interpretation of the excess density compensation seen in some island communities, since if isoclines are sufficiently concave (curved towards the origin) then this phenomenon is expected to be the rule rather than the exception.These observations do not depend upon any specific model, but in order to determine the shape of isoclines in natural communities a link must be made between the biological processes of competition and isocline shape. To this end three types of single-level competition model are distinguished (additive, multiplicative and temporal resource models) depending upon how gains from resources interact in determining individual fitness.The models are based upon resource availability functions (RAFs), which are decreasing functions of the level of competition and determine the availability of each resource to each species. Provided that the argument of these functions is always a weighted sum of the number of competitors then in the case of the additive resource model the shape of the RAFs determines directly the shape of the isoclines. For the multiplicative model, the shape of the logarithm of the RAFs adopts this role.Analysis of a special case of the additive resource model suggests that concave isoclines are likely to predominate, and that the degree of concavity is of an order which minimizes the tendency of total numbers to increase with species number. In some circumstances, involving “scramble”-type competition and habitat selection, the expected concavity is sufficient to cause a decrease. In any event the expected occurrence of concave competition isoclines predicts a much higher incidence of excess density compensation (due to carrying capacity differences) than expected from any model having linear isoclines.The effect of shifting from an additive to a multiplicative resource model is to make the existence of purely concave isoclines less probable and to raise the possibility of purely convex isoclines. On the other hand, shifting from an additive resource model to a temporal resource model apparently has no such simple interpretation and specific predictions must await further analysis.