Expansion of the pseudo-autosomal region and ongoing recombination suppression in the Silene latifolia sex chromosomes

Rapid and dynamic evolution of a giant Y chromosome in Silene latifolia

Some plants have massive sex-linked regions. To test hypotheses about their evolution, we sequenced the genome of Silene latifolia, in which giant heteromorphic sex chromosomes were first discovered in 1923. It has long been known that the Y chromosome consists mainly of a male-specific region that does not recombine with the X chromosome and carries the sex-determining genes and genes with other male functions. However, only with a whole Y chromosome assembly can candidate genes be validated experimentally and their locations determined and related to the suppression of recombination. We describe the genomic changes as the ancestral chromosome evolved into the current XY pair, testing ideas about the evolution of large nonrecombining regions and the mechanisms that created the present recombination pattern.

The Silene latifolia genome and its giant Y chromosome

In many species with sex chromosomes, the Y is a tiny chromosome. However, the dioecious plant Silene latifolia has a giant ~550-megabase Y chromosome, which has remained unsequenced so far. We used a long- and short-read hybrid approach to obtain a high-quality male genome. Comparative analysis of the sex chromosomes with their homologs in outgroups showed that the Y is highly rearranged and degenerated. Recombination suppression between X and Y extended in several steps and triggered a massive accumulation of repeats on the Y as well as in the nonrecombining pericentromeric region of the X, leading to giant sex chromosomes. Using sex phenotype mutants, we identified candidate sex-determining genes on the Y in locations consistent with their favoring recombination suppression events 11 and 5 million years ago.

Why should we study plant sex chromosomes?

Understanding plant sex chromosomes involves studying interactions between developmental and physiological genetics, genome evolution, and evolutionary ecology. We focus on areas of overlap between these. Ideas about how species with separate sexes (dioecious species, in plant terminology) can evolve are even more relevant to plants than to most animal taxa because dioecy has evolved many times from ancestral functionally hermaphroditic populations, often recently. One aim of studying plant sex chromosomes is to discover how separate males and females evolved from ancestors with no such genetic sex-determining polymorphism, and the diversity in the genetic control of maleness vs femaleness. Different systems share some interesting features, and their differences help to understand why completely sex-linked regions may evolve. In some dioecious plants, the sex-determining genome regions are physically small. In others, regions without crossing over have evolved sometimes extensive regions with properties very similar to those of the familiar animal sex chromosomes. The differences also affect the evolutionary changes possible when the environment (or pollination environment, for angiosperms) changes, as dioecy is an ecologically risky strategy for sessile organisms. Dioecious plants have repeatedly reverted to cosexuality, and hermaphroditic strains of fruit crops such as papaya and grapes are desired by plant breeders. Sex-linked regions are predicted to become enriched in genes with sex differences in expression, especially when higher expression benefits one sex function but harms the other. Such trade-offs may be important for understanding other plant developmental and physiological processes and have direct applications in plant breeding.

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.

Evolutionary History of the Poecilia picta Sex Chromosomes

The degree of divergence between the sex chromosomes is not always proportional to their age. In poeciliids, four closely related species all exhibit a male heterogametic sex chromosome system on the same linkage group, yet show a remarkable diversity in X and Y divergence. In Poecilia reticulata and P. wingei, the sex chromosomes remain homomorphic, yet P. picta and P. parae have a highly degraded Y chromosome. To test alternative theories about the origin of their sex chromosomes, we used a combination of pedigrees and RNA-seq data from P. picta families in conjunction with DNA-seq data collected from P. reticulata, P. wingei, P. parae, and P. picta. Phylogenetic clustering analysis of X and Y orthologs, identified through segregation patterns, and their orthologous sequences in closely related species demonstrates a similar time of origin for both the P. picta and P. reticulata sex chromosomes. We next used k-mer analysis to identify shared ancestral Y sequence across all four species, suggesting a single origin to the sex chromosome system in this group. Together, our results provide key insights into the origin and evolution of the poeciliid Y chromosome and illustrate that the rate of sex chromosome divergence is often highly heterogenous, even over relatively short evolutionary time frames.

Heterochiasmy and Sex Chromosome Evolution in Silene

