The effect of Wolbachia versus genetic incompatibilities on reinforcement and speciation

Wolbachia is a widespread group of intracellular bacteria commonly found in arthropods. In many insect species, Wolbachia induce a cytoplasmic mating incompatibility (CI). If different Wolbachia infections occur in the same host species, bidirectional CI is often induced. Bidirectional CI acts as a postzygotic isolation mechanism if parapatric host populations are infected with different Wolbachia strains. Therefore, it has been suggested that Wolbachia could promote speciation in their hosts. In this article we investigate theoretically whether Wolbachia-induced bidirectional CI selects for premating isolation and therefore reinforces genetic divergence between parapatric host populations. To achieve this we combined models for Wolbachia dynamics with a well-studied reinforcement model. This new model allows us to compare the effect of bidirectional CI on the evolution of female mating preferences with a situation in which postzygotic isolation is caused by nuclear genetic incompatibilities (NI). We distinguish between nuclear incompatibilities caused by two loci with epistatic interactions, and a single locus with incompatibility among heterozygotes in the diploid phase. Our main findings are: (1) bidirectional CI and single locus NI select for premating isolation with a higher speed and for a wider parameter range than epistatic NI; (2) under certain parameter values, runaway sexual selection leads to the increase of an introduced female preference allele and fixation of its preferred male trait allele in both populations, whereas under others it leads to divergence in the two populations in preference and trait alleles; and (3) bidirectional CI and single locus NI can stably persist up to migration rates that are two times higher than seen for epistatic NI. The latter finding is important because the speed with which mutants at the preference locus spread increases exponentially with the migration rate. In summary, our results show that bidirectional CI selects for rapid premating isolation and so generally support the view that Wolbachia can promote speciation in their hosts.     

A single-locus model of speciation

The crucial phase of speciation is argued to be the evolution of mating cross-incompatibility (prezygotic incompatibility) between the genotypes distinguishing the prospective species populations. Based on this idea, a single-locus model of speciation is presented, which is shown to be biologically plausible and may help to settle the controversy as to the biological significance of single-locus modes of speciation. The model involves three alleles, two of which characterize in homozygous state the prospective species populations and in heterozygous state their hybrids. The third allele represents a mutant which is equivalent to one of the first two alleles with the exception that it inhibits mating with carriers of the third allele. This third allele is fixed in one population and immigrates into a second population which contains the mutant inhibiting matings with members of the former population. Migration in the reverse direction does not occur. Proceeding from a widely applicable concept of fitness and mating preference it is shown that postzygotic incompatibility (hybrid or heterozygote disadvantage) alone suffices to trigger evolutionary replacement of the extant mating relations in the population receiving immigrants by any arbitrary degree of prezygotic incompatibility. This corroborates Wallace's hypothesis and emphasizes the potential biological relevance of speciation by reinforcement (parapatric speciation) at single gene loci.     

Effective population size may limit the power of laboratory experiments to demonstrate sympatric and parapatric speciation

Laboratory experiments designed to elucidate the mechanisms of sympatric and parapatric speciation may have been handicapped by too small population sizes, although this possibility has seldom been discussed. In this paper we review the published records of sympatric and parapatric speciation experiments to test the relative importance of selection intensity applied, duration of experiment and effective population size. Our results show that among these factors only effective population size has had a general effect on the generation of assortative mating. Reduced interbreeding is less likely to develop in small populations where the selection process often seems to have been opposed by inbreeding depression or loss of genetic variation. This study demonstrates that the experimental evidence frequently used as an argument against sympatric and parapatric speciation models is not as strong as previously believed.     

Speciation in peripheral populations: effects of drift load and mating systems

Speciation in peripheral populations has long been considered one of the most plausible scenarios for speciation with gene flow. In this study, however we identify two additional problems of peripatric speciation, as compared to the parapatric case, that may impede the completion of the speciation process for most parameter regions. First, with (predominantly) unidirectional gene flow, there is no selection pressure to evolve assortative mating on the continent. We discuss the implications of this for different mating schemes. Second, genetic load can build up in small populations. This can lead to extinction of the peripheral species, or generate selection pressure for lower assortative mating to avoid inbreeding. In this case, either a stable equilibrium with intermediate assortment evolves or there is cycling between phases of hybridization and phases of complete isolation.     

