Maintenance of Quantitative Genetic Variance Under Partial Self-Fertilization, with Implications for Evolution of Selfing

Polygenic selection to a changing optimum under self-fertilisation

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

Weak response to selection on stigma-anther distance in a primarily selfing population of yellow monkeyflower

Stebbins hypothesized that selfing lineages are evolutionary dead ends because they lack adaptive potential. While selfing populations often possess limited nucleotide variability compared with closely related outcrossers, reductions in the genetic variability of quantitative characters remain unclear, especially for key traits determining selfing rates. Yellow monkeyflower (Mimulus guttatus) populations generally outcross and maintain extensive quantitative genetic variation in floral traits. Here, we study the Joy Road population (Bodega Bay, CA, USA) of M. guttatus, where individuals exhibit stigma-anther distances (SAD) typical of primarily selfing monkeyflowers. We show that this population is closely related to nearby conspecifics on the Pacific Coast with a modest 33% reduction in genome-wide variation compared with a more highly outcrossing population. A five-generation artificial selection experiment challenged the hypothesis that the Joy Road population harbours comparatively low evolutionary potential in stigma-anther distance, a critical determinant of selfing rate in Mimulus. Artificial selection generated a weak phenotypic response, with low realized heritabilities (0.020-0.028) falling 84% below those measured for floral characters in more highly outcrossing M. guttatus. These results demonstrate substantial declines in evolutionary potential with a transition toward selfing. Whether these findings explain infrequent reversals to outcrossing or general limits on adaptation in selfers requires further investigation.

Early-acting inbreeding depression can evolve as an inbreeding avoidance mechanism

Despite the potential for mechanical, developmental and/or chemical mechanisms to prevent self-fertilization, incidental self-fertilization is inevitable in many predominantly outcrossing species. In such cases, inbreeding can compromise individual fitness. Unquestionably, much of this inbreeding depression is maladaptive. However, we show that when reproductive compensation allows for the replacement of inviable embryos lost early in development, selection can favour deleterious recessive variants that induce 'self-sacrificial' death of inbred embryos. Our theoretical results provide numerous testable predictions which could challenge the assumption that inbreeding depression is always maladaptive. Our work is applicable any species that cannot fully avoid inbreeding, exhibits substantial inbreeding depression, and has the potential to compensate embryos lost early in development. In addition to its general applicability, our theory suggests that self-sacrificial variants might be responsible for the remarkably low realized selfing rates of gymnosperms with high primary selfing rates, as gymnosperms exhibit strong inbreeding depression, have effective reproductive compensation mechanisms, and cannot evolve chemical self-incompatibility.

Gametic selection favours polyandry and selfing

Competition among pollen or sperm (gametic selection) can cause evolution. Mating systems shape the intensity of gametic selection by determining the competitors involved, which can in turn cause the mating system itself to evolve. We model the bidirectional relationship between gametic selection and mating systems, focusing on variation in female mating frequency (monandry-polyandry) and self-fertilisation (selfing-outcrossing). First, we find that monandry and selfing both reduce the efficiency of gametic selection in removing deleterious alleles. This means that selfing can increase mutation load, in contrast to cases without gametic selection where selfing purges deleterious mutations and decreases mutation load. Second, we explore how mating systems evolve via their effect on gametic selection. By manipulating gametic selection, polyandry can evolve to increase the fitness of the offspring produced. However, this indirect advantage of post-copulatory sexual selection is weak and is likely to be overwhelmed by any direct fitness effects of mating systems. Nevertheless, gametic selection can be potentially decisive for selfing evolution because it significantly reduces inbreeding depression, which favours selfing. Thus, the presence of gametic selection could be a key factor driving selfing evolution.

The infinitesimal model with dominance

The classical infinitesimal model is a simple and robust model for the inheritance of quantitative traits. In this model, a quantitative trait is expressed as the sum of a genetic and an environmental component, and the genetic component of offspring traits within a family follows a normal distribution around the average of the parents' trait values, and has a variance that is independent of the parental traits. In previous work, we showed that when trait values are determined by the sum of a large number of additive Mendelian factors, each of small effect, one can justify the infinitesimal model as a limit of Mendelian inheritance. In this paper, we show that this result extends to include dominance. We define the model in terms of classical quantities of quantitative genetics, before justifying it as a limit of Mendelian inheritance as the number, M, of underlying loci tends to infinity. As in the additive case, the multivariate normal distribution of trait values across the pedigree can be expressed in terms of variance components in an ancestral population and probabilities of identity by descent determined by the pedigree. Now, with just first-order dominance effects, we require two-, three-, and four-way identities. We also show that, even if we condition on parental trait values, the "shared" and "residual" components of trait values within each family will be asymptotically normally distributed as the number of loci tends to infinity, with an error of order 1/M. We illustrate our results with some numerical examples.

Can heterosis and inbreeding depression explain the maintenance of outcrossing in a cleistogamous perennial?

PREMISE

What maintains mixed mating is an evolutionary enigma. Cleistogamy-the production of both potentially outcrossing chasmogamous and obligately selfing cleistogamous flowers on the same individual plant-is an excellent system to study the costs of selfing. Inbreeding depression can prevent the evolution of greater selfing within populations, and heterosis in crosses between populations may further tip the balance in favor of outcrossing. Few empirical estimates of inbreeding depression and heterosis in the same system exist for cleistogamous species.

