Estimation of deleterious genomic mutation parameters in natural populations by accounting for variable mutation effects across loci

The genomes of all organisms are subject to continuous bombardment of deleterious genomic mutations (DGM). Our ability to accurately estimate various parameters of DGM has profound significance in population and evolutionary genetics. The Deng-Lynch method can estimate the parameters of DGM in natural selfing and outcrossing populations. This method assumes constant fitness effects of DGM and hence is biased under variable fitness effects of DGM. Here, we develop a statistical method to estimate DGM parameters by considering variable mutation effects across loci. Under variable mutation effects, the mean fitness and genetic variance for fitness of parental and progeny generations across selfing/outcrossing in outcrossing/selfing populations and the covariance between mean fitness of parents and that of their progeny are functions of DGM parameters: the genomic mutation rate U, average homozygous effect s, average dominance coefficient h, and covariance of selection and dominance coefficients cov(h, s). The DGM parameters can be estimated by the algorithms we developed herein, which may yield improved estimation of DGM parameters over the Deng-Lynch method as demonstrated by our simulation studies. Importantly, this method is the first one to characterize cov(h, s) for DGM.     

Estimation of parameters of deleterious mutations in partial selfing or partial outcrossing populations and in nonequilibrium populations

The Deng-Lynch method was developed to estimate the rate and effects of deleterious genomic mutations (DGM) in natural populations under the assumption that populations are either completely outcrossing or completely selfing and that populations are at mutation-selection (M-S) balance. However, in many plant and animal populations, selfing or outcrossing is often incomplete in that a proportion of populations undergo inbreeding while the rest are outcrossing. In addition, the degrees of deviation of populations from M-S balance are often not known. Through computer simulations, we investigated the robustness and the applicability of the Deng-Lynch method under different degrees of partial selfing or partial outcrossing and for nonequilibrium populations approaching M-S balance at different stages. The investigation was implemented under constant, variable, and epistatic mutation effects. We found that, generally, the estimation by the Deng-Lynch method is fairly robust if the selfing rate (S) is <0.10 in outcrossing populations and if S > 0.8 in selfing populations. The estimation may be unbiased under partial selfing with variable and epistatic mutation effects in predominantly outcrossing populations. The estimation is fairly robust in nonequilibrium populations at different stages approaching M-S balance. The dynamics of populations approaching M-S balance under various parameters are also studied. Under mutation and selection, populations approach balance at a rapid pace. Generally, it takes 400-2000 generations to reach M-S balance even when starting from homogeneous individuals free of DGM. Our investigation here provides a basis for characterizing DGM in partial selfing or outcrossing populations and for nonequilibrium populations.     

Mutation accumulation in a selfing population: consequences of different mutation rates between selfers and outcrossers

Currently existing theories predict that because deleterious mutations accumulate at a higher rate, selfing populations suffer from more intense genetic degradation relative to outcrossing populations. This prediction may not always be true when we consider a potential difference in deleterious mutation rate between selfers and outcrossers. By analyzing the evolutionary stability of selfing and outcrossing in an infinite population, we found that the genome-wide deleterious mutation rate would be lower in selfing than in outcrossing organisms. When this difference in mutation rate was included in simulations, we found that in a small population, mutations accumulated more slowly under selfing rather than outcrossing. This result suggests that under frequent and intense bottlenecks, a selfing population may have a lower risk of genetic extinction than an outcrossing population.     

The relative importance of factors determining genetic drift: mating system, spatial genetic structure, habitat and census size in Arabidopsis lyrata

? The mating system, dispersal and census size are predicted to determine the magnitude of genetic drift, but little is known about their relative importance in nature. ? We estimated the contributions of several population-level features to genetic drift in 18 populations of Arabidopsis lyrata. The factors were outcrossing rate, within-population spatial genetic structure, census size and substrate type. The expected heterozygosity (H(E)) at 10 microsatellite loci was taken to reflect the effective population size (N(e)) and the strength of genetic drift. ? The mating system explained most of the variation in H(E) (60%), followed by substrate (10%), genetic structure (9%) and census size (6%). The most outcrossing population had a +0.32 higher predicted H(E) than the most selfing population; the estimated N(e) of selfing populations was less than half that of outcrossing populations. Rocky outcrops supported populations with a +0.14 higher H(E) than did sandy substrates. The most structured population had a +0.24 higher H(E) than the least structured population, and the largest population had a +0.18 higher H(E) than the smallest population. ? This study illustrates the importance of outcrossing, genetic structure and the physical environment--together with census size--in maintaining H(E), and suggests that multiple population-level characteristics influence N(e) and the action of genetic drift.     

