Genetic degeneration of old and young Y chromosomes in the flowering plant Rumex hastatulus

Effect of heterozygous deletions on phenotypic changes and dosage compensation in Arabidopsis thaliana

Heterozygous deletions, which include a large number of genes, are often caused by the induction of mutations. The induction of gene dosage compensation should be considered when assessing the effects of heterozygous deletions on phenotypic changes. This mechanism is known to balance the expression levels of genes with different copy numbers in sex chromosomes, but it is also known to operate in autosomes. In the present study, 12 Arabidopsis thaliana BC1 mutants with heterozygous deletions were produced by crossing wild-type Col-0 plants with mutants induced by heavy ion beams. The sizes of the deletions ranged from 50.9 kb to 2.03 Mb, and the number of deleted genes ranged from 8 to 92. Nine of the 12 BC1 mutants showed phenotypic changes in fresh weight 14 days after cultivation or during the flowering period. RNA-sequencing (RNA-seq) analyses of 14-day-old leaves, 40-day-old leaves, and flower buds showed that dosage compensation did not occur in any stage or tissue tested. These results indicate that heterozygous deletions cause phenotypic changes owing to the absence of dosage compensation.

A model of Hill-Robertson interference caused by purifying selection in a nonrecombining genome

A new approach to modeling the effects of Hill-Robertson interference on levels of adaptation and patterns of variability in a nonrecombining genome or genomic region is described. The model assumes a set of L diallelic sites subject to reversible mutations between beneficial and deleterious alleles, with the same selection coefficient at each site. The assumption of reversibility allows the system to reach a stable statistical equilibrium with respect to the frequencies of deleterious mutations, in contrast to many previous models that assume irreversible mutations to deleterious alleles. The model is therefore appropriate for understanding the long-term properties of nonrecombining genomes such as Y chromosomes, and is applicable to haploid genomes or to diploid genomes when there is intermediate dominance with respect to the effects of mutations on fitness. Approximations are derived for the equilibrium frequencies of deleterious mutations, the effective population size that controls the fixation probabilities of mutations at sites under selection, the nucleotide site diversity at neutral sites located within the nonrecombining region, and the site frequency spectrum for segregating neutral variants. The approximations take into account the effects of linkage disequilibrium on the genetic variance at sites under selection. Comparisons with published and new computer simulation results show that the approximations are sufficiently accurate to be useful, and can provide insights into a wider range of parameter sets than is accessible by simulation. The relevance of the findings to data on nonrecombining genome regions is discussed.

Phylogenomics resolves key relationships in Rumex and uncovers a dynamic history of independently evolving sex chromosomes

Sex chromosomes have evolved independently many times across eukaryotes. Despite a considerable body of literature on sex chromosome evolution, the causes and consequences of variation in their formation, degeneration, and turnover remain poorly understood. Chromosomal rearrangements are thought to play an important role in these processes by promoting or extending the suppression of recombination on sex chromosomes. Sex chromosome variation may also contribute to barriers to gene flow, limiting introgression among species. Comparative approaches in groups with sexual system variation can be valuable for understanding these questions. Rumex is a diverse genus of flowering plants harboring significant sexual system and karyotypic variation, including hermaphroditic and dioecious clades with XY (and XYY) sex chromosomes. Previous disagreement in the phylogenetic relationships among key species has rendered the history of sex chromosome evolution uncertain. Resolving this history is important for investigating the interplay of chromosomal rearrangements, introgression, and sex chromosome evolution in the genus. Here, we use new transcriptome assemblies from 11 species representing major clades in the genus, along with a whole-genome assembly generated for a key hermaphroditic species. Using phylogenomic approaches, we find evidence for the independent evolution of sex chromosomes across two major clades, and introgression from unsampled lineages likely predating the formation of sex chromosomes in the genus. Comparative genomic approaches revealed high rates of chromosomal rearrangement, especially in dioecious species, with evidence for a complex origin of the sex chromosomes through multiple chromosomal fusions. However, we found no evidence of elevated rates of fusion on the sex chromosomes in comparison with autosomes, providing no support for an adaptive hypothesis of sex chromosome expansion due to sexually antagonistic selection. Overall, our results highlight a complex history of karyotypic evolution in Rumex, raising questions about the role that chromosomal rearrangements might play in the evolution of large heteromorphic sex chromosomes.

Pan-genome analysis of 13 Spinacia accessions reveals structural variations associated with sex chromosome evolution and domestication traits in spinach

