Findzx: an automated pipeline for detecting and visualising sex chromosomes using whole-genome sequencing data.

BACKGROUND: Sex chromosomes have evolved numerous times, as revealed by recent genomic studies. However, large gaps in our knowledge of sex chromosome diversity across the tree of life remain. Filling these gaps, through the study of novel species, is crucial for improved understanding of why and how sex chromosomes evolve. Characterization of sex chromosomes in already well-studied organisms is also important to avoid misinterpretations of population genomic patterns caused by undetected sex chromosome variation. RESULTS: Here we present findZX, an automated Snakemake-based computational pipeline for detecting and visualizing sex chromosomes through differences in genome coverage and heterozygosity between any number of males and females. A main feature of the pipeline is the option to perform a genome coordinate liftover to a reference genome of another species. This allows users to inspect sex-linked regions over larger contiguous chromosome regions, while also providing important between-species synteny information. To demonstrate its effectiveness, we applied findZX to publicly available genomic data from species belonging to widely different taxonomic groups (mammals, birds, reptiles, and fish), with sex chromosome systems of different ages, sizes, and levels of differentiation. We also demonstrate that the liftover method is robust over large phylogenetic distances (> 80 million years of evolution). CONCLUSIONS: With findZX we provide an easy-to-use and highly effective tool for identification of sex chromosomes. The pipeline is compatible with both Linux and MacOS systems, and scalable to suit different computational platforms.

Genomic trajectories of a near-extinction event in the Chatham Island black robin.

BACKGROUND: Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation. RESULTS: We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations. CONCLUSION: Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

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.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome.

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree-based linkage mapping. This species has extended Z and W chromosomes ("neo-sex chromosomes") formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex-linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A-part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein-coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A-derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z-linked, recombination patterns.

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

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

A novel neo-sex chromosome in Sylvietta brachyura (Macrosphenidae) adds to the extraordinary avian sex chromosome diversity among Sylvioidea songbirds.

We report the discovery of a novel neo-sex chromosome in an African warbler, Sylvietta brachyura (northern crombec; Macrosphenidae). This species is part of the Sylvioidea superfamily, where four separate autosome-sex chromosome translocation events have previously been discovered via comparative genomics of 11 of the 22 families in this clade. Our discovery here resulted from analyses of genomic data of single species-representatives from three additional Sylvioidea families (Macrosphenidae, Pycnonotidae and Leiothrichidae). In all three species, we confirmed the translocation of a part of chromosome 4A to the sex chromosomes, which originated basally in Sylvioidea. In S. brachyura, we found that a part of chromosome 8 has been translocated to the sex chromosomes, forming a unique neo-sex chromosome in this lineage. Furthermore, the non-recombining part of 4A in S. brachyura is smaller than in other Sylvioidea species, which suggests that recombination continued along this region after the fusion event in the Sylvioidea ancestor. These findings reveal additional sex chromosome diversity among the Sylvioidea, where five separate translocation events are now confirmed.

Why Do Sex Chromosomes Stop Recombining?

It is commonly assumed that sex chromosomes evolve recombination suppression because selection favours linkage between sex-determining and sexually antagonistic genes. However, although the role of sexual antagonism during sex chromosome evolution has attained strong support from theory, experimental and observational evidence is rare or equivocal. Here, we highlight alternative, often neglected, hypotheses for recombination suppression on sex chromosomes, which invoke meiotic drive, heterozygote advantage, and genetic drift, respectively. We contrast the hypotheses, the situations when they are likely to be of importance, and outline why it is surprisingly difficult to test them. Lastly, we discuss future research directions (including modelling, population genomics, comparative approaches, and experiments) to disentangle the different hypotheses of sex chromosome evolution.

Whole-genome analysis across 10 songbird families within Sylvioidea reveals a novel autosome-sex chromosome fusion.

