Chromosome Evolution Research Articles

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

Outcome of Allogeneic Hematopoietic Stem Cell Transplantation in a Child with Myelodysplastic Neoplasm with Complex Karyotype and ETV6 Variant.

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is the only treatment with the potential for cure in patients with myelodysplastic neoplasm (MDS). Nevertheless, disease relapse is the main cause of treatment failure. Due to the rarity of childhood MDS (cMDS), there are few studies showing the impact of cytogenetic alterations and genetic variants on outcomes of allo-HSCT, mainly focusing on clonal chromosomal evolution (CCE) during relapse post-transplant. Here, we describe a 3-year-old boy with cMDS-IB, who evolved to MDS/AML. At diagnosis, the patient showed a complex karyotype with chromosomal abnormalities not yet reported together in the same clone, and likely a pathogenic variant in the ETV6 gene. The patient was treated with HLA-matched sibling HSCT, but unfortunately relapsed after HSCT with clonal cytogenetic evolution. This unfavorable outcome suggests that the treatment of pediatric MDS associated with predictive factors of poor prognosis, such as complex karyotype and ETV6 variant remains a challenge.

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process.

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

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

BackgroundDifferent 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).ResultsFrom 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.ConclusionThe 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.

Two novel PCR-based assays for sexing of Silene latifolia and Silene dioica plants

Silene latifolia and S. dioica are model systems in studies of plant reproduction, chromosome evolution and sexual dimorphism, but sexing of plants based on morphology is only possible from flowering stage onwards. Both species show homogametic females (XX) and heterogametic males (XY).•Here we developed two assays (primer pairs ss816 and ss441) for molecular sexing of S. latifolia and S. dioica, targeting length polymorphisms between the X and Y-linked copies of the spermidine synthase gene.The two assays were successful in identifying known (flowering-stage) males and females from UK and Spanish populations, with an error rate of 3.1% (ss816; successful for both species) and 0% (ss441, only successful for S. latifolia). Our assays therefore represent novel tools for rapid, robust and simple determination of the genotypic sex of S. latifolia and S. dioica.

Repetitive DNAs and chromosome evolution in Megaleporinus obtusidens and M. reinhardti (Characiformes: Anostomidae).

The high dynamism of repetitive DNAs is a major driver of chromosome evolution. In particular, the accumulation of repetitive DNA sequences has been reported as part of the differentiation of sex-specific chromosomes. In turn, the fish species of the genus Megaleporinus are a monophyletic clade in which the presence of differentiated ZZ/ZW sex chromosomes represents a synapomorphic condition, thus serving as a suitable model to evaluate the dynamic evolution of repetitive DNA classes. Therefore, transposable elements (TEs) and in tandem repeats were isolated and located on chromosomes of Megaleporinus obtusidens and M. reinhardti to infer their role in chromosome differentiation with emphasis on sex chromosome systems. Despite the conserved karyotype features of both species, the location of repetitive sequences - Rex 1, Rex 3, (TTAGGG)n, (GATA)n, (GA)n, (CA)n, and (A)n - varied both intra and interspecifically, being mainly accumulated in Z and W chromosomes. The physical mapping of repetitive sequences confirmed the remarkable dynamics of repetitive DNA classes on sex chromosomes that might have promoted chromosome diversification and reproductive isolation in Megaleporinus species.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BackgroundWhole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.ResultsKaryotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.ConclusionChromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

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 (<10MYA) and neo (<200,000 years) sex chromosomes in the XYY cytotype of the dioecious plant Rumex hastatulus and a hermaphroditic outgroup R. 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.

Remarkable chromosomes and karyotypes: A top 10 list.

Chromosomes and karyotypes are particularly rich in oddities and extremes. Described below are 10 remarkable chromosomes and karyotypes sprinkled throughout the tree of life. These include variants in chromosome number, structure, and dynamics both natural and engineered. This versatility highlights the robustness and tolerance of the mitotic and meiotic machinery to dramatic changes in chromosome and karyotype architecture. These examples also illustrate that the robustness comes at a cost, enabling the evolution of chromosomes that subvert mitosis and meiosis.

Gene transfer agents: The ambiguous role of selfless viruses in genetic exchange and bacterial evolution.