The evolution of a non-recombining sex-specific region is a key step in sex chromosome evolution. Suppression of recombination between the (proto-) X- and Y-chromosomes in male meiosis creates a non-recombining Y-linked region (NRY), while the X-chromosome continues to recombine in females. Lack of recombination in the NRY defines its main properties—genetic degeneration and accumulation of repetitive DNA, making X and Y chromosomes very different from each other. How and why recombination suppression on sex chromosomes evolves remains controversial. A strong difference in recombination rates between the sexes (heterochiasmy) can facilitate or even cause recombination suppression. In the extreme case—complete lack of recombination in the heterogametic sex (achiasmy)—the entire sex-specific chromosome is automatically non-recombining. In this study, I analyse sex-specific recombination rates in a dioecious plant Silene latifolia (Caryophyllaceae), which evolved separate sexes and sex chromosomes ~11 million years ago. I reconstruct high-density RNAseq-based genetic maps including over five thousand genic markers for the two sexes separately. The comparison of the male and female maps reveals only modest heterochiasmy across the genome, with the exception of the sex chromosomes, where recombination is suppressed in males. This indicates that heterochiasmy likely played only a minor, if any, role in NRY evolution in S. latifolia, as recombination suppression is specific to NRY rather than to the entire genome in males. Other mechanisms such as structural rearrangements and/or epigenetic modifications were likely involved, and comparative genome analysis and genetic mapping in multiple Silene species will help to shed light on the mechanism(s) of recombination suppression that led to the evolution of sex chromosomes.

Evolution of plant sex and molecular mechanisms underlying plants sex separation

Unlike animals, plants exhibit more complexity of sexual morphs. The genetic mechanism underlying plant sex is a hot research topic in plant biology. In recent decades, advanced theories have been put forth on plant sex determination, but experimental proof is scarce. In recent years, vast achievements have been made to reveal the genetic mechanisms underlying sex separation of plants at the molecular level. Although the sex determination mechanisms have been clarified only in a limited number of plant species thus far, the discoveries offer us an opportunity to understand the genetic mechanisms triggering the separation of plant sexes. This paper reviewed the different aspects of the advanced studies on plant sex evolution and the molecular mechanisms underlying plant sex separation.

Comparison of methodological approaches to the study of young sex chromosomes: A case study in Poecilia

Studies of sex chromosome systems at early stages of divergence are key to understanding the initial process and underlying causes of recombination suppression. However, identifying signatures of divergence in homomorphic sex chromosomes can be challenging due to high levels of sequence similarity between the X and the Y. Variations in methodological precision and underlying data can make all the difference between detecting subtle divergence patterns or missing them entirely. Recent efforts to test for X-Y sequence differentiation in the guppy have led to contradictory results. Here, we apply different analytical methodologies to the same data set to test for the accuracy of different approaches in identifying patterns of sex chromosome divergence in the guppy. Our comparative analysis reveals that the most substantial source of variation in the results of the different analyses lies in the reference genome used. Analyses using custom-made genome assemblies for the focal population or species successfully recover a signal of divergence across different methodological approaches. By contrast, using the distantly related Xiphophorus reference genome results in variable patterns, due to both sequence evolution and structural variations on the sex chromosomes between the guppy and Xiphophorus. Changes in mapping and filtering parameters can additionally introduce noise and obscure the signal. Our results illustrate how analytical differences can alter perceived results and we highlight best practices for the study of nascent sex chromosomes.

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.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops

In sharp contrast to birds and mammals, in numerous cold-blooded vertebrates, sex chromosomes have been described as homomorphic. This sex chromosome homomorphy has been suggested to result from the high turnovers often observed across deeply diverged clades. However, little is known about the tempo and mode of sex chromosome evolution among the most closely related species. Here, we examined the evolution of sex chromosome among nine species of the torrent frog genus Amolops. We analyzed male and female GBS and RAD-seq from 182 individuals and performed PCR verification for 176 individuals. We identified signatures of sex chromosomes involving two pairs of chromosomes. We found that sex-chromosome homomorphy results from both turnover and X–Y recombination in the Amolops species, which simultaneously exhibits heterogeneous evolution on homologous and non-homologous sex chromosomes. A low turnover rate of non-homologous sex chromosomes exists in these torrent frogs. The ongoing X–Y recombination in homologous sex chromosomes will act as an indispensable force in preventing sex chromosomes from differentiating.

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.

Gene-Level, but Not Chromosome-Wide, Divergence between a Very Young House Fly Proto-Y Chromosome and Its Homologous Proto-X Chromosome

