PLOIDY AND EVOLUTION BY SEXUAL SELECTION: A COMPARISON OF HAPLOID AND DIPLOID FEMALE CHOICE MODELS NEAR FIXATION EQUILIBRIA

We compare the stability properties of haploid and diploid models of Fisherian sexual selection (with male contribution limited to sperm) by examining both models at equilibria for which a male trait is fixed or absent. Haploid and diploid two locus diallelic models share the property that the stability of such fixation equilibria is determined by the relationship between the harmonic mean of relative preference values for the common male trait, weighted by the frequency of the preferences, and the relative viability associated with the common male trait. When diploid females with heterozygotic-based preferences express preference strengths intermediate between homozygote-based preferences, then boundary equilibria of haploid and diploid models share many stability properties. However, even with intermediate heterozygote preferences, haploid and diploid models do differ: (1) for a particular frequency of the preference allele, both fixation boundaries can be stable for the diploid model, and (2) with over- or underdominance at the preference locus (a possibility precluded in the haploid model), a fixation boundary in the diploid model may show two switches in its stability state for increasing frequencies of one of the preference alleles. These differences are due not just to the impossibility of dominance in haploid models, but also to the larger number of diploid genotypes.     

DYNAMICS OF SEXUAL SELECTION IN DIPLOID POPULATIONS

We describe results for a diploid, two-locus model for the evolution of a female mating preference directed at an attractive male trait that is subject to viability and/or fertility selection. Using computer simulation, we studied a large, random sample of parameter values, assuming additivity of alleles at the preference locus and partial dominance at the trait locus. Simulation results were classifiable into nine types of parameter sets, each differing in equilibria, evolutionary trajectories, and rates of evolution. For many parameters, evolutionary trajectories converged on curves within the allelic frequency plane and subsequently evolved along the curves toward fixation. Neutrally stable curves of equilibria did not occur in Fisherian models that assume only viability and sexual selection unless there is complete dominance at the trait locus. The Fisherian models also exhibited oscillation of allelic frequencies and unique polymorphic equilibria. “Sexy son” models in which attractive males had reduced fertility were much less likely to lead to increase in traits and preferences than were the Fisherian models. However, if less fertile males had increased viability, trait polymorphisms and fixation of rare “sexy” alleles occurred. In general, the behavior of the diploid model was much more complex than that of analogous haploid or polygenic models.     

ON EVOLUTION UNDER SEXUAL AND VIABILITY SELECTION: A TWO-LOCUS DIPLOID MODEL

A two-locus diploid model of sexual selection is presented in which the two loci govern, respectively, a trait limited in expression in one sex (generally male) and the mating preferences of the other sex (generally female). The viability of a male depends on its genotype at the trait locus. In contrast, all females are equally viable and all individuals are equally fertile with respect to the two loci. Near fixation at both loci, evolution at the mating locus is neutral and hence a new mating preference allele will increase only through random genetic drift or through a correlated response to the increase of a new advantageous trait allele. If, however, a polymorphism is already maintained at the trait locus through overdominance in fitness then the increase of a rare preference allele depends only on the recombination rate between the loci and not on the new preference scheme.     

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.     

Natural Selection and Y-Linked Polymorphism

Several population genetic models allowing natural selection to act on Y-linked polymorphism are examined. The first incorporates pleiotropic effects of a Y-linked locus, such that viability, segregation distortion, fecundity and sexual selection can all be determined by one locus. It is shown that no set of selection parameters can maintain a stable Y-linked polymorphism. Interaction with the X chromosome is allowed in the next model, with viabilities determined by both X- and Y-linked factors. This model allows four fixation equilibria, two equilibria with X polymorphism and a unique point with both X- and Y-linked polymorphism. Stability analysis shows that the complete polymorphism is never stable. When X- and Y-linked loci influence meiotic drive, it is possible to have all fixation equilibria simultaneously unstable, and yet there is no stable interior equilibrium. Only when viability and meiotic drive are jointly affected by both X- and Y-linked genes can a Y-linked polymorphism be maintained. Unusual dynamics, including stable limit cycles, are generated by this model. Numerical simulations show that only a very small portion of the parameter space admits Y polymorphism, a result that is relevant to the interpretation of levels of Y-DNA sequence variation in natural populations.     

General analysis of frequency-dependent sexual selection at a multi-allelic locus

A general model is analyzed in which arbitrarily frequency-dependent selection acts on one sex of a diploid population with several alleles at one locus, as a result of viability or mating-success differences. The existence of boundary and polymorphic equilibria is examined, and conditions for local stability, internal and external, are obtained. The status of Hardy-Weinberg approximations in studying stability and approach to equilibria is also considered. The general principles are then applied to two specific models: one where genotypes fall into two phenotypic classes; and one with a hierarchy of dominance where viability and sexual selection are opposed. In the latter case it is found that, of all the equilibria present, there is one and only one which could possibly be stable: the existence of a unique globally stable equilibrium might then be inferred.     

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.     

