Patterns of molecular evolution in dioecious and non-dioecious Silene

In situ radiation explains the frequency of dioecious palms on islands

BACKGROUND AND AIMS

Dioecy has evolved up to 5000 times in angiosperms, despite the potentially high intrinsic costs to unisexuality. Dioecy prevents inbreeding, which is especially relevant on isolated islands when gene pools are small. Dioecy is also associated with certain dispersal traits, such as fruit size and type. However, the influence of dioecy on other life history traits and island distribution remains poorly understood. Here, we test the effect of dioecy on palm (Arecaceae) speciation rates, fruit size and frequency on islands.

METHODS

We used phylogenetic comparative methods to estimate the ancestral state of the sexual system and its impact on speciation rates and fruit size. Frequency of sexual systems, effect of insularity on the probability of being dioecious, and phylogenetic clustering of island dioecious vs. mainland species were inferred. Lastly, we determined the interplay of insularity and sexual system on speciation rates.

KEY RESULTS

Palms repeatedly evolved different sexual systems (dioecy, monoecy and polygamy) from a hermaphrodite origin. Differences in speciation rates and fruit size among the different sexual systems were not identified. An effect of islands on the probability of the palms being dioecious was also not found. However, we found a high frequency and phylogenetic clustering of dioecious palms on islands, which were not correlated with higher speciation rates.

CONCLUSIONS

The high frequency and phylogenetic clustering may be the result of in situ radiation and suggest an 'island effect' for dioecious palms, which was not explained by differential speciation rates. This island effect also cannot be attributed to long-distance dispersal due to the lack of fruit size difference among sexual systems, and particularly because palm dispersal to islands is highly constrained by the interaction between the sizes of fruit and frugivores. Taken together, we suggest that trait flexibility in sexual system evolution and the in situ radiation of dioecious lineages are the underlying causes of the outstanding distribution of palms on islands.

Dioecy Is Associated with High Genetic Diversity and Adaptation Rates in the Plant Genus Silene

About 15,000 angiosperm species (∼6%) have separate sexes, a phenomenon known as dioecy. Why dioecious taxa are so rare is still an open question. Early work reported lower species richness in dioecious compared with nondioecious sister clades, raising the hypothesis that dioecy may be an evolutionary dead-end. This hypothesis has been recently challenged by macroevolutionary analyses that detected no or even positive effect of dioecy on diversification. However, the possible genetic consequences of dioecy at the population level, which could drive the long-term fate of dioecious lineages, have not been tested so far. Here, we used a population genomics approach in the Silene genus to look for possible effects of dioecy, especially for potential evidence of evolutionary handicaps of dioecy underlying the dead-end hypothesis. We collected individual-based RNA-seq data from several populations in 13 closely related species with different sexual systems: seven dioecious, three hermaphroditic, and three gynodioecious species. We show that dioecy is associated with increased genetic diversity, as well as higher selection efficacy both against deleterious mutations and for beneficial mutations. The results hold after controlling for phylogenetic inertia, differences in species census population sizes and geographic ranges. We conclude that dioecious Silene species neither show signs of increased mutational load nor genetic evidence for extinction risk. We discuss these observations in the light of the possible demographic differences between dioecious and self-compatible hermaphroditic species and how this could be related to alternatives to the dead-end hypothesis to explain the rarity of dioecy.

Modeling a trait-dependent diversification process coupled with molecular evolution on a random species tree

Understanding the evolution of binary traits, which affects the birth and survival of species and also the rate of molecular evolution, remains challenging. In this work, we present a probabilistic modeling framework for binary trait, random species trees, in which the number of species and their traits are represented by an asymmetric, two-type, continuous time Markov branching process. The model involves a number of different parameters describing both character and molecular evolution on the so-called 'reduced' tree, consisting of only extant species at the time of observation. We expand our model by considering the impact of binary traits on dN/dS, the normalized ratio of nonsynonymous to synonymous substitutions. We also develop mechanisms which enable us to understand the substitution rates on a phylogenetic tree with regards to the observed traits. The properties obtained from the model are illustrated with a phylogeny of outcrossing and selfing plant species, which allows us to investigate not only the branching tree rates, but also the molecular rates and the intensity of selection.

SEX-DETector: A Probabilistic Approach to Study Sex Chromosomes in Non-Model Organisms

We propose a probabilistic framework to infer autosomal and sex-linked genes from RNA-seq data of a cross for any sex chromosome type (XY, ZW, and UV). Sex chromosomes (especially the non-recombining and repeat-dense Y, W, U, and V) are notoriously difficult to sequence. Strategies have been developed to obtain partially assembled sex chromosome sequences. Most of them remain difficult to apply to numerous non-model organisms, either because they require a reference genome, or because they are designed for evolutionarily old systems. Sequencing a cross (parents and progeny) by RNA-seq to study the segregation of alleles and infer sex-linked genes is a cost-efficient strategy, which also provides expression level estimates. However, the lack of a proper statistical framework has limited a broader application of this approach. Tests on empirical Silene data show that our method identifies 20-35% more sex-linked genes than existing pipelines, while making reliable inferences for downstream analyses. Approximately 12 individuals are needed for optimal results based on simulations. For species with an unknown sex-determination system, the method can assess the presence and type (XY vs. ZW) of sex chromosomes through a model comparison strategy. The method is particularly well optimized for sex chromosomes of young or intermediate age, which are expected in thousands of yet unstudied lineages. Any organisms, including non-model ones for which nothing is known a priori, that can be bred in the lab, are suitable for our method. SEX-DETector and its implementation in a Galaxy workflow are made freely available.

