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

The role of promiscuous molecular recognition in the evolution of RNase-based self-incompatibility in plants

How do biological networks evolve and expand? We study these questions in the context of the plant collaborative-non-self recognition self-incompatibility system. Self-incompatibility evolved to avoid self-fertilization among hermaphroditic plants. It relies on specific molecular recognition between highly diverse proteins of two families: female and male determinants, such that the combination of genes an individual possesses determines its mating partners. Though highly polymorphic, previous models struggled to pinpoint the evolutionary trajectories by which new specificities evolved. Here, we construct a novel theoretical framework, that crucially affords interaction promiscuity and multiple distinct partners per protein, as is seen in empirical findings disregarded by previous models. We demonstrate spontaneous self-organization of the population into distinct "classes" with full between-class compatibility and a dynamic long-term balance between class emergence and decay. Our work highlights the importance of molecular recognition promiscuity to network evolvability. Promiscuity was found in additional systems suggesting that our framework could be more broadly applicable.

Revisiting the number of self-incompatibility alleles in finite populations: From old models to new results

Under gametophytic self-incompatibility (GSI), plants are heterozygous at the self-incompatibility locus (S-locus) and can only be fertilized by pollen with a different allele at that locus. The last century has seen a heated debate about the correct way of modelling the allele diversity in a GSI population that was never formally resolved. Starting from an individual-based model, we derive the deterministic dynamics as proposed by Fisher (The genetical theory of natural selection - A complete, Variorum edition, Oxford University Press, 1958) and compute the stationary S-allele frequency distribution. We find that the stationary distribution proposed by Wright (Evolution, 18, 609, 1964) is close to our theoretical prediction, in line with earlier numerical confirmation. Additionally, we approximate the invasion probability of a new S-allele, which scales inversely with the number of resident S-alleles. Lastly, we use the stationary allele frequency distribution to estimate the population size of a plant population from an empirically obtained allele frequency spectrum, which complements the existing estimator of the number of S-alleles. Our expression of the stationary distribution resolves the long-standing debate about the correct approximation of the number of S-alleles and paves the way for new statistical developments for the estimation of the plant population size based on S-allele frequencies.

Parthenogenesis and the Evolution of Anisogamy

Recently, it was pointed out that classic models for the evolution of anisogamy do not take into account the possibility of parthenogenetic reproduction, even though sex is facultative in many relevant taxa (e.g., algae) that harbour both anisogamous and isogamous species. Here, we complement this recent analysis with an approach where we assume that the relationship between progeny size and its survival may differ between parthenogenetically and sexually produced progeny, favouring either the former or the latter. We show that previous findings that parthenogenesis can stabilise isogamy relative to the obligate sex case, extend to our scenarios. We additionally investigate two different ways for one mating type to take over the entire population. First, parthenogenesis can lead to biased sex ratios that are sufficiently extreme that one type can displace the other, leading to de facto asexuality for the remaining type that now lacks partners to fuse with. This process involves positive feedback: microgametes, being numerous, lack opportunities for syngamy, and should they proliferate parthenogenetically, the next generation makes this asexual route even more prominent for microgametes. Second, we consider mutations to strict asexuality in producers of micro- or macrogametes, and show that the prospects of asexual invasion depend strongly on the mating type in which the mutation arises. Perhaps most interestingly, we also find scenarios in which parthenogens have an intrinsic survival advantage yet facultatively sexual isogamous populations are robust to the invasion of asexuals, despite us assuming no genetic benefits of recombination. Here, equal contribution from both mating types to zygotes that are sufficiently well provisioned can outweigh the additional costs associated with syngamy.

Understanding evolutionary and ecological dynamics using a continuum limit

Continuum limits in the form of stochastic differential equations are typically used in theoretical population genetics to account for genetic drift or more generally, inherent randomness of the model. In evolutionary game theory and theoretical ecology, however, this method is used less frequently to study demographic stochasticity. Here, we review the use of continuum limits in ecology and evolution. Starting with an individual-based model, we derive a large population size limit, a (stochastic) differential equation which is called continuum limit. By example of the Wright-Fisher diffusion, we outline how to compute the stationary distribution, the fixation probability of a certain type, and the mean extinction time using the continuum limit. In the context of the logistic growth equation, we approximate the quasi-stationary distribution in a finite population.

Switching environments, synchronous sex, and the evolution of mating types

While facultative sex is common in sexually reproducing species, for reasons of tractability most mathematical models assume that such sex is asynchronous in the population. In this paper, we develop a model of switching environments to instead capture the effect of an entire population transitioning synchronously between sexual and asexual modes of reproduction. We use this model to investigate the evolution of the number of self-incompatible mating types in finite populations, which empirically can range from two to thousands. When environmental switching is fast, we recover the results of earlier studies that implicitly assumed populations were engaged in asynchronous sexual reproduction. However when the environment switches slowly, we see deviations from previous asynchronous theory, including a lower number of mating types at equilibrium and bimodality in the stationary distribution of mating types. We provide analytic approximations for both the fast and slow switching regimes, as well as a numerical scheme based on the Kolmogorov equations for the system to quickly evaluate the model dynamics at intermediate parameters. Our approach exploits properties of integer partitions in number theory. We also demonstrate how additional biological processes such as selective sweeps can be accounted for in this switching environment framework, showing that beneficial mutations can further erode mating type diversity in synchronous facultatively sexual populations.

Fitness differences suppress the number of mating types in evolving isogamous species

Sexual reproduction is not always synonymous with the existence of two morphologically different sexes; isogamous species produce sex cells of equal size, typically falling into multiple distinct self-incompatible classes, termed mating types. A long-standing open question in evolutionary biology is: what governs the number of these mating types across species? Simple theoretical arguments imply an advantage to rare types, suggesting the number of types should grow consistently; however, empirical observations are very different. While some isogamous species exhibit thousands of mating types, such species are exceedingly rare, and most have fewer than 10. In this paper, we present a mathematical analysis to quantify the role of fitness variation-characterized by different mortality rates-in determining the number of mating types emerging in simple evolutionary models. We predict that the number of mating types decreases as the variance of mortality increases.

Invasion and Extinction Dynamics of Mating Types Under Facultative Sexual Reproduction

In sexually reproducing isogamous species, syngamy between gametes is generally not indiscriminate, but rather restricted to occurring between complementary self-incompatible mating types. A longstanding question regards the evolutionary pressures that control the number of mating types observed in natural populations, which ranges from two to many thousands. Here, we describe a population genetic null model of this reproductive system, and derive expressions for the stationary probability distribution of the number of mating types, the establishment probability of a newly arising mating type, and the mean time to extinction of a resident type. Our results yield that the average rate of sexual reproduction in a population correlates positively with the expected number of mating types observed. We further show that the low number of mating types predicted in the rare-sex regime is primarily driven by low invasion probabilities of new mating type alleles, with established resident alleles being very stable over long evolutionary periods. Moreover, our model naturally exhibits varying selection strength dependent on the number of resident mating types. This results in higher extinction and lower invasion rates for an increasing number of residents.