Dynamics of drift, gene flow, and selection during speciation in Silene

The role of recombination landscape in species hybridisation and speciation

It is now well recognised that closely related species can hybridize and exchange genetic material, which may promote or oppose adaptation and speciation. In some cases, interspecific hybridisation is very common, making it surprising that species identity is preserved despite active gene exchange. The genomes of most eukaryotic species are highly heterogeneous with regard to gene density, abundance of repetitive DNA, chromatin compactisation etc, which can make certain genomic regions more prone or more resistant to introgression of genetic material from other species. Heterogeneity in local recombination rate underpins many of the observed patterns across the genome (e.g. actively recombining regions are typically gene rich and depleted for repetitive DNA) and it can strongly affect the permeability of genomic regions to interspecific introgression. The larger the region lacking recombination, the higher the chance for the presence of species incompatibility gene(s) in that region, making the entire non- or rarely recombining block impermeable to interspecific introgression. Large plant genomes tend to have highly heterogeneous recombination landscape, with recombination frequently occurring at the ends of the chromosomes and central regions lacking recombination. In this paper we review the relationship between recombination and introgression in plants and argue that large rarely recombining regions likely play a major role in preserving species identity in actively hybridising plant species.

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F_ST) and sequence divergence (d_XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for speciation in the face of gene flow.

Importance of Pollinator-Mediated Interspecific Pollen Transfer for Angiosperm Evolution

Understanding how pollen moves between species is critical to understanding speciation, diversification, and evolution of flowering plants. For co-flowering species that share pollinators, competition through interspecific pollen transfer (IPT) can profoundly impact floral evolution, decreasing female fitness via heterospecific pollen deposition on stigmas and male fitness via pollen misplacement during visits to heterospecific flowers. The pollination literature demonstrates that such reproductive interference frequently selects for reproductive character displacement in floral traits linked to pollinator attraction, pollen placement, and mating systems and has also revealed that IPT between given pairs of species is typically asymmetric. More recent work is starting to elucidate its importance to the speciation process, clarifying the link between IPT and current and historical patterns of hybridization, the evolution of phenotypic novelty through adaptive introgression, and the rise of reproductive isolation. Our review aims to stimulate further research on IPT as a ubiquitous mechanism that plays a central role in angiosperm diversification.

A coalescent-based estimator of genetic drift, and acoustic divergence in the Pteronotus parnellii species complex

Determining the processes responsible for phenotypic variation is one of the central tasks of evolutionary biology. While the importance of acoustic traits for foraging and communication in echolocating mammals suggests adaptation, the seldom-tested null hypothesis to explain trait divergence is genetic drift. Here we derive FST values from multi-locus coalescent isolation-with-migration models, and couple them with estimates of quantitative trait divergence, or PST, to test drift as the evolutionary process responsible for phenotypic divergence in island populations of the Pteronotus parnellii species complex. Compared to traditional comparisons of PST to FST, the migration-based estimates of FST are unidirectional instead of bidirectional, simultaneously integrate variation among loci and individuals, and posterior densities of PST and FST can be compared directly. We found the evolution of higher call frequencies is inconsistent with genetic drift for the Hispaniolan population, despite many generations of isolation from its Puerto Rican counterpart. While the Hispaniolan population displays dimorphism in call frequencies, the higher frequency of the females is incompatible with sexual selection. Instead, cultural drift toward higher frequencies among Hispaniolan females might explain the divergence. By integrating Bayesian coalescent and trait analyses, this study demonstrates a powerful approach to testing genetic drift as the default evolutionary mechanism of trait differentiation between populations.

Ecological divergence plays an important role in strong but complex reproductive isolation in campions (Silene)

New species arise through the evolution of reproductive barriers between formerly interbreeding lineages. Yet, comprehensive assessments of potential reproductive barriers, which are needed to make inferences on processes driving speciation, are only available for a limited number of systems. In this study, we estimated individual and cumulative strengths of seven prezygotic and six postzygotic reproductive barriers between the recently diverged taxa Silene dioica (L.) Clairv. and S. latifolia Poiret using both published and new data. A combination of multiple partial reproductive barriers resulted in near-complete reproductive isolation between S. dioica and S. latifolia, consistent with earlier estimates of gene flow between the taxa. Extrinsic barriers associated with adaptive ecological divergence were most important, while intrinsic postzygotic barriers had moderate individual strength but contributed only little to total reproductive isolation. These findings are in line with ecological divergence as driver of speciation. We further found extensive variation in extrinsic reproductive isolation, ranging from sites with very strong selection against migrants and hybrids to intermediate sites where substantial hybridization is possible. This situation may allow for, or even promote, heterogeneous genetic divergence.

Genetic architecture of traits associated with reproductive barriers in Silene: Coupling, sex chromosomes and variation

The evolution of reproductive barriers and their underlying genetic architecture is of central importance for the formation of new species. Reproductive barriers can be controlled either by few large-effect loci suggesting strong selection on key traits, or by many small-effect loci, consistent with gradual divergence or with selection on polygenic or multiple traits. Genetic coupling between reproductive barrier loci further promotes divergence, particularly divergence with ongoing gene flow. In this study, we investigated the genetic architectures of ten morphological, phenological and life history traits associated with reproductive barriers between the hybridizing sister species Silene dioica and S. latifolia; both are dioecious with XY-sex determination. We used quantitative trait locus (QTL) mapping in two reciprocal F₂ crosses. One to six QTLs per trait, including nine major QTLs (PVE > 20%), were detected on 11 of the 12 linkage groups. We found strong evidence for coupling of QTLs for uncorrelated traits and for an important role of sex chromosomes in the genetic architectures of reproductive barrier traits. Unexpectedly, QTLs detected in the two F₂ crosses differed largely, despite limited phenotypic differences between them and sufficient statistical power. The widely dispersed genetic architectures of traits associated with reproductive barriers suggest gradual divergence or multifarious selection. Coupling of the underlying QTLs likely promoted divergence with gene flow in this system. The low congruence of QTLs between the two crosses further points to variable and possibly redundant genetic architectures of traits associated with reproductive barriers, with important implications for the evolutionary dynamics of divergence and speciation.

Isolation and gene flow in a speciation continuum in newts

Because reproductive isolation often evolves gradually, differentiating lineages may retain the potential for genetic exchange for prolonged periods, providing an opportunity to quantify and to understand the fundamental role of gene flow during speciation. Here we delimit evolutionary lineages, reconstruct the phylogeny and infer gene flow in newts of the Lissotriton vulgaris species complex based on 74 nuclear markers sampled from 127 localities. We demonstrate that distinct lineages along the speciation continuum in newts exchange nontrivial amounts of genes, affecting their evolutionary trajectories. By integrating a wide array of methods, we delimit nine evolutionary lineages and show that two principal factors have driven their genetic differentiation: time since the last common ancestor determining levels of shared ancestral polymorphism, and shifts in geographic distributions determining the extent of secondary contact. Post-divergence gene flow, indicative of evolutionary non-independence, has been most extensive in Central Europe, while four southern European lineages have acquired the population-genetic hallmarks of independent species (L. graecus, L. kosswigi, L. lantzi, L. schmidtleri). We obtained strong statistical support for widespread mtDNA introgression following secondary contact, previously suggested by discordance between mtDNA phylogeny and morphology. Our study reveals long-term evolutionary persistence of evolutionary lineages that may periodically exchange genes with one another: although some of these lineages may become extinct or fuse, others will acquire complete reproductive isolation and will carry signatures of this complex history in their genomes.

Differential adaptation drives ecological speciation in campions (Silene): evidence from a multi-site transplant experiment

In order to investigate the role of differential adaptation for the evolution of reproductive barriers, we conducted a multi-site transplant experiment with the dioecious sister species Silene dioica and S. latifolia and their hybrids. Crosses within species as well as reciprocal first-generation (F₁ ) and second-generation (F₂ ) interspecific hybrids were transplanted into six sites, three within each species' habitat. Survival and flowering were recorded over 4 yr. At all transplant sites, the local species outperformed the foreign species, reciprocal F₁ hybrids performed intermediately and F₂ hybrids underperformed in comparison to F₁ hybrids (hybrid breakdown). Females generally had slightly higher cumulative fitness than males in both within- and between-species crosses and we thus found little evidence for Haldane's rule acting on field performance. The strength of selection against F₁ and F₂ hybrids as well as hybrid breakdown increased with increasing strength of habitat adaptation (i.e. the relative fitness difference between the local and the foreign species) across sites. Our results suggest that differential habitat adaptation led to ecologically dependent post-zygotic reproductive barriers and drives divergence and speciation in this Silene system.

Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex-linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within-species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex-determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation-based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex-linked genes, which requires closer linkage with the fully sex-linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.

Recombination changes at the boundaries of fully and partially sex-linked regions between closely related Silene species pairs

The establishment of a region of suppressed recombination is a critical change during sex chromosome evolution, leading to such properties as Y (and W) chromosome genetic degeneration, accumulation of repetitive sequences and heteromorphism. Although chromosome inversions can cause large regions to have suppressed recombination, and inversions are sometimes involved in sex chromosome evolution, gradual expansion of the non-recombining region could potentially sometimes occur. We here test whether closer linkage has recently evolved between the sex-determining region and several genes that are partially sex-linked in Silene latifolia, using Silene dioica, a closely related dioecious plants whose XY sex chromosome system is inherited from a common ancestor. The S. latifolia pseudoautosomal region (PAR) includes several genes extremely closely linked to the fully Y-linked region. These genes were added to an ancestral PAR of the sex chromosome pair in two distinct events probably involving translocations of autosomal genome regions causing multiple genes to become partially sex-linked. Close linkage with the PAR boundary must have evolved since these additions, because some genes added in both events now show almost complete sex linkage in S. latifolia. We compared diversity patterns of five such S. latifolia PAR boundary genes with their orthologues in S. dioica, including all three regions of the PAR (one gene that was in the ancestral PAR and two from each of the added regions). The results suggest recent recombination suppression in S. latifolia, since its split from S. dioica.

Patterns of gene flow and selection across multiple species of Acrocephalus warblers: footprints of parallel selection on the Z chromosome

Background

Understanding the mechanisms and selective forces leading to adaptive radiations and origin of biodiversity is a major goal of evolutionary biology. Acrocephalus warblers are small passerines that underwent an adaptive radiation in the last approximately 10 million years that gave rise to 37 extant species, many of which still hybridize in nature. Acrocephalus warblers have served as model organisms for a wide variety of ecological and behavioral studies, yet our knowledge of mechanisms and selective forces driving their radiation is limited. Here we studied patterns of interspecific gene flow and selection across three European Acrocephalus warblers to get a first insight into mechanisms of radiation of this avian group.

Results

We analyzed nucleotide variation at eight nuclear loci in three hybridizing Acrocephalus species with overlapping breeding ranges in Europe. Using an isolation-with-migration model for multiple populations, we found evidence for unidirectional gene flow from A. scirpaceus to A. palustris and from A. palustris to A. dumetorum. Gene flow was higher between genetically more closely related A. scirpaceus and A. palustris than between ecologically more similar A. palustris and A. dumetorum, suggesting that gradual accumulation of intrinsic barriers rather than divergent ecological selection are more efficient in restricting interspecific gene flow in Acrocephalus warblers. Although levels of genetic differentiation between different species pairs were in general not correlated, we found signatures of apparently independent instances of positive selection at the same two Z-linked loci in multiple species.

Conclusions

Our study brings the first evidence that gene flow occurred during Acrocephalus radiation and not only between sister species. Interspecific gene flow could thus be an important source of genetic variation in individual Acrocephalus species and could have accelerated adaptive evolution and speciation rate in this avian group by creating novel genetic combinations and new phenotypes. Independent instances of positive selection at the same loci in multiple species indicate an interesting possibility that the same loci might have contributed to reproductive isolation in several speciation events.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0692-2) contains supplementary material, which is available to authorized users.