METHODS

We investigate the potential costs of selfing by quantifying inbreeding depression and heterosis in three populations of the cleistogamous perennial Ruellia humilis Nutt (Acanthaceae). We performed three types of hand-pollinations-self, outcross-within, and outcross-between populations-and measured seed number, germination, total flower production, and estimated cumulative fitness for the resulting progeny in a greenhouse experiment.

RESULTS

We found moderate inbreeding depression for cumulative fitness (<30%) in two populations, but outbreeding depression for crosses within a third population (-26%). For between-population crosses, there was weak to modest heterosis (11-47%) in two of the population combinations, but modest to strong outbreeding depression (-21 to -71%) in the other four combinations.

CONCLUSIONS

Neither inbreeding depression nor heterosis was of sufficient magnitude to explain the continued production of chasmogamous flowers given the relative energetic advantage of cleistogamous flowers previously estimated for these populations. Outbreeding depression either within or between populations makes the maintenance of chasmogamous flowers even harder to explain. More information is needed on the genetic basis of cleistogamy to resolve this conundrum.

Influence of the mutation load on the genomic composition of hybrids between outcrossing and self-fertilizing species

Hybridization is a natural process whereby two diverging evolutionary lineages reproduce and create offspring of mixed ancestry. Differences in mating systems (e.g., self-fertilization and outcrossing) are expected to affect the direction and extent of hybridization and introgression in hybrid zones. Among other factors, selfers and outcrossers are expected to differ in their mutation loads. This has been studied both theoretically and empirically; however, conflicting predictions have been made on the effects mutation loads of parental species with different mating systems can have on the genomic composition of hybrids. Here, we develop a multi-locus, selective model to study how the different mutation load built up in selfers and outcrossers as a result of selective interference and homozygosity impact the long-term genetic composition of hybrid populations. Notably, our results emphasize that genes from the parental population with lesser mutation load get rapidly overrepresented in hybrid genomes, regardless of the hybrids own mating system. When recombination tends to be more important than mutation, outcrossers' genomes tend to be of higher quality and prevail. When recombination rates are low, however, selfers' genomes may reach higher quality than outcrossers' genomes and prevail in the hybrids. Taken together, these results provide concrete insights into one of the multiple factors influencing hybrid genome ancestry and introgression patterns in hybrid zones containing species with different mating systems.

Evolvability and trait function predict phenotypic divergence of plant populations

Understanding the causes and limits of population divergence in phenotypic traits is a fundamental aim of evolutionary biology, with the potential to yield predictions of adaptation to environmental change. Reciprocal transplant experiments and the evaluation of optimality models suggest that local adaptation is common but not universal, and some studies suggest that trait divergence is highly constrained by genetic variances and covariances of complex phenotypes. We analyze a large database of population divergence in plants and evaluate whether evolutionary divergence scales positively with standing genetic variation within populations (evolvability), as expected if genetic constraints are evolutionarily important. We further evaluate differences in divergence and evolvability-divergence relationships between reproductive and vegetative traits and between selfing, mixed-mating, and outcrossing species, as these factors are expected to influence both patterns of selection and evolutionary potentials. Evolutionary divergence scaled positively with evolvability. Furthermore, trait divergence was greater for vegetative traits than for floral (reproductive) traits, but largely independent of the mating system. Jointly, these factors explained ~40% of the variance in evolutionary divergence. The consistency of the evolvability-divergence relationships across diverse species suggests substantial predictability of trait divergence. The results are also consistent with genetic constraints playing a role in evolutionary divergence.

The potential for floral evolution in response to competing selection pressures following the loss of hawkmoth pollination in Ruellia humilis

PREMISE

In the absence of hawkmoth pollinators, chasmogamous (CH) flowers of Ruellia humilis self-pollinate by two secondary mechanisms. Other floral visitors might exert selection on CH floral traits to restore outcrossing, but at the same time preferential predation of CH seeds generates selection to increase the allocation of resources to cleistogamous (CL) flowers.

METHODS

To assess the potential for an evolutionary response to these competing selection pressures, we estimated additive genetic variances ( σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ ) and covariances for 14 reproductive traits and three fitness components in a Missouri population lacking hawkmoth pollinators.

RESULTS

We found significant σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ for all 11 floral traits and two measures of resource allocation to CL flowers, indicating the potential for a short-term response to selection on most reproductive traits. Selection generated by seed predators is predicted to increase the percentage of CL flowers by 0.24% per generation, and mean stigma-anther separation is predicted to decrease as a correlated response, increasing the fraction of plants that engage in prior selfing. However, the initial response to this selection is opposed by strong directional dominance.

CONCLUSIONS

The predicted evolutionary decrease in the number of CH flowers available for potential outcrossing, combined with the apparent preclusion of potential diurnal pollinators by the pollen-harvesting activities of sweat bees, suggest that 100% cleistogamy is the likely outcome of evolution in the absence of hawkmoths. However, rare mutations with large effects, such as delaying budbreak until after sunrise, could provide pathways for the restoration of outcrossing that are not reachable by gradual quantitative-genetic evolution.