X and Y chromosomes are usually derived from a pair of homologous autosomes, which then diverge from each other over time. Although Y-specific features have been characterized in sex chromosomes of various ages, the earliest stages of Y chromosome evolution remain elusive. In particular, we do not know whether early stages of Y chromosome evolution consist of changes to individual genes or happen via chromosome-scale divergence from the X. To address this question, we quantified divergence between young proto-X and proto-Y chromosomes in the house fly, Musca domestica. We compared proto-sex chromosome sequence and gene expression between genotypic (XY) and sex-reversed (XX) males. We find evidence for sequence divergence between genes on the proto-X and proto-Y, including five genes with mitochondrial functions. There is also an excess of genes with divergent expression between the proto-X and proto-Y, but the number of genes is small. This suggests that individual proto-Y genes, but not the entire proto-Y chromosome, have diverged from the proto-X. We identified one gene, encoding an axonemal dynein assembly factor (which functions in sperm motility), that has higher expression in XY males than XX males because of a disproportionate contribution of the proto-Y allele to gene expression. The upregulation of the proto-Y allele may be favored in males because of this gene's function in spermatogenesis. The evolutionary divergence between proto-X and proto-Y copies of this gene, as well as the mitochondrial genes, is consistent with selection in males affecting the evolution of individual genes during early Y chromosome evolution.

Divergence and Remarkable Diversity of the Y Chromosome in Guppies

The guppy sex chromosomes show an extraordinary diversity in divergence across populations and closely related species. In order to understand the dynamics of the guppy Y chromosome, we used linked-read sequencing to assess Y chromosome evolution and diversity across upstream and downstream population pairs that vary in predator and food abundance in three replicate watersheds. Based on our population-specific genome assemblies, we first confirmed and extended earlier reports of two strata on the guppy sex chromosomes. Stratum I shows significant accumulation of male-specific sequence, consistent with Y divergence, and predates the colonization of Trinidad. In contrast, Stratum II shows divergence from the X, but no Y-specific sequence, and this divergence is greater in three replicate upstream populations compared with their downstream pair. Despite longstanding assumptions that sex chromosome recombination suppression is achieved through inversions, we find no evidence of inversions associated with either Stratum I or Stratum II. Instead, we observe a remarkable diversity in Y chromosome haplotypes within each population, even in the ancestral Stratum I. This diversity is likely due to gradual mechanisms of recombination suppression, which, unlike an inversion, allow for the maintenance of multiple haplotypes. In addition, we show that this Y diversity is dominated by low-frequency haplotypes segregating in the population, suggesting a link between haplotype diversity and female preference for rare Y-linked color variation. Our results reveal the complex interplay between recombination suppression and Y chromosome divergence at the earliest stages of sex chromosome divergence.

Histone post-translational modifications in Silene latifolia X and Y chromosomes suggest a mammal-like dosage compensation system

Silene latifolia is a model organism to study evolutionary young heteromorphic sex chromosome evolution in plants. Previous research indicates a Y-allele gene degeneration and a dosage compensation system already operating. Here, we propose an epigenetic approach based on analysis of several histone post-translational modifications (PTMs) to find the first epigenetic hints of the X:Y sex chromosome system regulation in S. latifolia. Through chromatin immunoprecipitation we interrogated six genes from X and Y alleles. Several histone PTMS linked to DNA methylation and transcriptional repression (H3K27me3, H3K23me, H3K9me2 and H3K9me3) and to transcriptional activation (H3K4me3 and H4K5, 8, 12, 16ac) were used. DNA enrichment (Immunoprecipitated DNA/input DNA) was analyzed and showed three main results: (i) promoters of the Y allele are associated with heterochromatin marks, (ii) promoters of the X allele in males are associated with activation of transcription marks and finally, (iii) promoters of X alleles in females are associated with active and repressive marks. Our finding indicates a transcription activation of X allele and transcription repression of Y allele in males. In females we found a possible differential regulation (up X₁, down X₂) of each female X allele. These results agree with the mammal-like epigenetic dosage compensation regulation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Guppy Y Chromosome Integrity Maintained by Incomplete Recombination Suppression

The loss of recombination triggers divergence between the sex chromosomes and promotes degeneration of the sex-limited chromosome. Several livebearers within the genus Poecilia share a male-heterogametic sex chromosome system that is roughly 20 Myr old, with extreme variation in the degree of Y chromosome divergence. In Poecilia picta, the Y is highly degenerate and associated with complete X chromosome dosage compensation. In contrast, although recombination is restricted across almost the entire length of the sex chromosomes in Poecilia reticulata and Poecilia wingei, divergence between the X chromosome and the Y chromosome is very low. This clade therefore offers a unique opportunity to study the forces that accelerate or hinder sex chromosome divergence. We used RNA-seq data from multiple families of both P. reticulata and P. wingei, the species with low levels of sex chromosome divergence, to differentiate X and Y coding sequences based on sex-limited SNP inheritance. Phylogenetic tree analyses reveal that occasional recombination has persisted between the sex chromosomes for much of their length, as X- and Y-linked sequences cluster by species instead of by gametolog. This incomplete recombination suppression maintains the extensive homomorphy observed in these systems. In addition, we see differences between the previously identified strata in the phylogenetic clustering of X-Y orthologs, with those that cluster by chromosome located in the older stratum, the region previously associated with the sex-determining locus. However, recombination arrest appears to have expanded throughout the sex chromosomes more gradually instead of through a stepwise process associated with inversions.

