Steps in the evolution of heteromorphic sex chromosomes

Genetic drift drives faster-Z evolution in the salmon louse Lepeophtheirus salmonis

How sex chromosomes evolve compared to autosomes remains an unresolved question in population genetics. Most studies focus on only a handful of taxa, resulting in uncertainty over whether observed patterns reflect general processes or idiosyncrasies in particular clades. For example, in female heterogametic (ZW) systems, bird Z chromosomes tend to evolve quickly but not adaptively, while in Lepidopterans they evolve adaptively, but not always quickly. To understand how these observations fit into broader evolutionary patterns, we explore Z chromosome evolution outside of these two well-studied clades. We utilize a publicly available genome, gene expression, population, and outgroup data in the salmon louse Lepeophtheirus salmonis, an important agricultural pest copepod. We find that the Z chromosome is faster evolving than autosomes, but that this effect is driven by increased drift rather than adaptive evolution. Due to high rates of female reproductive failure, the Z chromosome exhibits a slightly lower effective population size than the autosomes which is nonetheless to decrease efficiency of hemizygous selection acting on the Z. These results highlight the usefulness of organismal life history in calibrating population genetic expectations and demonstrate the value of the ever-expanding wealth of publicly available data to help resolve outstanding evolutionary questions.

Satellite DNAs and the evolution of the multiple X1X2Y sex chromosomes in the wolf fish Hoplias malabaricus (Teleostei; Characiformes)

Multiple sex chromosomes usually arise from chromosomal rearrangements which involve ancestral sex chromosomes. There is a fundamental condition to be met for their long-term fixation: the meiosis must function, leading to the stability of the emerged system, mainly concerning the segregation of the sex multivalent. Here, we sought to analyze the degree of differentiation and meiotic pairing properties in the selected fish multiple sex chromosome system present in the wolf-fish Hoplias malabaricus (HMA). This species complex encompasses seven known karyotype forms (karyomorphs) where the karyomorph C (HMA-C) exhibits a nascent XY sex chromosomes from which the multiple X1X2Y system evolved in karyomorph HMA-D via a Y-autosome fusion. We combined genomic and cytogenetic approaches to analyze the satellite DNA (satDNA) content in the genome of HMA-D karyomorph and to investigate its potential contribution to X1X2Y sex chromosome differentiation. We revealed 56 satDNA monomers of which the majority was AT-rich and with repeat units longer than 100 bp. Seven out of 18 satDNA families chosen for chromosomal mapping by fluorescence in situ hybridization (FISH) formed detectable accumulation in at least one of the three sex chromosomes (X1, X2 and neo-Y). Nine satDNA monomers showed only two hybridization signals limited to HMA-D autosomes, and the two remaining ones provided no visible FISH signals. Out of seven satDNAs located on the HMA-D sex chromosomes, five mapped also to XY chromosomes of HMA-C. We showed that after the autosome-Y fusion event, the neo-Y chromosome has not substantially accumulated or eliminated satDNA sequences except for minor changes in the centromere-proximal region. Finally, based on the obtained FISHpatterns, we speculate on the possible contribution of satDNA to sex trivalent pairing and segregation.

Sex chromosome turnover and biodiversity in fishes

The impact of sex chromosomes and their turnover in speciation remains a subject of ongoing debate in the field of evolutionary biology. Fishes are the largest group of vertebrates, and they exhibit unparalleled sexual plasticity, as well as diverse sex-determining (SD) genes, sex chromosomes, and sex determination mechanisms. This diversity is hypothesized to be associated with the frequent turnover of sex chromosomes in fishes. Although it is evident that amh and amhr2 are repeatedly and independently recruited as SD genes, their relationship with the rapid turnover of sex chromosomes and the biodiversity of fishes remains unknown. We summarize the canonical models of sex chromosome turnover and highlight the vital roles of gene mutation and hybridization with empirical evidence. We revisit Haldane's rule and the large X-effect and propose the hypothesis that sex chromosomes accelerate speciation by multiplying genotypes via hybridization. By integrating recent findings on the turnover of SD genes, sex chromosomes, and sex determination systems in fish species, this review provides insights into the relationship between sex chromosome evolution and biodiversity in fishes.

The role of satellite DNAs in the chromosomal rearrangements and the evolution of the rare XY1Y2 sex system in Harttia (Siluriformes: Loricariidae)

The underlying processes behind the formation, evolution, and long-term maintenance of multiple sex chromosomes have been largely neglected. Among vertebrates, fishes represent the group with the highest diversity of multiple sex chromosome systems and, with six instances, the Neotropical fish genus Harttia stands out by presenting the most remarkable diversity. However, although the origin mechanism of their sex chromosome systems is well discussed, little is known about the importance of some repetitive DNA classes in the differentiation of multiple systems. In this work, by employing a combination of cytogenetic and genomic procedures, we evaluated the satellite DNA composition of H. carvalhoi with a focus on their role in the evolution, structure, and differentiation process of the rare XY1Y2 multiple-sex chromosome system. The genome of H. carvalhoi contains a total of 28 satellite DNA families, with the A + T content ranging between 38.1% and 68.1% and the predominant presence of long satellites. The in situ hybridization experiments detected 15 satellite DNAs with positive hybridization signals mainly on centromeric and pericentromeric regions of almost all chromosomes or clustered on a few pairs. Five of them presented clusters on X, Y1, and/or Y2 sex chromosomes which were therefore selected for comparative hybridization in the other three congeneric species. We found several conserved satellites accumulated on sex chromosomes and also in regions that were involved in chromosomal rearrangements. Our results provide a new contribution of satellitome studies in multiple sex chromosome systems in fishes and represent the first satellitome study for a Siluriformes species.

