Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution

Out with the old, in with the new: Meiotic driving of sex chromosome evolution

Chromosomal regions with meiotic drivers exhibit biased transmission (> 50 %) over their competing homologous chromosomal region. These regions often have two prominent genetic features: suppressed meiotic crossing over and rapidly evolving multicopy gene families. Heteromorphic sex chromosomes (e.g., XY) often share these two genetic features with chromosomal regions exhibiting meiotic drive. Here, we discuss parallels between meiotic drive and sex chromosome evolution, how the divergence of heteromorphic sex chromosomes can be influenced by meiotic drive, experimental approaches to study meiotic drive on sex chromosomes, and meiotic drive in traditional and non-traditional model organisms with high-quality genome assemblies. The newly available diversity of high-quality sex chromosome sequences allows us to revisit conventional models of sex chromosome evolution through the lens of meiotic drive.

Transforming growth factor-β (TGF-β): A master signal pathway in teleost sex determination

Sex determination and differentiation in fish has always been a hot topic in genetic breeding of aquatic animals. With the advances in next-generation sequencing (NGS) in recent years, sex chromosomes and sex determining genes can be efficiently identified in teleosts. To date, master sex determination genes have been elucidated in 114 species, of which 72 species have sex determination genes belonging to TGF-β superfamily. TGF-β is the only signaling pathway that the largest proportion of components, which including ligands (amhy, gsdfy, gdf6), receptors (amhr, bmpr), and regulator (id2bby), have opportunity recognized as a sex determination gene. In this review, we focus on the recent studies about teleost sex-determination genes within TGF-β superfamily and propose several hypotheses on how these genes regulate sex determination process. Differing from other reviews, our review specifically devotes significant attention to all members of the TGF-β signal pathway, not solely the sex determination genes within the TGF-β superfamily. However, the functions of the paralogous genes of TGF superfamily are still needed ongoing research. Further studies are required to more accurately interpret the molecular mechanism of TGF-β superfamily sex determination genes.

Improved assembly of the Pungitius pungitius reference genome

The nine-spined stickleback (Pungitius pungitius) has been increasingly used as a model system in studies of local adaptation and sex chromosome evolution but its current reference genome assembly is far from perfect, lacking distinct sex chromosomes. We generated an improved assembly of the nine-spined stickleback reference genome (98.3% BUSCO completeness) with the aid of linked-read mapping. While the new assembly (v8) was of similar size as the earlier version (v7), we were able to assign 4.4 times more contigs to the linkage groups and improve the contiguity of the genome. Moreover, the new assembly contains a ∼22.8 Mb Y-linked scaffold (LG22) consisting mainly of previously assigned X-contigs, putative Y-contigs, putative centromere contigs, and highly repetitive elements. The male individual showed an even mapping depth on LG12 (pseudo X chromosome) and LG22 (Y-linked scaffold) in the segregating sites, suggesting near-pure X and Y representation in the v8 assembly. A total of 26,803 genes were annotated, and about 33% of the assembly was found to consist of repetitive elements. The high proportion of repetitive elements in LG22 (53.10%) suggests it can be difficult to assemble the complete sequence of the species' Y chromosome. Nevertheless, the new assembly is a significant improvement over the previous version and should provide a valuable resource for genomic studies of stickleback fishes.

Sex chromosome cycle as a mechanism of stable sex determination

Recent advances in DNA sequencing technology have enabled the precise decoding of genomes in non-model organisms, providing a basis for unraveling the patterns and mechanisms of sex chromosome evolution. Studies of different species have yielded conflicting results regarding the traditional theory that sex chromosomes evolve from autosomes via the accumulation of deleterious mutations and degeneration of the Y (or W) chromosome. The concept of the 'sex chromosome cycle,' emerging from this context, posits that at any stage of the cycle (i.e., differentiation, degeneration, or loss), sex chromosome turnover can occur while maintaining stable sex determination. Thus, understanding the mechanisms that drive both the persistence and turnover of sex chromosomes at each stage of the cycle is crucial. In this review, we integrate recent findings on the mechanisms underlying maintenance and turnover, with a special focus on several organisms having unique sex chromosomes. Our review suggests that the diversity of sex chromosomes in the maintenance of stable sex determination is underappreciated and emphasizes the need for more research on the sex chromosome cycle.