The evolution of conditional dispersal and reproductive isolation along environmental gradients

Dispersal modulates gene flow throughout a population's spatial range. Gene flow affects adaptation at local spatial scales, and consequently impacts the evolution of reproductive isolation. A recent theoretical investigation has demonstrated that local adaptation along an environmental gradient, facilitated by the evolution of limited dispersal, can lead to parapatric speciation even in the absence of assortative mating. This and other studies assumed unconditional dispersal, so individuals start dispersing without regard to local environmental conditions. However, many species disperse conditionally; their propensity to disperse is contingent upon environmental cues, such as the degree of local crowding or the availability of suitable mates. Here, we use an individual-based model in continuous space to investigate by numerical simulation the relationship between the evolution of threshold-based conditional dispersal and parapatric speciation driven by frequency-dependent competition along environmental gradients. We find that, as with unconditional dispersal, parapatric speciation occurs under a broad range of conditions when reproduction is asexual, and under a more restricted range of conditions when reproduction is sexual. In both the asexual and sexual cases, the evolution of conditional dispersal is strongly influenced by the slope of the environmental gradient: shallow environmental gradients result in low dispersal thresholds and high dispersal distances, while steep environmental gradients result in high dispersal thresholds and low dispersal distances. The latter, however, remain higher than under unconditional dispersal, thus undermining isolation by distance, and hindering speciation in sexual populations. Consequently, the speciation of sexual populations under conditional dispersal is triggered by a steeper gradient than under unconditional dispersal. Enhancing the disruptiveness of frequency-dependent selection, more box-shaped competition kernels dramatically lower the speciation-enabling slope of the environmental gradient.     

Reinforcement of genetic coherence in a two-locus model

Background  

In order to maintain populations as units of reproduction and thus enable anagenetic evolution, genetic factors must exist which prevent continuing reproductive separation or enhance reproductive contact. This evolutionary principle is called genetic coherence and it marks the often ignored counterpart of cladistic evolution. Possibilities of the evolution of genetic coherence are studied with the help of a two-locus model with two alleles at each locus. The locus at which viability selection takes place is also the one that controls the fusion of gametes. The second locus acts on the first by modifying the control of the fusion probabilities. It thus acts as a mating modifier whereas the first locus plays the role of the object of selection and mating. Genetic coherence is enhanced by modifications which confer higher probabilities of fusion to heterotypic gametic combinations (resulting in heterozygous zygotes) at the object locus.     

Bidirectional cytoplasmic incompatibility and the stable coexistence of two Wolbachia strains in parapatric host populations

Wolbachiaare intracellular bacteria which are very widely distributed among arthropods. In many insect species Wolbachiaare known to induce cytoplasmic mating incompatibility (CI). It has been suggested that Wolbachiacould promote speciation in their hosts if parapatric host populations are infected with two different Wolbachiastrains causing bidirectional mating incompatibilities. A necessary condition for this speciation scenario to work is that the two Wolbachiastrains can stably coexist. The following study investigates this problem analysing a mathematical model with two host populations and migration between them. We show that the stability of bidirectional CI can be fully described in terms of a critical migration rate which is defined as the highest migration below which a stable coexistence of two Wolbachiastrains is possible. For some special cases we could derive analytical solutions for the critical migration rate; for the general case estimations of the critical migration rate are given. Our main finding is that bidirectional CI can stably persist in the face of high migration and can be as high as over 15% per generation for CI levels observed in nature. These results have implications for the potential of Wolbachiato promote genetic divergence and speciation in their hosts.     

A spatially explicit individual-based model of reinforcement in hybrid zones