Sex Chromosome Evolution: So Many Exceptions to the Rules

Genomic analysis of many nonmodel species has uncovered an incredible diversity of sex chromosome systems, making it possible to empirically test the rich body of evolutionary theory that describes each stage of sex chromosome evolution. Classic theory predicts that sex chromosomes originate from a pair of homologous autosomes and recombination between them is suppressed via inversions to resolve sexual conflict. The resulting degradation of the Y chromosome gene content creates the need for dosage compensation in the heterogametic sex. Sex chromosome theory also implies a linear process, starting from sex chromosome origin and progressing to heteromorphism. Despite many convergent genomic patterns exhibited by independently evolved sex chromosome systems, and many case studies supporting these theoretical predictions, emerging data provide numerous interesting exceptions to these long-standing theories, and suggest that the remarkable diversity of sex chromosomes is matched by a similar diversity in their evolution. For example, it is clear that sex chromosome pairs are not always derived from homologous autosomes. In addition, both the cause and the mechanism of recombination suppression between sex chromosome pairs remain unclear, and it may be that the spread of recombination suppression is a more gradual process than previously thought. It is also clear that dosage compensation can be achieved in many ways, and displays a range of efficacy in different systems. Finally, the remarkable turnover of sex chromosomes in many systems, as well as variation in the rate of sex chromosome divergence, suggest that assumptions about the inevitable linearity of sex chromosome evolution are not always empirically supported, and the drivers of the birth-death cycle of sex chromosome evolution remain to be elucidated. Here, we concentrate on how the diversity in sex chromosomes across taxa highlights an equal diversity in each stage of sex chromosome evolution.

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.

The Genomic Complexity of a Large Inversion in Great Tits

Chromosome inversions have clear effects on genome evolution and have been associated with speciation, adaptation, and the evolution of the sex chromosomes. In birds, these inversions may play an important role in hybridization of species and disassortative mating. We identified a large (≈64 Mb) inversion polymorphism in the great tit (Parus major) that encompasses almost 1,000 genes and more than 90% of Chromosome 1A. The inversion occurs at a low frequency in a set of over 2,300 genotyped great tits in the Netherlands with only 5% of the birds being heterozygous for the inversion. In an additional analysis of 29 resequenced birds from across Europe, we found two heterozygotes. The likely inversion breakpoints show considerable genomic complexity, including multiple copy number variable segments. We identified different haplotypes for the inversion, which differ in the degree of recombination in the center of the chromosome. Overall, this remarkable genetic variant is widespread among distinct great tit populations and future studies of the inversion haplotype, including how it affects the fitness of carriers, may help to understand the mechanisms that maintain it.

The Formation of Sex Chromosomes in Silene latifolia and S. dioica Was Accompanied by Multiple Chromosomal Rearrangements

The genus Silene includes a plethora of dioecious and gynodioecious species. Two species, Silene latifolia (white campion) and Silene dioica (red campion), are dioecious plants, having heteromorphic sex chromosomes with an XX/XY sex determination system. The X and Y chromosomes differ mainly in size, DNA content and posttranslational histone modifications. Although it is generally assumed that the sex chromosomes evolved from a single pair of autosomes, it is difficult to distinguish the ancestral pair of chromosomes in related gynodioecious and hermaphroditic plants. We designed an oligo painting probe enriched for X-linked scaffolds from currently available genomic data and used this probe on metaphase chromosomes of S. latifolia (2n = 24, XY), S. dioica (2n = 24, XY), and two gynodioecious species, S. vulgaris (2n = 24) and S. maritima (2n = 24). The X chromosome-specific oligo probe produces a signal specifically on the X and Y chromosomes in S. latifolia and S. dioica, mainly in the subtelomeric regions. Surprisingly, in S. vulgaris and S. maritima, the probe hybridized to three pairs of autosomes labeling their p-arms. This distribution suggests that sex chromosome evolution was accompanied by extensive chromosomal rearrangements in studied dioecious plants.

