Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances

Population Genomics Reveals Demographic History and Climate Adaptation in Japanese Arabidopsis halleri

Climate oscillations in the Quaternary forced species to major latitudinal or altitudinal range shifts. It has been suggested that adaptation concomitant with range shifts plays key roles in species responses during climate oscillations, but the role of selection for local adaptation to climatic changes remains largely unexplored. Here, we investigated population structure, demographic history and signatures of climate-driven selection based on genome-wide polymorphism data of 141 Japanese Arabidopsis halleri individuals, with European ones as outgroups. Coalescent-based analyses suggested a genetic differentiation between Japanese subpopulations since the Last Glacial Period (LGP), which would have contributed to shaping the current pattern of population structure. Population demographic analysis revealed the population size fluctuations in the LGP, which were particularly prominent since the subpopulations started to diverge (∼50, 000 years ago). The ecological niche modeling predicted the geographic or distribution range shifts from southern coastal regions to northern coastal and mountainous areas, possibly in association with the population size fluctuations. Through genome-wide association analyses of bioclimatic variables and selection scans, we investigated whether climate-associated loci are enriched in the extreme tails of selection scans, and demonstrated the prevailing signatures of selection, particularly toward a warmer climate in southern subpopulations and a drier environment in northern subpopulations, which may have taken place during or after the LGP. Our study highlights the importance of integrating climate associations, selection scans and population demographic analyses for identifying genomic signatures of population-specific adaptation, which would also help us predict the evolutionary responses to future climate changes.

Characterisation of the Gillenia S-locus provides insight into evolution of the nonself-recognition self-incompatibility system in apple

Self-incompatibility (SI) in plants has evolved independently multiple times and S-RNase-based gametophytic self-incompatibility (GSI) is most common. The Rosaceae family possesses both self-recognition (Prunus) and nonself-recognition (Malus) GSI systems, and the latter is widespread in flowering plants. Gillenia trifoliata is a Rosaceae species related to Prunus and Malus, providing utility for understanding SI evolution. Gillenia is sister taxon to Malus, but unlike Malus, has not undergone polyploidisation. In addition, the common ancestor of Gillenia and Prunus is close to the origin of the subfamily. Using a highly contiguous Gillenia genome, orthologous regions to both Malus and Prunus S-loci were identified. Only the Prunus-like S-locus was highly polymorphic and had signatures of a functional S-locus including positive selection of the S-RNase. This suggests a self-recognition system controls SI in Gillenia, and the common ancestors of Gillenia and Prunus, and Gillenia and the apple tribe, likely had a self-recognition SI system. Comparative genomics between Gillenia and Malus suggest apple lost the self-recognition mechanism, and a nonself-recognition mechanism evolved independently from a rudimentary locus with at least one male S-determinant. Repetitive sequences in the Malus-like S-locus in Gillenia may facilitate illegitimate recombination, suggesting putative mechanisms of evolution of nonself-recognition S-loci.

Rapid detection of RNase-based self-incompatibility in Lysimachia monelli (Primulaceae)

PREMISE

Primroses famously employ a system that simultaneously expresses distyly and filters out self-pollen. Other species in the Primulaceae family, including Lysimachia monelli (blue pimpernel), also express self-incompatibility (SI), but involving a system with distinct features and an unknown molecular genetic basis.

METHODS

We utilize a candidate-based transcriptome sequencing (RNA-seq) approach, relying on candidate T2/S-RNase Class III and S-linked F-box-motif-containing genes and harnessing the unusual evolutionary and genetic features of SI, to examine whether an RNase-based mechanism underlies SI in L. monelli. We term this approach "SI detection with RNA-seq" (SIDR).

RESULTS

The results of sequencing, crossing, population genetics, and molecular evolutionary features each support a causal association linking the recovered genotypes with SI phenotypes. The finding of RNase-based SI in Primulaceae (Ericales) all but cements the long-held view that this mechanism was present in the ancestral pentapetal eudicot, whose descendants now comprise two-thirds of angiosperms. It also significantly narrows the plausible maximum age for the heterostyly evolution within the family.

CONCLUSIONS

SIDR is powerful, flexible, inexpensive, and most critically enables work in often-neglected species. It may be used with or without candidate genes to close enormous gaps in understanding the genetic basis of SI and the history of breeding system evolution.

Exploring the Allelic Diversity of the Self-Incompatibility Gene Across Natural Populations in Petunia (Solanaceae)

