Origin and diversification dynamics of self-incompatibility haplotypes

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.

The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

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.

Widespread coexistence of self-compatible and self-incompatible phenotypes in a diallelic self-incompatibility system in Ligustrum vulgare (Oleaceae)

The breakdown of self-incompatibility (SI) in angiosperms is one of the most commonly observed evolutionary transitions. While multiple examples of SI breakdown have been documented in natural populations, there is strikingly little evidence of stable within-population polymorphism with both inbreeding (self-compatible) and outcrossing (self-incompatible) individuals. This absence of breeding system polymorphism corroborates theoretical expectations that predict that in/outbreeding polymorphism is possible only under very restricted conditions. However, theory also predicts that a diallelic sporophytic SI system should facilitate the maintenance of such polymorphism. We tested this prediction by studying the breeding system of Ligustrum vulgare L., an insect-pollinated hermaphroditic species of the Oleaceae family. Using stigma tests with controlled pollination and paternity assignment of open-pollinated progenies, we confirmed the existence of two self-incompatibility groups in this species. We also demonstrated the occurrence of self-compatible individuals in different populations of Western Europe arising from a mutation affecting the functioning of the pollen component of SI. Our results show that the observed low frequency of self-compatible individuals in natural populations is compatible with theoretical predictions only if inbreeding depression is very high.

Non-self recognition-based self-incompatibility can alternatively promote or prevent introgression

Self-incompatibility alleles (S-alleles), which prevent self-fertilisation in plants, have historically been expected to benefit from negative frequency-dependent selection and invade when introduced to a new population through gene flow. However, the most taxonomically widespread form of self-incompatibility, the ribonuclease-based system ancestral to the core eudicots, functions through collaborative non-self recognition, which can affect both short-term patterns of gene flow and the long-term process of S-allele diversification. We analysed a model of S-allele evolution in two populations connected by migration, focussing on comparisons among the fates of S-alleles initially unique to each population and those shared among populations. We found that both shared and unique S-alleles from the population with more unique S-alleles were usually fitter compared with S-alleles from the population with fewer S-alleles. Resident S-alleles often became extinct and were replaced by migrant S-alleles, although this outcome could be averted by pollen limitation or biased migration. Collaborative non-self recognition will usually either result in the whole-sale replacement of S-alleles from one population with those from another or else disfavour introgression of S-alleles altogether.

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.

The loss of self-incompatibility in a range expansion

It is commonly observed that plant species' range margins are enriched for increased selfing rates and, in otherwise self-incompatible species, for self-compatibility (SC). This has often been attributed to a response to selection under mate and/or pollinator limitation. However, range expansion can also cause reduced inbreeding depression, and this could facilitate the evolution of selfing in the absence of mate or pollinator limitation. Here, we explore this idea using spatially explicit individual-based simulations of a range expansion, in which inbreeding depression, variation in self-incompatibility (SI), and mate availability evolve. Under a wide range of conditions, the simulated range expansion brought about the evolution of selfing after the loss of SI in range-marginal populations. Under conditions of high recombination between the self-incompatibility locus (S-locus) and viability loci, SC remained marginal in the expanded metapopulation and could not invade the range core, which remained self-incompatible. In contrast, under low recombination and migration rates, SC was frequently able to displace SI in the range core by maintaining its association with a genomic background with purged genetic load. We conclude that the evolution of inbreeding depression during a range expansion promotes the evolution of SC at range margins, especially under high rates of recombination.‬.

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.

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.

Asymmetrical diversification of the receptor-ligand interaction controlling self-incompatibility in Arabidopsis

How two-component genetic systems accumulate evolutionary novelty and diversify in the course of evolution is a fundamental problem in evolutionary systems biology. In the Brassicaceae, self-incompatibility (SI) is a spectacular example of a diversified allelic series in which numerous highly diverged receptor-ligand combinations are segregating in natural populations. However, the evolutionary mechanisms by which new SI specificities arise have remained elusive. Using in planta ancestral protein reconstruction, we demonstrate that two allelic variants segregating as distinct receptor-ligand combinations diverged through an asymmetrical process whereby one variant has retained the same recognition specificity as their (now extinct) putative ancestor, while the other has functionally diverged and now represents a novel specificity no longer recognized by the ancestor. Examination of the structural determinants of the shift in binding specificity suggests that qualitative rather than quantitative changes of the interaction are an important source of evolutionary novelty in this highly diversified receptor-ligand system.