Polyploidization: Consequences of genome doubling on the evolutionary potential of populations

Whole-genome duplication is common in plants and is considered to have a broad range of effects on individuals' phenotypes and genomes and to be an important driver of plant adaptation and speciation. Despite their increased capacity to cope with challenging environments, polyploid lineages are generally as prone to extinction, and sometimes more prone, than their diploid progenitors. Although several explanations have been proposed to explain the short- and long-term disadvantages of polyploidy on the survival probability of populations, the consequences of whole-genome doubling on the heritable variance remain poorly studied. Whole-genome doubling can have major effects not only on the genetics, but also on the ecology and life history of the populations. Modifications of other properties of populations can reverse the effects of polyploidization per se on heritable variance. In this synthesis, I summarize the empirical and theoretical knowledge about the multifarious consequences of genome doubling on the heritable variance of quantitative traits and on the evolutionary potential of polyploid populations compared to their diploid progenitors. I propose several ways to decipher the consequences of whole-genome doubling on survival probability and to study the further consequences of shifting the ecological niche and life-history traits of a population. I also highlight some practical considerations for comparing the heritable variance of a trait among different cytotypes. Such investigations appear to be timely and necessary to understand more about the paradoxical aspects of polyploidization and to understand the evolutionary potential of polyploid lineages in a global warming context.

Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift

Resurrection studies are a useful tool to measure how phenotypic traits have changed in populations through time. If these trait modifications correlate with the environmental changes that occurred during the time period, it suggests that the phenotypic changes could be a response to selection. Selfing, through its reduction of effective size, could challenge the ability of a population to adapt to environmental changes. Here, we used a resurrection study to test for adaptation in a selfing population of Medicago truncatula, by comparing the genetic composition and flowering times across 22 generations. We found evidence for evolution toward earlier flowering times by about two days and a peculiar genetic structure, typical of highly selfing populations, where some multilocus genotypes (MLGs) are persistent through time. We used the change in frequency of the MLGs through time as a multilocus fitness measure and built a selection gradient that suggests evolution toward earlier flowering times. Yet, a simulation model revealed that the observed change in flowering time could be explained by drift alone, provided the effective size of the population is small enough (<150). These analyses suffer from the difficulty to estimate the effective size in a highly selfing population, where effective recombination is severely reduced.

Evolvability: A Quantitative-Genetics Perspective

The concept of evolvability emerged in the early 1990s and soon became fashionable as a label for different streams of research in evolutionary biology. In evolutionary quantitative genetics, evolvability is defined as the ability of a population to respond to directional selection. This differs from other fields by treating evolvability as a property of populations rather than organisms or lineages and in being focused on quantification and short-term prediction rather than on macroevolution. While the term evolvability is new to quantitative genetics, many of the associated ideas and research questions have been with the field from its inception as biometry. Recent research on evolvability is more than a relabeling of old questions, however. New operational measures of evolvability have opened possibilities for understanding adaptation to rapid environmental change, assessing genetic constraints, and linking micro- and macroevolution.

Fitness consequences of hybridization in a predominantly selfing species: insights into the role of dominance and epistatic incompatibilities

Studying the consequences of hybridization on plant performance is insightful to understand the adaptive potential of populations, notably at local scales. Due to reduced effective recombination, predominantly selfing species are organized in highly homozygous multi-locus-genotypes (or lines) that accumulate genetic differentiation both among- and within-populations. This high level of homozygosity facilitates the dissection of the genetic basis of hybrid performance in highly selfing species, which gives insights into the mechanisms of reproductive isolation between lines. Here, we explored the fitness consequences of hybridization events between natural inbred lines of the predominantly selfing species Medicago truncatula, at both within- and among-populations scales. We found that hybridization has opposite effects pending on studied fitness proxies, with dry mass showing heterosis, and seed production showing outbreeding depression. Although we found significant patterns of heterosis and outbreeding depression, they did not differ significantly for within- compared to among-population crosses. Family-based analyses allowed us to determine that hybrid differentiation was mostly due to dominance and epistasis. Dominance and/or dominant epistatic interactions increased dry mass, while decreasing seed production, and recessive epistatic interactions mostly had a positive effect on both fitness proxies. Our results illustrate how genetic incompatibilities can accumulate at a very local scale among multi-locus-genotypes, and how non-additive genetic effects contribute to heterosis and outbreeding depression.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

Genetics of quantitative traits with dominance under stabilizing and directional selection in partially selfing species

Recurrent self-fertilization is thought to lead to reduced adaptive potential by decreasing the genetic diversity of populations, thus leading selfing lineages down an evolutionary "blind alley." Although well supported theoretically, empirical support for reduced adaptability in selfing species is limited. One limitation of classical theoretical models is that they assume pure additivity of the fitness-related traits that are under stabilizing selection, despite ample evidence that quantitative traits are subject to dominance. Here, we relax this assumption and explore the effect of dominance on a fitness-related trait under stabilizing selection for populations that differ in selfing rates. By decomposing the genetic variance into additional components specific to inbred populations, we show that dominance components can explain a substantial part of the genetic variance of inbred populations. We also show that ignoring these components leads to an upward bias in the predicted response to selection. Finally, we show that when considering the effect of dominance, the short-term evolutionary potential of populations remains comparable across the entire gradient in outcrossing rates, and genetic associations can even make selfing populations more evolvable on the longer term, reconciling theoretical, and empirical results.