Immediate Dosage Compensation Is Triggered by the Deletion of Y-Linked Genes in Silene latifolia

The loss of functional genes from non-recombining sex-specific chromosomes [1, 2], such as the Y chromosomes in mammals [3] or W chromosomes in birds [4], should result in an imbalance of gene products for sex-linked genes [5]. Different chromosome-wide systems that rebalance gene expression are known to operate in organisms with relatively old sex chromosomes [6]; e.g., Drosophila overexpress X-linked genes in males [7], while mammals shut down one of the X chromosomes in females [8]. It is not known how long it takes for a chromosome-wide dosage compensation system to evolve. To shed light on the early evolution of dosage compensation, we constructed a high-density Y-deletion map and used deletion mutants to manipulate gene dose and analyze gene expression in white campion (Silene latifolia), which evolved dioecy and sex chromosomes only 11 million years ago [9]. We demonstrate that immediate dosage compensation can be triggered by deletions in a large portion of the p arm of the Y chromosome. Our results indicate that dosage compensation in S. latifolia does not have to evolve gene by gene because a system to upregulate gene expression is already operating on part of the X chromosome, which likely represents an intermediate step in the evolution of a chromosome-wide dosage compensation system in this species.

Evolution of sex determination and heterogamety changes in section Otites of the genus Silene

Switches in heterogamety are known to occur in both animals and plants. Although plant sex determination systems probably often evolved more recently than those in several well-studied animals, including mammals, and have had less time for switches to occur, we previously detected a switch in heterogamety in the plant genus Silene: section Otites has both female and male heterogamety, whereas S. latifolia and its close relatives, in a different section of the genus, Melandrium (subgenus Behenantha), all have male heterogamety. Here we analyse the evolution of sex chromosomes in section Otites, which is estimated to have evolved only about 0.55 MYA. Our study confirms female heterogamety in S. otites and newly reveals female heterogamety in S. borysthenica. Sequence analyses and genetic mapping show that the sex-linked regions of these two species are the same, but the region in S. colpophylla, a close relative with male heterogamety, is different. The sex chromosome pairs of S. colpophylla and S. otites each correspond to an autosome of the other species, and both differ from the XY pair in S. latifolia. Silene section Otites species are suitable for detailed studies of the events involved in such changes, and our phylogenetic analysis suggests a possible change from female to male heterogamety within this section. Our analyses suggest a possibility that has so far not been considered, change in heterogamety through hybridization, in which a male-determining chromosome from one species is introgressed into another one, and over-rides its previous sex-determining system.

A Genetic Map of Ostrich Z Chromosome and the Role of Inversions in Avian Sex Chromosome Evolution

Recombination arrest is a necessary step for the evolution of distinct sex chromosomes. Structural changes, such as inversions, may represent the mechanistic basis for recombination suppression and comparisons of the structural organization of chromosomes as given by chromosome-level assemblies offer the possibility to infer inversions across species at some detail. In birds, deduction of the process of sex chromosome evolution has been hampered by the lack of a validated chromosome-level assembly from a representative of one of the two basal clades of modern birds, Paleognathae. We therefore developed a high-density genetic linkage map of the ostrich Z chromosome and used this to correct an existing assembly, including correction of a large chimeric superscaffold and the order and orientation of other superscaffolds. We identified the pseudoautosomal region as a 52 Mb segment (≈60% of the Z chromosome) where recombination occurred in both sexes. By comparing the order and location of genes on the ostrich Z chromosome with that of six bird species from the other major clade of birds (Neognathae), and of reptilian outgroup species, 25 Z-linked inversions were inferred in the avian lineages. We defined Z chromosome organization in an early avian ancestor and identified inversions spanning the candidate sex-determining DMRT1 gene in this ancestor, which could potentially have triggered the onset of avian sex chromosome evolution. We conclude that avian sex chromosome evolution has been characterized by a complex process of probably both Z-linked and W-linked inversions (and/or other processes). This study illustrates the need for validated chromosome-level assemblies for inference of genome evolution.

Role of recombination and faithfulness to partner in sex chromosome degeneration