Mating system variation in hybrid zones: facilitation, barriers and asymmetries to gene flow

Plant mating systems play a key role in structuring genetic variation both within and between species. In hybrid zones, the outcomes and dynamics of hybridization are usually interpreted as the balance between gene flow and selection against hybrids. Yet, mating systems can introduce selective forces that alter these expectations; with diverse outcomes for the level and direction of gene flow depending on variation in outcrossing and whether the mating systems of the species pair are the same or divergent. We present a survey of hybridization in 133 species pairs from 41 plant families and examine how patterns of hybridization vary with mating system. We examine if hybrid zone mode, level of gene flow, asymmetries in gene flow and the frequency of reproductive isolating barriers vary in relation to mating system/s of the species pair. We combine these results with a simulation model and examples from the literature to address two general themes: (1) the two-way interaction between introgression and the evolution of reproductive systems, and (2) how mating system can facilitate or restrict interspecific gene flow. We conclude that examining mating system with hybridization provides unique opportunities to understand divergence and the processes underlying reproductive isolation.

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.

Evolutionary Pathways for the Generation of New Self-Incompatibility Haplotypes in a Nonself-Recognition System

Self-incompatibility (SI) is a genetically based recognition system that functions to prevent self-fertilization and mating among related plants. An enduring puzzle in SI is how the high diversity observed in nature arises and is maintained. Based on the underlying recognition mechanism, SI can be classified into two main groups: self-recognition (SR) and nonself-recognition (NSR). Most work has focused on diversification within SR systems despite expected differences between the two groups in the evolutionary pathways and outcomes of diversification. Here, we use a deterministic population genetic model and stochastic simulations to investigate how novel S-haplotypes evolve in a gametophytic NSR [SRNase/S Locus F-box (SLF)] SI system. For this model, the pathways for diversification involve either the maintenance or breakdown of SI and can vary in the order of mutations of the female (SRNase) and male (SLF) components. We show analytically that diversification can occur with high inbreeding depression and self-pollination, but this varies with evolutionary pathway and level of completeness (which determines the number of potential mating partners in the population), and, in general, is more likely for lower haplotype number. The conditions for diversification are broader in stochastic simulations of finite population size. However, the number of haplotypes observed under high inbreeding and moderate-to-high self-pollination is less than that commonly observed in nature. Diversification was observed through pathways that maintain SI as well as through self-compatible intermediates. Yet the lifespan of diversified haplotypes was sensitive to their level of completeness. By examining diversification in a NSR SI system, this model extends our understanding of the evolution and maintenance of haplotype diversity observed in a recognition system common in flowering plants.

Patterns of Polymorphism at the Self-Incompatibility Locus in 1,083 Arabidopsis thaliana Genomes

Although the transition to selfing in the model plant Arabidopsis thaliana involved the loss of the self-incompatibility (SI) system, it clearly did not occur due to the fixation of a single inactivating mutation at the locus determining the specificities of SI (the S-locus). At least three groups of divergent haplotypes (haplogroups), corresponding to ancient functional S-alleles, have been maintained at this locus, and extensive functional studies have shown that all three carry distinct inactivating mutations. However, the historical process of loss of SI is not well understood, in particular its relation with the last glaciation. Here, we took advantage of recently published genomic resequencing data in 1,083 Arabidopsis thaliana accessions that we combined with BAC sequencing to obtain polymorphism information for the whole S-locus region at a species-wide scale. The accessions differed by several major rearrangements including large deletions and interhaplogroup recombinations, forming a set of haplogroups that are widely distributed throughout the native range and largely overlap geographically. “Relict” A. thaliana accessions that directly derive from glacial refugia are polymorphic at the S-locus, suggesting that the three haplogroups were already present when glacial refugia from the last Ice Age became isolated. Interhaplogroup recombinant haplotypes were highly frequent, and detailed analysis of recombination breakpoints suggested multiple independent origins. These findings suggest that the complete loss of SI in A. thaliana involved independent self-compatible mutants that arose prior to the last Ice Age, and experienced further rearrangements during postglacial colonization.