Self-incompatibility (SI) is a genetic mechanism to prevent self-fertilization and thereby promote outcrossing in hermaphroditic plant species through discrimination of self and nonself-pollen by pistils. In many SI systems, recognition between pollen and pistils is controlled by a single multiallelic locus (called the S-locus), in which multiple alleles (called S-alleles) are segregating. Because of the extreme level of polymorphism of the S-locus, identification of S-alleles has been a major issue in many SI studies for decades. Here, we report an RNA-seq-based method to explore allelic diversity of the S-locus by employing the long-read sequencing technology of the Oxford Nanopore MinION and applied it for the gametophytic SI system of Petunia (Solanaceae), in which the female determinant is a secreted ribonuclease called S-RNase that inhibits the elongation of self-pollen tubes by degrading RNA. We developed a method to identify S-alleles by the search of S-RNase sequences, using the previously reported sequences as queries, and found in total 62 types of S-RNase including 45 novel types. We validated this method through Sanger sequencing and crossing experiments, confirming the sequencing accuracy and SI phenotypes corresponding to genotypes. Then, using the obtained sequence data together with polymerase chain reaction-based genotyping in a larger sample set of 187 plants, we investigated the diversity, frequency, and the level of shared polymorphism of S-alleles across populations and species. The method and the dataset obtained in Petunia will be an important basis for further studying the evolution of S-RNase-based gametophytic SI systems in natural populations.

Heterozygote advantage can explain the extraordinary diversity of immune genes

The majority of highly polymorphic genes are related to immune functions and with over 100 alleles within a population, genes of the major histocompatibility complex (MHC) are the most polymorphic loci in vertebrates. How such extraordinary polymorphism arose and is maintained is controversial. One possibility is heterozygote advantage (HA), which can in principle maintain any number of alleles, but biologically explicit models based on this mechanism have so far failed to reliably predict the coexistence of significantly more than 10 alleles. We here present an eco-evolutionary model showing that evolution can result in the emergence and maintenance of more than 100 alleles under HA if the following two assumptions are fulfilled: first, pathogens are lethal in the absence of an appropriate immune defence; second, the effect of pathogens depends on host condition, with hosts in poorer condition being affected more strongly. Thus, our results show that HA can be a more potent force in explaining the extraordinary polymorphism found at MHC loci than currently recognised.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Temporal challenges in detecting balancing selection from population genomic data

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

What does effective population size tell us about loss of allelic variation?

There are two primary measures of the amount of genetic variation in a population at a locus: heterozygosity and the number of alleles. Effective population size (N e) provides both an expectation of the amount of heterozygosity in a population at drift-mutation equilibrium and the rate of loss of heterozygosity because of genetic drift. In contrast, the number of alleles in a population at drift-mutation equilibrium is a function of both N e and census size (N C). In addition, populations with the same N e can lose allelic variation at very different rates. Allelic variation is generally much more sensitive to bottlenecks than heterozygosity. Expressions used to adjust for the effects of violations of the ideal population on N e do not provide good predictions of the loss of allelic variation. These effects are much greater for loci with many alleles, which are often important for adaptation. We show that there is a linear relationship between the reduction of N C and the corresponding reduction of the expected number of alleles at drift-mutation equilibrium. This makes it possible to predict the expected effect of a bottleneck on allelic variation. Heterozygosity provides good estimates of the rate of adaptive change in the short-term, but allelic variation provides important information about long-term adaptive change. The guideline of long-term N e being greater than 500 is often used as a primary genetic metric for evaluating conservation status. We recommend that this guideline be expanded to take into account allelic variation as well as heterozygosity.

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

Absence of long-term balancing selection on variation in EuMYB3, an R2R3-MYB gene responsible for the anther-color polymorphism in Erythronium umbilicatum

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

Gametophytic self-incompatibility in Maxillariinae orchids

Gametophytic self-incompatibility (GSI) has been mainly described in species-rich clades within Orchidaceae. We report GSI for a genus within Maxillariinae, one of the most conspicuous and diverse subtribes of neotropical orchids. Here, we describe the reproductive system of Brasiliorchis picta, B. phoenicanthera, and B. porphyrostele. Anatomical studies of fruit development showed that pollen tubes of aborted, self-pollinated flowers degenerate half-way in the stylar channel and never reach the ovules. Spontaneous self-pollination and emasculation set no fruits for none of the sampled species, thus supporting the hypothesis that these three species studied rely on the agency of pollinators and pollinator-mediated cross-pollination to set fruit. Fruit set from cross-pollinations ranged from 33.4 to 77.5%. One self-pollinated fruit of B. porphyrostele developed to completion. All other fruits aborted between 10 and 21 days after pollination. These data support previous evidence of variable strength GSI being exhibited in orchid species. Additional studies of self-incompatibility systems are needed to evaluate their role in species diversification and evolution of reproductive strategies in Maxillariinae and to allow for effective conservation strategies of threatened orchid species.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence

Self-incompatibility in Prunus species is governed by a single locus consisting of two highly multi-allelic and tightly linked genes, one coding for an F-box protein-i.e., SFB in Prunus- controlling the pollen specificity and one coding for an S-RNase gene controlling the pistil specificity. Genotyping the allelic combination in a fruit tree species is an essential procedure both for cross-based breeding and for establishing pollination requirements. Gel-based PCR techniques using primer pairs designed from conserved regions and spanning polymorphic intronic regions are traditionally used for this task. However, with the great advance of massive sequencing techniques and the lowering of sequencing costs, new genotyping-by-sequencing procedures are emerging. The alignment of resequenced individuals to reference genomes, commonly used for polymorphism detection, yields little or no coverage in the S-locus region due to high polymorphism between different alleles within the same species, and cannot be used for this purpose. Using the available sequences of Japanese plum S-loci concatenated in a rosary-like structure as synthetic reference sequence, we describe a procedure to accurately genotype resequenced individuals that allowed the analysis of the S-genotype in 88 Japanese plum cultivars, 74 of them are reported for the first time. In addition to unraveling two new S-alleles from published reference genomes, we identified at least two S-alleles in 74 cultivars. According to their S-allele composition, they were assigned to 22 incompatibility groups, including nine new incompatibility groups reported here for the first time (XXVII-XXXV).

Large-scale fungal strain sequencing unravels the molecular diversity in mating loci maintained by long-term balancing selection

Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of strains. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci define mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty strains of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. The sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. 17,550 mating types were predicted based on the allelic classes. In vitro crosses allowed us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. We conclude that the genetic diversity of mating loci in Trichaptum is due to long-term balancing selection, with limited recombination and duplication activity. The large number of sequenced strains highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Urbanization increases floral specialization of pollinators

Understanding how urbanization alters functional interactions among pollinators and plants is critically important given increasing anthropogenic land use and declines in pollinator populations. Pollinators often exhibit short-term specialization and visit plants of the same species during one foraging trip. This facilitates plant receipt of conspecific pollen-pollen on a pollinator that is the same species as the plant on which the pollinator was foraging. Conspecific pollen receipt facilitates plant reproductive success and is thus important to plant and pollinator persistence. We investigated how urbanization affects short-term specialization of insect pollinators by examining pollen loads on insects' bodies and identifying the number and species of pollen grains on insects caught in urban habitat fragments and natural areas. We assessed possible drivers of differences between urban and natural areas, including frequency dependence in foraging, species richness and diversity of the plant and pollinator communities, floral abundance, and the presence of invasive plant species. Pollinators were more specialized in urban fragments than in natural areas, despite no differences in the species richness of plant communities across site types. These differences were likely driven by higher specialization of common pollinators, which were more abundant in urban sites. In addition, pollinators preferred to forage on invasive plants at urban sites and native plants at natural sites. Our findings reveal indirect effects of urbanization on pollinator fidelity to individual plant species and have implications for the maintenance of plant species diversity in small habitat fragments. Higher preference of pollinators for invasive plants at urban sites suggests that native species may receive fewer visits by pollinators. Therefore, native plant species diversity may decline in urban sites without continued augmentation of urban flora or removal of invasive species.

Frequency dependence and the predictability of evolution in a changing environment

Frequency‐dependent (FD) selection, whereby fitness and selection depend on the genetic or phenotypic composition of the population, arises in numerous ecological contexts (competition, mate choice, crypsis, mimicry, etc.) and can strongly impact evolutionary dynamics. In particular, negative frequency‐dependent selection (NFDS) is well known for its ability to potentially maintain stable polymorphisms, but it has also been invoked as a source of persistent, predictable frequency fluctuations. However, the conditions under which such fluctuations persist are not entirely clear. In particular, previous work rarely considered that FD is unlikely to be the sole driver of evolutionary dynamics when it occurs, because most environments are not static but instead change dynamically over time. Here, we investigate how FD interacts with a temporally fluctuating environment to shape the dynamics of population genetic change. We show that a simple metric introduced by Lewontin, the slope of frequency change against frequency near equilibrium, works as a key criterion for distinguishing microevolutionary outcomes, even in a changing environment. When this slope D is between 0 and –2 (consistent with the empirical examples we review), substantial fluctuations would not persist on their own in a large population occupying a constant environment, but they can still be maintained indefinitely as quasi‐cycles fueled by environmental noise or genetic drift. However, such moderate NFDS buffers and temporally shifts evolutionary responses to periodic environments (e.g., seasonality). Stronger FD, with slope D < –2, can produce self‐sustained cycles that may overwhelm responses to a changing environment, or even chaos that fundamentally limits predictability. This diversity of expected outcomes, together with the empirical evidence for both FD and environment‐dependent selection, suggests that the interplay of internal dynamics with external forcing should be investigated more systematically to reach a better understanding and prediction of evolution.

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.

Mates Matter: Gametophyte Kinship Recognition and Inbreeding in the Giant Kelp, Macrocystis pyrifera (Laminariales, Phaeophyceae)

