Avian W and mammalian Y chromosomes convergently retained dosage-sensitive regulators

The Human Developing Cerebral Cortex Is Characterized by an Elevated De Novo Expression of Long Noncoding RNAs in Excitatory Neurons

The outstanding human cognitive capacities are computed in the cerebral cortex, a mammalian-specific brain region and the place of massive biological innovation. Long noncoding RNAs have emerged as gene regulatory elements with higher evolutionary turnover than mRNAs. The many long noncoding RNAs identified in neural tissues make them candidates for molecular sources of cerebral cortex evolution and disease. Here, we characterized the genomic and cellular shifts that occurred during the evolution of the long noncoding RNA repertoire expressed in the developing cerebral cortex and explored putative roles for these long noncoding RNAs in the evolution of the human brain. Using transcriptomics and comparative genomics, we comprehensively annotated the cortical transcriptomes of humans, rhesus macaques, mice, and chickens and classified human cortical long noncoding RNAs into evolutionary groups as a function of their predicted minimal ages. Long noncoding RNA evolutionary groups showed differences in expression levels, splicing efficiencies, transposable element contents, genomic distributions, and transcription factor binding to their promoters. Furthermore, older long noncoding RNAs showed preferential expression in germinative zones, outer radial glial cells, and cortical inhibitory (GABAergic) neurons. In comparison, younger long noncoding RNAs showed preferential expression in cortical excitatory (glutamatergic) neurons, were enriched in primate and human-specific gene co-expression modules, and were dysregulated in neurodevelopmental disorders. These results suggest different evolutionary routes for older and younger cortical long noncoding RNAs, highlighting old long noncoding RNAs as a possible source of molecular evolution of conserved developmental programs; conversely, we propose that the de novo expression of primate- and human-specific young long noncoding RNAs is a putative source of molecular evolution and dysfunction of cortical excitatory neurons, warranting further investigation.

Inverse and Proportional Trans Modulation of Gene Expression in Human Aneuploidies

Genomic imbalance in aneuploidy is often detrimental to organisms. To gain insight into the molecular basis of aneuploidies in humans, we analyzed transcriptome data from several autosomal and sex chromosome aneuploidies. The results showed that in human aneuploid cells, genes located on unvaried chromosomes are inversely or proportionally trans-modulated, while a subset of genes on the varied chromosomes are compensated. Less genome-wide modulation is found for sex chromosome aneuploidy compared with autosomal aneuploidy due to X inactivation and the retention of dosage sensitive regulators on both sex chromosomes to limit the effective dosage change. We also found that lncRNA and mRNA can have different responses to aneuploidy. Furthermore, we analyzed the relationship between dosage-sensitive transcription factors and their targets, which illustrated the modulations and indicates genomic imbalance is related to stoichiometric changes in components of gene regulatory complexes.In summary, this study demonstrates the existence of trans-acting effects and compensation mechanisms in human aneuploidies and contributes to our understanding of gene expression regulation in unbalanced genomes and disease states.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution

BACKGROUND

Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected.

METHODS

We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes.

RESULTS

Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host.

CONCLUSIONS

Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Satellite DNAs, heterochromatin, and sex chromosomes of the wattled jacana (Charadriiformes; Jacanidae): a species with highly rearranged karyotype

Charadriiformes, which comprises shorebirds and their relatives, is one of the most diverse avian orders, with over 390 species showing a wide range of karyotypes. Here, we isolated and characterized the whole collection of satellite DNAs (satDNAs) at both molecular and cytogenetic levels of one of its representative species, named the wattled jacana (Jacana jacana), a species that contains a typical ZZ/ZW sex chromosome system and a highly rearranged karyotype. In addition, we also investigate the in situ location of telomeric and microsatellite repeats. A small catalog of 11 satDNAs was identified that typically accumulated on microchromosomes and on the W chromosome. The latter also showed a significant accumulation of telomeric signals, being (GA)10 the only microsatellite with positive hybridization signals among all the 16 tested ones. These current findings contribute to our understanding of the genomic organization of repetitive DNAs in a bird species with high degree of chromosomal reorganization contrary to the majority of bird species that have stable karyotypes.

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.

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

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

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba)

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

