The evolutionary maintenance of ancient recombining sex chromosomes in the ostrich.

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

The association between haemosporidian infection and non-breeding moult location in great reed warblers revisited by combining feather stable isotope profiles and geolocator data.

Stable isotope analysis provides valuable insights into the ecology of long-distance migratory birds during periods spent away from a specific study site. In a previous study, Swedish great reed warblers (Acrocephalus arundinaceus) infected with haemosporidian parasites differed in feather isotope ratios compared to non-infected birds, suggesting that infected and non-infected birds spent the non-breeding season in different locations or habitats. Here, we use a novel dataset comprising geolocator data, isotopes, and haemosporidian infection status of 92 individuals from four Eurasian populations to investigate whether parasite transmission varies with geography or habitats. We found that the probability of harbouring Plasmodium and Leucocytozoon parasites was higher in birds moulting in the eastern region of the non-breeding grounds. However, no geographic pattern occurred for Haemoproteus infections or overall infection status. In contrast to the previous study, we did not find any relationship between feather isotope ratios and overall haemosporidian infection for the entire current dataset. Plasmodium-infected birds had lower feather δ15N values indicating that they occupied more mesic habitats. Leucocytozoon-infected birds had higher feather δ34S values suggesting more coastal sites or wetlands with anoxic sulphate reduction. As the composition and prevalence of haemosporidian parasites differed between the old and the current dataset, we suggest that the differences might be a consequence of temporal dynamics of haemosporidian parasites. Our results emphasize the importance of replicating studies conducted on a single population over a restricted time period, as the patterns can become more complex for data from wider geographical areas and different time periods.

RecView: an interactive R application for locating recombination positions using pedigree data.

