RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region

A hemizygous supergene controls homomorphic and heteromorphic self-incompatibility systems in Oleaceae

Self-incompatibility (SI) has evolved independently multiple times and prevents self-fertilization in hermaphrodite angiosperms. Several groups of Oleaceae such as jasmines exhibit distylous flowers, with two compatibility groups each associated with a specific floral morph.1 Other Oleaceae species in the olive tribe have two compatibility groups without associated morphological variation.2,3,4,5 The genetic basis of both homomorphic and dimorphic SI systems in Oleaceae is unknown. By comparing genomic sequences of three olive subspecies (Olea europaea) belonging to the two compatibility groups, we first locate the genetic determinants of SI within a 700-kb hemizygous region present only in one compatibility group. We then demonstrate that the homologous hemizygous region also controls distyly in jasmine. Phylogenetic analyses support a common origin of both systems, following a segmental genomic duplication in a common ancestor. Examination of the gene content of the hemizygous region in different jasmine and olive species suggests that the mechanisms determining compatibility groups and floral phenotypes (whether homomorphic or dimorphic) in Oleaceae rely on the presence/absence of two genes involved in gibberellin and brassinosteroid regulation.

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

The evolutionary maintenance of ancient recombining sex chromosomes in the ostrich

Sex chromosomes have evolved repeatedly across the tree of life and often exhibit extreme size dimorphism due to genetic degeneration of the sex-limited chromosome (e.g. the W chromosome of some birds and Y chromosome of mammals). However, in some lineages, ancient sex-limited chromosomes have escaped degeneration. Here, we study the evolutionary maintenance of sex chromosomes in the ostrich (Struthio camelus), where the W remains 65% the size of the Z chromosome, despite being more than 100 million years old. Using genome-wide resequencing data, we show that the population scaled recombination rate of the pseudoautosomal region (PAR) is higher than similar sized autosomes and is correlated with pedigree-based recombination rate in the heterogametic females, but not homogametic males. Genetic variation within the sex-linked region (SLR) (π = 0.001) was significantly lower than in the PAR, consistent with recombination cessation. Conversely, genetic variation across the PAR (π = 0.0016) was similar to that of autosomes and dependent on local recombination rates, GC content and to a lesser extent, gene density. In particular, the region close to the SLR was as genetically diverse as autosomes, likely due to high recombination rates around the PAR boundary restricting genetic linkage with the SLR to only ~50Kb. The potential for alleles with antagonistic fitness effects in males and females to drive chromosome degeneration is therefore limited. While some regions of the PAR had divergent male-female allele frequencies, suggestive of sexually antagonistic alleles, coalescent simulations showed this was broadly consistent with neutral genetic processes. Our results indicate that the degeneration of the large and ancient sex chromosomes of the ostrich may have been slowed by high recombination in the female PAR, reducing the scope for the accumulation of sexually antagonistic variation to generate selection for recombination cessation.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

A CLAVATA3-like Gene Acts as a Gynoecium Suppression Function in White Campion

How do separate sexes originate and evolve? Plants provide many opportunities to address this question as they have diverse mating systems and separate sexes (dioecy) that evolved many times independently. The classic "two-factor" model for evolution of separate sexes proposes that males and females can evolve from hermaphrodites via the spread of male and female sterility mutations that turn hermaphrodites into females and males, respectively. This widely accepted model was inspired by early genetic work in dioecious white campion (Silene latifolia) that revealed the presence of two sex-determining factors on the Y-chromosome, though the actual genes remained unknown. Here, we report identification and functional analysis of the putative sex-determining gene in S. latifolia, corresponding to the gynoecium suppression factor (GSF). We demonstrate that GSF likely corresponds to a Y-linked CLV3-like gene that is specifically expressed in early male flower buds and encodes the protein that suppresses gynoecium development in S. latifolia. Interestingly, GSFY has a dysfunctional X-linked homolog (GSFX) and their synonymous divergence (dS = 17.9%) is consistent with the age of sex chromosomes in this species. We propose that female development in S. latifolia is controlled via the WUSCHEL-CLAVATA feedback loop, with the X-linked WUSCHEL-like and Y-linked CLV3-like genes, respectively. Evolution of dioecy in the S. latifolia ancestor likely involved inclusion of ancestral GSFY into the nonrecombining region on the nascent Y-chromosome and GSFX loss of function, which resulted in disbalance of the WUSCHEL-CLAVATA feedback loop between the sexes and ensured gynoecium suppression in males.

Searching for signatures of sexually antagonistic selection on stickleback sex chromosomes