Compensation of gene dosage on the mammalian X

Changes in gene dosage can have tremendous evolutionary potential (e.g. whole-genome duplications), but without compensatory mechanisms, they can also lead to gene dysregulation and pathologies. Sex chromosomes are a paradigmatic example of naturally occurring gene dosage differences and their compensation. In species with chromosome-based sex determination, individuals within the same population necessarily show 'natural' differences in gene dosage for the sex chromosomes. In this Review, we focus on the mammalian X chromosome and discuss recent new insights into the dosage-compensation mechanisms that evolved along with the emergence of sex chromosomes, namely X-inactivation and X-upregulation. We also discuss the evolution of the genetic loci and molecular players involved, as well as the regulatory diversity and potentially different requirements for dosage compensation across mammalian species.

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.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes1

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Although difficult to test empirically, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering" hypothesis which suggests that expansion of the sex-linked region (SLR) reduces the homozygous expression of deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of partially recessive deleterious mutational variation for the evolution of otherwise neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. Both autosomal and SLR-expanding inversions face a race against time: lightly-loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, after which their chances of fixing become negligible. In contrast, initially unloaded inversions eventually become neutral as their deleterious load reaches the same equilibrium as non-inverted haplotypes. Despite the differences in inheritance and indirect selection, SLR-expanding inversions exhibit similar evolutionary dynamics to autosomal inversions over many biologically plausible parameter conditions. Differences emerge when the population average mutation load is quite high; in this case large autosomal inversions that are lucky enough to be mutation-free can rise to intermediate to high frequencies where selection in homozygotes becomes important (Y-linked inversions never appear as homozygous karyotypes); conditions requiring either high mutation rates, highly recessive deleterious mutations, weak selection, or a combination thereof.

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.

Carica papaya L. sex chromosome review and physical mapping of the serk 2, svp-like and mdar 4 sequences

Physical mapping evidences the chromosome organization and structure. Despite the data about plant cytogenomics, physical mapping has been conducted from single-copy and/or low-copy genes for few species. Carica papaya cytogenomics has been accomplished from BAC-FISH and repeatome sequences. We aimed to map the serk 2, svp-like and mdar 4 sequences in C. papaya. The sequences were amplified and the amplicons sequenced, showing similarity in relation to serk 2, svp-like and mdar 4 genes. Carica papaya diploidy was confirmed and the mitotic chromosomes characterized. The chromosome 1 exhibited the secondary constriction pericentromeric to the centromere of the long arm. So, we concluded that it is the sex chromosomes. serk 2 was mapped in the long arm interstitial portion of the sex chromosomes, and the interphase nuclei showed two fluorescence signals. Considering these results and the sequencing data from the C. papaya sex chromosomes, svp-like and mdar 4 genes were mapped in the interstitial region of the sex chromosome long arm. Both sequences showed only one fluorescence signal in the interphase nuclei. The procedure adopted here can be reproduced for other single-copy and/or low-copy genes, allowing the construction of cytogenetic maps. In addition, we revisited the cytogenomics data about C. papaya sex chromosomes, presenting a revised point of view about the structure and evolution to these chromosomes.

Gap-free X and Y chromosome assemblies of Salix arbutifolia reveal an evolutionary change from male to female heterogamety in willows, without a change in the position of the sex-determining locus

In the Vetrix clade of Salix, a genus of woody flowering plants, sex determination involves chromosome 15, but an XY system has changed to a ZW system. We studied the detailed genetic changes involved. We used genome sequencing, with chromosome conformation capture (Hi-C) and PacBio HiFi reads to assemble chromosome level gap-free X and Y of Salix arbutifolia, and distinguished the haplotypes in the 15X- and 15Y-linked regions, to study the evolutionary history of the sex-linked regions (SLRs). Our sequencing revealed heteromorphism of the X and Y haplotypes of the SLR, with the X-linked region being considerably larger than the corresponding Y region, mainly due to accumulated repetitive sequences and gene duplications. The phylogenies of single-copy orthogroups within the SLRs indicate that S. arbutifolia and Salix purpurea share an ancestral SLR within a repeat-rich region near the chromosome 15 centromere. During the change in heterogamety, the X-linked region changed to a W-linked one, while the Z was derived from the Y.

Initial Steps of XY Sex Chromosome Differentiation in the Armored Catfish Hypostomus albopunctatus (Siluriformes: Loricariidae) Revealed by Heterochromatin Accumulation

In fish species, heterochromatinization is one process that could trigger sex chromosome differentiation. The present article describes a nascent XX/XY sex chromosome system evidenced by heterochromatin accumulation and microsatellite (GATA)8 in Hypostomus albopunctatus from two populations of the Paraná River basin. The specimens of H. albopunctatus from the Campo and Bossi Rivers share the same karyotype. The species exhibits 74 chromosomes (8m+14sm +16st +36a, fundamental number = 112). The C-banding technique suggests male heterogamety in H. albopunctatus, where the Y-chromosome is morphologically like the X-chromosome but differs from it for having long arms that are entirely heterochromatic. Double fluorescence in situ hybridization (FISH) with 18S and 5S rDNA probes confirmed the Ag-nucleolus organizer region sites in a single pair for both populations, and minor rDNA clusters showed interpopulational variation. FISH with the microsatellite (GATA)8 probe showed a dispersed pattern in the karyotype, accumulating these sequences of sex chromosomes of both populations. FISH with microsatellite (CGC)10 probe showed interpopulational variation. The absence of differentiated sex chromosomes in H. albopunctatus is described previously, and a new variant is documented herein where XY chromosomes can be seen in an early stage of differentiation.

