The population genetics of sporophytic self-incompatibility in Senecio squalidus L. (Asteraceae): avoidance of mating constraints imposed by low S-allele number

Genomic changes and stabilization following homoploid hybrid speciation of the Oxford ragwort Senecio squalidus

Oxford ragwort (Senecio squalidus) is one of only two homoploid hybrid species known to have originated very recently, so it is a unique model for determining genomic changes and stabilization following homoploid hybrid speciation. Here, we provide a chromosome-level genome assembly of S. squalidus with 95% of the assembly contained in the 10 longest scaffolds, corresponding to its haploid chromosome number. We annotated 30,249 protein-coding genes and estimated that ∼62% of the genome consists of repetitive elements. We then characterized genome-wide patterns of linkage disequilibrium, polymorphism, and divergence in S. squalidus and its two parental species, finding that (1) linkage disequilibrium is highly heterogeneous, with a region on chromosome 4 showing increased values across all three species but especially in S. squalidus; (2) regions harboring genetic incompatibilities between the two parental species tend to be large, show reduced recombination, and have lower polymorphism in S. squalidus; (3) the two parental species have an unequal contribution (70:30) to the genome of S. squalidus, with long blocks of parent-specific ancestry supporting a very rapid stabilization of the hybrid lineage after hybrid formation; and (4) genomic regions with major parent ancestry exhibit an overrepresentation of loci with evidence for divergent selection occurring between the two parental species on Mount Etna. Our results show that both genetic incompatibilities and natural selection play a role in determining genome-wide reorganization following hybrid speciation and that patterns associated with homoploid hybrid speciation-typically seen in much older systems-can evolve very quickly following hybridization.

Genome-wide identification and characterization of SRLK gene family reveal their roles in self-incompatibility of Erigeron breviscapus

Self-incompatibility (SI) is a reproductive protection mechanism that plants acquired during evolution to prevent self-recession. As the female determinant of SI specificity, SRK has been shown to be the only recognized gene on the stigma and plays important roles in SI response. Asteraceae is the largest family of dicotyledonous plants, many of which exhibit self-incompatibility. However, systematic studies on SRK gene family in Asteraceae are still limited due to lack of high-quality genomic data. In this study, we performed the first systematic genome-wide identification of S-locus receptor like kinases (SRLKs) in the self-incompatible Asteraceae species, Erigeron breviscapus, which is also a widely used perennial medicinal plant endemic to China.52 SRLK genes were identified in the E. breviscapus genome. Structural analysis revealed that the EbSRLK proteins in E. breviscapus are conserved. SRLK proteins from E. breviscapus and other SI plants are clustered into 7 clades, and the majority of the EbSRLK proteins are distributed in Clade I. Chromosomal and duplication analyses indicate that 65% of the EbSRLK genes belong to tandem repeats and could be divided into six tandem gene clusters. Gene expression patterns obtained in E. breviscapus multiple-tissue RNA-Seq data revealed differential temporal and spatial features of EbSRLK genes. Among these, two EbSRLK genes having high expression levels in tongue flowers were cloned. Subcellular localization assay demonstrated that both of their fused proteins are localized on the plasma membrane. All these results indicated that EbSRLK genes possibly involved in SI response in E. breviscapus. This comprehensive genome-wide study of the SRLK gene family in E. breviscapus provides valuable information for understanding the mechanism of SSI in Asteraceae.

Molecular cloning and yeast two-hybrid provide new evidence for unique sporophytic self-incompatibility system of Corylus

