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

The fourspine stickleback (Apeltes quadracus) has an XY sex chromosome system with polymorphic inversions on both X and Y chromosomes

Teleost fish are well-known for possessing a diversity of sex chromosomes and for undergoing frequent turnovers of these sex chromosomes. However, previous studies have mainly focused on variation between species, while comparatively little attention has been given to sex chromosome polymorphisms within species, which may capture early stages of sex chromosome changes. To better understand the evolution of sex chromosomes, we used the fourspine stickleback (Apeltes quadracus) as a model organism. Previous cytogenetic studies suggested that females of this species possessed a ZW heteromorphic sex chromosome system. However, genetic crosses and our whole-genome sequencing of three geographically distinct wild populations revealed that A. quadracus has an XY sex chromosome on chromosome 23. This chromosome has not previously been identified as a sex chromosome in any other stickleback species, indicating a recent sex chromosome turnover. We also identified two genes - rxfp2a and zar1l - as novel candidate sex determination genes. Notably, we observed inversions on both the X and Y chromosomes in different populations, resulting in distinctive patterns of differentiation between the X and Y chromosomes across populations. The new sex chromosome and intraspecies inversion polymorphisms observed in A. quadracus provide an excellent system for future work assessing the relative fitness effects of the inversions, which will enable testing theoretical models about the drivers of sex chromosome evolution and turnover.

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.

Sex Chromosome Evolution: Hallmarks and Question Marks

Sex chromosomes are widespread in species with separate sexes. They have evolved many times independently and display a truly remarkable diversity. New sequencing technologies and methodological developments have allowed the field of molecular evolution to explore this diversity in a large number of model and nonmodel organisms, broadening our vision on the mechanisms involved in their evolution. Diverse studies have allowed us to better capture the common evolutionary routes that shape sex chromosomes; however, we still mostly fail to explain why sex chromosomes are so diverse. We review over half a century of theoretical and empirical work on sex chromosome evolution and highlight pending questions on their origins, turnovers, rearrangements, degeneration, dosage compensation, gene content, and rates of evolution. We also report recent theoretical progress on our understanding of the ultimate reasons for sex chromosomes' existence.

Sexy ways: approaches to studying plant sex chromosomes

Sex chromosomes have evolved in many plant species with separate sexes. Current plant research is shifting from examining the structure of sex chromosomes to exploring their functional aspects. New studies are progressively unveiling the specific genetic and epigenetic mechanisms responsible for shaping distinct sexes in plants. While the fundamental methods of molecular biology and genomics are generally employed for the analysis of sex chromosomes, it is often necessary to modify classical procedures not only to simplify and expedite analyses but sometimes to make them possible at all. In this review, we demonstrate how, at the level of structural and functional genetics, cytogenetics, and bioinformatics, it is essential to adapt established procedures for sex chromosome analysis.

Evidence for feminized genetic males in a flea beetle using newly identified X-linked markers

The equilibrium of sex ratios in sexually reproducing species is often disrupted by various environmental and genetic factors, including endosymbionts like Wolbachia. In this study, we explore the highly female-biased sex ratio observed in the flea beetle, Altica lythri, and its underlying mechanisms. Ancient hybridization events between Altica species have led to mitochondrial DNA introgression, resulting in distinct mitochondrial haplotypes that go along with different Wolbachia infections (HT1-wLytA1, HT1*- uninfected, HT2-wLytA2, and HT3-wLytB). Notably, beetles with some haplotypes exclusively produce female offspring, suggesting potential Wolbachia-induced phenomena such as feminization of genetic males. However, the observed female bias could also be a consequence of the ancient hybridization resulting in nuclear-cytoplasmic conflicts between introgressed mtDNA and nuclear genes. Through transcriptomic analysis and the program SEX-DETector, we established markers for genotypic sex differentiation for A. lythri, enabling genetic sexing via qPCR. Our findings suggest that feminization of genetic males is contributing to the skewed sex ratios, highlighting the intricate dynamics of sex determination and reproductive strategies in this flea beetle. This study provides valuable insights into the dynamics of genetic conflicts, endosymbionts, and sex ratios, revealing the novel phenomenon of genetic male feminization in the flea beetle A. lythri.

Using a Handful of Transcriptomes to Detect Sex-Linked Markers and Develop Molecular Sexing Assays in a Species with Homomorphic Sex Chromosomes