The evolution of sex along an environmental gradient

Although temporally changing environments generally favor sex and recombination, the effects of spatial environmental heterogeneity have been less explored. In this article, we use a classical model of adaptation along with an environmental gradient to study the selective forces acting on reproductive mode evolution in the central and marginal parts of the distribution range of a species. The model considers a polygenic trait under stabilizing selection (the optimal trait value changing across space) and includes a demographic component imposing range limits. The results show that in the central part of the range (where populations are well adapted), recombination tends to increase the mean fitness of offspring in regimes where drift is sufficiently strong (generating a benefit for sex), while it has the opposite effect when the effect of drift stays negligible. However, these effects remain weak and are easily overwhelmed by slight intrinsic fitness differences between sexuals and asexuals. In agreement with previous results, asexuality may be favored in marginal populations, as it can preserve adaptation to extreme conditions. However, a substantial advantage of asexuality is possible only in conditions maintaining a strong maladaptation of sexuals at range limits (high effective environmental gradient, weak selection at loci coding for the trait).

Addressing Darwin's dilemma: Can pseudo-overdominance explain persistent inbreeding depression and load?

Darwin spent years investigating the effects of self-fertilization, concluding that "nature abhors perpetual self-fertilization." Given that selection purges inbred populations of strongly deleterious mutations and drift fixes mild mutations, why does inbreeding depression (ID) persist in highly inbred taxa and why do no purely selfing taxa exist? Background selection, associations and interference among loci, and drift within small inbred populations all limit selection while often increasing fixation. These mechanisms help to explain why more inbred populations in most species consistently show more fixed load. This drift load is manifest in the considerable heterosis regularly observed in between-population crosses. Such heterosis results in subsequent high ID, suggesting a mechanism by which small populations could retain variation and inbreeding load. Multiple deleterious recessive mutations linked in repulsion generate pseudo-overdominance. Many tightly linked load loci could generate a balanced segregating load high enough to sustain ID over many generations. Such pseudo-overdominance blocks (or "PODs") are more likely to occur in regions of low recombination. They should also result in clear genetic signatures including genomic hotspots of heterozygosity; distinct haplotypes supporting alleles at intermediate frequency; and high linkage disequilibrium in and around POD regions. Simulation and empirical studies tend to support these predictions. Additional simulations and comparative genomic analyses should explore POD dynamics in greater detail to resolve whether PODs exist in sufficient strength and number to account for why ID and load persist within inbred lineages.

Maintenance of high inbreeding depression in selfing populations: Two-stage effect of early- and late-acting mutations

High estimates of inbreeding depression have been obtained in many plant populations with high selfing rates. However, deleterious mutations might be purged from such populations as a result of selfing. I developed a simulation model assuming the presence of mutations at two sets of loci, namely, early- and late-acting loci, and the selective abortion of embryos coupled with ovule overproduction. In the model, early-acting loci are expressed during embryo initiation, and less vigorous embryos are aborted. Late-acting loci are expressed after selective abortion ends; the surviving embryos (seeds) compete, and some of them form the next generation. If mutations are allowed to occur in both early- and late-acting loci, both types increase in frequency in populations with high selfing rates. However, this phenomenon does not occur if mutations occur only in the early- or only in the late-acting loci. Consistent results are observed even if the total number of loci in which mutations are allowed to occur is the same for simulations with both early- and late-acting loci, only early-acting loci, or only late-acting loci, indicating that the presence of both types of loci is the causal factor. Thus, the two-stage effect, or occurrence of both early- and late-acting mutations, promotes the maintenance of these mutations in populations with high selfing rates.

Identification of the First Oomycete Mating-type Locus Sequence in the Grapevine Downy Mildew Pathogen, Plasmopara viticola

Mating types are self-incompatibility systems that promote outcrossing in plants, fungi, and oomycetes. Mating-type genes have been widely studied in plants and fungi but have yet to be identified in oomycetes, eukaryotic organisms closely related to brown algae that cause many destructive animal and plant diseases. We identified the mating-type locus of Plasmopara viticola, the oomycete responsible for grapevine downy mildew, one of the most damaging grapevine diseases worldwide. Using a genome-wide association approach, we identified a 570-kb repeat-rich non-recombining region controlling mating types, with two highly divergent alleles. We showed that one mating type was homozygous, whereas the other was heterozygous at this locus. The mating-type locus encompassed 40 genes, including one encoding a putative hormone receptor. Functional studies will, however, be required to validate the function of these genes and find the actual determinants of mating type. Our findings have fundamental implications for our understanding of the evolution of mating types, as they reveal a unique determinism involving an asymmetry of heterozygosity, as in sex chromosomes and unlike other mating-type systems. This identification of the mating-type locus in such an economically important crop pathogen also has applied implications, as outcrossing facilitates rapid evolution and resistance to harsh environmental conditions.

Hidden genetic variance contributes to increase the short-term adaptive potential of selfing populations