Gene transfer agents (GTAs) are genetic elements derived from ancestral bacteriophages that have become domesticated by the host. GTAs are present in diverse prokaryotic organisms, where they can facilitate horizontal gene transfer under certain conditions. Unlike typical bacteriophages, GTAs do not exhibit any preference for the replication or transfer of the genes encoding them; instead, they exhibit a remarkable capacity to package chromosomal, and sometimes extrachromosomal, DNA into virus-like capsids and disseminate it to neighboring cells. Because GTAs resemble defective prophages, identification of novel GTAs is not trivial. The detection of candidates relies on the genetic similarity to known GTAs, which has been fruitful in α-proteobacterial lineages but challenging in more distant bacteria. Here we consider several fundamental questions: What is the true prevalence of GTAs in prokaryote genomes? Given there are high costs for GTA production, what advantage do GTAs provide to the bacterial host to justify their maintenance? How is the bacterial chromosome recognized and processed for inclusion in GTA particles? This article highlights the challenges in comprehensively understanding GTAs' prevalence, function and DNA packaging method. Going forward, broad study of atypical GTAs and use of ecologically relevant conditions are required to uncover their true impact on bacterial chromosome evolution.

Absence of chromosome axis protein recruitment prevents meiotic recombination chromosome-wide in the budding yeast Lachancea kluyveri

Meiotic recombination shows broad variations across species and along chromosomes and is often suppressed at and around genomic regions determining sexual compatibility such as mating type loci in fungi. Here, we show that the absence of Spo11-DSBs and meiotic recombination on Lakl0C-left, the chromosome arm containing the sex locus of the Lachancea kluyveri budding yeast, results from the absence of recruitment of the two chromosome axis proteins Red1 and Hop1, essential for proper Spo11-DSBs formation. Furthermore, cytological observation of spread pachytene meiotic chromosomes reveals that Lakl0C-left does not undergo synapsis. However, we show that the behavior of Lakl0C-left is independent of its particularly early replication timing and is not accompanied by any peculiar chromosome structure as detectable by Hi-C in this yet poorly studied yeast. Finally, we observed an accumulation of heterozygous mutations on Lakl0C-left and a sexual dimorphism of the haploid meiotic offspring, supporting a direct effect of this absence of meiotic recombination on L. kluyveri genome evolution and fitness. Because suppression of meiotic recombination on sex chromosomes is widely observed across eukaryotes, the mechanism for recombination suppression described here may apply to other species, with the potential to impact sex chromosome evolution.

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots.

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Genomic insights into the chromosomal elongation in a family of Collembola.

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

The genomic study of repetitive elements in Solea senegalensis reveals multiple impacts of transposable elements in the evolution and architecture of Pleuronectiformes chromosomes

Pleuronectiformes are flatfishes with high commercial value and a prominent example of successful marine adaptation through chromosomal evolution. Hence, the aim of this study was to analyze the 14 relative abundance of repetitive elements (satellite DNA and transposable elements (TE)) in the 15 genome of 10 fish species (8 flatfish) delving into the study of the species of special relevance, 16 Senegalese sole, Solea senegalensis. The results showed differences in the abundance of repetitive elements, with S. senegalensis exhibiting the highest frequency and coverage of these elements reaching the 40% of the genome and not at random distribution. It is noteworthy the presence of relevant peaks of Helitrons in centromeric/pericentromeric positions mainly in the bi-armed chromosomes 1, 2, 4, 6, 7, and 9. The position of the centromeres of this species determined through the genomic localization of the family of satellite DNA PvuII, and other repetitive sequences was obtained de novo. This allowed us to know the genomic position of the centromeres in 19 out of the 21 chromosomes of S. senegalensis. Helitrons showed an accumulation of tandem copies mainly in the pericentromeric positions of chromosomes 1 and 2, occupying a region, in the first case, of 600Kb of tandem repeats. That has only been previously described in mammals and plants. Divergence and copy number studies indicated the presence of active families in the species’ genome and the existence of two important events of transposon activity (burst) in the genome of S. senegalensis, mainly accentuated in Helitrons. The results showed that only the families of DNA transposons exhibited a landscape with symmetrical bell-shaped distribution. The phylogenetic analysis of Helitron families revealed the presence of two large groups of families and the presence of four groups of sequences with heterogeneous distribution among chromosomes. Finally, the phylogenomic analysis of 8615 sequences belonging to Helitron insertions from 5 families of flatfish and two external species, allowed to classify the copies into nine groups of sequences with different levels of divergence and clusters, including some branches with distant phylogenetically species. The implications of this study will help to expand the knowledge of chromosome structure and evolution of these species.

Anatomy of a mega-radiation: Biogeography and niche evolution in Astragalus.