Intralocus sexually antagonistic selection occurs when an allele is beneficial to one sex but detrimental to the other. This form of selection is thought to be key to the evolution of sex chromosomes but is hard to detect. Here we perform an analysis of phased young sex chromosomes to look for signals of sexually antagonistic selection in the Japan Sea stickleback (Gasterosteus nipponicus). Phasing allows us to date the suppression of recombination on the sex chromosome and provides unprecedented resolution to identify sexually antagonistic selection in the recombining region of the chromosome. We identify four windows with elevated divergence between the X and Y in the recombining region, all in or very near genes associated with phenotypes potentially under sexually antagonistic selection in humans. We are unable, however, to rule out the alternative hypothesis that the peaks of divergence result from demographic effects. Thus, although sexually antagonistic selection is a key hypothesis for the formation of supergenes on sex chromosomes, it remains challenging to detect.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

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.

Molecular Sex Identification in the Hardy Rubber Tree (Eucommia ulmoides Oliver) via ddRAD Markers

Eucommia ulmoides, also known as the industrially and medicinally important hardy rubber tree, is the sole species of Eucommiaceae. Nevertheless, its dioecious property hinders sex recognition by traditional morphological observation at very early developmental stages, thus inhibiting breeding and economic cropping. In this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) was applied to screen sex-linked molecular markers for sex identification and investigation of the sex determination system in 20 male and female E. ulmoides individual plants, respectively. In consequence, five candidate male-specific loci but no female-specific loci were predicated among the 183,752 male and 147,122 female catalogue loci by bioinformatics analysis. Subsequent PCR (polymerase chain reaction) amplification and Sanger sequencing examinations were performed on another 24 individuals, 12 for each sex, from a separate population. One ideal sex-linked locus, MSL4, was identified among the five putative male-specific loci that were found using ddRAD data. MSL4 is 479 bp in length and highly conserved in all the male individuals, suggesting its feature of being stable and repeatable. Our results also indicated that the sex of E. ulmoides is likely determined genetically. In short, this study provides a consistent and reproducible ddRAD marker (MSL4) that is able to discriminate male from female seedlings in E. ulmoides, which will be valuable for rapid breeding practice and better commercial production of this economically important tree.

The Location of the Pseudoautosomal Boundary in Silene latifolia

Y-chromosomes contain a non-recombining region (NRY), and in many organisms it was shown that the NRY expanded over time. How and why the NRY expands remains unclear. Young sex chromosomes, where NRY expansion occurred recently or is on-going, offer an opportunity to study the causes of this process. Here, we used the plant Silene latifolia, where sex chromosomes evolved ~11 million years ago, to study the location of the boundary between the NRY and the recombining pseudoautosomal region (PAR). The previous work devoted to the NRY/PAR boundary in S. latifolia was based on a handful of genes with locations approximately known from the genetic map. Here, we report the analysis of 86 pseudoautosomal and sex-linked genes adjacent to the S. latifolia NRY/PAR boundary to establish the location of the boundary more precisely. We take advantage of the dense genetic map and polymorphism data from wild populations to identify 20 partially sex-linked genes located in the “fuzzy boundary”, that rarely recombines in male meiosis. Genes proximal to this fuzzy boundary show no evidence of recombination in males, while the genes distal to this partially-sex-linked region are actively recombining in males. Our results provide a more accurate location for the PAR boundary in S. latifolia, which will help to elucidate the causes of PAR boundary shifts leading to NRY expansion over time.

Evolutionary potential for genomic islands of sexual divergence on recombining sex chromosomes

Differentiated sex chromosomes are thought to develop through the accumulation of polymorphisms at loci subject to opposing selection between males and females, and/or between haploids and diploids. As sex chromosomes differentiate, reduced recombination becomes favored between selected loci and the sex-determining region, strengthening genetic associations between alleles favored in a sex and the corresponding sex chromosome. Here a model is analyzed to explore whether polymorphism at one sexually or ploidally antagonistic locus facilitates the spread of rare alleles at other loci experiencing antagonistic selection, promoting further differentiation of the sex chromosomes. It is found that antagonistic polymorphisms can spread and capture other such loci, building 'genomic islands' of differentiation on sex chromosomes, but the conditions are very restrictive, requiring the loci to be strongly selected, tightly linked and distant from the sex-determining region. Epistatic interactions can facilitate the promotion of polymorphism among selected loci, but only if preferentially favoring heterozygotes. Although these results apply to any taxa, plants provide a fertile ground for testing these and related theories given the recurrent evolutionary transitions to dioecy, which provide multiple opportunities to track the early evolution of sex chromosomes.