Sex determination in mammals is strongly linked to sex chromosomes. In most cases, females possess two copies of X chromosome while males have one X and one Y chromosome. It is assumed that these chromosomes originated from a pair of homologous autosomes, which diverged when recombination between them was suppressed. However, it is still debated why the sex chromosomes stopped recombining and how this process spread out over most part of the chromosomes. To study this problem, we developed a simulation model, in which the recombination rate between the sex chromosomes can freely evolve. We found that the suppression of recombination between the X and Y is spontaneous and proceeds very quickly during the evolution of population, which leads to the degeneration of the Y in males. Interestingly, the degeneration happens only when mating pairs are unfaithful. This evolutionary strategy purifies the X chromosome from defective alleles and leads to the larger number of females than males in the population. In consequence, the reproductive potential of the whole population increases. Our results imply that both the suppression of recombination and the degeneration of Y chromosome may be associated with reproductive strategy and favoured in polygamous populations with faithless mating partners.

The Mutation Rate and the Age of the Sex Chromosomes in Silene latifolia

Many aspects of sex chromosome evolution are common to both plants and animals [1], but the process of Y chromosome degeneration, where genes on the Y become non-functional over time, may be much slower in plants due to purifying selection against deleterious mutations in the haploid gametophyte [2, 3]. Testing for differences in Y degeneration between the kingdoms has been hindered by the absence of accurate age estimates for plant sex chromosomes. Here, we used genome resequencing to estimate the spontaneous mutation rate and the age of the sex chromosomes in white campion (Silene latifolia). Screening of single nucleotide polymorphisms (SNPs) in parents and 10 F₁ progeny identified 39 de novo mutations and yielded a rate of 7.31 × 10^-9 (95% confidence interval: 5.20 × 10^-9 - 8.00 × 10^-9) mutations per site per haploid genome per generation. Applying this mutation rate to the synonymous divergence between homologous X- and Y-linked genes (gametologs) gave age estimates of 11.00 and 6.32 million years for the old and young strata, respectively. Based on SNP segregation patterns, we inferred which genes were Y-linked and found that at least 47% are already dysfunctional. Applying our new estimates for the age of the sex chromosomes indicates that the rate of Y degeneration in S. latifolia is nearly 2-fold slower when compared to animal sex chromosomes of a similar age. Our revised estimates support Y degeneration taking place more slowly in plants, a discrepancy that may be explained by differences in the life cycles of animals and plants.

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.

Agrobacterium rhizogenes-mediated transformation of a dioecious plant model Silene latifolia

Silene latifolia serves as a model species to study dioecy, the evolution of sex chromosomes, dosage compensation and sex-determination systems in plants. Currently, no protocol for genetic transformation is available for this species, mainly because S. latifolia is considered recalcitrant to in vitro regeneration and infection with Agrobacterium tumefaciens. Using cytokinins and their synthetic derivatives, we markedly improved the efficiency of regeneration. Several agrobacterial strains were tested for their ability to deliver DNA into S. latifolia tissues leading to transient and stable expression of the GUS reporter. The use of Agrobacterium rhizogenes strains resulted in the highest transformation efficiency (up to 4.7% of stable transformants) in hairy root cultures. Phenotypic and genotypic analyses of the T1 generation suggested that the majority of transformation events contain a small number of independent T-DNA insertions and the transgenes are transmitted to the progeny in a Mendelian pattern of inheritance. In short, we report an efficient and reproducible protocol for leaf disc transformation and subsequent plant regeneration in S. latifolia, based on the unique combination of infection with A. rhizogenes and plant regeneration from hairy root cultures using synthetic cytokinins. A protocol for the transient transformation of S.latifolia protoplasts was also developed and applied to demonstrate the possibility of targeted mutagenesis of the sex linked gene SlAP3 by TALENs and CRISPR/Cas9.

The Evolution of Sex Chromosomes and Dosage Compensation in Plants

Plant sex chromosomes can be vastly different from those of the few historical animal model organisms from which most of our understanding of sex chromosome evolution is derived. Recently, we have seen several advancements from studies on green algae, brown algae, and land plants that are providing a broader understanding of the variable ways in which sex chromosomes can evolve in distant eukaryotic groups. Plant sex-determining genes are being identified and, as expected, are completely different from those in animals. Species with varying levels of differentiation between the X and Y have been found in plants, and these are hypothesized to be representing different stages of sex chromosome evolution. However, we are also finding that sex chromosomes can remain morphologically unchanged over extended periods of time. Where degeneration of the Y occurs, it appears to proceed similarly in plants and animals. Dosage compensation (a phenomenon that compensates for the consequent loss of expression from the Y) has now been documented in a plant system, its mechanism, however, remains unknown. Research has also begun on the role of sex chromosomes in sexual conflict resolution, and it appears that sex-biased genes evolve similarly in plants and animals, although the functions of these genes remain poorly studied. Because the difficulty in obtaining sex chromosome sequences is increasingly being overcome by methodological developments, there is great potential for further discovery within the field of plant sex chromosome evolution.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Sex-oriented research on dioecious crops of Indian subcontinent: an updated review

