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.     

Joint effect of selection, linkage and partial inbreeding on additive genetic variance in an infinite population

Under the inifinitesimal model of gene effects, selection reduces the additive genetic variance by inducing negative linkage disequilibrium among selected genes. If the selected genes are linked, the decay of linkage disequilibrium is delayed, and the reduction of additive genetic variance is enhanced. Inbreeding in an infinite population also alters the additive genetic variance through the generation of positive association among genes within a locus. In the present study, the joint effect of selection, linkage and partial inbreeding (partial selfing or partial full-sib mating) on the additive genetic variance was modeled. The recurrence relations of the additive genetic variance between successive generations and the prediction equation of the asymptotic additive genetic variance were derived. Numerical computation showed that although partially inbred populations initially maintain larger genetic variances, the accumulated effect of selection overrides the effect of inbreeding. Stochastic simulation was carried out to check the precision of prediction, showing that the obtained equations give a satisfactory prediction during initial generations. However, the predicted values always overestimate the simulated values, especially in later generations. Based on these results, possible extensions and perspectives of the assumed model were discussed.     

Effects of Partial Inbreeding on Fixation Rates and Variation of Mutant Genes

Diffusion methods were used to investigate the fixation probability, average time until fixation and extinction, and cumulative heterozygosity and genetic variance for single mutant genes in finite populations with partial inbreeding. The critical parameters in the approximation are the coefficient of inbreeding due to nonrandom mating (F) and the effective population size (Ne), which also depends on F and the variance of family size. For large Ns, the fixation probability (u) is u = 2(Ne/N)s (F + h - Fh), where N is the population census, s is the coefficient of selection of the mutant homozygote and h is the coefficient of dominance. For Poisson family size (independent Poisson distributions of selfed and nonselfed offspring with partial selfing, and independent Poisson distributions of male and female numbers with partial sib mating), Ne = N/(1 + F), and the time until fixation is approximately equal to Ne/N times the time to fixation with random mating, but this relation does not hold, however, for other distributions of family size. The cumulative nonadditive variance until fixation or loss for dominant genes is reduced with increasing F while for recessive genes it is increased with intermediate values of F. The average time until extinction of deleterious mutations is reduced by increasing F. This reduction, when expressed as a proportion, is approximately independent of the initial gene frequency as well as the selective disadvantage if this is large.     

Effect of linkage on the control of inbreeding in selection programmes

Selection and mating methods for controlling inbreeding in selection programmes are based on relationships obtained from pedigrees. The efficiency of these methods has always been tested by studies using genetic models of independent loci. However, under linkage the rate of inbreeding obtained from pedigrees can be different from the probability of identity by descent of genes. We simulated a quantitative trait under artificial selection controlled by a large number of genes spread on genome regions of different sizes. A method to control inbreeding based on minimising the average coancestry of selected individuals with a restriction in the loss of selection response, and a mating procedure to control inbreeding were applied. These methods, that use coancestry relationships, were not effective in controlling inbreeding when the genome sizes were smaller than five morgans or so. However, for larger genome sizes the methods were sufficiently efficient. For very tight linkage, methods that utilise molecular information from markers should be used. We finally discuss the effects of the selection of individual major genes on the neutral variability of adjacent genome regions.     

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.     

Prediction of rates of inbreeding in selected populations

A method is presented for the prediction of rate of inbreeding for populations with discrete generations. The matrix of Wright's numerator relationships is partitioned into 'contribution' matrices which describe the contribution of the Mendelian sampling of genes of ancestors in a given generation to the relationship between individuals in later generations. These contributions stabilize with time and the value to which they stabilize is shown to be related to the asymptotic rate of inbreeding and therefore also the effective population size, Ne approximately 2N/(mu 2r + sigma 2r), where N is the number of individuals per generation and mu r and sigma 2r are the mean and variance of long-term relationships or long-term contributions. These stabilized values are then predicted using a recursive equation via the concept of selective advantage for populations with hierarchical mating structures undergoing mass selection. Account is taken of the change in genetic parameters as a consequence of selection and also the increasing 'competitiveness' of contemporaries as selection proceeds. Examples are given and predicted rates of inbreeding are compared to those calculated in simulations. For populations of 20 males and 20, 40, 100 or 200 females the rate of inbreeding was found to increase by as much as 75% over the rate of inbreeding in an unselected population depending on mating ratio, selection intensity and heritability of the selected trait. The prediction presented here estimated the rate of inbreeding usually within 5% of that calculated from simulation.     

The effect of linkage and population size on inbreeding depression due to mutational load.