Abstract An individual-based model consisting of two dioecious populations in a two-dimensional environmental grid was constructed. Each population began with, and never exceeded, 1000 individuals; extinction was allowed. Genomes consisting of 30 biallelic loci for male sexual advertisement call, female mate preference, and population origin were constructed, and lineages of each individual in the starting populations were followed for 2000 generations. Type and level of hybrid disadvantage, initial population distribution, patchiness of environmental resources, and level of mate choice were varied. Persistence of bimodal hybrid zones was nonexistent at low levels of hybrid disadvantage and universal at high levels of hybrid disadvantage, with a narrow threshold in which persistence was unpredictable. Persistence occurred at lower levels of hybrid disadvantage when populations were initially parapatric rather than sympatric, and environments were patchy rather than homogeneous. Increased divergence in mating systems occurred when hybrid disadvantage was high, hybrids were infertile, populations were initially parapatric, and increased female choice was allowed. Mating system divergence was much higher in interacting populations compared with noninteracting populations, indicating that reinforcement caused most of the observed divergence. When hybrids were infertile, reinforcement contributed to speciation, because under hybrid infertility the probability of persistence at low levels of hybrid disadvantage was positively related to mate choice. The results agree with previous one-dimensional spatial models in finding that population persistence is more likely in parapatric and patchy population distributions. In addition, the results show that hybrid infertility may facilitate the process of reinforcement and speciation.     

How phenotypic matching based on neutral mating cues enables speciation in locally adapted populations

Maynard Smith's (American Naturalist, 1966, 100, 637) suggestion that in some cases a prerequisite for speciation is the existence of local ecological adaptations has not received much attention to date. Here, we test the hypothesis using a model like that of Maynard Smith but differing in the way animals disperse between niches. In previous studies, males disperse randomly between niches but females stay put in their natal niche. As a first step toward generalizing the model, we here analyze the case that equal proportions of the two sexes disperse between niches before breeding. Supporting Maynard Smith's (1966) hypothesis, we find that once local adaptations are established, a neutral mating cue at an independent locus can rapidly enable speciation in populations with a suitable mechanism for phenotype matching. We find that stable ecological polymorphisms are relatively insensitive to the strength of selection, but depend crucially on the extent of dispersal between niches, with a threshold of ~5% if population sizes in two niches are equal. At higher levels of dispersal, ecological differentiation is lost. These results contrast with those of earlier studies and shed light on why parapatric speciation is limited by the extent of gene flow. Our testable model provides a candidate explanation for the rapid speciation rates, diversity of appearance and occurrence of “species flocks” observed among some African cichlids and neotropical birds and may also have implications for the occurrence of punctuational change on phylogenies.     

Coevolution of self-fertilization and inbreeding depression. II. Symmetric overdominance in viability

We describe the evolutionary dynamics of a modifier of selfing coevolving with a locus subject to symmetric overdominance in viability under general levels of reduction in pollination success as a consequence of self-fertilization (pollen discounting). Simple models of the evolution of breeding systems that represent inbreeding depression as a constant parameter do not admit the possibility of stable mixed mating systems involving both inbreeding and random mating. Contrary to this expectation, we find that coevolution between a modifier of selfing and a single overdominant locus situated anywhere in the genome can generate evolutionarily attracting mixed mating systems. Two forms of association between the modifier locus and the viability locus promote the evolution of outcrossing. The favored heterozygous genotype at the viability locus develops positive associations with modifier alleles that enhance outcrossing and with the heterozygous genotype at the modifier locus. Associations between outcrossing and high viability evolve immediately upon the introduction of a rare modifier allele, even in the absence of linkage.     

Segregation and the evolution of sex under overdominant selection

The segregation of alleles disrupts genetic associations at overdominant loci, causing a sexual population to experience a lower mean fitness compared to an asexual population. To investigate whether circumstances promoting increased sex exist within a population with heterozygote advantage, a model is constructed that monitors the frequency of alleles at a modifier locus that changes the relative allocation to sexual and asexual reproduction. The frequency of these modifier alleles changes over time as a correlated response to the dynamics at a fitness locus under overdominant selection. Increased sex can be favored in partially sexual populations that inbreed to some extent. This surprising finding results from the fact that inbred populations have an excess of homozygous individuals, for whom sex is always favorable. The conditions promoting increased levels of sex depend on the selection pressure against the homozygotes, the extent of sex and inbreeding in the population, and the dominance of the invading modifier allele.     

