Complex evolutionary trajectories of sex chromosomes across bird taxa

Functional innovation through new genes as a general evolutionary process

In the past decade, our understanding of how new genes originate in diverse organisms has advanced substantially, and more than a dozen molecular mechanisms for generating initial gene structures were identified, in addition to gene duplication. These new genes have been found to integrate into and modify pre-existing gene networks primarily through mutation and selection, revealing new patterns and rules with stable origination rates across various organisms. This progress has challenged the prevailing belief that new proteins evolve from pre-existing genes, as new genes may arise de novo from noncoding DNA sequences in many organisms, with high rates observed in flowering plants. New genes have important roles in phenotypic and functional evolution across diverse biological processes and structures, with detectable fitness effects of sexual conflict genes that can shape species divergence. Such knowledge of new genes can be of translational value in agriculture and medicine.

Hierarchical architecture of neo-sex chromosomes and accelerated adaptive evolution in tortricid moths

Sex chromosomes can expand through fusion with autosomes, thereby acquiring unique evolutionary patterns. In butterflies and moths (Lepidoptera), these sex chromosome-autosome (SA) fusions occur relatively frequently, suggesting possible evolutionary advantages. Here, we investigated how SA fusion affects chromosome features and molecular evolution in leafroller moths (Lepidoptera: Tortricidae). Phylogenomic analysis showed that Tortricidae diverged ∼124 million years ago, accompanied by an SA fusion between the Merian elements M(20 + 17) and MZ. In contrast to partial autosomal fusions, the fused neo-Z Chromosome developed a hierarchical architecture, in which the three elements exhibit heterogeneous sequence features and evolutionary patterns. Specifically, the M17 part had a distinct base composition and chromatin domains. Unlike M20 and MZ, M17 was expressed at the same levels as autosomes in both sexes, compensating for the lost gene dosage in females. Concurrently, the SA fusion drove M17 as an evolutionary hotspot, accelerating the evolution of several genes related to ecological adaptation (e.g., ABCCs) and facilitating the divergence of closely related species, whereas the undercompensated M20 did not show such an effect. Thus, accelerated evolution under a novel pattern of dosage compensation may have favored the adaptive radiation of this group. This study demonstrates the association between a karyotype variant and adaptive evolution and explains the recurrent SA fusion in the Lepidoptera.

Application of Minimally Invasive Oral Swab Samples for qPCR-Based Sexing in Neognathae Birds

Birds are inherently social creatures that rely on pairing to enhance their well-being. Since many bird species lack obvious physical differences between females and males, sex identification is essential for ensuring their welfare. Additionally, early determination of the sexes of birds is crucial for their breeders, especially considering that most companion birds do not display clear sexual characteristics. Molecular genetic sexing has been demonstrated to be the most reliable method for determining the sexes of monomorphic birds. The objective of the present study was to demonstrate rapid, effective, and precise identification of sex in birds through quantitative real-time PCR (qPCR) using samples obtained via a minimally invasive technique (oral swabs). This qPCR method assesses variations in gene copy numbers within conserved Z-specific genes such as CHRNA6, DDX4, VPS13A, LPAR1, and TMEM161B, which are absent from the W chromosome. A total of 34 samples were included in this study from the following 17 bird species: domestic pigeon (Columba livia domestica), domestic chicken (Gallus gallus domesticus), domestic goose (Anser anser f domesticus), domestic duck (Anas platyrhynchos domesticus), Mute swan (Cygnus olor), Budgerigar (Melopsittacus undulatus), Lovebird (Agapornis roseicollis), Cockatiel (Nymphicus hollandicus), Red-rumped parrot (Psephotus haematonotus), Rose-ringed parakeet (Psittacula krameri), African grey parrot (Psittacus erithacus), domestic Canary (Serinus canaria forma domestica), Goldfinch (Carduelis carduelis major), Gouldian Finch (Chloebia gouldiae), Red Siskin (Carduelis cucullata), Australian Zebra Finch (Taeniopygia castanotis), and Common buzzard (Buteo buteo). The results proved that the CHRNA6, DDX4, VPS13A, LPAR1, and TMEM161B genes can reveal the sexes in the Neognath birds tested.

Transposable elements shape the landscape of heterozygous structural variation in a bird genome