Structural variations (SVs) are major genetic variants that can be involved in the origin, adaptation and domestication of species. However, the identification and characterization of SVs in Spinacia species are rare due to the lack of a pan-genome. Here, we report eight chromosome-scale assemblies of cultivated spinach and its two wild species. After integration with five existing assemblies, we constructed a comprehensive Spinacia pan-genome and identified 193 661 pan-SVs, which were genotyped in 452 Spinacia accessions. Our pan-SVs enabled genome-wide association study identified signals associated with sex and clarified the evolutionary direction of spinach. Most sex-linked SVs (86%) were biased to occur on the Y chromosome during the evolution of the sex-linked region, resulting in reduced Y-linked gene expression. The frequency of pan-SVs among Spinacia accessions further illustrated the contribution of these SVs to domestication, such as bolting time and seed dormancy. Furthermore, compared with SNPs, pan-SVs act as efficient variants in genomic selection (GS) because of their ability to capture missing heritability information and higher prediction accuracy. Overall, this study provides a valuable resource for spinach genomics and highlights the potential utility of pan-SV in crop improvement and breeding programmes.

Allopolyploidization from two dioecious ancestors leads to recurrent evolution of sex chromosomes

Polyploidization presents an unusual challenge for species with sex chromosomes, as it can lead to complex combinations of sex chromosomes that disrupt reproductive development. This is particularly true for allopolyploidization between species with different sex chromosome systems. Here, we assemble haplotype-resolved chromosome-level genomes of a female allotetraploid weeping willow (Salix babylonica) and a male diploid S. dunnii. We show that weeping willow arose from crosses between a female ancestor from the Salix-clade, which has XY sex chromosomes on chromosome 7, and a male ancestor from the Vetrix-clade, which has ancestral XY sex chromosomes on chromosome 15. We find that weeping willow has one pair of sex chromosomes, ZW on chromosome 15, that derived from the ancestral XY sex chromosomes in the male ancestor of the Vetrix-clade. Moreover, the ancestral 7X chromosomes from the female ancestor of the Salix-clade have reverted to autosomal inheritance. Duplicated intact ARR17-like genes on the four homologous chromosomes 19 likely have contributed to the maintenance of dioecy during polyploidization and sex chromosome turnover. Taken together, our results suggest the rapid evolution and reversion of sex chromosomes following allopolyploidization in weeping willow.

The rate of W chromosome degeneration across multiple avian neo-sex chromosomes

When sex chromosomes evolve recombination suppression, the sex-limited chromosome (Y/W) commonly degenerate by losing functional genes. The rate of Y/W degeneration is believed to slow down over time as the most essential genes are maintained by purifying selection, but supporting data are scarce especially for ZW systems. Here, we study W degeneration in Sylvioidea songbirds where multiple autosomal translocations to the sex chromosomes, and multiple recombination suppression events causing separate evolutionary strata, have occurred during the last ~ 28.1-4.5 million years (Myr). We show that the translocated regions have maintained 68.3-97.7% of their original gene content, compared to only 4.2% on the much older ancestral W chromosome. By mapping W gene losses onto a dated phylogeny, we estimate an average gene loss rate of 1.0% per Myr, with only moderate variation between four independent lineages. Consistent with previous studies, evolutionarily constrained and haploinsufficient genes were preferentially maintained on W. However, the gene loss rate did not show any consistent association with strata age or with the number of W genes at strata formation. Our study provides a unique account on the pace of W gene loss and reinforces the significance of purifying selection in maintaining essential genes on sex chromosomes.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

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.

Phased Assembly of Neo-Sex Chromosomes Reveals Extensive Y Degeneration and Rapid Genome Evolution in Rumex hastatulus

Y chromosomes are thought to undergo progressive degeneration due to stepwise loss of recombination and subsequent reduction in selection efficiency. However, the timescales and evolutionary forces driving degeneration remain unclear. To investigate the evolution of sex chromosomes on multiple timescales, we generated a high-quality phased genome assembly of the massive older (<10 MYA) and neo (<200,000 yr) sex chromosomes in the XYY cytotype of the dioecious plant Rumex hastatulus and a hermaphroditic outgroup Rumex salicifolius. Our assemblies, supported by fluorescence in situ hybridization, confirmed that the neo-sex chromosomes were formed by two key events: an X-autosome fusion and a reciprocal translocation between the homologous autosome and the Y chromosome. The enormous sex-linked regions of the X (296 Mb) and two Y chromosomes (503 Mb) both evolved from large repeat-rich genomic regions with low recombination; however, the complete loss of recombination on the Y still led to over 30% gene loss and major rearrangements. In the older sex-linked region, there has been a significant increase in transposable element abundance, even into and near genes. In the neo-sex-linked regions, we observed evidence of extensive rearrangements without gene degeneration and loss. Overall, we inferred significant degeneration during the first 10 million years of Y chromosome evolution but not on very short timescales. Our results indicate that even when sex chromosomes emerge from repetitive regions of already-low recombination, the complete loss of recombination on the Y chromosome still leads to a substantial increase in repetitive element content and gene degeneration.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

Evolution of a ZW sex chromosome system in willows