Using a stochastic model of a finite population in which there is mutation to partially recessive detrimental alleles at many loci, we study the effects of population size and linkage between the loci on the population mean fitness and inbreeding depression values. Although linkage between the selected loci decreases the amount of inbreeding depression, neither population size nor recombination rate have strong effects on these quantities, unless extremely small values are assumed. We also investigate how partial linkage between the loci that determine fitness affects the invasion of populations by alleles at a modifier locus that controls the selfing rate. In most of the cases studied, the direction of selection on modifiers was consistent with that found in our previous deterministic calculations. However, there was some evidence that linkage between the modifier locus and the selected loci makes outcrossing less likely to evolve; more losses of alleles promoting outcrossing occurred in runs with linkage than in runs with free recombination. We also studied the fate of neutral alleles introduced into populations carrying detrimental mutations. The times to loss of neutral alleles introduced at low frequency were shorter than those predicted for alleles in the absence of selected loci, taking into account the reduction of the effective population size due to inbreeding. Previous studies have been confined to outbreeding populations, and to alleles at frequencies close to one-half, and have found an effect in the opposite direction. It therefore appears that associations between neutral and selected loci may produce effects that differ according to the initial frequencies of the neutral alleles.     

THE EVOLUTION OF ASSORTATIVE MATING AND SELFING WITH IN- AND OUTBREEDING DEPRESSION

The effect of inbreeding and outbreeding depression on the evolution of assortment are often considered separately. For instance, inbreeding depression is usually thought to shape selfing rates whereas outbreeding depression is commonly thought to affect the evolution of assortative mating. In this article, we consider the evolution of assortment in a context of local adaptation and we show that it is a typical situation in which both effects act simultaneously to shape the degree of selfing or assortative mating. More specifically, we show that selection on a modifier of mating can be partitioned into three distinct effects: a transmission advantage, an association to heterozygosity (proportional to inbreeding depression), and an association to beneficial alleles (proportional to outbreeding depression), so that random mating may evolve even with strong local adaptation. In addition, we show that it is necessary to carefully delimit the conditions for polymorphism at local adaptation locus to study the evolution of assortment. In particular, the range of parameters most favorable to the maintenance of polymorphism corresponds to situations favoring less assortment.     

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.     

Small effective population sizes in a widespread selfing species, Lymnaea truncatula (Gastropoda: Pulmonata)

We present here a spatial and temporal population genetic survey of a common freshwater snail, also a predominantly selfing species, Lymnaea truncatula. The rate of genetic diversity loss was quantified by estimating the effective size (Ne) of the snail populations, using two different methods. A temporal survey allowed estimation of a variance effective size of the populations, and a spatial survey allowed the estimation of an inbreeding effective size, from two-locus identity disequilibria estimates. Both methods were consistent and provided low Ne values. Drift due to (i) high amounts of selfing and (ii) fluctuations in population sizes because of temporary habitats, and also selection coupled to genome-wide linkage disequilibria, could explain such reductions in Ne. The loss of genetic diversity appears to be counterbalanced only very partially by low apparent rates of gene flow.     

Effective Size of Populations under Selection

Equations to approximate the effective size (N(e)) of populations under continued selection are obtained that include the possibility of partial full-sib mating and other systems such as assortative mating. The general equation for the case of equal number of sexes and constant number of breeding individuals (N) is N(e) = 4N/[2(1 - α(I)) + (S(k)(2) + 4Q(2)C(2)) (1 + α(I) + 2α(O))], where S(k)(2) is the variance of family size due to sampling without selection, C(2) is the variance of selective advantages among families (the squared coefficient of variation of the expected number of offspring per family), α(I) is the deviation from Hardy-Weinberg proportions, α(O) is the correlation between genes of male and female parents, and Q(2) is the term accounting for the cumulative effect of selection on an inherited trait. This is obtained as Q = 2/[2 - G(1 + r)], where G is the remaining proportion of genetic variance in selected individuals and r is the correlation of the expected selective values of male and female parents. The method is also extended to the general case of different numbers of male and female parents. The predictive value of the formulae is tested under a model of truncation selection with the infinitesimal model of gene effects, where C(2) and G are a function of the selection intensity, the heritability and the intraclass correlation of sibs. Under random mating r = α(I) = -1/(N - 1) and α(O) = 0. Under partial full-sib mating with an average proportion β of full-sib matings per generation, r & β and α(O) & α(I) & β/ (4 - 3β). The prediction equation is compared to other approximations based on the long-term contributions of ancestors to descendants. Finally, based on the approach followed, a system of mating (compensatory mating) is proposed to reduce rates of inbreeding without loss of response in selection programs in which selected individuals from the largest families are mated to those from the smallest families.     

RESPONSE TO SELECTION IN PARTIALLY SELF-FERTILIZING POPULATIONS. I. SELECTION ON A SINGLE TRAIT

