Limits to Adaptation in Partially Selfing Species

Pervasive G × E interactions shape adaptive trajectories and the exploration of the phenotypic space in artificial selection experiments

Quantitative genetics models have shown that long-term selection responses depend on initial variance and mutational influx. Understanding limits of selection requires quantifying the role of mutational variance. However, correlative responses to selection on nonfocal traits can perturb the selection response on the focal trait; and generations are often confounded with selection environments so that genotype by environment (G×E) interactions are ignored. The Saclay divergent selection experiments (DSEs) on maize flowering time were used to track the fate of individual mutations combining genotyping data and phenotyping data from yearly measurements (DSEYM) and common garden experiments (DSECG) with four objectives: (1) to quantify the relative contribution of standing and mutational variance to the selection response, (2) to estimate genotypic mutation effects, (3) to study the impact of G×E interactions in the selection response, and (4) to analyze how trait correlations modulate the exploration of the phenotypic space. We validated experimentally the expected enrichment of fixed beneficial mutations with an average effect of +0.278 and +0.299 days to flowering, depending on the genetic background. Fixation of unfavorable mutations reached up to 25% of incoming mutations, a genetic load possibly due to antagonistic pleiotropy, whereby mutations fixed in the selection environment (DSEYM) turned to be unfavorable in the evaluation environment (DSECG). Global patterns of trait correlations were conserved across genetic backgrounds but exhibited temporal patterns. Traits weakly or uncorrelated with flowering time triggered stochastic exploration of the phenotypic space, owing to microenvironment-specific fixation of standing variants and pleiotropic mutational input.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

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.

Population Rescue through an Increase in the Selfing Rate under Pollen Limitation: Plasticity versus Evolution

AbstractIncreased rates of self-fertilization offer reproductive assurance when plant populations experience pollen limitation, but self-fertilization may reduce fitness by exposing deleterious mutations. If an environmental change responsible for pollen limitation also induces plastic mating system shifts toward self-pollination, the reproductive assurance benefit and inbreeding depression cost of increased self-fertilization occur immediately, while the benefit and cost happen more gradually when increased self-fertilization occur through evolution. I built eco-evolutionary models to explore the demographic and genetic conditions in which higher self-fertilization by plasticity and/or evolution rescues populations, following deficits due to a sudden onset of pollen limitation. Rescue is most likely under an intermediate level of selfing rate increase, either through plasticity or evolution, and this critical level of selfing rate increase is higher under stronger pollen limitation. Generally, rescue is more likely through plasticity than through evolution. Under weak pollen limitation, rescue by enhanced self-fertilization may mainly occur through purging of deleterious mutations rather than reproductive assurance. The selfing rate increase conferring the highest rescue probability is lower when the initial population size is smaller. This article shows the importance of plasticity during plant population rescue and offers insights for future studies of the evolution of mating system plasticity.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Drivers of strong isolation and small effective population size at a leading range edge of a widespread plant

Climate change has influenced species distributions worldwide with upward elevational shifts observed in many systems. Leading range edge populations, like those at upper elevation limits, are crucial for climate change responses but can exhibit low genetic diversity due to founder effects, isolation, or limited outbreeding. These factors can hamper local adaptation at range limits. Using the widespread herb, Argentina anserina, we measured ecological attributes (population density on the landscape, area of population occupancy, and plant and flower density) spanning a 1000 m elevation gradient, with high elevation populations at the range limit. We measured vegetative clonal potential in the greenhouse for populations spanning the gradient. We combined these data with a ddRAD-seq dataset to test the hypotheses that high elevation populations would exhibit ecological and genomic signatures of leading range edge populations. We found that population density on the landscape declined towards the high elevation limit, as is expected towards range edges. However, plant density was elevated within edge populations. In the greenhouse, high elevation plants exhibited stronger clonal potential than low elevation plants, likely explaining increased plant density in the field. Phylogeographic analysis supported more recent colonization of high elevation populations which were also more genetically isolated, had more extreme heterozygote excess and had smaller effective population size than low. Results support that colonization of high elevations was likely accompanied by increased asexuality, contributing to a decline in effective population size. Despite high plant density in leading edge populations, their small effective size, isolation and clonality could constrain adaptive potential.

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.