Inbreeding, the mating between genetically related individuals, often results in reduced survival and fecundity of offspring, relative to outcrossing. Yet, high inbreeding rates are commonly observed in seaweeds, suggesting compensatory reproductive traits may affect the costs and benefits of the mating system. We experimentally manipulated inbreeding levels in controlled crossing experiments, using gametophytes from 19 populations of Macrocystis pyrifera along its Eastern Pacific coastal distribution (EPC). The objective was to investigate the effects of male-female kinship on female fecundity and fertility, to estimate inbreeding depression in the F1 progeny, and to assess the variability of these effects among different regions and habitats of the EPC. Results revealed that the presence and kinship of males had a significant effect on fecundity and fertility of female gametophytes. Females left alone or in the presence of sibling males express the highest gametophyte size, number, and size of oogonia, suggesting they were able to sense the presence and the identity of their mates before gamete contact. The opposite trend was observed for the production of embryos per female gametes, indicating higher costs of selfing and parthenogenesis than outcrossing on fertility. However, the increased fecundity compensated for the reduced fertility, leading to a stable overall reproductive output. Inbreeding also affected morphological traits of juvenile sporophytes, but not their heatwave tolerance. The male-female kinship effect was stronger in high-latitude populations, suggesting that females from low-latitude marginal populations might have evolved to mate with any male gamete to guarantee reproductive success.

Studying models of balancing selection using phase-type theory

Balancing selection (BLS) is the evolutionary force that maintains high levels of genetic variability in many important genes. To further our understanding of its evolutionary significance, we analyze models with BLS acting on a biallelic locus: an equilibrium model with long-term BLS, a model with long-term BLS and recent changes in population size, and a model of recent BLS. Using phase-type theory, a mathematical tool for analyzing continuous time Markov chains with an absorbing state, we examine how BLS affects polymorphism patterns in linked neutral regions, as summarized by nucleotide diversity, the expected number of segregating sites, the site frequency spectrum, and the level of linkage disequilibrium (LD). Long-term BLS affects polymorphism patterns in a relatively small genomic neighborhood, and such selection targets are easier to detect when the equilibrium frequencies of the selected variants are close to 50%, or when there has been a population size reduction. For a new mutation subject to BLS, its initial increase in frequency in the population causes linked neutral regions to have reduced diversity, an excess of both high and low frequency derived variants, and elevated LD with the selected locus. These patterns are similar to those produced by selective sweeps, but the effects of recent BLS are weaker. Nonetheless, compared to selective sweeps, nonequilibrium polymorphism and LD patterns persist for a much longer period under recent BLS, which may increase the chance of detecting such selection targets. An R package for analyzing these models, among others (e.g., isolation with migration), is available.

Confirmation of a Gametophytic Self-Incompatibility in Oryza longistaminata

Oryza longistaminata, a wild species of African origin, has been reported to exhibit self-incompatibility (SI). However, the genetic pattern of its SI remained unknown. In this study, we conducted self-pollination and reciprocal cross-pollination experiments to verify that O. longistaminata is a strictly self-incompatible species. The staining of pollen with aniline blue following self-pollination revealed that although pollen could germinate on the stigma, the pollen tube was unable to enter the style to complete pollination, thereby resulting in gametophytic self-incompatibility (GSI). LpSDUF247, a S-locus male determinant in the gametophytic SI system of perennial ryegrass, is predicted to encode a DUF247 protein. On the basic of chromosome alignment with LpSDUF247, we identified OlSS1 and OlSS2 as Self-Incompatibility Stamen candidate genes in O. longistaminata. Chromosome segment analysis revealed that the Self-Incompatibility Pistil candidate gene of O. longistaminata (OlSP) is a polymorphic gene located in a region flanking OlSS1. OlSS1 was expressed mainly in the stamens, whereas OlSS2 was expressed in both the stamens and pistils. OlSP was specifically highly expressed in the pistils, as revealed by RT-PCR and qRT-PCR analyses. Collectively, our observations indicate the occurrence of GSI in O. longistaminata and that this process is potentially controlled by OlSS1, OlSS2, and OlSP. These findings provide further insights into the genetic mechanisms underlying self-compatibility in plants.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Population Genomics of the "Arcanum" Species Group in Wild Tomatoes: Evidence for Separate Origins of Two Self-Compatible Lineages

Given their diverse mating systems and recent divergence, wild tomatoes (Solanum section Lycopersicon) have become an attractive model system to study ecological divergence, the build-up of reproductive barriers, and the causes and consequences of the breakdown of self-incompatibility. Here we report on a lesser-studied group of species known as the "Arcanum" group, comprising the nominal species Solanum arcanum, Solanum chmielewskii, and Solanum neorickii. The latter two taxa are self-compatible but are thought to self-fertilize at different rates, given their distinct manifestations of the morphological "selfing syndrome." Based on experimental crossings and transcriptome sequencing of a total of 39 different genotypes from as many accessions representing each species' geographic range, we provide compelling evidence for deep genealogical divisions within S. arcanum; only the self-incompatible lineage known as "var. marañón" has close genealogical ties to the two self-compatible species. Moreover, there is evidence under multiple inference schemes for different geographic subsets of S. arcanum var. marañón being closest to S. chmielewskii and S. neorickii, respectively. To broadly characterize the population-genomic consequences of these recent mating-system transitions and their associated speciation events, we fit demographic models indicating strong reductions in effective population size, congruent with reduced nucleotide and S-locus diversity in the two independently derived self-compatible species.