Evolution of Young Sex Chromosomes in Two Dioecious Sister Plant Species with Distinct Sex Determination Systems

In the last decade, progress has been made in methods to identify the sex determination system in plants. This gives the opportunity to study sex chromosomes that arose independently at different phylogenetic scales, and thus allows the discovery and the understanding of early stages of sex chromosome evolution. In the genus Silene, sex chromosomes have evolved independently in at least two clades from a nondioecious ancestor, the Melandrium and Otites sections. In the latter, sex chromosomes could be younger than in the section Melandrium, based on phylogenetic studies and as no heteromorphic sex chromosomes have been detected. This section might also exhibit lability in sex determination, because male heterogamy and female heterogamy have been suggested to occur.

In this study, we investigated the sex determination system of two dioecious species in the section Otites (Silene otites and its close relative Silene pseudotites). Applying the new probabilistic method SEX-DETector on RNA-seq data from cross-controlled progenies, we inferred their most likely sex determination system and a list of putative autosomal and sex-linked contigs. We showed that the two phylogenetically close species differed in their sex determination system (XY versus ZW) with sex chromosomes that derived from two different pairs of autosomes. We built a genetic map of the sex chromosomes and showed that both pairs exhibited a large region with lack of recombination. However, the sex-limited chromosomes exhibited no strong degeneration. Finally, using the “ancestral” autosomal expression of sex-linked orthologs of nondioecious S. nutans, we found a slight signature of dosage compensation in the heterogametic females of S. otites.

Evolution of recombination rates between sex chromosomes

In species with genetic sex-determination, the chromosomes carrying the sex-determining genes have often evolved non-recombining regions and subsequently evolved the full set of characteristics denoted by the term 'sex chromosomes'. These include size differences, creating chromosomal heteromorphism, and loss of gene functions from one member of the chromosome pair. Such characteristics and changes have been widely reviewed, and underlie molecular genetic approaches that can detect sex chromosome regions. This review deals mainly with the evolution of new non-recombining regions, focusing on how certain evolutionary situations select for suppressed recombination (rather than the proximate mechanisms causing suppressed recombination between sex chromosomes). Particularly important is the likely involvement of sexually antagonistic polymorphisms in genome regions closely linked to sex-determining loci. These may be responsible for the evolutionary strata of sex chromosomes that have repeatedly formed by recombination suppression evolving across large genome regions. More studies of recently evolved non-recombining sex-determining regions should help to test this hypothesis empirically, and may provide evidence about whether other situations can sometimes lead to sex-linked regions evolving. Similarities with other non-recombining genome regions are discussed briefly, to illustrate common features of the different cases, though no general properties apply to all of them.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Size and Content of the Sex-Determining Region of the Y Chromosome in Dioecious Mercurialis annua, a Plant with Homomorphic Sex Chromosomes

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua, a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

A Comparison of Selective Pressures in Plant X-Linked and Autosomal Genes

Selection is expected to work differently in autosomal and X-linked genes because of their ploidy difference and the exposure of recessive X-linked mutations to haploid selection in males. However, it is not clear whether these expectations apply to recently evolved sex chromosomes, where many genes retain functional X- and Y-linked gametologs. We took advantage of the recently evolved sex chromosomes in the plant Silene latifolia and its closely related species to compare the selective pressures between hemizygous and non-hemizygous X-linked genes as well as between X-linked genes and autosomal genes. Our analysis, based on over 1000 genes, demonstrated that, similar to animals, X-linked genes in Silene evolve significantly faster than autosomal genes—the so-called faster-X effect. Contrary to expectations, faster-X divergence was detectable only for non-hemizygous X-linked genes. Our phylogeny-based analyses of selection revealed no evidence for faster adaptation in X-linked genes compared to autosomal genes. On the other hand, partial relaxation of purifying selection was apparent on the X-chromosome compared to the autosomes, consistent with a smaller genetic diversity in S. latifolia X-linked genes (π_x = 0.016; π_aut = 0.023). Thus, the faster-X divergence in S. latifolia appears to be a consequence of the smaller effective population size rather than of a faster adaptive evolution on the X-chromosome. We argue that this may be a general feature of “young” sex chromosomes, where the majority of X-linked genes are not hemizygous, preventing haploid selection in heterogametic sex.

The slowdown of Y chromosome expansion in dioecious Silene latifolia due to DNA loss and male-specific silencing of retrotransposons

