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

Cell size and selection for stress-induced cell fusion in unicellular eukaryotes

In unicellular organisms, sexual reproduction typically begins with the fusion of two cells (plasmogamy) followed by the fusion of their two haploid nuclei (karyogamy) and finally meiosis. Most work on the evolution of sexual reproduction focuses on the benefits of the genetic recombination that takes place during meiosis. However, the selection pressures that may have driven the early evolution of binary cell fusion, which sets the stage for the evolution of karyogamy by bringing nuclei together in the same cell, have seen less attention. In this paper we develop a model for the coevolution of cell size and binary cell fusion rate. The model assumes that larger cells experience a survival advantage from their larger cytoplasmic volume. We find that under favourable environmental conditions, populations can evolve to produce larger cells that undergo obligate binary cell fission. However, under challenging environmental conditions, populations can evolve to subsequently produce smaller cells under binary cell fission that nevertheless retain a survival advantage by fusing with other cells. The model thus parsimoniously recaptures the empirical observation that sexual reproduction is typically triggered by adverse environmental conditions in many unicellular eukaryotes and draws conceptual links to the literature on the evolution of multicellularity.

Genomic determination of breeding systems and trans-specific evolution of HD MAT genes in suilloid fungi

Studying the signatures of evolution can help to understand genetic processes. Here, we demonstrate how the existence of balancing selection can be used to identify the breeding systems of fungi from genomic data. The breeding systems of fungi are controlled by self-incompatibility loci that determine mating types between potential mating partners, resulting in strong balancing selection at the loci. Within the fungal phylum Basidiomycota, two such self-incompatibility loci, namely HD MAT locus and P/R MAT locus, control mating types of gametes. Loss of function at one or both MAT loci results in different breeding systems and relaxes the MAT locus from balancing selection. By investigating the signatures of balancing selection at MAT loci, one can infer a species' breeding system without culture-based studies. Nevertheless, the extreme sequence divergence among MAT alleles imposes challenges for retrieving full variants from both alleles when using the conventional read-mapping method. Therefore, we employed a combination of read-mapping and local de novo assembly to construct haplotypes of HD MAT alleles from genomes in suilloid fungi (genera Suillus and Rhizopogon). Genealogy and pairwise divergence of HD MAT alleles showed that the origins of mating types predate the split between these two closely related genera. High sequence divergence, trans-specific polymorphism, and the deeply diverging genealogy confirm the long-term functionality and multiallelic status of HD MAT locus in suilloid fungi. This work highlights a genomics approach to studying breeding systems regardless of the culturability of organisms based on the interplay between evolution and genetics.

Origin and persistence of polymorphism in loci targeted by disassortative preference: a general model

The emergence and persistence of polymorphism within populations generally requires specific regimes of natural or sexual selection. Here, we develop a unified theoretical framework to explore how polymorphism at targeted loci can be generated and maintained by either disassortative mating choice or balancing selection due to, for example, heterozygote advantage. To this aim, we model the dynamics of alleles at a single locus A in a population of haploid individuals, where reproductive success depends on the combination of alleles carried by the parents at locus A. Our theoretical study of the model confirms that the conditions for the persistence of a given level of allelic polymorphism depend on the relative reproductive advantages among pairs of individuals. Interestingly, equilibria with unbalanced allelic frequencies were shown to emerge from successive introduction of mutants. We then investigate the role of the function linking allelic divergence to reproductive advantage on the evolutionary fate of alleles within the population. Our results highlight the significance of the shape of this function for both the number of alleles maintained and their level of genetic divergence. Large number of alleles are maintained with substantial replacement of alleles, when disassortative advantage slowly increases with allelic differentiation . In contrast, few highly differentiated alleles are predicted to be maintained when genetic differentiation has a strong effect on disassortative advantage. These opposite effects predicted by our model explain how disassortative mate choice may lead to various levels of allelic differentiation and polymorphism, and shed light on the effect of mate preferences on the persistence of balanced and unbalanced polymorphism in natural population.

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.

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.

Drivers of Mating Type Composition in Tetrahymena thermophila

Sex offers advantages even in primarily asexual species. Some ciliates appear to utilize such reproductive strategy with many mating types. However, the factors determining the composition of mating types in the unicellular ciliate Tetrahymena thermophila are poorly understood, and this is further complicated by non-Mendelian determination of mating type in the offspring. We therefore developed a novel population genetics model to predict how various factors influence the dynamics of mating type composition, including natural selection. The model predicted either the coexistence of all seven mating types or fixation of a single mating type in a population, depending on parameter combinations, irrespective of natural selection. To understand what factor(s) may be more influential and to test the validity of theoretical prediction, five replicate populations were maintained in laboratory such that several factors could be controlled or measured. Whole-genome sequencing was used to identify newly arising mutations and determine mating type composition. Strikingly, all populations were found to be driven by strong selection on newly arising beneficial mutations to fixation of their carrying mating types, and the trajectories of speed to fixation agreed well with our theoretical predictions. This study illustrates the evolutionary strategies that T. thermophila can utilize to optimize population fitness.