Differences in floral traits and flower visitation rates in mating systems in Prunella vulgaris (Lamiaceae)

夏枯草交配系统对花特征和访花频率差异的影响 植物花特征和传粉者的访问次数与交配系统类型密切相关。唇形科植物夏枯草(Prunella vulgaris)存 在两种植株类型,分别为柱头伸出花冠和柱头在花冠内部的植株,而且两种植株的比例在不同种群中存在差异。本研究选择柱头伸出花冠外植株占绝大多数、柱头伸出花冠外植株占多数和柱头在花冠内部植株占多数的3个种群,通过比较每个种群中两种植株类型的开花物候、花形态特征、昆虫访问频率、自交能力、传粉者对结实的贡献以及近交衰退的水平,以检验花特征和传粉者访问次数与交配系统类型的关系。研究结果表明,与柱头在花冠内部的植株相比,柱头伸出花冠外的植株具有更大和更多的花,产生更多的花粉和花蜜,具有更高的访花频率,并主要通过异交产生种子。在种群水平,柱头伸出花冠外的植株占多数种群的访花频率显著高于柱头在花冠内部植株占多数的访花频率。柱头在花冠内部的植株比柱头伸出花冠外的植株具有更强的自动自交能力,在传粉者缺乏时为其提供了繁殖保障,但繁殖保障和异交率在不同种群中差异不显著,表明较低的昆虫访问能够满足夏枯草的授粉需求以产生种子,这可能与夏枯草较少的胚珠数量(每朵花仅有4个胚珠)有关。柱头在花冠内部植株的近交衰退水平低于柱头伸出花冠外植株的近交衰退水平,但两种植株类型的近交衰退水平均低于0.5,说明近交衰退不足以阻止该物种中自交的进化。综上所述,柱头在花冠内部的植株能够通过自交为夏枯草提供繁殖保障,而柱头伸出花冠外的植株能够利用昆虫传粉确保异交,表明混合交配系统在该物种中是一个稳定的状态。     

Demographic signatures accompanying the evolution of selfing in Leavenworthia alabamica

The evolution of selfing from outcrossing is a common transition, yet little is known about the mutations and selective factors that promote this shift. In the mustard family, single-locus self-incompatibility (SI) enforces outcrossing. In this study, we test whether mutations causing self-compatibility (SC) are linked to the self-incompatibility locus (S-locus) in Leavenworthia alabamica, a species where two selfing races (a2 and a4) co-occur with outcrossing populations. We also infer the ecological circumstances associated with origins of selfing using molecular sequence data. Genealogical reconstruction of the Lal2 locus, the putative ortholog of the SRK locus, showed that both selfing races are fixed for one of two different S-linked Lal2 sequences, whereas outcrossing populations harbor many S-alleles. Hybrid crosses demonstrated that S-linked mutations cause SC in each selfing race. These results strongly suggest two origins of selfing in this species, a result supported by population admixture analysis of 16 microsatellite loci and by a population tree built from eight nuclear loci. One selfing race (a4) shows signs of a severe population bottleneck, suggesting that reproductive assurance might have caused the evolution of selfing in this case. In contrast, the population size of race a2 cannot be distinguished from that of outcrossing populations after correcting for differences in selfing rates. Coalescent-based analyses suggest a relatively old origin of selfing in the a4 race (～150 ka ago), whereas selfing evolved recently in the a2 race (～12-48 ka ago). These results imply that S-locus mutations have triggered two recent shifts to selfing in L. alabamica, but that these transitions are not always associated with a severe population bottleneck, suggesting that factors other than reproductive assurance may play a role in its evolution.     

Small reductions in corolla size and pollen: ovule ratio,but no changes in flower shape in selfing populations of the North American <Emphasis Type="Italic">Arabidopsis lyrata</Emphasis>

The shift from outcrossing to selfing is often accompanied by striking changes in floral morphology towards a “selfing syndrome”, which is characterized by flowers with reduction in size, pollen: ovule (P/O) ratio, and herkogamy. This study aims to test whether such changes have occurred in the North American Arabidopsis lyrata, which is of particular interest because of the relatively recent transitions to selfing in this system. Flower size, flower shape, herkogamy levels, P/O ratio, and floral integration of six self-incompatible (outcrossing) and six self-compatible (selfing) populations of A. lyrata were measured in a common environment using conventional and geometric morphometrics methods. Although selfers had on average 9.2% smaller corollas, 8.4% longer pistils, and 21.5% lower P/O ratios than outcrossers, there were no differences in shape, floral integration, and herkogamy between outcrossing and selfing populations. Moreover, most variation in floral traits was explained by population genetic background rather than by mating system. We conclude that selfing populations in A. lyrata have not evolved a selfing syndrome.     