Molecular and Evolutionary Characterization of Pollen S Determinant (SFB Alleles) in Four Diploid and Hexaploid Plum Species (Prunus spp.)

In more than 60 families of angiosperms, the self- and cross-fertilization is avoided through a complex widespread genetic system called self-incompatibility (SI). One of the major puzzling issues concerning the SI is the evolution of this system in species with complex polyploid genomes. Among plums, one of the first fruits species to attract human interest, polyploid species represent enormous genetic potential, which can be exploited in breeding programs. However, molecular studies in these species are very scarce due to the complexity of their genome. In order to study the SFB gene [the male component of gametophytic self-incompatibility system (GSI)] in plum species, 36 plum accessions belonging to diploid and hexaploid species were used. A total of 19 different alleles were identified; 1 of them was revealed after analyzing sequences. Peptide sequence analysis allowed identifying the five domains features of the SFB gene. Polymorphism analysis showed a subtle difference between domesticated and open pollinated Tunisian accessions and suggested a probable influence of the ploidy level. Divergence analysis between studied sequences showed that a new specificity may appear after 5.3% of divergence at synonymous sites between pairs of sequences in Prunus insititia, 6% in Prunus cerasifera, 8% and 9% in Prunus domestica and Prunus salicina respectively. Furthermore, sites under positive selection, the ones more likely to be responsible for specificity determination, were identified. A positive and significant Pearson correlation was found between the divergence between sequences, divergence time, fixed substitutions (MK test), and PSS number. These results supported the model assuming that functionally distinct proteins have arisen not as a result of chance fixation of neutral variants, but rather as a result of positive Darwinian selection. Further, the role that plays recombination can not be ruled out, since a rate of 0.08 recombination event per polymorphic sites was identified.

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.

Base-Pairing Requirements for Small RNA-Mediated Gene Silencing of Recessive Self-Incompatibility Alleles in Arabidopsis halleri

Small noncoding RNAs are central regulators of genome activity and stability. Their regulatory function typically involves sequence similarity with their target sites, but understanding the criteria by which they specifically recognize and regulate their targets across the genome remains a major challenge in the field, especially in the face of the diversity of silencing pathways involved. The dominance hierarchy among self-incompatibility alleles in Brassicaceae is controlled by interactions between a highly diversified set of small noncoding RNAs produced by dominant S-alleles and their corresponding target sites on recessive S-alleles. By controlled crosses, we created numerous heterozygous combinations of S-alleles in Arabidopsis halleri and developed an real-time quantitative PCR assay to compare allele-specific transcript levels for the pollen determinant of self-incompatibility (SCR). This provides the unique opportunity to evaluate the precise base-pairing requirements for effective transcriptional regulation of this target gene. We found strong transcriptional silencing of recessive SCR alleles in all heterozygote combinations examined. A simple threshold model of base pairing for the small RNA-target interaction captures most of the variation in SCR transcript levels. For a subset of S-alleles, we also measured allele-specific transcript levels of the determinant of pistil specificity (SRK), and found sharply distinct expression dynamics throughout flower development between SCR and SRK In contrast to SCR, both SRK alleles were expressed at similar levels in the heterozygote genotypes examined, suggesting no transcriptional control of dominance for this gene. We discuss the implications for the evolutionary processes associated with the origin and maintenance of the dominance hierarchy among self-incompatibility alleles.

Negative frequency dependent selection contributes to the maintenance of a global polymorphism in mitochondrial DNA

BACKGROUND

Understanding the forces that maintain diversity across a range of scales is at the very heart of biology. Frequency-dependent processes are generally recognized as the most central process for the maintenance of ecological diversity. The same is, however, not generally true for genetic diversity. Negative frequency dependent selection, where rare genotypes have an advantage, is often regarded as a relatively weak force in maintaining genetic variation in life history traits because recombination disassociates alleles across many genes. Yet, many regions of the genome show low rates of recombination and genetic variation in such regions (i.e., supergenes) may in theory be upheld by frequency dependent selection.

RESULTS

We studied what is essentially a ubiquitous life history supergene (i.e., mitochondrial DNA) in the fruit fly Drosophila subobscura, showing sympatric polymorphism with two main mtDNA genotypes co-occurring in populations world-wide. Using an experimental evolution approach involving manipulations of genotype starting frequencies, we show that negative frequency dependent selection indeed acts to maintain genetic variation in this region. Moreover, the strength of selection was affected by food resource conditions.