Background

The rise and fall of the Y chromosome was demonstrated in animals but plants often possess the large evolutionarily young Y chromosome that is thought has expanded recently. Break-even points dividing expansion and shrinkage phase of plant Y chromosome evolution are still to be determined. To assess the size dynamics of the Y chromosome, we studied intraspecific genome size variation and genome composition of male and female individuals in a dioecious plant Silene latifolia, a well-established model for sex-chromosomes evolution.

Results

Our genome size data are the first to demonstrate that regardless of intraspecific genome size variation, Y chromosome has retained its size in S. latifolia. Bioinformatics study of genome composition showed that constancy of Y chromosome size was caused by Y chromosome DNA loss and the female-specific proliferation of recently active dominant retrotransposons. We show that several families of retrotransposons have contributed to genome size variation but not to Y chromosome size change.

Conclusions

Our results suggest that the large Y chromosome of S. latifolia has slowed down or stopped its expansion. Female-specific proliferation of retrotransposons, enlarging the genome with exception of the Y chromosome, was probably caused by silencing of highly active retrotransposons in males and represents an adaptive mechanism to suppress degenerative processes in the haploid stage. Sex specific silencing of transposons might be widespread in plants but hidden in traditional hermaphroditic model plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4547-7) contains supplementary material, which is available to authorized users.

Identification of the Sex-Biased Gene Expression and Putative Sex-Associated Genes in Eucommia ulmoides Oliver Using Comparative Transcriptome Analyses

Eucommia ulmoides is a model representative of the dioecious plants with sex differentiation at initiation. Nevertheless, the genetic mechanisms of sexual dimorphism and sex determination in E. ulmoides remain poorly understood. In this study de novo transcriptome sequencing on Illumina platform generated >45 billion high-quality bases from fresh leaves of six male and female individuals of E. ulmoides. A total of 148,595 unigenes with an average length of 801 base-pairs (bp) were assembled. Through comparative transcriptome analyses, 116 differentially expressed genes (DEGs) between the males and the females were detected, including 73 male-biased genes and 43 female-biased genes. Of these DEGs, three female-biased genes were annotated to be related with the sexually dimorphic gutta content in E. ulmoides. One male-biased DEG was identified as putative MADS box gene APETALA3, a B class floral organ identity gene in the flowering plants. SNPs calling analyses further confirmed that the APETALA3-like gene was probably involved in the sex determination in E. ulmoides. Four other male-biased DEGs were potential sex-associated genes as well with segregated SNPs in accord with sex type. In addition, the SNPs density was 1.02 per kilobase (kb) in the expressed genes of E. ulmoides, implying a relatively high genetic diversity.

Evidence for a genetic sex determination in Cnidaria, the Mediterranean red coral (Corallium rubrum)

Sexual reproduction is widespread among eukaryotes, and the sex-determining processes vary greatly among species. While genetic sex determination (GSD) has been intensively described in bilaterian species, no example has yet been recorded among non-bilaterians. However, the quasi-ubiquitous repartition of GSD among multicellular species suggests that similar evolutionary forces can promote this system, and that these forces could occur also in non-bilaterians. Studying sex determination across the range of Metazoan diversity is indeed important to understand better the evolution of this mechanism and its lability. We tested the existence of sex-linked genes in the gonochoric red coral (Corallium rubrum, Cnidaria) using restriction site-associated DNA sequencing. We analysed 27 461 single nucleotide polymorphisms (SNPs) in 354 individuals from 12 populations including 53 that were morphologically sexed. We found a strong association between the allele frequencies of 472 SNPs and the sex of individuals, suggesting an XX/XY sex-determination system. This result was confirmed by the identification of 435 male-specific loci. An independent test confirmed that the amplification of these loci enabled us to identify males with absolute certainty. This is the first demonstration of a GSD system among non-bilaterian species and a new example of its convergence in multicellular eukaryotes.

Evidence for a genetic sex determination in Cnidaria, the Mediterranean red coral (Corallium rubrum)