THE EVOLUTION OF SELF-FERTILIZATION AND INBREEDING DEPRESSION IN PLANTS. II. EMPIRICAL OBSERVATIONS

A bimodal distribution of outcrossing rates was observed for natural plant populations, with more primarily selfing and primarily outcrossing species, and fewer species with intermediate outcrossing rate than expected by chance. We suggest that this distribution results from selection for the maintenance of outcrossing in historically large, outcrossing populations with substantial inbreeding depression, and from selection for selfing when increased inbreeding, due to pollinator failure or population bottlenecks, reduces the level of inbreeding depression. Few species or populations are fixed at complete selfing or complete outcrossing. A low level of selfing in primarily outcrossing species is unlikely to be selectively advantageous, but will not reduce inbreeding depression to the level where selfing is selectively favored, particularly if accompanied by reproductive compensation. Similarly, occasional outcrossing in primarily selfing species is unlikely to regularly provide sufficient heterosis to maintain selection for outcrossing through individual selection. Genetic, morphological and ecological constraints may limit the potential for outcrossing rates in selfers to be reduced below some minimum level.     

The evolutionary response of mating system to heterosis

Isolation allows populations to diverge and to fix different alleles. Deleterious alleles that reach locally high frequencies contribute to genetic load, especially in inbred or selfing populations, in which selection is relaxed. In the event of secondary contact, the recessive portion of the genetic load is masked in the hybrid offspring, producing heterosis. This advantage, only attainable through outcrossing, should favour evolution of greater outcrossing even if inbreeding depression has been purged from the contributing populations. Why, then, are selfing‐to‐outcrossing transitions not more common? To evaluate the evolutionary response of mating system to heterosis, we model two monomorphic populations of entirely selfing individuals, introduce a modifier allele that increases the rate of outcrossing and investigate whether the heterosis among populations is sufficient for the modifier to invade and fix. We find that the outcrossing mutation invades for many parameter choices, but it rarely fixes unless populations harbour extremely large unique fixed genetic loads. Reversions to outcrossing become more likely as the load becomes more polygenic, or when the modifier appears on a rare background, such as by dispersal of an outcrossing genotype into a selfing population. More often, the outcrossing mutation instead rises to moderate frequency, which allows recombination in hybrids to produce superior haplotypes that can spread without the mutation's further assistance. The transience of heterosis can therefore explain why secondary contact does not commonly yield selfing‐to‐outcrossing transitions.     

Limited phenological and pollinator-mediated isolation among selfing and outcrossing Arabidopsis lyrata populations

Transitions from outcrossing to selfing have been a frequent evolutionary shift in plants and clearly play a role in species divergence. However, many questions remain about the initial mechanistic basis of reproductive isolation during the evolution of selfing. For instance, how important are pre-zygotic pre-pollination mechanisms (e.g. changes in phenology and pollinator visitation) in maintaining reproductive isolation between newly arisen selfing populations and their outcrossing ancestors? To test whether changes in phenology and pollinator visitation isolate selfing populations of Arabidopsis lyrata from outcrossing populations, we conducted a common garden experiment with plants from selfing and outcrossing populations as well as their between-population hybrids. Specifically, we asked whether there was isolation between outcrossing and selfing plants and their between-population hybrids through differences in (1) the timing or intensity of flowering; and/or (2) pollinator visitation. We found that phenology largely overlapped between plants from outcrossing and selfing populations. There were also no differences in pollinator preference related to mating system. Additionally, pollinators preferred to visit flowers on the same plant rather than exploring nearby plants, creating a large opportunity for self-fertilization. Overall, this suggests that pre-zygotic pre-pollination mechanisms do not strongly reproductively isolate plants from selfing and outcrossing populations of Arabidopsis lyrata.     

THE EVOLUTION OF SELF-FERTILIZATION AND INBREEDING DEPRESSION IN PLANTS. I. GENETIC MODELS