The Corylus genus contains several important nut producing species and exhibits sporophytic self-incompatibility (SSI). However, the underlying molecular mechanisms of SSI in Corylus remain largely unknown. To clarify whether Corylus and Brassica share the same SSI molecular mechanism. We cloned ChaTHL1/2, ChaMLPK, ChaARC1, ChaEX70A1 genes from Ping'ou hybrid hazelnut using RACE techniques and tested the interaction between the ChaARC1 and ChaSRK1/2. We also examined the pistil-pollen interactions using scanning electron microscopy. We found no differences in the stigma surface within 1 h after compatible or incompatible pollination. Compatible pollen tubes penetrated the stigma surface, while incompatible pollen did not penetrate the stigma 4 h after pollination. Bioinformatics analysis revealed that ChaTHL1/2, ChaMLPK, ChaARC1 and ChaEX70A1 have corresponding functional domains. Quantitative real-time PCR (qRT-PCR) analysis showed that ChaTHL1/2, ChaMLPK, ChaARC1 and ChaEX70A1 were not regularly expressed in compatible or incompatible pollination. Furthermore, the expression patterns of ARC1, THL1/2, MLPK and Exo70A1 were quite distinct between Corylus and Brassica. According to yeast two-hybrid assays, ChaSRK1/2 did not interact with ChaARC1, confirming that the SRK-ARC1 signalling pathway implicated in the SSI response of Brassica was not conserved in Corylus. These results further reinforce the conclusion that, notwithstanding the similarity of the genetic basis, the SSI mechanism of Corylus does not conform in many respects with that of Brassica. Our findings could be helpful to better explore the potential mechanism of SSI system in Corylus.

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.‬.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

Population genetics of self-incompatibility in a clade of relict cliff-dwelling plant species

The mating systems of species in small or fragmented populations impact upon their persistence. Small self-incompatible (SI) populations risk losing S allele diversity, responsible for the SI response, by drift thereby limiting mate availability and leading to population decline or SI system breakdown. But populations of relict and/or endemic species have resisted these demographic conditions over long periods suggesting their mating systems have adapted. To address a lack of empirical data on this topic, we studied the SI systems of three relict cliff-dwelling species of Sonchus section Pustulati (Asteraceae): S. masguindalii, S. fragilis and S. pustulatus in the western Mediterranean region. We performed controlled pollinations within and between individuals to measure index of SI (ISI) expression and identify S alleles in multiple population samples. Sonchus masguindalii and S. pustulatus showed strong SI (ISI = 0.6-1.0) compared to S. fragilis (ISI = 0.1-0.7). Just five S alleles were estimated for Spanish S. pustulatus and a moderate 11-15 S alleles for Moroccan S. pustulatus and S. fragilis, respectively. The fact that autonomous fruit set was generally improved by active self-pollination in self-compatible S. fragilis suggests that individuals with weak SI can show a wide range of outcrossing levels dependent on the degree of self or outcross pollen that pollinators bear. We conclude that frequent S allele dominance interactions that mask the incompatibility interactions of recessive S alleles leading to higher mate availability and partial breakdown of SI leading to mixed mating, both contribute to reproductive resilience in this group.

Transcriptomic comparison of the self-pollinated and cross-pollinated flowers of Erigeron breviscapus to analyze candidate self-incompatibility-associated genes

BACKGROUND

Self-incompatibility (SI) is a widespread and important mating system that promotes outcrossing in plants. Erigeron breviscapus, a medicinal herb used widely in traditional Chinese medicine, is a self-incompatible species of Asteraceae. However, the genetic characteristics of SI responses in E. breviscapus remain largely unknown. To understand the possible mechanisms of E. breviscapus in response to SI, we performed a comparative transcriptomic analysis with capitulum of E. breviscapus after self- and cross-pollination, which may provide valuable information for analyzing the candidate SI-associated genes of E. breviscapus.

METHODS

Using a high-throughput next-generation sequencing (Illumina) approach, the transcriptionexpression profiling of the different genes of E. breviscapus were obtained, some results were verified by quantitative real time PCR (qRT-PCR).

RESULTS