Diversity of sexual systems within different lineages of the genus Silene

Species and populations can be categorized by their sexual systems, depending on the spatial distribution of female and male reproductive structures within and among plants. Although a high diversity of sexual systems exists in Silene, their relative frequency at the genus and infrageneric level is unknown. Here, we carried out an extensive literature search for direct or indirect descriptions of sexual systems in Silene species. We found descriptions of sexual systems for 98 Silene species, where 63 and 35 correspond to the phylogenetically supported subgenera Silene and Behenantha, respectively. Hermaphroditism was the commonest sexual system (58.2 %), followed by dioecy (14.3 %), gynodioecy (13.3 %) and gynodioecy-gynomonoecy (i.e. hermaphroditic, female and gynomonoecious plants coexisting in the same population; 12.2 %). The presence of these sexual systems in both subgenera suggests their multiple origins. In 17 species, the description of sexual systems varied, and in most cases these differences corresponded to variations within or among populations. Interestingly, the poorly studied gynodioecy-gynomonoecy sexual system showed similar frequency to dioecy and gynodioecy in both subgenera. In addition, the incidence of gynodioecy-gynomonoecy was analysed in the species of section Psammophilae (Silene littorea, S. psammitis, S. adscendens and S. cambessedesii), in a survey of 26 populations across the distribution area of the species. The four species showed gynomonoecy-gynodioecy in most populations. Hermaphrodites were the most frequent morph, with a low number of females and gynomonoecious plants in all populations. The frequency of sexual morphs varied significantly among the studied populations but not among species. Female plants generally produced smaller numbers of flowers than hermaphroditic or gynomonoecious plants, and the percentages of female flowers per population were low. All these findings suggest that the gynodioecious-gynomonoecious sexual system in section Psammophilae is closer to hermaphroditism or gynomonoecy than gynodioecy.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

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

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

Standard sister clade comparison fails when testing derived character States

Comparing species richness in sister clades that differ in a character state is one of the ways to study factors influencing diversification. While most of its applications have focussed on traits that increase diversification, some have been used to study the association of a trait with lower species richness, e.g., the occurrence of dioecy in flowering plants. We show here, using simulations and an analytical model, that the null expectation of equal species richness that is generally used in sister clade comparisons is wrong in the case of a derived trait occurring independently from speciation: one should expect fewer species in the clade with the derived character state when there is no difference in diversification rates. This is due to the waiting time for the derived state to appear, which causes it to occur more often on longer branches. This has the important implication that the probability for a clade to possess the derived state depends on the tree geometry, and thus on species richness: species-poorer clades are more likely to possess the derived state. We develop a statistical test for sister clade comparisons to study the effect of a derived character state. Applying it to a data set of dioecious clades, we find that we cannot confirm earlier work that concluded that dioecy decreases diversification; on the contrary, it seems to be associated to higher species richness than expected. [angiosperms; dioecy; diversification; sister clades; species richness.].

Evolution of sex determination systems with heterogametic males and females in silene

The plant genus Silene has become a model for evolutionary studies of sex chromosomes and sex-determining mechanisms. A recent study performed in Silene colpophylla showed that dioecy and the sex chromosomes in this species evolved independently from those in Silene latifolia, the most widely studied dioecious Silene species. The results of this study show that the sex-determining system in Silene otites, a species related to S. colpophylla, is based on female heterogamety, a sex determination system that is unique among the Silene species studied to date. Our phylogenetic data support the placing of S. otites and S. colpophylla in the subsection Otites and the analysis of ancestral states suggests that the most recent common ancestor of S. otites and S. colpophylla was most probably dioecious. These observations imply that a switch from XX/XY sex determination to a ZZ/ZW system (or vice versa) occurred in the subsection Otites. This is the first report of two different types of heterogamety within one plant genus of this mostly nondioecious plant family.

Disentangling the effects of mating systems and mutation rates on cytoplasmic [correction of cytoplamic] diversity in gynodioecious Silene nutans and dioecious Silene otites

Many flowering plant species exhibit a variety of distinct sexual morphs, the two most common cases being the co-occurrence of females and males (dioecy) or the co-occurrence of hermaphrodites and females (gynodioecy). In this study, we compared DNA sequence variability of the three genomes (nuclear, mitochondrial and chloroplastic) of a gynodioecious species, Silene nutans, with that of a closely related dioecious species, Silene otites. In the light of theoretical models, we expect cytoplasmic diversity to differ between the two species due to the selective dynamics that acts on cytoplasmic genomes in gynodioecious species: under an epidemic scenario, the gynodioecious species is expected to exhibit lower cytoplasmic diversity than the dioecious species, while the opposite is expected in the case of balancing selection maintaining sterility cytoplasms in the gynodioecious species. We found no difference between the species for nuclear gene diversity, but, for the cytoplasmic loci, the gynodioecious S. nutans had more haplotypes, and higher nucleotide diversity, than the dioecious relative, S. otites, even though the latter has a relatively high rate of mitochondrial synonymous substitutions, and therefore presumably a higher mutation rate. Therefore, as the mitochondrial mutation rate cannot account for the higher cytoplasmic diversity found in S. nutans, our findings support the hypothesis that gynodioecy in S. nutans has been maintained by balancing selection rather than by epidemic-like dynamics.