Maintenance of Species Boundaries Despite Ongoing Gene Flow in Ragworts

The role of hybridization between diversifying species has been the focus of a huge amount of recent evolutionary research. While gene flow can prevent speciation or initiate species collapse, it can also generate new hybrid species. Similarly, while adaptive divergence can be wiped out by gene flow, new adaptive variation can be introduced via introgression. The relative frequency of these outcomes, and indeed the frequency of hybridization and introgression in general are largely unknown. One group of closely-related species with several documented cases of hybridization is the Mediterranean ragwort (genus: Senecio) species-complex. Examples of both polyploid and homoploid hybrid speciation are known in the clade, although their evolutionary relationships and the general frequency of introgressive hybridization among them remain unknown. Using a whole genome gene-space dataset comprising eight Senecio species we fully resolve the phylogeny of these species for the first time despite phylogenetic incongruence across the genome. Using a D-statistic approach, we demonstrate previously unknown cases of introgressive hybridization between multiple pairs of taxa across the species tree. This is an important step in establishing these species as a study system for diversification with gene flow, and suggests that introgressive hybridization may be a widespread and important process in plant evolution.

The large-X effect in plants: increased species divergence and reduced gene flow on the Silene X-chromosome

The disproportionately large involvement of the X-chromosome in the isolation of closely related species (the large-X effect) has been reported for many animals, where X-linked genes are mostly hemizygous in the heterogametic sex. The expression of deleterious recessive mutations is thought to drive the frequent involvement of the X-chromosome in hybrid sterility, as well as to reduce interspecific gene flow for X-linked genes. Here, we evaluate the role of the X-chromosome in the speciation of two closely related plant species - the white and red campions (Silene latifolia and S. dioica) - that hybridize widely across Europe. The two species evolved separate sexes and sex chromosomes relatively recently (~10(7)  years), and unlike most animal species, most X-linked genes have intact Y-linked homologs. We demonstrate that the X-linked genes show a very small and insignificant amount of interspecific gene flow, while gene flow involving autosomal loci is significant and sufficient to homogenize the gene pools of the two species. These findings are consistent with the hypothesis of the large-X effect in Silene and comprise the first report of this effect in plants. Nonhemizygosity of many X-linked genes in Silene males indicates that exposure of recessive mutations to selection may not be essential for the occurrence of the large-X effect. Several possible causes of the large-X effect in Silene are discussed.