The genetic basis of selfing rate evolution

Evolution of selfing is common in plant populations, but the genetic basis of selfing rate evolution remains unclear. Although the effects of genetic properties on fixation for mating-unrelated alleles have been investigated, loci that modify the selfing rate (selfing modifiers) differ from mating-unrelated loci in several aspects. Using population genetic models, I investigate the genetic basis of selfing rate evolution. For mating-unrelated alleles, selfing promotes fixation only for recessive mutations, but for selfing modifiers, because the selection coefficient depends on the background selfing rate, selfing can promote fixation even for dominant modifiers. For mating-unrelated alleles, the fixation probability from standing variation is independent of dominance and decreases with an increased background selfing rate. However, for selfing modifiers, the fixation probability peaks at an intermediate selfing rate and when alleles are recessive, because a change of its selection coefficient necessarily involves a change of the inbreeding coefficient, because both depend on the level of inbreeding depression. Furthermore, evolution of selfing involving multiple modifier loci is more likely when selfing is controlled by few large-effect rather than many slight-effect modifiers. I discuss how these characteristics of selfing modifiers have implications for the unidirectional transition from outcrossing to selfing and other empirical patterns.

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.

Establishment of local adaptation in partly self-fertilizing populations

Populations often inhabit multiple ecological patches and thus experience divergent selection, which can lead to local adaptation if migration is not strong enough to swamp locally adapted alleles. Conditions for the establishment of a locally advantageous allele have been studied in randomly mating populations. However, many species reproduce, at least partially, through self-fertilization, and how selfing affects local adaptation remains unclear and debated. Using a two-patch branching process formalism, we obtained a closed-form approximation under weak selection for the probability of establishment of a locally advantageous allele (P) for arbitrary selfing rate and dominance level, where selection is allowed to act on viability or fecundity, and migration can occur via seed or pollen dispersal. This solution is compared to diffusion approximation and used to investigate the consequences of a shift in a mating system on P, and the establishment of protected polymorphism. We find that selfing can either increase or decrease P, depending on the patterns of dominance in the two patches, and has conflicting effects on local adaptation. Globally, selfing favors local adaptation when locally advantageous alleles are (partially) recessive, when selection between patches is asymmetrical and when migration occurs through pollen rather than seed dispersal. These results establish a rigorous theoretical background to study heterogeneous selection and local adaptation in partially selfing species.

Balancing selection in self-fertilizing populations

Self-fertilization commonly occurs in hermaphroditic species, either occasionally or as the main reproductive mode. It strongly affects the genetic functioning of a population by increasing homozygosity and genetic drift and reducing the effectiveness of recombination. Balancing selection is a form of selection that maintains polymorphism, which has been extensively studied in outcrossing species. Yet, despite recent developments, the analysis of balancing selection in partially selfing species is limited to specific cases and a general treatment is still lacking. In particular, it is unclear whether selfing globally reduced the efficacy of balancing selection as in the well-known case of overdominance. I provide a unifying framework, quantify how selfing affects the maintenance of polymorphism and the efficacy of the different form of balancing selection, and show that they can be classified into two main categories: overdominance-like selection (including true overdominance, selection variable in space and time, and antagonistic selection), which is strongly affected by selfing, and negative frequency dependent selection, which is barely affected by selfing, even at multiple loci. I also provide simple analytical results for all cases under the assumption of weak selection. This framework provides theoretical background to analyze the genomic signature of balancing selection in partially selfing species. It also sheds new light on the evolution of selfing species, including the evolution of selfing syndrome, the interaction with pathogens, and the evolutionary fate of selfing lineages.