An effective strategy for assembling the sex-limited chromosome

BACKGROUND

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

FINDINGS

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

CONCLUSIONS

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

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

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

Gene Expression Profiling Reveals Potential Players of Sex Determination and Asymmetrical Development in Chicken Embryo Gonads

Despite the notable progress made in recent years, the understanding of the genetic control of gonadal sex differentiation and asymmetrical ovariogenesis in chicken during embryonic development remains incomplete. This study aimed to identify potential key genes and speculate about the mechanisms associated with ovary and testis development via an analysis of the results of PacBio and Illumina transcriptome sequencing of embryonic chicken gonads at the initiation of sexual differentiation (E4.5, E5.5, and E6.5). PacBio sequencing detected 328 and 233 significantly up-regulated transcript isoforms in females and males at E4.5, respectively. Illumina sequencing detected 95, 296 and 445 DEGs at E4.5, E5.5, and E6.5, respectively. Moreover, both sexes showed asymmetrical expression in gonads, and more DEGs were detected on the left side. There were 12 DEGs involved in cell proliferation shared between males and females in the left gonads. GO analysis suggested that coagulation pathways may be involved in the degradation of the right gonad in females and that blood oxygen transport pathways may be involved in preventing the degradation of the right gonad in males. These results provide a comprehensive expression profile of chicken embryo gonads at the initiation of sexual differentiation, which can serve as a theoretical basis for further understanding the mechanism of bird sex determination and its evolutionary process.

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5 Collapse

Key Words Collapse

MESH Headings

• Animals
• Female
• Male
• Chromatin/genetics
• Crows/genetics
• Epigenesis, Genetic
• Methylation
• Dosage Compensation, Genetic
• Sex Chromosomes

Collapse

Grants

• HORIZON EUROPE European Research Council
• Swedish Research Council
• Knut och Alice Wallenbergs Stiftelse

Collapse

Affiliation(s)

Ana Catalán Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Justin Merondun Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Ulrich Knief Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany Evolutionary Biology & Ecology,Faculty of Biology, University of Freiburg, Freiburg, Germany
Jochen B. W. Wolf Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany

Collapse

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

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

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

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.

Preliminary Study on Expression and Function of the Chicken W Chromosome Gene MIER3 in Embryonic Gonads

The sex chromosomes of birds are designated Z and W. The male is homogamous (ZZ), and the female is heterogamous (ZW). The chicken W chromosome is a degenerate version of the Z chromosome and harbors only 28 protein-coding genes. We studied the expression pattern of the W chromosome gene MIER3 (showing differential expression during gonadogenesis) in chicken embryonic gonads and its potential role in gonadal development. The W copy of MIER3 (MIER3-W) shows a gonad-biased expression in chicken embryonic tissues which was different from its Z copy. The overall expression of MIER3-W and MIER3-Z mRNA and protein is correlated with the gonadal phenotype being higher in female gonads than in male gonads or female-to-male sex-reversed gonads. Chicken MIER3 protein is highly expressed in the nucleus, with relatively lower expression in the cytoplasm. Overexpression of MIER3-W in male gonad cells suggested its effect on the GnRH signaling pathway, cell proliferation, and cell apoptosis. MIER3 expression is associated with the gonadal phenotype. MIER3 may promote female gonadal development by regulating EGR1 and αGSU genes. These findings enrich our knowledge of chicken W chromosome genes and support a more systematic and in-depth understanding of gonadal development in chickens.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

A new haplotype-resolved turkey genome to enable turkey genetics and genomics research

BACKGROUND

The domesticated turkey (Meleagris gallopavo) is a species of significant agricultural importance and is the second largest contributor, behind broiler chickens, to world poultry meat production. The previous genome is of draft quality and partly based on the chicken (Gallus gallus) genome. A high-quality reference genome of M. gallopavo is essential for turkey genomics and genetics research and the breeding industry.

RESULTS

