Evolution of mating types in finite populations: The precarious advantage of being rare

Sexually reproducing populations with self‐incompatibility bear the cost of limiting potential mates to individuals of a different type. Rare mating types escape this cost since they are unlikely to encounter incompatible partners, leading to the deterministic prediction of continuous invasion by new mutants and an ever‐increasing number of types. However, rare types are also at an increased risk of being lost by random drift. Calculating the number of mating types that a population can maintain requires consideration of both the deterministic advantages and the stochastic risks. By comparing the relative importance of selection and drift, we show that a population of size N can maintain a maximum of approximately N^1/3 mating types for intermediate population sizes, whereas for large N, we derive a formal estimate. Although the number of mating types in a population is quite stable, the rare‐type advantage promotes turnover of types. We derive explicit formulas for both the invasion and turnover probabilities in finite populations.     

The divergence history of the perennial plant Linaria cavanillesii confirms a recent loss of self‐incompatibility

Many angiosperms prevent inbreeding through a self‐incompatibility (SI) system, but the loss of SI has been frequent in their evolutionary history. The loss of SI may often lead to an increase in the selfing rate, with the purging of inbreeding depression and the ultimate evolution of a selfing syndrome, where plants have smaller flowers with reduced pollen and nectar production. In this study, we used approximate Bayesian computation (ABC) to estimate the timing of divergence between populations of the plant Linaria cavanillesii that differ in SI status and in which SI is associated with low inbreeding depression but not with a transition to full selfing or a selfing syndrome. Our analysis suggests that the mixed‐mating self‐compatible (SC) population may have begun to diverge from the SI populations around 2810 generation ago, a period perhaps too short for the evolution of a selfing syndrome. We conjecture that the SC population of L. cavanillesii is at an intermediate stage of transition between outcrossing and selfing.     

Inbreeding depression is high in a self‐incompatible perennial herb population but absent in a self‐compatible population showing mixed mating

High inbreeding depression is thought to be one of the major factors preventing evolutionary transitions in hermaphroditic plants from self‐incompatibility (SI) and outcrossing toward self‐compatibility (SC) and selfing. However, when selfing does evolve, inbreeding depression can be quickly purged, allowing the evolution of complete self‐fertilization. In contrast, populations that show intermediate selfing rates (a mixed‐mating system) typically show levels of inbreeding depression similar to those in outcrossing species, suggesting that selection against inbreeding might be responsible for preventing the transition toward complete self‐fertilization. By implication, crosses among populations should reveal patterns of heterosis for mixed‐mating populations that are similar to those expected for outcrossing populations. Using hand‐pollination crosses, we compared levels of inbreeding depression and heterosis between populations of Linaria cavanillesii (Plantaginaceae), a perennial herb showing contrasting mating systems. The SI population showed high inbreeding depression, whereas the SC population displaying mixed mating showed no inbreeding depression. In contrast, we found that heterosis based on between‐population crosses was similar for SI and SC populations. Our results are consistent with the rapid purging of inbreeding depression in the derived SC population, despite the persistence of mixed mating. However, the maintenance of outcrossing after a transition to SC is inconsistent with the prediction that populations that have purged their inbreeding depression should evolve toward complete selfing, suggesting that the transition to SC in L. cavanillesii has been recent. SC in L. cavanillesii thus exemplifies a situation in which the mating system is likely not at an equilibrium with inbreeding depression.     

Breaking linkage between mating compatibility factors: Tetrapolarity in Microbotryum