Avian genomes exhibit compact organization and remarkable chromosomal stability. However, the extent and mechanisms by which structural variation in avian genomes differ from those in other vertebrate lineages are poorly explored. This study generated a diploid genome assembly for the golden pheasant ( Chrysolophus pictus), a species distinguished by the vibrant plumage of males. Each haploid genome assembly included complete chromosomal models, incorporating all microchromosomes. Analysis revealed extensive tandem amplification of immune-related genes across the smallest microchromosomes (dot chromosomes), with an average copy number of 54. Structural variation between the haploid genomes was primarily shaped by large insertions and deletions (indels), with minimal contributions from inversions or duplications. Approximately 28% of these large indels were associated with recent insertions of transposable elements, despite their typically low activity in bird genomes. Evidence for significant effects of transposable elements on gene expression was minimal. Evolutionary strata on the sex chromosomes were identified, along with a drastic rearrangement of the W chromosome. These analyses of the high-quality diploid genome of the golden pheasant provide valuable insights into the evolutionary patterns of structural variation in avian genomes.

Chromosome-level echidna genome illuminates evolution of multiple sex chromosome system in monotremes

BACKGROUND

A thorough analysis of genome evolution is fundamental for biodiversity understanding. The iconic monotremes (platypus and echidna) feature extraordinary biology. However, they also exhibit rearrangements in several chromosomes, especially in the sex chromosome chain. Therefore, the lack of a chromosome-level echidna genome has limited insights into genome evolution in monotremes, in particular the multiple sex chromosomes complex.

RESULTS

Here, we present a new long reads-based chromosome-level short-beaked echidna (Tachyglossus aculeatus) genome, which allowed the inference of chromosomal rearrangements in the monotreme ancestor (2n = 64) and each extant species. Analysis of the more complete sex chromosomes uncovered homology between 1 Y chromosome and multiple X chromosomes, suggesting that it is the ancestral X that has undergone reciprocal translocation with ancestral autosomes to form the complex. We also identified dozens of ampliconic genes on the sex chromosomes, with several ancestral ones expressed during male meiosis, suggesting selective constraints in pairing the multiple sex chromosomes.

CONCLUSION

The new echidna genome provides an important basis for further study of the unique biology and conservation of this species.

Sex Chromosome Turnovers and Stability in Snakes

For a long time, snakes were presented as a textbook example of a group with gradual differentiation of homologous ZZ/ZW sex chromosomes. However, recent advances revealed that the ZZ/ZW sex chromosomes characterize only caenophidian snakes and certain species of boas and pythons have nonhomologous XX/XY sex chromosomes. We used genome coverage analysis in four non-caenophidian species to identify their sex chromosomes, and we examined the homology of sex chromosomes across phylogenetically informative snake lineages. We identified sex chromosomes for the first time in 13 species of non-caenophidian snakes, providing much deeper insights into the evolutionary history of snake sex chromosomes. The evolution of sex chromosomes in snakes is more complex than previously thought. Snakes may have had ancestral XX/XY sex chromosomes, which are still present in a blind snake and some boas, and there were several transitions to derived XX/XY sex chromosomes with different gene content and two or even three transitions to ZZ/ZW sex chromosomes. However, we discuss more alternative scenarios. In any case, we document that (1) some genomic regions were likely repeatedly co-opted as sex chromosomes in phylogenetically distant lineages, even with opposite types of heterogamety; (2) snake lineages differ greatly in the rate of differentiation of sex chromosomes; (3) snakes likely originally possessed sex chromosomes prone to turnovers. The sex chromosomes became evolutionarily highly stable once their differentiation progressed in the megadiverse caenophidian snakes. Snakes thus provide an ideal system for studying the evolutionary factors that drive unequal rates of differentiation, turnovers and stability of sex chromosomes.

Karyotype evolution of suliformes and description of a ♂Z1Z1Z2Z2/♀Z1Z2W multiple sex chromosome system in boobies (Sula spp.)