BACKGROUND: Recombination reshuffles alleles at linked loci, allowing genes to evolve independently and consequently enhancing the efficiency of selection. This makes quantifying recombination along chromosomes an important goal for understanding how selection and drift are acting on genes and chromosomes. RESULTS: We present RecView, an interactive R application and its homonymous R package, to facilitate locating recombination positions along chromosomes or scaffolds using whole-genome genotype data of a three-generation pedigree. RecView analyses and plots the grandparent-of-origin of all informative alleles along each chromosome of the offspring in the pedigree, and infers recombination positions with either of two built-in algorithms: one based on change in the proportion of the alleles with specific grandparent-of-origin, and one on the degree of continuity of alleles with the same grandparent-of-origin. RecView handles multiple offspring and chromosomes simultaneously, and all putative recombination positions are reported in base pairs together with an estimated precision based on the local density of informative alleles. We demonstrate RecView using genotype data of a passerine bird with an available reference genome, the great reed warbler (Acrocephalus arundinaceus), and show that recombination events can be located to specific positions. CONCLUSIONS: RecView is an easy-to-use and highly effective application for locating recombination positions with high precision. RecView is available on GitHub ( https://github.com/HKyleZhang/RecView.git ).

Animal tracing with sulfur isotopes: Spatial segregation and climate variability in Africa likely contribute to population trends of a migratory songbird.

Climatic conditions affect animals but range-wide impacts at the population level remain largely unknown, especially in migratory species. However, studying climate-population relationships is still challenging in small migrants due to a lack of efficient and cost-effective geographic tracking method. Spatial distribution patterns of environmental stable isotopes (so called 'isoscapes') generally overcome these limitations but none of the currently available isoscapes provide a substantial longitudinal gradient in species-rich sub-Saharan Africa. In this region, sulphur (δ34 S) has not been sufficiently explored on a larger scale. We developed a δ34 S isoscape to trace animal origins in sub-Saharan Africa by coupling known-origin samples from tracked migratory birds with continental remotely sensed environmental data building on environment-δ34 S relationships using a flexible machine learning technique. Furthermore, we link population-specific nonbreeding grounds with interannual climatic variation that might translate to breeding population trends. The predicted δ34 S isotopic map featured east-west and coast-to-inland isotopic gradients and was applied to predict nonbreeding grounds of three breeding populations of Eurasian Reed Warblers Acrocephalus scirpaceus with two distinct migratory phenotypes. Breeding populations as well as migratory phenotypes exhibited large-scale segregation within the African nonbreeding range. These regions also differed substantially in the interannual climatic variation, with higher interannual variability in the eastern part of the range during 2001-2012. Over the same period, the eastern European breeding population seemed to have experienced a more steep decline in population size. The link between migratory patterns and large-scale climatic variability appears important to better understand population trajectories in many declining migratory animals. We believe animal tracing using sulphur isotopes will facilitate these efforts and offers manifold ecological and forensic applications in the biodiversity hotspot of sub-Saharan Africa.

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.

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.

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.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies.

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

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.

Latitudinal clines in sexual selection, sexual size dimorphism and sex-specific genetic dispersal during a poleward range expansion.

Range expansions can be shaped by sex differences in behaviours and other phenotypic traits affecting dispersal and reproduction. Here, we investigate sex differences in morphology, behaviour and genomic population differentiation along a climate-mediated range expansion in the common bluetail damselfly (Ischnura elegans) in northern Europe. We sampled 65 sites along a 583-km gradient spanning the I. elegans range in Sweden and quantified latitudinal gradients in site relative abundance, sex ratio and sex-specific shifts in body size and mating status (a measure of sexual selection). Using single nucleotide polymorphism (SNP) data for 426 individuals from 25 sites, we further investigated sex-specific landscape and climatic effects on neutral genetic connectivity and migration patterns. We found evidence for sex differences associated with the I. elegans range expansion, namely (a) increased male body size with latitude, but no latitudinal effect on female body size, resulting in reduced sexual dimorphism towards the range limit, (b) a steeper decline in male genetic similarity with increasing geographic distance than in females, (c) male-biased genetic migration propensity and (d) a latitudinal cline in migration distance (increasing migratory distances towards the range margin), which was stronger in males. Cooler mean annual temperatures towards the range limit were associated with increased resistance to gene flow in both sexes. Sex ratios became increasingly male biased towards the range limit, and there was evidence for a changed sexual selection regime shifting from favouring larger males in the south to favouring smaller males in the north. Our findings suggest sex-specific spatial phenotype sorting at the range limit, where larger males disperse more under higher landscape resistance associated with cooler climates. The combination of latitudinal gradients in sex-biased dispersal, increasing male body size and (reduced) sexual size dimorphism should have emergent consequences for sexual selection dynamics and the mating system at the expanding range front. Our study illustrates the importance of considering sex differences in the study of range expansions driven by ongoing climate change.

Genome assembly, sex-biased gene expression and dosage compensation in the damselfly Ischnura elegans.

The evolution of sex chromosomes, and patterns of sex-biased gene expression and dosage compensation, are poorly known among early winged insects such as odonates. We assembled and annotated the genome of Ischnura elegans (blue-tailed damselfly), which, like other odonates, has a male-hemigametic sex-determining system (X0 males, XX females). By identifying X-linked genes in I. elegans and their orthologs in other insect genomes, we found homologies between the X chromosome in odonates and chromosomes of other orders, including the X chromosome in Coleoptera. Next, we showed balanced expression of X-linked genes between sexes in adult I. elegans, i.e. evidence of dosage compensation. Finally, among the genes in the sex-determining pathway only fruitless was found to be X-linked, while only doublesex showed sex-biased expression. This study reveals partly conserved sex chromosome synteny and independent evolution of dosage compensation among insect orders separated by several hundred million years of evolutionary history.

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.

Detecting Polygenic Evolution: Problems, Pitfalls, and Promises.

Unraveling the genetic basis of organismal form and function remains one of the major goals of evolutionary biology. Theory has long supported a model of polygenic evolution in which quantitative traits are underpinned by many genes of small effect, but empirical methods have lacked the power to detect causative loci when effect sizes are small or moderate. We (i) review traditional approaches used for identifying the molecular basis of phenotypic traits, to highlight the inherent problems and pitfalls that bias them towards the detection of large-effect loci. We then (ii) outline the promises of recent statistical frameworks to detect polygenic signatures of trait evolution, and discuss some of the first studies in evolutionary biology employing these approaches. Lastly, we (iii) outline future directions and point to areas that still need development.

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 origin of the world's smallest flightless bird, the Inaccessible Island Rail Atlantisia rogersi (Aves: Rallidae).

Rails (Aves: Rallidae) are renowned for their extreme dispersal capability, which has given rise to numerous island lineages. Many insular species lost the ability to fly as a response to release from predator pressure-a feature causing rapid extinction when humans subsequently introduced mammals. The world's smallest extant flightless bird, the Inaccessible Island Rail Atlantisia rogersi, is endemic to Inaccessible Island, Tristan da Cunha archipelago, in the central South Atlantic Ocean. It is placed in a monotypic genus, but its taxonomic affinity, as well as geographic origin, are disputed. Contrary to its suggested Old World origin, we demonstrate that the Inaccessible Island Rail is nested within the mainly South American 'Laterallus clade' and that it colonized ≥3 million-year-old Inaccessible Island from South America c. 1.5 million years ago. The taxonomy of rails has traditionally been based on morphology, and convergent evolution has caused many cases of misclassification. We suggest a re-classification within the 'Laterallus clade' and call for extended coverage of taxon sampling for DNA sequencing.

Evidence for sexual conflict over major histocompatibility complex diversity in a wild songbird.

Sex differences in parasite load and immune responses are found across a wide range of animals, with females generally having lower parasite loads and stronger immune responses than males. Intrigued by these general patterns, we investigated if there was any sign of sex-specific selection on an essential component of adaptive immunity that is known to affect fitness, the major histocompatibility complex class I (MHC-I) genes, in a 20-year study of great reed warblers. Our analyses on fitness related to MHC-I diversity showed a highly significant interaction between MHC-I diversity and sex, where males with higher, and females with lower, MHC-I diversity were more successful in recruiting offspring. Importantly, mean MHC-I diversity did not differ between males and females, and consequently neither sex reached its MHC-I fitness optimum. Thus, there is an unresolved genetic sexual conflict over MHC-I diversity in great reed warblers. Selection from pathogens is known to maintain MHC diversity, but previous theory ignores that the immune environments are considerably different in males and females. Our results suggest that sexually antagonistic selection is an important, previously neglected, force in the evolution of vertebrate adaptive immunity, and have implications for evolutionary understanding of costs of immune responses and autoimmune diseases.

Signatures of local adaptation along environmental gradients in a range-expanding damselfly (Ischnura elegans).

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range-expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site-associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V-ATPase) and visual processes (long-wavelength-sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate-mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.