Evolution and regulation of animal sex chromosomes

Animal sex chromosomes typically carry the upstream sex-determining gene that triggers testis or ovary development and, in some species, are regulated by global dosage compensation in response to functional decay of the Y chromosome. Despite the importance of these pathways, they exhibit striking differences across species, raising fundamental questions regarding the mechanisms underlying their evolutionary turnover. Recent studies of non-model organisms, including insects, reptiles and teleosts, have yielded a broad view of the diversity of sex chromosomes that challenges established theories. Moreover, continued studies in model organisms with recently developed technologies have characterized the dynamics of sex determination and dosage compensation in three-dimensional nuclear space and at single-cell resolution. Here, we synthesize recent insights into sex chromosomes from a variety of species to review their evolutionary dynamics with respect to the canonical model, as well as their diverse mechanisms of regulation.

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.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

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.

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.

Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

A phylogenetics-based nomenclature system for steroid receptors in teleost fishes

Teleost fishes have emerged as tractable models for studying the neuroendocrine regulation of social behavior via molecular genetic techniques, such as CRISPR/Cas9 gene editing. Moreover, teleosts provide an opportunity to investigate the evolution of steroid receptors and their functions, as species within this lineage possess novel steroid receptor paralogs that resulted from a teleost-specific whole genome duplication. Although teleost fishes have grown in popularity as models for behavioral neuroendocrinology, there is not a consistent nomenclature system for steroid receptors and their genes, which may impede a clear understanding of steroid receptor paralogs and their functions. Here, we used a phylogenetic approach to assess the relatedness of protein sequences encoding steroid receptor paralogs in 18 species from 12 different orders of the Infraclass Teleostei. While most similarly named sequences grouped based on the established phylogeny of the teleost lineage, our analysis revealed several inconsistencies in the nomenclature of steroid receptor paralogs, particularly for sequences encoding estrogen receptor beta (ERβ). Based on our results, we propose a nomenclature system for teleosts in which Greek symbols refer to proteins and numbers refer to genes encoding different subtypes of steroid receptors within the five major groups of this nuclear receptor subfamily. Collectively, our results bridge a critical gap by providing a cohesive naming system for steroid receptors in teleost fishes, which will serve to improve communication, promote collaboration, and enhance our understanding of the evolution and function of steroid receptors across vertebrates.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

A distinct neurogenomic response to a trade-off between social challenge and opportunity in male sticklebacks (Gasterosteus aculeatus)

Animals frequently make adaptive decisions about what to prioritize when faced with multiple, competing demands simultaneously. However, the proximate mechanisms of decision-making in the face of competing demands are not well understood. We explored this question using brain transcriptomics in a classic model system: threespined sticklebacks, where males face conflict between courtship and territorial defence. We characterized the behaviour and brain gene expression profiles of males confronted by a trade-off between courtship and territorial defence by comparing them to males not confronted by this trade-off. When faced with the trade-off, males behaviourally prioritized defence over courtship, and this decision was reflected in their brain gene expression profiles. A distinct set of genes and biological processes was recruited in the brain when males faced a trade-off and these responses were largely non-overlapping across two brain regions. Combined, these results raise new questions about the interplay between the neural and molecular mechanisms involved in decision-making.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

Evaluating the role of sexual antagonism in the evolution of sex chromosomes: new data from fish

The recent increase in available molecular and genomic data for diverse taxa helps to shed new light on long-standing theories. Research into sex chromosome evolution has particularly benefited from a growing number of studies of fish, motivated by their highly diverse mechanisms of sex determination. Sexual antagonism is regularly cited as an influential force in sex chromosome emergence; however, this so far proves difficult to demonstrate. In this review, we highlight recent developments in the investigation of sexual antagonism in sex chromosome research in fish. We find strong emphasis placed on study-organism specific genomic features and patterns of recombination, rather than evidence for a comprehensive role of sexual antagonism. In this light, we discuss the alternative models of sex chromosome evolution. We conclude that fish represents a key resource for further research, provided attention is given to species-specific effects while simultaneously integrating comparative studies across taxa for a vital and comprehensive understanding of sex chromosome evolution and investigation of proposed models.