Inferring the Demographic History of Inbred Species from Genome-Wide SNP Frequency Data

Demographic inference using the site frequency spectrum (SFS) is a common way to understand historical events affecting genetic variation. However, most methods for estimating demography from the SFS assume random mating within populations, precluding these types of analyses in inbred populations. To address this issue, we developed a model for the expected SFS that includes inbreeding by parameterizing individual genotypes using beta-binomial distributions. We then take the convolution of these genotype probabilities to calculate the expected frequency of biallelic variants in the population. Using simulations, we evaluated the model's ability to coestimate demography and inbreeding using one- and two-population models across a range of inbreeding levels. We also applied our method to two empirical examples, American pumas (Puma concolor) and domesticated cabbage (Brassica oleracea var. capitata), inferring models both with and without inbreeding to compare parameter estimates and model fit. Our simulations showed that we are able to accurately coestimate demographic parameters and inbreeding even for highly inbred populations (F = 0.9). In contrast, failing to include inbreeding generally resulted in inaccurate parameter estimates in simulated data and led to poor model fit in our empirical analyses. These results show that inbreeding can have a strong effect on demographic inference, a pattern that was especially noticeable for parameters involving changes in population size. Given the importance of these estimates for informing practices in conservation, agriculture, and elsewhere, our method provides an important advancement for accurately estimating the demographic histories of these species.

Environmental fluctuations dampen the effects of clonal reproduction on evolutionary rescue

Evolutionary rescue occurs when genetic change allows a population to persist in response to an environmental change that would otherwise have led to extinction. Most studies of evolutionary rescue assume that species have either fully clonal or fully sexual reproduction; however, many species have partially clonal reproductive strategies in which they reproduce both clonally and sexually. Furthermore, the few evolutionary rescue studies that have evaluated partially clonal reproduction did not consider fluctuations in the environment, which are nearly ubiquitous in nature. Here, we use individual-based simulations to investigate how environmental fluctuations (either uncorrelated or positively autocorrelated) influence the effect of clonality on evolutionary rescue. We show that, for moderate magnitudes of environmental fluctuations, as was found in the absence of fluctuations, increasing the degree of clonality increases the probability of population persistence in response to an abrupt environmental change, but decreases persistence in response to a continuous, directional environmental change. However, with large magnitudes of fluctuations, both the benefits of clonality following a step change and the detrimental effects of clonality following a continuous, directional change are generally reduced; in fact, in the latter scenario, increasing clonality can even become beneficial if environmental fluctuations are autocorrelated. We also show that increased generational overlap dampens the effects of environmental fluctuations. Overall, we demonstrate that understanding the evolutionary rescue of partially clonal organisms requires not only knowledge of the species life history and the type of environmental change, but also an understanding of the magnitude and autocorrelation of environmental fluctuations.

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.

How and When Does Outcrossing Occur in the Predominantly Selfing Species Medicago truncatula?

Empirical studies on natural populations of Medicago truncatula revealed selfing rates higher than 80%, but never up to 100%. Similarly, several studies of predominantly selfing species show variability in the level of residual outcrossing between populations and also between temporal samples of the same population. However, these studies measure global selfing rates at the scale of the population and we do not know whether there is intra-population variation and how outcrossing events are distributed, between genotypes, plants, flowers, or seeds. Theoretical studies predict the maintenance of residual outcrossing in highly selfing species due to environmental (e.g., pollen biology) and/or genetic determinants and decompositions of the variation in outcrossing rate using experimental data can be very informative to test these hypotheses. Here, we focus on one natural population of M. truncatula in order to describe precisely its mating system. In particular, we investigated the determinants of the selfing rate by testing for seasonal variations (environmental determinism) and variations between genotypes (genetic determinism). We measured selfing rates in maternal progenies from plants collected widely across a natural population. For each plant, we collected pods from flowers produced at the beginning and at the end of the flowering season to test for a seasonal variation in the outcrossing rate. For each collected offspring, we also estimated the likelihood that it was issued from a self-fertilization event and assessed the genetic component of variation of this mating system measure. We found a significant, albeit small, increase in outcrossing rate in progenies collected at the end [t _m = 0.137 (SD = 0.025)] compared to those collected at the beginning [t _m = 0.083 (0.016)] of the flowering season. A significant between genotypes variation in selfing rate was also detected, resulting in a heritability of 9% for the rate of residual outcrossing. Altogether, our work shows that despite a predominantly selfing reproductive mode, M. truncatula displays variation in residual outcrossing rate, and that this trait is likely under a complex determinism combining environmental and genetic factors. We discuss the evolutionary implications of our results for the population.