Our comprehension of avian karyotypes still needs to be improved, especially for Suliform birds. To enhance understanding of chromosomal evolution in this order, we conducted conventional and molecular cytogenetic analysis in five species, named Sula dactylatra, Sula leucogaster, Sula sula (Sulidae), Fregata magnificens (Fregatidae), and Nannopterum brasilianum (Phalacrocoracidae). The diploid chromosome number for S. dactylatra and S. leucogaster was established as 2n = 76 in males, and 2n = 75 in females, but S. sula displayed a karyotype of 2n = 76 chromosomes in males. The disparity in diploid chromosome numbers between male and female Sula is due to a multiple sex chromosome system of the Z1Z1Z2Z2/Z1Z2W type. We propose that the emergence of this multiple-sex chromosome system resulted from a Robertsonian translocation involving the W chromosome and the smallest microchromosome. Fregata magnificens exhibited a diploid number 76 (2n = 76), while N. brasilianum displayed a diploid number of 74 (2n = 74) in both sexes. The ribosomal cluster was located in one microchromosome pair in S. dactylatra, S. leucogaster, S. sula, and F. magnificens and in four pairs in N. brasilianum. Our findings provide evidence of a conserved multiple-sex chromosome system within the Sula genus, shedding light on the high karyotype diversity in Suliformes.

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.

Orthopteran Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and Evolution of Supergenes

The early evolution of sex chromosomes has remained obscure for more than a century. The Vandiemenella viatica species group of morabine grasshoppers is highly suited for studying the early stages of sex chromosome divergence and degeneration of the Y chromosome. This stems from the fact that neo-XY sex chromosomes have independently evolved multiple times by X-autosome fusions with different autosomes. Here, we generated new chromosome-level assemblies for two chromosomal races representing karyotypes with and without neo-sex chromosomes (P24XY and P24X0), and sequence data of a third chromosomal race with a different neo-XY chromosome system (P25XY). Interestingly, these two neo-XY chromosomal races are formed by different X-autosome fusions (involving chr1 and chrB, respectively), and we found that both neo-Y chromosomes have partly ceased to recombine with their neo-X counterpart. We show that the neo-XY chromosomes have diverged through accumulation of SNPs and structural mutations, and that many neo-Y-linked genes have degenerated since recombination ceased. However, the non-recombining regions of neo-Y chromosomes host non-degenerated genes crucial for sex determination, such as sex-lethal and transformer, alongside genes associated with spermatogenesis, fertility, and reproduction, illustrating their integrative role as a masculinizing supergene. Contrary to expectations, the neo-Y chromosomes showed (slightly) lower density of transposable elements (TEs) compared to other genomic regions. The study reveals the unique dynamics of young sex chromosomes, with evolution of recombination suppression and pronounced decay of (some) neo-sex chromosome genes, and provides a compelling case illustrating how chromosomal fusions and post-fusion mutational processes contribute to the evolution of supergenes.

Genome assemblies for Chromidotilapia guntheri (Teleostei: Cichlidae) identify a novel candidate gene for vertebrate sex determination, RIN3

Advances in genome sequencing have greatly accelerated the identification of sex chromosomes in a variety of species. Many of these species have experienced structural rearrangements that reduce recombination between the sex chromosomes, allowing the accumulation of sequence differences over many megabases. Identification of the genes that are responsible for sex determination within these sometimes large regions has proved difficult. Here, we identify an XY sex chromosome system on LG19 in the West African cichlid fish Chromidotilapia guntheri in which the region of differentiation extends over less than 400 kb. We develop high-quality male and female genome assemblies for this species, which confirm the absence of structural variants, and which facilitate the annotation of genes in the region. The peak of differentiation lies within rin3, which has experienced several debilitating mutations on the Y chromosome. We suggest two hypotheses about how these mutations might disrupt endocytosis, leading to Mendelian effects on sexual development.

Genomic Anatomy of Homozygous XX Females and YY Males Reveals Early Evolutionary Trajectory of Sex-determining Gene and Sex Chromosomes in Silurus Fishes