Transitions in the heterogamety of sex chromosomes (e.g., XY to ZW or vice versa) fundamentally alter the genetic basis of sex determination, however the details of these changes have been studied in only a few cases. In an XY to ZW transition, the X is likely to give rise to the W because they both carry feminizing genes and the X is expected to harbour less genetic load than the Y. Here, using a new reference genome for Salix exigua, we trace the X, Y, Z, and W sex determination regions during the homologous transition from an XY system to a ZW system in willow (Salix). We show that both the W and the Z arose from the Y chromosome. We find that the new Z chromosome shares multiple homologous putative masculinizing factors with the ancestral Y, whereas the new W lost these masculinizing factors and gained feminizing factors. The origination of both the W and Z from the Y was permitted by an unexpectedly low genetic load on the Y and this indicates that the origins of sex chromosomes during homologous transitions may be more flexible than previously considered.

A plastome phylogeny of Rumex (Polygonaceae) illuminates the divergent evolutionary histories of docks and sorrels

The genus Rumex L. (Polygonaceae) provides a unique system for investigating the evolutionary development of sex determination and molecular rate evolution. Historically, Rumex has been divided, both taxonomically and colloquially into two groups: 'docks' and 'sorrels'. A well-resolved phylogeny can help evaluate a genetic basis for this division. Here we present a plastome phylogeny for 34 species of Rumex, inferred using maximum likelihood criteria. The historical 'docks' (Rumex subgenus Rumex) were resolved as monophyletic. The historical 'sorrels' (Rumex subgenera Acetosa and Acetosella) were resolved together, though not monophyletic due to the inclusion of R. bucephalophorus (Rumex subgenus Platypodium). Emex is supported as its own subgenus within Rumex, instead of resolved as sister taxa. We found remarkably low nucleotide diversity among the docks, consistent with recent diversification in that group, especially as compared to the sorrels. Fossil calibration of the phylogeny suggested that the common ancestor for Rumex (including Emex) has origins in the lower Miocene (22.13 MYA). The sorrels appear to have subsequently diversified at a relatively constant rate. The origin of the docks, however, was placed in the upper Miocene, but with most speciation occurring in the Plio-Pleistocene.

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.

A new phylogeny of Rumex (Polygonaceae) adds evolutionary context to the diversity of reproductive systems present in the genus

Rumex is one of about 50 genera in the knotweed family, Polygonaceae. The genus comprises about 200 species with bisexual, or more commonly, unisexual flowers, with the species displaying monoecious, dioecious, synoecious (hermaphroditic) or polygamous reproductive systems. Some of the dioecious species have heteromorphic sex chromosomes, which is rare amongst angiosperms. We here present a plastid phylogeny of 67 species, representing all four subgenera. For this study, we used three chloroplast markers, rbcL, trnH-psbA, trnL-F and dense taxon sampling to reconstruct the most comprehensive molecular phylogeny of Rumex to date. The reconstructed phylogeny for this work resolves six major clades and one large grade in Rumexsubg.Rumex. In addition, the species with known dioecious reproductive systems are resolved within a broader clade we term "the dioecious clade". These results suggest that the species with divergent reproductive systems are more closely related to each other than to other species comprising the rest of the Rumex genus.

Physiological aspects of sex differences and Haldane's rule in Rumex hastatulus

Haldane's rule (HR, impairment of fertility and/or viability of interracial hybrids) seems to be one of few generalizations in evolutionary biology. The validity of HR has been confirmed in animals, and more recently in some dioecious plants (Silene and Rumex). Dioecious Rumex hastatulus has two races differing in the sex chromosome system: Texas (T) and North Carolina (NC), and T × NC males showed both reduced pollen fertility and rarity-two classical symptoms of Haldane's rule (HR). The reduced fertility of these plants has a simple mechanistic explanation, but the reason for their rarity was not elucidated. Here, we measured selected physiological parameters related to the antioxidant defense system in parental races and reciprocal hybrids of R. hastatulus. We showed that the X-autosome configurations, as well as asymmetries associated with Y chromosomes and cytoplasm, could modulate this system in hybrids. The levels and quantitative patterns of the measured parameters distinguish the T × NC hybrid from the other analyzed forms. Our observations suggest that the rarity of T × NC males is caused postzygotically and most likely related to the higher level of oxidative stress induced by the chromosomal incompatibilities. It is the first report on the physiological aspects of HR in plants.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.

Are plant and animal sex chromosomes really all that different?