The unusual S locus of Leavenworthia is composed of two sets of paralogous loci

The Leavenworthia self-incompatibility locus (S locus) consists of paralogs (Lal2, SCRL) of the canonical Brassicaceae S locus genes (SRK, SCR), and is situated in a genomic position that differs from the ancestral one in the Brassicaceae. Unexpectedly, in a small number of Leavenworthia alabamica plants examined, sequences closely resembling exon 1 of SRK have been found, but the function of these has remained unclear. BAC cloning and expression analyses were employed to characterize these SRK-like sequences. An SRK-positive Bacterial Artificial Chromosome clone was found to contain complete SRK and SCR sequences located close by one another in the derived genomic position of the Leavenworthia S locus, and in place of the more typical Lal2 and SCRL sequences. These sequences are expressed in stigmas and anthers, respectively, and crossing data show that the SRK/SCR haplotype is functional in self-incompatibility. Population surveys indicate that < 5% of Leavenworthia S loci possess such alleles. An ancestral translocation or recombination event involving SRK/SCR and Lal2/SCRL likely occurred, together with neofunctionalization of Lal2/SCRL, and both haplotype groups now function as Leavenworthia S locus alleles. These findings suggest that S locus alleles can have distinctly different evolutionary origins.

Evolutionary history of two pollen self-incompatibility factors reveals alternate routes to self-compatibility within Solanum

PREMISE OF THE STUDY

Self-incompatibility (SI) prevents self-fertilization and reduces inbreeding. While SI is common in plants, transitions to self-compatibility (SC) occur frequently. Little is known about the genetic changes and evolutionary steps underlying these shifts.

METHODS

In the Solanaceae, SI is gametophytic, with specificity determined by S-RNases in the pistil and S-locus F-box proteins (SLFs) in pollen. We examined the role of two pollen factors, Cullin1 (CUL1) and SLF-23, in SI → SC transitions in wild tomato species from the Arcanum species group (Solanum arcanum, S. neorickii, and S. chmielewskii). Pollen compatibility was assessed on tester lines that reject pollen lacking functional SLF-23 or CUL1. Complementation tests, gene sequencing, and phylogenetic analyses were used to characterize both functional and nonfunctional alleles.

KEY RESULTS

We found evidence for multiple independent SI → SC transitions. In S. arcanum and S. chmielewskii, SC is caused by loss of pistil S-RNase activity, while in S. neorickii SC is associated with expression of a functional SLF-23 that recognizes the S₉ type S-RNase expressed in its pistils. Interestingly, we found identical deletion mutations in CUL1 exon 7 of S. chmielewskii as previously seen in S. habrochaites.

CONCLUSIONS

Mating system transitions in the Arcanum group have occurred via both pistil loss-of-function and pollen gain-of-function SC mutations. Mutations common to S. chmielewskii and S. habrochaites must have arisen in a common ancestor, possibly to the entire tomato clade, then became fixed in different lineages after loss of pistil-side SI function.

Elucidation of the genetic architecture of self-incompatibility in olive: Evolutionary consequences and perspectives for orchard management

The olive (Olea europaea L.) is a typical important perennial crop species for which the genetic determination and even functionality of self‐incompatibility (SI) are still largely unresolved. It is still not known whether SI is under gametophytic or sporophytic genetic control, yet fruit production in orchards depends critically on successful ovule fertilization. We studied the genetic determination of SI in olive in light of recent discoveries in other genera of the Oleaceae family. Using intra‐ and interspecific stigma tests on 89 genotypes representative of species‐wide olive diversity and the compatibility/incompatibility reactions of progeny plants from controlled crosses, we confirmed that O. europaea shares the same homomorphic diallelic self‐incompatibility (DSI) system as the one recently identified in Phillyrea angustifolia and Fraxinus ornus. SI is sporophytic in olive. The incompatibility response differs between the two SI groups in terms of how far pollen tubes grow before growth is arrested within stigma tissues. As a consequence of this DSI system, the chance of cross‐incompatibility between pairs of varieties in an orchard is high (50%) and fruit production may be limited by the availability of compatible pollen. The discovery of the DSI system in O. europaea will undoubtedly offer opportunities to optimize fruit production.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Overcoming self-incompatibility in grasses: a pathway to hybrid breeding