Standing genetic variation is considered a major contributor to the adaptive potential of species. The low heritable genetic variation observed in self-fertilizing populations has led to the hypothesis that species with this mating system would be less likely to adapt. However, a non-negligible amount of cryptic genetic variation for polygenic traits, accumulated through negative linkage disequilibrium, could prove to be an important source of standing variation in self-fertilizing species. To test this hypothesis, we simulated populations under stabilizing selection subjected to an environmental change. We demonstrate that, when the mutation rate is high (but realistic), selfing populations are better able to store genetic variance than outcrossing populations through genetic associations, notably due to the reduced effective recombination rate associated with predominant selfing. Following an environmental shift, this diversity can be partially remobilized, which increases the additive variance and adaptive potential of predominantly (but not completely) selfing populations. In such conditions, despite initially lower observed genetic variance, selfing populations adapt as readily as outcrossing ones within a few generations. For low mutation rates, purifying selection impedes the storage of diversity through genetic associations, in which case, as previously predicted, the lower genetic variance of selfing populations results in lower adaptability compared to their outcrossing counterparts. The population size and the mutation rate are the main parameters to consider, as they are the best predictors of the amount of stored diversity in selfing populations. Our results and their impact on our knowledge of adaptation under high selfing rates are discussed.

Epistasis, inbreeding depression, and the evolution of self-fertilization

Inbreeding depression resulting from partially recessive deleterious alleles is thought to be the main genetic factor preventing self-fertilizing mutants from spreading in outcrossing hermaphroditic populations. However, deleterious alleles may also generate an advantage to selfers in terms of more efficient purging, while the effects of epistasis among those alleles on inbreeding depression and mating system evolution remain little explored. In this article, we use a general model of selection to disentangle the effects of different forms of epistasis (additive-by-additive, additive-by-dominance, and dominance-by-dominance) on inbreeding depression and on the strength of selection for selfing. Models with fixed epistasis across loci, and models of stabilizing selection acting on quantitative traits (generating distributions of epistasis) are considered as special cases. Besides its effects on inbreeding depression, epistasis may increase the purging advantage associated with selfing (when it is negative on average), while the variance in epistasis favors selfing through the generation of linkage disequilibria that increase mean fitness. Approximations for the strengths of these effects are derived, and compared with individual-based simulation results.

Bayesian inference of natural selection from spatiotemporal phenotypic data

Spatiotemporal variations of natural selection may influence the evolution of various features of organisms such as local adaptation or specialisation. This article develops a method for inferring how selection varies between locations and between generations from phenotypic data. It is assumed that generations are non-overlapping and that individuals reproduce by selfing or asexually. A quantitative genetics model taking account of the effects of stabilising natural selection, the environment and mutation on phenotypic means and variances is developed. Explicit results on the evolution of populations are derived and used to develop a Bayesian inference method. The latter is applied to simulated data and to data from a wheat participatory plant breeding programme. It has some ability to infer evolutionary parameters, but estimates may be sensitive to prior distributions, for example when phenotypic time series are short and when environmental effects are large. In such cases, sensitivity to prior distributions may be reported or more data may be collected.

Inbreeding depression is difficult to purge in self-incompatible populations of Leavenworthia alabamica

The extent to which inbreeding depression can be purged is a major determinant of mating system evolution and is important to conservation and crop improvement. Studies of inbreeding depression purging have not been conducted in self-incompatible plants before. An experimental ('ancestral') treatment was first created from self-incompatible plants of Leavenworthia alabamica. Lines derived from this population were maintained by self-pollination for three generations in the attempt to create a 'purged' population with fewer recessive, deleterious mutations of large effect. Fitness components and the frequency of malformed phenotypes were monitored in progeny derived from selfing and outcrossing in the ancestral and purged treatments. Fitness component means and inbreeding depression were largely unchanged by three generations of forced self-pollination, and there was no reduction in the frequency of plants exhibiting malformed phenotypes. Our findings indicate that inbreeding depression in this species is largely a result of mutations of mild effect, consistent with the observation that self-incompatibility is maintained in most populations of L. alabamica, despite the presence of genetic variants with weaker self-incompatibility. Moreover, although population theory suggests that deleterious mutations of large effect should be sheltered from selection in the region of self-incompatibility locus, our results do not support this prediction.

How relatedness between mates influences reproductive success: An experimental analysis of self-fertilization and biparental inbreeding in a marine bryozoan

Kin associations increase the potential for inbreeding. The potential for inbreeding does not, however, make inbreeding inevitable. Numerous factors influence whether inbreeding preference, avoidance, or tolerance evolves, and, in hermaphrodites where both self-fertilization and biparental inbreeding are possible, it remains particularly difficult to predict how selection acts on the overall inbreeding strategy, and to distinguish the type of inbreeding when making inferences from genetic markers. Therefore, we undertook an empirical analysis on an understudied type of mating system (spermcast mating in the marine bryozoan, Bugula neritina) that provides numerous opportunities for inbreeding preference, avoidance, and tolerance. We created experimental crosses, containing three generations from two populations to estimate how parental reproductive success varies across parental relatedness, ranging from self, siblings, and nonsiblings from within the same population. We found that the production of viable selfed offspring was extremely rare (only one colony produced three selfed offspring) and biparental inbreeding more common. Paternity analysis using 16 microsatellite markers confirmed outcrossing. The production of juveniles was lower for sib mating compared with nonsib mating. We found little evidence for consistent inbreeding, in terms of nonrandom mating, in adult samples collected from three populations, using multiple population genetic inferences. Our results suggest several testable hypotheses that potentially explain the overall mating and dispersal strategy in this species, including early inbreeding depression, inbreeding avoidance through cryptic mate choice, and differential dispersal distances of sperm and larvae.