Sex chromosomes display remarkable diversity and variability among vertebrates. Compared with research on the X/Y and Z/W chromosomes, which have long evolutionary histories in mammals and birds, studies on the sex chromosomes at early evolutionary stages are limited. Here, we precisely assembled the genomes of homozygous XX female and YY male Lanzhou catfish (Silurus lanzhouensis) derived from an artificial gynogenetic family and a self-fertilized family, respectively. Chromosome 24 (Chr24) was identified as the sex chromosome based on resequencing data. Comparative analysis of the X and Y chromosomes showed an approximate 320 kb Y-specific region with a Y-specific duplicate of anti-Mullerian hormone type II receptor (amhr2y), which is consistent with findings in 2 other Silurus species but on different chromosomes (Chr24 of Silurus meridionalis and Chr5 of Silurus asotus). Deficiency of amhr2y resulted in male-to-female sex reversal, indicating that amhr2y plays a male-determining role in S. lanzhouensis. Phylogenetic analysis and comparative genomics revealed that the common sex-determining gene amhr2y was initially translocated to Chr24 of the Silurus ancestor along with the expansion of transposable elements. Chr24 was maintained as the sex chromosome in S. meridionalis and S. lanzhouensis, whereas a sex-determining region transition triggered sex chromosome turnover from Chr24 to Chr5 in S. asotus. Additionally, gene duplication, translocation, and degeneration were observed in the Y-specific regions of Silurus species. These findings present a clear case for the early evolutionary trajectory of sex chromosomes, including sex-determining gene origin, repeat sequence expansion, gene gathering and degeneration in sex-determining region, and sex chromosome turnover.

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.

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

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

Why do sex chromosomes progressively lose recombination?

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

Y-chromosome Degeneration due to Speciation and Founder Effect

The Y chromosome in the XY sex-determination system is often shorter than its X counterpart, a condition attributed to degeneration after Y recombination ceases. Contrary to the traditional view of continuous, gradual degeneration, our study reveals stabilization within large mating populations. In these populations, we demonstrate that both mutant and active alleles on the Y chromosome can reach equilibrium through a mutation-selection balance. However, the emergence of a new species, particularly through the founder effect, can disrupt this equilibrium. Specifically, if the male founders of a new species carry only a mutant allele for a particular Y-linked gene, this allele becomes fixed, leading to the loss of the corresponding active gene on the Y chromosome. Our findings suggest that the rate of Y-chromosome degeneration may be linked to the frequency of speciation events associated with single-male founder events.

Molecular sexing of birds using quantitative PCR (qPCR) of sex-linked genes and logistic regression models

The ability to sex individuals is an important component of many behavioural and ecological investigations and provides information for demographic models used in conservation and species management. However, many birds are difficult to sex using morphological characters or traditional molecular sexing methods. In this study, we developed probabilistic models for sexing birds using quantitative PCR (qPCR) data. First, we quantified distributions of gene copy numbers at a set of six sex-linked genes, including the sex-determining gene DMRT1, for individuals across 17 species and seven orders of birds (n = 150). Using these data, we built predictive logistic models for sex identification and tested their performance with independent samples from 51 species and 13 orders (n = 209). Models using the two loci most highly correlated with sex had greater accuracy than models using the full set of sex-linked loci, across all taxonomic levels of analysis. Sex identification was highly accurate when individuals to be assigned were of species used in model building. Our analytical approach was widely applicable across diverse neognath bird lineages spanning millions of years of evolutionary divergence. Unlike previous methods, our probabilistic framework incorporates uncertainty around qPCR measurements as well as biological variation within species into decision-making rules. We anticipate that this method will be useful for sexing birds, including those of high conservation concern and/or subsistence value, that have proven difficult to sex using traditional approaches. Additionally, the general analytical framework presented in this paper may also be applicable to other organisms with sex chromosomes.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Evolution and expression patterns of the neo-sex chromosomes of the crested ibis

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

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.

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

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.

A high-quality genome assembly highlights the evolutionary history of the great bustard (Otis tarda, Otidiformes)

Conservation genomics often relies on non-invasive methods to obtain DNA fragments which limit the power of multi-omic analyses for threatened species. Here, we report multi-omic analyses based on a well-preserved great bustard individual (Otis tarda, Otidiformes) that was found dead in the mountainous region in Gansu, China. We generate a near-complete genome assembly containing only 18 gaps scattering in 8 out of the 40 assembled chromosomes. We characterize the DNA methylation landscape which is correlated with GC content and gene expression. Our phylogenomic analysis suggests Otidiformes and Musophagiformes are sister groups that diverged from each other 46.3 million years ago. The genetic diversity of great bustard is found the lowest among the four available Otidiformes genomes, possibly due to population declines during past glacial periods. As one of the heaviest migratory birds, great bustard possesses several expanded gene families related to cardiac contraction, actin contraction, calcium ion signaling transduction, as well as positively selected genes enriched for metabolism. Finally, we identify an extremely young evolutionary stratum on the sex chromosome, a rare case among birds. Together, our study provides insights into the conservation genomics, adaption and chromosome evolution of the great bustard.

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.