Allogamous grasses exhibit an effective two-locus gametophytic self-incompatibility (SI) system, limiting the range of breeding techniques applicable for cultivar development. Current breeding methods based on populations are characterized by comparably low genetic gains for important traits such as biomass yield. To implement more efficient breeding schemes, the overall understanding of the SI system is crucial as are the mechanisms involved in the breakdown of SI. Self-fertile variants in outcrossing grasses have been studied, and the current level of knowledge includes approximate gene locations, linked molecular markers and first hypotheses on their mode of action. Environmental conditions increasing seed set upon self-pollination have also been described. Even though some strategies were proposed to take advantage of self-fertility, there have, so far, not been changes in the methods applied in cultivar development for allogamous grasses. In this review, we describe the current knowledge about self-fertility in allogamous grasses and outline strategies to incorporate this trait for implementation in synthetic and hybrid breeding schemes.

Mutations in two pollen self-incompatibility factors in geographically marginal populations of Solanum habrochaites impact mating system transitions and reproductive isolation

PREMISE OF THE STUDY

Self-incompatibility (SI) is a mechanism that prevents inbreeding in many plant species. The mutational breakdown of SI occurs frequently, yet relatively little is known about the evolutionary steps involved in the progressive loss of pistil and pollen SI function.

METHODS

In Solanaceae, SI is the S-RNase-based gametophytic type. We used SI and SC populations of the wild tomato species Solanum habrochaites to study natural variation for two pollen SI factors: a Cullin1 (CUL1) protein and an S-locus F-box protein (SLF-23). Pollen compatibility was assessed on an allotriploid tester line encoding an S-RNase recognized by SLF-23. Both pollen factors are required for compatibility on this tester line. Complementation tests and gene sequencing were used to identify mutations in CUL1 or SLF-23.

KEY RESULTS

We detected loss-of-function mutations in CUL1 and/or SLF-23 in SC populations collected near the northern and southern geographic margins of this taxon's natural range. Nonmarginal SC and all SI accessions expressed mostly functional alleles of these pollen factors. Comparison of the CUL1 sequences identified several shared deletion mutations present in both northern and southern margin SC accessions.

CONCLUSIONS

Loss-of-function mutations in CUL1 and SLF-23 likely became fixed relatively late during SI to SC transitions, after loss of pistil SI function. Mutations in CUL1 establish unilateral incompatibility with SI populations and strengthen reproductive isolation. Point mutations common to northern and southern SC biotypes likely derive from shared ancestral variants found in more central SI populations.

Helitron-like transposons contributed to the mating system transition from out-crossing to self-fertilizing in polyploid Brassica napus L

The mating system transition in polyploid Brassica napus (AACC) from out-crossing to selfing is a typical trait to differentiate it from their diploid progenitors. Elucidating the mechanism of mating system transition has profound consequences for understanding the speciation and evolution in B. napus. Functional complementation experiment has shown that the insertion of 3.6 kb into the promoter of self-incompatibility male determining gene, BnSP11-1 leads to its loss of function in B. napus. The inserted fragment was found to be a non-autonomous Helitron transposon. Further analysis showed that the inserted 3.6 kb non-autonomous Helitron transposon was widely distributed in B. napus accessions which contain the S haplotype BnS-1. Through promoter deletion analysis, an enhancer and a putative cis-regulatory element (TTCTA) that were required for spatio-temporal specific expression of BnSP11-1 were identified, and both might be disrupted by the insertion of Helitron transposon. We suggested that the insertion of Helitron transposons in the promoter of BnSP11-1 gene had altered the mating system and might facilitated the speciation of B. napus. Our findings have profound consequences for understanding the self-compatibility in B. napus as well as for the trait variations during evolutionary process of plant polyploidization.

How Have Self-Incompatibility Haplotypes Diversified? Generation of New Haplotypes during the Evolution of Self-Incompatibility from Self-Compatibility

