Species-specific dynamics may cause deviations from general biogeographical predictions - evidence from a population genomics study of a New Guinean endemic passerine bird family (Melampittidae).

The family Melampittidae is endemic to New Guinea and consists of two monotypic genera: Melampitta lugubris (Lesser Melampitta) and Megalampitta gigantea (Greater Melampitta). Both Melampitta species have scattered and disconnected distributions across New Guinea in the central mountain range and in some of the outlying ranges. While M. lugubris is common and found in most montane regions of the island, M. gigantaea is elusive and known from only six localities in isolated pockets on New Guinea with very specific habitats of limestone and sinkholes. In this project, we apply museomics to determine the population structure and demographic history of these two species. We re-sequenced the genomes of all seven known M. gigantaea samples housed in museum collections as well as 24 M. lugubris samples from across its distribution. By comparing population structure between the two species, we investigate to what extent habitat dependence, such as in M. gigantaea, may affect population connectivity. Phylogenetic and population genomic analyses, as well as acoustic variation revealed that M. gigantaea consists of a single population in contrast to M. lugubris that shows much stronger population structure across the island. We suggest a recent collapse of M. gigantaea into its fragmented habitats as an explanation to its unexpected low diversity and lack of population structure. The deep genetic divergences between the M. lugubris populations on the Vogelkop region, in the western central range and the eastern central range, respectively, suggests that these three populations should be elevated to full species level. This work sheds new light on the mechanisms that have shaped the intriguing distribution of the two species within this family and is a prime example of the importance of museum collections for genomic studies of poorly known and rare species.

The formation of the Indo-Pacific montane avifauna.

The processes generating the earth's montane biodiversity remain a matter of debate. Two contrasting hypotheses have been advanced to explain how montane populations form: via direct colonization from other mountains, or, alternatively, via upslope range shifts from adjacent lowland areas. We seek to reconcile these apparently conflicting hypotheses by asking whether a species' ancestral geographic origin determines its mode of mountain colonization. Island-dwelling passerine birds at the faunal crossroads between Eurasia and Australo-Papua provide an ideal study system. We recover the phylogenetic relationships of the region's montane species and reconstruct their ancestral geographic ranges, elevational ranges, and migratory behavior. We also perform genomic population studies of three super-dispersive montane species/clades with broad island distributions. Eurasian-origin species populated archipelagos via direct colonization between mountains. This mode of colonization appears related to ancestral adaptations to cold and seasonal climates, specifically short-distance migration. Australo-Papuan-origin mountain populations, by contrast, evolved from lowland ancestors, and highland distribution mostly precludes their further colonization of island mountains. Our study explains much of the distributional variation within a complex biological system, and provides a synthesis of two seemingly discordant hypotheses for montane community formation.

Introgression underlies phylogenetic uncertainty but not parallel plumage evolution in a recent songbird radiation.

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails", a group of five species in the songbird genus Motacilla: one species, Motacilla alba, shows wide intra-specific plumage variation, while the four others form two pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the two species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome, to quantify the extent of discordance across genomic regions, and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investigation of genome-wide patterns of gene tree heterogeneity to help understanding the mechanisms underlying phenotypic evolution.

Megafaunal extinctions, not climate change, may explain Holocene genetic diversity declines in Numenius shorebirds.

Understanding the relative contributions of historical and anthropogenic factors to declines in genetic diversity is important for informing conservation action. Using genome-wide DNA of fresh and historic specimens, including that of two species widely thought to be extinct, we investigated fluctuations in genetic diversity and present the first complete phylogenomic tree for all nine species of the threatened shorebird genus Numenius, known as whimbrels and curlews. Most species faced sharp declines in effective population size, a proxy for genetic diversity, soon after the Last Glacial Maximum (around 20,000 years ago). These declines occurred prior to the Anthropocene and in spite of an increase in the breeding area predicted by environmental niche modeling, suggesting that they were not caused by climatic or recent anthropogenic factors. Crucially, these genetic diversity declines coincide with mass extinctions of mammalian megafauna in the Northern Hemisphere. Among other factors, the demise of ecosystem-engineering megafauna which maintained open habitats may have been detrimental for grassland and tundra-breeding Numenius shorebirds. Our work suggests that the impact of historical factors such as megafaunal extinction may have had wider repercussions on present-day population dynamics of open habitat biota than previously appreciated.