Revisiting santa rosalia to unfold a degeneracy of classic models of speciation

Abstract Many classic models of speciation incorporate assortative mating based on mating groups, such as plants with different flowering times, and they investigate whether an ecological trait under disruptive natural selection becomes genetically associated with the selectively neutral mating trait. It is well known that this genetic association is potently destroyed by recombination. In this note, we point out a more fundamental difficulty: if a "knife-edge" symmetry assumption of previous models is violated, then the mating trait is no longer neutral and sexual selection eliminates the polymorphism in the mating locus. This result strengthens the growing consensus that magic traits are the more likely route to nonallopatric speciation. We expand the model assuming also ecological selection on the mating trait and investigate the conditions for natural selection to overcome sexual selection and maintain mating polymorphism; we find that the combination of natural and sexual selection can cause also bistability of allele frequencies.     

Selective sweeps in multilocus models of quantitative traits

We study the trajectory of an allele that affects a polygenic trait selected toward a phenotypic optimum. Furthermore, conditioning on this trajectory we analyze the effect of the selected mutation on linked neutral variation. We examine the well-characterized two-locus two-allele model but we also provide results for diallelic models with up to eight loci. First, when the optimum phenotype is that of the double heterozygote in a two-locus model, and there is no dominance or epistasis of effects on the trait, the trajectories of selected mutations rarely reach fixation; instead, a polymorphic equilibrium at both loci is approached. Whether a polymorphic equilibrium is reached (rather than fixation at both loci) depends on the intensity of selection and the relative distances to the optimum of the homozygotes at each locus. Furthermore, if both loci have similar effects on the trait, fixation of an allele at a given locus is less likely when it starts at low frequency and the other locus is polymorphic (with alleles at intermediate frequencies). Weaker selection increases the probability of fixation of the studied allele, as the polymorphic equilibrium is less stable in this case. When we do not require the double heterozygote to be at the optimum we find that the polymorphic equilibrium is more difficult to reach, and fixation becomes more likely. Second, increasing the number of loci decreases the probability of fixation, because adaptation to the optimum is possible by various combinations of alleles. Summaries of the genealogy (height, total length, and imbalance) and of sequence polymorphism (number of polymorphisms, frequency spectrum, and haplotype structure) next to a selected locus depend on the frequency that the selected mutation approaches at equilibrium. We conclude that multilocus response to selection may in some cases prevent selective sweeps from being completed, as described in previous studies, but that conditions causing this to happen strongly depend on the genetic architecture of the trait, and that fixation of selected mutations is likely in many instances.     

On the logic of Fisherian sexual selection*

In Fisher's model of sexual selection, a female preference for a male trait spreads together with the trait because their genetic bases become correlated. This can be interpreted as a “greenbeard” system: a preference gene, by inducing a female to mate with a trait-bearing male, favors itself because the male is disproportionately likely also to carry the preference gene. Here, we use this logic to argue that Fisherian sexual selection in diploids proceeds via two channels: (i) trait-bearing males are disproportionately the product of matings between preference-bearing mothers and trait-bearing fathers, and thus trait and preference genes are correlated “in trans”; (ii) trait and preference genes come into gametic phase disequilibrium, and thus are correlated “in cis.” Gametic phase disequilibrium is generated by three distinct mechanisms that we identify. The trans channel does not operate when sexual selection is restricted to the haploid phase, and therefore represents a fundamental difference between haploid and diploid models of sexual selection. We show that the cis and trans channels contribute equally to the spread of the preference when recombination between the preference and trait loci is free, but that the trans channel is substantially more important when linkage is tight.     

Conditions for the existence of stationary densities for some two-dimensional diffusion processes with applications in population biology

Conditions are derived that we conjecture are necessary and sufficient for the existence of stationary densities for a class of two-dimensional diffusion processes. The derivation of the conditions rests on the assumption that a two-dimensional stationary density (which can be viewed as a stable “internal equilibrium”) exists if and only if all “boundary equilibria” are unstable in the sense that small perturbations lead to moving away from the boundaries with high probability. For the models considered, the boundary equilibria are one-dimensional stationary densities and equilibrium points. To demonstrate the usefulness of the conditions, three random environment models are analyzed: a three-allele selection model, a two-species competition model, and a two-locus selection model. Several of the results obtained have been verified by alternate methods.     

Sexually antagonistic co‐evolution: a model and an empirical test