Evidence for asymmetrical hybridization despite pre- and post-pollination reproductive barriers between two Silene species

Interspecific hybridization is widespread among plants; nevertheless, pre- and post-zygotic isolating mechanisms may maintain species integrity for interfertile species in sympatry despite some gene flow. Interspecific hybridization and potential isolating barriers were evaluated between co-flowering Silene asclepiadea and Silene yunnanensis in an alpine community in southwest China. We investigated morphological and molecular (nuclear microsatellites and chloroplast gene sequence) variation in sympatric populations of S. asclepiadea and S. yunnanensis. Additionally, we analyzed pollinator behaviour and compared reproductive success between the putative hybrids and their parental species. Both the molecular and morphological data indicate that there were putative natural hybrids in the field, with S. asclepiadae the ovule parent and S. yunnanensis the pollen parent. Bumblebees were the primary visitors to S. asclepiadae and putative hybrids, while butterflies were the primary visitors to S. yunnanensis Pollen production and viability were significantly lower in putative hybrids than the parental species. The direction of hybridization is quite asymmetric from S. yunnanensis to S. asclepiadea Protandry combined with later peak flowering of S. yunnanensis, and pollinator preference may have contributed to the asymmetric pattern of hybridization, but putative hybrids were rare. Our results thus suggest that despite gene flow, S. asclepiadea and S. yunnanensis can maintain species boundaries, perhaps as a result of floral isolation and low fecundity of the hybrids.