Single-Cell RNA Sequencing Reveals Microevolution of the Stickleback Immune System

The risk and severity of pathogen infections in humans, livestock, or wild organisms depend on host immune function, which can vary between closely related host populations or even among individuals. This immune variation can entail between-population differences in immune gene coding sequences, copy number, or expression. In recent years, many studies have focused on population divergence in immunity using whole-tissue transcriptomics. But, whole-tissue transcriptomics cannot distinguish between evolved differences in gene regulation within cells, versus changes in cell composition within the focal tissue. Here, we leverage single-cell transcriptomic approaches to document signatures of microevolution of immune system structure in a natural system, the three-spined stickleback (Gasterosteus aculeatus). We sampled nine adult fish from three populations with variability in resistance to a cestode parasite, Schistocephalus solidus, to create the first comprehensive immune cell atlas for G. aculeatus. Eight broad immune cell types, corresponding to major vertebrate immune cells, were identified. We were also able to document significant variation in both abundance and expression profiles of the individual immune cell types among the three populations of fish. Furthermore, we demonstrate that identified cell type markers can be used to reinterpret traditional transcriptomic data: we reevaluate previously published whole-tissue transcriptome data from a quantitative genetic experimental infection study to gain better resolution relating infection outcomes to inferred cell type variation. Our combined study demonstrates the power of single-cell sequencing to not only document evolutionary phenomena (i.e., microevolution of immune cells) but also increase the power of traditional transcriptomic data sets.

Molecular parallelism in the evolution of a master sex-determining role for the anti-Mullerian hormone receptor 2 gene (amhr2) in Midas cichlids

The evolution of sex chromosomes and their differentiation from autosomes is a major event during genome evolution that happened many times in several lineages. The repeated evolution and lability of sex-determination mechanisms in fishes makes this a well-suited system to test for general patterns in evolution. According to current theory, differentiation is triggered by the suppression of recombination following the evolution of a new master sex-determining gene. However, the molecular mechanisms that establish recombination suppression are known from few examples, owing to the intrinsic difficulties of assembling sex-determining regions (SDRs). The development of forward-genetics and long-read sequencing have generated a wealth of data questioning central aspects of the current theory. Here, we demonstrate that sex in Midas cichlids is determined by an XY system, and identify and assemble the SDR by combining forward-genetics, long-read sequencing and optical mapping. We show how long-reads aid in the detection of artefacts in genotype-phenotype mapping that arise from incomplete genome assemblies. The male-specific region is restricted to a 100-kb segment on chromosome 4 that harbours transposable elements and a Y-specific duplicate of the anti-Mullerian receptor 2 gene, which has evolved master sex-determining functions repeatedly. Our data suggest that amhr2Y originated by an interchromosomal translocation from chromosome 20 to 4 pre-dating the split of Midas and Flier cichlids. In the latter, it is pseudogenized and translocated to another chromosome. Duplication of anti-Mullerian genes is a common route to establishing new sex determiners, highlighting the role of molecular parallelism in the evolution of sex determination.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole-genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h²  = 0.42-0.65), indicating adaptive potential for all traits. Environment-association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation - supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

Sex chromosome differentiation via changes in the Y chromosome repeat landscape in African annual killifishes Nothobranchius furzeri and N. kadleci