CONCLUSIONS

Our work provides novel experimental support for the view that balancing selection through negative frequency dependency acts to maintain genetic variation in life history genes. We suggest that the emergence of negative frequency dependent selection on mtDNA is symptomatic of the fundamental link between ecological processes related to resource use and the maintenance of genetic variation.

Evolution of self-compatibility by a mutant Sm-RNase in citrus

Self-incompatibility (SI) is an important mechanism that prevents self-fertilization and inbreeding in flowering plants. The most widespread SI system utilizes S ribonucleases (S-RNases) and S-locus F-boxes (SLFs) as S determinants. In citrus, SI is ancestral, and Citrus maxima (pummelo) is self-incompatible, while Citrus reticulata (mandarin) and its hybrids are self-compatible (SC). Here, we identify nine highly polymorphic pistil-specific, developmentally expressed S-RNases from pummelo that segregate with S haplotypes in a gametophytic manner and cluster with authentic S-RNases. We provide evidence that these S-RNases function as the female S determinants in citrus. Moreover, we show that each S-RNase is linked to approximately nine SLFs. In an analysis of 117 citrus SLF and SFL-like (SLFL) genes, we reveal that they cluster into 12 types and that the S-RNases and intra-haplotypic SLF and SLFL genes co-evolved. Our data support the notion that citrus have a S locus comprising a S-RNase and several SLFs that fit the non-self-recognition model. We identify a predominant single nucleotide mutation, S_m-RNase, in SC citrus, which provides a 'natural' loss of function. We show that SI-SC transitions due to the S_m-RNase initially arose in mandarin, spreading to its hybrids and became fixed. Identification of an evolutionarily distant new genus utilizing the S-RNase-based SI system, >100 million years separated from the nearest S-RNase family, is a milestone for evolutionary comparative studies.

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.

Genetic basis and timing of a major mating system shift in Capsella

A crucial step in the transition from outcrossing to self-fertilization is the loss of genetic self-incompatibility (SI). In the Brassicaceae, SI involves the interaction of female and male specificity components, encoded by the genes SRK and SCR at the self-incompatibility locus (S-locus). Theory predicts that S-linked mutations, and especially dominant mutations in SCR, are likely to contribute to loss of SI. However, few studies have investigated the contribution of dominant mutations to loss of SI in wild plant species. Here, we investigate the genetic basis of loss of SI in the self-fertilizing crucifer species Capsella orientalis, by combining genetic mapping, long-read sequencing of complete S-haplotypes, gene expression analyses and controlled crosses. We show that loss of SI in C. orientalis occurred < 2.6 Mya and maps as a dominant trait to the S-locus. We identify a fixed frameshift deletion in the male specificity gene SCR and confirm loss of male SI specificity. We further identify an S-linked small RNA that is predicted to cause dominance of self-compatibility. Our results agree with predictions on the contribution of dominant S-linked mutations to loss of SI, and thus provide new insights into the molecular basis of mating system transitions.

Nonlinear frequency-dependent selection promotes long-term coexistence between bacteria species

Negative frequency-dependent selection (NFDS) is an important mechanism for species coexistence and for the maintenance of genetic polymorphism. Long-term coexistence nevertheless requires NFDS interactions to be resilient to further evolution of the interacting species or genotypes. For closely related genotypes, NFDS interactions have been shown to be preserved through successive rounds of evolution in coexisting lineages. On the contrary, the evolution of NFDS interactions between distantly related species has received less attention. Here, we tracked the co-evolution of Escherichia coli and Citrobacter freundii that initially differ in their ecological characteristics. We showed that these two bacterial species engaged in an NFDS interaction particularly resilient to further evolution: despite a very strong asymmetric rate of adaptation, their coexistence was maintained owing to an NFDS pattern where fitness increases steeply as the frequency decreases towards zero. Using a model, we showed how and why such NFDS pattern can emerge. These findings provide a robust explanation for the long-term maintenance of species at very low frequencies.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Increased fluctuation in a butterfly metapopulation leads to diploid males and decline of a hyperparasitoid

Climate change can increase spatial synchrony of population dynamics, leading to large-scale fluctuation that destabilizes communities. High trophic level species such as parasitoids are disproportionally affected because they depend on unstable resources. Most parasitoid wasps have complementary sex determination, producing sterile males when inbred, which can theoretically lead to population extinction via the diploid male vortex (DMV). We examined this process empirically using a hyperparasitoid population inhabiting a spatially structured host population in a large fragmented landscape. Over four years of high host butterfly metapopulation fluctuation, diploid male production by the wasp increased, and effective population size declined precipitously. Our multitrophic spatially structured model shows that host population fluctuation can cause local extinctions of the hyperparasitoid because of the DMV. However, regionally it persists because spatial structure allows for efficient local genetic rescue via balancing selection for rare alleles carried by immigrants. This is, to our knowledge, the first empirically based study of the possibility of the DMV in a natural host–parasitoid system.