The opportunity for outcrossing varies across the geographic range of the primarily selfing Clarkia xantiana ssp. parviflora

PREMISE

The timing of self-fertilization has potentially important consequences for the trajectory of mating system evolution, the opportunity for outcrossing, and the maintenance of genetic variation in populations. For primarily selfing taxa, it remains poorly understood as to how floral variation influences the opportunity for outcrossing and whether those attributes vary among populations across geographic ranges.

METHODS

We examined variation in floral traits (herkogamy, protandry, flower size, stigma stage at anthesis, timing of stigma receptivity) in seven populations of Clarkia xantiana ssp. parviflora, a primarily selfing taxon, spanning from the western to eastern margins of its distribution. We also performed experimental emasculations and pollinations (followed by stigma severing) to quantify the extent of opportunities for outcrossing across flower development.

RESULTS

There was marked among-population variation in all floral traits, particularly between far eastern and western populations. Emasculation experiments showed that the eastern populations had minimal autonomous selfing, but western populations had high rates of selfing within 24 h after anthesis. Population variation in autofertility was significantly predicted by floral trait variation, especially protandry and petal size.

CONCLUSIONS

Greater protandry both extended the time over which outcrossing could potentially occur and reduced the probability of autonomous selfing, suggesting that there may be a tradeoff that results in fitness loss when pollinator visits are not common. The east-west pattern of differentiation in some floral traits parallels that of postglacial range expansion, suggesting that selection on the mating system may have been strong in the process of range expansion.

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.

Evolution of self-incompatibility in the Brassicaceae: Lessons from a textbook example of natural selection

Self-incompatibility (SI) is a self-recognition genetic system enforcing outcrossing in hermaphroditic flowering plants and results in one of the arguably best understood forms of natural (balancing) selection maintaining genetic variation over long evolutionary times. A rich theoretical and empirical population genetics literature has considerably clarified how the distribution of SI phenotypes translates into fitness differences among individuals by a combination of inbreeding avoidance and rare-allele advantage. At the same time, the molecular mechanisms by which self-pollen is specifically recognized and rejected have been described in exquisite details in several model organisms, such that the genotype-to-phenotype map is also pretty well understood, notably in the Brassicaceae. Here, we review recent advances in these two fronts and illustrate how the joint availability of detailed characterization of genotype-to-phenotype and phenotype-to-fitness maps on a single genetic system (plant self-incompatibility) provides the opportunity to understand the evolutionary process in a unique perspective, bringing novel insight on general questions about the emergence, maintenance, and diversification of a complex genetic system.

Selective Sweeps Under Dominance and Inbreeding

A major research goal in evolutionary genetics is to uncover loci experiencing positive selection. One approach involves finding 'selective sweeps' patterns, which can either be 'hard sweeps' formed by de novo mutation, or 'soft sweeps' arising from recurrent mutation or existing standing variation. Existing theory generally assumes outcrossing populations, and it is unclear how dominance affects soft sweeps. We consider how arbitrary dominance and inbreeding via self-fertilization affect hard and soft sweep signatures. With increased self-fertilization, they are maintained over longer map distances due to reduced effective recombination and faster beneficial allele fixation times. Dominance can affect sweep patterns in outcrossers if the derived variant originates from either a single novel allele, or from recurrent mutation. These models highlight the challenges in distinguishing hard and soft sweeps, and propose methods to differentiate between scenarios.

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.