About 20,000 years ago, the Earth was a much colder world roamed by giant mastodons, gigantic elks, woolly mammoths and sabre-tooth tigers. Yet these imposing creatures were living on borrowed time: by the start of the Holocene, around 10,000 years later, many animals over 45kg had vanished across the Northern Hemisphere, closing the book on what is known as the Quaternary extinction event. As large grazers disappeared, the landscape likely changed too. Where open tundra and grasslands may have once dominated, woodlands and shrubs probably took over, creating ripple effects for surviving species. These extinction events took place in a changing world, with glaciers starting to retreat about 20,000 years ago and human populations colonizing an increasing share of this planet's land area. In fact, since the end of this last glacial maximum, ecosystems have been reshaped by a succession and a combination of climatic, historical and human-driven forces. This makes it difficult for scientists to disentangle the relative contribution of these factors on the lives of animals. Tan et al. decided to explore this question by reconstructing how effective population sizes changed over the past 20,000 years for nine species of curlews and whimbrels. These shorebirds, which together comprise the genus Numenius, breed slowly and nest in open environments such as moorlands or tundra. Many are currently under threat. Fluctuations in the numbers of breeding individuals affect the genetic diversity of a species, and these events leave tell-tale genetic signatures that can be uncovered through DNA analyses. Tan et al. had enough fresh and museum samples to infer these changes for five Numenius species, revealing that genetic diversity brutally dropped soon after the last glacial period ended. At the time, humans were yet to make significant changes on their environment and a warming world should have supported population growth. Tan et al. suggest that, instead, this sharp decline is linked to the late Quaternary extinctions of large mammals: with the demise of grazing animals which could keep woodlands at bay, the shorebirds lost their open nesting grounds. This event has left its mark in the genome of existing species, with these birds still exhibiting a low level of genetic diversity that may put them at further risk for extinction.

Patterns of phylogenetic diversification in the Dollarbird (Eurystomus orientalis) and Azure Roller (Eurystomus azureus) complex.

Genetic isolation and morphological differentiation are two important factors in the speciation process that not always act in concert. A rapid morphological change in a lineage can hide its close relationship to another lineage, while slight morphological differentiation between two taxa can give the appearance of a closer relationship than is actually the case. The Dollarbird (Eurystomus orientalis) and the Azure Roller (Eurystomus azureus) is such an example. Today the Dollarbird and the Azure Roller are unanimously considered to constitute two distinct species, but in a recent genetic study it has been shown that the latter taxon, despite being larger and having a distinctly different coloration, is phylogenetically nested within the former. Its precise placement within this complex has not been determined, however. In this study, we investigate the phylogenetic relationships within the Dollarbird/Azure Roller complex. We estimate divergence times and infer phylogenetic relationships using sequence data from 6,475 genome-wide intronic regions, as well as complete mitochondrial genomes, using both concatenation and multispecies coalescence approaches. We find that within the Dollarbird/Azure Roller complex there are several examples of discrepancies between genetic and morphological differentiation. The Dollarbird is currently divided into between nine to twelve subspecies. Some of these subspecies are poorly differentiated, whereas others are morphologically more clearly discernable. Our data suggest that the complex consist of at least seven distinct genetic lineages that do not entirely match the morphological variation within the group. For instance, our results show that the subspecies solomonensis from the Solomon Islands, despite being morphologically very similar to its geographically closest neighbors, in fact is a highly distinct lineage that became isolated more than 700,000 years ago. In contrast, the morphologically distinct Azure Roller, which is currently treated as a distinct species, is nested within the Dollarbird and forms a slightly younger lineage than solomonensis and is the sister group to a clade with Australian and New Guinean Dollarbirds. Our results also show a deep genetic split within the Dollarbirds on the Asian mainland. This stands in contrast to the apparent clinal morphological variation reported for the birds on the Asian mainland. We also find support for the presence of a genetically distinct clade in the Wallacea region. The birds from the Wallacea region has previously been recognized as a distinct subspecies, connectens, but is currently placed in synonymy of other subspecies. Our results are thus at odds with the current division of the Dollarbird/Azure Roller complex into two species. Given that the species status of azureus is undisputed, the apparent genetic isolation of solomonensis and its clear separation from the other lineages suggests that this taxon also warrants species status. Based on the genetic and morphological variation observed within the Dollarbird/Azure Roller complex there is little doubt that even more taxa should regarded as species, but this require further examination.