To understand the biology of a species, it is often crucial to be able to differentiate males and females. However, many species lack easily identifiable sexually dimorphic traits. In those that possess sex chromosomes, molecular sexing offers a good alternative, and molecular sexing assays can be developed through the comparison of male and female genomic sequences. However, in many nonmodel species, sex chromosomes are poorly differentiated, and identifying sex-linked sequences and developing sexing assays can be challenging. In this study, we highlight a simple transcriptome-based procedure for the detection of sex-linked markers suitable for the development of sexing assays that circumvents limitations of more commonly used approaches. We apply it to the spotted snow skink Carinascincus ocellatus, a viviparous lizard with homomorphic XY chromosomes that has environmentally induced sex reversal. With transcriptomes from three males and three females alone, we identify thousands of putative Y-linked sequences. We confirm linkage through alignment of assembled transcripts to a distantly related lizard genome and readily design multiple single locus polymerase chain reaction primers to sex C. ocellatus and related species. Our approach also facilitates valuable comparisons of sex determining systems on a broad taxonomic scale.

An effective strategy for assembling the sex-limited chromosome

BACKGROUND

Most currently available reference genomes lack the sequence map of sex-limited (such as Y and W) chromosomes, which results in incomplete assemblies that hinder further research on sex chromosomes. Recent advancements in long-read sequencing and population sequencing have provided the opportunity to assemble sex-limited chromosomes without the traditional complicated experimental efforts.

FINDINGS

We introduce the first computational method, Sorting long Reads of Y or other sex-limited chromosome (SRY), which achieves improved assembly results compared to flow sorting. Specifically, SRY outperforms in the heterochromatic region and demonstrates comparable performance in other regions. Furthermore, SRY enhances the capabilities of the hybrid assembly software, resulting in improved continuity and accuracy.

CONCLUSIONS

Our method enables true complete genome assembly and facilitates downstream research of sex-limited chromosomes.

Can a Y Chromosome Degenerate in an Evolutionary Instant? A Commentary on Fong et al. 2023

It is well known that the Y chromosomes of Drosophila and mammals and the W chromosomes of birds carry only small fractions of the genes carried by the homologous X or Z chromosomes, and this "genetic degeneration" is associated with loss of recombination between the sex chromosome pair. However, it is still not known how much evolutionary time is needed to reach such nearly complete degeneration. The XY pair of species in a group of closely related poecilid fish is homologous but has been found to have either nondegenerated or completely degenerated Y chromosomes. We evaluate evidence described in a recent paper and show that the available data cast doubt on the view that degeneration has been extraordinarily rapid in the latter (Micropoecilia species).

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.

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

The evolution of huge Y chromosomes in Coccinia grandis and its sister, Coccinia schimperi

Microscopically dimorphic sex chromosomes in plants are rare, reducing our ability to study them. One difficulty has been the paucity of cultivatable species pairs for cytogenetic, genomic and experimental work. Here, we study the newly recognized sisters Coccinia grandis and Coccinia schimperi, both with large Y chromosomes as we here show for Co. schimperi. We built genetic maps for male and female Co. grandis using a full-sibling family, inferred gene sex-linkage, and, with Co. schimperi transcriptome data, tested whether X- and Y-alleles group by species or by sex. Most sex-linked genes for which we could include outgroups grouped the X- and Y-alleles by species, but some 10% instead grouped the two species' X-alleles. There was no relationship between XY synonymous-site divergences in these genes and gene position on the non-recombining part of the X, suggesting recombination arrest shortly before or after species divergence, here dated to about 3.6 Ma. Coccinia grandis and Co. schimperi are the species pair with the most heteromorphic sex chromosomes in vascular plants (the condition in their sister remains unknown), and future work could use them to study mechanisms of Y chromosome enlargement and parallel degeneration, or to test Haldane's rule about lower hybrid fitness in the heterogametic sex. 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'.

Labile sex expression in angiosperm species with sex chromosomes

Here, we review the literature on sexual lability in dioecious angiosperm species with well-studied sex chromosomes. We distinguish three types of departures from strict dioecy, concerning either a minority of flowers in some individuals (leakiness) or the entire individual, which can constantly be bisexual or change sex. We found that for only four of the 22 species studied, reports of lability are lacking. The occurrence of lability is only weakly related to sex chromosome characteristics (number of sex-linked genes, age of the non-recombining region). These results contradict the naive idea that lability is an indication of the absence or the recent evolution of sex chromosomes, and thereby contribute to a growing consensus that sex chromosomes do not necessarily fix sex determination once and for all. We discuss some implications of these findings for the evolution of sex chromosomes, and suggest that more species with well-characterized lability should be studied with genomic data and tools. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Findzx: an automated pipeline for detecting and visualising sex chromosomes using whole-genome sequencing data