Linkage of genes determining separate self‐incompatibility mechanisms is a general expectation of sexual eukaryotes that helps to resolve conflicts between reproductive assurance and recombination. However, in some organisms, multiple loci are required to be heterozygous in offspring while segregating independently in meiosis. This condition, termed “tetrapolarity” in basidiomycete fungi, originated in the ancestor to that phylum, and there have been multiple reports of subsequent transitions to “bipolarity” (i.e., linkage of separate mating factors). In the genus Microbotryum, we present the first report of the breaking of linkage between two haploid self‐incompatibility factors and derivation of a tetrapolar breeding system. This breaking of linkage is associated with major alteration of genome structure, with the compatibility factors residing on separate mating‐type chromosome pairs, reduced in size but retaining the structural dimorphism characteristic for regions of recombination suppression. The challenge to reproductive assurance from unlinked compatibility factors may be overcome by the automictic mating system in Microbotryum (i.e., mating among products of the same meiosis). As a curious outcome, this linkage transition and its effects upon outcrossing compatibility rates may reinforce automixis as a mating system. These observations contribute to understanding mating systems and linkage as fundamental principles of sexual life cycles, with potential impacts on conventional wisdom regarding mating‐type evolution.     

Disentangling the effects of breakdown of self-incompatibility and transition to selfing in North American Arabidopsis lyrata

A breakdown of self‐incompatibility (SI) followed by a shift to selfing is commonly observed in the evolution of flowering plants. Both are expected to reduce the levels of heterozygosity and genetic diversity. However, breakdown of SI should most strongly affect the region of the SI locus (S‐locus) because of the relaxation of balancing selection that operates on a functional S‐locus, and a potential selective sweep. In contrast, a transition to selfing should affect the whole genome. We set out to disentangle the effects of breakdown of SI and transition to selfing on the level and distribution of genetic diversity in North American populations of Arabidopsis lyrata. Specifically, we compared sequence diversity of loci linked and unlinked to the S‐locus for populations ranging from complete selfing to fully outcrossing. Regardless of linkage to the S‐locus, heterozygosity and genetic diversity increased with population outcrossing rate. High heterozygosity of self‐compatible individuals in outcrossing populations suggests that SI is not the only factor preventing the evolution of self‐fertilization in those populations. There was a strong loss of diversity in selfing populations, which was more pronounced at the S‐locus. In addition, selfing populations showed an accumulation of derived mutations at the S‐locus. Our results provide evidence that beyond the genome‐wide consequences of the population bottleneck associated with the shift to selfing, the S‐locus of A. lyrata shows a specific signal either reflecting the relaxation of balancing selection or positive selection.     

The evolution of self-compatible and self-incompatible populations in a hermaphroditic perennial, <Emphasis Type="Italic">Trillium camschatcense</Emphasis> (Melanthiaceae)

The evolution of selfing from outcrossing ancestors is known to have occurred repeatedly in angiosperms. Theoretical studies have argued that the transition from outcrossing to selfing is accomplished more easily than the reverse case, and phylogenetic analyses involving self-compatible (SC) and self-incompatible (SI) species has basically supported this assumption. The evolutionary direction of self-compatibility and self-incompatibility was examined in Trillium camschatcense, which contains geographically widespread SC populations, and restricted SI populations. Ecological surveys have revealed that the SC populations were suitable for outcrossing, and selfing in these populations did not confer any fitness advantage. Since reproductive fitness indicates the possibility of an evolutionary shift from self-compatibility to self-incompatibility, the phylogenetic relationships of SI and SC populations of T. camschatcense were investigated based on cpDNA variations and nuclear DNA microsatellite polymorphisms. Although phylogenetic analyses did not provide credible evidence to determine evolutionary direction, the SI populations turned out to be monophyletic with extremely low genetic differentiation. Based on these results, we proposed two possible scenarios for the evolutionary backgrounds of SI and SC populations in T. camschatcense. The plausibility of each scenario was evaluated based on the reproductive and geographical features of the mating systems. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Models of Evolution of Reproductive Isolation