By adopting the trio-binning approach, we were able to assemble a high-quality chromosome-level F1 assembly and 2 parental haplotype assemblies, leveraging long-read technologies and genome-wide chromatin interaction data (Hi-C). From a total of 40 chromosomes (2n = 80), we captured 35 chromosomes in a single scaffold, showing much improved genome completeness and continuity compared to the old assembly build. The 3 assemblies are of higher quality than the previous draft quality assembly and comparable to the chicken assemblies (GRCg7) shown by the largest contig N50 (26.6 Mb) and comparable BUSCO gene set completeness scores (96-97%). Comparative analyses confirm a previously identified large inversion of around 19 Mbp on the Z chromosome not found in other Galliformes. Structural variation between the parent haplotypes was identified, which poses potential new target genes for breeding.

CONCLUSIONS

We contribute a new high-quality turkey genome at the chromosome level, benefiting turkey genetics and other avian genomics research as well as the turkey breeding industry.

The contribution of sex chromosome conflict to disrupted spermatogenesis in hybrid house mice

Incompatibilities on the sex chromosomes are important in the evolution of hybrid male sterility, but the evolutionary forces underlying this phenomenon are unclear. House mice (Mus musculus) lineages have provided powerful models for understanding the genetic basis of hybrid male sterility. X chromosome-autosome interactions cause strong incompatibilities in M. musculus F1 hybrids, but variation in sterility phenotypes suggests a more complex genetic basis. In addition, XY chromosome conflict has resulted in rapid expansions of ampliconic genes with dosage-dependent expression that is essential to spermatogenesis. Here, we evaluated the contribution of XY lineage mismatch to male fertility and stage-specific gene expression in hybrid mice. We performed backcrosses between two house mouse subspecies to generate reciprocal Y-introgression strains and used these strains to test the effects of XY mismatch in hybrids. Our transcriptome analyses of sorted spermatid cells revealed widespread overexpression of the X chromosome in sterile F1 hybrids independent of Y chromosome subspecies origin. Thus, postmeiotic overexpression of the X chromosome in sterile F1 mouse hybrids is likely a downstream consequence of disrupted meiotic X-inactivation rather than XY gene copy number imbalance. Y chromosome introgression did result in subfertility phenotypes and disrupted expression of several autosomal genes in mice with an otherwise nonhybrid genomic background, suggesting that Y-linked incompatibilities contribute to reproductive barriers, but likely not as a direct consequence of XY conflict. Collectively, these findings suggest that rapid sex chromosome gene family evolution driven by genomic conflict has not resulted in strong male reproductive barriers between these subspecies of house mice.

Blood transcriptome comparison between sexes and their function in 4-week Rhode Island red chickens

Sex is a major biological factor in the development and physiology of a sexual reproductive organism, and its role in the growing process is needed to be investigated in various species. We compare blood transcriptome between 5 males and 5 females in 4-week-old Rhode Island Red chickens and perform functional annotation of differentially expressed genes (DEGs). The results are as follows. 141 and 109 DEGs were located in autosomes and sex chromosomes, respectively. The gene ontology (GO) terms are significantly (p < 0.05) enriched, which were limb development, inner ear development, positive regulation of dendrite development, the KEGG pathway the TGF-beta signaling pathway, and melanogenesis (p < 0.05). These pathways are related to morphological maintenance and growth of the tissues. In addition, the SMAD2W and the BMP5 were involved in the TGF-beta signaling pathway, and both play an important role in maintaining tissue development. The major DEGs related to the development of neurons and synapses include the up-regulated NRN1, GDF10, SLC1A1, BMP5, NBEA, and NRXN1. Also, 7 DEGs were validated using RT-qPCR with high correlation (r ² = 0.74). In conclusion, the differential expression of blood tissue in the early growing chicken was enriched in TGF-beta signaling and related to the development of neurons and synapses including SMAD2W and BMP5. These results suggest that blood in the early growing stage is differentially affected in tissue development, nervous system, and pigmentation by sex. For future research, experimental characterization of DEGs and a holistic investigation of various tissues and growth stages will be required.

The Gallus gallus RJF reference genome reveals an MHCY haplotype organized in gene blocks that contain 107 loci including 45 specialized, polymorphic MHC class I loci, 41 C-type lectin-like loci, and other loci amid hundreds of transposable elements