Models reveal that sexually antagonistic co‐evolution exaggerates female resistance and male persistence traits. Here we adapt an established model by including directional sexual selection acting against persistence. We find similar equilibria to previous models showing that sexually antagonistic co‐evolution can be limited by counteracting sexual, as well as, natural selection. We tested the model using empirical data for the seaweed fly, Coelopa ursina, in which body size acts as a persistence and a resistance trait. Our model can generate continuous co‐evolutionary cycles and stable equilibria, however, all simulations using empirically derived parameter estimates reach stable equilibria. Thus, stable equilibria might be more common in nature than continuous co‐evolutionary cycles, suggesting that sexual conflict is unlikely to promote speciation. The model predicts male biased sexual size dimorphism for C. ursina, comparable with empirically observed values. Male persistence is shown to be more sensitive than female resistance to changes in model parameters.     

Models of sexual selection on a quantitative genetic trait when preference is acquired by sexual imprinting

The evolution of a quantitative genetic trait under stabilizing viability selection and sexual selection is modeled for a polygynous species in which female mating preferences are acquired by sexual imprinting on the parents and by exposure to the surviving population at large. Stabilizing viability selection acts equally on both sexes in the case of a sexually monomorphic trait and on males only in the case of a dimorphic trait. A genetically fixed sensory or perceptual bias defines the origin of the scale on which the trait is measured, and the possibility is incorporated that female preferences may deviate asymmetrically from the familiar-either toward or away from this origin. When viability selection is strong relative to sexual selection, the models predict that the mean trait value will evolve to the viability optimum. With intermediate ratios of the strength of viability to sexual selection, a stable equilibrium can occur on either side of this viability optimum, depending on the direction of asymmetry in female preferences. When viability selection is relatively weak and certain other conditions are also satisfied, runaway selection is predicted.     

Identifying signatures of sexual selection using genomewide selection components analysis

Sexual selection must affect the genome for it to have an evolutionary impact, yet signatures of selection remain elusive. Here we use an individual‐based model to investigate the utility of genome‐wide selection components analysis, which compares allele frequencies of individuals at different life history stages within a single population to detect selection without requiring a priori knowledge of traits under selection. We modeled a diploid, sexually reproducing population and introduced strong mate choice on a quantitative trait to simulate sexual selection. Genome‐wide allele frequencies in adults and offspring were compared using weighted F_ST values. The average number of outlier peaks (i.e., those with significantly large F_ST values) with a quantitative trait locus in close proximity (“real” peaks) represented correct diagnoses of loci under selection, whereas peaks above the F_ST significance threshold without a quantitative trait locus reflected spurious peaks. We found that, even with moderate sample sizes, signatures of strong sexual selection were detectable, but larger sample sizes improved detection rates. The model was better able to detect selection with more neutral markers, and when quantitative trait loci and neutral markers were distributed across multiple chromosomes. Although environmental variation decreased detection rates, the identification of real peaks nevertheless remained feasible. We also found that detection rates can be improved by sampling multiple populations experiencing similar selection regimes. In short, genome‐wide selection components analysis is a challenging but feasible approach for the identification of regions of the genome under selection.     

Stochasticity in Sexual Selection Enables Divergence: Implications for Moth Pheromone Evolution

Sexual selection has long been hypothesized to lead to allopatric speciation, and one possible mechanism for this is that its interaction with stochasticity, which perturbs the trait and preference equilibria, can result in different traits being preferred in different populations. Here we specifically examine the role that stochastic changes in sexual selection strength plays in the shift of predominance between pairs of preferences and traits within a single population. We first create a single-locus null model of stochasticity during frequency dependent selection and then compare it to a two-locus population genetic model with stochastic strengths of female preferences for male traits. We find some interesting differences between the two models, primarily that in the two-locus sexual selection model shifts between preference and trait regimes occur more often with both weak and strong preferences, compared to intermediate preference strengths. We discuss the implications of our results for the evolution of pheromone blends and male responses during speciation in moths, a case that seems to match the assumptions of our model.     

EVOLUTION OF DOMINANCE UNDER FREQUENCY-DEPENDENT INTRASPECIFIC COMPETITION IN AN ASSORTATIVELY MATING POPULATION

