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.     

Modifiers of mutation rate: a general reduction principle

A deterministic two-locus population genetic model with random mating is studied. The first locus, with two alleles, is subject to mutation and arbitrary viability selection. The second locus, with an arbitrary number of alleles, controls the mutation at the first locus. A class of viability-analogous Hardy-Weinberg equilibria is analyzed in which the selected gene and the modifier locus are in linkage equilibrium. It is shown that at these equilibria a reduction principle for the success of new mutation-modifying alleles is valid. A new allele at the modifier locus succeeds if its marginal average mutation rate is less than the mean mutation rate of the resident modifier allele evaluated at the equilibrium. Internal stability properties of these equilibria are also described.     

A single locus model of selection in permanent translocation heterozygotes

A model for evolution at a single locus in permanent translocation heterozygotes is described. It is also applicable to other permanent structural heterozygotes that possess the mating systems discussed. Recombination occurs between the locus and the chromosomes, which are of two types. The mating system includes selfing and random mating. When recombination is rare, selection will result in almost complete fixation on the single most fit genotype present in a population, regardless of the frequency of selfing. This provides a possible explanation for the ecological success of permanent translocation heterozygotes in some groups of organisms, like Oenothera. It is not necessary to postulate that their success is the result of hybrid vigor. Furthermore, a complete lack of recombination is not necessary for explaining the observed association of alleles with particular segmental arrangements. A small amount of recombination is consistent with the observation that different segmental arrangements often carry different alleles.     

EVALUATING A SIMPLE APPROXIMATION TO MODELING THE JOINT EVOLUTION OF SELF‐FERTILIZATION AND INBREEDING DEPRESSION

A comprehensive understanding of plant mating system evolution requires detailed genetic models for both the mating system and inbreeding depression, which are often intractable. A simple approximation assuming that the mating system evolves by small infrequent mutational steps has been proposed. We examine its accuracy by comparing the evolutionarily stable selfing rates it predicts to those obtained from an explicit genetic model of the selfing rate, when inbreeding depression is caused by partly recessive deleterious mutations at many loci. Both models also include pollen limitation and pollen discounting. The approximation produces reasonably accurate predictions with a low or moderate genomic mutation rate to deleterious alleles, on the order of U = 0.02–0.2. However, for high mutation rates, the predictions of the full genetic model differ substantially from those of the approximation, especially with nearly recessive lethal alleles. This occurs because when a modifier allele affecting the selfing rate is rare, homozygous modifiers are produced mainly by selfing, which enhances the opportunity for purging nearly recessive lethals and increases the marginal fitness of the allele modifying the selfing rate. Our results confirm that explicit genetic models of selfing rate and inbreeding depression are required to understand mating system evolution.     

On the Origin of Meiotic Reproduction: A Genetic Modifier Model

We study the conditions under which a rare allele that modifies the relative rates of meiotic reproduction and apomixis increases in a population in which meiotic reproduction entails selfing as well as random outcrossing. A distinct locus, at which mutation maintains alleles that are lethal in homozygous form, determines viability. We find that low viability of carriers of the lethal alleles, high rates of selfing, dominance of the introduced modifier allele, and lower rates of recombination promote the evolution of meiosis. Meiotic reproduction can evolve even in the absence of linkage between the modifier and the viability locus. The adaptive value of meiotic reproduction depends on the relative viabilities of offspring derived by meiosis and by apomixis, and on associations between the modifier and the viability locus. Meiotic reproduction, particularly under selfing, generates more diverse offspring, including those with very high and very low viability. Elimination of offspring with low viability generates positive associations between enhancers of meiotic reproduction and high viability. In addition, partial selfing generates positive associations in heterozygosity (identity disequilibrium) between the modifier and the viability locus, even in the absence of linkage. The two kinds of associations together can compensate for initial reductions in mean offspring viability under meiotic reproduction.     

Effect of selection against deleterious mutations on the decline in heterozygosity at neutral loci in closely inbreeding populations.

Transition matrices for selfing and full-sib mating were derived to investigate the effect of selection against deleterious mutations on the process of inbreeding at a linked neutral locus. Selection was allowed to act within lines only (selection type I) or equally within and between lines (type II). For selfing lines under selection type I, inbreeding is always retarded, the retardation being determined by the recombination fraction between the neutral and selected loci and the inbreeding depression from the selected locus, irrespective of the selection coefficient (s) and dominance coefficient (h) of the mutant allele. For selfing under selection type II or full-sib mating under both selection types, inbreeding is delayed by weak selection (small s and sh), due to the associative overdominance created at the neutral locus, and accelerated by strong selection, due to the elevated differential contributions between alternative alleles at the neutral locus within individuals and between lines (for selection type II). For multiple fitness loci under selection, stochastic simulations were run for populations with selfing, full-sib mating, and random mating, using empirical estimates of mutation parameters and inbreeding load in Drosophila. The simulations results are in general compatible with empirical observations.     