The evolutionary maintenance of ancient recombining sex chromosomes in the ostrich

Sex chromosomes have evolved repeatedly across the tree of life and often exhibit extreme size dimorphism due to genetic degeneration of the sex-limited chromosome (e.g. the W chromosome of some birds and Y chromosome of mammals). However, in some lineages, ancient sex-limited chromosomes have escaped degeneration. Here, we study the evolutionary maintenance of sex chromosomes in the ostrich (Struthio camelus), where the W remains 65% the size of the Z chromosome, despite being more than 100 million years old. Using genome-wide resequencing data, we show that the population scaled recombination rate of the pseudoautosomal region (PAR) is higher than similar sized autosomes and is correlated with pedigree-based recombination rate in the heterogametic females, but not homogametic males. Genetic variation within the sex-linked region (SLR) (π = 0.001) was significantly lower than in the PAR, consistent with recombination cessation. Conversely, genetic variation across the PAR (π = 0.0016) was similar to that of autosomes and dependent on local recombination rates, GC content and to a lesser extent, gene density. In particular, the region close to the SLR was as genetically diverse as autosomes, likely due to high recombination rates around the PAR boundary restricting genetic linkage with the SLR to only ~50Kb. The potential for alleles with antagonistic fitness effects in males and females to drive chromosome degeneration is therefore limited. While some regions of the PAR had divergent male-female allele frequencies, suggestive of sexually antagonistic alleles, coalescent simulations showed this was broadly consistent with neutral genetic processes. Our results indicate that the degeneration of the large and ancient sex chromosomes of the ostrich may have been slowed by high recombination in the female PAR, reducing the scope for the accumulation of sexually antagonistic variation to generate selection for recombination cessation.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Micro Germline-Restricted Chromosome in Blue Tits: Evidence for Meiotic Functions

The germline-restricted chromosome (GRC) is likely present in all songbird species but differs widely in size and gene content. This extra chromosome has been described as either a microchromosome with only limited basic gene content or a macrochromosome with enriched gene functions related to female gonad and embryo development. Here, we assembled, annotated, and characterized the first micro-GRC in the blue tit (Cyanistes caeruleus) using high-fidelity long-read sequencing data. Although some genes on the blue tit GRC show signals of pseudogenization, others potentially have important functions, either currently or in the past. We highlight the GRC gene paralog BMP15, which is among the highest expressed GRC genes both in blue tits and in zebra finches (Taeniopygia guttata) and is known to play a role in oocyte and follicular maturation in other vertebrates. The GRC genes of the blue tit are further enriched for functions related to the synaptonemal complex. We found a similar functional enrichment when analyzing published data on GRC genes from two nightingale species (Luscinia spp.). We hypothesize that these genes play a role in maintaining standard maternal inheritance or in recombining maternal and paternal GRCs during potential episodes of biparental inheritance.

Heterochiasmy and Sex Chromosome Evolution in Silene

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

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

Genetic sex determination, sex chromosome size and sex-specific lifespans across tetrapods

Sex differences in lifespan are ubiquitous across the tree of life and exhibit broad taxonomic patterns that remain a puzzle, such as males living longer than females in birds and vice versa in mammals. The prevailing unguarded X hypothesis explains sex differences in lifespan by differential expression of recessive mutations on the X or Z chromosome of the heterogametic sex, but has only received indirect support to date. An alternative hypothesis is that the accumulation of deleterious mutations and repetitive elements on the Y or W chromosome might lower the survival of the heterogametic sex ('toxic Y' hypothesis). Here, we use a new database to report lower survival of the heterogametic relative to the homogametic sex across 136 species of birds, mammals, reptiles and amphibians, as expected if sex chromosomes shape sex-specific lifespans, and consistent with previous findings. We also found that the relative sizes of both the X and the Y chromosomes in mammals (but not the Z or the W chromosomes in birds) are associated with sex differences in lifespan, as predicted by the unguarded X and the 'toxic Y'. Furthermore, we report that the relative size of the Y is negatively associated with male lifespan in mammals, so that small Y size correlates with increased male lifespan. In theory, toxic Y effects are expected to be particularly strong in mammals, and we did not find similar effects in birds. Our results confirm the role of sex chromosomes in explaining sex differences in lifespan across tetrapods and further suggest that, at least in mammals, 'toxic Y' effects may play an important part in this role.