The amounts of inbreeding depression upon selfing and of heterosis upon outcrossing determine the strength of selection on the selfing rate in a population when this evolves polygenically by small steps. Genetic models are constructed which allow inbreeding depression to change with the mean selfing rate in a population by incorporating both mutation to recessive and partially dominant lethal and sublethal alleles at many loci and mutation in quantitative characters under stabilizing selection. The models help to explain observations of high inbreeding depression (> 50%) upon selfing in primarily outcrossing populations, as well as considerable heterosis upon outcrossing in primarily selfing populations. Predominant selfing and predominant outcrossing are found to be alternative stable states of the mating system in most plant populations. Which of these stable states a species approaches depends on the history of its population structure and the magnitude of effect of genes influencing the selfing rate.     

Mating system and the evolution of quantitative traits: an experimental study of Mimulus guttatus

Abstract The mating system of a population profoundly influences its evolution. Inbreeding alters the balance of evolutionary forces that determine the amount of genetic variation within a population. It redistributes that variation among individuals, altering heritabilities and genetic correlations. Inbreeding even changes the basic relationships between these genetic statistics and response to selection. If populations differing only in mating system are exposed to the same selection pressures, will they respond in qualitatively different ways? Here, we address this question by imposing selection on an index of two negatively correlated traits (flower size and development rate) within experimental populations that reproduce entirely by outcrossing, entirely by self‐fertilizing, or by a mixture of outcrossing and selfing. Entirely selfing populations responded mainly by evolving larger flowers whereas outcrossing populations also evolved more rapid development. Divergence occurred despite an equivalent selection regime and no direct effect of mating system on fitness. The study provides an experimental demonstration of how the interaction of selection, genetic drift, and mating system can produce dramatic short‐term changes in trait means, variances, and covariances.     

Evolution of the Selfing Syndrome in Arabis alpina (Brassicaceae)

Introduction

The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.

Methods

Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.

Results and Discussion

Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.

Conclusion

We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.     

Variation in Outcrossing Rates in Eichhornia paniculata: Temporal Changes in Populations of Contrasting Style Morph Structure

Abstract Surveys of mating-system parameters in populations of the annual, self-compatible, tristylous, emergent aquatic, Eichhornia paniculata (Pontederiaceae) from N.E. Brazil and Jamaica have indicated that the species exhibits a wide range of outcrossing rates. To investigate whether temporal variation in outcrossing rate was also a feature of populations, open-pollinated families were sampled from five populations of contrasting style morph structure from N.E. Brazil over three consecutive years (1987–1989). Multilocus estimates of outcrossing rate ( t ) were obtained from assays of isozyme polymorphisms using starch gel electrophoresis. There was significant variation both among populations and between years in the frequency of outcrossing. Outcrossing in three tristylous populations was high ( t > 0.80), with relatively small fluctuations occurring over the three-year sampling period. In contrast, in a dimorphic and monomorphic population considerable self-fertilization occurred and the frequency of outcrossing declined significantly from 1987 to 1989 in both populations. In the dimorphic population, increased selfing was associated with a marked reduction in population size and an increase in the frequency of selfing variants of the mid-styled morph. The significance of temporal variation in outcrossing frequency in plant populations is discussed in relation to its effect on population genetic structure and recent models of mating-system evolution.     

Fine-scale estimation of outcrossing in western redcedar with microsatellite assay of bulked DNA

Western redcedar (Thuja plicata, Cupressaceae) is a self-fertile conifer with a mixed mating system and significant variation for outcrossing among populations. In this paper, we conducted a fine-scale study of mating system variation to identify correlates of outcrossing in natural populations. We examined variation for outcrossing within and among individual trees, and describe a new method to estimate outcrossing using bulked DNA samples. Bulking (assaying DNA tissues from several individuals simultaneously) increases the experimental power without increasing the experimental effort. We sampled 80 trees from four natural populations in southwestern British Columbia. From each tree, we sampled from up to six crown positions (three heights and inner vs outer branches). From each position, two samples of three seedlings each were bulked before DNA extractions. Using four microsatellite loci, we obtained outcrossing rates for each tree and for each of the six crown positions. We found individual tree selfing rates to increase with tree height in all four populations, but selfing rates did not differ among crown positions. The higher selfing rate of larger trees is probably due to their greater proportional contribution to local pollen clouds. Individual tree outcrossing rates ranged from 22 to 100% and the population outcrossing rates from 66 to 78%. Missed alleles due to bulking and the estimation method used both cause a downward bias in outcrossing rates, so that these estimates are probably lower than the actual outcrossing rates. Nevertheless, the trends we observed are not affected by systematic biases of estimation.     

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.     

The similarity index and DNA fingerprinting