An Unexplored Side of Regeneration Niche: Seed Quantity and Quality Are Determined by the Effect of Temperature on Pollen Performance

In 1977, Peter Grubb introduced the regeneration niche concept, which assumes that a plant species cannot persist if the environmental conditions are only suitable for adult plant growth and survival, but not for seed production, dispersal, germination, and seedling establishment. During the last decade, this concept has received considerable research attention as it helps to better understand community assembly, population dynamics, and plant responses to environmental changes. Yet, in its present form, it focuses too much on the post-fertilization stages of plant sexual reproduction, neglecting the fact that the environment can operate as a constraint at many points in the chain of processes necessary for successful regeneration. In this review, we draw the attention of the plant ecology research community to the pre-fertilization stages of plant sexual reproduction, an almost ignored but important aspect of the regeneration niche, and their potential consequences for successful seed production. Particularly, we focus on how temperature affects pollen performance and determines plant reproduction success by playing an important role in the temporal and spatial variations in seed quality and quantity. We also review the pollen adaptations to temperature stresses at different levels of plant organization and discuss the plasticity of the performance of pollen under changing temperature conditions. The reviewed literature demonstrates that pre-fertilization stages of seed production, particularly the extreme sensitivity of male gametophyte performance to temperature, are the key determinants of a species' regeneration niche. Thus, we suggest that previous views stating that the regeneration niche begins with the production of seeds should be modified to include the preceding stages. Lastly, we identify several gaps in pollen-related studies revealing a framework of opportunities for future research, particularly how these findings could be used in the field of plant biology and ecology.

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.

Negative Frequency-Dependent Selection Is Frequently Confounding

Persistent genetic variation within populations presents an evolutionary problem, as natural selection and genetic drift tend to erode genetic diversity. Models of balancing selection were developed to account for the maintenance of genetic variation observed in natural populations. Negative frequency-dependent selection is a powerful type of balancing selection that maintains many natural polymorphisms, but it is also commonly misinterpreted. This review aims to clarify the processes underlying negative frequency-dependent selection, describe classes of polymorphisms that can and cannot result from these processes, and discuss the empirical data needed to accurately identify processes that generate or maintain diversity in nature. Finally, the importance of accurately describing the processes affecting genetic diversity within populations as it relates to research progress is considered.

Rarity and persistence

Rarity is a population characteristic that is usually associated with a high risk of extinction. We argue here, however, that chronically rare species (those with low population densities over many generations across their entire ranges) may have individual-level traits that make populations more resistant to extinction. The major obstacle to persistence at low density is successful fertilisation (union between egg and sperm), and chronically rare species are more likely to survive when (1) fertilisation occurs inside or close to an adult, (2) mate choice involves long-distance signals, (3) adults or their surrogate gamete dispersers are highly mobile, or (4) the two sexes are combined in a single individual. In contrast, external fertilisation and wind- or water-driven passive dispersal of gametes, or sluggish or sedentary adult life habits in the absence of gamete vectors, appear to be incompatible with sustained rarity. We suggest that the documented increase in frequency of these traits among marine genera over geological time could explain observed secular decreases in rates of background extinction. Unanswered questions remain about how common chronic rarity actually is, which traits are consistently associated with chronic rarity, and how chronically rare species are distributed among taxa, and among the world's ecosystems and regions.

Balancing selection on the number of repeats in the ribosomal intergenic spacer present in naturally occurring yellow perch (Perca flavescens) populations

The ribosomal intergenic spacer (IGS), responsible for the rate of transcription of rRNA genes, is associated with the growth and fecundity of individuals. A previous study of IGS length variants in a yellow perch (Perca flavescens) population revealed the presence of two predominant alleles differing by 1 kb due to variation in the number of repeat units. This study aims to assess whether length variation of IGS is the result of selection in natural populations. Length variation of IGS and 11 neutral microsatellite loci were assessed in geographically distant yellow perch populations. Most populations displayed the very same IGS alleles; they did not differ in frequencies among populations and the F_ST was not significantly different from zero. In contrast, diversity at microsatellite loci was high and differed among populations (F_ST = 0.18). Selection test based on F_ST identified IGS as a significant outlier from neutral expectations for population differentiation. Heterozygote excess was also detected in one specific cohort, suggesting temporal variation in the selection regime. While the exact mechanism remains to be specified, together the results of this study support the contention that balancing selection is acting to maintain two distinct IGS alleles in natural fish populations.

Plant mating systems: self-incompatibility and evolutionary transitions to self-fertility in the mustard family