Mathematical models are presented for the evolution of postmating and premating reproductive isolation. In the case of postmating isolation it is assumed that hybrid sterility or inviability is caused by incompatibility of alleles at one or two loci, and evolution of reproductive isolation occurs by random fixation of different incompatibility alleles in different populations. Mutations are assumed to occur following either the stepwise mutation model or the infinite-allele model. Computer simulations by using It?'s stochastic differential equations have shown that in the model used the reproductive isolation mechanism evolves faster in small populations than in large populations when the mutation rate remains the same. In populations of a given size it evolves faster when the number of loci involved is large than when this is small. In general, however, evolution of isolation mechanisms is a very slow process, and it would take thousands to millions of generations if the mutation rate is of the order of 10(-5) per generation. Since gene substitution occurs as a stochastic process, the time required for the establishment of reproductive isolation has a large variance. Although the average time of evolution of isolation mechanisms is very long, substitution of incompatibility genes in a population occurs rather quickly once it starts. The intrapopulational fertility or viability is always very high. In the model of premating isolation it is assumed that mating preference or compatibility is determined by male- and female-limited characters, each of which is controlled by a single locus with multiple alleles, and mating occurs only when the male and female characters are compatible with each other. Computer simulations have shown that the dynamics of evolution of premating isolation mechanism is very similar to that of postmating isolation mechanism, and the mean and variance of the time required for establishment of premating isolation are very large. Theoretical predictions obtained from the present study about the speed of evolution of reproductive isolation are consistent with empirical data available from vertebrate organisms.     

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.     

Breakdown of gametophytic self-incompatibility in subdivided populations

In many hermaphroditic flowering plants, self-fertilization is prevented by self-incompatibility (SI), often controlled by a single locus, the S-locus. In single isolated populations, the maintenance of SI depends chiefly on inbreeding depression and the number of SI alleles at the S-locus. In subdivided populations, however, population subdivision has complicated effects on both the number of SI alleles and the level of inbreeding depression, rendering the maintenance of SI difficult to predict. Here, we explore the conditions for the invasion of a self-compatible mutant in a structured population. We find that the maintenance of SI is strongly compromised when a population becomes subdivided. We show that this effect is mainly caused by the decrease in the local diversity of SI alleles rather than by a change in the dynamics of inbreeding depression. Strikingly, we also find that the diversity of SI alleles at the whole population level is a poor predictor of the maintenance of SI. We discuss the implications of our results for the interpretation of empirical data on the loss of SI in natural populations.     

GENETIC ARCHITECTURE OF INBREEDING DEPRESSION AND THE MAINTENANCE OF GAMETOPHYTIC SELF‐INCOMPATIBILITY

Gametophytic self‐incompatibility (GSI) is a widespread genetic system, which enables hermaphroditic plants to avoid self‐fertilization and mating with close relatives. Inbreeding depression is thought to be the major force maintaining SI; however, inbreeding depression is a dynamical variable that depends in particular on the mating system. In this article we use multilocus, individual‐based simulations to examine the coevolution of SI and inbreeding depression within finite populations. We focus on the conditions for the maintenance of SI when self‐compatible (SC) mutants are introduced in the population by recurrent mutation, and compare simulation results with predictions from an analytical model treating inbreeding depression as a fixed parameter (thereby neglecting effects of purging within the SC subpopulation). In agreement with previous models, we observe that the maintenance of SI is associated with high inbreeding depression and is facilitated by high rates of self‐pollination. Purging of deleterious mutations by SC mutants has little effect on the spread of those mutants as long as most deleterious alleles have weak fitness effects: in this case, the genetic architecture of inbreeding depression has little effect on the maintenance of SI. By contrast, purging may greatly enhance the spread of SC mutants when deleterious alleles have strong fitness effects.     

Evolutionary dynamics of mating system shifts in Arabidopsis lyrata

Outcrossing is the prevalent mode of reproduction in plants and animals despite its substantial costs, while selfing and mixed mating occur at much lower frequency. Comparative research on plants has demonstrated the lability of self‐incompatibility, but there is little information about the transition on a within‐species level from self‐incompatibility to predominant selfing. We studied variation in mating system among 18 populations of Arabidopsis lyrata within a phylogenetic context to shed light on the evolution of selfing. Realized and potential mating systems were assessed by genetic analysis with microsatellite markers and hand‐self‐pollinations on 30 plants from each population. The fraction of self‐incompatible plants in a population was highly correlated with the outcrossing rate, showing that the spread of self‐compatibility is accompanied by or soon followed by an increase in the rate of selfing. The four predominantly selfing populations (outcrossing rates < 0.25) fell into more than one phylogenetic cluster, suggesting that the transition to selfing occurred more than once independently. Hence, A. lyrata offers an opportunity for the comparative analysis of outcrossing as a predominant mode of reproduction in plants and of the causes of the shift to selfing.     