Sexual reproduction is widespread among eukaryotes, and the sex-determining processes vary greatly among species. While genetic sex determination (GSD) has been intensively described in bilaterian species, no example has yet been recorded among non-bilaterians. However, the quasi-ubiquitous repartition of GSD among multicellular species suggests that similar evolutionary forces can promote this system, and that these forces could occur also in non-bilaterians. Studying sex determination across the range of Metazoan diversity is indeed important to understand better the evolution of this mechanism and its lability. We tested the existence of sex-linked genes in the gonochoric red coral (Corallium rubrum, Cnidaria) using restriction site-associated DNA sequencing. We analysed 27 461 single nucleotide polymorphisms (SNPs) in 354 individuals from 12 populations including 53 that were morphologically sexed. We found a strong association between the allele frequencies of 472 SNPs and the sex of individuals, suggesting an XX/XY sex-determination system. This result was confirmed by the identification of 435 male-specific loci. An independent test confirmed that the amplification of these loci enabled us to identify males with absolute certainty. This is the first demonstration of a GSD system among non-bilaterian species and a new example of its convergence in multicellular eukaryotes.

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.

Homomorphic ZW chromosomes in a wild strawberry show distinctive recombination heterogeneity but a small sex-determining region

Recombination in ancient, heteromorphic sex chromosomes is typically suppressed at the sex-determining region (SDR) and proportionally elevated in the pseudoautosomal region (PAR). However, little is known about recombination dynamics of young, homomorphic plant sex chromosomes. We examine male and female function in crosses and unrelated samples of the dioecious octoploid strawberry Fragaria chiloensis in order to map the small and recently evolved SDR controlling both traits and to examine recombination patterns on the incipient ZW chromosome. The SDR of this ZW system is located within a 280 kb window, in which the maternal recombination rate is lower than the paternal one. In contrast to the SDR, the maternal PAR recombination rate is much higher than the rates of the paternal PAR or autosomes, culminating in an elevated chromosome-wide rate. W-specific divergence is elevated within the SDR and a single polymorphism is observed in high species-wide linkage disequilibrium with sex. Selection for recombination suppression within the small SDR may be weak, but fluctuating sex ratios could favor elevated recombination in the PAR to remove deleterious mutations on the W. The recombination dynamics of this nascent sex chromosome with a modestly diverged SDR may be typical of other dioecious plants.

Comparison of Spinach Sex Chromosomes with Sugar Beet Autosomes Reveals Extensive Synteny and Low Recombination at the Male-Determining Locus

Spinach (Spinacia oleracea, 2n = 12) and sugar beet (Beta vulgaris, 2n = 18) are important crop members of the family Chenopodiaceae ss Sugar beet has a basic chromosome number of 9 and a cosexual breeding system, as do most members of the Chenopodiaceae ss. family. By contrast, spinach has a basic chromosome number of 6 and, although certain cultivars and genotypes produce monoecious plants, is considered to be a dioecious species. The loci determining male and monoecious sexual expression were mapped to different loci on the spinach sex chromosomes. In this study, a linkage map with 46 mapped protein-coding sequences was constructed for the spinach sex chromosomes. Comparison of the linkage map with a reference genome sequence of sugar beet revealed that the spinach sex chromosomes exhibited extensive synteny with sugar beet chromosomes 4 and 9. Tightly linked protein-coding genes linked to the male-determining locus in spinach corresponded to genes located in or around the putative pericentromeric and centromeric regions of sugar beet chromosomes 4 and 9, supporting the observation that recombination rates were low in the vicinity of the male-determining locus. The locus for monoecism was confined to a chromosomal segment corresponding to a region of approximately 1.7Mb on sugar beet chromosome 9, which may facilitate future positional cloning of the locus.

Can the experimental evolution programme help us elucidate the genetic basis of adaptation in nature?

There have been a variety of approaches taken to try to characterize and identify the genetic basis of adaptation in nature, spanning theoretical models, experimental evolution studies and direct tests of natural populations. Theoretical models can provide formalized and detailed hypotheses regarding evolutionary processes and patterns, from which experimental evolution studies can then provide important proofs of concepts and characterize what is biologically reasonable. Genetic and genomic data from natural populations then allow for the identification of the particular factors that have and continue to play an important role in shaping adaptive evolution in the natural world. Further to this, experimental evolution studies allow for tests of theories that may be difficult or impossible to test in natural populations for logistical and methodological reasons and can even generate new insights, suggesting further refinement of existing theories. However, as experimental evolution studies often take place in a very particular set of controlled conditions--that is simple environments, a small range of usually asexual species, relatively short timescales--the question remains as to how applicable these experimental results are to natural populations. In this review, we discuss important insights coming from experimental evolution, focusing on four key topics tied to the evolutionary genetics of adaptation, and within those topics, we discuss the extent to which the experimental work compliments and informs natural population studies. We finish by making suggestions for future work in particular a need for natural population genomic time series data, as well as the necessity for studies that combine both experimental evolution and natural population approaches.

Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1698-7) contains supplementary material, which is available to authorized users.