We study the evolution of higher levels of dominance as a response to negative frequency-dependent selection. In contrast to previous studies, we focus on the effect of assortative mating on the evolution of dominance under frequency-dependent intraspecific competition. We analyze a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under a mixture of frequency-independent stabilizing selection, density-dependent selection, and frequency-dependent selection caused by intraspecific competition for a continuum of resources. The second (modifier) locus determines the degree of dominance at the trait level. Additionally, the population mates assortatively with respect to similarities in the ecological trait. Our analysis shows that the parameter region in which dominance can be established decreases if small levels of assortment are introduced. In addition, the degree of dominance that can be established also decreases. In contrast, if assortment is intermediate, sexual selection for extreme types can be established, which leads to evolution of higher levels of dominance than under random mating. For modifiers with large effects, intermediate levels of assortative mating are most favorable for the evolution of dominance. For large modifiers, the speed of fixation can even be higher for intermediate levels of assortative mating than for random mating.     

Maintenance of Genetic Variability under Strong Stabilizing Selection: A Two-Locus Model

We study a two locus model with additive contributions to the phenotype to explore the relationship between stabilizing selection and recombination. We show that if the double heterozygote has the optimum phenotype and the contributions of the loci to the trait are different, then any symmetric stabilizing selection fitness function can maintain genetic variability provided selection is sufficiently strong relative to linkage. We present results of a detailed analysis of the quadratic fitness function which show that selection need not be extremely strong relative to recombination for the polymorphic equilibria to be stable. At these polymorphic equilibria the mean value of the trait, in general, is not equal to the optimum phenotype, there exists a large level of negative linkage disequilibrium which ``hides' additive genetic variance, and different equilibria can be stable simultaneously. We analyze dependence of different characteristics of these equilibria on the location of optimum phenotype, on the difference in allelic effect, and on the strength of selection relative to recombination. Our overall result that stabilizing selection does not necessarily eliminate genetic variability is compatible with some experimental results where the lines subject to strong stabilizing selection did not have significant reductions in genetic variability.     

The effects of sexual selection on trait divergence in a peripheral population with gene flow

The unique aspects of speciation and divergence in peripheral populations have long sparked much research. Unidirectional migration, received by some peripheral populations, can hinder the evolution of distinct differences from their founding populations. Here, we explore the effects that sexual selection, long hypothesized to drive the divergence of distinct traits used in mate choice, can play in the evolution of such traits in a partially isolated peripheral population. Using population genetic continent‐island models, we show that with phenotype matching, sexual selection increases the frequency of an island‐specific mating trait only when female preferences are of intermediate strength. We identify regions of preference strength for which sexual selection can instead cause an island‐specific trait to be lost, even when it would have otherwise been maintained at migration‐selection balance. When there are instead separate preference and trait loci, we find that sexual selection can lead to low trait frequencies or trait loss when female preferences are weak to intermediate, but that sexual selection can increase trait frequencies when preferences are strong. We also show that novel preference strengths almost universally cannot increase, under either mating mechanism, precluding the evolution of premating isolation in peripheral populations at the early stages of species divergence.     

Evolution of sexual preferences in quantitative characters

An analysis of equilibria and dynamics of the means, variances, and covariances of female mating preference for a quantitative male secondary sexual character following a Gaussian model is presented. For many combinations of viability and sexual selection parameters the evolving Gaussian distribution of phenotypes can diverge. The results on the cases of convergence and their limiting forms suggest some reinterpretations of Fisher's "runaway" process of sexual selection. One possibility is to interpret Fisher's postulated "initial advantage not due to female preference" as a shift in viability selection where runaway evolution occurs if the mean preferred trait evolves beyond its new viability optimum (due to sexual selection). This definition is contrasted with situations in which the new viability optimum is undershot. The quantitative and qualitative conclusions differ from models that approximate genetic covariance evolution involving a constant covariance.     

Analysis of partial assortative mating and sexual selection models for a polygamous species involving a trait based on levels of heterozygosity

The consequences of preferential mating in the presence of partial assortative and sexual selection mechanisms are ascertained for a two-allele one-locus trait involving two phenotype classes C₁ = {all homozygotes} and C₂ = {heterozygotes}. Relevant biological cases may include Burley (1977, Proc. Nat. Acad. Sci. USA74, 3476–3479), Wilbur et al. (1978, Evolution32, 264–270), and Singh and Zouros (1978, Evolution32, 342–353). When the preference rate for the heterozygote exceeds that for homozygotes, it is established that the unique stable state is the central Hardy-Weinberg equilibrium. The rate of approach is faster with sexual selection than for the corresponding model of assortative mating. When the preference rates favor the homozygotes then in this symmetric model of sexual selection two asymmetric Hardy-Weinberg polymorphisms can evolve, and which succeeds depends on initial conditions. The models are also analyzed with natural selection acting on phenotypes superimposed on assortative mating. In this case we can have up to three coexisting stable states involving both fixation alternatives and a central polymorphism. The corresponding model with sexual selection maintains either the central equilibrium as in assortative mating or two asymmetric polymorphic equilibria.