Evolution of bird sex chromosomes: a cytogenomic approach in Palaeognathae species

BACKGROUND

Different patterns of sex chromosome differentiation are seen in Palaeognathae birds, a lineage that includes the ratites (Struthioniformes, Rheiformes, Apterygiformes, Casuariiformes, and the sister group Tinamiformes). While some Tinamiform species have well-differentiated W chromosomes, both Z and W of all the flightless ratites are still morphologically undifferentiated. Here, we conducted a comprehensive analysis of the ZW differentiation in birds using a combination of cytogenetic, genomic, and bioinformatic approaches. The whole set of satDNAs from the emu (Dromaius novaehollandiae) was described and characterized. Furthermore, we examined the in situ locations of these satDNAs alongside several microsatellite repeats and carried out Comparative Genomic Hybridizations in two related species: the greater rhea (Rhea americana) and the tataupa tinamou (Crypturellus tataupa).

RESULTS

From the 24 satDNA families identified (which represent the greatest diversity of satDNAs ever uncovered in any bird species), only three of them were found to accumulate on the emu's sex chromosomes, with no discernible accumulation observed on the W chromosome. The W chromosomes of both the greater rhea and the emu did not exhibit a significant buildup of either C-positive heterochromatin or repetitive DNAs, indicating their large undifferentiation both at morphological and molecular levels. In contrast, the tataupa tinamou has a highly differentiated W chromosome that accumulates several DNA repeats.

CONCLUSION

The findings provide new information on the architecture of the avian genome and an inside look at the starting points of sex chromosome differentiation in birds.

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.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.

Cephalopod Sex Determination and its Ancient Evolutionary Origin Revealed by Chromosome-level Assembly of the California Two-Spot Octopus

Sex chromosomes are critical elements of sexual reproduction in many animal and plant taxa, however they show incredible diversity and rapid turnover even within clades. Here, using a chromosome-level assembly generated with long read sequencing, we report the first evidence for genetic sex determination in cephalopods. We have uncovered a sex chromosome in California two-spot octopus (Octopus bimaculoides) in which males/females show ZZ/ZO karyotypes respectively. We show that the octopus Z chromosome is an evolutionary outlier with respect to divergence and repetitive element content as compared to other chromosomes and that it is present in all coleoid cephalopods that we have examined. Our results suggest that the cephalopod Z chromosome originated between 455 and 248 million years ago and has been conserved to the present, making it the among the oldest conserved animal sex chromosomes known.

Turnover of sex chromosomes in the Lake Tanganyika cichlid tribe Tropheini (Teleostei: Cichlidae)

Sex chromosome replacement is frequent in many vertebrate clades, including fish, frogs, and lizards. In order to understand the mechanisms responsible for sex chromosome turnover and the early stages of sex chromosome divergence, it is necessary to study lineages with recently evolved sex chromosomes. Here we examine sex chromosome evolution in a group of African cichlid fishes (tribe Tropheini) which began to diverge from one another less than 4 MYA. We have evidence for a previously unknown sex chromosome system, and preliminary indications of several additional systems not previously reported in this group. We find a high frequency of sex chromosome turnover and estimate a minimum of 14 turnovers in this tribe. We date the origin of the most common sex determining system in this tribe (XY-LG5/19) near the base of one of two major sub-clades of this tribe, about 3.4 MY ago. Finally, we observe variation in the size of one sex-determining region that suggests independent evolution of evolutionary strata in species with a shared sex-determination system. Our results illuminate the rapid rate of sex chromosome turnover in the tribe Tropheini and set the stage for further studies of the dynamics of sex chromosome evolution in this group.

X Chromosome-Specific Repeats in Non-Domestic Bovidae

Repetitive sequences form a substantial and still enigmatic part of the mammalian genome. We isolated repetitive DNA blocks of the X chromosomes of three species of the family Bovidae: Kobus defassa (KDEXr sequence), Bos taurus (BTAXr sequence) and Antilope cervicapra (ACEXr sequence). The copy numbers of the isolated sequences were assessed using qPCR, and their chromosomal localisations were analysed using FISH in ten bovid tribes and in outgroup species. Besides their localisation on the X chromosome, their presence was also revealed on the Y chromosome and autosomes in several species. The KDEXr sequence abundant in most Bovidae species also occurs in distant taxa (Perissodactyla and Carnivora) and seems to be evolutionarily older than BTAXr and ACEXr. The ACEXr sequence, visible only in several Antilopini species using FISH, is probably the youngest, and arised in an ancestor common to Bovidae and Cervidae. All three repetitive sequences analysed in this study are interspersed among gene-rich regions on the X chromosomes, apparently preventing the crossing-over in their close vicinity. This study demonstrates that repetitive sequences on the X chromosomes have undergone a fast evolution, and their variation among related species can be beneficial for evolutionary studies.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

A super-pangenome of the North American wild grape species

BACKGROUND

Capturing the genetic diversity of wild relatives is crucial for improving crops because wild species are valuable sources of agronomic traits that are essential to enhance the sustainability and adaptability of domesticated cultivars. Genetic diversity across a genus can be captured in super-pangenomes, which provide a framework for interpreting genomic variations.

RESULTS