I developed a gametophytic self-incompatibility (SI) model to study the conditions leading to diversification in SI haplotypes. In the model, the SI system is assumed to be incomplete, and the pollen expressing a given specificity is not fully rejected by the pistils expressing the same specificity. I also assumed that mutations can occur that enhance the rejection of pollen by pistils with the same haplotype variant and reduce rejection by pistils with other variants in the same haplotype. I found that if such mutations occur, the new haplotypes (mutant variants) can stably coexist with the ancestral haplotype in which the mutant arose. This is because pollen bearing the new haplotype is most strongly rejected by pistils bearing the same new haplotype among the pistils in the population; hence, negative frequency-dependent selection prevents their fixation. I also performed simulations and found that the nearly complete SI system evolves from completely self-compatible populations and that SI haplotypes can increase to about 40-50 within a few thousand generations. On the basis of my findings, I propose that diversification of SI haplotypes occurred during the evolution of SI from self-compatibility.

Evidence for the long-term maintenance of a rare self-incompatibility system in Oleaceae

A rare homomorphic diallelic self-incompatibility (DSI) system discovered in Phillyrea angustifolia (family Oleaceae, subtribe Oleinae) can promote the transition from hermaphroditism to androdioecy. If widespread and stable in Oleaceae, DSI may explain the exceptionally high rate of androdioecious species reported in this plant family. Here, we set out to determine whether DSI occurs in another Oleaceae lineage. We tested for DSI in subtribe Fraxininae, a lineage that diverged from subtribe Oleinae c. 40 million yr ago. We explored the compatibility relationships in Fraxinus ornus using 81 hermaphrodites and 25 males from one natural stand and two naturalized populations using intra- and interspecific stigma tests performed on F. ornus and P. angustifolia testers. We uncovered a DSI system with hermaphrodites belonging to one of two self-incompatibility (SI) groups and males compatible with both groups, making for a truly androdioecious reproductive system. The two human-founded populations contained only one of the two SI groups. Our results provide evidence for the evolutionary persistence of DSI. We discuss how its stability over time may have affected transitions to other sexual systems, such as dioecy.

Selection of sporophytic and gametophytic self-incompatibility in the absence of a superlocus

Self-incompatibility (SI) is a complex trait that enforces outcrossing in plant populations. SI generally involves tight linkage of genes coding for the proteins that underlie self-pollen detection and pollen identity specification. Here, we develop two-locus genetic models to address the question of whether sporophytic SI (SSI) and gametophytic SI (GSI) can invade populations of self-compatible plants when there is no linkage or weak linkage of the underlying pollen detection and identity genes (i.e., no S-locus supergene). The models assume that SI evolves as a result of exaptation of genes formerly involved in functions other than SI. Model analysis reveals that SSI and GSI can invade populations even when the underlying genes are loosely linked, provided that inbreeding depression and selfing rate are sufficiently high. Reducing recombination between these genes makes conditions for invasion more lenient. These results can help account for multiple, independent evolution of SI systems as seems to have occurred in the angiosperms.

The joint evolution and maintenance of self-incompatibility with gynodioecy or androdioecy

Mating systems show two kinds of frequent transitions: from hermaphroditism to dioecy, gynodioecy or androdioecy, or from self-incompatibility (SI) to self-compatibility (SC). While models have mostly investigated these two kinds of transitions as independent, empirical observations suggest that, to some extent, they can evolve jointly. Here, we study the joint evolution and maintenance of SI and androdioecy or SI and gynodioecy by the means of phenotypic models. Our models focus on three parameters: the unisexuals׳ advantage relative to that of the hermaphrodites due to resource reallocation, inbreeding depression and the selfing rate. We assume no pollen limitation or discounting. We show that SI helps the maintenance of androdioecy, but favors the loss of gynodioecy, and also that androdioecy facilitates the maintenance of SI, whereas gynodioecy does not affect it. We finally investigate how gynodioecy and androdioecy may affect the diversification of SI groups, especially considering an evolutionary pathway through SC intermediates. We show that while androdioecy prevents the increase of the number of SI groups, under certain conditions of inbreeding depression and selfing rates, gynodioecy allows it.

Gene duplication and genetic exchange drive the evolution of S-RNase-based self-incompatibility in Petunia

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16-20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16-20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations.

Initial invasion of gametophytic self-incompatibility alleles in the absence of tight linkage between pollen and pistil S alleles