Homomorphic sex chromosomes and their turnover are common in teleosts. We investigated the evolution of nascent sex chromosomes in several populations of two sister species of African annual killifishes, Nothobranchius furzeri and N. kadleci, focusing on their under-studied repetitive landscape. We combined bioinformatic analyses of the repeatome with molecular cytogenetic techniques, including comparative genomic hybridization, fluorescence in situ hybridization with satellite sequences, ribosomal RNA genes (rDNA) and bacterial artificial chromosomes (BACs), and immunostaining of SYCP3 and MLH1 proteins to mark lateral elements of synaptonemal complexes and recombination sites, respectively. Both species share the same heteromorphic XY sex chromosome system, which thus evolved prior to their divergence. This was corroborated by sequence analysis of a putative master sex determining (MSD) gene gdf6Y in both species. Based on their divergence, differentiation of the XY sex chromosome pair started approximately 2 million years ago. In all populations, the gdf6Y gene mapped within a region rich in satellite DNA on the Y chromosome long arms. Despite their heteromorphism, X and Y chromosomes mostly pair regularly in meiosis, implying synaptic adjustment. In N. kadleci, Y-linked paracentric inversions like those previously reported in N. furzeri were detected. An inversion involving the MSD gene may suppress occasional recombination in the region, which we otherwise evidenced in the N. furzeri population MZCS-121 of the Limpopo clade lacking this inversion. Y chromosome centromeric repeats were reduced compared with the X chromosome and autosomes, which points to a role of relaxed meiotic drive in shaping the Y chromosome repeat landscape. We speculate that the recombination rate between sex chromosomes was reduced due to heterochiasmy. The observed differences between the repeat accumulations on the X and Y chromosomes probably result from high repeat turnover and may not relate closely to the divergence inferred from earlier SNP analyses.

Meiotic pairing and double-strand break formation along the heteromorphic threespine stickleback sex chromosomes

Double-strand break repair during meiosis is normally achieved using the homologous chromosome as a repair template. Heteromorphic sex chromosomes share little sequence homology, presenting unique challenges to the repair of double-strand breaks. Our understanding of how heteromorphic sex chromosomes behave during meiosis has been focused on ancient sex chromosomes, where the X and Y differ markedly in overall structure and gene content. It remains unclear how more recently evolved sex chromosomes that share considerably more sequence homology with one another pair and form double-strand breaks. One possibility is barriers to pairing evolve rapidly. Alternatively, recently evolved sex chromosomes may exhibit pairing and double-strand break repair that more closely resembles that of their autosomal ancestors. Here, we use the recently evolved X and Y chromosomes of the threespine stickleback fish (Gasterosteus aculeatus) to study patterns of pairing and double-stranded break formation using molecular cytogenetics. We found that the sex chromosomes of threespine stickleback fish did not pair exclusively in the pseudoautosomal region. Instead, the chromosomes fully paired in a non-homologous fashion. To achieve this, the X chromosome underwent synaptic adjustment during pachytene to match the axis length of the Y chromosome. Double-strand break formation and repair rate also matched that of the autosomes. Our results highlight that recently evolved sex chromosomes exhibit meiotic behavior that is reminiscent of autosomes and argues for further work to identify the homologous templates that are used to repair double-strand breaks on the X and Y chromosomes.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

Sexually mediated phenotypic variation within and between sexes as a continuum structured by ecology: The mosaic nature of skeletal variation across body regions in Threespine stickleback (Gasterosteus aculeatus L.)

Ecological character displacement between the sexes, and sexual selection, integrate into a convergent set of factors that produce sexual variation. Ecologically modulated, sexually mediated variation within and between sexes may be a major contributor to the amount of total variation that selection can act on in species. Threespine stickleback (Gasterosteus aculeatus) display rapid adaptive responses and sexual variation in many phenotypic traits. We examined phenotypic variation in the skull, pectoral and pelvic girdles of threespine stickleback from two freshwater and two coastal marine sites on the Sunshine Coast of British Columbia, Canada, using an approach that avoids a priori assumptions about bimodal patterns of variation. We quantified shape and size of the cranial, pectoral and pelvic regions of sticklebacks in marine and freshwater habitats using 3D geometric morphometrics and an index of sexually mediated variation. We show that the expression of phenotypic variation is structured in part by the effects of both habitat marine vs freshwater and the effects of individual sites within each habitat. Relative size exerts variable influence, and patterns of phenotypic variation associated with sex vary among body regions. This fine-grained quantification of sexually mediated variation in the context of habitat difference and different anatomical structures indicates a complex relationship between genetically inferred sex and environmental factors, demonstrating that the interplay between shared genetic background and sexually mediated, ecologically based selective pressures structures the phenotypic expression of complex traits.

Searching for signatures of sexually antagonistic selection on stickleback sex chromosomes

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

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

Assessment of phylogenetic approaches to study the timing of recombination cessation on sex chromosomes