Sex chromosomes in plants have often been contrasted with those in animals with the goal of identifying key differences that can be used to elucidate fundamental evolutionary properties. For example, the often homomorphic sex chromosomes in plants have been compared to the highly divergent systems in some animal model systems, such as birds, Drosophila and therian mammals, with many hypotheses offered to explain the apparent dissimilarities, including the younger age of plant sex chromosomes, the lesser prevalence of sexual dimorphism, or the greater extent of haploid selection. Furthermore, many plant sex chromosomes lack complete sex chromosome dosage compensation observed in some animals, including therian mammals, Drosophila, some poeciliids, and Anolis, and plant dosage compensation, where it exists, appears to be incomplete. Even the canonical theoretical models of sex chromosome formation differ somewhat between plants and animals. However, the highly divergent sex chromosomes observed in some animal groups are actually the exception, not the norm, and many animal clades are far more similar to plants in their sex chromosome patterns. This begs the question of how different are plant and animal sex chromosomes, and which of the many unique properties of plants would be expected to affect sex chromosome evolution differently than animals? In fact, plant and animal sex chromosomes exhibit more similarities than differences, and it is not at all clear that they differ in terms of sexual conflict, dosage compensation, or even degree of divergence. Overall, the largest difference between these two groups is the greater potential for haploid selection in plants compared to animals. This may act to accelerate the expansion of the non-recombining region at the same time that it maintains gene function within it. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Recombination landscape dimorphism and sex chromosome evolution in the dioecious plant Rumex hastatulus

There is growing evidence from diverse taxa for sex differences in the genomic landscape of recombination, but the causes and consequences of these differences remain poorly understood. Strong recombination landscape dimorphism between the sexes could have important implications for the dynamics of sex chromosome evolution because low recombination in the heterogametic sex can favour the spread of sexually antagonistic alleles. Here, we present a sex-specific linkage map and revised genome assembly of Rumex hastatulus and provide the first evidence and characterization of sex differences in recombination landscape in a dioecious plant. We present data on significant sex differences in recombination, with regions of very low recombination in males covering over half of the genome. This pattern is evident on both sex chromosomes and autosomes, suggesting that pre-existing differences in recombination may have contributed to sex chromosome formation and divergence. Our analysis of segregation distortion suggests that haploid selection due to pollen competition occurs disproportionately in regions with low male recombination. We hypothesize that sex differences in the recombination landscape have contributed to the formation of a large heteromorphic pair of sex chromosomes in R. hastatulus, but more comparative analyses of recombination will be important to investigate this hypothesis further. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Some thoughts about the words we use for thinking about sex chromosome evolution

Sex chromosomes are familiar to most biologists since they first learned about genetics. However, research over the past 100 years has revealed that different organisms have evolved sex-determining systems independently. The differences in the ages of systems, and in how they evolved, both affect whether sex chromosomes have evolved. However, the diversity means that the terminology used tends to emphasize either the similarities or the differences, sometimes causing misunderstandings. In this article, I discuss some concepts where special care is needed with terminology. The following four terms regularly create problems: ‘sex chromosome’, ‘master sex-determining gene’, ‘evolutionary strata’ and ‘genetic degeneration’. There is no generally correct or wrong use of these words, but efforts are necessary to make clear how they are to be understood in specific situations. I briefly outline some widely accepted ideas about sex chromosomes, and then discuss these ‘problem terms’, highlighting some examples where careful use of the words helps bring to light current uncertainties and interesting questions for future work.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Dosage compensation evolution in plants: theories, controversies and mechanisms

In a minority of flowering plants, separate sexes are genetically determined by sex chromosomes. The Y chromosome has a non-recombining region that degenerates, causing a reduced expression of Y genes. In some species, the lower Y expression is accompanied by dosage compensation (DC), a mechanism that re-equalizes male and female expression and/or brings XY male expression back to its ancestral level. Here, we review work on DC in plants, which started as early as the late 1960s with cytological approaches. The use of transcriptomics fired a controversy as to whether DC existed in plants. Further work revealed that various plants exhibit partial DC, including a few species with young and homomorphic sex chromosomes. We are starting to understand the mechanisms responsible for DC in some plants, but in most species, we lack the data to differentiate between global and gene-by-gene DC. Also, it is unknown why some species evolve many dosage compensated genes while others do not. Finally, the forces that drive DC evolution remain mysterious, both in plants and animals. We review the multiple evolutionary theories that have been proposed to explain DC patterns in eukaryotes with XY or ZW sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Recurrent allopolyploidization, Y-chromosome introgression and the evolution of sexual systems in the plant genus Mercurialis

The plant genus Mercurialis includes dioecious, monoecious and androdioecious species (where males coexist with hermaphrodites). Its diversification involved reticulate evolution via hybridization and polyploidization. The Y chromosome of the diploid species Mercurialis annua shows only mild signs of degeneration. We used sequence variation at a Y-linked locus in several species and at multiple autosomal and pseudoautosomal loci to investigate the origin and evolution of the Y chromosome across the genus. Our study provides evidence for further cases of allopolyploid speciation. It also reveals that all lineages with separate sexes (with one possible exception) share the same ancestral Y chromosome. Surprisingly, males in androdioecious populations of hexaploid M. annua carry a Y chromosome that is not derived from either of its two putative progenitor lineages but from a more distantly related perennial dioecious lineage via introgression. These results throw new light on the evolution of sexual systems and polyploidy in Mercurialis and secure it as a promising model for further study of plant sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