Partial Selfing Can Reduce Genetic Loads While Maintaining Diversity During Experimental Evolution

Partial selfing, whereby self- and cross- fertilization occur in populations at intermediate frequencies, is generally thought to be evolutionarily unstable. Yet, it is found in natural populations. This could be explained if populations with partial selfing are able to reduce genetic loads and the possibility for inbreeding depression while keeping genetic diversity that may be important for future adaptation. To address this hypothesis, we compare the experimental evolution of Caenorhabditis elegans populations under partial selfing, exclusive selfing or predominant outcrossing, while they adapt to osmotically challenging conditions. We find that the ancestral genetic load, as measured by the risk of extinction upon inbreeding by selfing, is maintained as long as outcrossing is the main reproductive mode, but becomes reduced otherwise. Analysis of genome-wide single-nucleotide polymorphisms (SNPs) during experimental evolution and among the inbred lines that survived enforced inbreeding indicates that populations with predominant outcrossing or partial selfing maintained more genetic diversity than expected with neutrality or purifying selection. We discuss the conditions under which this could be explained by the presence of recessive deleterious alleles and/or overdominant loci. Taken together, our observations suggest that populations evolving under partial selfing can gain some of the benefits of eliminating unlinked deleterious recessive alleles and also the benefits of maintaining genetic diversity at partially dominant or overdominant loci that become associated due to variance of inbreeding levels.

Effect of partial selfing and polygenic selection on establishment in a new habitat

This article analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the inbreeding history model to characterize mutation‐selection balance in a large, partially selfing source population under selection involving multiple nonidentical loci. I then use individual‐based simulations to study the eco‐evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long‐term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection is discussed.

How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species

Most theoretical works predict that selfing should reduce the level of additive genetic variance available for quantitative traits within natural populations. Despite a growing number of quantitative genetic studies undertaken during the last two decades, this prediction is still not well supported empirically. To resolve this issue and confirm or reject theoretical predictions, we reviewed quantitative trait heritability estimates from natural plant populations with different rates of self-fertilization and carried out a meta-analysis. In accordance with models of polygenic traits under stabilizing selection, we found that the fraction of additive genetic variance is negatively correlated with the selfing rate. Although the mating system explains a moderate fraction of the variance, the mean reduction of narrow-sense heritability values between strictly allogamous and predominantly selfing populations is strong, around 60%. Because some nonadditive components of genetic variance become selectable under inbreeding, we determine whether self-fertilization affects the relative contribution of these components to genetic variance by comparing narrow-sense heritability estimates from outcrossing populations with broad-sense heritability estimated in autogamous populations. Results suggest that these nonadditive components of variance may restore some genetic variance in predominantly selfing populations; it remains, however, uncertain how these nonadditive components will contribute to adaptation.

Structure of multilocus genetic diversity in predominantly selfing populations

Predominantly selfing populations are expected to have reduced effective population sizes due to nonrandom sampling of gametes, demographic stochasticity (bottlenecks or extinction-recolonization), and large scale hitchhiking (reduced effective recombination). Thus, they are expected to display low genetic diversity, which was confirmed by empirical studies. The structure of genetic diversity in predominantly selfing species is dramatically different from outcrossing ones, with populations often dominated by one or a few multilocus genotypes (MLGs) coexisting with several rare genotypes. Therefore, multilocus diversity indices are relevant to describe diversity in selfing populations. Here, we use simulations to provide analytical expectations for multilocus indices and examine whether selfing alone can be responsible for the high-frequency MLGs persistent through time in the absence of selection. We then examine how combining single and multilocus indices of diversity may be insightful to distinguish the effects of selfing, population size, and more complex demographic events (bottlenecks, migration, admixture, or extinction-recolonization). Finally, we examine how temporal changes in MLG frequencies can be insightful to understand the evolutionary trajectory of a given population. We show that combinations of selfing and small demographic sizes can result in high-frequency MLGs, as observed in natural populations. We also show how different demographic scenarios can be distinguished by the parallel analysis of single and multilocus indices of diversity, and we emphasize the importance of temporal data for the study of predominantly selfing populations. Finally, the comparison of our simulations with empirical data on populations of Medicago truncatula confirms the pertinence of our simulation framework.

Sexual selection and inbreeding: Two efficient ways to limit the accumulation of deleterious mutations