The evolution of mate choice in a fluctuating environment

This paper analyzes the evolutionary dynamics of a locus controlling the degree of female mating preference in a temporally fluctuating environment. Preference for mating with males with respect to their genotypes at a locus that is subject to temporally varying natural selection pressure is considered first. With weak selection and free recombination between the choice locus and the selected locus, preference for mating with heterozygotes appears to be favored. With strong selection, preference for homozygous mates may be favored. In each case, choice alleles may increase from very low initial frequencies to near fixation, in contrast to previous models of mate choice in varying environments. Linkages between the two loci has complex effects on the strength and direction of selection for mate choice. Preference for mating with males with the currently fitter genotypes at the locus under natural selection is also modelled. Provided that the environmental period is not too short, a rare allele conferring such preference may be favored and spread to fixation. Strong natural selection, tight linkage and a short environmental period may produce polymorphism for the level of mate choice.     

INBREEDING DEPRESSION WITH HETEROZYGOTE ADVANTAGE AND ITS EFFECT ON SELECTION FOR MODIFIERS CHANGING THE OUTCROSSING RATE

The equilibrium level of inbreeding depression in populations with different selfing rates is studied for models with symmetrical or asymmetrical heterozygous advantage at several loci with partial linkage. As for the case of a single locus, the inbreeding depression caused by loci with heterozygous advantage can be higher for partially selfing populations than for complete outcrossing. The spread of modifier alleles at another locus that affects the selfing rate is studied. The stability of outcrossing populations to invasion by alleles that give increased selfing is found to depend on levels of inbreeding depression being greater than one-half, in accordance with earlier models that assumed a fixed level of inbreeding depression. However, in partially selfing populations the spread of such alleles can be checked by smaller levels of inbreeding depression than one-half, so that they do not always spread to fixation. This is interpreted as being due to associations between the genotypes at the modifier locus and the selected loci, together with increasing inbreeding depression as selfing increases, and does not occur if the inbreeding depression is due to mutation-selection balance.     

Hitchhiking: A Comparison of Linkage and Partial Selfing

Genetic hitchhiking occurs when alleles at unselected loci are changed in frequency because of an association with alleles at a selected locus. This association may be mediated either by linkage or partial selfing (inbreeding) and can affect the gene frequency and gametic disequilibrium at the neutral loci. Hitchhiking from partial selfing (unlinked loci) occurs more quickly than linkage hitchhiking and generally has a greater effect. In addition, partial-selfing hitchhiking can cause increases or changes in sign in gametic disequilibrium between neutral loci. The effects of the two types of hitchhiking with different levels of dominance, zygotic frequencies and number of selected loci are also examined. The general conditions for linkage and partial-selfing hitchhiking are outlined and the implications of hitchhiking are discussed for marker or electrophoretic loci.     

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.     

Clinal genetic variation and the 'rare allele phenomenon' in random mating populations of <Emphasis Type="Italic">Urophora cardui</Emphasis> (Diptera: Tephritidae)

In the present study we investigate a contact zone between two population groups of the tephritid fly Urophora cardui. We investigate scenarios that may have produced the genetic differentiation of the two groups, and we describe the 'rare allele phenomenon' from the contact zone. The rare allele phenomenon refers to alleles that are found at high frequency in contact zones but are rare or lacking outside the contact zone. The phenomenon is often observed in hybrid zones between subspecies of limited reproductive compatibility, but seldom in populations with random mating. Clinal genetic variation was observed at three loci in the contact zone. Three alleles at the locus Aat showed steep clines, between 20--70 km wide. A rare Aat-A allele occurred at high frequency in the centre of the contact zone. Two further loci, Hk and Pgd, showed less steep clinal genetic variation, the transition being in and slightly south of the centre of the Aat cline. Populations showed Hardy--Weinberg proportions and there was no evidence for linkage dis-equlibrium. These findings suggest random mating and gradual introgression between the population systems, which may originate from at least two range expansions. Aat's steep clines and rare allele may indicate selection on Aat alleles, although we presently can not quantify any agents. Because U. cardui experiences random mating in the contact zone with no apparent 'hybrid' incompatibility, mating experiments offer the possibility for future enquiries about the genetic basis of the rare allele phenomenon.     

Host mating system and the prevalence of disease in a plant population