Here we report the sequencing, assembly, and annotation of nine wild North American grape genomes, which are phased and scaffolded at chromosome scale. We generate a reference-unbiased super-pangenome using pairwise whole-genome alignment methods, revealing the extent of the genomic diversity among wild grape species from sequence to gene level. The pangenome graph captures genomic variation between haplotypes within a species and across the different species, and it accurately assesses the similarity of hybrids to their parents. The species selected to build the pangenome are a great representation of the genus, as illustrated by capturing known allelic variants in the sex-determining region and for Pierce's disease resistance loci. Using pangenome-wide association analysis, we demonstrate the utility of the super-pangenome by effectively mapping short reads from genus-wide samples and identifying loci associated with salt tolerance in natural populations of grapes.

CONCLUSIONS

This study highlights how a reference-unbiased super-pangenome can reveal the genetic basis of adaptive traits from wild relatives and accelerate crop breeding research.

Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes.

Social antagonism facilitates supergene expansion in ants

Antagonistic selection has long been considered a major driver of the formation and expansion of sex chromosomes. For example, sexually antagonistic variation on an autosome can select for suppressed recombination between that autosome and the sex chromosome, leading to a neo-sex chromosome. Autosomal supergenes, chromosomal regions containing tightly linked variants affecting the same complex trait, share similarities with sex chromosomes, raising the possibility that sex chromosome evolution models can explain the evolution of genome structure and recombination in other contexts. We tested this premise in a Formica ant species, wherein we identified four supergene haplotypes on chromosome 3 underlying colony social organization and sex ratio. We discovered a novel rearranged supergene variant (9r) on chromosome 9 underlying queen miniaturization. The 9r is in strong linkage disequilibrium with one chromosome 3 haplotype (P2) found in multi-queen (polygyne) colonies. We suggest that queen miniaturization is strongly disfavored in the single-queen (monogyne) background and is thus socially antagonistic. As such, divergent selection experienced by ants living in alternative social "environments" (monogyne and polygyne) may have contributed to the emergence of a genetic polymorphism on chromosome 9 and associated queen-size dimorphism. Consequently, an ancestral polygyne-associated haplotype may have expanded to include the polymorphism on chromosome 9, resulting in a larger region of suppressed recombination spanning two chromosomes. This process is analogous to the formation of neo-sex chromosomes and consistent with models of expanding regions of suppressed recombination. We propose that miniaturized queens, 16%-20% smaller than queens without 9r, could be incipient intraspecific social parasites.

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos)

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.

Testing immediate dosage compensation in Drosophila miranda via irradiation with heavy-ion beams

Many organisms with heteromorphic sex chromosomes possess a mechanism of dosage compensation (DC) in which X-linked genes are upregulated in males to mitigate the dosage imbalance between sexes and between chromosomes. However, how quickly the DC is established during evolution remains unknown. In this study, by irradiating Drosophila miranda male flies, which carry young sex chromosomes (the so-called neo-sex chromosomes), with heavy-ion beams, we induced deletions in the neo-Y chromosome to mimic the condition of Y-chromosome degeneration, in which functional neo-Y-linked genes are nonfunctionalized; furthermore, we tested whether their neo-X-linked gametologs were immediately upregulated. Because the males that received 2-Gy iron-ion beam irradiation exhibited lower fertility, we sequenced the genomes and transcriptomes of six F1 males derived from these males. Our pipeline identified 82 neo-Y-linked genes in which deletions were predicted in the F1 males. Only three of them showed a one-to-one gametologous relationship with the neo-X-linked genes. The candidate deletions in these three genes occurred in UTRs and did not seriously affect their expression levels. These observations indirectly suggest that DC was unlikely to have operated on the neo-X-linked genes immediately after the pseudogenization of their neo-Y-linked gametologs in D. miranda. Therefore, the dosage imbalance caused by deletions in the neo-Y-linked genes without paralogs may not have effectively been compensated, and individuals with such deletions could have exhibited lethality. Future studies on sex chromosomes at different ages will further reveal the relationship between the age of sex chromosomes and the stringency of DC.

MicrosatNavigator: exploring nonrandom distribution and lineage-specificity of microsatellite repeat motifs on vertebrate sex chromosomes across 186 whole genomes

Microsatellites are short tandem DNA repeats, ubiquitous in genomes. They are believed to be under selection pressure, considering their high distribution and abundance beyond chance or random accumulation. However, limited analysis of microsatellites in single taxonomic groups makes it challenging to understand their evolutionary significance across taxonomic boundaries. Despite abundant genomic information, microsatellites have been studied in limited contexts and within a few species, warranting an unbiased examination of their genome-wide distribution in distinct versus closely related-clades. Large-scale comparisons have revealed relevant trends, especially in vertebrates. Here, "MicrosatNavigator", a new tool that allows quick and reliable investigation of perfect microsatellites in DNA sequences, was developed. This tool can identify microsatellites across the entire genome sequences. Using this tool, microsatellite repeat motifs were identified in the genome sequences of 186 vertebrates. A significant positive correlation was noted between the abundance, density, length, and GC bias of microsatellites and specific lineages. The (AC)n motif is the most prevalent in vertebrate genomes, showing distinct patterns in closely related species. Longer microsatellites were observed on sex chromosomes in birds and mammals but not on autosomes. Microsatellites on sex chromosomes of non-fish vertebrates have the lowest GC content, whereas high-GC microsatellites (≥ 50 M% GC) are preferred in bony and cartilaginous fishes. Thus, similar selective forces and mutational processes may constrain GC-rich microsatellites to different clades. These findings should facilitate investigations into the roles of microsatellites in sex chromosome differentiation and provide candidate microsatellites for functional analysis across the vertebrate evolutionary spectrum.