The contributions of Nettie Stevens to the field of sex chromosome biology

The early 1900s delivered many foundational discoveries in genetics, including re-discovery of Mendel's research and the chromosomal theory of inheritance. Following these insights, many focused their research on whether the development of separate sexes had a chromosomal basis or if instead it was caused by environmental factors. It is Dr Nettie M. Stevens' Studies in spermatogenesis (1905) that provided the unequivocal evidence that the inheritance of the Y chromosome initiated male development in mealworms. This result established that sex is indeed a Mendelian trait with a genetic basis and that the sex chromosomes play a critical role. In Part II of Studies in spermatogenesis (1906), an XY pair was identified in dozens of additional species, further validating the function of sex chromosomes. Since this formative work, a wealth of studies in animals and plants have examined the genetic basis of sex. The goal of this review is to shine a light again on Stevens’ Studies in spermatogenesis and the lasting impact of this work. We additionally focus on key findings in plant systems over the last century and open questions that are best answered, as in Stevens' work, by synthesizing across many systems.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Evolution of the sex-determining region in Ginkgo biloba

Sex chromosomes or sex-determining regions (SDR) have been discovered in many dioecious plant species, including the iconic 'living fossil' Ginkgo biloba, though the location and size of the SDR in G. biloba remain contradictory. Here we resolve these controversies and analyse the evolution of the SDR in this species. Based on transcriptome sequencing data from four genetic crosses we reconstruct male- and female-specific genetic maps and locate the SDR to the middle of chromosome 2. Integration of the genetic maps with the genome sequence reveals that recombination in and around the SDR is suppressed in a region of about 50 Mb in both males and females. However, occasional recombination does occur except a small, less than 5 Mb long region that does not recombine in males. Based on synonymous divergence between homologous X- and Y-linked genes in this region, we infer that the Ginkgo SDR is fairly old-at least of Cretaceous origin. The analysis of substitution rates and gene expression reveals only slight Y-degeneration. These results are consistent with findings in other dioecious plants with homomorphic sex chromosomes, where the SDR is typically small and evolves in a region with pre-existing reduced recombination, surrounded by long actively recombining pseudoautosomal regions. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Y recombination arrest and degeneration in the absence of sexual dimorphism

Current theory proposes that degenerated sex chromosomes-such as the mammalian Y-evolve through three steps: (i) recombination arrest, linking male-beneficial alleles to the Y chromosome; (ii) Y degeneration, resulting from the inefficacy of selection in the absence of recombination; and (iii) dosage compensation, correcting the resulting low expression of X-linked genes in males. We investigate a model of sex chromosome evolution that incorporates the coevolution of cis and trans regulators of gene expression. We show that the early emergence of dosage compensation favors the maintenance of Y-linked inversions by creating sex-antagonistic regulatory effects. This is followed by degeneration of these nonrecombining inversions caused by regulatory divergence between the X and Y chromosomes. In contrast to current theory, the whole process occurs without any selective pressure related to sexual dimorphism.

Heterogeneous Histories of Recombination Suppression on Stickleback Sex Chromosomes

How consistent are the evolutionary trajectories of sex chromosomes shortly after they form? Insights into the evolution of recombination, differentiation, and degeneration can be provided by comparing closely related species with homologous sex chromosomes. The sex chromosomes of the threespine stickleback (Gasterosteus aculeatus) and its sister species, the Japan Sea stickleback (G. nipponicus), have been well characterized. Little is known, however, about the sex chromosomes of their congener, the blackspotted stickleback (G. wheatlandi). We used pedigrees to obtain experimentally phased whole genome sequences from blackspotted stickleback X and Y chromosomes. Using multispecies gene trees and analysis of shared duplications, we demonstrate that Chromosome 19 is the ancestral sex chromosome and that its oldest stratum evolved in the common ancestor of the genus. After the blackspotted lineage diverged, its sex chromosomes experienced independent and more extensive recombination suppression, greater X-Y differentiation, and a much higher rate of Y degeneration than the other two species. These patterns may result from a smaller effective population size in the blackspotted stickleback. A recent fusion between the ancestral blackspotted stickleback Y chromosome and Chromosome 12, which produced a neo-X and neo-Y, may have been favored by the very small size of the recombining region on the ancestral sex chromosome. We identify six strata on the ancestral and neo-sex chromosomes where recombination between the X and Y ceased at different times. These results confirm that sex chromosomes can evolve large differences within and between species over short evolutionary timescales.

How did the guppy Y chromosome evolve?