The evolution of sex chromosomes is hypothesized to be punctuated by consecutive recombination cessation events, forming "evolutionary strata" that ceased to recombine at different time points. The demarcation of evolutionary strata is often assessed by estimates of the timing of recombination cessation (t_RC ) along the sex chromosomes, commonly inferred from the level of synonymous divergence or with species phylogenies at gametologous (X-Y or Z-W) sequence data. However, drift and selection affect sequences unpredictably and introduce uncertainty when inferring t_RC . Here, we assess two alternative phylogenetic approaches to estimate t_RC ; (i) the expected likelihood weight (ELW) approach that finds the most likely topology among a set of hypothetical topologies and (ii) the BEAST approach that estimates t_RC with specified calibration priors on a reference species topology. By using Z and W gametologs of an old and a young evolutionary stratum on the neo-sex chromosome of Sylvioidea songbirds, we show that the ELW and BEAST approaches yield similar t_RC estimates, and that both outperform two frequently applied approaches utilizing synonymous substitution rates (dS) and maximum likelihood (ML) trees, respectively. Moreover, we demonstrate that both ELW and BEAST provide more precise t_RC estimates when sequences of multiple species are included in the analyses.

Repetitive DNA Sequences in the Human Y Chromosome and Male Infertility

The male-specific Y chromosome, which is well known for its diverse and complex repetitive sequences, has different sizes, genome structures, contents and evolutionary trajectories from other chromosomes and is of great significance for testis development and function. The large number of repetitive sequences and palindrome structure of the Y chromosome play an important role in maintaining the stability of male sex determining genes, although they can also cause non-allelic homologous recombination within the chromosome. Deletion of certain Y chromosome sequences will lead to spermatogenesis disorders and male infertility. And Y chromosome genes are also involved in the occurrence of reproductive system cancers and can increase the susceptibility of other tumors. In addition, the Y chromosome has very special value in the personal identification and parentage testing of male-related cases in forensic medicine because of its unique paternal genetic characteristics. In view of the extremely high frequency and complexity of gene rearrangements and the limitations of sequencing technology, the analysis of Y chromosome sequences and the study of Y-gene function still have many unsolved problems. This article will introduce the structure and repetitive sequence of the Y chromosome, summarize the correlation between Y chromosome various sequence deletions and male infertility for understanding the repetitive sequence of Y chromosome more systematically, in order to provide research motivation for further explore of the molecules mechanism of Y-deletion and male infertility and theoretical foundations for the transformation of basic research into applications in clinical medicine and forensic medicine.

Identification of a candidate sex determination gene in Culaea inconstans suggests convergent recruitment of an Amh duplicate in two lineages of stickleback

Sex chromosomes vary greatly in their age and levels of differentiation across the tree of life. This variation is largely due to the rates of sex chromosome turnover in different lineages; however, we still lack an explanation for why sex chromosomes are so conserved in some lineages (e.g. mammals, birds) but so labile in others (e.g. teleosts, amphibians). To identify general mechanisms driving transitions in sex determination systems or forces which favour their conservation, we first require empirical data on sex chromosome systems from multiple lineages. Stickleback fishes are a valuable model lineage for the study of sex chromosome evolution due to variation in sex chromosome systems between closely-related species. Here, we identify the sex chromosome and a strong candidate for the master sex determination gene in the brook stickleback, Culaea inconstans. Using whole-genome sequencing of wild-caught samples and a lab cross, we identify AmhY, a male specific duplication of the gene Amh, as the candidate master sex determination gene. AmhY resides on Chromosome 20 in C. inconstans and is likely a recent duplication, as both AmhY and the sex-linked region of Chromosome 20 show little sequence divergence. Importantly, this duplicate AmhY represents the second independent duplication and recruitment of Amh as the sex determination gene in stickleback and the eighth example known across teleosts. We discuss this convergence in the context of sex chromosome turnovers and the role that the Amh/AmhrII pathway, which is crucial for sex determination, may play in the evolution of sex chromosomes in teleosts.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish

SignificanceAlthough turnover of sex chromosomes is very common in many vertebrate lineages, the transition process is still elusive. We studied the sex-determining region (SDR) of 12 congeneric fish species of Takifugu fish that compose an important model for the study of genomics and sex determination. We found that while nine species retained their ancestral SDR, three species had acquired derived SDRs. Although the derived SDRs resided in three different chromosomes, they harbored a shared supergene flanked by two putative transposable elements. The results highlight the underestimated role of a mobile supergene in turnover of sex chromosomes in vertebrates.