BACKGROUND

Sex chromosomes have evolved numerous times, as revealed by recent genomic studies. However, large gaps in our knowledge of sex chromosome diversity across the tree of life remain. Filling these gaps, through the study of novel species, is crucial for improved understanding of why and how sex chromosomes evolve. Characterization of sex chromosomes in already well-studied organisms is also important to avoid misinterpretations of population genomic patterns caused by undetected sex chromosome variation.

RESULTS

Here we present findZX, an automated Snakemake-based computational pipeline for detecting and visualizing sex chromosomes through differences in genome coverage and heterozygosity between any number of males and females. A main feature of the pipeline is the option to perform a genome coordinate liftover to a reference genome of another species. This allows users to inspect sex-linked regions over larger contiguous chromosome regions, while also providing important between-species synteny information. To demonstrate its effectiveness, we applied findZX to publicly available genomic data from species belonging to widely different taxonomic groups (mammals, birds, reptiles, and fish), with sex chromosome systems of different ages, sizes, and levels of differentiation. We also demonstrate that the liftover method is robust over large phylogenetic distances (> 80 million years of evolution).

CONCLUSIONS

With findZX we provide an easy-to-use and highly effective tool for identification of sex chromosomes. The pipeline is compatible with both Linux and MacOS systems, and scalable to suit different computational platforms.

Evolution of the canonical sex chromosomes of the guppy and its relatives

The sex chromosomes of the guppy, Poecilia reticulata, and its close relatives are of particular interest: they are much younger than the highly degenerate sex chromosomes of model systems such as humans and Drosophila melanogaster, and they carry many of the genes responsible for the males' dramatic coloration. Over the last decade, several studies have analyzed these sex chromosomes using a variety of approaches including sequencing genomes and transcriptomes, cytology, and linkage mapping. Conflicting conclusions have emerged, in particular concerning the history of the sex chromosomes and the evolution of suppressed recombination between the X and Y. Here, we address these controversies by reviewing the evidence and reanalyzing data. We find no evidence of a nonrecombining sex-determining region or evolutionary strata in P. reticulata. Furthermore, we find that the data most strongly support the hypothesis that the sex-determining regions of 2 close relatives of the guppy, Poecilia wingei and Micropoecilia picta, evolved independently after their lineages diverged. We identify possible causes of conflicting results in previous studies and suggest best practices going forward.

A derived ZW chromosome system in Amborella trichopoda, representing the sister lineage to all other extant flowering plants

The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

Characterization of a Sex-Determining Region and Its Genomic Context via Statistical Estimates of Haplotype Frequencies in Daughters and Sons Sequenced in Pools

Sex chromosomes are generally derived from a pair of autosomes that have acquired a locus controlling sex. Sex chromosomes may evolve reduced recombination around this locus and undergo a long process of molecular divergence. At that point, the original loci controlling sex may be difficult to pinpoint. This difficulty has affected many model species from mammals to birds to flies, which present highly diverged sex chromosomes. Identifying sex-controlling loci is easier in species with molecularly similar sex chromosomes. Here we aimed at pinpointing the sex-determining region (SDR) of Armadillidium vulgare, a terrestrial isopod with female heterogamety (ZW females and ZZ males) and whose sex chromosomes appear to show low genetic divergence. To locate the SDR, we assessed single-nucleotide polymorphism (SNP) allele frequencies in F1 daughters and sons sequenced in pools (pool-seq) in several families. We developed a Bayesian method that uses the SNP genotypes of individually sequenced parents and pool-seq data from F1 siblings to estimate the genetic distance between a given genomic region (contig) and the SDR. This allowed us to assign more than 43 Mb of contigs to sex chromosomes, and to demonstrate extensive recombination and very low divergence between these chromosomes. By taking advantage of multiple F1 families, we delineated a very short genomic region (∼65 kb) that presented no evidence of recombination with the SDR. In this short genomic region, the comparison of sequencing depths between sexes highlighted female-specific genes that have undergone recent duplication, and which may be involved in sex determination in A. vulgare.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F₁ family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