MHCY is a second major histocompatibility complex-like gene region in chickens originally identified by the presence of major histocompatibility complex class I-like and class II-like gene sequences. Up to now, the MHCY gene region has been poorly represented in genomic sequence data. A high density of repetitive sequence and multiple members of several gene families prevented the accurate assembly of short-read sequence data for MHCY. Identified here by single-molecule real-time sequencing sequencing of BAC clones for the Gallus gallus Red Jungle Fowl reference genome are 107 MHCY region genes (45 major histocompatibility complex class I-like, 41 c-type-lectin-like, 8 major histocompatibility complex class IIβ, 8 LENG9-like, 4 zinc finger protein loci, and a single only zinc finger-like locus) located amid hundreds of retroelements within 4 contigs representing the region. Sequences obtained for nearby ribosomal RNA genes have allowed MHCY to be precisely mapped with respect to the nucleolar organizer region. Gene sequences provide insights into the unusual structure of the MHCY class I molecules. The MHCY class I loci are polymorphic and group into 22 types based on predicted amino acid sequences. Some MHCY class I loci are full-length major histocompatibility complex class I genes. Others with altered gene structure are considered gene candidates. The amino acid side chains at many of the polymorphic positions in MHCY class I are directed away rather than into the antigen-binding groove as is typical of peptide-binding major histocompatibility complex class I molecules. Identical and nearly identical blocks of genomic sequence contribute to the observed multiplicity of identical MHCY genes and the large size (&gt;639 kb) of the Red Jungle Fowl MHCY haplotype. Multiple points of hybridization observed in fluorescence in situ hybridization suggest that the Red Jungle Fowl MHCY haplotype is made up of linked, but physically separated genomic segments. The unusual gene content, the evidence of highly similar duplicated segments, and additional evidence of variation in haplotype size distinguish polymorphic MHCY from classical polymorphic major histocompatibility complex regions.

Prediction of sex-determination mechanisms in avian primordial germ cells using RNA-seq analysis

In birds, sex is determined through cell-autonomous mechanisms and various factors, such as the dosage of DMRT1. While the sex-determination mechanism in gonads is well known, the mechanism in germ cells remains unclear. In this study, we explored the gene expression profiles of male and female primordial germ cells (PGCs) during embryogenesis in chickens to predict the mechanism underlying sex determination. Male and female PGCs were isolated from blood and gonads with a purity > 96% using flow cytometry and analyzed using RNA-seq. Prior to settlement in the gonads, female circulating PGCs (cPGCs) obtained from blood displayed sex-biased expression. Gonadal PGCs (gPGCs) also exhibited sex-biased expression, and the number of female-biased genes detected was higher than that of male-biased genes. The female-biased genes in gPGCs were enriched in some metabolic processes. To reveal the mechanisms underlying the transcriptional regulation of female-biased genes in gPGCs, we performed stimulation tests. Retinoic acid stimulation of cultured gPGCs derived from male embryos resulted in the upregulation of several female-biased genes. Overall, our results suggest that sex determination in avian PGCs involves aspects of both cell-autonomous and somatic-cell regulation. Moreover, it appears that sex determination occurs earlier in females than in males.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

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

Sex-limited chromosomes and non-reproductive traits

Sex chromosomes are typically viewed as having originated from a pair of autosomes, and differentiated as the sex-limited chromosome (e.g. Y) has degenerated by losing most genes through cessation of recombination. While often thought that degenerated sex-limited chromosomes primarily affect traits involved in sex determination and sex cell production, accumulating evidence suggests they also influence traits not sex-limited or directly involved in reproduction. Here, we provide an overview of the effects of sex-limited chromosomes on non-reproductive traits in XY, ZW or UV sex determination systems, and discuss evolutionary processes maintaining variation at sex-limited chromosomes and molecular mechanisms affecting non-reproductive traits.

Fragile, unfaithful and persistent Ys-on how meiosis can shape sex chromosome evolution