A number of dioecious species are grown across India and some of those plants play a crucial role in the agro-based economy of the country. The diagnosis of sex is very difficult in the dioecious plant prior flowering wherein sex identification at the seedling stage is of great importance to breeders as well as farmers for crop improvement or production purpose. A comprehensive approach of sex determination comprising morphological, biochemical, cytological and molecular attributes is a must required for gender differentiation in dioecious plant species. In the present review, we highlighted the economical, medicinal as well as industrial importance of most of the dioecious species extensively grown in Indian subcontinent. In addition to that, the cytogenetic, genetic as well as molecular information in connection to their sex determination were critically discussed in this review.

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.

Fruitful analysis of sex chromosomes reveals X-treme genetic diversity

A new study on sex chromosome evolution in papaya helps to illuminate sex chromosome biology, including deviations from expected trajectories.

Please see related Research article: https://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-1095-9

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.

How to make a sex chromosome

Sex chromosomes can evolve once recombination is halted between a homologous pair of chromosomes. Owing to detailed studies using key model systems, we have a nuanced understanding and a rich review literature of what happens to sex chromosomes once recombination is arrested. However, three broad questions remain unanswered. First, why do sex chromosomes stop recombining in the first place? Second, how is recombination halted? Finally, why does the spread of recombination suppression, and therefore the rate of sex chromosome divergence, vary so substantially across clades? In this review, we consider each of these three questions in turn to address fundamental questions in the field, summarize our current understanding, and highlight important areas for future work.

Sex chromosome evolution begins when recombination between a homologous pair of chromosomes is halted. Here, Wright et al. review our current understanding of the causes and mechanisms of recombination suppression between incipient sex chromosomes and suggest future directions for the field.

Plant Sex Chromosomes

Although individuals in most flowering plant species, and in many haploid plants, have both sex functions, dioecious species-in which individuals have either male or female functions only-are scattered across many taxonomic groups, and many species have genetic sex determination. Among these, some have visibly heteromorphic sex chromosomes, and molecular genetic studies are starting to uncover sex-linked markers in others, showing that they too have fully sex-linked regions that are either too small or are located in chromosomes that are too small to be cytologically detectable from lack of pairing, lack of visible crossovers, or accumulation of heterochromatin. Detailed study is revealing that, like animal sex chromosomes, plant sex-linked regions show evidence for accumulation of repetitive sequences and genetic degeneration. Estimating when recombination stopped confirms the view that many plants have young sex-linked regions, making plants of great interest for studying the timescale of these changes.

Deciphering evolutionary strata on plant sex chromosomes and fungal mating-type chromosomes through compositional segmentation

Sex chromosomes have evolved from a pair of homologous autosomes which differentiated into sex determination systems, such as XY or ZW system, as a consequence of successive recombination suppression between the gametologous chromosomes. Identifying the regions of recombination suppression, namely, the "evolutionary strata", is central to understanding the history and dynamics of sex chromosome evolution. Evolution of sex chromosomes as a consequence of serial recombination suppressions is well-studied for mammals and birds, but not for plants, although 48 dioecious plants have already been reported. Only two plants Silene latifolia and papaya have been studied until now for the presence of evolutionary strata on their X chromosomes, made possible by the sequencing of sex-linked genes on both the X and Y chromosomes, which is a requirement of all current methods that determine stratum structure based on the comparison of gametologous sex chromosomes. To circumvent this limitation and detect strata even if only the sequence of sex chromosome in the homogametic sex (i.e. X or Z chromosome) is available, we have developed an integrated segmentation and clustering method. In application to gene sequences on the papaya X chromosome and protein-coding sequences on the S. latifolia X chromosome, our method could decipher all known evolutionary strata, as reported by previous studies. Our method, after validating on known strata on the papaya and S. latifolia X chromosome, was applied to the chromosome 19 of Populus trichocarpa, an incipient sex chromosome, deciphering two, yet unknown, evolutionary strata. In addition, we applied this approach to the recently sequenced sex chromosome V of the brown alga Ectocarpus sp. that has a haploid sex determination system (UV system) recovering the sex determining and pseudoautosomal regions, and then to the mating-type chromosomes of an anther-smut fungus Microbotryum lychnidis-dioicae predicting five strata in the non-recombining region of both the chromosomes.

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.

A new physical mapping approach refines the sex-determining gene positions on the Silene latifolia Y-chromosome