In homomorphic self-incompatibility (SI) systems of plants, the loci controlling the pollen and pistil types are tightly linked, and this prevents the generation of compatible combinations of alleles expressing pollen and pistil types, which would result in self-fertilization. We modeled the initial invasion of the first pollen and pistil alleles in gametophytic SI to determine whether these alleles can stably coexist in a population without tight linkage. We assume pollen and pistil loci each carry an incompatibility allele S and an allele without an incompatibility function N. We assume that pollen with an S allele are incompatible with pistils carrying S alleles, whereas other crosses are compatible. Ovules in pistils carrying an S allele suffer viability costs because recognition consumes resources. We found that the cost of carrying a pistil S allele allows pollen and pistil S alleles to coexist in a stable equilibrium if linkage is partial. This occurs because parents that carry pistil S alleles but are homozygous for pollen N alleles cannot avoid self-fertilization; however, they suffer viability costs. Hence, pollen N alleles are selected again. When pollen and pistil S alleles can coexist in a polymorphic equilibrium, selection will favor tighter linkage.

Trait transitions in explicit ecological and genomic contexts: plant mating systems as case studies

Plants are astonishingly diverse in how they reproduce sexually, and the study of plant mating systems provides some of the most compelling cases of parallel and independent evolutionary transitions. In this chapter, we review how the massive amount of genomic data being produced is allowing long-standing predictions from ecological and evolutionary theory to be put to test. After a review of theoretical predictions about the importance of considering the genomic architecture of the mating system, we focus on a set of recent discoveries on how the mating system is controlled in a variety of model and non-model species. In parallel, genomic approaches have revealed the complex interaction between the evolution of genes controlling mating systems and genome evolution, both genome-wide and in the mating system control region. In several cases, major transitions in the mating system can be clearly associated with important ecological changes, hence illuminating an important interplay between ecological and genomic approaches. We also list a number of major unsolved questions that remain for the field, and highlight foreseeable conceptual developments that are likely to play a major role in our understanding of how plant mating systems evolve in Nature.

Effects of pollen availability and the mutation bias on the fixation of mutations disabling the male specificity of self-incompatibility

The evolution of self-compatibility (SC) by the loss of self-incompatibility (SI) is regarded as one of the most frequent transitions in flowering plants. SI systems are generally characterized by specific interactions between the male and female specificity genes encoded at the S-locus. Recent empirical studies have revealed that the evolution of SC is often driven by male SC-conferring mutations at the S-locus rather than by female mutations. In this study, using a forward simulation model, we compared the fixation probabilities of male vs. female SC-conferring mutations at the S-locus. We explicitly considered the effects of pollen availability in the population and bias in the occurrence of SC-conferring mutations on the male and female specificity genes. We found that male SC-conferring mutations were indeed more likely to be fixed than were female SC-conferring mutations in a wide range of parameters. This pattern was particularly strong when pollen availability was relatively high. Under such a condition, even if the occurrence of mutations was biased strongly towards the female specificity gene, male SC-conferring mutations were much more often fixed. Our study demonstrates that fixation probabilities of those two types of mutation vary strongly depending on ecological and genetic conditions, although both types result in the same evolutionary consequence-the loss of SI.

Transitions from reproductive systems governed by two self-incompatible loci to one in fungi

Self-incompatibility (SI), a reproductive system broadly present in plants, chordates, fungi, and protists, might be controlled by one or several multiallelic loci. How a transition in the number of SI loci can occur and the consequences of such events for the population's genetics and dynamics have not been studied theoretically. Here, we provide analytical descriptions of two transition mechanisms: linkage of the two SI loci (scenario 1) and the loss of function of one incompatibility gene within a mating type of a population with two SI loci (scenario 2). We show that invasion of populations by the new mating type form depends on whether the fitness of the new type is lowered, and on the allelic diversity of the SI loci and the recombination between the two SI loci in the starting population. Moreover, under scenario 1, it also depends on the frequency of the SI alleles that became linked. We demonstrate that, following invasion, complete transitions in the reproductive system occurs under scenario 2 and is predicted only for small populations under scenario 1. Interestingly, such events are associated with a drastic reduction in mating type number.