Madagascar Leaf-Tail Geckos (Uroplatus spp.) Share Independently Evolved Differentiated ZZ/ZW Sex Chromosomes

Geckos are an excellent group to study the evolution of sex determination, as they possess a remarkable variability ranging from a complete absence of sex chromosomes to highly differentiated sex chromosomes. We explored sex determination in the Madagascar leaf-tail geckos of the genus Uroplatus. The cytogenetic analyses revealed highly heterochromatic W chromosomes in all three examined species (Uroplatus henkeli, U. alluaudi, U. sikorae). The comparative gene coverage analysis between sexes in U. henkeli uncovered an extensive Z-specific region, with a gene content shared with the chicken chromosomes 8, 20, 26 and 28. The genomic region homologous to chicken chromosome 28 has been independently co-opted for the role of sex chromosomes in several vertebrate lineages, including monitors, beaded lizards and monotremes, perhaps because it contains the amh gene, whose homologs were repeatedly recruited as a sex-determining locus. We demonstrate that all tested species of leaf-tail geckos share homologous sex chromosomes despite the differences in shape and size of their W chromosomes, which are not homologous to the sex chromosomes of other closely related genera. The rather old (at least 40 million years), highly differentiated sex chromosomes of Uroplatus geckos can serve as a great system to study the convergence of sex chromosomes evolved from the same genomic region.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

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.

ZZ/ZW Sex Chromosomes in the Madagascar Girdled Lizard, Zonosaurus madagascariensis (Squamata: Gerrhosauridae)

Scincoidea, the reptilian clade that includes girdled lizards, night lizards, plated lizards and skinks, are considered as a lineage with diversity in sex-determining systems. Recently, the hypothesis on the variability in sex determination in skinks and even more the absence of sex chromosomes in some of them has been rivalling. Homologous, evolutionary stable XX/XY sex chromosomes were documented to be widespread across skinks. However, sex determination in the other scincoidean families is highly understudied. ZZ/ZW sex chromosomes have been identified only in night lizards and a single species of plated lizards. It seems that although there are different sex chromosome systems among scincoidean lineages, they share one common trait: they are mostly poorly differentiated and often undetectable by cytogenetic methods. Here, we report one of the exceptions, demonstrating for the first time ZZ/ZW sex chromosomes in the plated lizard Zonosaurus madagascariensis. Its sex chromosomes are morphologically similar, but the W is clearly detectable by comparative genomic hybridization (CGH), suggesting that the Z and W chromosomes highly differ in sequence content. Our findings confirm the presence of female heterogamety in plated lizards and provides novel insights to expand our understanding of sex chromosome evolution in scincoidean lizards.

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.

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.

New insights into Xenopus sex chromosome genomics from the Marsabit clawed frog X. borealis

In many groups, sex chromosomes change frequently but the drivers of their rapid evolution are varied and often poorly characterized. With an aim of further understanding sex chromosome turnover, we investigated the polymorphic sex chromosomes of the Marsabit clawed frog, Xenopus borealis, using genomic data and a new chromosome-scale genome assembly. We confirmed previous findings that 54.1 Mb of chromosome 8L is sex-linked in animals from east Kenya and a laboratory strain, but most (or all) of this region is not sex-linked in natural populations from west Kenya. Previous work suggests possible degeneration of the Z chromosomes in the east population because many sex-linked transcripts of this female heterogametic population have female-biased expression, and we therefore expected this chromosome to not be present in the west population. In contrast, our simulations support a model where most or all of the sex-linked portion of the Z chromosome from the east acquired autosomal segregation in the west, and where much genetic variation specific to the large sex-linked portion of the W chromosome from the east is not present in the west. These recent changes are consistent with the hot-potato model, wherein sex chromosome turnover is favoured by natural selection if it purges a (minimally) degenerate sex-specific sex chromosome, but counterintuitively suggest natural selection failed to purge a Z chromosome that has signs of more advanced and possibly more ancient regulatory degeneration. These findings highlight complex evolutionary dynamics of young, rapidly evolving Xenopus sex chromosomes and set the stage for mechanistic work aimed at pinpointing additional sex-determining genes in this group.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

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.

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.