The sex chromosome pairs of many species do not undergo genetic recombination, unlike the autosomes. It has been proposed that the suppressed recombination results from natural selection favouring close linkage between sex-determining genes and mutations on this chromosome with advantages in one sex, but disadvantages in the other (these are called sexually antagonistic mutations). No example of such selection leading to suppressed recombination has been described, but populations of the guppy display sexually antagonistic mutations (affecting male coloration), and would be expected to evolve suppressed recombination. In extant close relatives of the guppy, the Y chromosomes have suppressed recombination, and have lost all the genes present on the X (this is called genetic degeneration). However, the guppy Y occasionally recombines with its X, despite carrying sexually antagonistic mutations. We describe evidence that a new Y evolved recently in the guppy, from an X chromosome like that in these relatives, replacing the old, degenerated Y, and explaining why the guppy pair still recombine. The male coloration factors probably arose after the new Y evolved, and have already evolved expression that is confined to males, a different way to avoid the conflict between the sexes.

We report new findings concerning the long-studied the guppy XY pair, which has remained somewhat mystifying. We show that it can be understood as a case of a recent sex chromosome turnover event in which an older, highly degenerated Y chromosome was lost, and creation of a new sex chromosome from the ancestral X. This chromosome acquired a male-determining factor, possibly by a mutation in (or a duplication of) a previously X-linked gene, or (less likely) by movement of an ancestral Y-linked maleness factor onto the X. We relate the findings to theoretical models of such events, and argue that the proposed change was free from factors thought to impede such turnovers. The change resulted in the intriguing situation where the X chromosome is old and the Y is much younger, and we discuss some other species where a similar change seems likely to have occurred.

The timing of genetic degeneration of sex chromosomes

Genetic degeneration is an extraordinary feature of sex chromosomes, with the loss of functions of Y-linked genes in species with XY systems, and W-linked genes in ZW systems, eventually affecting almost all genes. Although degeneration is familiar to most biologists, important aspects are not yet well understood, including how quickly a Y or W chromosome can become completely degenerated. I review the current understanding of the time-course of degeneration. Degeneration starts after crossing over between the sex chromosome pair stops, and theoretical models predict an initially fast degeneration rate and a later much slower one. It has become possible to estimate the two quantities that the models suggest are the most important in determining degeneration rates-the size of the sex-linked region, and the time when recombination became suppressed (which can be estimated using Y-X or W-Z sequence divergence). However, quantifying degeneration is still difficult. I review evidence on gene losses (based on coverage analysis) or loss of function (by classifying coding sequences into functional alleles and pseudogenes). I also review evidence about whether small genome regions degenerate, or only large ones, whether selective constraints on the genes in a sex-linked region also strongly affect degeneration rates, and about how long it takes before all (or almost all) genes are lost. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Schistosome W-linked genes inform temporal dynamics of sex chromosome evolution and suggest candidate for sex determination

Schistosomes, the human parasites responsible for snail fever, are female-heterogametic. Different parts of their ZW sex chromosomes have stopped recombining in distinct lineages, creating "evolutionary strata" of various ages. While the Z-chromosome is well characterized at the genomic and molecular level, the W-chromosome has remained largely unstudied from an evolutionary perspective, as only a few W-linked genes have been detected outside of the model species Schistosoma mansoni. Here, we characterize the gene content and evolution of the W-chromosomes of S. mansoni and of the divergent species S. japonicum. We use a combined RNA/DNA k-mer based pipeline to assemble around one hundred candidate W-specific transcripts in each of the species. About half of them map to known protein coding genes, the majority homologous to S. mansoni Z-linked genes. We perform an extended analysis of the evolutionary strata present in the two species (including characterizing a previously undetected young stratum in S. japonicum) to infer patterns of sequence and expression evolution of W-linked genes at different time points after recombination was lost. W-linked genes show evidence of degeneration, including high rates of protein evolution and reduced expression. Most are found in young lineage-specific strata, with only a few high expression ancestral W-genes remaining, consistent with the progressive erosion of non-recombining regions. Among these, the splicing factor U2AF2 stands out as a promising candidate for primary sex determination, opening new avenues for understanding the molecular basis of the reproductive biology of this group.

Transcriptional regulation of dosage compensation in Carica papaya