Male Differentiation in the Marine Copepod Oithona nana Reveals the Development of a New Nervous Ganglion and Lin12-Notch-Repeat Protein-Associated Proteolysis

Simple Summary

Copepods are tiny crustaceans and the most abundant animals on Earth; they also play a crucial role in the marine food chain. Among copepods, Oithona is one of most ecologically successful and is known for its differential behavior between males and females. The males adopt the strategy “live fast, die young”: they are constantly in motion to search for females, more vulnerable to predators, feed less, and have a higher mortality rate. In our study, we found the presence of a new male-specific ganglion in Oithona nana probably involved in female cues sensing. We also demonstrate the potential role of new Lin-12 Notch Repeat proteins in the development of this ganglion by interacting with proteins involved in the development of the nervous system. Thanks to our findings, we propose that the “live fast, die young” strategy of the O. nana males is optimized by the explosion of these Lin-12 Notch Repeat proteins in the male proteome involved in the development of the male-specific olfactory ganglion to increase female cue sensing and mating.

Abstract

Copepods are among the most numerous animals, and they play an essential role in the marine trophic web and biogeochemical cycles. The genus Oithona is described as having the highest density of copepods. The Oithona male paradox describes the activity states of males, which are obliged to alternate between immobile and mobile phases for ambush feeding and mate searching, respectively, while the female is less mobile and feeds less. To characterize the molecular basis of this sexual dimorphism, we combined immunofluorescence, genomics, transcriptomics, and protein–protein interaction approaches and revealed the presence of a male-specific nervous ganglion. Transcriptomic analysis showed male-specific enrichment for nervous system development-related transcripts. Twenty-seven Lin12-Notch Repeat domain-containing protein coding genes (LDPGs) of the 75 LDPGs identified in the genome were specifically expressed in males. Furthermore, some LDPGs coded for proteins with predicted proteolytic activity, and proteases-associated transcripts showed a male-specific enrichment. Using yeast double–hybrid assays, we constructed a protein–protein interaction network involving two LDPs with proteases, extracellular matrix proteins, and neurogenesis-related proteins. We also hypothesized possible roles of the LDPGs in the development of the lateral ganglia through helping in extracellular matrix lysis, neurites growth guidance, and synapses genesis.

Detecting sex-linked genes using genotyped individuals sampled in natural populations

We propose a method, SDpop, able to infer sex-linkage caused by recombination suppression typical of sex chromosomes. The method is based on the modeling of the allele and genotype frequencies of individuals of known sex in natural populations. It is implemented in a hierarchical probabilistic framework, accounting for different sources of error. It allows statistical testing for the presence or absence of sex chromosomes, and detection of sex-linked genes based on the posterior probabilities in the model. Furthermore, for gametologous sequences, the haplotype and level of nucleotide polymorphism of each copy can be inferred, as well as the divergence between them. We test the method using simulated data, as well as data from both a relatively recent and an old sex chromosome system (the plant Silene latifolia and humans) and show that, for most cases, robust predictions are obtained with 5 to 10 individuals per sex.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility

Whole‐genome sequencing of non‐model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data‐handling. Three such regions are the major histocompatibility complex (MHC), sex‐determining regions (SDRs) and the plant self‐incompatibility locus (S‐locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short‐read sequences against multiple alternative references improves genotyping comprehensiveness at the S‐locus thereby contributing to more accurate assessments of allelic frequencies. Long‐read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S‐locus haplotypes. These advances are adding value to short‐read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

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

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

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication

Background

A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized.

Results

We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined.

Conclusions

We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

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.

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.

Cytonuclear Genetic Incompatibilities in Plant Speciation

Due to the endosymbiotic origin of organelles, a pattern of coevolution and coadaptation between organellar and nuclear genomes is required for proper cell function. In this review, we focus on the impact of cytonuclear interaction on the reproductive isolation of plant species. We give examples of cases where species exhibit barriers to reproduction which involve plastid-nuclear or mito-nuclear genetic incompatibilities, and describe the evolutionary processes at play. We also discuss potential mechanisms of hybrid fitness recovery such as paternal leakage. Finally, we point out the possible interplay between plant mating systems and cytonuclear coevolution, and its consequence on plant speciation.

Sex differences in spiders: from phenotype to genomics

Sexual reproduction is pervasive in animals and has led to the evolution of sexual dimorphism. In most animals, males and females show marked differences in primary and secondary sexual traits. The formation of sex-specific organs and eventually sex-specific behaviors is defined during the development of an organism. Sex determination processes have been extensively studied in a few well-established model organisms. While some key molecular regulators are conserved across animals, the initiation of sex determination is highly diverse. To reveal the mechanisms underlying the development of sexual dimorphism and to identify the evolutionary forces driving the evolution of different sexes, sex determination mechanisms must thus be studied in detail in many different animal species beyond the typical model systems. In this perspective article, we argue that spiders represent an excellent group of animals in which to study sex determination mechanisms. We show that spiders are sexually dimorphic in various morphological, behavioral, and life history traits. The availability of an increasing number of genomic and transcriptomic resources and functional tools provides a great starting point to scrutinize the extensive sexual dimorphism present in spiders on a mechanistic level. We provide an overview of the current knowledge of sex determination in spiders and propose approaches to reveal the molecular and genetic underpinnings of sexual dimorphism in these exciting animals.

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.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data

Background

Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference.

Results

Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3).