Population genomics of the island thrush elucidates one of earth's great archipelagic radiations.

Tropical islands are renowned as natural laboratories for evolutionary study. Lineage radiations across tropical archipelagos are ideal systems for investigating how colonization, speciation, and extinction processes shape biodiversity patterns. The expansion of the island thrush across the Indo-Pacific represents one of the largest yet most perplexing island radiations of any songbird species. The island thrush exhibits a complex mosaic of pronounced plumage variation across its range and is arguably the world's most polytypic bird. It is a sedentary species largely restricted to mountain forests, yet it has colonized a vast island region spanning a quarter of the globe. We conducted a comprehensive sampling of island thrush populations and obtained genome-wide SNP data, which we used to reconstruct its phylogeny, population structure, gene flow, and demographic history. The island thrush evolved from migratory Palearctic ancestors and radiated explosively across the Indo-Pacific during the Pleistocene, with numerous instances of gene flow between populations. Its bewildering plumage variation masks a biogeographically intuitive stepping stone colonization path from the Philippines through the Greater Sundas, Wallacea, and New Guinea to Polynesia. The island thrush's success in colonizing Indo-Pacific mountains can be understood in light of its ancestral mobility and adaptation to cool climates; however, shifts in elevational range, degree of plumage variation and apparent dispersal rates in the eastern part of its range raise further intriguing questions about its biology.

Museomics reveals the phylogenetic position of the extinct Moroccan trout Salmo pallaryi.

The authors used museomics to reconstruct the mitochondrial genome from two individuals of the Moroccan, endemic and extinct trout, Salmo pallaryi. They further obtained partial data from 21 nuclear genes previously used for trout phylogenetic analyses. Phylogenetic analyses, including publicly available data from the mitochondrial control region and the cytochrome b gene, and the 21 nuclear genes, place S. pallaryi among other North African trouts. mtDNA places S. pallaryi close to Salmo macrostigma within a single North African clade. Although the nuclear coverage of the genome was low, both specimens were independently positioned as sisters to one of two distantly related North African clades, viz. the Atlas clade with the Dades trout, Salmo multipunctatus. Phylogenetic discordance between mtDNA and nuclear DNA phylogenies is briefly discussed. As several specimens that were extracted failed to produce DNA of sufficient quality, the authors discuss potential reasons for the failure. They suggest that museum specimens in poor physical condition may be better for DNA extraction compared to better-preserved ones, possibly related to the innovation of formalin as a fixative before ethanol storage in the early 20th century.

Satellite DNA evolution in Corvoidea inferred from short and long reads.

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

A 14,000-year-old genome sheds light on the evolution and extinction of a Pleistocene vulture.

The New World Vulture [Coragyps] occidentalis (L. Miller, 1909) is one of many species that were extinct by the end of the Pleistocene. To understand its evolutionary history we sequenced the genome of a 14,000 year old [Coragyps] occidentalis found associated with megaherbivores in the Peruvian Andes. occidentalis has been viewed as the ancestor, or possibly sister, to the extant Black Vulture Coragyps atratus, but genomic data shows occidentalis to be deeply nested within the South American clade of atratus. Coragyps atratus inhabits lowlands, but the fossil record indicates that occidentalis mostly occupied high elevations. Our results suggest that occidentalis evolved from a population of atratus in southwestern South America that colonized the High Andes 300 to 400 kya. The morphological and morphometric differences between occidentalis and atratus may thus be explained by ecological diversification following from the natural selection imposed by this new and extreme, high elevation environment. The sudden evolution of a population with significantly larger body size and different anatomical proportions than atratus thus constitutes an example of punctuated evolution.

A guide to avian museomics: Insights gained from resequencing hundreds of avian study skins.