Population genetics of sporophytic self-incompatibility in Senecio squalidus L. (Asteraceae) II: a spatial autocorrelation approach to determining mating behaviour in the presence of low S allele diversity

We recently estimated that as few as six S alleles represent the extent of S locus diversity in a British population of the self-incompatible (SI) coloniser Senecio squalidus (Oxford Ragwort). Despite the predicted constraints to mating imposed by such a low number of S alleles, S. squalidus maintains a strong sporophytic self-incompatibility (SSI) system and there is no evidence for a breakdown of SSI or any obvious negative reproductive consequences for this highly successful coloniser. The present paper assesses mating behaviour in an Oxford S. squalidus population through observations of its effect on spatial patterns of genetic diversity and thus the extent to which it is responsible for ameliorating the potentially detrimental reproductive consequences of low S allele diversity in British S. squalidus. A spatial autocorrelation (SA) treatment of S locus and allozyme polymorphism data for four loci indicates that mating events regularly occur at all the distance classes examined from 60 to 480 m throughout the entire sample population. Less SA is observed for S locus data than for allozyme data in accordance with the hypothesis that SSI and low diversity at the S locus are driving these large-scale mating events. The limited population structure at small distances of 60 m and less observed for SA analysis of the Me-2 locus and by F-statistics for all the allozyme data, is evidence of some local relatedness due to limited seed and pollen dispersal in S. squalidus. However, the overall impression of mating dynamics in this S. squalidus population is that of ample potential mating opportunities with many individuals at large population scales, indicating that reproductive success is not seriously affected by few S alleles available for mating interactions.     

RECONSTRUCTING ORIGINS OF LOSS OF SELF‐INCOMPATIBILITY AND SELFING IN NORTH AMERICAN ARABIDOPSIS LYRATA: A POPULATION GENETIC CONTEXT

Theoretical and empirical comparisons of molecular diversity in selfing and outcrossing plants have primarily focused on long‐term consequences of differences in mating system (between species). However, improving our understanding of the causes of mating system evolution requires ecological and genetic studies of the early stages of mating system transition. Here, we examine nuclear and chloroplast DNA sequences and microsatellite variation in a large sample of populations of Arabidopsis lyrata from the Great Lakes region of Eastern North American that show intra‐ and interpopulation variation in the degree of self‐incompatibility and realized outcrossing rates. Populations show strong geographic clustering irrespective of mating system, suggesting that selfing either evolved multiple times or has spread to multiple genetic backgrounds. Diversity is reduced in selfing populations, but not to the extent of the severe loss of variation expected if selfing evolved due to selection for reproductive assurance in connection with strong founder events. The spread of self‐compatibility in this region may have been favored as colonization bottlenecks following glaciation or migration from Europe reduced standing levels of inbreeding depression. However, our results do not suggest a single transition to selfing in this system, as has been suggested for some other species in the Brassicaceae.     

A one‐locus model of androdioecy with two homomorphic self‐incompatibility groups: expected vs. observed male frequencies

Androdioecy, the occurrence of males and hermaphrodites in a single population, is a rare breeding system because the conditions for maintenance of males are restrictive. In the androdioecious shrub Phillyrea angustifolia, high male frequencies are observed in some populations. The species has a sporophytic self‐incompatibility (SI) system with two self‐incompatibility groups, which ensures that two groups of hermaphrodites can each mate only with the other group, whereas males can fertilize hermaphrodites of both groups. Here, we analyse a population genetic model to investigate the dynamics of such an androdioecious species, assuming that self‐incompatibility and sex phenotypes are determined by a single locus. Our model confirms a previous prediction that a slight reproductive advantage of males relative to hermaphrodites allows the maintenance of males at high equilibrium frequencies. The model predicts different equilibria between hermaphrodites of the two SI groups and males, depending on the male advantage, the initial composition of the population and the population size, whose effect is studied through stochastic simulations. Although the model can generate high male frequencies, observed frequencies are considerably higher than the model predicts. We finally discuss how this model may help explain the large male frequency variation observed in other androdioecious species of Oleaceae: some species show only androdioecious populations, as P. angustifolia, whereas others show populations either completely hermaphrodite or androdioecious.     