Theory and empirical data showed that two processes can boost selection against deleterious mutations, thus facilitating the purging of the mutation load: inbreeding, by exposing recessive deleterious alleles to selection in homozygous form, and sexual selection, by enhancing the relative reproductive success of males with small mutation loads. These processes tend to be mutually exclusive because sexual selection is reduced under mating systems that promote inbreeding, such as self‐fertilization in hermaphrodites. We estimated the relative efficiency of inbreeding and sexual selection at purging the genetic load, using 50 generations of experimental evolution, in a hermaphroditic snail (Physa acuta). To this end, we generated lines that were exposed to various intensities of inbreeding, sexual selection (on the male function), and nonsexual selection (on the female function). We measured how these regimes affected the mutation load, quantified through the survival of outcrossed and selfed juveniles. We found that juvenile survival strongly decreased in outbred lines with reduced male selection, but not when female selection was relaxed, showing that male‐specific sexual selection does purge deleterious mutations. However, in lines exposed to inbreeding, where sexual selection was also relaxed, survival did not decrease, and even increased for self‐fertilized juveniles, showing that purging through inbreeding can compensate for the absence of sexual selection. Our results point to the further question of whether a mixed strategy combining the advantages of both mechanisms of genetic purging could be evolutionary stable.

Heterogeneity in local density allows a positive evolutionary relationship between self-fertilisation and dispersal

