Case studies and mathematical models of ecological speciation. 2. Palms on an oceanic island

A recent study of a pair of sympatric species of palms on the Lord Howe Island is viewed as providing probably one of the most convincing examples of sympatric speciation to date. Here we describe and study a stochastic, individual-based, explicit genetic model tailored for this palms system. Overall, our results show that relatively rapid (<50,000 generations) colonization of a new ecological niche, and sympatric or parapatric speciation via local adaptation and divergence in flowering periods are theoretically plausible if (i) the number of loci controlling the ecological and flowering period traits is small; (ii) the strength of selection for local adaptation is intermediate; and (iii) an acceleration of flowering by a direct environmental effect associated with the new ecological niche is present. We discuss patterns and time-scales of ecological speciation identified by our model, and we highlight important parameters and features that need to be studied empirically in order to provide information that can be used to improve the biological realism and power of mathematical models of ecological speciation.     

Modelling sympatric speciation by means of biologically plausible mechanistic processes as exemplified by threespine stickleback species pairs

We investigate the plausibility of sympatric speciation through a modelling study. We built up a series of models with increasing complexity while focussing on questioning the realism of model assumptions by checking them critically against a particular biological system, namely the sympatric benthic and limnetic species of threespine stickleback in British Columbia, Canada. These are morphologically adapted to their feeding habits: each performs better in its respective habitat than do hybrids with intermediate morphology. Ecological character displacement through disruptive selection and competition, and reinforcement through mating preferences may have caused their divergence. Our model assumptions include continuous morphological trait(s) instead of a dimorphic trait, and mating preferences based on the same trait(s) as selected for in food competition. Initially, morphology is intermediate. We apply disruptive selection against intermediates, frequency-dependent resource competition, and one of two alternative mating preference mechanisms. Firstly, preference is based on similarity where mating preference may result from “imprinting” on conspecifics encountered in their preferred foraging habitat. Here, speciation occurs easily—ecological hybrid inferiority is not necessary. Hybrid inferiority reinforces the stringency of assortative mating. Secondly, individual preferences exist for different trait values. Here, speciation occurs when linkage disequilibrium between trait and preference develops, and some hybrid inferiority is required. Finally, if the morphology subject to disruptive selection, frequency-dependent competition, and mate choice, is coded for by two loci, linkage disequilibrium between the two loci is required for speciation. Speciation and reinforcement of stringency of choosiness are possible in this case too, but rarely. Results demonstrate the contingency of speciation, with the same starting point not necessarily producing the same outcome. The study resulted in flagging issues where models often lack in biological realism and issues where more empirical studies could inform on whether assumptions are likely valid.     

Host choice promotes reproductive isolation between host races of the larch budmoth Zeiraphera diniana

The chances for sympatric speciation are improved if ecological divergence leads to assortative mating as a by-product. This effect is known in parasites that find mates using host cues, but studies of larch- and pine-feeding races of the larch budmoth (Zeiraphera diniana, Lepidoptera: Tortricidae) suggest it may also occur when mate attraction is via sex pheromones that are independent of habitat. We have previously shown that females releasing pheromones on or near their own host attract more males of their own race than if placed on the alternative host. This host effect would enhance assortative mating provided adults preferentially alight on their native hosts. Here we investigate alighting preferences in natural mixed forest using a novel likelihood analysis of genotypic clusters based on three semidiagnostic allozyme loci. Both larch and pine females show a realized alighting preference for their own host of 86%. The equivalent preferences of males were 79% for the larch race and 85% for the pine race. These preferences are also detectable in small-scale laboratory experiments, where alighting preferences of larch and pine races towards their own hosts were, respectively, 67 and 66% in females and 69 and 63% in males. Pure larch race moths reared in the laboratory had alighting choice similar to moths from natural populations, while hybrids were intermediate, showing that alighting preferences were heritable and approximately additive. The field estimates of alighting preference, coupled with earlier work on mate choice, yield an estimated rate of natural hybridization between sympatric host races of 2.2-3.8% per generation. Divergent alighting choice enhances pheromone-mediated assortative mating today, and is likely to have been an important cause of assortative mating during initial divergence in host use. Because resources are normally 'coarse-grained' in space and time, assortative mating due to ecological divergence may be a more important catalyst of sympatric speciation than generally realized.     

Comparing the consequences of natural selection,adaptive phenotypic plasticity,and matching habitat choice for phenotype–environment matching,population genetic structure,and reproductive isolation in meta‐populations