Mating system impacts the genetic architecture of adaptation to heterogeneous environments

Self-fertilisation has consequences for variation across the genome as it reduces effective population size, effect recombination rates and pollen flow, with implications for local adaptation. We conducted simulations of divergent stabilising selection on a quantitative trait with drift, pollen flow, mutation, recombination and different outcrossing rates. We quantified trait divergence and the genetic architecture of adaptation. We conducted an F_ST outlier analysis to identify candidate loci and quantified the impact of mating system on detectability. Selfing promoted trait divergence mainly through reductions in pollen flow. Moreover, trait architecture became more diffuse with selfing. Average effect size of trait loci was lower, while the number of loci, and their clustering distance increased. The genetic architecture of selfers was also more diffuse than outcrossers for equivalent migration rates. However, when deleterious alleles were included, architectures became more concentrated in selfers, likely to be because of reductions in population size caused by mutational meltdown and impacts of background selection on N_e . Our simulations demonstrate that mating system has important impacts on adaptive divergence of traits and the genetic landscape underlying that divergence. Selfing has a significant effect on detectability of regions of the genome important for adaptation because of neutral divergence and diffuse trait architecture.

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.

Evolution: Selfing Takes Species Down Stebbins's Blind Alley

Shifts from outcrossing to selfing have occurred thousands of times across the tree of life. By reducing the size of the gene pool, selfing should limit adaptive potential. A refreshing empirical experiment with snails supports this long-standing hypothesis.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

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.

The interaction between sex-specific selection and local adaptation in species without separate sexes

Local adaptation in hermaphrodite species can be based on a variety of fitness components, including survival, as well as both female and male sex-functions within individuals. When selection via female and male fitness components varies spatially (e.g. due to environmental heterogeneity), local adaptation will depend, in part, on variation in selection through each fitness component, and the extent to which genetic trade-offs between sex-functions maintain genetic variation necessary for adaptation. Local adaptation will also depend on the hermaphrodite mating system because self-fertilization alters several key factors influencing selection and the maintenance of genetic variance underlying trade-offs between the sex-functions (sexually antagonistic polymorphism). As a first step to guide intuition regarding sex-specific adaptation in hermaphrodites, we develop a simple theoretical model incorporating the essential features of hermaphrodite mating and adaptation in a spatially heterogeneous environment, and explore the interaction between sex-specific selection, self-fertilization and local adaptation. Our results suggest that opportunities for sex-specific local adaptation in hermaphrodites depend strongly on the extent of self-fertilization and inbreeding depression. Using our model as a conceptual framework, we provide a broad overview of the literature on sex-specific selection and local adaptation in hermaphroditic plants and animals, emphasizing promising future directions in light of our theoretical predictions.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Mutation Rate Evolution in Partially Selfing and Partially Asexual Organisms

Different factors can influence the evolution of the mutation rate of a species: costs associated with DNA replication fidelity, indirect selection caused by the mutations produced (that should generally favor lower mutation rates, given that most mutations affecting fitness are deleterious), and genetic drift, which may render selection acting on weak mutators inefficient. In this paper, we use a two-locus model to compute the strength of indirect selection acting on a modifier locus that affects the mutation rate toward a deleterious allele at a second, linked, locus, in a population undergoing partial selfing or partial clonality. The results show that uniparental reproduction increases the effect of indirect selection for lower mutation rates. Extrapolating to the case of a whole genome with many deleterious alleles, and introducing a direct cost to DNA replication fidelity, the results can be used to compute the evolutionarily stable mutation rate, U In the absence of mutational bias toward higher U, the analytical prediction fits well with individual-based, multilocus simulation results. When such a bias is added into the simulations, however, genetic drift may lead to the maintenance of higher mutation rates, and this effect may be amplified in highly selfing or highly clonal populations due to their reduced effective population size.

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.