Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Understanding the chromosomal evolution in cuckoos (Aves, Cuculiformes): a journey through unusual rearrangements

The Cuculiformes are a family of over 150 species that live in a range of habitats, such as forests, savannas, and deserts. Here, bacterial artificial chromosome (BAC) probes (75 from chicken and 14 from zebra finch macrochromosomes 1-10 +ZW and for microchromosomes 11-28 (except 16)) were used to investigate chromosome homologies between chicken and the squirrel cuckoo (Piaya cayana). In addition, repetitive DNA probes were applied to characterize the chromosome organization and to explore the role of these sequences in the karyotype evolution of P. cayana. We also applied BAC probes for chicken chromosome 17 and Z to the guira cuckoo (Guira guira) to test whether this species has an unusual Robertsonian translocation between a microchromosome and the Z chromosome, recently described in the smooth-billed ani (Crotophaga ani). Our results revealed extensive chromosome reorganization with inter- and intrachromosomal rearrangements in P. cayana, including a conspicuous chromosome size and heterochromatin polymorphism on chromosome pair 20. Furthermore, we confirmed that the Z-autosome Robertsonian translocation found in C. ani is also found in G. guira, not P. cayana. These findings suggest that this translocation occurred prior to the divergence between C. ani and G. guira, but after the divergence with P. cayana.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Whole-genome phylogeography of the blue-faced honeyeater (Entomyzon cyanotis) and discovery and characterization of a neo-Z chromosome

Whole-genome surveys of genetic diversity and geographic variation often yield unexpected discoveries of novel structural variation, which long-read DNA sequencing can help clarify. Here, we report on whole-genome phylogeography of a bird exhibiting classic vicariant geographies across Australia and New Guinea, the blue-faced honeyeater (Entomyzon cyanotis), and the discovery and characterization of a novel neo-Z chromosome by long-read sequencing. Using short-read genome-wide SNPs, we inferred population divergence events within E. cyanotis across the Carpentarian and other biogeographic barriers during the Pleistocene (~0.3-1.7 Ma). Evidence for introgression between nonsister populations supports a hypothesis of reticulate evolution around a triad of dynamic barriers around Pleistocene Lake Carpentaria between Australia and New Guinea. During this phylogeographic survey, we discovered a large (134 Mbp) neo-Z chromosome and we explored its diversity, divergence and introgression landscape. We show that, as in some sylvioid passerine birds, a fusion occurred between chromosome 5 and the Z chromosome to form a neo-Z chromosome; and in E. cyanotis, the ancestral pseudoautosomal region (PAR) appears nonrecombinant between Z and W, along with most of the fused chromosome 5. The added recombination-suppressed portion of the neo-Z (~37.2 Mbp) displays reduced diversity and faster population genetic differentiation compared with the ancestral-Z. Yet, the new PAR (~17.4 Mbp) shows elevated diversity and reduced differentiation compared to autosomes, potentially resulting from introgression. In our case, long-read sequencing helped clarify the genomic landscape of population divergence on autosomes and sex chromosomes in a species where prior knowledge of genome structure was still incomplete.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

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.

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.

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.

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots

The karyotype of most birds has remained considerably stable during more than 100 million years’ evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

Parrots have undergone substantial karyotype evolution compared to most other birds. Here, Huang et al. analyze chromosome-level genome assemblies for four parrot species and elucidate the complex evolutionary history of parrot chromosomes.

Neo-sex chromosome evolution and phenotypic differentiation across an elevational gradient in horned larks (Eremophila Alpestris)

Genetic structure and phenotypic variation among populations is affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD dataset of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment.

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.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.

Long-term stability of sex chromosome gene content allows accurate qPCR-based molecular sexing across birds

Embryos, juveniles, and even adults of many bird species lack pronounced external sexually dimorphic characteristics. Accurate identification of sex is crucial for research (e.g., developmental, population, and evolutionary studies), management of wildlife species, and captive breeding programmes for both conservation and poultry. An accurate molecular sexing method applicable across the entire bird radiation is theoretically possible thanks to the long-term stability of their ZZ/ZW sex chromosomes, but current methods are not applicable in a wide range of bird lineages. Here, we developed a novel molecular sexing method based on the comparison of gene copy number variation by quantitative real-time PCR (qPCR) in conserved Z-specific genes (CHRNA6, DDX4, LPAR1, TMEM161B, VPS13A), i.e. genes linked to Z but absent from W chromosomes. We tested the method across three paleognath and 70 neognath species covering the avian phylogeny. In addition, we designed primers for four Z-specific genes (DOCK8, FUT10, PIGG and PSD3) for qPCR-based molecular sexing in three paleognath species. We have demonstrated that the genes DOCK8, FUT10, PIGG and PSD3 can identify sex in paleognath birds and the genes CHRNA6, DDX4, TMEM161B, and VPS13A can reveal sex in neognath birds. The gene LPAR1 can be used to accurately identify sex in both paleognath and neognath species. Along with outlining a novel method of practical importance for molecular sexing in birds, our study also documents in detail the conservation of sex chromosomes across the avian phylogeny.

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single-island endemic bird, but no evidence of MHC-based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC-I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC-I genotype/diversity on mate choice or survival. However, we demonstrate that MHC-I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2-8 linked MHC-I loci. Within-locus homozygosity is high, contributing to low population-wide diversity. Conversely, each individual had comparably many alleles, 6-16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within-individual MHC-I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species' adaptive potential.

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.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.