Organisms commonly experience significant spatiotemporal variation in their environments. In response to such heterogeneity, different mechanisms may act that enhance ecological performance locally. However, depending on the nature of the mechanism involved, the consequences for populations may differ greatly. Building on a previous model that investigated the conditions under which different adaptive mechanisms (co)evolve, this study compares the ecological and evolutionary population consequences of three very different responses to environmental heterogeneity: matching habitat choice (directed gene flow), adaptive plasticity (associated with random gene flow), and divergent natural selection. Using individual‐based simulations, we show that matching habitat choice can have a greater adaptive potential than plasticity or natural selection: it allows for local adaptation while protecting genetic polymorphism despite global mating or strong environmental changes. Our simulations further reveal that increasing environmental fluctuations and unpredictability generally favor the emergence of specialist genotypes but that matching habitat choice is better at preventing local maladaptation by individuals. This confirms that matching habitat choice can speed up the genetic divergence among populations, cause indirect assortative mating via spatial clustering, and hence even facilitate sympatric speciation. This study highlights the potential importance of directed dispersal in local adaptation and speciation, stresses the difficulty of deriving its operation from nonexperimental observational data alone, and helps define a set of ecological conditions which should favor its emergence and subsequent detection in nature.     

Sympatric Speciation: When Is It Possible in Bacteria?

According to theory, sympatric speciation in sexual eukaryotes is favored when relatively few loci in the genome are sufficient for reproductive isolation and adaptation to different niches. Here we show a similar result for clonally reproducing bacteria, but which comes about for different reasons. In simulated microbial populations, there is an evolutionary tradeoff between early and late stages of niche adaptation, which is resolved when relatively few loci are required for adaptation. At early stages, recombination accelerates adaptation to new niches (ecological speciation) by combining multiple adaptive alleles into a single genome. Later on, without assortative mating or other barriers to gene flow, recombination generates unfit intermediate genotypes and homogenizes incipient species. The solution to this tradeoff may be simply to reduce the number of loci required for speciation, or to reduce recombination between species over time. Both solutions appear to be relevant in natural microbial populations, allowing them to diverge into ecological species under similar constraints as sexual eukaryotes, despite differences in their life histories.     

Sensory drive speciation and patterns of variation at selectively neutral genes

Speciation by sensory drive can occur if divergent adaptation of sensory systems causes rapid evolution of mating traits and the resulting development of assortative mating. Previous theoretical studies have shown that sensory drive can cause rapid divergent adaptive evolution from one to two phenotypes. In this study, we examined two topics: the possibility of adaptive radiation by sensory drive from one to more than two phenotypes and the relationships of patterns of variation at selectively neutral genes to levels of viability selection, habitat and mating preferences and migration. We conducted individual-based simulations assuming a sensory trait and a mating trait controlled by a small number of loci. We found that adaptive radiation is possible when the number of loci controlling the sensory trait is small; the levels of viability selection, habitat and mating preferences are intermediate; and the emigration rate is high. We also found that emigration rates as well as the levels of habitat and mating preferences are related to F _ST values at neutral loci, but F _ST proved to be insensitive to a small change in the number of loci controlling the mating trait. This suggests that an estimation of the past population history is possible without an accurate genetic model.     

Inheritance of female mating preference in a sympatric sibling species pair of Lake Victoria cichlids: implications for speciation

Female mate choice has often been proposed to play an important role in cases of rapid speciation, in particular in the explosively evolved haplochromine cichlid species flocks of the Great Lakes of East Africa. Little, if anything, is known in cichlid radiations about the heritability of female mating preferences. Entirely sympatric distribution, large ecological overlap and conspicuous differences in male nuptial coloration, and female preferences for these, make the sister species Pundamilia pundamilia and P. nyererei from Lake Victoria an ideally suited species pair to test assumptions on the genetics of mating preferences made in models of sympatric speciation. Female mate choice is necessary and sufficient to maintain reproductive isolation between these species, and it is perhaps not unlikely therefore, that female mate choice has been important during speciation. A prerequisite for this, which had remained untested in African cichlid fish, is that variation in female mating preferences is heritable. We investigated mating preferences of females of these sister species and their hybrids to test this assumption of most sympatric speciation models, and to further test the assumption of some models of sympatric speciation by sexual selection that female preference is a single-gene trait. We find that the differences in female mating preferences between the sister species are heritable, possibly with quite high heritabilities, and that few but probably more than one genetic loci contribute to this behavioural speciation trait with no apparent dominance. We discuss these results in the light of speciation models and the debate about the explosive radiation of cichlid fishes in Lake Victoria.     