Sex chromosomes in birds have long been considered to be extremely stable. However, this notion has lately been challenged by findings of independent autosome-sex chromosome fusions within songbirds, several of which occur within a single clade, the superfamily Sylvioidea. To understand what ecological and evolutionary processes drive changes in sex chromosome systems, we need complete descriptions of sex chromosome diversity across taxonomic groups. Here, we characterize the sex chromosome systems across Sylvioidea using whole-genome data of species representatives of 10 different families, including two published and eight new genomes. We describe a novel fusion in the family Cisticolidae (represented by Cisticola juncidis) involving a part of chromosome 4. We also confirm the previously identified fusion between chromosome Z and a part of chromosome 4A in all 10 families and show that fusions involving parts of chromosomes 3 and 5 are not found outside the families where they were first discovered (Alaudidae and Panuridae). These findings add to the complexity of the sex chromosome system in Sylvioidea, where four independent autosome-sex chromosome fusions have now been identified.

Repeated sex chromosome evolution in vertebrates supported by expanded avian sex chromosomes.

Sex chromosomes have evolved from the same autosomes multiple times across vertebrates, suggesting that selection for recombination suppression has acted repeatedly and independently on certain genetic backgrounds. Here, we perform comparative genomics of a bird clade (larks and their sister lineage; Alaudidae and Panuridae) where multiple autosome-sex chromosome fusions appear to have formed expanded sex chromosomes. We detected the largest known avian sex chromosome (195.3 Mbp) and show that it originates from fusions between parts of four avian chromosomes: Z, 3, 4A and 5. Within these four chromosomes, we found evidence of five evolutionary strata where recombination had been suppressed at different time points, and show that stratum age explained the divergence rate of Z-W gametologs. Next, we analysed chromosome content and found that chromosome 3 was significantly enriched for genes with predicted sex-related functions. Finally, we demonstrate extensive homology to sex chromosomes in other vertebrate lineages: chromosomes Z, 3, 4A and 5 have independently evolved into sex chromosomes in fish (Z), turtles (Z, 5), lizards (Z, 4A), mammals (Z, 4A) and frogs (Z, 3, 4A, 5). Our results provide insights into and support for repeated evolution of sex chromosomes in vertebrates.

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

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

The genomic architecture of the passerine MHC region: High repeat content and contrasting evolutionary histories of single copy and tandemly duplicated MHC genes.

The major histocompatibility complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assembled and annotated the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5 Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologues among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

Biomolecular analyses reveal the age, sex and species identity of a near-intact Pleistocene bird carcass.

Ancient remains found in permafrost represent a rare opportunity to study past ecosystems. Here, we present an exceptionally well-preserved ancient bird carcass found in the Siberian permafrost, along with a radiocarbon date and a reconstruction of its complete mitochondrial genome. The carcass was radiocarbon dated to approximately 44-49 ka BP, and was genetically identified as a female horned lark. This is a species that usually inhabits open habitat, such as the steppe environment that existed in Siberia at the time. This near-intact carcass highlights the potential of permafrost remains for evolutionary studies that combine both morphology and ancient nucleic acids.

Insights into Avian Incomplete Dosage Compensation: Sex-Biased Gene Expression Coevolves with Sex Chromosome Degeneration in the Common Whitethroat.

Non-recombining sex chromosomes (Y and W) accumulate deleterious mutations and degenerate. This poses a problem for the heterogametic sex (XY males; ZW females) because a single functional gene copy often implies less gene expression and a potential imbalance of crucial expression networks. Mammals counteract this by dosage compensation, resulting in equal sex chromosome expression in males and females, whereas birds show incomplete dosage compensation with significantly lower expression in females (ZW). Here, we study the evolution of Z and W sequence divergence and sex-specific gene expression in the common whitethroat (Sylvia communis), a species within the Sylvioidea clade where a neo-sex chromosome has been formed by a fusion between an autosome and the ancestral sex chromosome. In line with data from other birds, females had lower expression than males at the majority of sex-linked genes. Results from the neo-sex chromosome region showed that W gametologs have diverged functionally to a higher extent than their Z counterparts, and that the female-to-male expression ratio correlated negatively with the degree of functional divergence of these gametologs. We find it most likely that sex-linked genes are being suppressed in females as a response to W chromosome degradation, rather than that these genes experience relaxed selection, and thus diverge more, by having low female expression. Overall, our data of this unique avian neo-sex chromosome system suggest that incomplete dosage compensation evolves, at least partly, through gradual accumulation of deleterious mutations at the W chromosome and declining female gene expression.

