PHYLOGENY OF LEAVENWORTHIA S-ALLELES SUGGESTS UNIDIRECTIONAL MATING SYSTEM EVOLUTION AND ENHANCED POSITIVE SELECTION FOLLOWING AN ANCIENT POPULATION BOTTLENECK

Self-incompatibility and the genetic architecture of inbreeding depression

Inbreeding depression plays a fundamental role in evolution. To help detect and characterize the loci that underlie inbreeding depression, we used bud pollination and salt treatments to circumvent self-incompatibility (SI) in plants from populations of Leavenworthia alabamica and produced families of progeny that were then genotyped at genetically mapped single-nucleotide polymorphism (SNP) loci. Using Bayesian inference, the segregation patterns for each SNP were used to explore support for different dominance and selection coefficients at linked viability loci in different genomic regions. There was support for several partially recessive viability loci in one of the populations, and one such locus mapped to the genomic region of the novel SI locus in L. alabamica. These results are consistent with earlier findings that showed purging of inbreeding depression for germination rate in L. alabamica. They are also consistent with expectations from evolutionary genetic theory that recessive, deleterious alleles linked to loci under balancing selection can be sheltered from selection.

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.

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.

Genetic variation for pseudo-self-compatibility in self-incompatible populations of Leavenworthia alabamica (Brassicaceae)

Self-incompatibility (SI) promotes outcrossing, but transitions to self-compatibility (SC) are frequent. Population genetic theory describing the breakdown of SI to SC suggests that, under most conditions, populations should be composed of either SI or SC individuals. Under a narrow range of conditions, theory suggests that SI may persist alongside reduced expression of SI (pseudo-SI, PSI) in mixed-mating populations. We studied genetic variation for PSI segregating in four SI populations of Leavenworthia alabamica by measurement of the heritability of pollen tube number after self-pollination. We tested for the role of the S-locus in this variation by sequencing seven S-alleles from plants with high pseudo-SC (PSC) and testing for the co-segregation of these alleles with PSC. We found a continuous distribution of PSC in all populations and 90% of plants exhibited PSC. The heritability ranged from 0.39 to 0.57. All seven S-alleles from plants with high PSC exhibited trans-specific polymorphism, and no stop codons were observed within the c. 600-bp region sequenced. One of these S-alleles was directly associated with the inheritance of PSC. We conclude that heritable variation in PSC is largely a result of genetic variation in the signaling cascade downstream of the S-locus reaction, together with the presence of one leaky S-allele.

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.

Recent selection for self-compatibility in a population of Leavenworthia alabamica

The evolution of self-compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self-incompatibility (SI) to self-fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S-locus. Despite permitting inbreeding, mutations at the S-locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self- and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S-locus, we collected sequence data from a ∼690 Kb region surrounding the S-locus, as well as from regions not linked to the S-locus. To test for recent positive selection acting at the S-locus, we examined polymorphism and the site-frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.

Evolution of Selfing: Recurrent Patterns in Molecular Adaptation

Selfing has evolved in animals, fungi, and plants, and since Darwin's pioneering study, it is considered one of the most frequent evolutionary trends in flowering plants. Generally, the evolution of selfing is characterized by a loss of self-incompatibility, the selfing syndrome, and changes in genome-wide polymorphism patterns. Recent interdisciplinary studies involving molecular functional experiments, genome-wide data, experimental evolution, and evolutionary ecology using Arabidopsis thaliana, Caenorhabditis elegans, and other species show that the evolution of selfing is not merely a degradation of outcrossing traits but a model for studying the recurrent patterns underlying adaptive molecular evolution. For example, in wild Arabidopsis relatives, self-compatibility evolved from mutations in the male specificity gene, S-LOCUS CYSTEINE-RICH PROTEIN/S-LOCUS PROTEIN 11 (SCR/SP11), rather than the female specificity gene, S-LOCUS RECEPTOR KINASE (SRK), supporting the theoretical prediction of sexual asymmetry. Prevalence of dominant self-compatible mutations is consistent with Haldane's sieve, which acts against recessive adaptive mutations. Time estimates based on genome-wide polymorphisms and self-incompatibility genes generally support the recent origin of selfing.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

BACKGROUND

The transition from outcrossing to selfing has long been portrayed as an 'evolutionary dead end' because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a 'rostellum') co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.

RESULTS

The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4-3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.

CONCLUSIONS

The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the 'rostellum') is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.

Convergent evolution at the gametophytic self-incompatibility system in Malus and Prunus