Sex chromosome evolution results in the disparity in gene content between heterogametic sex chromosomes and creates the need for dosage compensation to counteract the effects of gene dose imbalance of sex chromosomes in males and females. It is not known at which stage of sex chromosome evolution dosage compensation would evolve. We used global gene expression profiling in male and female papayas to assess gene expression patterns of sex-linked genes on the papaya sex chromosomes. By analyzing expression ratios of sex-linked genes to autosomal genes and sex-linked genes in males relative to females, our results showed that dosage compensation was regulated on a gene-by-gene level rather than whole sex-linked region in papaya. Seven genes on the papaya X chromosome exhibited dosage compensation. We further compared gene expression ratios in the two evolutionary strata. Y alleles in the older evolutionary stratum showed reduced expression compared to X alleles, while Y alleles in the younger evolutionary stratum showed elevated expression compared to X alleles. Reduced expression of Y alleles in the older evolutionary stratum might be caused by accumulation of deleterious mutations in regulatory regions or transposable element-mediated methylation spreading. Most X-hemizygous genes exhibited either no or very low expression, suggesting that gene silencing might play a role in maintaining transcriptional balance between females and males.

When and how do sex-linked regions become sex chromosomes?

The attention given to heteromorphism and genetic degeneration of "classical sex chromosomes" (Y chromosomes in XY systems, and the W in ZW systems that were studied first and are best described) has perhaps created the impression that the absence of recombination between sex chromosomes is inevitable. I here argue that continued recombination is often to be expected, that absence of recombination is surprising and demands further study, and that the involvement of selection in reduced recombination is not yet well understood. Despite a long history of investigations of sex chromosome pairs, there is a need for more quantitative approaches to studying sex-linked regions. I describe a scheme to help understand the relationships between different properties of sex-linked regions. Specifically, I focus on their sizes (differentiating between small regions and extensive fully sex-linked ones), the times when they evolved, and their differentiation, and review studies using DNA sequencing in nonmodel organisms that are providing information about the processes causing these properties.

Widespread Recombination Suppression Facilitates Plant Sex Chromosome Evolution

Classical models suggest that recombination rates on sex chromosomes evolve in a stepwise manner to localize sexually antagonistic variants in the sex in which they are beneficial, thereby lowering rates of recombination between X and Y chromosomes. However, it is also possible that sex chromosome formation occurs in regions with preexisting recombination suppression. To evaluate these possibilities, we constructed linkage maps and a chromosome-scale genome assembly for the dioecious plant Rumex hastatulus. This species has a polymorphic karyotype with a young neo-sex chromosome, resulting from a Robertsonian fusion between the X chromosome and an autosome, in part of its geographic range. We identified the shared and neo-sex chromosomes using comparative genetic maps of the two cytotypes. We found that sex-linked regions of both the ancestral and the neo-sex chromosomes are embedded in large regions of low recombination. Furthermore, our comparison of the recombination landscape of the neo-sex chromosome to its autosomal homolog indicates that low recombination rates mainly preceded sex linkage. These patterns are not unique to the sex chromosomes; all chromosomes were characterized by massive regions of suppressed recombination spanning most of each chromosome. This represents an extreme case of the periphery-biased recombination seen in other systems with large chromosomes. Across all chromosomes, gene and repetitive sequence density correlated with recombination rate, with patterns of variation differing by repetitive element type. Our findings suggest that ancestrally low rates of recombination may facilitate the formation and subsequent evolution of heteromorphic sex chromosomes.

Male-specific Y-chromosomal regions in waterhemp (Amaranthus tuberculatus) and Palmer amaranth (Amaranthus palmeri)

Amaranthus tuberculatus and Amaranthus palmeri are agronomically important weed species, both with stable dioecious reproductive systems. An understanding of the genetic basis of sex determination may lead to new methods of managing these troublesome weeds. Previous research identified genomic sequences associated with maleness in each species. Male-specific sequences were used to identify genomic regions in both species that are believed to contain sex-determining genes, i.e. the male-specific Y (MSY) region. These regions were compared to understand if sex determination is controlled via the same physiological pathway and if dioecy evolved independently. A contiguously assembled candidate MSY region identified in Amaranthus palmeri is approximately 1.3 Mb with 121 predicted gene models. In Amaranthus tuberculatus, several contigs, with combined length of 4.6 Mb and with 147 gene models, were identified as belonging to the MSY region. Synteny was not detected between the two species' candidate MSY regions but they shared two predicted genes. With lists of candidate genes for sex determination containing fewer than 200 in each species, future research can address whether sex determination is controlled via similar physiological pathways and whether dioecy has indeed evolved independently in these species.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa

BACKGROUND AND AIMS

Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes.

METHODS

We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons.

KEY RESULTS

We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes.

CONCLUSIONS

The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Evolutionary Genomics of Plant Gametophytic Selection

It has long been recognized that natural selection during the haploid gametophytic phase of the plant life cycle may have widespread importance for rates of evolution and the maintenance of genetic variation. Recent theoretical advances have further highlighted the significance of gametophytic selection for diverse evolutionary processes. Genomic approaches offer exciting opportunities to address key questions about the extent and effects of gametophytic selection on plant evolution and adaptation. Here, we review the progress and prospects for integrating functional and evolutionary genomics to test theoretical predictions, and to examine the importance of gametophytic selection on genetic diversity and rates of evolution. There is growing evidence that selection during the gametophyte phase of the plant life cycle has important effects on both gene and genome evolution and is likely to have important pleiotropic effects on the sporophyte. We discuss the opportunities to integrate comparative population genomics, genome-wide association studies, and experimental approaches to further distinguish how differential selection in the two phases of the plant life cycle contributes to genetic diversity and adaptive evolution.