Contrasting results from GWAS and QTL mapping on wing length in great reed warblers.

A major goal in evolutionary biology is to understand the genetic basis of adaptive traits. In migratory birds, wing morphology is such a trait. Our previous work on the great reed warbler (Acrocephalus arundinaceus) shows that wing length is highly heritable and under sexually antagonistic selection. Moreover, a quantitative trait locus (QTL) mapping analysis detected a pronounced QTL for wing length on chromosome 2, suggesting that wing morphology is partly controlled by genes with large effects. Here, we re-evaluate the genetic basis of wing length in great reed warblers using a genomewide association study (GWAS) approach based on restriction site-associated DNA sequencing (RADseq) data. We use GWAS models that account for relatedness between individuals and include covariates (sex, age and tarsus length). The resulting association landscape was flat with no peaks on chromosome 2 or elsewhere, which is in line with expectations for polygenic traits. Analysis of the distribution of p-values did not reveal biases, and the inflation factor was low. Effect sizes were however not uniformly distributed on some chromosomes, and the Z chromosome had weaker associations than autosomes. The level of linkage disequilibrium (LD) in the population decayed to background levels within c. 1 kbp. There could be several reasons to why our QTL study and GWAS gave contrasting results including differences in how associations are modelled (cosegregation in pedigree vs. LD associations), how covariates are accounted for in the models, type of marker used (multi- vs. biallelic), difference in power or a combination of these. Our study highlights that the genetic architecture even of highly heritable traits is difficult to characterize in wild populations.

Genetic differences between willow warbler migratory phenotypes are few and cluster in large haplotype blocks.

It is well established that differences in migratory behavior between populations of songbirds have a genetic basis but the actual genes underlying these traits remains largely unknown. In an attempt to identify such candidate genes we de novo assembled the genome of the willow warbler Phylloscopus trochilus, and used whole-genome resequencing and a SNP array to associate genomic variation with migratory phenotypes across two migratory divides around the Baltic Sea that separate SW migrating P. t. trochilus wintering in western Africa and SSE migrating P. t. acredula wintering in eastern and southern Africa. We found that the genomes of the two migratory phenotypes lack clear differences except for three highly differentiated regions located on chromosomes 1, 3, and 5 (containing 146, 135, and 53 genes, respectively). Within each migratory phenotype we found virtually no differences in allele frequencies for thousands of SNPs, even when comparing geographically distant populations breeding in Scandinavia and Far East Russia (>6000 km). In each of the three differentiated regions, multidimensional scaling-based clustering of SNP genotypes from more than 1100 individuals demonstrates the presence of distinct haplotype clusters that are associated with each migratory phenotype. In turn, this suggests that recombination is absent or rare between haplotypes, which could be explained by inversion polymorphisms. Whereas SNP alleles on chromosome 3 correlate with breeding altitude and latitude, the allele distribution within the regions on chromosomes 1 and 5 perfectly matches the geographical distribution of the migratory phenotypes. The most differentiated 10 kb windows and missense mutations within these differentiated regions are associated with genes involved in fatty acid synthesis, possibly representing physiological adaptations to the different migratory strategies. The ∼200 genes in these regions, of which several lack described function, will direct future experimental and comparative studies in the search for genes that underlie important migratory traits.