Biological specimens in natural history collections constitute a massive repository of genetic information. Many specimens have been collected in areas in which they no longer exist or in areas where present-day collecting is not possible. There are also specimens in collections representing populations or species that have gone extinct. Furthermore, species or populations may have been sampled throughout an extensive time period, which is particularly valuable for studies of genetic change through time. With the advent of high-throughput sequencing, natural history museum resources have become accessible for genomic research. Consequently, these unique resources are increasingly being used across many fields of natural history. In this paper, we summarize our experiences of resequencing hundreds of genomes from historical avian museum specimens. We publish the protocols we have used and discuss the entire workflow from sampling and laboratory procedures, to the bioinformatic processing of historical specimen data.

Speciation and population divergence in a mutualistic seed dispersing bird.

Bird-mediated seed dispersal is crucial for the regeneration and viability of ecosystems, often resulting in complex mutualistic species networks. Yet, how this mutualism drives the evolution of seed dispersing birds is still poorly understood. In the present study we combine whole genome re-sequencing analyses and morphometric data to assess the evolutionary processes that shaped the diversification of the Eurasian nutcracker (Nucifraga), a seed disperser known for its mutualism with pines (Pinus). Our results show that the divergence and phylogeographic patterns of nutcrackers resemble those of other non-mutualistic passerine birds and suggest that their early diversification was shaped by similar biogeographic and climatic processes. The limited variation in foraging traits indicates that local adaptation to pines likely played a minor role. Our study shows that close mutualistic relationships between bird and plant species might not necessarily act as a primary driver of evolution and diversification in resource-specialized birds.

The formation of avian montane diversity across barriers and along elevational gradients.

Tropical mountains harbor exceptional concentrations of Earth's biodiversity. In topographically complex landscapes, montane species typically inhabit multiple mountainous regions, but are absent in intervening lowland environments. Here we report a comparative analysis of genome-wide DNA polymorphism data for population pairs from eighteen Indo-Pacific bird species from the Moluccan islands of Buru and Seram and from across the island of New Guinea. We test how barrier strength and relative elevational distribution predict population differentiation, rates of historical gene flow, and changes in effective population sizes through time. We find population differentiation to be consistently and positively correlated with barrier strength and a species' altitudinal floor. Additionally, we find that Pleistocene climate oscillations have had a dramatic influence on the demographics of all species but were most pronounced in regions of smaller geographic area. Surprisingly, even the most divergent taxon pairs at the highest elevations experience gene flow across barriers, implying that dispersal between montane regions is important for the formation of montane assemblages.

Utilizing museomics to trace the complex history and species boundaries in an avian-study system of conservation concern.

A taxonomic classification that accurately captures evolutionary history is essential for conservation. Genomics provides powerful tools for delimiting species and understanding their evolutionary relationships. This allows for a more accurate and detailed view on conservation status compared with other, traditionally used, methods. However, from a practical and ethical perspective, gathering sufficient samples for endangered taxa may be difficult. Here, we use museum specimens to trace the evolutionary history and species boundaries in an Asian oriole clade. The endangered silver oriole has long been recognized as a distinct species based on its unique coloration, but a recent study suggested that it might be nested within the maroon oriole-species complex. To evaluate species designation, population connectivity, and the corresponding conservation implications, we assembled a de novo genome and used whole-genome resequencing of historical specimens. Our results show that the silver orioles form a monophyletic lineage within the maroon oriole complex and that maroon and silver forms continued to interbreed after initial divergence, but do not show signs of recent gene flow. Using a genome scan, we identified genes that may form the basis for color divergence and act as reproductive barriers. Taken together, our results confirm the species status of the silver oriole and highlight that taxonomic revision of the maroon forms is urgently needed. Our study demonstrates how genomics and Natural History Collections (NHC) can be utilized to shed light on the taxonomy and evolutionary history of natural populations and how such insights can directly benefit conservation practitioners when assessing wild populations.

Novel genome reveals susceptibility of popular gamebird, the red-legged partridge (Alectoris rufa, Phasianidae), to climate change.

We produced a high-quality de novo genome assembly of the red-legged partridge A. rufa, the first reference genome of its genus, by utilising novel 10× Chromium technology. The estimated genome size was 1.19 Gb with an overall genome heterozygosity of 0.0022; no runs of homozygosity were observed. In total, 21,589 protein coding genes were identified and assigned to 16,772 orthologs. Of these, 201 emerged as unique to Alectoris and were enriched for positive regulation of epithelial cell migration, viral genome integration and maturation. Using PSMC analysis, we inferred a major demographic decline commencing ~140,000 years ago, consistent with forest expansion and reduction of open habitats during the Eemian interglacial. Present-day populations exhibit the historically lowest genetic diversity. Besides implications for management and conservation, this genome also promises key insights into the physiology of these birds with a view to improving poultry husbandry practices.