Characterization of the male-specific region containing the candidate sex-determining gene in Amur catfish (Silurus asotus) using third-generation- and pool-sequencing data

Amur catfish (Silurus asotus) is an ecologically and economically important fish species in Asia. Here, we assembled the female and male Amur catfish genomes, with genome sizes of 757.15 and 755.44 Mb, respectively, at the chromosome level using nanopore and Hi-C technologies. Consistent with the known diploid chromosome count, both genomes contained 29 chromosome-size scaffolds covering 98.80 and 98.73 % of the complete haplotypic assembly with scaffold N50 of 28.87 and 27.29 Mb, respectively. The female (n = 40) and male (n = 40) pools were re-sequenced. Comparative analysis of sequencing and re-sequencing data from both sexes confirmed the presence of an XX/XY sex determination system in Amur catfish and revealed Chr5 as the sex chromosome containing an approximately 400 kb Y-specific region (MSY). Gene annotation revealed a male-specific duplicate of amhr2, namely amhr2y, in MSY, which is male-specific in different wild populations and expressed only in the testes. Amur catfish shared partially syntenic MSY and amhr2y genes with the southern catfish (S. meridionalis, Chr24), which were located on different chromosomes. High sequence divergence between amhr2y and amhr2 and high sequence similarity with amhr2y were observed in both species. These results indicate the common origin of the sex-determining (SD) gene and transition of amhr2y in the two Silurus species. Accumulation of repetitive elements in the MSY of both species may be the main driver of the transition of amhr2y. Overall, our study provides valuable catfish genomic resources. Moreover, determination of amhr2y as the candidate SD gene in Amur catfish provides another example of amhr2 as the SD gene in fish.

Integrating Genomic and Chromosomal Data: A Cytogenetic Study of Transancistrus santarosensis (Loricariidae: Hypostominae) with Characterization of a ZZ/ZW Sex Chromosome System

The plecos (Loricariidae) fish represent a great model for cytogenetic investigations due to their variety of karyotypes, including diploid and polyploid genomes, and different types of sex chromosomes. In this study we investigate Transancistrus santarosensis a rare loricariid endemic to Ecuador, integrating cytogenetic methods with specimens' molecular identification by mtDNA, to describe the the species karyotype. We aim to verify whether sex chromosomes are cytologically identifiable and if they are associated with the accumulation of repetitive sequences present in other species of the family. The analysis of the karyotype (2n = 54 chromosomes) excludes recent centric fusion and pericentromeric inversion and suggests the presence of a ZZ/ZW sex chromosome system at an early stage of differentiation: the W chromosome is degenerated but is not characterized by the presence of differential sex-specific repetitive DNAs. Data indicate that although T. santarosensis has retained the ancestral diploid number of Loricariidae, it accumulated heterochromatin and shows non-syntenic ribosomal genes localization, chromosomal traits considered apomorphic in the family.

Eighty million years of rapid evolution of the primate Y chromosome

The Y chromosome usually plays a critical role in determining male sex and comprises sequence classes that have experienced unique evolutionary trajectories. Here we generated 19 new primate sex chromosome assemblies, analysed them with 10 existing assemblies and report rapid evolution of the Y chromosome across primates. The pseudoautosomal boundary has shifted at least six times during primate evolution, leading to the formation of a Simiiformes-specific evolutionary stratum and to the independent start of young strata in Catarrhini and Platyrrhini. Different primate lineages experienced different rates of gene loss and structural and chromatin change on their Y chromosomes. Selection on several Y-linked genes has contributed to the evolution of male developmental traits across the primates. Additionally, lineage-specific expansions of ampliconic regions have further increased the diversification of the structure and gene composition of the Y chromosome. Overall, our comprehensive analysis has broadened our knowledge of the evolution of the primate Y chromosome.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

Structural variation of a sex-linked region confers monoecy and implicates GATA15 as a master regulator of sex in Salix purpurea

The Salicaceae, including Populus and Salix, are dioecious perennials that utilize different sex determination systems. This family provides a useful system to better understand the evolution of dioecy and sex chromosomes. Here, a rare monoecious genotype of Salix purpurea, 94003, was self- and cross-pollinated and progeny sex ratios were used to test hypotheses on possible mechanisms of sex determination. To delimit genomic regions associated with monoecious expression, the 94003 genome sequence was assembled and DNA- and RNA-Seq of progeny inflorescences was performed. Based on alignments of progeny shotgun DNA sequences to the haplotype-resolved monoecious 94003 genome assembly and reference male and female genomes, a 1.15 Mb sex-linked region on Chr15W was confirmed to be absent in monecious plants. Inheritance of this structural variation is responsible for the loss of a male-suppressing function in what would otherwise be genetic females (ZW), resulting in monoecy (ZW^H or WW^H ), or lethality, if homozygous (W^H W^H ). We present a refined, two-gene sex determination model for Salix purpurea, mediated by ARR17 and GATA15 that is different from the single-gene ARR17-mediated system in the related genus Populus.

Giant Sex Chromosomes in Omophoita Species (Coleoptera, Chrysomelidae): Structural and Evolutionary Relationships Revealed by Zoo-FISH and Comparative Genomic Hybridization (CGH)