Conclusions

Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

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

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

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

Transitions in sex determination and sex chromosomes across vertebrate species

Despite the prevalence of sexual reproduction across eukaryotes, there is a remarkable diversity of sex-determination mechanisms. The underlying causes of this diversity remain unclear, and it is unknown whether there are convergent trends in the directionality of turnover in sex-determination mechanisms. We used the recently assembled Tree of Sex database to assess patterns in the evolution of sex-determination systems in the remarkably diverse vertebrate clades of teleost fish, squamate reptiles and amphibians. Contrary to theoretical predictions, we find no evidence that the evolution of separate sexes is irreversible, as transitions from separate sexes to hermaphroditism occur at higher rates than the reverse in fish. We also find that transitions from environmental sex determination to genetic sex determination occur at higher rates than the reverse in both squamates and fish, suggesting that genetic sex determination is more stable. However, our data are not consistent with the hypothesis that heteromorphic sex chromosomes are an "evolutionary trap." Rather, we find similar transition rates between homomorphic and heteromorphic sex chromosomes in both fish and amphibians, and to environmental sex determination from heteromorphic vs. homomorphic sex chromosome systems in fish. Finally, we find that transitions between male and female heterogamety occur at similar rates in amphibians and squamates, while transitions to male heterogamety occur at higher rates in fish. Together, these results provide the most comprehensive view to date of the evolution of vertebrate sex determination in a phylogenetic context, providing new insight into long-standing questions about the evolution of sexual reproduction.

Genomic imprinting mediates dosage compensation in a young plant XY system

Sex chromosomes have repeatedly evolved from a pair of autosomes. Consequently, X and Y chromosomes initially have similar gene content, but ongoing Y degeneration leads to reduced expression and eventual loss of Y genes¹. The resulting imbalance in gene expression between Y genes and the rest of the genome is expected to reduce male fitness, especially when protein networks have components from both autosomes and sex chromosomes. A diverse set of dosage compensating mechanisms that alleviates these negative effects has been described in animals^2-4. However, the early steps in the evolution of dosage compensation remain unknown, and dosage compensation is poorly understood in plants⁵. Here, we describe a dosage compensation mechanism in the evolutionarily young XY sex determination system of the plant Silene latifolia. Genomic imprinting results in higher expression from the maternal X chromosome in both males and females. This compensates for reduced Y expression in males, but results in X overexpression in females and may be detrimental. It could represent a transient early stage in the evolution of dosage compensation. Our finding has striking resemblance to the first stage proposed by Ohno⁶ for the evolution of X inactivation in mammals.

Has adaptation occurred in males and females since separate sexes evolved in the plant Silene latifolia?

The evolution of separate sexes may involve changed expression of many genes, as each sex adapts to its new state. Evidence is accumulating for sex differences in expression even in organisms that have recently evolved separate sexes from hermaphrodite or monoecious (cosexual) ancestors, such as some dioecious flowering plants. We describe evidence that a dioecious plant species with recently evolved dioecy, Silene latifolia, has undergone adaptive changes that improve functioning in females, in addition to changes that are probably pleiotropic effects of male sterility. The results suggest pervasive adaptations as soon as males and females evolve from their cosexual ancestor.

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.

Evolution of Sex Chromosome Dosage Compensation in Animals: A Beautiful Theory, Undermined by Facts and Bedeviled by Details