S-RNase-based gametophytic self-incompatibility (GSI) has evolved once before the split of the Asteridae and Rosidae. This conclusion is based on the phylogenetic history of the S-RNase that determines pistil specificity. In Rosaceae, molecular characterizations of Prunus species, and species from the tribe Pyreae (i.e., Malus, Pyrus, Sorbus) revealed different numbers of genes determining S-pollen specificity. In Prunus only one pistil and pollen gene determine GSI, while in Pyreae there is one pistil but multiple pollen genes, implying different specificity recognition mechanisms. It is thus conceivable that within Rosaceae the genes involved in GSI in the two lineages are not orthologous but possibly paralogous. To address this hypothesis we characterised the S-RNase lineage and S-pollen lineage genes present in the genomes of five Rosaceae species from three genera: M. × domestica (apple, self-incompatible (SI); tribe Pyreae), P. persica (peach, self-compatible (SC); Amygdaleae), P. mume (mei, SI; Amygdaleae), Fragaria vesca (strawberry, SC; Potentilleae), and F. nipponica (mori-ichigo, SI; Potentilleae). Phylogenetic analyses revealed that the Malus and Prunus S-RNase and S-pollen genes belong to distinct gene lineages, and that only Prunus S-RNase and SFB-lineage genes are present in Fragaria. Thus, S-RNase based GSI system of Malus evolved independently from the ancestral system of Rosaceae. Using expression patterns based on RNA-seq data, the ancestral S-RNase lineage gene is inferred to be expressed in pistils only, while the ancestral S-pollen lineage gene is inferred to be expressed in tissues other than pollen.

Reproductive character displacement shapes a spatially structured petal color polymorphism in Leavenworthia stylosa

Character displacement is a potentially important process driving trait evolution and species diversification. Floral traits may experience character displacement in response to pollinator-mediated competition (ecological character displacement) or the risk of forming hybrids with reduced fitness (reproductive character displacement). We test these and alternative hypotheses to explain a yellow-white petal color polymorphism in Leavenworthia stylosa, where yellow morphs are spatially associated with a white-petaled congener (Leavenworthia exigua) that produces hybrids with complete pollen sterility. A reciprocal transplant experiment found limited evidence of local adaptation of yellow color morphs via increased survival and seed set. Pollinator observations revealed that Leavenworthia attract various pollinators that generally favor white petals and exhibit color constancy. Pollen limitation experiments showed that yellow petals do not alleviate competition for pollination. Interspecific pollinator movements were infrequent and low hybridization rates (∼0.40-0.85%) were found in each morph, with natural rates likely being lower. Regardless, hybridization rates were significantly higher in white morphs of L. stylosa, yielding a small selection coefficient of s = 0.0042 against this phenotype in sympatry with L. exigua. These results provide support for RCD as a mechanism contributing to the pattern of petal color polymorphism in L. stylosa.

Self-incompatibility in Brassicaceae: identification and characterization of SRK-like sequences linked to the S-locus in the tribe Biscutelleae

Self-incompatibility (SI) is a genetic system that prevents self-fertilization in many Angiosperms. Although plants from the Brassicaceae family present an apparently unique SI system that is ancestral to the family, investigations at the S-locus responsible for SI have been mostly limited to two distinct lineages (Brassica and Arabidopsis-Capsella, respectively). Here, we investigated SI in a third deep-branching lineage of Brassicaceae: the tribe Biscutelleae. By coupling sequencing of the SI gene responsible for pollen recognition (SRK) with phenotypic analyses based on controlled pollinations, we identified 20 SRK-like sequences functionally linked to 13 S-haplotypes in 21 individuals of Biscutella neustriaca and 220 seedlings. We found two genetic and phylogenetic features of SI in Biscutelleae that depart from patterns observed in the reference Arabidopsis clade: (1) SRK-like sequences cluster into two main phylogenetic lineages interspersed within the many SRK lineages of Arabidopsis; and (2) some SRK-like sequences are transmitted by linked pairs, suggesting local duplication within the S-locus. Strikingly, these features also were observed in the Brassica clade but probably evolved independently, as the two main SRK clusters in Biscutella are distinct from those in Brassica. In the light of our results and of what has been previously observed in other Brassicaceae, we discuss the ecological and evolutionary implications on SI plant populations of the high diversity and the complex dominance relationships we found at the S-locus in Biscutelleae.

Limited genomic consequences of mixed mating in the recently derived sister species pair, Collinsia concolor and Collinsia parryi

Highly selfing species often show reduced effective population sizes and reduced selection efficacy. Whether mixed mating species, which produce both self and outcross progeny, show similar patterns of diversity and selection remains less clear. Examination of patterns of molecular evolution and levels of diversity in species with mixed mating systems can be particularly useful for investigating the relative importance of linked selection and demographic effects on diversity and the efficacy of selection, as the effects of linked selection should be minimal in mixed mating populations, although severe bottlenecks tied to founder events could still be frequent. To begin to address this gap, we assembled and analysed the transcriptomes of individuals from a recently diverged mixed mating sister species pair in the self-compatible genus, Collinsia. The de novo assembly of 52 and 37 Mbp C. concolor and C. parryi transcriptomes resulted in ~40 000 and ~55 000 contigs, respectively, both with an average contig size ~945. We observed a high ratio of shared polymorphisms to fixed differences in the species pair and minimal differences between species in the ratio of synonymous to replacement substitutions or codon usage bias implying comparable effective population sizes throughout species divergence. Our results suggest that differences in effective population size and selection efficacy in mixed mating taxa shortly after their divergence may be minimal and are likely influenced by fluctuating mating systems and population sizes.

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