DNA-fingerprint similarity is being used increasingly to make inferences about levels of genetic variation within and between natural populations. It is shown that the similarity index--the average fraction of shared restriction fragments--provides upwardly biased estimates of population homozygosity but nearly unbiased estimates of the average identity-in-state for random pairs of individuals. A method is suggested for partitioning the DNA-fingerprint dissimilarity into within- and between-population components. Some simple expressions are given for the sampling variances of these estimators.     

BREEDING SYSTEM EVOLUTION IN LEAVENWORTHIA: BREEDING SYSTEM VARIATION AND REPRODUCTIVE SUCCESS IN NATURAL POPULATIONS OF LEAVENWORTHIA CRASSA (CRUCIFERAE)

Populations of Leavenworthia crassa (Cruciferae) studied for 3 years exhibited among- and within-population genetic variation for a suite of floral and reproductive traits (flower width, petal length, anther position, ability to set seeds in the absence of pollinators, time to first flowering) associated with breeding system. We used electrophoretic markers to show that a population with small, monomorphically colored flowers with introrse anthers had a significantly lower outcrossing rate (t = 0.03) than did a population with larger, polymorphically colored flowers with extrorse anthers (t = 0.33). In the more-outcrossing population the correlation between higher maternal plant outcrossing rate and the suite of six traits approached significance (P < 0.067), with greater petal size, greater flower width, and reduced ability to set seeds in the absence of pollinators contributing significantly. Plants in selfing populations had a generally higher reproductive success, with a higher number of flowers per plant, a smaller proportion of unfertilized ovules, a smaller proportion of fertilized ovules aborted, a higher rate of fruit set, and overall a larger number of seeds matured than did plants from the more outcrossing populations. Pollen limitation did not appear to account for lower reproductive success in outcrossing populations. Resource limitation did not differ substantially between populations. However, within-ovary patterns of fertilization, abortion, and seed weight were significantly less random in outcrossing populations than in selfing populations, suggesting that differential gamete and embryo success may be responsible for lower reproductive success in outcrossing populations.     

Inferring recent outcrossing rates using multilocus individual heterozygosity: application to evolving wheat populations

Using multilocus individual heterozygosity, a method is developed to estimate the outcrossing rates of a population over a few previous generations. Considering that individuals originate either from outcrossing or from n successive selfing generations from an outbred ancestor, a maximum-likelihood (ML) estimator is described that gives estimates of past outcrossing rates in terms of proportions of individuals with different n values. Heterozygosities at several unlinked codominant loci are used to assign n values to each individual. This method also allows a test of whether populations are in inbreeding equilibrium. The estimator's reliability was checked using simulations for different mating histories. We show that this ML estimator can provide estimates of outcrossing rates for the final generation outcrossing rate (t(0)) and a mean of the preceding rates (t(p)) and can detect major temporal variation in the mating system. The method is most efficient for low to intermediate outcrossing levels. Applied to nine populations of wheat, this method gave estimates of t(0) and t(p). These estimates confirmed the absence of outcrossing t(0) = 0 in the two populations subjected to manual selfing. For free-mating wheat populations, it detected lower final generation outcrossing rates t(0) = 0-0.06 than those expected from global heterozygosity t = 0.02-0.09. This estimator appears to be a new and efficient way to describe the multilocus heterozygosity of a population, complementary to Fis and progeny analysis approaches.     

Postcopulatory sexual selection reduces genetic diversity in experimental populations of Caenorhabditis elegans

Postcopulatory sexual selection affects the evolution of numerous features ranging from mating behavior to seminal fluid toxicity to the size of gametes. In an earlier study of the effect of sperm competition risk on sperm size evolution, experimental populations of the nematode Caenorhabditis elegans were maintained either by outcrossing (sperm competition present) or by selfing (no sperm competition), and after 60 generations, significantly larger sperm had evolved in the outcrossing populations. To determine the effects of this selection on population genetic variation, we assessed genetic diversity in a large number of loci using random amplification of polymorphic DNA-PCR. Nearly 80% of the alleles present in parental strain populations persisted in the 6 experimental populations after the 60 generations and, despite a 2.2-fold difference in expected heterozygosity, the resulting levels of genetic variation were equivalent between the outcrossing and selfing experimental populations. By inference, we conclude that genetic hitchhiking due to sexual selection in the experimental populations dramatically reduced genetic diversity. We use the levels of variation in the selfing populations as a control for the effects of drift, and estimate the strength of sexual selection to be strong in obligatorily outcrossing populations. Although sequential hermaphrodites like C. elegans probably experience little sexual selection in nature, these data suggest that sexual selection can profoundly affect diversity in outcrossing taxa.