After assembly, 63,485 gene models were obtained (average gene size 882 bp; N50 = 1485 bp), among which 38,540 unigenes (60.70% of total gene models) were annotated by comparisons with four public databases (Nr, Swiss-Prot, KEGG and COG): 38,338 unigenes (60.38% of total gene models) showed high homology with sequences in the Nr database. Differentially expressed genes were identified among the three cDNA libraries (non-, self- and cross-pollinated capitulum of E. breviscapus), and approximately 230 genes might be associated with SI responses. Several these genes were upregulated in self-pollinated capitulum but downregulated in cross-pollinated capitulum, such as SRLK (SRK-like) and its downstream signal factor, MLPK. qRT-PCR confirmed that the expression patterns of EbSRLK1 and EbSRLK3 genes were not closely related to SI of E. breviscapus.

CONCLUSIONS

This work represents the first large-scale analysis of gene expression in the self-pollinated and cross-pollinated flowers of E. breviscapus. A larger number of notable genes potentially involved in SI responses showed differential expression, including genes playing crucial roles in cell-cell communication, signal transduction and the pollination process. We thus hypothesized that those genes showing differential expression and encoding critical regulators of SI responses, such as MLPK, ARC1, CaM, Exo70A1, MAP, SF21 and Nod, might affect SI responses in E. breviscapus. Taken together, our study provides a pool of SI-related genes in E. breviscapus and offers a valuable resource for elucidating the mechanisms of SI in Asteraceae.

Incest versus abstinence: reproductive trade-offs between mate limitation and progeny fitness in a self-incompatible invasive plant

Plant mating systems represent an evolutionary and ecological trade-off between reproductive assurance through selfing and maximizing progeny fitness through outbreeding. However, many plants with sporophytic self-incompatibility systems exhibit dominance interactions at the S-locus that allow biparental inbreeding, thereby facilitating mating between individuals that share alleles at the S-locus. We investigated this trade-off by estimating mate availability and biparental inbreeding depression in wild radish from five different populations across Australia. We found dominance interactions among S-alleles increased mate availability relative to estimates based on individuals that did not share S-alleles. Twelve of the sixteen fitness variables were significantly reduced by inbreeding. For all the three life-history phases evaluated, self-fertilized offspring suffered a greater than 50% reduction in fitness, while full-sib and half-sib offspring suffered a less than 50% reduction in fitness. Theory indicates that fitness costs greater than 50% can result in an evolutionary trajectory toward a stable state of self-incompatibility (SI). This study suggests that dominance interactions at the S-locus provide a possible third stable state between SI and SC where biparental inbreeding increases mate availability with relatively minor fitness costs. This strategy allows weeds to establish in new environments while maintaining a functional SI system.

A review of the allozyme data set for the Canarian endemic flora: causes of the high genetic diversity levels and implications for conservation

Background and Aims Allozyme and reproductive data sets for the Canarian flora are updated in order to assess how the present levels and structuring of genetic variation have been influenced by the abiotic island traits and by phylogenetically determined biotic traits of the corresponding taxa; and in order to suggest conservation guidelines. Methods Kruskal-Wallis tests are conducted to assess the relationships of 27 variables with genetic diversity (estimated by A, P, Ho and He) and structuring (GST) of 123 taxa representing 309 populations and 16 families. Multiple linear regression analyses (MLRAs) are carried out to determine the relative influence of the less correlated significant abiotic and biotic factors on the genetic diversity levels. Key Results and Conclusions The interactions between biotic features of the colonizing taxa and the abiotic island features drive plant diversification in the Canarian flora. However, the lower weight of closeness to the mainland than of (respectively) high basic chromosome number, partial or total self-incompatibility and polyploidy in the MLRAs indicates substantial phylogenetic constraint; the importance of a high chromosome number is feasibly due to the generation of a larger number of linkage groups, which increase gametic and genotypic diversity. Genetic structure is also more influenced by biotic factors (long-range seed dispersal, basic chromosome number and partial or total self-incompatibility) than by distance to the mainland. Conservation-wise, genetic structure estimates (FST/GST) only reflect endangerment under intensive population sampling designs, and neutral genetic variation levels do not directly relate to threat status or to small population sizes. Habitat protection is emphasized, but the results suggest the need for urgent implementation of elementary reproductive studies in all cases, and for ex situ conservation measures for the most endangered taxa, even without prior studies. In non-endangered endemics, multidisciplinary research is needed before suggesting case-specific conservation strategies. The molecular information relevant for conservation should be conserved in a standardized format to facilitate further insight.