The evolution of widespread recombination suppression on the dwarf hamster (Phodopus) X chromosome

The X chromosome of therian mammals shows strong conservation among distantly related species, limiting insights into the distinct selective processes that have shaped sex chromosome evolution. We constructed a chromosome-scale de novo genome assembly for the Siberian dwarf hamster (Phodopus sungorus), a species reported to show extensive recombination suppression across an entire arm of the X chromosome. Combining a physical genome assembly based on shotgun and long-range proximity ligation sequencing with a dense genetic map, we detected widespread suppression of female recombination across ∼65% of the Phodopus X chromosome. This region of suppressed recombination likely corresponds to the Xp arm, which has previously been shown to be highly heterochromatic. Using additional sequencing data from two closely related species (P. campbelli and P. roborovskii), we show that recombination suppression on Xp appears to be independent of major structural rearrangements. The suppressed Xp arm was enriched for several transposable element families and de-enriched for genes primarily expressed in placenta, but otherwise showed similar gene densities, expression patterns, and rates of molecular evolution when compared to the recombinant Xq arm. Phodopus Xp gene content and order was also broadly conserved relative to the more distantly related rat X chromosome. These data suggest that widespread suppression of recombination has likely evolved through the transient induction of facultative heterochromatin on the Phodopus Xp arm without major changes in chromosome structure or genetic content. Thus, substantial changes in the recombination landscape have so far had relatively subtle influences on patterns of X-linked molecular evolution in these species.

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 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’.

Identification of sex chromosome and sex-determining gene of southern catfish (Silurus meridionalis) based on XX, XY and YY genome sequencing

Teleosts are important models to study sex chromosomes and sex-determining (SD) genes because they present a variety of sex determination systems. Here, we used Nanopore and Hi-C technologies to generate a high-contiguity chromosome-level genome assembly of a YY southern catfish (Silurus meridionalis). The assembly is 750.0 Mb long, with contig N50 of 15.96 Mb and scaffold N50 of 27.22 Mb. We also sequenced and assembled an XY male genome with a size of 727.2 Mb and contig N50 of 13.69 Mb. We identified a candidate SD gene through comparisons to our previous assembly of an XX individual. By resequencing male and female pools, we characterized a 2.38 Mb sex-determining region (SDR) on Chr24. Analysis of read coverage and comparison of the X and Y chromosome sequences showed a Y specific insertion (approx. 500 kb) in the SDR which contained a male-specific duplicate of amhr2 (named amhr2y). amhr2y and amhr2 shared high-nucleotide identity (81.0%) in the coding region but extremely low identity in the promotor and intron regions. The exclusive expression in the male gonadal primordium and loss-of-function inducing male to female sex reversal confirmed the role of amhr2y in male sex determination. Our study provides a new example of amhr2 as the SD gene in fish and sheds light on the convergent evolution of the duplication of AMH/AMHR2 pathway members underlying the evolution of sex determination in different fish lineages.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Unique structure and positive selection promote the rapid divergence of Drosophila Y chromosomes

Y chromosomes across diverse species convergently evolve a gene-poor, heterochromatic organization enriched for duplicated genes, LTR retrotransposons, and satellite DNA. Sexual antagonism and a loss of recombination play major roles in the degeneration of young Y chromosomes. However, the processes shaping the evolution of mature, already degenerated Y chromosomes are less well-understood. Because Y chromosomes evolve rapidly, comparisons between closely related species are particularly useful. We generated de novo long-read assemblies complemented with cytological validation to reveal Y chromosome organization in three closely related species of the Drosophila simulans complex, which diverged only 250,000 years ago and share >98% sequence identity. We find these Y chromosomes are divergent in their organization and repetitive DNA composition and discover new Y-linked gene families whose evolution is driven by both positive selection and gene conversion. These Y chromosomes are also enriched for large deletions, suggesting that the repair of double-strand breaks on Y chromosomes may be biased toward microhomology-mediated end joining over canonical non-homologous end-joining. We propose that this repair mechanism contributes to the convergent evolution of Y chromosome organization across organisms.