The relationship of nucleotide polymorphism, recombination rate and selection in wild tomato species

We analyzed the effects of mating system and recombination rate on single nucleotide polymorphisms using 14 single-copy nuclear loci from single populations of five species of wild tomatoes (Solanum section Lycopersicon). The taxa investigated comprise two self-compatible (SC) and three self-incompatible (SI) species. The observed reduction in nucleotide diversity in the SC populations compared to the SI populations is much stronger than expected under the neutral effects of the mating system on effective population size. Importantly, outgroup sequences available for 11 of the 14 loci yield strong positive correlations between silent nucleotide diversity and silent divergence, indicative of marked among-locus differences in mutation rates and/or selective constraints. Furthermore, using a physical estimate of local recombination rates, we find that silent nucleotide diversity (but not divergence) is positively correlated with recombination rate in two of the SI species. However, this correlation is not nearly as strong as in other well-characterized species (in particular, Drosophila). We propose that nucleotide diversity in Lycopersicon is dominated mainly by differences in neutral mutation rates and/or selective constraints among loci, demographic processes (such as population subdivision), and background selection. In addition, we hypothesize that the soil seed bank plays an important role in the maintenance of the large genetic diversity in the SI species (in particular L. peruvianum).     

On the evolutionary costs of self-incompatibility: incomplete reproductive compensation due to pollen limitation

Pollen limitation affects plants with diverse reproductive systems and ecologies. In self-incompatible (SI) species, pollen limitation may preclude full reproductive compensation for prezygotic rejection of pollen. We present a model designed to explore the effects of incomplete reproductive compensation on evolutionary changes at a modifier locus that regulates the level of SI expression. Our results indicate that incomplete reproductive compensation greatly increases the evolutionary costs of SI, particularly in populations with low S-allele diversity. The evolutionary fate of modifiers of SI expression depends on the rate at which they are transmitted to future generations as well as the effects of SI on offspring number and quality. Partial SI expression can represent a stable condition rather than an evolutionarily transient state between full expression and full suppression. This unanticipated result provides the first theoretical support for the evolutionary stability of such mixed mating systems, the existence of which has recently been documented.     

Mixed mating in natural populations of the chestnut blight fungus, Cryphonectria parasitica

As in plants, fungi exhibit wide variation in reproductive strategies and mating systems. Although most sexually reproducing fungi are either predominantly outcrossing or predominantly selfing, there are some notable exceptions. The haploid, ascomycete chestnut blight pathogen, Cryphonectria parasitica, has previously been shown to have a mixed mating system in one population in USA. In this report, we show that both selfing and outcrossing occur in 10 additional populations of C. parasitica sampled from Japan, Italy, Switzerland and USA. Progeny arrays from each population were assayed for segregation at vegetative incompatibility (vic) and DNA fingerprinting loci. Outcrossing rates (t(m)) were estimated as the proportion of progeny arrays showing segregation at one or more loci, corrected by the probability of nondetection of outcrossing (alpha). Estimates of t(m) varied from 0.74 to 0.97, with the lowest rates consistently detected in USA populations (0.74-0.78). Five populations (four in USA and one in Italy) had t(m) significantly less than 1, supporting the conclusion that these populations exhibit mixed mating. The underlying causes of variation in outcrossing rates among populations of C. parasitica are not known, but we speculate that--as in plants--outcrossing is a function of ecological, demographic and genetic factors.     

Costs of selfing prevent the spread of a self‐compatibility mutation that causes reproductive assurance