Socially mediated speciation

Abstract.— We employ a simple model to show that social selection can lead to prezygotic reproductive isolation. The evolution of social discrimination causes the congealing of phenotypically similar individuals into different, spatially distinct tribes. However, tribal formation is only obtained for certain types of social behavior: altruistic and selfish acts can produce tribes, whereas spiteful and mutualistic behaviors never do. Moreover, reduced hybrid fitness at tribal borders leads to the selection of mating preferences, which then spread to the core areas of the respective tribes. Unlike models of resource competition, our model generates reproductive isolation in an ecologically homogeneous environment. We elaborate on how altruistic acts can lead to reproductive isolation, but also predict that certain types of competition can lead to the speciation effect. Our theory provides a framework for how individual-level interactions mold lineage diversification, with parapatric speciation as a possible end product.     

Assortative mating and the cost of inbreeding — A simulation approach

Assortative mating is an important factor in the process of speciation. Models of speciation frequently deal with small founder populations often with mating preferences based on ecological traits or habitat preferences. Small populations, on the other hand might suffer from inbreeding. However, few studies have explored the combined effects of assortative mating and inbreeding in such populations. Can they speciate, or are they doomed to eventually go extinct? With this simulation we show that assortative mating based on similarities increases the possibility for change in a population, as long as the population does not suffer from inbreeding depression. Inbred populations seem not to be able to cope with strong assortative mating, as this is likely to elevate the level of inbreeding, increasing the risks of inbreeding depression and as a result decreasing population mean fitness. This in turn hinders the possibility of change, and instead might drive the population to extinction.     

The genetics and ecology of sympatric speciation: A case study

Mating occurs on the larval host plant in allRhagoletis species (Diptera: Tephritidae). We show how this attribute, when coupled with certain differences in other biological traits, strongly influences the mode of speciation. In species of thesuavis species group, host shifts have never occurred during speciation, and larvae feed in the husks of any walnut species(Juglans spp.), which are highly toxic. Taxa are allopatric or parapatric and exhibit deep phylogenetic nodes suggesting relatively ancient speciation events. Traits responsible for species and mate recognition, particularly in parapatric species, are morphologically distinct and strongly sexually dimorphic. All aspects of their biology, genetics and distribution are consistent with a slow rate of allopatric speciation followed by morphological divergence in secondary contact. In contrast, speciation in thepomonella species group has always involved a shift to a new, usually unrelated, non-toxic host, and all taxa within these groups are sympatric, monophagous and morphologically indistinguishable from one another. Phylogenetic nodes are very shallow, indicating recent sympatric speciation. Sympatric divergence is promoted by genetic variation which allows a portion of the original species to shift to a new habitat or host. Evidence suggests that changes in a few key loci responsible for host selection and fitness on a new host may initiate host shifts. By exploiting different habitats, competition for resources between diverging populations is reduced or avoided. We provide evidence that in phytophagous and parasitic insects sufficient intrinsic barriers to gene flow can evolve between sister populations as they adapt to different habitats or hosts to allow each population to establish independent evolutionary lineages in sympatry.     

Towards a theory of the evolution of modifier genes

The main findings of a study of the evolution of modifier gene frequencies in models of deterministic population genetics are presented. A wide variety of random mating systems are subject to selection with modifiers operating, in different cases, on mutation rates, migration between subpopulations, and linkage between other loci. In all these instances, the modifier frequencies evolve in such a way as to maximize the mean fitness of the population at equilibrium. This is remarkable since, the modifier genes are selectively neutral in the sense that they do not affect the fitness of their individual carriers. In nonrandom mating systems, the mean fitness concept is not well-defined, and there does not appear to be such a simple principle governing the evolution of modifier frequencies. In assortative mating systems, modifiers favoring reduced assortment propensities tend to increase. In contrast, for selfing-outcrossing systems, modifiers favoring increased selfing tend to increase.     

Quantitative traits and mode of speciation in Martinique anoles