Flowering plants have evolved diverse mechanisms that promote outcrossing. The most widespread of these outbreeding devices are self-incompatibility systems, the highly selective prefertilization mating barriers that prevent self-fertilization by disrupting pollen-pistil interactions. Despite the advantages of outcrossing, loss of self-incompatibility has occurred repeatedly in many plant families. In the mustard family, the highly polymorphic receptors and ligands that mediate the recognition and inhibition of self-pollen in self-incompatibility have been characterized and the 3D structure of the receptor-ligand complex has been solved. Sequence analyses and empirical studies in self-incompatible and self-compatible species are elucidating the genetic basis of switches from the outcrossing to selfing modes of mating and beginning to provide clues to the diversification of the self recognition repertoire.

Genome Sequencing Reveals the Origin of the Allotetraploid Arabidopsis suecica

Polyploidy is an example of instantaneous speciation when it involves the formation of a new cytotype that is incompatible with the parental species. Because new polyploid individuals are likely to be rare, establishment of a new species is unlikely unless polyploids are able to reproduce through self-fertilization (selfing), or asexually. Conversely, selfing (or asexuality) makes it possible for polyploid species to originate from a single individual-a bona fide speciation event. The extent to which this happens is not known. Here, we consider the origin of Arabidopsis suecica, a selfing allopolyploid between Arabidopsis thaliana and Arabidopsis arenosa, which has hitherto been considered to be an example of a unique origin. Based on whole-genome re-sequencing of 15 natural A. suecica accessions, we identify ubiquitous shared polymorphism with the parental species, and hence conclusively reject a unique origin in favor of multiple founding individuals. We further estimate that the species originated after the last glacial maximum in Eastern Europe or central Eurasia (rather than Sweden, as the name might suggest). Finally, annotation of the self-incompatibility loci in A. suecica revealed that both loci carry non-functional alleles. The locus inherited from the selfing A. thaliana is fixed for an ancestral non-functional allele, whereas the locus inherited from the outcrossing A. arenosa is fixed for a novel loss-of-function allele. Furthermore, the allele inherited from A. thaliana is predicted to transcriptionally silence the allele inherited from A. arenosa, suggesting that loss of self-incompatibility may have been instantaneous.

Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures

Explaining the maintenance of adaptive diversity within populations is a long-standing goal in evolutionary biology, with important implications for conservation, medicine, and agriculture. Adaptation often leads to the fixation of beneficial alleles, and therefore it erodes local diversity so that understanding the coexistence of multiple adaptive phenotypes requires deciphering the ecological mechanisms that determine their respective benefits. Here, we show how antagonistic frequency-dependent selection (FDS), generated by natural and sexual selection acting on the same trait, maintains mimicry polymorphism in the toxic butterfly Heliconius numata Positive FDS imposed by predators on mimetic signals favors the fixation of the most abundant and best-protected wing-pattern morph, thereby limiting polymorphism. However, by using mate-choice experiments, we reveal disassortative mate preferences of the different wing-pattern morphs. The resulting negative FDS on wing-pattern alleles is consistent with the excess of heterozygote genotypes at the supergene locus controlling wing-pattern variation in natural populations of H. numata The combined effect of positive and negative FDS on visual signals is sufficient to maintain a diversity of morphs displaying accurate mimicry with other local prey, although some of the forms only provide moderate protection against predators. Our findings help understand how alternative adaptive phenotypes can be maintained within populations and emphasize the need to investigate interactions between selective pressures in other cases of puzzling adaptive polymorphism.

Identification, genealogical structure and population genetics of S-alleles in Malus sieversii, the wild ancestor of domesticated apple

The self-incompatibility (SI) gene that is specifically expressed in pistils encodes the SI-associated ribonuclease (S-RNase), functioning as the female-specificity determinant of a gametophytic SI system. Despite extensive surveys in Malus domestica, the S-alleles have not been fully investigated for Malus sieversii, the primary wild ancestor of the domesticated apple. Here we screened the M. sieversii S-alleles via PCR amplification and sequencing, and identified 14 distinct alleles in this species. By contrast, nearly 40 are present in its close wild relative, Malus sylvestris. We further sequenced 8 nuclear genes to provide a neutral reference, and investigated the evolution of S-alleles via genealogical and population genetic analyses. Both shared ancestral polymorphism and an excess of non-synonymous substitution were detected in the S-RNases of the tribe Maleae in Rosaceae, indicating the action of long-term balancing selection. Approximate Bayesian Computations based on the reference neutral loci revealed a severe bottleneck in four of the six studied M. sieversii populations, suggesting that the low number of S-alleles found in this species is mainly the result of diversity loss due to a drastic population contraction. Such a bottleneck may lead to ambiguous footprints of ongoing balancing selection detected at the S-locus. This study not only elucidates the constituents and number of S-alleles in M. sieversii but also illustrates the potential utility of S-allele number shifts in demographic inference for self-incompatible plant species.