In flowering plants, shifts from outcrossing to partial or complete self‐fertilization have occurred independently thousands of times, yet the underlying adaptive processes are difficult to discern. Selfing's ability to provide reproductive assurance when pollination is uncertain is an oft‐cited ecological explanation for its evolution, but this benefit may be outweighed by costs diminishing its selective advantage over outcrossing. We directly studied the fitness effects of a self‐compatibility mutation that was backcrossed into a self‐incompatible (SI) population of Leavenworthia alabamica, illuminating the direction and magnitude of selection on the mating‐system modifier. In array experiments conducted in two years, self‐compatible (SC) plants produced 17–26% more seed, but this advantage was counteracted by extensive seed discounting—the replacement of high‐quality outcrossed seeds by selfed seeds. Using a simple model and simulations, we demonstrate that SC mutations with these attributes rarely spread to high frequency in natural populations, unless inbreeding depression falls below a threshold value (0.57 ≤ δ_threshold ≤ 0.70) in SI populations. A combination of heavy seed discounting and inbreeding depression likely explains why outcrossing adaptations such as self‐incompatibility are maintained generally, despite persistent input of selfing mutations, and frequent limits on outcross seed production in nature.     

Individual and Temporal Variation in Outcrossing Rates and Pollen Flow Patterns in Ceiba pentandra (Malvaceae: Bombacoidea)

Ceiba pentandra is a tropical tree with high rates of selfing in some populations. In mixed‐mating species, variation in selfing is due to changes in adult density or variability of incompatibility systems. The effect of spatial isolation and phenology on selfing rates and pollen flow distances was analyzed using microsatellites in a fragmented population of Ceiba pentandra, in southern Costa Rica. Adult trees within a heterogeneous landscape were classified as grouped or isolated. We compared selfing rates at the individual level, between isolation conditions and 2 yr (2007, 2009), which differed in the number of flowering individuals. Mixed mating was estimated in both years (t_m = 0.624–0.759). Trees mated predominantly by outcrossing, while only a few trees reproduced through selfing. Spatial isolation did not significantly affect outcrossing rates. The progeny of grouped trees was mostly sired by near‐neighbors (<1 km) and by long‐distance pollen flow events in isolated trees. A reduction in the number of flowering individuals in 2009 reduced near‐neighbor matings, increased selfing in grouped trees, and decreased the number of unsampled sires in the progeny. Comparing selfing rates on individuals that flowered in both reproductive periods suggests a flexible mating system. Variation in self‐fertilization rates in this population appears to depend on variation of individual traits, such as genetic variability in self‐incompatibility genes, but it is independent of landscape heterogeneity. In contrast, pollen flow distances depend on local tree density as bats concentrate their foraging between near individuals to maximize energy efficiency.     

Selective acquisition of novel mating type and vegetative incompatibility genes via interspecies gene transfer in the globally invading eukaryote Ophiostoma novo-ulmi

The Dutch elm disease fungus Ophiostoma novo-ulmi, which has destroyed billions of elm trees worldwide, originally invaded Europe as a series of clonal populations with a single mating type (MAT-2) and a single vegetative incompatibility (vic) type. The populations then rapidly became diverse with the appearance of the MAT-1 type and many vegetative incompatibility types. Here, we have investigated the mechanism using isolates from sites in Portugal at which the rapid evolution of O. novo-ulmi populations from clonality to heterogeneity was well established. We show by genetic mapping of vic and MAT loci with AFLP markers and by sequence analysis of MAT loci that this diversification was due to selective acquisition by O. novo-ulmi of the MAT-1 and vic loci from another species, Ophiostoma ulmi. A global survey showed that interspecies transfer of the MAT-1 locus occurred on many occasions as O. novo-ulmi spread across the world. We discuss the possibility that fixation of the MAT-1 and vic loci occurred in response to spread of deleterious viruses in the originally clonal populations. The process demonstrates the potential of interspecies gene transfer for facilitating rapid adaptation of invasive organisms to a new environment.