Astragalus (Fabaceae), with more than 3000 species, represents a globally successful radiation of morphologically highly similar species predominant across the northern hemisphere. It has attracted attention from systematists and biogeographers, who have asked what factors might be behind the extraordinary diversity of this important arid-adapted clade and what sets it apart from close relatives with far less species richness. Here, for the first time using extensive phylogenetic sampling, we asked whether (1) Astragalus is uniquely characterized by bursts of radiation or whether diversification instead is uniform and no different from closely related taxa. Then we tested whether the species diversity of Astragalus is attributable specifically to its predilection for (2) cold and arid habitats, (3) particular soils, or to (4) chromosome evolution. Finally, we tested (5) whether Astragalus originated in central Asia as proposed and (6) whether niche evolutionary shifts were subsequently associated with the colonization of other continents. Our results point to the importance of heterogeneity in the diversification of Astragalus, with upshifts associated with the earliest divergences but not strongly tied to any abiotic factor or biogeographic regionalization tested here. The only potential correlate with diversification we identified was chromosome number. Biogeographic shifts have a strong association with the abiotic environment and highlight the importance of central Asia as a biogeographic gateway. Our investigation shows the importance of phylogenetic and evolutionary studies of logistically challenging "mega-radiations." Our findings reject any simple key innovation behind high diversity and underline the often nuanced, multifactorial processes leading to species-rich clades.

Heteromorphic ZZ/ZW sex chromosomes sharing gene content with mammalian XX/XY are conserved in Madagascan chameleons of the genus Furcifer

Chameleons are well-known lizards with unique morphology and physiology, but their sex determination has remained poorly studied. Madagascan chameleons of the genus Furcifer have cytogenetically distinct Z and W sex chromosomes and occasionally Z1Z1Z2Z2/Z1Z2W multiple neo-sex chromosomes. To identify the gene content of their sex chromosomes, we microdissected and sequenced the sex chromosomes of F. oustaleti (ZZ/ZW) and F. pardalis (Z1Z1Z2Z2/Z1Z2W). In addition, we sequenced the genomes of a male and a female of F. lateralis (ZZ/ZW) and F. pardalis and performed a comparative coverage analysis between the sexes. Despite the notable heteromorphy and distinctiveness in heterochromatin content, the Z and W sex chromosomes share approximately 90% of their gene content. This finding demonstrates poor correlation of the degree of differentiation of sex chromosomes at the cytogenetic and gene level. The test of homology based on the comparison of gene copy number variation revealed that female heterogamety with differentiated sex chromosomes remained stable in the genus Furcifer for at least 20 million years. These chameleons co-opted for the role of sex chromosomes the same genomic region as viviparous mammals, lacertids and geckos of the genus Paroedura, which makes these groups excellent model for studies of convergent and divergent evolution of sex chromosomes.

Turnovers of Sex-Determining Mutation in the Golden Pompano and Related Species Provide Insights into Microevolution of Undifferentiated Sex Chromosome.

The suppression of recombination is considered a hallmark of sex chromosome evolution. However, previous research has identified undifferentiated sex chromosomes and sex determination by single SNP in the greater amberjack (Seriola dumerili). We observed the same phenomena in the golden pompano (Trachinotus ovatus) of the same family Carangidae and discovered a different sex-determining SNP within the same gene Hsd17b1. We propose an evolutionary model elucidating the turnover of sex-determining mutations by highlighting the contrasting dynamics between purifying selection, responsible for maintaining W-linked Hsd17b1, and neutral evolution, which drives Z-linked Hsd17b1. Additionally, sporadic loss-of-function mutations in W-linked Hsd17b1 contribute to the conversion of W chromosomes into Z chromosomes. This model was directly supported by simulations, closely related species, and indirectly by zebrafish mutants. These findings shed new light on the early stages of sex chromosome evolution.

Phylogeography and ecological niche modelling of the South American cricetid rodent Graomys griseoflavus, with insights into its chromosomal evolution

Abstract The rodent Graomys griseoflavus has a wide geographical distribution in the Arid Diagonal of South America, showing variation in the diploid number (2n = 33–38) caused by three different Robertsonian (Rb) translocations. Two contrasting hypotheses have been proposed to explain the origin and evolution of this polymorphism: (i) unique and sequential events; or (ii) multiple and independent origins of Rb translocations. Here, we combined phylogeography and ecological niche modelling to elucidate the evolutionary history of G. griseoflavus and help to contrast these hypotheses. The results indicated a demographic increase that would have begun ~150 kya from the High Monte ecoregion. The spatial distributions of these populations were coincident with the areas of highest habitat suitability. The palaeoclimatic projections suggested the presence of two main glacial refugia coincident with the main genetic clusters. Graomys griseoflavus underwent two pulses of southward expansion at ~220 and ~210 kya. More recently (~20 kya), there was a secondary contact between the southern populations (Low Monte ecoregion) expanding northwards and the northern populations (High Monte ecoregion) expanding into marginal areas of the Chaco ecoregion. Combined with the cytogenetic evidence, our results support the multiple and independent origins of Rb translocations.

Chromosomal evolution, environmental heterogeneity, and migration drive spatial patterns of species richness in Calochortus (Liliaceae)

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness in Calochortus (Liliaceae, 74 spp.). Calochortus occupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades-inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)-began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence.

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.