Muller's ratchet of the Y chromosome with gene conversion

Muller's ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller's ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller's ratchet to explore the evolution of the Y chromosome. The model assumes a nonrecombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.

Roles of anti-Müllerian hormone and its duplicates in sex determination and germ cell proliferation of Nile tilapia

Duplicates of amh are crucial for fish sex determination and differentiation. In Nile tilapia, unlike in other teleosts, amh is located on X chromosome. The Y chromosome amh (amhΔ-y) is mutated with 5 bp insertion and 233 bp deletion in the coding sequence, and tandem duplicate of amh on Y chromosome (amhy) has been identified as the sex determiner. However, the expression of amh, amhΔ-y, and amhy, their roles in germ cell proliferation and the molecular mechanism of how amhy determines sex is still unclear. In this study, expression and functions of each duplicate were analyzed. Sex reversal occurred only when amhy was mutated as revealed by single, double, and triple mutation of the 3 duplicates in XY fish. Homozygous mutation of amhy in YY fish also resulted in sex reversal. Earlier and higher expression of amhy/Amhy was observed in XY gonads compared with amh/Amh during sex determination. Amhy could inhibit the transcription of cyp19a1a through Amhr2/Smads signaling. Loss of cyp19a1a rescued the sex reversal phenotype in XY fish with amhy mutation. Interestingly, mutation of both amh and amhy in XY fish or homozygous mutation of amhy in YY fish resulted in infertile females with significantly increased germ cell proliferation. Taken together, these results indicated that up-regulation of amhy during the critical period of sex determination makes it the sex-determining gene, and it functions through repressing cyp19a1a expression via Amhr2/Smads signaling pathway. Amh retained its function in controlling germ cell proliferation as reported in other teleosts, while amhΔ-y was nonfunctionalized.

Chromosomal fusions facilitate adaptation to divergent environments in threespine stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Reference genome of Lumpfish Cyclopterus lumpus Linnaeus provides evidence of male heterogametic sex determination through the AMH pathway

Teleosts exhibit extensive diversity of sex determination (SD) systems and mechanisms, providing the opportunity to study the evolution of sex determination and sex chromosomes. Here we sequenced the genome of the Common Lumpfish (Cyclopterus lumpus Linnaeus), a species of increasing importance to aquaculture, and identified the SD region and master SD locus using a 70K SNP array and tissue-specific expression data. The chromosome-level assembly identified 25 diploid chromosomes with a total size of 572.89 Mb, a scaffold N50 of 23.86 Mb, and genome annotation predicted 21,480 protein-coding genes. Genome wide association analysis located a highly sex-associated region on chromosome 13, suggesting that anti-Müllerian hormone (AMH) is the putative SD factor. Linkage disequilibrium and heterozygosity across chromosome 13 support a proto-XX/XY system, with an absence of widespread chromosome divergence between sexes. We identified three copies of AMH in the Lumpfish primary and alternate haplotype assemblies localized in the SD region. Comparison to sequences from other teleosts suggested a monophyletic relationship and conservation within the Cottioidei. One AMH copy showed similarity to AMH/AMHY in a related species and was also the only copy with expression in testis tissue, suggesting this copy may be the functional copy of AMH in Lumpfish. The two other copies arranged in tandem inverted duplication were highly similar, suggesting a recent duplication event. This study provides a resource for the study of early sex chromosome evolution and novel genomic resources that benefits Lumpfish conservation management and aquaculture.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

Against the mainstream: exceptional evolutionary stability of ZW sex chromosomes across the fish families Triportheidae and Gasteropelecidae (Teleostei: Characiformes)