Differences in style length confer prezygotic isolation between two dioecious species of Silene in sympatry

One fundamental signature of reinforcement is elevated prezygotic reproductive isolation between related species in sympatry relative to allopatry. However, this alone is inadequate evidence for reinforcement, as traits conferring reproductive isolation can occur as a by-product of other forces. We conducted crosses between Silene latifolia and S. diclinis, two closely related dioecious flowering plant species. Crosses with S. latifolia mothers from sympatry exhibited lower seed set than mothers from five allopatric populations when S. diclinis was the father. However, two other allopatric populations also exhibited low seed set. A significant interaction between style length and sire species revealed that seed set declined as style length increased when interspecific, but not intraspecific, fathers where used. Moreover, by varying the distance pollen tubes had to traverse, we found interspecific pollen placement close to the ovary resulted in seed set in both long- and short-styled S. latifolia mothers. Our results reveal that the long styles of S. latifolia in sympatry with S. diclinis contribute to the prevention of hybrid formation. We argue that forces other than reinforcing selection are likely to be responsible for the differences in style length seen in sympatry.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae)

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Models of selection, isolation, and gene flow in speciation

Many marine ecologists aspire to use genetic data to understand how selection and demographic history shape the evolution of diverging populations as they become reproductively isolated species. I propose combining two types of genetic analysis focused on this key early stage of the speciation process to identify the selective agents directly responsible for population divergence. Isolation-with-migration (IM) models can be used to characterize reproductive isolation between populations (low gene flow), while codon models can be used to characterize selection for population differences at the molecular level (especially positive selection for high rates of amino acid substitution). Accessible transcriptome sequencing methods can generate the large quantities of data needed for both types of analysis. I highlight recent examples (including our work on fertilization genes in sea stars) in which this confluence of interest, models, and data has led to taxonomically broad advances in understanding marine speciation at the molecular level. I also highlight new models that incorporate both demography and selection: simulations based on these theoretical advances suggest that polymorphisms shared among individuals (a key source of information in IM models) may lead to false-positive evidence of selection (in codon models), especially during the early stages of population divergence and speciation that are most in need of study. The false-positive problem may be resolved through a combination of model improvements plus experiments that document the phenotypic and fitness effects of specific polymorphisms for which codon models and IM models indicate selection and reproductive isolation (such as genes that mediate sperm-egg compatibility at fertilization).