Case studies and mathematical models of ecological speciation. 3: Ecotype formation in a Swedish snail

Formation of partially reproductively isolated ecotypes in the rough periwinkle, Littorina saxatilis , may be a case of incipient nonallopatric ecological speciation. To better understand the dynamics of ecotype formation, its timescale, driving forces and evolutionary consequences, we developed a spatially explicit, individual-based model incorporating relevant ecological, spatial and mate selection data for Swedish L. saxatilis . We explore the impact of bounded hybrid superiority, ecological scenarios and mate selection systems on ecotype formation, gene flow and the evolution of prezygotic isolation. Our model shows that ecotypes are expected to form rapidly in parapatry under conditions applicable to Swedish L. saxatilis and may proceed to speciation. However, evolution of nonrandom mating had complex behaviour. Ecotype evolution was inhibited by pre-existing mating preferences, but facilitated by the evolution of novel preferences. While in many scenarios positive assortative mating reduced gene flow between ecotypes, in others negative assortative mating arose, preferences were lost after ecotype formation, preferences were confined to one ecotype or the ancestral ecotype became extinct through sexual selection. Bounded hybrid superiority (as observed in nature) enhanced ecotype formation but increased gene flow. Our results highlight that ecotype formation and speciation are distinct processes: factors that contribute to ecotype formation can be detrimental to speciation and vice versa. The complex interactions observed between local adaptation and nonrandom mating imply that generalization from data is unreliable without quantitative theory for speciation.     

Adaptive,but not condition‐dependent,body shape differences contribute to assortative mating preferences during ecological speciation

Assortative mating is critical for reproductive isolation during speciation; however, the mechanisms underlying mating preferences are often unknown. Assortative mating can be mediated through preferences for condition‐dependent and adaptive (“magic”) traits, but rigorously testing these hypotheses has been impeded by trait covariation in living organisms. We used computer‐generated models to examine the role of body shape in producing association preferences between fish populations undergoing ecological speciation in different habitat types. We demonstrate that body shape can serve as an adaptive trait (variation in head size between populations) and a condition‐dependent signal (variation in abdominal distention among individuals). Female preferences for stimuli varying in only one aspect of body shape uncovered evidence for body shape as a magic trait across population pairs, but no evidence for body shape serving as a condition‐dependent signal. Evolution of preferences only in females from one habitat type as well as stronger preferences in sympatric nonsulfidic as opposed to allopatric nonsulfidic populations suggests that reinforcement may have played a role in producing the observed patterns.     

CASE STUDIES AND MATHEMATICAL MODELS OF ECOLOGICAL SPECIATION. 4. HYBRID SPECIATION IN BUTTERFLIES IN A JUNGLE

We build a spatial individual-based multilocus model of homoploid hybrid speciation tailored for a tentative case of hybrid origin of Heliconius heurippa from H. melpomene and H. cydno in South America. Our model attempts to account for empirical patterns and data on genetic incompatibility, mating preferences and selection by predation (both based on coloration patterns), habitat preference, and local adaptation for all three Heliconius species. Using this model, we study the likelihood of recombinational speciation and identify the effects of various ecological and genetic parameters on the dynamics, patterns, and consequences of hybrid ecological speciation. Overall, our model supports the possibility of hybrid origin of H. heurippa under certain conditions. The most plausible scenario would include hybridization between H. melpomene and H. cydno in an area geographically isolated from the rest of both parental species with subsequent long-lasting geographic isolation of the new hybrid species, followed by changes in the species ranges, the secondary contact, and disappearance of H. melpomene -type ecomorph in the hybrid species. However, much more work (both empirical and theoretical) is necessary to be able to make more definite conclusions on the importance of homoploid hybrid speciation in animals.     

Adaptive sympatric speciation of polychromatic "roundfin" sailfin silverside fish in Lake Matano (Sulawesi)

The significance of sympatric speciation is one of the most controversial topics in evolutionary biology. Theory suggests that different factors can lead to speciation in full geographical contact, including selection and nonrandom mating. Strict criteria have been established for assessing sympatric speciation, which have been met in only a very few cases. Here, we investigate differentiation among sympatric morphospecies and color morphs of "roundfin" sailfin silversides (Telmatherinidae), small freshwater fish endemic to ancient Lake Matano in Central Sulawesi (Indonesia). Morphospecies are distinct according to body shape (geometric morphometrics), population structure (population-level amplified fragment length polymorphism [AFLP] markers), ecology, and mating behavior (habitat transects, stomach contents). Explorative genome scans based on AFLPs indicate that divergent selection affects only 1.3-4.2% of the analyzed loci, suggesting an early stage of speciation. Transect data demonstrate strong assortative mating and adaptive niche differentiation. However, we find no restrictions in gene flow among the conspicuous male color morphs. In summary, our data are consistent with a sympatric mode of divergence among three morphospecies under conditions effectively ruling out allopatric scenarios. Substantial, but incomplete, reproductive isolation suggests an early stage of speciation, most likely due to ecological selection pressure.     