Draft genome assemblies of four manakins

Manakins are a family of small suboscine passerine birds characterized by their elaborate courtship displays, non-monogamous mating system, and sexual dimorphism. This family has served as a good model for the study of sexual selection. Here we present genome assemblies of four manakin species, including Cryptopipo holochlora, Dixiphia pipra (also known as Pseudopipra pipra), Machaeropterus deliciosus and Masius chrysopterus, generated by Single-tube Long Fragment Read (stLFR) technology. The assembled genome sizes ranged from 1.10 Gb to 1.19 Gb, with average scaffold N50 of 29 Mb and contig N50 of 169 Kb. On average, 12,055 protein-coding genes were annotated in the genomes, and 9.79% of the genomes were annotated as repetitive elements. We further identified 75 Mb of Z-linked sequences in manakins, containing 585 to 751 genes and an ~600 Kb pseudoautosomal region (PAR). One notable finding from these Z-linked sequences is that a possible Z-to-autosome/PAR reversal could have occurred in M. chrysopterus. These de novo genomes will contribute to a deeper understanding of evolutionary history and sexual selection in manakins.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The genome is the complete DNA sequence of an individual. It is a crucial foundation for many studies in medicine, agriculture, and conservation biology. Advances in genetics have made it possible to rapidly sequence, or read out, the genome of many organisms. For closely related species, scientists can then do detailed comparisons, revealing similar genes with a shared past or a common role, but comparing more distantly related organisms remains difficult.

One major challenge is that genes are often lost or duplicated over evolutionary time. One way to be more confident is to look at ‘synteny’, or how genes are organized or ordered within the genome. In some groups of species, synteny persists across millions of years of evolution. Combining sequence similarity with gene order could make comparisons between distantly related species more robust.

To do this, Lovell et al. developed GENESPACE, a software that links similarities between DNA sequences to the order of genes in a genome. This allows researchers to visualize and explore related DNA sequences and determine whether genes have been lost or duplicated. To demonstrate the value of GENESPACE, Lovell et al. explored evolution in vertebrates and flowering plants. The software was able to highlight the shared sequences between unique sex chromosomes in birds and mammals, and it was able to track the positions of genes important in the evolution of grass crops including maize, wheat, and rice.

Exploring the genetic code in this way could lead to a better understanding of the evolution of important sections of the genome. It might also allow scientists to find target genes for applications like crop improvement. Lovell et al. have designed the GENESPACE software to be easy for other scientists to use, allowing them to make graphics and perform analyses with few programming skills.

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.

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.

Chromosome-level genome assembly defines female-biased genes associated with sex determination and differentiation in the human blood fluke Schistosoma japonicum

Schistosomiasis is a neglected tropical disease of humans caused by blood flukes of the genus Schistosoma, the only dioecious parasitic flatworm. Although aspects of sex determination, differentiation and reproduction have been studied in some Schistosoma species, almost nothing is known for Schistosoma japonicum, the causative agent of schistosomiasis japonica. This mainly reflects the lack of high-quality genomic and transcriptomic resources for this species. As current genomes for S. japonicum are highly fragmented, we assembled and report a chromosome-level reference genome (seven autosomes, the Z-chromosome and partial W-chromosome), achieving a substantially enhanced gene annotation. Utilizing this genome, we discovered that the sex chromosomes of S. japonicum and its congener S. mansoni independently suppressed recombination during evolution, forming five and two evolutionary strata, respectively. By exploring the W-chromosome and sex-specific transcriptomes, we identified 35 W-linked genes and 257 female-preferentially transcribed genes (FTGs) from our chromosomal assembly and uncovered a signature for sex determination and differentiation in S. japonicum. These FTGs clustering within autosomes or the Z-chromosome exhibit a highly dynamic transcription profile during the pairing of female and male schistosomula, thereby representing a critical phase for the maturation of the female worms and suggesting distinct layers of regulatory control of gene transcription at this development stage. Collectively, these data provide a valuable resource for further functional genomic characterization of S. japonicum, shed light on the evolution of sex chromosomes in this highly virulent human blood fluke, and provide a pathway to identify novel targets for development of intervention tools against schistosomiasis.

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.