Sex-linked inheritance is a stark exception to Mendel's Laws of Heredity. Here we discuss how the evolution of heteromorphic sex chromosomes (mainly the Y) has been shaped by the intricacies of the meiotic programme. We propose that persistence of Y chromosomes in distantly related mammalian phylogroups can be explained in the context of pseudoautosomal region (PAR) size, meiotic pairing strategies, and the presence of Y-borne executioner genes that regulate meiotic sex chromosome inactivation. We hypothesise that variation in PAR size can be an important driver for the evolution of recombination frequencies genome wide, imposing constraints on Y fate. If small PAR size compromises XY segregation during male meiosis, the stress of producing aneuploid gametes could drive function away from the Y (i.e., a fragile Y). The Y chromosome can avoid fragility either by acquiring an achiasmatic meiotic XY pairing strategy to reduce aneuploid gamete production, or gain meiotic executioner protection (a persistent Y). Persistent Ys will then be under strong pressure to maintain high recombination rates in the PAR (and subsequently genome wide), as improper segregation has fatal consequences for germ cells. In the event that executioner protection is lost, the Y chromosome can be maintained in the population by either PAR rejuvenation (extension by addition of autosome material) or gaining achiasmatic meiotic pairing, the alternative is Y loss. Under this dynamic cyclic evolutionary scenario, understanding the meiotic programme in vertebrate and invertebrate species will be crucial to further understand the plasticity of the rise and fall of heteromorphic sex chromosomes.

A high-quality assembly reveals genomic characteristics, phylogenetic status, and causal genes for leucism plumage of Indian peafowl

Background

The dazzling phenotypic characteristics of male Indian peafowl (Pavo cristatus) are attractive both to the female of the species and to humans. However, little is known about the evolution of the phenotype and phylogeny of these birds at the whole-genome level. So far, there are no reports regarding the genetic mechanism of the formation of leucism plumage in this variant of Indian peafowl.

Results

A draft genome of Indian peafowl was assembled, with a genome size of 1.05 Gb (the sequencing depth is 362×), and contig and scaffold N50 were up to 6.2 and 11.4 Mb, respectively. Compared with other birds, Indian peafowl showed changes in terms of metabolism, immunity, and skeletal and feather development, which provided a novel insight into the phenotypic evolution of peafowl, such as the large body size and feather morphologies. Moreover, we determined that the phylogeny of Indian peafowl was more closely linked to turkey than chicken. Specifically, we first identified that PMEL was a potential causal gene leading to the formation of the leucism plumage variant in Indian peafowl.

Conclusions

This study provides an Indian peafowl genome of high quality, as well as a novel understanding of phenotypic evolution and phylogeny of Indian peafowl. These results provide a valuable reference for the study of avian genome evolution. Furthermore, the discovery of the genetic mechanism for the development of leucism plumage is both a breakthrough in the exploration of peafowl plumage and also offers clues and directions for further investigations of the avian plumage coloration and artificial breeding in peafowl.

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

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

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

Sexually dimorphic expression of a chicken sex chromosome gene (VCP) reflects differences in gonadal development between males and females

The chicken has a Z-W sex chromosome system, in which the males are the homogametic sex (ZZ) and the females the heterogametic sex (ZW). The smaller W chromosome is generally considered to be a highly degraded copy of the Z chromosome that retains around 28-30 homologous protein-coding genes' These Z-W homologues are thought to have important, but undefined, roles in development, and here we explore the role of one of these genes, VCP (Valosin Containing Protein) in gonadogenesis. We established RNA expression levels of both Z and W VCP homologues, the levels of VCP protein, and the cellular localization of VCP protein in male and female embryonic gonads during development. We also assessed the effects of female-to-male sex-reversal on VCP expression in developing gonads. The results showed that both VCP RNA and protein are expressed at higher levels in female than male gonads, and the expression levels of VCP protein and VCP-Z transcript, but not VCP-W transcript, are decreased in female-to-male sex reversed gonads. In addition, the spatial expression of VCP protein differs between male and female embryonic gonads: in testes, VCP protein is mainly confined to the medullary sex cords, while in ovaries, VCP protein is expressed throughout the medulla and at higher levels in the cortex. The results suggest that sexually dimorphic expression of chicken VCP reflects differences in gonadal morphology between sexes.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Sex-related differences in aging rate are associated with sex chromosome system in amphibians