A modified susceptible-infected-recovered (SIR) host-pathogen model is used to determine the influence of plant mating system on the outcome of a host-pathogen interaction. Unlike previous models describing how interactions between mating system and pathogen infection affect individual fitness, this model considers the potential consequences of varying mating systems on the prevalence of resistance alleles and disease within the population. If a single allele for disease resistance is sufficient to confer complete resistance in an individual and if both homozygote and heterozygote resistant individuals have the same mean birth and death rates, then, for any parameter set, the selfing rate does not affect the proportions of resistant, susceptible or infected individuals at equilibrium. If homozygote and heterozygote individual birth rates differ, however, the mating system can make a difference in these proportions. In that case, depending on other parameters, increased selfing can either increase or decrease the rate of infection in the population. Results from this model also predict higher frequencies of resistance alleles in predominantly selfing compared to predominantly outcrossing populations for most model conditions. In populations that have higher selfing rates, the resistance alleles are concentrated in homozygotes, whereas in more outcrossing populations, there are more resistant heterozygotes.     

The Hill-Robertson effect and the evolution of recombination

In finite populations, genetic drift generates interference between selected loci, causing advantageous alleles to be found more often on different chromosomes than on the same chromosome, which reduces the rate of adaptation. This "Hill-Robertson effect" generates indirect selection to increase recombination rates. We present a new method to quantify the strength of this selection. Our model represents a new beneficial allele (A) entering a population as a single copy, while another beneficial allele (B) is sweeping at another locus. A third locus affects the recombination rate between selected loci. Using a branching process model, we calculate the probability distribution of the number of copies of A on the different genetic backgrounds, after it is established but while it is still rare. Then, we use a deterministic model to express the change in frequency of the recombination modifier, due to hitchhiking, as A goes to fixation. We show that this method can give good estimates of selection for recombination. Moreover, it shows that recombination is selected through two different effects: it increases the fixation probability of new alleles, and it accelerates selective sweeps. The relative importance of these two effects depends on the relative times of occurrence of the beneficial alleles.     

Consequences of genetic linkage for the maintenance of sexually antagonistic polymorphism in hermaphrodites

When selection differs between males and females, pleiotropic effects among genes expressed by both sexes can result in sexually antagonistic selection (SA), where beneficial alleles for one sex are deleterious for the other. For hermaphrodites, alleles with opposing fitness effects through each sex function represent analogous genetic constraints on fitness. Recent theory based on single‐locus models predicts that the maintenance of SA genetic variation should be greatly reduced in partially selfing populations. However, selfing also reduces the effective rate of recombination, which should facilitate selection on linked allelic combinations and expand opportunities for balancing selection in a multilocus context. Here, I develop a two‐locus model of SA selection for simultaneous hermaphrodites, and explore the joint influence of linkage, self‐fertilization, and dominance on the maintainance of SA polymorphism. I find that the effective reduction in recombination caused by selfing significantly expands the parameter space where SA polymorphism can be maintained relative to single‐locus models. In particular, linkage facilitates the invasion of male‐beneficial alleles, partially compensating for the “female‐bias” in the net direction of selection created by selfing. I discuss the implications of accounting for linkage among SA loci for the maintenance of SA genetic variation and mixed mating systems in hermaphrodites.     

On the evolution of fluctuating segregation distortion

Previous mathematical models of the genetic control by one locus of the segregation at another have all concluded that alleles causing departures from Mendelian segregation should succeed. In this study the segregation ratios induced at the major locus by the modifier locus fluctuate cyclically. It is shown that if initially there is Mendelian segregation and if the rare modifying allele induces symmetric fluctuation about the Mendelian ratios it cannot succeed. It is further proven that if initially there are symmetric fluctuations about Mendelian segregation then an allele reducing the amplitude of the fluctuation will succeed.     

Linkage analysis of "necessary" disease loci versus "susceptibility" loci.

The association of some diseases with specific alleles of certain genetic markers has been difficult to explain. Several explanations have been proposed for the phenomenon of association, e.g. the existence of multiple, interacting genes (epistasis) or a disease locus in linkage disequilibrium with the marker locus. One might suppose that when marker data from families with associated diseases are analyzed for linkage, the existence of the association would assure that linkage will be found, and found at a tight recombination fraction. In fact, however, linkage analyses of some diseases associated with HLA, as well as diseases associated with alleles at other loci located throughout the genome, show significant evidence against linkage, and others show loose linkage, to the puzzlement of many researchers. In part, the puzzlement arises because linkage analysis is ideal for looking for loci that are necessary, even if not sufficient, for disease expression but may be much less useful for finding loci that are neither necessary nor sufficient for disease expression (so-called susceptibility loci). This work explores what happens when one looks for linkage to susceptibility loci. A susceptibility locus in this case means that the allele increases risk but is neither necessary nor sufficient for disease expression. It might be either an allele at the marker locus itself that is increasing susceptibility or an allele at a locus in linkage disequilibrium with the marker. This work uses computer simulation to examine how linkage analyses behave when confronted with data from such a model.(ABSTRACT TRUNCATED AT 250 WORDS)     