Sporophytic self-incompatibility in Senecio squalidus (Asteraceae): S allele dominance interactions and modifiers of cross-compatibility and selfing rates

Understanding genetic mechanisms of self-incompatibility (SI) and how they evolve is central to understanding the mating behaviour of most outbreeding angiosperms. Sporophytic SI (SSI) is controlled by a single multi-allelic locus, S, which is expressed in the diploid (sporophyte) plant to determine the SI phenotype of its haploid (gametophyte) pollen. This allows complex patterns of independent S allele dominance interactions in male (pollen) and female (pistil) reproductive tissues. Senecio squalidus is a useful model for studying the genetic regulation and evolution of SSI because of its population history as an alien invasive species in the UK. S. squalidus maintains a small number of S alleles (7-11) with a high frequency of dominance interactions. Some S. squalidus individuals also show partial selfing and/or greater levels of cross-compatibility than expected under SSI. We previously speculated that these might be adaptations to invasiveness. Here we describe a detailed characterization of the regulation of SSI in S. squalidus. Controlled crosses were used to determine the S allele dominance hierarchy of six S alleles and effects of modifiers on cross-compatibility and partial selfing. Complex dominance interactions among S alleles were found with at least three levels of dominance and tissue-specific codominance. Evidence for S gene modifiers that increase selfing and/or cross-compatibility was also found. These empirical findings are discussed in the context of theoretical predictions for maintenance of S allele dominance interactions, and the role of modifier loci in the evolution of SI.

Genetic diversity and fitness in small populations of partially asexual, self-incompatible plants

How self-incompatibility systems are maintained in plant populations is still a debated issue. Theoretical models predict that self-incompatibility systems break down according to the intensity of inbreeding depression and number of S-alleles. Other studies have explored the function of asexual reproduction in the maintenance of self-incompatibility. However, the population genetics of partially asexual, self-incompatible populations are poorly understood and previous studies have failed to consider all possible effects of asexual reproduction or could only speculate on those effects. In this study, we investigated how partial asexuality may affect genetic diversity at the S-locus and fitness in small self-incompatible populations. A genetic model including an S-locus and a viability locus was developed to perform forward simulations of the evolution of populations of various sizes. Drift combined with partial asexuality produced a decrease in the number of alleles at the S-locus. In addition, an excess of heterozygotes was present in the population, causing an increase in mutation load. This heterozygote excess was enhanced by the self-incompatibility system in small populations. In addition, in highly asexual populations, individuals produced asexually had some fitness advantages over individuals produced sexually, because sexual reproduction produces homozygotes of the deleterious allele, contrary to asexual reproduction. Our results suggest that future research on the function of asexuality for the maintenance of self-incompatibility will need to (1) account for whole-genome fitness (mutation load generated by asexuality, self-incompatibility and drift) and (2) acknowledge that the maintenance of self-incompatibility may not be independent of the maintenance of sex itself.

Effect of variation in self-incompatibility on pollen limitation and inbreeding depression in Flourensia cernua (Asteraceae) scrubs of contrasting density

BACKGROUND AND AIMS

Selection may favour a partial or complete loss of self-incompatibility (SI) if it increases the reproductive output of individuals in the presence of low mate availability. The reproductive output of individuals varying in their strength of SI may also be affected by population density via its affect on the spatial structuring and number of S-alleles in populations. Modifiers increasing levels of self-compatibility can be selected when self-compatible individuals receive reproductive compensation by, for example, increasing seed set and/or when they become associated with high fitness genotypes.

METHODS