Self-fertilization is a common form of reproduction in plants and it has important implications for quantitative trait evolution. Here, I present a model of selection on quantitative traits that can accommodate any level of self-fertilization. The “structured linear model” (SLM) predicts the evolution of the mean phenotype as a function of three distinct quantities: the mean additive genetic value, the directional dominance, and the mean inbreeding coefficient. Stochastic simulations of truncation selection demonstrate the accuracy of the SLM in predicting changes in the mean and variance of a quantitative trait over the full range of selfing rates. They also illustrate how complex interactions between selection and mating system determine the population distribution of inbreeding coefficients and also the amount of linkage disequilibrium. Changes in the genetic variance due to linkage disequilibria, which are commonly referred to as the “Bulmer effect,” are greatly magnified by selfing. This complicates the relationship between selfing rate and response to selection. Like the random mating theory, the parameters of the SLM can be estimated from phenotypic data.     

Postglacial ecotype formation under outcrossing and self‐fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate‐related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole‐genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate‐related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate‐related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome‐wide linkage disequilibrium.     

Functional pleiotropy and mating system evolution in plants: frequency-independent mating

Mutations that alter the morphology of floral displays (e.g., flower size) or plant development can change multiple functions simultaneously, such as pollen export and selfing rate. Given the effect of these various traits on fitness, pleiotropy may alter the evolution of both mating systems and floral displays, two characters with high diversity among angiosperms. The influence of viability selection on mating system evolution has not been studied theoretically. We model plant mating system evolution when a single locus simultaneously affects the selfing rate, pollen export, and viability. We assume frequency-independent mating, so our model characterizes prior selfing. Pleiotropy between increased viability and selfing rate reduces opportunities for the evolution of pure outcrossing, can favor complete selfing despite high inbreeding depression, and notably, can cause the evolution of mixed mating despite very high inbreeding depression. These results highlight the importance of pleiotropy for mating system evolution and suggest that selection by nonpollinating agents may help explain mixed mating, particularly in species with very high inbreeding depression.     

Genetic differentiation of neutral markers and quantitative traits in predominantly selfing metapopulations: confronting theory and experiments with Arabidopsis thaliana

The comparison of the genetic differentiation of quantitative traits (QST) and molecular markers (FST) can inform on the strength and spatial heterogeneity of selection in natural populations, provided that markers behave neutrally. However, selection may influence the behaviour of markers in selfing species with strong linkage disequilibria among loci, therefore invalidating this test of detection of selection. We address this issue by monitoring the genetic differentiation of five microsatellite loci (FST) and nine quantitative traits (QST) in experimental metapopulations of the predominantly selfing species Arabidopsis thaliana, that evolved during eight generations. Metapopulations differed with respect to population size and selection heterogeneity. In large populations, the genetic differentiation of neutral microsatellites was much larger under heterogeneous selection than under uniform selection. Using simulations, we show that this influence of selection heterogeneity on FST can be attributable to initial linkage disequilibria among loci, creating stronger genetic differentiation of QTL than expected under a simple additive model with no initial linkage. We found no significant differences between FST and QST regardless of selection heterogeneity, despite a demonstrated effect of selection on QST values. Additional data are required to validate the role of mating system and linkage disequilibria in the joint evolution of neutral and selected genetic differentiation, but our results suggest that FST/QST comparisons can be conservative tests to detect selection in selfing species.     

INBREEDING DEPRESSION IN A SELF-COMPATIBLE,ANDRODIOECIOUS CRUSTACEAN,EULIMNADIA TEXANA

The observation that offspring produced by the mating of close relatives are often less fit than those produced by matings between unrelated individuals (i.e., inbreeding depression) has commonly been explained in terms of the increased probability of expressing deleterious recessive alleles among inbred offspring (the partial dominance model). This model predicts that inbreeding depression should be limited in regularly inbreeding populations because the deleterious alleles that cause inbreeding depression (i.e., the genetic load) should be purged by regularly exposing these alleles to natural selection. We indirectly test the partial dominance model using four highly inbred populations of an androdioecious crustacean, the clam shrimp Eulimnadia texana. These shrimp are comprised of males and hermaphrodites, the latter capable of either self-fertilizing or mating with a male (i.e., outcrossing between hermaphrodites is impossible). Hermaphrodites are further subdivided into monogenics (produced via self-fertilization) and amphigenics (produced via self-fertilization or outcrossing). Electrophoretic evidence suggests significant differences in heterozygosity among populations, but that selfing rates were not statistically different (average s = 0.67). Additional electrophoretic analyses reveal that three previously described sex-linked loci (Fum, Idh-1, and Idh-2) are all tightly linked to each other, with crossing over on the order of 1% per generation. Although selfing rates are clearly high, we present evidence that early inbreeding depression (hatching rates, juvenile survival, and age at sexual maturity) exists in all four populations. For all of these factors, inbreeding depression was inferred by the positive correlation of multilocus heterozygosity and fitness. Cumulative inbreeding depression (8) is between 0.41 and 0.47 across all populations, which appears to be too low to limit the effects of purging via identity disequilibrium. Instead, we suggest that the maintenance of inbreeding depression in these populations is due to the observed linkage group, which we suggest contains a large number of genes including many related to fitness. Segregation of such a large linkage group would explain our observations of the predominance of amphigenic hermaphrodites in our field samples and of survival differences between monogenics and amphigenics within selfed clutches. We propose that a modified form of the overdominance model for inbreeding depression operating at the level of linkage groups maintains the observed levels of inbreeding depression in these populations even in the face of high rates of selfing.     