Teleost fishes exhibit a breath-taking diversity of sex determination and differentiation mechanisms. They encompass at least nine sex chromosome systems with often low degree of differentiation, high rate of inter- and intra-specific variability, and frequent turnovers. Nevertheless, several mainly female heterogametic systems at an advanced stage of genetic differentiation and high evolutionary stability have been also found across teleosts, especially among Neotropical characiforms. In this study, we aim to characterize the ZZ/ZW sex chromosome system in representatives of the Triportheidae family (Triportheus auritus, Agoniates halecinus, and the basal-most species Lignobrycon myersi) and its sister clade Gasteropelecidae (Carnegiella strigata, Gasteropelecus levis, and Thoracocharax stellatus). We applied both conventional and molecular cytogenetic approaches including chromosomal mapping of 5S and 18S ribosomal DNA clusters, cross-species chromosome painting (Zoo-FISH) with sex chromosome-derived probes and comparative genomic hybridization (CGH). We identified the ZW sex chromosome system for the first time in A. halecinus and G. levis and also in C. strigata formerly reported to lack sex chromosomes. We also brought evidence for possible mechanisms underlying the sex chromosome differentiation, including inversions, repetitive DNA accumulation, and exchange of genetic material. Our Zoo-FISH experiments further strongly indicated that the ZW sex chromosomes of Triportheidae and Gasteropelecidae are homeologous, suggesting their origin before the split of these lineages (approx. 40-70 million years ago). Such extent of sex chromosome stability is almost exceptional in teleosts, and hence, these lineages afford a special opportunity to scrutinize unique evolutionary forces and pressures shaping sex chromosome evolution in fishes and vertebrates in general.

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.

Adding New Pieces to the Puzzle of Karyotype Evolution in Harttia (Siluriformes, Loricariidae): Investigation of Amazonian Species

A remarkable morphological diversity and karyotype variability can be observed in the Neotropical armored catfish genus Harttia. These fishes offer a useful model to explore both the evolution of karyotypes and sex chromosomes, since many species possess male-heterogametic sex chromosome systems and a high rate of karyotype repatterning. Based on the karyotype organization, the chromosomal distribution of several repetitive DNA classes, and the rough estimates of genomic divergences at the intraspecific and interspecific levels via Comparative Genomic Hybridization, we identified shared diploid chromosome numbers (2n = 54) but different karyotype compositions in H. dissidens (20m + 26sm + 8a) and Harttia sp. 3 (16m + 18sm + 14st + 6a), and different 2n in H. guianensis (2n = 58; 20m + 26sm + 2st + 10a). All species further displayed similar patterns of chromosomal distribution concerning constitutive heterochromatin, 18S ribosomal DNA (rDNA) sites, and most of the surveyed microsatellite motifs. Furthermore, differences in the distribution of 5S rDNA sites and a subset of microsatellite sequences were identified. Heteromorphic sex chromosomes were lacking in H. dissidens and H. guianensis at the scale of our analysis. However, one single chromosome pair in Harttia sp. 3 males presented a remarkable accumulation of male genome-derived probe after CGH, pointing to a tentative region of early sex chromosome differentiation. Thus, our data support already previously outlined evidence that Harttia is a vital model for the investigation of teleost karyotype and sex chromosome dynamics.

Expanding the classical paradigm: what we have learnt from vertebrates about sex chromosome evolution

Until recently, the field of sex chromosome evolution has been dominated by the canonical unidirectional scenario, first developed by Muller in 1918. This model postulates that sex chromosomes emerge from autosomes by acquiring a sex-determining locus. Recombination reduction then expands outwards from this locus, to maintain its linkage with sexually antagonistic/advantageous alleles, resulting in Y or W degeneration and potentially culminating in their disappearance. Based mostly on empirical vertebrate research, we challenge and expand each conceptual step of this canonical model and present observations by numerous experts in two parts of a theme issue of Phil. Trans. R. Soc. B. We suggest that greater theoretical and empirical insights into the events at the origins of sex-determining genes (rewiring of the gonadal differentiation networks), and a better understanding of the evolutionary forces responsible for recombination suppression are required. Among others, crucial questions are: Why do sex chromosome differentiation rates and the evolution of gene dose regulatory mechanisms between male versus female heterogametic systems not follow earlier theory? Why do several lineages not have sex chromosomes? And: What are the consequences of the presence of (differentiated) sex chromosomes for individual fitness, evolvability, hybridization and diversification? We conclude that the classical scenario appears too reductionistic. Instead of being unidirectional, we show that sex chromosome evolution is more complex than previously anticipated and principally forms networks, interconnected to potentially endless outcomes with restarts, deletions and additions of new genomic material. 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 II)'.