Complex population structure of the Atlantic puffin revealed by whole genome analyses.

The factors underlying gene flow and genomic population structure in vagile seabirds are notoriously difficult to understand due to their complex ecology with diverse dispersal barriers and extensive periods at sea. Yet, such understanding is vital for conservation management of seabirds that are globally declining at alarming rates. Here, we elucidate the population structure of the Atlantic puffin (Fratercula arctica) by assembling its reference genome and analyzing genome-wide resequencing data of 72 individuals from 12 colonies. We identify four large, genetically distinct clusters, observe isolation-by-distance between colonies within these clusters, and obtain evidence for a secondary contact zone. These observations disagree with the current taxonomy, and show that a complex set of contemporary biotic factors impede gene flow over different spatial scales. Our results highlight the power of whole genome data to reveal unexpected population structure in vagile marine seabirds and its value for seabird taxonomy, evolution and conservation.

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)'.

Genomic signatures of rapid adaptive divergence in a tropical montane species.

Mountain regions contain extraordinary biodiversity. The environmental heterogeneity and glacial cycles often accelerate speciation and adaptation of montane species, but how these processes influence the genomic differentiation of these species is largely unknown. Using a novel chromosome-level genome and population genomic comparisons, we study allopatric divergence and selection in an iconic bird living in a tropical mountain region in New Guinea, Archbold's bowerbird (Amblyornis papuensis). Our results show that the two populations inhabiting the eastern and western Central Range became isolated ca 11 800 years ago, probably because the suitable habitats for this cold-tolerating bird decreased when the climate got warmer. Our genomic scans detect that genes in highly divergent genomic regions are over-represented in developmental processes, which is probably associated with the observed differences in body size between the populations. Overall, our results suggest that environmental differences between the eastern and western Central Range probably drive adaptive divergence between them.

Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise.

Genome assemblies are currently being produced at an impressive rate by consortia and individual laboratories. The low costs and increasing efficiency of sequencing technologies now enable assembling genomes at unprecedented quality and contiguity. However, the difficulty in assembling repeat-rich and GC-rich regions (genomic "dark matter") limits insights into the evolution of genome structure and regulatory networks. Here, we compare the efficiency of currently available sequencing technologies (short/linked/long reads and proximity ligation maps) and combinations thereof in assembling genomic dark matter. By adopting different de novo assembly strategies, we compare individual draft assemblies to a curated multiplatform reference assembly and identify the genomic features that cause gaps within each assembly. We show that a multiplatform assembly implementing long-read, linked-read and proximity sequencing technologies performs best at recovering transposable elements, multicopy MHC genes, GC-rich microchromosomes and the repeat-rich W chromosome. Telomere-to-telomere assemblies are not a reality yet for most organisms, but by leveraging technology choice it is now possible to minimize genome assembly gaps for downstream analysis. We provide a roadmap to tailor sequencing projects for optimized completeness of both the coding and noncoding parts of nonmodel genomes.

Phylogenomics of white-eyes, a 'great speciator', reveals Indonesian archipelago as the center of lineage diversity.

Archipelagoes serve as important 'natural laboratories' which facilitate the study of island radiations and contribute to the understanding of evolutionary processes. The white-eye genus Zosterops is a classical example of a 'great speciator', comprising c. 100 species from across the Old World, most of them insular. We achieved an extensive geographic DNA sampling of Zosterops by using historical specimens and recently collected samples. Using over 700 genome-wide loci in conjunction with coalescent species tree methods and gene flow detection approaches, we untangled the reticulated evolutionary history of Zosterops, which comprises three main clades centered in Indo-Africa, Asia, and Australasia, respectively. Genetic introgression between species permeates the Zosterops phylogeny, regardless of how distantly related species are. Crucially, we identified the Indonesian archipelago, and specifically Borneo, as the major center of diversity and the only area where all three main clades overlap, attesting to the evolutionary importance of this region.