The beetles of the subtribe Oedionychina (Chrysomelidae, Alticinae) are the only ones that have the atypical giant and achiasmatic sex chromosomes, which are substantially larger than the autosomes. Previous cytogenetic analyses suggest a large accumulation of repetitive DNA in the sex chromosomes. In this study, we examined the similarity of X and Y chromosomes in four Omophoita species and compared genomic differentiation to better understand the evolutionary process and the giant sex chromosomes origin. Intraspecific genomic comparation using male and female genomes of O. octoguttata and interespecific analyses using genomic DNA of O. octoguttata, O. sexnotata, O. magniguttis, and O. personata were performed. In addition, whole chromosome painting (WCP) experiments were performed with X and Y chromosome probes of O. octogutatta. CGH analysis revealed great genomic similarity between the sexes and a sex-specific region on the Y chromosome, and interspecific analysis revealed a genomic divergence between species. In contrast, WCP results revealed that the sex chromosomes of O. octoguttata have high intra- and interspecific similarity with the studied species. Our data support a common origin under the canonical evolution of the sex chromosomes in this group, as they have high genomic similarity between them.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

Recurrent neo-sex chromosome evolution in kiwifruit

Sex chromosome evolution is thought to be tightly associated with the acquisition and maintenance of sexual dimorphisms. Plant sex chromosomes have evolved independently in many lineages^1,2 and can provide a powerful comparative framework to study this. We assembled and annotated genome sequences of three kiwifruit species (genus Actinidia) and uncovered recurrent sex chromosome turnovers in multiple lineages. Specifically, we observed structural evolution of the neo-Y chromosomes, which was driven via rapid bursts of transposable element insertions. Surprisingly, sexual dimorphisms were conserved in the different species studied, despite the fact that the partially sex-linked genes differ between them. Using gene editing in kiwifruit, we demonstrated that one of the two Y-chromosome-encoded sex-determining genes, Shy Girl, shows pleiotropic effects that can explain the conserved sexual dimorphisms. These plant sex chromosomes therefore maintain sexual dimorphisms through the conservation of a single gene, without a process involving interactions between separate sex-determining genes and genes for sexually dimorphic traits.

Occurrence of Sex Chromosomes in Fish of the Genus Ancistrus with a New Description of Multiple Sex Chromosomes in the Ecuadorian Endemic Ancistrus clementinae (Loricariidae)

Ancistrus Kner, 1854, is the most diverse genus among the Ancistrini (Loricariidae) with 70 valid species showing a wide geographic distribution and great taxonomic and systematic complexity. To date, about 40 Ancistrus taxa have been karyotyped, all from Brazil and Argentina, but the statistic is uncertain because 30 of these reports deal with samples that have not yet been identified at the species level. This study provides the first cytogenetic description of the bristlenose catfish, Ancistrus clementinae Rendahl, 1937, a species endemic to Ecuador, aiming to verify whether a sex chromosome system is identifiable in the species and, if so, which, and if its differentiation is associated with the presence of repetitive sequences reported for other species of the family. We associated the karyotype analysis with the COI molecular identification of the specimens. Karyotype analysis suggested the presence of a ♂ZZ/♀ZW₁W₂ sex chromosome system, never detected before in Ancistrus, with both W₁W₂ chromosomes enriched with heterochromatic blocks and 18S rDNA, in addition to GC-rich repeats (W₂). No differences were observed between males and females in the distribution of 5S rDNA or telomeric repeats. Cytogenetic data here obtained confirm the huge karyotype diversity of Ancistrus, both in chromosome number and sex-determination systems.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Genomic basis of Y-linked dwarfism in cichlids pursuing alternative reproductive tactics

Sexually antagonistic selection, which favours different optima in males and females, is predicted to play an important role in the evolution of sex chromosomes. Body size is a sexually antagonistic trait in the shell-brooding cichlid fish Lamprologous callipterus, as "bourgeois" males must be large enough to carry empty snail shells to build nests whereas females must be small enough to fit into shells for breeding. In this species, there is also a second male morph: smaller "dwarf" males employ an alternative reproductive strategy by wriggling past spawning females into shells to fertilize eggs. L. callipterus male morphology is passed strictly from father to son, suggesting Y-linkage. However, sex chromosomes had not been previously identified in this species, and the genomic basis of size dimorphism was unknown. Here we used whole-genome sequencing to identify a 2.4-Mb sex-linked region on scaffold_23 with reduced coverage and single nucleotide polymorphism density in both male morphs compared to females. Within this sex region, distinct Y-haplotypes delineate the two male morphs, and candidate genes for body size (GHRHR, a known dwarfism gene) and sex determination (ADCYAP1R1) are in high linkage disequilibrium. Because differences in body size between females and males are under strong selection in L. callipterus, we hypothesize that sexual antagonism over body size initiated early events in sex chromosome evolution, followed by Y divergence to give rise to bourgeois and dwarf male reproductive strategies. Our results are consistent with the hypothesis that sexually antagonistic traits should be linked to young sex chromosomes.

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.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

Sex chromosomes in the tribe Cyprichromini (Teleostei: Cichlidae) of Lake Tanganyika