Sex chromosomes are particularly interesting regions of the genome for both molecular genetics and evolutionary studies; yet, for most species, we lack basic information, such as the gene order along the chromosome. Because they lack recombination, Y-linked genes cannot be mapped genetically, leaving physical mapping as the only option for establishing the extent of synteny and homology with the X chromosome. Here, we developed a novel and general method for deletion mapping of non-recombining regions by solving “the travelling salesman problem”, and evaluate its accuracy using simulated datasets. Unlike the existing radiation hybrid approach, this method allows us to combine deletion mutants from different experiments and sources. We applied our method to a set of newly generated deletion mutants in the dioecious plant Silene latifolia and refined the locations of the sex-determining loci on its Y chromosome map.

Rapid Y degeneration and dosage compensation in plant sex chromosomes

The nonrecombining regions of animal Y chromosomes are known to undergo genetic degeneration, but previous work has failed to reveal large-scale gene degeneration on plant Y chromosomes. Here, we uncover rapid and extensive degeneration of Y-linked genes in a plant species, Silene latifolia, that evolved sex chromosomes de novo in the last 10 million years. Previous transcriptome-based studies of this species missed unexpressed, degenerate Y-linked genes. To identify sex-linked genes, regardless of their expression, we sequenced male and female genomes of S. latifolia and integrated the genomic contigs with a high-density genetic map. This revealed that 45% of Y-linked genes are not expressed, and 23% are interrupted by premature stop codons. This contrasts with X-linked genes, in which only 1.3% of genes contained stop codons and 4.3% of genes were not expressed in males. Loss of functional Y-linked genes is partly compensated for by gene-specific up-regulation of X-linked genes. Our results demonstrate that the rate of genetic degeneration of Y-linked genes in S. latifolia is as fast as in animals, and that the evolutionary trajectories of sex chromosomes are similar in the two kingdoms.

Plant contributions to our understanding of sex chromosome evolution

A minority of angiosperms have male and female flowers separated in distinct individuals (dioecy), and most dioecious plants do not have cytologically different (heteromorphic) sex chromosomes. Plants nevertheless have several advantages for the study of sex chromosome evolution, as genetic sex determination has evolved repeatedly and is often absent in close relatives. I review sex-determining regions in non-model plant species, which may help us to understand when and how (and, potentially, test hypotheses about why) recombination suppression evolves within young sex chromosomes. I emphasize high-throughput sequencing approaches that are increasingly being applied to plants to test for non-recombining regions. These data are particularly illuminating when combined with sequence data that allow phylogenetic analyses, and estimates of when these regions evolved. Together with comparative genetic mapping, this has revealed that sex-determining loci and sex-linked regions evolved independently in many plant lineages, sometimes in closely related dioecious species, and often within the past few million years. In reviewing recent progress, I suggest areas for future work, such as the use of phylogenies to allow the informed choice of outgroup species suitable for inferring the directions of changes, including testing whether Y chromosome-like regions are undergoing genetic degeneration, a predicted consequence of losing recombination.

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

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.

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.

Gene loss from a plant sex chromosome system

Sex chromosomes have evolved independently in numerous animal and plant lineages. After recombination becomes suppressed between two homologous sex chromosomes, genes on the non-recombining Y chromosomes (and W chromosomes in ZW systems) undergo genetic degeneration, losing functions retained by their X- or Z-linked homologs, changing their expression, and becoming lost [1, 2]. Adaptive changes may also occur, both on the non-recombining Y chromosome, to shut down expression of maladapted genes [3], and on the X chromosome (or the Z in ZW systems), which may evolve dosage compensation to increase low expression or compensate for poor protein function in the heterogametic sex [2, 4, 5]. Although empirical approaches to studying genetic degeneration have been developed for model species [3, 6], the onset and dynamics of these changes are still poorly understood, particularly in de novo evolving sex chromosomes. Sex chromosomes of some plants evolved much more recently than those of mammals, birds, and Drosophila [7-9], making them suitable for studying the early stages of genetic degeneration in de novo evolving sex chromosomes. In plants, haploid selection should oppose gene loss from Y chromosomes, but recent work on sex chromosomes of two plant species has estimated that Y-linked transcripts are lacking for 10%-30% of X-linked genes [10-12]. Here, we provide evidence that, in Silene latifolia, this largely involved losses of Y-linked genes, and not suppressed expression of Y-linked alleles, or gene additions to the X chromosome. Our results also suggest that chromosome-wide dosage compensation does not occur in this plant.