On the evolution of genetic incompatibility systems. IV. Modification of response to an existing antigen polymorphism under partial selfing

A 2-locus model of the evolution of self-incompatibility in a population practicing partial selfing is presented. An allele is introduced at a modifier locus which influences the strength of the rejection reaction expressed by the style in response to antigens recognized in pollen. Two causes of inbreeding depression are investigated. First, offspring viability depends solely on the source (self or non-self) of the fertilizing pollen. Second, offspring viability declines with the expression of recessive deleterious alleles, segregating at a third (disease) locus, which exhibit an imperfect association with antigen alleles. Evolutionary changes occurring at the disease locus are not considered in this study. The condition under which a modifier allele that intensifies the incompatibility reaction increases when rare depends upon the number of antigens, the frequency of recessive deleterious alleles at the disease locus, and the level of association between the antigen locus and the disease locus. It is the improvement of viability among offspring derived by outcrossing, rather than the prevention of self-fertilization, that may represent the primary evolutionary function of genetic incompatibility systems.     

On the Theory of Partially Inbreeding Finite Populations. I. Partial Selfing

Some stochastic theory is developed for monoecious populations of size N in which there are probabilities beta and 1 - beta of reproduction by selfing and by random mating. It is assumed that beta much greater than N-1. Expressions are derived for the inbreeding coefficient of one random individual and the coefficient of kinship of two random separate individuals at time t. The mean and between-lines variance of the fraction of copies of a locus that are identical in two random separate individuals in an equilibrium population are obtained under the assumption that there is an infinite number of possible alleles. It is found that the theory for random mating populations holds if the effective population number is Ne = N'/(1 + FIS), where FIS is the inbreeding coefficient at equilibrium when N is infinite and N' is the reciprocal of the probability that two gametes contributing to random separate adults come from the same parent. When there is a binomial distribution of successful gametes emanating from each adult, N' = N. An approximation to the probability that an allele A survives if it is originally present in one AA heterozygote is found to be 2(N'/N)(FISS1 + (1 - FIS)S2), where S1 and S2 are the selective advantages of AA and AA in comparison with AA. In the last section it is shown that if there is partial full sib mating and binomial offspring distributions Ne = N/(1 + 3FIS).     

Effects of population size and metapopulation dynamics on a mating-system polymorphism

The evolutionary dynamics of neutral alleles under the Wright-Fisher model are well understood. Similarly, the effect of population turnover on neutral genetic diversity in a metapopulation has attracted recent attention in theoretical studies. Here we present the results of computer simulations of a simple model that considers the effects of finite population size and metapopulation dynamics on a mating-system polymorphism involving selfing and outcrossing morphs. The details of the model are based on empirical data from dimorphic populations of the annual plant Eichhornia paniculata, but the results are also of relevance to species with density-dependent selfing rates in general. In our model, the prior selfing rate is determined by two alleles segregating at a single diploid locus. After prior selfing occurs, some remaining ovules are selfed through competing self-fertilisation in finite populations as a result of random mating among gametes. Fitness differences between the mating-system morphs were determined by inbreeding depression and pollen discounting in a context-dependent manner. Simulation results showed evidence of frequency dependence in the action of pollen discounting and inbreeding depression in finite populations. In particular, as a result of selfing in outcrossers through random mating among gametes, selfers experienced a "fixation bias" through drift, even when the mating-system locus was selectively neutral. In a metapopulation, high colony turnover generally favoured the fixation of the outcrossing morph, because inbreeding depression reduced opportunities for colony establishment by selfers through seed dispersal. Our results thus demonstrate that population size and metapopulation processes can lead to evolutionary dynamics involving pollen and seed dispersal that are not predicted for large populations with stable demography.     

Intratetrad mating and the evolution of linkage relationships

Mating among the immediate products of meiosis (intratetrad mating) is a common feature of many organisms with parthenogenesis or with mating-type determination in the haploid phase. Using a three-locus deterministic model we show that intratetrad mating, unlike other systems of mating, allows sheltering of deleterious recessive alleles even if there is only partial linkage between a mating locus and a load locus. Moreover, modifiers that reduce recombination between the load and mating-type locus will spread to fixation, even when there is no linkage disequilibrium between these loci in the population as a whole. This seeming contradiction to classical expectation is because partial linkage generates linkage disequilibrium among segregating loci within a tetrad, which then acts as the "mating unit."