Sex determining loci have been described on at least 12 of 22 chromosomes in East African cichlid fishes, indicating a high rate of sex chromosome turnover. To better understand the rates and patterns of sex chromosome replacement, we used new methods to characterize the sex chromosomes of the cichlid tribe Cyprichromini from Lake Tanganyika. Our k-mer based methods successfully identified sex-linked polymorphisms without the need for a reference genome. We confirm the three previously reported sex chromosomes in this group. We determined the polarity of the sex chromosome turnover on LG05 in Cyprichromis as ZW to XY. We identified a new ZW locus on LG04 in Paracyprichromis brieni. The LG15 XY locus in Paracyprichromis nigripinnis was not found in other Paracyprichromis species, and the sample of Paracyprichromis sp. "tembwe" is likely to be of hybrid origin. Although highly divergent sex chromosomes are thought to develop in a stepwise manner, we show two cases (LG05-ZW and LG05-XY) in which the region of differentiation encompasses most of the chromosome, but appears to have arisen in a single step. This study expands our understanding of sex chromosome evolution in the Cyprichromini, and indicates an even higher level of sex chromosome turnover than previously thought.

Evolutionary dynamics of pseudoautosomal region 1 in humans and great apes

BACKGROUND

The pseudoautosomal region 1 (PAR1) is a 2.7 Mb telomeric region of human sex chromosomes. PAR1 has a crucial role in ensuring proper segregation of sex chromosomes during male meiosis, exposing it to extreme recombination and mutation processes. We investigate PAR1 evolution using population genomic datasets of extant humans, eight populations of great apes, and two archaic human genome sequences.

RESULTS

We find that PAR1 is fast evolving and closer to evolutionary nucleotide equilibrium than autosomal telomeres. We detect a difference between substitution patterns and extant diversity in PAR1, mainly driven by the conflict between strong mutation and recombination-associated fixation bias at CpG sites. We detect excess C-to-G mutations in PAR1 of all great apes, specific to the mutagenic effect of male recombination. Despite recent evidence for Y chromosome introgression from humans into Neanderthals, we find that the Neanderthal PAR1 retained similarity to the Denisovan sequence. We find differences between substitution spectra of these archaics suggesting rapid evolution of PAR1 in recent hominin history. Frequency analysis of alleles segregating in females and males provided no evidence for recent sexual antagonism in this region. We study repeat content and double-strand break hotspot regions in PAR1 and find that they may play roles in ensuring the obligate X-Y recombination event during male meiosis.

CONCLUSIONS

Our study provides an unprecedented quantification of population genetic forces governing PAR1 biology across extant and extinct hominids. PAR1 evolutionary dynamics are predominantly governed by recombination processes with a strong impact on mutation patterns across all species.

ZW sex-chromosome evolution and contagious parthenogenesis in Artemia brine shrimp

Eurasian brine shrimp (genus Artemia) have closely related sexual and asexual lineages of parthenogenetic females, which produce rare males at low frequencies. Although they are known to have ZW chromosomes, these are not well characterized, and it is unclear whether they are shared across the clade. Furthermore, the underlying genetic architecture of the transmission of asexuality, which can occur when rare males mate with closely related sexual females, is not well understood. We produced a chromosome-level assembly for the sexual Eurasian species Artemia sinica and characterized in detail the pair of sex chromosomes of this species. We combined this new assembly with short-read genomic data for the sexual species Artemia sp. Kazakhstan and several asexual lineages of Artemia parthenogenetica, allowing us to perform an in-depth characterization of sex-chromosome evolution across the genus. We identified a small differentiated region of the ZW pair that is shared by all sexual and asexual lineages, supporting the shared ancestry of the sex chromosomes. We also inferred that recombination suppression has spread to larger sections of the chromosome independently in the American and Eurasian lineages. Finally, we took advantage of a rare male, which we backcrossed to sexual females, to explore the genetic basis of asexuality. Our results suggest that parthenogenesis is likely partly controlled by a locus on the Z chromosome, highlighting the interplay between sex determination and asexuality.

Repeated turnovers keep sex chromosomes young in willows

Background

Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively.

Results

Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the “deleterious mutation load” and “sexually antagonistic selection” theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration.

Conclusions

Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02769-w.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops

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

The Rattlesnake W Chromosome: A GC-Rich Retroelement Refugium with Retained Gene Function Across Ancient Evolutionary Strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

X chromosome agents of sexual differentiation

Understanding sex differences in physiology and disease requires the identification of the molecular agents that cause phenotypic sex differences. Two groups of such agents are genes located on the sex chromosomes, and gonadal hormones. The former have coherent linkage to chromosomes that form differently in the two sexes under the influence of genomic forces that are not related to reproductive function, whereas the latter have a direct or indirect relationship to reproduction. Evidence published in the past 5 years supports the identification of several agents of sexual differentiation encoded by the X chromosome in mice, including Kdm5c, Kdm6a, Ogt and Xist. These X chromosome agents have wide pleiotropic effects, potentially influencing sex differences in many different tissues, a characteristic shared with the gonadal hormones. The identification of X chromosome agents of sexual differentiation will facilitate understanding of complex intersecting gene pathways underlying sex differences in disease.

Differing associations between sex determination and sex‐linked inversions in two ecotypes of Littorina saxatilis