A combination of sexual and ecological divergence contributes to rearrangement spread during initial stages of speciation

Chromosomal rearrangements between sympatric species often contain multiple loci contributing to assortative mating, local adaptation and hybrid sterility. When and how these associations arise during the process of speciation remains a subject of debate. Here, we address the relative roles of local adaptation and assortative mating on the dynamics of rearrangement evolution by studying how a rearrangement covaries with sexual and ecological trait divergence within a species. Previously, a chromosomal rearrangement that suppresses recombination on the Z (sex) chromosome was identified in European corn borer moths (Ostrinia nubilalis). We further characterize this recombination suppressor and explore its association with variation in sex pheromone communication and seasonal ecological adaptation in pairs of populations that are divergent in one or both of these characteristics. Direct estimates of recombination suppression in pedigree mapping families indicated that more than 39% of the Z chromosome (encompassing up to ~10 megabases and ~300 genes) resides within a nonrecombining unit, including pheromone olfactory receptor genes and a major quantitative trait locus that contributes to ecotype differences (Pdd). Combining direct and indirect estimates of recombination suppression, we found that the rearrangement was occasionally present between sexually isolated strains (E vs. Z) and between divergent ecotypes (univoltine vs. bivoltine). However, it was only consistently present when populations differed in both sexual and ecological traits. Our results suggest that independent of the forces that drove the initial establishment of the rearrangement, a combination of sexual and ecological divergence is required for rearrangement spread during speciation.     

Multilocus models of sympatric speciation: Bush versus Rice versus Felsenstein

Abstract. In populations of phytophagous insects that use the host plant as a rendezvous for mating, divergence in host preference could lead to sympatric speciation. Speciation requires the elimination of "generalist" genotypes, that is, those with intermediate host preference. This could occur because such genotypes have an inherent fitness disadvantage, or because preference alleles become associated with alleles that are oppositely selected on the two hosts. Although the former mechanism has been shown to be plausible, the latter mechanism has not been studied in detail. I consider a multilocus model (the "Bush model") in which one set of biallelic loci affects host preference, and a second set affects viability on the hosts once chosen. Alleles that increase viability on one host decrease viability on the other, and all loci are assumed to be unlinked. With moderately strong selection on the viability loci, preference alleles rapidly become associated with viability alleles, and the population splits into two reproductively isolated host specialist populations. The conditions for speciation to occur in this model, as measured by the strength of selection required, are somewhat more stringent than in a model in which preference and viability are controlled by the same loci (one-trait model). In contrast, the conditions are much less stringent than in a model in which speciation requires buildup of associations between viability loci and loci controlling a host-independent assortative mating trait (canonical two-trait model). Moreover, in the one-trait model, and to a lesser extent the Bush model, the strength of selection needed to initiate speciation is only slightly greater than that needed to complete it. This indicates that documenting instances of sympatric species that are reproductively isolated only by host or habitat preference would provide evidence for the plausibility of sympatric speciation in nature.     

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F _ST outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F _ST partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F _ST differed between candidate and noncandidate loci for all partitions except in the case of microgeographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.     

Speciation through the learning of habitat features

Learning of environmental features can influence both mating behaviour and the location where young are produced. This may lead to speciation in three steps: (i) colonization of a new habitat, (ii) genetic divergence of the two groups by adaptation to the habitats, and (iii) a decrease of genetic mixing between the lineages (similar to reinforcement). In a previous paper we showed that steps (i) and (ii) occur readily for a wide range of fixed mating and habitat preferences. Here, we study whether this can ultimately lead to speciation through selective changes in these preferences. We show that this indeed occurs, and, furthermore, it is very general: for a large class of models there is selection toward producing young more frequently in the natal habitat. Once habitat preference is strong, there is selection toward stronger assortative mating. Even when steps (i) and (ii) initially fail, genetic divergence may succeed at a later evolutionary stage, after which a decrease of genetic mixing completes speciation. Our results show that speciation by the learning of habitat features is an extremely effective mechanism.     

A TEST FOR SEXUAL SELECTION ON HYBRIDS OF TWO SYMPATRIC STICKLEBACKS