Sex‐related differences in mortality are widespread in the animal kingdom. Although studies have shown that sex determination systems might drive lifespan evolution, sex chromosome influence on aging rates have not been investigated so far, likely due to an apparent lack of demographic data from clades including both XY (with heterogametic males) and ZW (heterogametic females) systems. Taking advantage of a unique collection of capture–recapture datasets in amphibians, a vertebrate group where XY and ZW systems have repeatedly evolved over the past 200 million years, we examined whether sex heterogamy can predict sex differences in aging rates and lifespans. We showed that the strength and direction of sex differences in aging rates (and not lifespan) differ between XY and ZW systems. Sex‐specific variation in aging rates was moderate within each system, but aging rates tended to be consistently higher in the heterogametic sex. This led to small but detectable effects of sex chromosome system on sex differences in aging rates in our models. Although preliminary, our results suggest that exposed recessive deleterious mutations on the X/Z chromosome (the “unguarded X/Z effect”) or repeat‐rich Y/W chromosome (the “toxic Y/W effect”) could accelerate aging in the heterogametic sex in some vertebrate clades.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

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

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

Genetic Regulation of Avian Testis Development

As in other vertebrates, avian testes are the site of spermatogenesis and androgen production. The paired testes of birds differentiate during embryogenesis, first marked by the development of pre-Sertoli cells in the gonadal primordium and their condensation into seminiferous cords. Germ cells become enclosed in these cords and enter mitotic arrest, while steroidogenic Leydig cells subsequently differentiate around the cords. This review describes our current understanding of avian testis development at the cell biology and genetic levels. Most of this knowledge has come from studies on the chicken embryo, though other species are increasingly being examined. In chicken, testis development is governed by the Z-chromosome-linked DMRT1 gene, which directly or indirectly activates the male factors, HEMGN, SOX9 and AMH. Recent single cell RNA-seq has defined cell lineage specification during chicken testis development, while comparative studies point to deep conservation of avian testis formation. Lastly, we identify areas of future research on the genetics of avian testis development.

The avian W chromosome is a refugium for endogenous retroviruses with likely effects on female-biased mutational load and genetic incompatibilities

It is a broadly observed pattern that the non-recombining regions of sex-limited chromosomes (Y and W) accumulate more repeats than the rest of the genome, even in species like birds with a low genome-wide repeat content. Here, we show that in birds with highly heteromorphic sex chromosomes, the W chromosome has a transposable element (TE) density of greater than 55% compared to the genome-wide density of less than 10%, and contains over half of all full-length (thus potentially active) endogenous retroviruses (ERVs) of the entire genome. Using RNA-seq and protein mass spectrometry data, we were able to detect signatures of female-specific ERV expression. We hypothesize that the avian W chromosome acts as a refugium for active ERVs, probably leading to female-biased mutational load that may influence female physiology similar to the 'toxic-Y' effect in Drosophila males. Furthermore, Haldane's rule predicts that the heterogametic sex has reduced fertility in hybrids. We propose that the excess of W-linked active ERVs over the rest of the genome may be an additional explanatory variable for Haldane's rule, with consequences for genetic incompatibilities between species through TE/repressor mismatches in hybrids. Together, our results suggest that the sequence content of female-specific W chromosomes can have effects far beyond sex determination and gene dosage. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

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

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

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

The timing of genetic degeneration of sex chromosomes

Genetic degeneration is an extraordinary feature of sex chromosomes, with the loss of functions of Y-linked genes in species with XY systems, and W-linked genes in ZW systems, eventually affecting almost all genes. Although degeneration is familiar to most biologists, important aspects are not yet well understood, including how quickly a Y or W chromosome can become completely degenerated. I review the current understanding of the time-course of degeneration. Degeneration starts after crossing over between the sex chromosome pair stops, and theoretical models predict an initially fast degeneration rate and a later much slower one. It has become possible to estimate the two quantities that the models suggest are the most important in determining degeneration rates-the size of the sex-linked region, and the time when recombination became suppressed (which can be estimated using Y-X or W-Z sequence divergence). However, quantifying degeneration is still difficult. I review evidence on gene losses (based on coverage analysis) or loss of function (by classifying coding sequences into functional alleles and pseudogenes). I also review evidence about whether small genome regions degenerate, or only large ones, whether selective constraints on the genes in a sex-linked region also strongly affect degeneration rates, and about how long it takes before all (or almost all) genes are lost. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Schistosome W-linked genes inform temporal dynamics of sex chromosome evolution and suggest candidate for sex determination