Despite empirical evidence for a positive relationship between dispersal and self-fertilization (selfing), theoretical work predicts that these traits should always be negatively correlated, and the Good Coloniser Syndrome of high dispersal and selfing (Cf. Baker's Law) should not evolve. Critically, previous work assumes that adult density is spatiotemporally homogeneous, so selfing results in identical offspring production for all patches, eliminating the benefit of dispersal for escaping from local resource competition. We investigate the joint evolution of dispersal and selfing in a demographically structured metapopulation model where local density is spatiotemporally heterogeneous due to extinction-recolonization dynamics. Selfing alleviates outcrossing failure due to low local density (an Allee effect) while dispersal alleviates competition through dispersal of propagules from high- to low-density patches. Because local density is spatiotemporally heterogeneous in our model, selfing does not eliminate heterogeneity in competition, so dispersal remains beneficial even under full selfing. Hence the Good Coloniser Syndrome is evolutionarily stable under a broad range of conditions, and both negative and positive relationships between dispersal and selfing are possible, depending on the environment. Our model thus accommodates positive empirical relationships between dispersal and selfing not predicted by previous theoretical work and provides additional explanations for negative relationships.

The Evolutionary Interplay between Adaptation and Self-Fertilization

Genome-wide surveys of nucleotide polymorphisms, obtained from next-generation sequencing, have uncovered numerous examples of adaptation in self-fertilizing organisms, especially regarding changes to climate, geography, and reproductive systems. Yet existing models for inferring attributes of adaptive mutations often assume idealized outcrossing populations, which risks mischaracterizing properties of these variants. Recent theoretical work is emphasizing how various aspects of self-fertilization affects adaptation, yet empirical data on these properties are lacking. We review theoretical and empirical studies demonstrating how self-fertilization alters the process of adaptation, illustrated using examples from current sequencing projects. We propose ideas for how future research can more accurately quantify aspects of adaptation in self-fertilizers, including incorporating the effects of standing variation, demographic history, and polygenic adaptation.

Analysis of large-scale next-generation sequencing datasets are finding more examples of adaptive evolution at the genomic level.

Advances in theoretical work has demonstrated how self-fertilisation affects different aspects of adaptation in these organisms, compared to outcrossers.

Current software and statistical methods do not take different mating systems into account, which risks mischaracterising the presence or strength of adaptive mutations from genome scans.

Development of new mathematical and statistical methods that explicitly consider self-fertilization and associated demographic effects will enable researchers to more accurately quantify adaptation in these organisms.

What can genome-wide association studies tell us about the evolutionary forces maintaining genetic variation for quantitative traits?

Contents 21 I. 21 II. 22 III. 24 IV. 25 V. 29 30 References 30 SUMMARY: Understanding the evolutionary forces that shape genetic variation within species has long been a goal of evolutionary biology. Integrating data for the genetic architecture of traits from genome-wide association mapping studies (GWAS) along with the development of new population genetic methods for identifying selection in sequence data may allow us to evaluate the roles of mutation-selection balance and balancing selection in shaping genetic variation at various scales. Here, we review the theoretical predictions for genetic architecture and additional signals of selection on genomic sequence for the loci that affect traits. Next, we review how plant GWAS have tested for the signatures of various selective scenarios. Limited evidence to date suggests that within-population variation is maintained primarily by mutation-selection balance while variation across the landscape is the result of local adaptation. However, there are a number of inherent biases in these interpretations. We highlight these challenges and suggest ways forward to further understanding of the maintenance of variation.

Inbreeding depression under mixed outcrossing, self-fertilization and sib-mating

BACKGROUND

Biparental inbreeding, mating between two relatives, occurs at a low frequency in many natural plant populations, which also often have substantial rates of self-fertilization. Although biparental inbreeding is likely to influence the dynamics of inbreeding depression and the evolution of selfing rates, it has received limited theoretical attention in comparison to selfing. The only previous model suggested that biparental inbreeding can favour the maintenance of stable intermediate selfing rates, but made unrealistic assumptions about the genetic basis of inbreeding depression. Here we extend a genetic model of inbreeding depression, describing nearly recessive lethal mutations at a very large number of loci, to incorporate sib-mating. We also include a constant component of inbreeding depression modelling the effects of mildly deleterious, nearly additive alleles. We analyze how observed rates of sib-mating influence the mean number of heterozygous lethals alleles and inbreeding depression in a population reproducing by a mixture of self-fertilization, sib-mating and outcrossing. We finally use the ensuing relationship between equilibrium inbreeding depression and population selfing rate to infer the evolutionarily stable selfing rates expected under such a mixed mating system.

RESULTS

We show that for a given rate of inbreeding, sib-mating is more efficient at purging inbreeding depression than selfing, because homozygosity of lethals increases more gradually through sib-mating than through selfing. Because sib-mating promotes the purging of inbreeding depression and the evolution of selfing, our genetic model of inbreeding depression also predicts that sib-mating is unlikely to maintain stable intermediate selfing rates.

CONCLUSIONS

Our results imply that even low rates of sib-mating affect plant mating system evolution, by facilitating the evolution of selfing via more efficient purging of inbreeding depression. Alternative mechanisms, such as pollination ecology, are necessary to explain stable mixed selfing and outcrossing.

Limits to Adaptation in Partially Selfing Species

In outcrossing populations, "Haldane's sieve" states that recessive beneficial alleles are less likely to fix than dominant ones, because they are less exposed to selection when rare. In contrast, selfing organisms are not subject to Haldane's sieve and are more likely to fix recessive types than outcrossers, as selfing rapidly creates homozygotes, increasing overall selection acting on mutations. However, longer homozygous tracts in selfers also reduce the ability of recombination to create new genotypes. It is unclear how these two effects influence overall adaptation rates in partially selfing organisms. Here, we calculate the fixation probability of beneficial alleles if there is an existing selective sweep in the population. We consider both the potential loss of the second beneficial mutation if it has a weaker advantage than the first one, and the possible replacement of the initial allele if the second mutant is fitter. Overall, loss of weaker adaptive alleles during a first selective sweep has a larger impact on preventing fixation of both mutations in highly selfing organisms. Furthermore, the presence of linked mutations has two opposing effects on Haldane's sieve. First, recessive mutants are disproportionally likely to be lost in outcrossers, so it is likelier that dominant mutations will fix. Second, with elevated rates of adaptive mutation, selective interference annuls the advantage in selfing organisms of not suffering from Haldane's sieve; outcrossing organisms are more able to fix weak beneficial mutations of any dominance value. Overall, weakened recombination effects can greatly limit adaptation in selfing organisms.

Lower selfing rates in metallicolous populations than in non-metallicolous populations of the pseudometallophyte Noccaea caerulescens (Brassicaceae) in Southern France

BACKGROUND AND AIMS

The pseudometallophyte Noccaea caerulescens is an excellent model to study evolutionary processes, as it grows both on normal and on heavy-metal-rich, toxic soils. The evolution and demography of populations are critically impacted by mating system and, yet, information about the N. caerulescens mating system is limited.

METHODS

Mean selfing rates were assessed using microsatellite loci and a robust estimation method (RMES) in five metallicolous and five non-metallicolous populations of N. caerulescens in Southern France, and this measure was replicated for two successive reproductive seasons. As a part of the study, the patterns of gene flow among populations were analysed. The mating system was then characterized at a fine spatial scale in three populations using the MLTR method on progeny arrays.

KEY RESULTS

The results confirm that N. caerulescens has a mixed mating system, with selfing rates ranging from 0·2 to 0·5. Selfing rates did not vary much among populations within ecotypes, but were lower in the metallicolous than in the non-metallicolous ecotype, in both seasons. Effective population size was also lower in non-metallicolous populations. Biparental inbreeding was null to moderate. Differentiation among populations was generally high, but neither ecotype nor isolation by distance explained it.

CONCLUSIONS

The consequences of higher selfing rates on adaptation are expected to be weak to moderate in non-metallicolous populations and they are expected to suffer less from inbreeding depression, compared to metallicolous populations.

Increased heterosis in selfing populations of a perennial forb

Quantifying the importance of random genetic drift in natural populations is central to understanding the potential limits to natural selection. One approach is to estimate the magnitude of heterosis, the increased fitness of progeny derived from crosses between populations relative to crosses within populations caused by the heterozygous masking of deleterious recessive or nearly recessive alleles that have been fixed by drift within populations. Self-fertilization is expected to reduce the effective population size by half relative to outcrossing, and population bottlenecks may be common during the transition to selfing. Therefore, chance fixation of deleterious alleles due to drift in selfing populations should increase heterosis between populations. Increased homozygosity due to fixation or loss of alleles should also decrease inbreeding depression within populations. Most populations of the perennial herb Arabidopsis lyrata ssp. lyrata are self-incompatible (SI), but several have evolved self-compatibility and are highly selfing. We quantified heterosis and inbreeding depression in two predominantly self-compatible (SC) and seven SI populations in a field common garden experiment within the species' native range and examined the correlation between these metrics to gauge the similarity in their genetic basis. We measured proportion germination in the lab, and survival and fecundity (flower and seed production) for 2 years in the field, and calculated estimates of cumulative fitness. We found 7.2-fold greater heterosis in SC compared with SI populations, despite substantial heterosis in SI populations (56 %). Inbreeding depression was >61 %, and not significantly different between SC and SI populations. There was no correlation between population estimates of heterosis and inbreeding depression, suggesting that they have somewhat different genetic bases. Combined with other sources of information, our results suggest a history of bottlenecks in all of these populations. The bottlenecks in SC populations may have been severe, but their strong inbreeding depression remains enigmatic.