Evidence for Dosage Compensation in Coccinia grandis, a Plant with a Highly Heteromorphic XY System

About 15,000 angiosperms are dioecious, but the mechanisms of sex determination in plants remain poorly understood. In particular, how Y chromosomes evolve and degenerate, and whether dosage compensation evolves as a response, are matters of debate. Here, we focus on Coccinia grandis, a dioecious cucurbit with the highest level of X/Y heteromorphy recorded so far. We identified sex-linked genes using RNA sequences from a cross and a model-based method termed SEX-DETector. Parents and F1 individuals were genotyped, and the transmission patterns of SNPs were then analyzed. In the >1300 sex-linked genes studied, maximum X-Y divergence was 0.13-0.17, and substantial Y degeneration is implied by an average Y/X expression ratio of 0.63 and an inferred gene loss on the Y of ~40%. We also found reduced Y gene expression being compensated by elevated expression of corresponding genes on the X and an excess of sex-biased genes on the sex chromosomes. Molecular evolution of sex-linked genes in C. grandis is thus comparable to that in Silene latifolia, another dioecious plant with a strongly heteromorphic XY system, and cucurbits are the fourth plant family in which dosage compensation is described, suggesting it might be common in plants.

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.

Sex Chromosome Degeneration by Regulatory Evolution

In many species, the Y (or W) sex chromosome is degenerate. Current theory proposes that this degeneration follows the arrest of recombination and results from the accumulation of deleterious mutations due to selective interference-the inefficacy of natural selection on non-recombining genomic regions. This theory requires very few assumptions, but it does not robustly predict fast erosion of the Y (or W) in large populations or the stepwise degeneration of several small non-recombining strata. We propose a new mechanism for Y/W erosion that works over faster timescales, in large populations, and for small non-recombining regions (down to a single sex-linked gene). The mechanism is based on the instability and divergence of cis-regulatory sequences in non-recombining genome regions, which become selectively haploidized to mask deleterious mutations on coding sequences. This haploidization is asymmetric, because cis-regulators on the X cannot be silenced (otherwise there would be no expression in females). This process causes rapid Y/W degeneration and simultaneous evolution of dosage compensation, provided that autosomal trans-regulatory sequences with sex-limited effects are available to compensate for cis-regulatory divergence. Although this "degeneration by regulatory evolution" does not require selective interference, both processes may act in concert to further accelerate Y degeneration.

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 Y Chromosome as a Battleground for Intragenomic Conflict

Y chromosomes are typically viewed as genetic wastelands with few intact genes. Recent genomic analyses in Drosophila, however, show that gene gain is prominent on young Y chromosomes. Meiosis- and RNAi-related genes often coamplify on recently formed X and Y chromosomes, are testis-expressed, and produce antisense transcripts and short RNAs. RNAi pathways are also involved in suppressing sex ratio drive in Drosophila. These observations paint a dynamic picture of sex chromosome differentiation, suggesting that rapidly evolving genomic battles over segregation are rampant on young sex chromosomes and utilize RNAi to defend the genome against selfish elements that manipulate fair meiosis. Recurrent sex chromosome drive can have profound ecological, evolutionary, and cellular impacts and account for unique features of sex chromosomes.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabissativa–derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabis sativa-derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

Ancestral and neo-sex chromosomes contribute to population divergence in a dioecious plant

Empirical evidence from several animal groups suggests sex chromosomes disproportionately contribute to reproductive isolation. This effect may be enhanced when sex chromosomes are associated with turnover of sex determination systems resulting from structural rearrangements to the chromosomes. We investigated these predictions in the dioecious plant Rumex hastatulus, which is composed of populations of two different sex chromosome cytotypes caused by an X-autosome fusion. Using population genomic analyses, we investigated the demographic history of R. hastatulus and explored the contributions of ancestral and neo-sex chromosomes to population genetic divergence. Our study revealed that the cytotypes represent genetically divergent populations with evidence for historical but not contemporary gene flow between them. In agreement with classical predictions, we found that the ancestral X chromosome was disproportionately divergent compared with the rest of the genome. Excess differentiation was also observed on the Y chromosome, even when we used measures of differentiation that control for differences in effective population size. Our estimates of the timing of the origin of neo-sex chromosomes in R. hastatulus are coincident with cessation of gene flow, suggesting that the chromosomal fusion event that gave rise to the origin of the XYY cytotype may have also contributed to reproductive isolation.