The effect of variation in the strength of SI and scrub density (low versus high) on seed set, seed germination and inbreeding depression in seed germination (delta(germ)) was investigated in the partially self-incompatible species Flourensia cernua by analysing data from self-, cross- and open-pollinated florets.

KEY RESULTS

Examination of 100 plants in both high and low scrub densities revealed that 51% of plants were strongly self-incompatible and 49 % varied from being self-incompatible to self-compatible. Seed set after hand cross-pollination was higher than after open-pollination for self-incompatible, partially self-incompatible and self-compatible plants but was uniformly low for strongly self-incompatible plants. Strongly self-incompatible and self-incompatible plants exhibited lower seed set, seed germination and multiplicative female fitness (floral display x seed set x seed germination) in open-pollinated florets compared with partially self-incompatible and self-compatible plants. Scrub density also had an effect on seed set and inbreeding depression: in low-density scrubs seed set was higher after open-pollination and delta(germ) was lower.

CONCLUSIONS

These data suggest that (a) plants suffered outcross pollen limitation, (b) female fitness in partially self-incompatible and self-compatible plants is enhanced by increased mate-compatibility and (c) plants in low-density scrubs received higher quality pollen via open-pollination than plants in high-density scrubs.

The evolution of self-incompatibility when mates are limiting

Self-incompatibility (SI) is a genetic barrier to inbreeding that is broadly distributed in angiosperms. In finite populations of SI plants, the loss of S-allele diversity can limit plant reproduction by reducing the availability of compatible mates. Many studies have shown that small or fragmented plant populations suffer from mate limitation. The advent of molecular typing of S-alleles in many species has paved the way to address quantitatively the importance of mate limitation, and to provide greater insight into why and how SI systems breakdown frequently in nature. In this review, we highlight the ecological factors that contribute to mate limitation in SI taxa, discuss their consequences for the evolution and functioning of SI, and propose new empirical research directions.

A tale of two continents: Baker's rule and the maintenance of self-incompatibility in Lycium (Solanaceae)

Over 50 years ago, Baker (1955, 1967) suggested that self-compatible species were more likely than self-incompatible species to establish new populations on oceanic islands. His logic was straightforward and rested on the assumption that colonization was infrequent; thus, mate limitation favored the establishment of self-fertilizing individuals. In support of Baker's rule, many authors have documented high frequencies of self-compatibility on islands, and recent work has solidified the generality of Baker's ideas. The genus Lycium (Solanaceae) has ca. 80 species distributed worldwide, and phylogenetic studies suggest that Lycium originated in South America and dispersed to the Old World a single time. Previous analyses of the S-RNase gene, which controls the stylar component of self-incompatibility, have shown that gametophytically controlled self-incompatibility is ancestral within the genus, making Lycium a good model for investigating Baker's assertions concerning reproductive assurance following oceanic dispersal. Lycium is also useful for investigations of reproductive evolution, given that species vary both in sexual expression and the presence of self-incompatibility. A model for the evolution of gender dimorphism suggests that polyploidy breaks down self-incompatibility, leading to the evolution of gender dimorphism, which arises as an alternative outcrossing mechanism. There is a perfect association of dimorphic gender expression, polyploidy, and self-compatibility (vs. cosexuality, diploidy, and self-incompatibility) among North American Lycium. Although the association between ploidy level and gender expression also holds for African Lycium, to date no studies of mating systems have been initiated in Old World species. Here, using controlled pollinations, we document strong self-incompatibility in two cosexual, diploid species of African Lycium. Further, we sequence the S-RNase gene in 15 individuals from five cosexual, diploid species of African Lycium and recover 24 putative alleles. Genealogical analyses indicate reduced trans-generic diversity of S-RNases in the Old World compared to the New World. We suggest that genetic diversity at this locus was reduced as a result of a founder event, but, despite the bottleneck, self-incompatibility was maintained in the Old World. Maximum-likelihood analyses of codon substitution patterns indicate that positive Darwinian selection has been relatively strong in the Old World, suggesting the rediversification of S-RNases following a bottleneck. The present data thus provide a dramatic exception to Baker's rule, in addition to supporting a key assumption of the Miller and Venable (2000) model, namely that self-incompatibility is associated with diploidy and cosexuality.