Sexual antagonism is a common hypothesis for driving the evolution of sex chromosomes, whereby recombination suppression is favored between sexually antagonistic loci and the sex‐determining locus to maintain beneficial combinations of alleles. This results in the formation of a sex‐determining region. Chromosomal inversions may contribute to recombination suppression but their precise role in sex chromosome evolution remains unclear. Because local adaptation is frequently facilitated through the suppression of recombination between adaptive loci by chromosomal inversions, there is potential for inversions that cover sex‐determining regions to be involved in local adaptation as well, particularly if habitat variation creates environment‐dependent sexual antagonism. With these processes in mind, we investigated sex determination in a well‐studied example of local adaptation within a species: the intertidal snail, Littorina saxatilis. Using SNP data from a Swedish hybrid zone, we find novel evidence for a female‐heterogametic sex determination system that is restricted to one ecotype. Our results suggest that four putative chromosomal inversions, two previously described and two newly discovered, span the putative sex chromosome pair. We determine their differing associations with sex, which suggest distinct strata of differing ages. The same inversions are found in the second ecotype but do not show any sex association. The striking disparity in inversion‐sex associations between ecotypes that are connected by gene flow across a habitat transition that is just a few meters wide indicates a difference in selective regime that has produced a distinct barrier to the spread of the newly discovered sex‐determining region between ecotypes. Such sex chromosome‐environment interactions have not previously been uncovered in L. saxatilis and are known in few other organisms. A combination of both sex‐specific selection and divergent natural selection is required to explain these highly unusual patterns.

Ecology and the evolution of sex chromosomes

Sex chromosomes are common features of animal genomes, often carrying a sex determination gene responsible for initiating the development of sexually dimorphic traits. The specific chromosome that serves as the sex chromosome differs across taxa as a result of fusions between sex chromosomes and autosomes, along with sex chromosome turnover-autosomes becoming sex chromosomes and sex chromosomes 'reverting' back to autosomes. In addition, the types of genes on sex chromosomes frequently differ from the autosomes, and genes on sex chromosomes often evolve faster than autosomal genes. Sex-specific selection pressures, such as sexual antagonism and sexual selection, are hypothesized to be responsible for sex chromosome turnovers, the unique gene content of sex chromosomes and the accelerated evolutionary rates of genes on sex chromosomes. Sex-specific selection has pronounced effects on sex chromosomes because their sex-biased inheritance can tilt the balance of selection in favour of one sex. Despite the general consensus that sex-specific selection affects sex chromosome evolution, most population genetic models are agnostic as to the specific sources of these sex-specific selection pressures, and many of the details about the effects of sex-specific selection remain unresolved. Here, I review the evidence that ecological factors, including variable selection across heterogeneous environments and conflicts between sexual and natural selection, can be important determinants of sex-specific selection pressures that shape sex chromosome evolution. I also explain how studying the ecology of sex chromosome evolution can help us understand important and unresolved aspects of both sex chromosome evolution and sex-specific selection.

Stepwise evolution of a butterfly supergene via duplication and inversion

Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning coloration in the African monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of more than 1 million bp underwent several segmental duplications at least 7.5 Ma. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree and suggest that structural polymorphisms have persisted for at least 4.1 Myr. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes.

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

Remarkable Divergence of the Sex-Linked Region between Two Wild Spinach Progenitors, Spinacia turkestanica and Spinacia tetrandra

The sex-linked region (SLR) plays an important role in determining the sex of a plant. The SLR of the Y chromosome, composed of a 14.1-Mb inversion and a 10-Mb Y-duplication region (YDR), was deciphered in Spinacia oleracea previously. However, our understanding of the SLR in its wild relatives, S. turkestanica and S. tetrandra, remains limited. In this study, we used 63 resequencing data from the three Spinacia species to infer the evolution of the SLR among the Spinacia species. In the SLR, all the cultivated spinach and S. turkestanica accessions were clustered into two distinct categories with both sexes, while the S. tetrandra accessions of both sexes were grouped. This suggests that S. oleracea shared a similar SLR with S. turkestanica, but not with S. tetrandra, which was further confirmed based on the population structure and principal component analysis. Furthermore, we identified 3910 fully sex-linked SNPs in S. oleracea and 92.82% of them were available in S. turkestanica, while none of the SNPs were adopted in S. tetrandra. Genome coverage in males and females supported the hypothesis that the YDR increasingly expanded during its evolution. Otherwise, we identified 13 sex-linked transposable element insertion polymorphisms within the inversion in both S. oleracea and S. turkestanica, demonstrating that the transposable element insertions might have occurred before the recombination suppression event of the inversion. The SLR was conserved compared with the pseudoautosomal region given that the genetic hitchhiking process occurred in the SLR during its evolution. Our findings will significantly advance our understanding of the characteristics and evolution of the SLR in Spinacia species.

Chromosome-scale assembly with a phased sex-determining region resolves features of early Z and W chromosome differentiation in a wild octoploid strawberry

When sex chromosomes stop recombining, they start to accumulate differences. The sex-limited chromosome (Y or W) especially is expected to degenerate via the loss of nucleotide sequence and the accumulation of repetitive sequences. However, how early signs of degeneration can be detected in a new sex chromosome is still unclear. The sex-determining region of the octoploid strawberries is young, small, and dynamic. Using PacBio HiFi reads, we obtained a chromosome-scale assembly of a female (ZW) Fragaria chiloensis plant carrying the youngest and largest of the known sex-determining region on the W in strawberries. We fully characterized the previously incomplete sex-determining region, confirming its gene content, genomic location, and evolutionary history. Resolution of gaps in the previous characterization of the sex-determining region added 10 kb of sequence including a noncanonical long terminal repeat-retrotransposon; whereas the Z sequence revealed a Harbinger transposable element adjoining the sex-determining region insertion site. Limited genetic differentiation of the sex chromosomes coupled with structural variation may indicate an early stage of W degeneration. The sex chromosomes have a similar percentage of repeats but differ in their repeat distribution. Differences in the pattern of repeats (transposable element polymorphism) apparently precede sex chromosome differentiation, thus potentially contributing to recombination cessation as opposed to being a consequence of it.