Schistosomes, the human parasites responsible for snail fever, are female-heterogametic. Different parts of their ZW sex chromosomes have stopped recombining in distinct lineages, creating "evolutionary strata" of various ages. While the Z-chromosome is well characterized at the genomic and molecular level, the W-chromosome has remained largely unstudied from an evolutionary perspective, as only a few W-linked genes have been detected outside of the model species Schistosoma mansoni. Here, we characterize the gene content and evolution of the W-chromosomes of S. mansoni and of the divergent species S. japonicum. We use a combined RNA/DNA k-mer based pipeline to assemble around one hundred candidate W-specific transcripts in each of the species. About half of them map to known protein coding genes, the majority homologous to S. mansoni Z-linked genes. We perform an extended analysis of the evolutionary strata present in the two species (including characterizing a previously undetected young stratum in S. japonicum) to infer patterns of sequence and expression evolution of W-linked genes at different time points after recombination was lost. W-linked genes show evidence of degeneration, including high rates of protein evolution and reduced expression. Most are found in young lineage-specific strata, with only a few high expression ancestral W-genes remaining, consistent with the progressive erosion of non-recombining regions. Among these, the splicing factor U2AF2 stands out as a promising candidate for primary sex determination, opening new avenues for understanding the molecular basis of the reproductive biology of this group.

Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae)

Heteromorphic W and Y sex chromosomes often experience gene loss and heterochromatinization, which is frequently viewed as their "degeneration". However, the evolutionary trajectories of the heterochromosomes are in fact more complex since they may not only lose but also acquire new sequences. Previously, we found that the heterochromatic W chromosome of a lizard Eremias velox (Lacertidae) is decondensed and thus transcriptionally active during the lampbrush stage. To determine possible sources of this transcription, we sequenced DNA from a microdissected W chromosome sample and a total female DNA sample and analyzed the results of reference-based and de novo assembly. We found a new repetitive sequence, consisting of fragments of an autosomal protein-coding gene ATF7IP2, several SINE elements, and sequences of unknown origin. This repetitive element is distributed across the whole length of the W chromosome, except the centromeric region. Since it retained only 3 out of 10 original ATF7IP2 exons, it remains unclear whether it is able to produce a protein product. Subsequent studies are required to test the presence of this element in other species of Lacertidae and possible functionality. Our results provide further evidence for the view of W and Y chromosomes as not just "degraded" copies of Z and X chromosomes but independent genomic segments in which novel genetic elements may arise.

Metagenomic Analysis of the Respiratory Microbiome of a Broiler Flock from Hatching to Processing

Elucidating the complex microbial interactions in biological environments requires the identification and characterization of not only the bacterial component but also the eukaryotic viruses, bacteriophage, and fungi. In a proof of concept experiment, next generation sequencing approaches, accompanied by the development of novel computational and bioinformatics tools, were utilized to examine the evolution of the microbial ecology of the avian trachea during the growth of a healthy commercial broiler flock. The flock was sampled weekly, beginning at placement and concluding at 49 days, the day before processing. Metagenomic sequencing of DNA and RNA was utilized to examine the bacteria, virus, bacteriophage, and fungal components during flock growth. The utility of using a metagenomic approach to study the avian respiratory virome was confirmed by detecting the dysbiosis in the avian respiratory virome of broiler chickens diagnosed with infection with infectious laryngotracheitis virus. This study provides the first comprehensive analysis of the ecology of the avian respiratory microbiome and demonstrates the feasibility for the use of this approach in future investigations of avian respiratory diseases.

Multi-Copy Gene Family Evolution on the Avian W Chromosome