Secondary evolution of a self-incompatibility locus in the Brassicaceae genus Leavenworthia

Self-incompatibility (SI) is the flowering plant reproductive system in which self pollen tube growth is inhibited, thereby preventing self-fertilization. SI has evolved independently in several different flowering plant lineages. In all Brassicaceae species in which the molecular basis of SI has been investigated in detail, the product of the S-locus receptor kinase (SRK) gene functions as receptor in the initial step of the self pollen-rejection pathway, while that of the S-locus cysteine-rich (SCR) gene functions as ligand. Here we examine the hypothesis that the S locus in the Brassicaceae genus Leavenworthia is paralogous with the S locus previously characterized in other members of the family. We also test the hypothesis that self-compatibility in this group is based on disruption of the pollen ligand-producing gene. Sequence analysis of the S-locus genes in Leavenworthia, phylogeny of S alleles, gene expression patterns, and comparative genomics analyses provide support for both hypotheses. Of special interest are two genes located in a non-S locus genomic region of Arabidopsis lyrata that exhibit domain structures, sequences, and phylogenetic histories similar to those of the S-locus genes in Leavenworthia, and that also share synteny with these genes. These A. lyrata genes resemble those comprising the A. lyrata S locus, but they do not function in self-recognition. Moreover, they appear to belong to a lineage that diverged from the ancestral Brassicaceae S-locus genes before allelic diversification at the S locus. We hypothesize that there has been neo-functionalization of these S-locus-like genes in the Leavenworthia lineage, resulting in evolution of a separate ligand-receptor system of SI. Our results also provide support for theoretical models that predict that the least constrained pathway to the evolution of self-compatibility is one involving loss of pollen gene function.

Secondary evolution of a self-incompatibility locus in the Brassicaceae genus Leavenworthia

Self-incompatibility (SI) is the flowering plant reproductive system in which self pollen tube growth is inhibited, thereby preventing self-fertilization. SI has evolved independently in several different flowering plant lineages. In all Brassicaceae species in which the molecular basis of SI has been investigated in detail, the product of the S-locus receptor kinase (SRK) gene functions as receptor in the initial step of the self pollen-rejection pathway, while that of the S-locus cysteine-rich (SCR) gene functions as ligand. Here we examine the hypothesis that the S locus in the Brassicaceae genus Leavenworthia is paralogous with the S locus previously characterized in other members of the family. We also test the hypothesis that self-compatibility in this group is based on disruption of the pollen ligand-producing gene. Sequence analysis of the S-locus genes in Leavenworthia, phylogeny of S alleles, gene expression patterns, and comparative genomics analyses provide support for both hypotheses. Of special interest are two genes located in a non-S locus genomic region of Arabidopsis lyrata that exhibit domain structures, sequences, and phylogenetic histories similar to those of the S-locus genes in Leavenworthia, and that also share synteny with these genes. These A. lyrata genes resemble those comprising the A. lyrata S locus, but they do not function in self-recognition. Moreover, they appear to belong to a lineage that diverged from the ancestral Brassicaceae S-locus genes before allelic diversification at the S locus. We hypothesize that there has been neo-functionalization of these S-locus-like genes in the Leavenworthia lineage, resulting in evolution of a separate ligand-receptor system of SI. Our results also provide support for theoretical models that predict that the least constrained pathway to the evolution of self-compatibility is one involving loss of pollen gene function.

Secondary evolution of a self-incompatibility locus in the Brassicaceae genus Leavenworthia

Self-incompatibility (SI) is the flowering plant reproductive system in which self pollen tube growth is inhibited, thereby preventing self-fertilization. SI has evolved independently in several different flowering plant lineages. In all Brassicaceae species in which the molecular basis of SI has been investigated in detail, the product of the S-locus receptor kinase (SRK) gene functions as receptor in the initial step of the self pollen-rejection pathway, while that of the S-locus cysteine-rich (SCR) gene functions as ligand. Here we examine the hypothesis that the S locus in the Brassicaceae genus Leavenworthia is paralogous with the S locus previously characterized in other members of the family. We also test the hypothesis that self-compatibility in this group is based on disruption of the pollen ligand-producing gene. Sequence analysis of the S-locus genes in Leavenworthia, phylogeny of S alleles, gene expression patterns, and comparative genomics analyses provide support for both hypotheses. Of special interest are two genes located in a non-S locus genomic region of Arabidopsis lyrata that exhibit domain structures, sequences, and phylogenetic histories similar to those of the S-locus genes in Leavenworthia, and that also share synteny with these genes. These A. lyrata genes resemble those comprising the A. lyrata S locus, but they do not function in self-recognition. Moreover, they appear to belong to a lineage that diverged from the ancestral Brassicaceae S-locus genes before allelic diversification at the S locus. We hypothesize that there has been neo-functionalization of these S-locus-like genes in the Leavenworthia lineage, resulting in evolution of a separate ligand-receptor system of SI. Our results also provide support for theoretical models that predict that the least constrained pathway to the evolution of self-compatibility is one involving loss of pollen gene function.