Many animals with genetic sex determination harbor heteromorphic sex chromosomes, where the heterogametic sex has half the gene dose of the homogametic sex. This imbalance, if reflected in the abundance of transcripts or proteins, has the potential to deleteriously disrupt interactions between X-linked and autosomal loci in the heterogametic sex. Classical theory predicts that molecular mechanisms will evolve to provide dosage compensation that recovers expression levels comparable to ancestral expression prior to sex chromosome divergence. Such dosage compensating mechanisms may also, secondarily, result in balanced sex-linked gene expression between males and females. However, numerous recent studies addressing sex chromosome dosage compensation (SCDC) in a diversity of animals have yielded a surprising array of patterns concerning dosage compensation in the heterogametic sex, as well as dosage balance between sexes. These results substantially contradict longstanding theory, catalyzing both novel perspectives and new approaches in dosage compensation research. In this review, we summarize the theory, analytical approaches, and recent results concerning evolutionary patterns of SCDC in animals. We also discuss methodological challenges and discrepancies encountered in this research, which often underlie conflicting results. Finally, we discuss what outstanding questions and opportunities exist for future research on SCDC.

De novo transcriptome assembly from flower buds of dioecious, gynomonoecious and chemically masculinized female Coccinia grandis reveals genes associated with sex expression and modification

BACKGROUND

Coccinia grandis (ivy gourd), is a dioecious member of Cucurbitaceae having heteromorphic sex chromosomes. Chromosome constitution of male and female plants of C. grandis is 22A + XY and 22A + XX respectively. Earlier we showed that a unique gynomonoecious form of C. grandis (22A + XX) also exists in nature bearing morphologically hermaphrodite flowers (GyM-H). Additionally, application of silver nitrate (AgNO₃) on female plants induces stamen development leading to the formation of morphologically hermaphrodite flowers (Ag-H) despite the absence of Y-chromosome. Due to the unavailability of genome sequence and the slow pace at which sex-linked genes are identified, sex expression and modification in C. grandis are not well understood.

RESULTS

We have carried out a comprehensive RNA-Seq study from early-staged male, female, GyM-H, and Ag-H as well as middle-staged male and GyM-H flower buds. A de novo transcriptome was assembled using Trinity and annotated by BLAST2GO and Trinotate pipelines. The assembled transcriptome consisted of 467,233 'Trinity Transcripts' clustering into 378,860 'Trinity Genes'. Female_Early_vs_Male_Early, Ag_Early_vs_Female_Early, and GyM-H_Middle_vs_Male_Middle comparisons exhibited 35,694, 3574, and 14,954 differentially expressed transcripts respectively. Further, qRT-PCR analysis of selected candidate genes validated digital gene expression profiling results. Interestingly, ethylene response-related genes were found to be upregulated in female buds compared to male buds. Also, we observed that AgNO₃ treatment suppressed ethylene responses in Ag-H flowers by downregulation of ethylene-responsive transcription factors leading to stamen development. Further, GO terms related to stamen development were enriched in early-staged male, GyM-H, and Ag-H buds compared to female buds supporting the fact that stamen growth gets arrested in female flowers.

CONCLUSIONS

Suppression of ethylene responses in both male and Ag-H compared to female buds suggests a probable role of ethylene in stamen suppression similar to monoecious cucurbits such as melon and cucumber. Also, pollen fertility associated GO terms were depleted in middle-staged GyM-H buds compared to male buds indicating the necessity of Y-chromosome for pollen fertility. Overall, this study would enable identification of new sex-biased genes for further investigation of stamen arrest, pollen fertility, and AgNO₃-mediated sex modification.

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.

Evolution of sex-biased gene expression in a dioecious plant

Separate sexes and sex-biased gene expression have repeatedly evolved in animals and plants, but the underlying changes in gene expression remain unknown. Here, we studied a pair of plant species, one in which separate sexes and sex chromosomes evolved recently and one which maintained hermaphrodite flowers resembling the ancestral state, to reconstruct expression changes associated with the evolution of dioecy. We found that sex-biased gene expression has evolved in autosomal and sex-linked genes in the dioecious species. Most expression changes relative to hermaphrodite flowers occurred in females rather than males, with higher and lower expression in females leading to female-biased and male-biased expression, respectively. Expression changes were more common in genes located on the sex chromosomes than the autosomes and led to feminization of the X chromosome and masculinization of the Y chromosome. Our results support a scenario in which sex-biased gene expression evolved during the evolution of dioecy to resolve intralocus sexual conflicts over the allocation of resources.