The sex chromosomes often follow unusual evolutionary trajectories. In particular, the sex-limited chromosomes frequently exhibit a small but unusual gene content in numerous species, where many genes have undergone massive gene amplification. The reasons for this remain elusive with a number of recent studies implicating meiotic drive, sperm competition, genetic drift, and gene conversion in the expansion of gene families. However, our understanding is primarily based on Y chromosome studies as few studies have systematically tested for copy number variation on W chromosomes. Here, we conduct a comprehensive investigation into the abundance, variability, and evolution of ampliconic genes on the avian W. First, we quantified gene copy number and variability across the duck W chromosome. We find a limited number of gene families as well as conservation in W-linked gene copy number across duck breeds, indicating that gene amplification may not be such a general feature of sex chromosome evolution as Y studies would initially suggest. Next, we investigated the evolution of HINTW, a prominent ampliconic gene family hypothesized to play a role in female reproduction and oogenesis. In particular, we investigated the factors driving the expansion of HINTW using contrasts between modern chicken and duck breeds selected for different female-specific selection regimes and their wild ancestors. Although we find the potential for selection related to fecundity in explaining small-scale gene amplification of HINTW in the chicken, purifying selection seems to be the dominant mode of evolution in the duck. Together, this challenges the assumption that HINTW is key for female fecundity across the avian phylogeny.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

A new duck genome reveals conserved and convergently evolved chromosome architectures of birds and mammals

BACKGROUND

Ducks have a typical avian karyotype that consists of macro- and microchromosomes, but a pair of much less differentiated ZW sex chromosomes compared to chickens. To elucidate the evolution of chromosome architectures between ducks and chickens, and between birds and mammals, we produced a nearly complete chromosomal assembly of a female Pekin duck by combining long-read sequencing and multiplatform scaffolding techniques.

RESULTS

A major improvement of genome assembly and annotation quality resulted from the successful resolution of lineage-specific propagated repeats that fragmented the previous Illumina-based assembly. We found that the duck topologically associated domains (TAD) are demarcated by putative binding sites of the insulator protein CTCF, housekeeping genes, or transitions of active/inactive chromatin compartments, indicating conserved mechanisms of spatial chromosome folding with mammals. There are extensive overlaps of TAD boundaries between duck and chicken, and also between the TAD boundaries and chromosome inversion breakpoints. This suggests strong natural selection pressure on maintaining regulatory domain integrity, or vulnerability of TAD boundaries to DNA double-strand breaks. The duck W chromosome retains 2.5-fold more genes relative to chicken. Similar to the independently evolved human Y chromosome, the duck W evolved massive dispersed palindromic structures, and a pattern of sequence divergence with the Z chromosome that reflects stepwise suppression of homologous recombination.

CONCLUSIONS

Our results provide novel insights into the conserved and convergently evolved chromosome features of birds and mammals, and also importantly add to the genomic resources for poultry studies.

An Investigation of ZZ/ZW and XX/XY Sex Determination Systems in North African Catfish (Clarias gariepinus, )

An investigation of sex-specific loci may provide important insights into fish sex determination strategies. This may be useful for biotechnological purposes, for example, to produce all-male or all-female fish for commercial breeding. The North African catfish species, Clarias gariepinus, has been widely adopted for aquaculture because its superior growth and disease resistance render the species suitable for hybridization with other catfish to improve the productivity and quality of fish meat. This species has either a ZZ/ZW or XX/XY sex determination system. Here, we investigate and characterize these systems using high-throughput genome complexity reduction sequencing as Diversity Arrays Technology. This approach was effective in identifying moderately sex-linked loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 30 perfectly sexed individuals of C. gariepinus. However, SNPs based markers were not found in this study. In total, 41 loci met the criteria for being moderately male-linked (with male vs. female ratios 80:20 and 70:30), while 25 loci were found to be moderately linked to female sex. No strictly male- or female-linked loci were detected. Seven moderately male-linked loci were partially homologous to some classes of transposable elements and three moderately male-linked loci were partially homologous to functional genes. Our data showed that the male heterogametic XX/XY sex determination system should co-exist with the ZZ/ZW system in C. gariepinus. Our finding of the co-existence of XX/XY and ZZ/ZW systems can be applied to benefit commercial breeding of this species in Thailand. This approach using moderately sex-linked loci provides a solid baseline for revealing sex determination mechanisms and identify potential sex determination regions in catfish, allowing further investigation of genetic improvements in breeding programs.