Change in frequency of a rare mutant allele: A general formula and applications to partial inbreeding models

 We deduce and prove a general formula to approximate the change in frequency of a mutant allele under weak selection, when this allele is introduced in small frequency into a population which was previously at a fixation state. We apply the formula to autosomal genes in partial selfing models and to autosomal as well as sex-linked genes in partial sib mating models. It is shown that the fate of a rare mutant allele depends not only on the selection parameters, the inbreeding coefficient and the reproductive values of the sexes in sex-differentiated populations, but also on coefficients of relatedness between mates. This is interpreted as a kin selection effect caused by inbreeding per se. Received: 3 December 2001 / Revised version: 10 April 2002 / Published online: 19 November 2002 Research supported in part by NSERC of Canada and FCAR of Québec. Mathematics Subject Classification (2000): Primary 60J80, Secondary 92D10, 92D25 Keywords or phrases: Adaptive topography – Partial selfing – Partial sib mating – Kin selection     

Effective Size of Nonrandom Mating Populations

Nonrandom mating whereby parents are related is expected to cause a reduction in effective population size because their gene frequencies are correlated and this will increase the genetic drift. The published equation for the variance effective size, Ne, which includes the possibility of nonrandom mating, does not take into account such a correlation, however. Further, previous equations to predict effective sizes in populations with partial sib mating are shown to be different, but also incorrect. In this paper, a corrected form of these equations is derived and checked by stochastic simulation. For the case of stable census number, N, and equal progeny distributions for each sex, the equation is [formula: see text], where Sk2 is the variance of family size and alpha is the departure from Hardy-Weinberg proportions. For a Poisson distribution of family size (Sk2 = 2), it reduces to Ne = N/(1 + alpha), as when inbreeding is due to selfing. When nonrandom mating occurs because there is a specified system of partial inbreeding every generation, alpha can be substituted by Wright's FIS statistic, to give the effective size as a function of the proportion of inbred mates.     

The Effect of Inbreeding on Competitive Male-Mating Ability in DROSOPHILA MELANOGASTER

The effect of full-sib inbreeding on competitive male-mating ability (CI♂) in Drosophila melanogaster was investigated in two experiments. In the first, five inbred lines (with reserves) were assessed up to 18 generations. Linear inbreeding depression, of 5.9% per 10% increase in homozygosity, was observed. In a second experiment, 21 inbred lines were tested after three generations of full-sib mating (without reserves), and the decline with inbreeding was more severe, the male competitive index (CI♂) decreasing by 10.7% per 10% increase in F. The difference between these results is attributed to natural selection acting on variation within the inbred lines in extent of homozygosity, which can arise because of the peculiarly strong influence of linkage in Drosophila. Furthermore, differentiation between the lines may have reflected this variation rather than the various effects of different alleles fixed.—These results imply that the genetic variation in male-mating ability is largely due to dominance (no epistasis was detected) and are consonant with the proposition that intermale sexual selection is a very important component of fitness in D. melanogaster . There was no evidence of a positive correlation between male body size and competitive mating ability.     

Joint evolution of differential seed dispersal and self‐fertilization

Differential seed dispersal, in which selfed and outcrossed seeds possess different dispersal propensities, represents a potentially important individual‐level association. A variety of traits can mediate differential seed dispersal, including inflorescence and seed size variation. However, how natural selection shapes such associations is poorly known. Here, we developed theoretical models for the evolution of mating system and differential seed dispersal in metapopulations, incorporating heterogeneous pollination, dispersal cost, cost of outcrossing and environment‐dependent inbreeding depression. We considered three models. In the ‘fixed dispersal model’, only selfing rate is allowed to evolve. In the ‘fixed selfing model’, in which selfing is fixed but differential seed dispersal can evolve, we showed that natural selection favours a higher, equal or lower dispersal rate for selfed seeds to that for outcrossed seeds. However, in the ‘joint evolution model’, in which selfing and dispersal can evolve together, evolution necessarily leads to higher or equal dispersal rate for selfed seeds compared to that for outcrossed. Further comparison revealed that outcrossed seed dispersal is selected against by the evolution of mixed mating or selfing, whereas the evolution of selfed seed dispersal undergoes independent processes. We discuss the adaptive significance and constraints for mating system/dispersal association.