Dominance relationships between self-incompatibility alleles controlled by DNA methylation

In crucifers, the pollen S-determinant gene, SP11, is sporophytically expressed in the anther tapetum, and the pollen self-incompatibility phenotype is determined by the dominance relationships between the two S-haplotypes it carries. We report here that 5' promoter sequences of recessive SP11 alleles are specifically methylated in the tapetum before the initiation of SP11 transcription. These results suggest that tissue-specific monoallelic de novo DNA methylation is involved in determining the dominance interactions that determine the cruciferous self-incompatibility phenotype.

Polyploidy and self-compatibility: is there an association?

•  Researchers have hypothesized that self-compatibility (SC) should be more common in polyploid taxa than their diploid counterparts because of selection for reproductive assurance and/or the expected decline in inbreeding depression associated with having 'extra' gene copies. Support for this view has come from an observed breakdown of self-incompatibility (SI) in some species with a gametophytic system (GSI). The purpose of this research was to assess the strength of this relationship across a wider array of SI systems. •  A large database, of diploid chromosome numbers, ploidy levels, and types of SI system, was assembled for angiosperm species and used to test for an association between ploidy and SC. •  No strong association was found between SC and polyploidy at the level of species or families, and there was no evidence that those having a functional SI system also had fewer polyploid taxa or that most polyploids experience a breakdown in SI. •  These results challenge the assumption that self-fertilization is strongly associated with polyploidy and suggest directions for further research on the evolution of polyploidy in relation to SI.

Pollen recognition and rejection during the sporophytic self-incompatibility response: Brassica and beyond

Many hermaphrodite flowering plants avoid self-fertilization through genetic systems of self-incompatibility (SI). SI allows a plant to recognize and to reject self or self-related pollen, thereby preserving its ovules for outcrossing. Genes situated at the S-locus encode the 'male' (pollen) and 'female' (pistil) recognition determinants of SI. In sporophytic SI (SSI) the male determinant is expressed in the diploid anther, therefore haploid pollen grains behave with a diploid S phenotype. In Brassica, the male and the female determinants of SSI have been identified as a peptide ligand and its cognate receptor, respectively, and recent studies have identified downstream signalling molecules involved in pollen rejection. It now needs to be established whether the Brassica mechanism is universal in species with SSI, or unique to the Brassicaceae.

The different mechanisms of sporophytic self-incompatibility

Flowering plants have evolved a multitude of mechanisms to avoid self-fertilization and promote outbreeding. Self-incompatibility (SI) is by far the most common of these, and is found in ca. 60% of flowering plants. SI is a genetically controlled pollen-pistil recognition system that provides a barrier to fertilization by self and self-related pollen in hermaphrodite (usually co-sexual) flowering plants. Two genetically distinct forms of SI can be recognized: gametophytic SI (GSI) and sporophytic SI (SSI), distinguished by how the incompatibility phenotype of the pollen is determined. GSI appears to be the most common mode of SI and can operate through at least three different mechanisms, two of which have been characterized extensively at a molecular level in the Solanaceae and Papaveraceae. Because molecular studies of SSI have been largely confined to species from the Brassicaceae, predominantly Brassica species, it is not yet known whether SSI, like GSI, can operate through different molecular mechanisms. Molecular studies of SSI are now being carried out on Ipomoea trifida (Convolvulaceae) and Senecio squalidus (Asteraceae) and are providing important preliminary data suggesting that SSI in these two families does not share the same molecular mechanism as that of the Brassicaceae. Here, what is currently known about the molecular regulation of SSI in the Brassicaceae is briefly reviewed, and the emerging data on SSI in I. trifida, and more especially in S. squalidus, are discussed.