We investigate extensive quantitative trait variation (dewlap hue, colour pattern, dorsum hue, body proportions and scalation) in the Martinique anole across eight transects representing nascent parapatric ecological speciation, nascent allopatric speciation and allopatric divergence without sufficient genetic structure to suggest speciation. Quantitative trait divergence can be extremely large between adjacent sets of populations, but with one exception that this is associated with difference in habitat rather than past allopatry. Nascent ecological speciation shows the greatest level of quantitative trait divergence across all character sets including those implicated in natural, as well as sexual selection. The sole example of nascent allopatric speciation is associated with fairly strong quantitative trait divergence among most character sets, but not the set most implicated in natural (rather than sexual) selection. The role of sexual selection in ecological speciation is discussed, both in terms of female choice with assortative mating and male–male competition with condition‐dependant sexual signals.     

Mimicry and the evolution of premating isolation in Heliconius melpomene Linnaeus

Ecological divergence can cause speciation if adaptive traits have pleiotropic effects on mate choice. In Heliconius butterflies, mimetic patterns play a role in mate detection between sister species, as well as signalling to predators. Here we show that male butterflies from four recently diverged parapatric populations of Heliconius melpomene are more likely to approach and court their own colour patterns as compared with those of other races. A few exceptions, where males were more attracted to patterns other than their own, suggest that some mimetic patterns are sub-optimal in mate choice. Genotype frequencies in hybrid zones between races of H. melpomene suggest that mating is random, so reinforcement is unlikely to have played a role in intra-specific divergence. In summary, co-evolved divergence of colour pattern and mate preference occurs rapidly and is likely the first step in Heliconius speciation.     

Frequency-dependent selection and the evolution of assortative mating

A long-standing goal in evolutionary biology is to identify the conditions that promote the evolution of reproductive isolation and speciation. The factors promoting sympatric speciation have been of particular interest, both because it is notoriously difficult to prove empirically and because theoretical models have generated conflicting results, depending on the assumptions made. Here, we analyze the conditions under which selection favors the evolution of assortative mating, thereby reducing gene flow between sympatric groups, using a general model of selection, which allows fitness to be frequency dependent. Our analytical results are based on a two-locus diploid model, with one locus altering the trait under selection and the other locus controlling the strength of assortment (a "one-allele" model). Examining both equilibrium and nonequilibrium scenarios, we demonstrate that whenever heterozygotes are less fit, on average, than homozygotes at the trait locus, indirect selection for assortative mating is generated. While costs of assortative mating hinder the evolution of reproductive isolation, they do not prevent it unless they are sufficiently great. Assortative mating that arises because individuals mate within groups (formed in time or space) is most conducive to the evolution of complete assortative mating from random mating. Assortative mating based on female preferences is more restrictive, because the resulting sexual selection can lead to loss of the trait polymorphism and cause the relative fitness of heterozygotes to rise above homozygotes, eliminating the force favoring assortment. When assortative mating is already prevalent, however, sexual selection can itself cause low heterozygous fitness, promoting the evolution of complete reproductive isolation (akin to "reinforcement") regardless of the form of natural selection.     

Diazinon resistance, fluctuating asymmetry and fitness in the Australian sheep blowfly, lucilia cuprina

Genetic evidence suggests that the evolution of resistance to the insecticide diazinon in Lucilia cuprina initially produced an increase in asymmetry. At that time resistant flies were presumed to be at a selective disadvantage in the absence of diazinon. Subsequent evolution in natural populations selected modifiers to ameliorate these effects. The fitness and fluctuating asymmetry levels of resistant flies are currently similar to those of susceptibles. Previous genetic analyses have shown the fitness modifier to co-segregate with the region of chromosome III marked by the white eyes, w, locus, unlinked to the diazinon resistance locus, Rop-1, on chromosome IV. This study maps the asymmetry modifier to the same region, shows, as in the case of the fitness modifier, its effect to be dominant and presents data consistent with the fitness/asymmetry modifier being the same gene (gene complex). These results suggest changes in fluctuating asymmetry reflect changes in fitness.