In this study we assessed whether sexual selection against hybrids contributes to reproductive isolation between two sympatric stickleback species. The species are recently diverged and possibly in the final stages of speciation. Our aim was to find whether mating discrimination of the parental species selects against F₁ hybrids, and what conditions are necessary for such sexual selection to operate. We used conservative no-choice laboratory trials with reproductively naive, lab-reared fish to measure female mating preferences. Females exhibited ranked preferences, preferring in order: conspecific, hybrid, then heterospecific males. However, intermediate attractiveness does not necessarily imply selection against hybrids: two-way ANOVAs suggested that limnetic, benthic, and hybrid males were statistically equivalent when averaged across females. Thus, this experiment found no evidence for a hybrid mating disadvantage. Our interpretation is that if sexual selection against hybrids is present in the wild, then some factor that biases encounter rates between hybrids and parental species (e.g., habitat selection) is necessary to produce it.     

Mechanisms of species divergence through visual adaptation and sexual selection: Perspectives from a cichlid model sys-tem

The theory of ecological speciation suggests that assortative mating evolves most easily when mating preferences are directly linked to ecological traits that are subject to divergent selection. Sensory adaptation can play a major role in this process, because selective mating is often mediated by sexual signals: bright colours, complex song, pheromone blends and so on. When divergent sensory adaptation affects the perception of such signals, mating patterns may change as an immediate consequence. Alternatively, mating preferences can diverge as a result of indirect effects: assortative mating may be promoted by selection against intermediate phenotypes that are maladapted to their (sensory) environment. For Lake Victoria cichlids, the visual envi-ronment constitutes an important selective force that is heterogeneous across geographical and water depth gradients. We investi-gate the direct and indirect effects of this heterogeneity on the evolution of female preferences for alternative male nuptial colours (red and blue) in the genus Pundamilia. Here, we review the current evidence for divergent sensory drive in this system, extract general principles, and discuss future perspectives.     

Refining the conditions for sympatric ecological speciation

Can speciation occur in a single population when different types of resources are available, in the absence of any geographical isolation, or any spatial or temporal variation in selection? The controversial topics of sympatric speciation and ecological speciation have already stimulated many theoretical studies, most of them agreeing on the fact that mechanisms generating disruptive selection, some level of assortment, and enough heterogeneity in the available resources, are critical for sympatric speciation to occur. Few studies, however, have combined the three factors and investigated their interactions. In this article, I analytically derive conditions for sympatric speciation in a general model where the distribution of resources can be uni‐ or bimodal, and where a parameter controls the range of resources that an individual can exploit. This approach bridges the gap between models of a unimodal continuum of resources and Levene‐type models with discrete resources. I then test these conditions against simulation results from a recently published article (Thibert‐Plante & Hendry, 2011, J. Evol. Biol. 24 : 2186–2196) and confirm that sympatric ecological speciation is favoured when (i) selection is disruptive (i.e. individuals with an intermediate trait are at a local fitness minimum), (ii) resources are differentiated enough and (iii) mating is assortative. I also discuss the role of mating preference functions and the need (or lack thereof) for bimodality in resource distributions for diversification.     

HABITAT AVOIDANCE: OVERLOOKING AN IMPORTANT ASPECT OF HOST‐SPECIFIC MATING AND SYMPATRIC SPECIATION?

Understanding speciation requires discerning how reproductive barriers to gene flow evolve between previously interbreeding populations. Models of sympatric speciation for phytophagous insects posit that reproductive isolation can evolve in the absence of geographic isolation as a consequence of an insect shifting and ecologically adapting to a new host plant. One important adaptation contributing to sympatric differentiation is host-specific mating. When organisms mate in preferred habitats, a system of positive assortative mating is established that facilitates sympatric divergence. Models of host fidelity generally assume that host choice is determined by the aggregate effect of alleles imparting positive preferences for different plant species. But negative effect genes for avoiding nonnatal plants may also influence host use. Previous studies have shown that apple and hawthorn-infesting races of Rhagoletis pomonella flies use volatile compounds emitted from the surface of fruit as key chemosensory cues to recognize and distinguish between their host plants. Here, we report results from field trials indicating that in addition to preferring the odor of their natal fruit, apple and hawthorn flies, and their undescribed sister species infesting flowering dogwood (Cornus florida), also avoid the odors of nonnatal fruit. We discuss the implications of nonnatal fruit avoidance for the evolutionary dynamics and genetics of sympatric speciation. Our findings reveal an underappreciated role for habitat avoidance as a potential postmating, as well as prezygotic, barrier to gene flow.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.