Maintenance of strong morphological differentiation despite ongoing natural hybridization between sympatric species of Lomatia (Proteaceae)

BACKGROUND AND AIMS

When species cohesion is maintained despite ongoing natural hybridization, many questions are raised about the evolutionary processes operating in the species complex. This study examined the extensive natural hybridization between the Australian native shrubs Lomatia myricoides and L. silaifolia (Proteaceae). These species exhibit striking differences in morphology and ecological preferences, exceeding those found in most studies of hybridization to date.

METHODS

Nuclear microsatellite markers (nSSRs), genotyping methods and morphometric analyses were used to uncover patterns of hybridization and the role of gene flow in morphological differentiation between sympatric species.

KEY RESULTS

The complexity of hybridization patterns differed markedly between sites, however, signals of introgression were present at all sites. One site provided evidence of a large hybrid swarm and the likely presence of multiple hybrid generations and backcrosses, another site a handful of early generational hybrids and a third site only traces of admixture from a past hybridization event. The presence of cryptic hybrids and a pattern of morphological bimodality amongst hybrids often disguised the extent of underlying genetic admixture.

CONCLUSIONS

Distinct parental habitats and phenotypes are expected to form barriers that contribute to the rapid reversion of hybrid populations to their parental character state, due to limited opportunities for hybrid/intermediate advantage. Furthermore, strong genomic filters may facilitate continued gene flow between species without the danger of assimilation. Stochastic fire events facilitate temporal phenological isolation between species and may partly explain the bi-directional and site-specific patterns of hybridization observed. Furthermore, the findings suggest that F1 hybrids are rare, and backcrosses may occur rapidly following these initial hybridization events.

Testing for the footprint of sexually antagonistic polymorphisms in the pseudoautosomal region of a plant sex chromosome pair

The existence of sexually antagonistic (SA) polymorphism is widely considered the most likely explanation for the evolution of suppressed recombination of sex chromosome pairs. This explanation is largely untested empirically, and no such polymorphisms have been identified, other than in fish, where no evidence directly implicates these genes in events causing loss of recombination. We tested for the presence of loci with SA polymorphism in the plant Silene latifolia, which is dioecious (with separate male and female individuals) and has a pair of highly heteromorphic sex chromosomes, with XY males. Suppressed recombination between much of the Y and X sex chromosomes evolved in several steps, and the results in Bergero et al. (2013) show that it is still ongoing in the recombining or pseudoautosomal, regions (PARs) of these chromosomes. We used molecular evolutionary approaches to test for the footprints of SA polymorphisms, based on sequence diversity levels in S. latifolia PAR genes identified by genetic mapping. Nucleotide diversity is high for at least four of six PAR genes identified, and our data suggest the existence of polymorphisms maintained by balancing selection in this genome region, since molecular evolutionary (HKA) tests exclude an elevated mutation rate, and other tests also suggest balancing selection. The presence of sexually antagonistic alleles at a locus or loci in the PAR is suggested by the very different X and Y chromosome allele frequencies for at least one PAR gene.

Patterns of molecular evolution in dioecious and non-dioecious Silene

Dioecy (i.e. having separate sexes) is a rather rare breeding system in flowering plants. Such rareness may result from a high probability of extinction in dioecious species because of less efficient dispersal and the costs of sexual selection, which are expected to harm dioecious species' survival on the long term. These handicaps should decrease the effective population size (Ne) of dioecious species, which in turn should reduce the efficacy of selection. Moreover, sexual selection in dioecious species is expected to specifically affect some genes, which will evolve under positive selection. The relative contribution of these effects is currently unknown and we tried to disentangle them by comparing sequence evolution between dioecious and non-dioecious species in the Silene genus (Caryophyllaceae), where dioecy has evolved at least twice. For the dioecious species in the section Melandrium, where dioecy is the oldest, we found a global reduction of purifying selection, while on some, male-biased genes, positive selection was found. For section Otites, where dioecy evolved more recently, we found no significant differences between dioecious and non-dioecious species. Our results are consistent with the view that dioecy is an evolutionary dead end in flowering plants, although other scenarios for explaining reduced Ne cannot be ruled out. Our results also show that contrasting forces act on the genomes of dioecious plants, and suggest that some time is required before the genome of such plants bears the footprints of dioecy.