Temporal Dynamics of Avian Populations during Pleistocene Revealed by Whole-Genome Sequences

Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana)

The nonrandom distribution of chromosomal characteristics and functional elements-genomic architecture-impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure-ordination, clustering, and phylogenetic methods-we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

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.

Whole genomes show contrasting trends of population size changes and genomic diversity for an Amazonian endemic passerine over the late quaternary

The "Amazon tipping point" is a global change scenario resulting in replacement of upland terra-firme forests by large-scale "savannization" of mostly southern and eastern Amazon. Reduced rainfall accompanying the Last Glacial Maximum (LGM) has been proposed to have acted as such a tipping point in the past, with the prediction that terra-firme inhabiting species should have experienced reductions in population size as drier habitats expanded. Here, we use whole-genomes of an Amazonian endemic organism (Scale-backed antbirds - Willisornis spp.) sampled from nine populations across the region to test this historical demography scenario. Populations from southeastern Amazonia and close to the Amazon-Cerrado ecotone exhibited a wide range of demographic patterns, while most of those from northern and western Amazonia experienced uniform expansions between 400 kya and 80-60 kya, with gradual declines toward 20 kya. Southeastern populations of Willisornis were the last to diversify and showed smaller heterozygosity and higher runs of homozygosity values than western and northern populations. These patterns support historical population declines throughout the Amazon that affected more strongly lineages in the southern and eastern areas, where historical "tipping point" conditions existed due to the widespread replacement of humid forest by drier and open vegetation during the LGM.

The role of biogeographical barriers on the historical dynamics of passerine birds with a circum-Amazonian distribution

Common distributional patterns have provided the foundations of our knowledge of Neotropical biogeography. A distinctive pattern is the "circum-Amazonian distribution", which surrounds Amazonia across the forested lowlands south and east of the basin, the Andean foothills, the Venezuelan Coastal Range, and the Tepuis. The underlying evolutionary and biogeographical mechanisms responsible for this widespread pattern of avian distribution have yet to be elucidated. Here, we test the effects of biogeographical barriers in four species in the passerine family Thamnophilidae by performing comparative demographic analyses of genome-scale data. Specifically, we used flanking regions of ultraconserved regions to estimate population historical parameters and genealogical trees and tested demographic models reflecting contrasting biogeographical scenarios explaining the circum-Amazonian distribution. We found that taxa with circum-Amazonian distribution have at least two main phylogeographical clusters: (1) Andes, often extending into Central America and the Tepuis; and (2) the remaining of their distribution. These clusters are connected through corridors along the Chaco-Cerrado and southeastern Amazonia, allowing gene flow between Andean and eastern South American populations. Demographic histories are consistent with Pleistocene climatic fluctuations having a strong influence on the diversification history of circum-Amazonian taxa, Refugia played a crucial role, enabling both phenotypic and genetic differentiation, yet maintaining substantial interconnectedness to keep considerable levels of gene flow during different dry/cool and warm/humid periods. Additionally, steep environmental gradients appear to play a critical role in maintaining both genetic and phenotypic structure.

Demographic responses of oceanic island birds to local and regional ecological disruptions revealed by whole-genome sequencing

Disentangling the effects of ecological disruptions operating at different spatial and temporal scales in shaping past species' demography is particularly important in the current context of rapid environmental changes driven by both local and regional factors. We argue that volcanic oceanic islands provide useful settings to study the influence of past ecological disruptions operating at local and regional scales on population demographic histories. We investigate potential drivers of past population dynamics for three closely related species of passerine birds from two volcanic oceanic islands, Reunion and Mauritius (Mascarene archipelago), with distinct volcanic history. Using ABC and PSMC inferences from complete genomes, we reconstructed the demographic history of the Reunion Grey White-eye (Zosterops borbonicus (Pennant, 1781)), the Reunion Olive White-eye (Z. olivaceus (Linnaeus, 1766)) and the Mauritius Grey White-eye (Z. mauritianus (Gmelin, 1789)) and searched for possible causes underlying similarities or differences between species living on the same or different islands. Both demographic inferences strongly support ancient and long-term expansions in all species. They also reveal different trajectories between species inhabiting different islands, but consistent demographic trajectories in species or populations from the same island. Species from Reunion appear to have experienced synchronous reductions in population size during the Last Glacial Maximum, a trend not seen in the Mauritian species. Overall, this study suggests that local events may have played a role in shaping population trajectories of these island species. It also highlights the potential of our conceptual framework to disentangle the effects of local and regional drivers on past species' demography and long-term population processes.

Hybrid de novo genome assembly of the sexually dimorphic Lady Amherst's pheasant

Pheasants are an important group of birds, valued for their economic benefit as poultry birds, game birds, and as ornamental species for their plumage. Lady Amherst's pheasant Chrysolophus amherstiae is an ornamental species, valued for its elaborate and beautiful plumage. In this study, we present a high-quality de novo hybrid genome assembly of C. amherstiae. Previous attempts to sequence the genome of this species resulted in draft-level assemblies, which are not available in the public domain. Using a combination of Illumina short reads and Oxford Nanopore's long-reads, we assembled a high-quality genome of N50 ~3.9 Mb and near complete BUSCO assessment. We observed a correlation between effective population size and past climatic conditions, with an increase in population size during the warm interglacial periods. We further observed significant fluctuations in genes involved with the immune system and visual perception. C. amherstiae is a highly dimorphic species, and significant fluctuations in gene families involved in immune response, visual perception, among others, suggesting a role of mate choice and sexual selection in the evolution and maintenance of exaggerated traits in the males.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

Habitat specialization predicts demographic response and vulnerability of floodplain birds in Amazonia

The annual flooding cycle of Amazonian rivers sustains the largest floodplains on Earth, which harbour a unique bird community. Recent studies suggest that habitat specialization drove different patterns of population structure and gene flow in floodplain birds. However, the lack of a direct estimate of habitat affinity prevents a proper test of its effects on population histories. In this work, we used occurrence data, satellite images and genomic data (ultra-conserved elements) from 24 bird species specialized on a variety of seasonally flooded environments to classify habitat affinities and test its influence on evolutionary histories of Amazonian floodplain birds. We demonstrate that birds with higher specialization in river islands and dynamic environments have gone through more recent demographic expansion and currently have less genetic diversity than floodplain generalist birds. Our results indicate that there is an intrinsic relationship between habitat affinity and environmental dynamics, influencing patterns of population structure, demographic history and genetic diversity. Within the floodplains, historical landscape changes have had more severe impacts on island specialists, making them more vulnerable to current and future anthropogenic changes, as those imposed by hydroelectric dams in the Amazon Basin.

Past volcanic activity predisposes an endemic threatened seabird to negative anthropogenic impacts

Humans are regularly cited as the main driver of current biodiversity extinction, but the impact of historic volcanic activity is often overlooked. Pre-human evidence of wildlife abundance and diversity are essential for disentangling anthropogenic impacts from natural events. Réunion Island, with its intense and well-documented volcanic activity, endemic biodiversity, long history of isolation and recent human colonization, provides an opportunity to disentangle these processes. We track past demographic changes of a critically endangered seabird, the Mascarene petrel Pseudobulweria aterrima, using genome-wide SNPs. Coalescent modeling suggested that a large ancestral population underwent a substantial population decline in two distinct phases, ca. 125,000 and 37,000 years ago, coinciding with periods of major eruptions of Piton des Neiges. Subsequently, the ancestral population was fragmented into the two known colonies, ca. 1500 years ago, following eruptions of Piton de la Fournaise. In the last century, both colonies declined significantly due to anthropogenic activities, and although the species was initially considered extinct, it was rediscovered in the 1970s. Our findings suggest that the current conservation status of wildlife on volcanic islands should be firstly assessed as a legacy of historic volcanic activity, and thereafter by the increasing anthropogenic impacts, which may ultimately drive species towards extinction.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

The genome of a globally invasive passerine, the common myna, Acridotheres tristis

In an era of global climate change, biodiversity conservation is receiving increased attention. Conservation efforts are greatly aided by genetic tools and approaches, which seek to understand patterns of genetic diversity and how they impact species health and their ability to persist under future climate regimes. Invasive species offer vital model systems in which to investigate questions regarding adaptive potential, with a particular focus on how changes in genetic diversity and effective population size interact with novel selection regimes. The common myna (Acridotheres tristis) is a globally invasive passerine and is an excellent model species for research both into the persistence of low-diversity populations and the mechanisms of biological invasion. To underpin research on the invasion genetics of this species, we present the genome assembly of the common myna. We describe the genomic landscape of this species, including genome wide allelic diversity, methylation, repeats, and recombination rate, as well as an examination of gene family evolution. Finally, we use demographic analysis to identify that some native regions underwent a dramatic population increase between the two most recent periods of glaciation, and reveal artefactual impacts of genetic bottlenecks on demographic analysis.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

Amazonian colonization from the Atlantic Forest: New perspectives on the connections of South American tropical forests

An open and dry vegetation belt separates Amazonia (AM) and the Atlantic Forest (AF). Evidence from palaeoclimatic and phylogenetic studies suggests past connections between these forests during cycles of increased humidity through the formation of forest corridors. The distinctive northern AF avifauna is known to have affinities both with AM and the southern AF. Still, the extent of how these two regions contributed to the assemblage of this avifauna remains poorly understood. Using historical demographic analyses and comparative phylogeography based on sub-genomic genetic sampling, we assessed how past connections between AM and AF led to shared vicariance and colonization events in four avian AF endemic taxa. Our results supported the occurrence of humid forest corridors promoting the contact between AF and AM populations and suggested two vicariant events and two colonization events from AF to AM. Population divergences were mostly non-synchronous and occurred multiple times during the Pleistocene. Historical gene flow was prevalent across study groups, supporting migration flows after the initial separation between AM and AF - a pattern previously unknown in birds between these regions. Idiosyncratic histories and divergent demographic syndromes suggest that organisms' responses to climate-driven habitat shifts broadly depend on their ecological attributes. This study strengthened our knowledge of past connections between AM and AF and provided demographic scenarios amenable for testing in other groups of co-distributed organisms.

Genetic Load and Adaptive Potential of a Recovered Avian Species that Narrowly Avoided Extinction

High genetic diversity is a good predictor of long-term population viability, yet some species persevere despite having low genetic diversity. Here we study the genomic erosion of the Seychelles paradise flycatcher (Terpsiphone corvina), a species that narrowly avoided extinction after having declined to 28 individuals in the 1960s. The species recovered unassisted to over 250 individuals in the 1990s and was downlisted from Critically Endangered to Vulnerable in the International Union for the Conservation of Nature Red List in 2020. By comparing historical, prebottleneck (130+ years old) and modern genomes, we uncovered a 10-fold loss of genetic diversity. Highly deleterious mutations were partly purged during the bottleneck, but mildly deleterious mutations accumulated. The genome shows signs of historical inbreeding during the bottleneck in the 1960s, but low levels of recent inbreeding after demographic recovery. Computer simulations suggest that the species long-term small Ne reduced the masked genetic load and made the species more resilient to inbreeding and extinction. However, the reduction in genetic diversity due to the chronically small Ne and the severe bottleneck is likely to have reduced the species adaptive potential to face environmental change, which together with a higher load, compromises its long-term population viability. Thus, small ancestral Ne offers short-term bottleneck resilience but hampers long-term adaptability to environmental shifts. In light of rapid global rates of population decline, our work shows that species can continue to suffer the effect of their decline even after recovery, highlighting the importance of considering genomic erosion and computer modeling in conservation assessments.

Redefining the Evolutionary History of the Rock Dove, Columba livia, Using Whole Genome Sequences

The domestic pigeon's exceptional phenotypic diversity was key in developing Darwin's Theory of Evolution and establishing the concept of artificial selection. However, unlike its domestic counterpart, its wild progenitor, the rock dove Columba livia has received considerably less attention. Therefore, questions regarding its domestication, evolution, taxonomy, and conservation status remain unresolved. We generated whole-genome sequencing data from 65 historical rock doves that represent all currently recognized subspecies and span the species' original geographic distribution. Our dataset includes 3 specimens from Darwin's collection, and the type specimens of 5 different taxa. We characterized their population structure, genomic diversity, and gene-flow patterns. Our results show the West African subspecies C. l. gymnocyclus is basal to rock doves and domestic pigeons, and suggests gene-flow between the rock dove's sister species C. rupestris, and the ancestor of rock doves after its split from West African populations. These genomes allowed us to propose a model for the evolution of the rock dove in light of the refugia theory. We propose that rock dove genetic diversity and introgression patterns derive from a history of allopatric cycles and dispersion waves during the Quaternary glacial and interglacial periods. To explore the rock dove domestication history, we combined our new dataset with available genomes from domestic pigeons. Our results point to at least 1 domestication event in the Levant that gave rise to all domestic breeds analysed in this study. Finally, we propose a species-level taxonomic arrangement to reflect the evolutionary history of the West African rock dove populations.

Evidence of Site-Specific and Male-Biased Germline Mutation Rate in a Wild Songbird

Germline mutations are the ultimate source of genetic variation and the raw material for organismal evolution. Despite their significance, the frequency and genomic locations of mutations, as well as potential sex bias, are yet to be widely investigated in most species. To address these gaps, we conducted whole-genome sequencing of 12 great reed warblers (Acrocephalus arundinaceus) in a pedigree spanning 3 generations to identify single-nucleotide de novo mutations (DNMs) and estimate the germline mutation rate. We detected 82 DNMs within the pedigree, primarily enriched at CpG sites but otherwise randomly located along the chromosomes. Furthermore, we observed a pronounced sex bias in DNM occurrence, with male warblers exhibiting three times more mutations than females. After correction for false negatives and adjusting for callable sites, we obtained a mutation rate of 7.16 × 10-9 mutations per site per generation (m/s/g) for the autosomes and 5.10 × 10-9 m/s/g for the Z chromosome. To demonstrate the utility of species-specific mutation rates, we applied our autosomal mutation rate in models reconstructing the demographic history of the great reed warbler. We uncovered signs of drastic population size reductions predating the last glacial period (LGP) and reduced gene flow between western and eastern populations during the LGP. In conclusion, our results provide one of the few direct estimates of the mutation rate in wild songbirds and evidence for male-driven mutations in accordance with theoretical expectations.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

Chromosomal-level reference genome of a wild North American mallard (Anas platyrhynchos)

The mallard (Anas platyrhynchos) is one of the most common, economically, and socially important birds around the world. Mallards were not only an important food source for early humans but eventually becoming intimately linked with people as they were domesticated over the last 2,000 years. To date, mallard genomes are largely reconstructed from samples of domestic or unknown genetic heritage. Here, we report the first high-quality genome assembly and annotation of a genetically vetted wild mallard from North America (NAwild_v1.0). The genome was assembled using a combination of shotgun libraries, proximity ligation Chicago, and Dovetail Hi-C libraries. The final assembly is ∼1.04 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome-level scaffolds, and with a N50/L50 of 79.1 Mb/4 scaffolds. We used a combination of gene prediction and similarity approaches to annotate a total of 23,584 functional genes, of which 19,242 were associated to GO terms. The genome assembly and the set of annotated genes yielded a 95.4% completeness score when compared with the BUSCO aves_odb10 dataset. Next, we aligned 3 previously published mallard genomes to ours, and demonstrate how runs of homozygosity and nucleotide diversity are substantially higher and lower, respectively, to ours and how these artificially changed genomes resulted in profoundly different and biased demographic histories. Our wild mallard assembly not only provides a valuable resource to shed light onto genome evolution, speciation, and other adaptive processes, but also helping with identifying functional genes that have been significantly altered during the domestication process.

No evidence of inbreeding depression despite a historical severe bottleneck in the endangered Bermuda petrel (Pterodroma cahow)

The Bermuda petrel Pterodroma cahow is an island endemic seabird that belongs to the Procellariiformes, one of the most endangered orders of birds. Historical records suggest a significant population size decline following human settlement in Bermuda, bringing the species to near extinction. Since the 1950s, the population has been recovering aided by the implementation of an ongoing conservation plan. However, it still faces several threats, and negative genetic effects resulting from that drastic decline are to be expected, including inbreeding and genetic drift. We studied genetic diversity and levels of inbreeding, and their effects on individual fitness and mating choice. We also tested for a genetic signature of the recent demographic bottleneck. For this, we analyzed variation in thousands of nuclear single-nucleotide polymorphisms derived from double digest restriction site-associated DNA sequencing and 1 mitochondrial gene (cytochrome oxidase I). The results revealed that the Bermuda petrel suffered a recent genetic bottleneck and shows low mitochondrial diversity compared with other petrel species. Conversely, nuclear diversity was similar to that of other endangered petrels. Inbreeding levels were not high overall, although some individuals were highly inbred. However, we found no evidence that individual inbreeding or relatedness between mates affected hatching success, or that mate choice is influenced by kinship in this very small population.

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.

Beringia as a high-latitude engine of avian speciation

Beringia is a biogeographically dynamic region that extends from northeastern Asia into northwestern North America. This region has affected avian divergence and speciation in three important ways: (i) by serving as a route for intercontinental colonisation between Asia and the Americas; (ii) by cyclically splitting (and often reuniting) populations, subspecies, and species between these continents; and (iii) by providing isolated refugia through glacial cycles. The effects of these processes can be seen in taxonomic splits of shallow to increasing depths and in the presence of regional endemics. We review the taxa involved in the latter two processes (splitting-reuniting and isolation), with a focus on three research topics: avian diversity, time estimates of the generation of that diversity, and the regions within Beringia that might have been especially important. We find that these processes have generated substantial amounts of avian diversity, including 49 pairs of avian subspecies or species whose breeding distributions largely replace one another across the divide between the Old World and the New World in Beringia, and 103 avian species and subspecies endemic to this region. Among endemics, about one in three is recognised as a full biological species. Endemic taxa in the orders Charadriiformes (shorebirds, alcids, gulls, and terns) and Passeriformes (perching birds) are particularly well represented, although they show very different levels of diversity through evolutionary time. Endemic Beringian Charadriiformes have a 1.31:1 ratio of species to subspecies. In Passeriformes, endemic taxa have a 0.09:1 species-to-subspecies ratio, suggesting that passerine (and thus terrestrial) endemism might be more prone to long-term extinction in this region, although such 'losses' could occur through their being reconnected with wider continental populations during favourable climatic cycles (e.g. subspecies reintegration with other populations). Genetic evidence suggests that most Beringian avian taxa originated over the past 3 million years, confirming the importance of Quaternary processes. There seems to be no obvious clustering in their formation through time, although there might be temporal gaps with lower rates of diversity generation. For at least 62 species, taxonomically undifferentiated populations occupy this region, providing ample potential for future evolutionary diversification.

A high-quality genome assembly highlights the evolutionary history of the great bustard (Otis tarda, Otidiformes)

Conservation genomics often relies on non-invasive methods to obtain DNA fragments which limit the power of multi-omic analyses for threatened species. Here, we report multi-omic analyses based on a well-preserved great bustard individual (Otis tarda, Otidiformes) that was found dead in the mountainous region in Gansu, China. We generate a near-complete genome assembly containing only 18 gaps scattering in 8 out of the 40 assembled chromosomes. We characterize the DNA methylation landscape which is correlated with GC content and gene expression. Our phylogenomic analysis suggests Otidiformes and Musophagiformes are sister groups that diverged from each other 46.3 million years ago. The genetic diversity of great bustard is found the lowest among the four available Otidiformes genomes, possibly due to population declines during past glacial periods. As one of the heaviest migratory birds, great bustard possesses several expanded gene families related to cardiac contraction, actin contraction, calcium ion signaling transduction, as well as positively selected genes enriched for metabolism. Finally, we identify an extremely young evolutionary stratum on the sex chromosome, a rare case among birds. Together, our study provides insights into the conservation genomics, adaption and chromosome evolution of the great bustard.

Species-specific traits mediate avian demographic responses under past climate change

Anticipating species' responses to environmental change is a pressing mission in biodiversity conservation. Despite decades of research investigating how climate change may affect population sizes, historical context is lacking, and the traits that mediate demographic sensitivity to changing climate remain elusive. We use whole-genome sequence data to reconstruct the demographic histories of 263 bird species over the past million years and identify networks of interacting morphological and life history traits associated with changes in effective population size (N_e) in response to climate warming and cooling. Our results identify direct and indirect effects of key traits representing dispersal, reproduction and survival on long-term demographic responses to climate change, thereby highlighting traits most likely to influence population responses to ongoing climate warming.

Small population of the largest water strider after the late Pleistocene and the implications for its conservation

Climate oscillation and its synergistic impacts on habitat fragmentation have been identified as threatening the survival of some extant species. However, the mechanisms by which semi-aquatic insects impacted by such events remain poorly understood. Herein, we studied the largest water strider in the world, Gigantometra gigas, to explore the effect of these two factors on its evolutionary history. The sequences of mitogenomic and nrDNA cluster were utilized to reconstruct phylogenetic relationship among G. gigas populations and its demographic history. Mitochondrial genes were separately reconstructed topologies of that populations and detected remarkable differences. We found that G. gigas populations conform to the isolation-by-distance model, and decline occurred at about 120 ka, which was probably influenced by the climate change during the late Pleistocene and eventually maintained a small effective population size (Ne) around 85,717. The populations in Guangdong Province of China are worthy of note in that they exhibit low genetic diversity, a small Ne around 18,899 individuals, and occupy an area with little suitable future habitat for G. gigas. This work recommends that conservation efforts are implemented to ensure the long-term survival of small G. gigas populations, and notes that further evaluation of their extinction risk under the impacts of human activities is required.

Mitochondrial genome analysis, phylogeny and divergence time evaluation of Strixaluco (Aves, Strigiformes, Strigidae)

Background Prior research has shown that the European peninsulas were the main sources of Strixaluco colonisation of Northern Europe during the late glacial period. However, the phylogenetic relationship and the divergence time between S.aluco from Leigong Mountain Nature Reserve, Guizhou Province, China and the Strigiformes from overseas remains unclear. The mitochondrial genome structure of birds is a covalent double-chain loop structure that is highly conserved and, thus, suitable for phylogenetic analysis. This study examined the phylogenetic relationship and divergence time of Strix using the whole mitochondrial genome of S.aluco. New information In this study, the complete mitochondrial genome of Strixaluco, with a total length of 18,632 bp, is reported for the first time. A total of 37 genes were found, including 22 tRNAs, two rRNAs, 13 protein-coding genes and two non-coding control regions. Certain species of Tytoninae were used as out-group and PhyloSuite software was applied to build the ML-tree and BI-tree of Strigiformes. Finally, the divergence time tree was constructed using BEAST 2.6.7 software and the age of Miosurniadiurna fossil-bearing sediments (6.0-9.5 Ma) was set as internal correction point. The common ancestor of Strix was confirmed to have diverged during the Pleistocene (2.58-0.01 Ma). The combined action of the dramatic uplift of the Qinling Mountains in the Middle Pleistocene and the climate oscillation of the Pleistocene caused Strix divergence between the northern and southern parts of mainland China. The isolation of glacial-interglacial rotation and glacier refuge was the main reason for the divergence of Strixuralensis and S.aluco from their common ancestor during this period. This study provides a reference for the evolutionary history of S.aluco.

Evolution of the germline mutation rate across vertebrates

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself¹. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies². Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis³. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

Changes in the functional diversity of modern bird species over the last million years

Despite evidence of declining biosphere integrity, we currently lack understanding of how the functional diversity associated with changes in abundance among ecological communities has varied over time and before widespread human disturbances. We combine morphological, ecological, and life-history trait data for >260 extant bird species with genomic-based estimates of changing effective population size (N_e) to quantify demographic-based shifts in avian functional diversity over the past million years and under pre-anthropogenic climate warming. We show that functional diversity was relatively stable over this period, but underwent significant changes in some key areas of trait space due to changing species abundances. Our results suggest that patterns of population decline over the Pleistocene have been concentrated in particular regions of trait space associated with extreme reproductive strategies and low dispersal ability, consistent with an overall erosion of functional diversity. Further, species most sensitive to climate warming occupied a relatively narrow region of functional space, indicating that the largest potential population increases and decreases under climate change will occur among species with relatively similar trait sets. Overall, our results identify fluctuations in functional space of extant species over evolutionary timescales and represent the demographic-based vulnerability of different regions of functional space among these taxa. The integration of paleodemographic dynamics with functional trait data enhances our ability to quantify losses of biosphere integrity before anthropogenic disturbances and attribute contemporary biodiversity loss to different drivers over time.

Genomic architecture underlying morphological and physiological adaptation to high elevation in a songbird

Organisms often acquire physiological and morphological modifications to conquer ecological challenges when colonizing new environments which lead to their adaptive evolution. However, deciphering the genomic mechanism of ecological adaptation is difficult because ecological environments are often too complex for straightforward interpretation. Thus, we examined the adaptation of a widespread songbird-the rufous-capped babbler (Cyanoderma ruficeps)-to a relatively simple system: distinct environments across elevational gradients on the mountainous island of Taiwan. We focused on the genomic sequences of 43 birds from five populations to show that the Taiwan group split from its sister group in mainland China around 1-2 million years ago (Ma) and colonized the montane habitats of Taiwan at least twice around 0.03-0.22 Ma. The montane and lowland Taiwan populations diverged with gene flow between them, suggesting strong selection associated with different elevations. We found that the montane babblers had smaller beaks than the lowland ones, consistent with Allen's rule, and identified candidate genes-COL9A1 and SOX11-underlying the beak size changes. We also found that altitudinally divergent mutations were mostly located in noncoding regions and tended to accumulate in chromosomal inversions and autosomes. The altitudinally divergent mutations might regulate genes related to haematopoietic, metabolic, immune, auditory and vision functions, as well as cerebrum morphology and plumage development. The results reveal the genomic bases of morphological and physiological adaptation in this species to the low temperature, hypoxia and high UV light environment at high elevation. These findings improve our understanding of how ecological adaptation drives population divergence from the perspective of genomic architecture.

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.

Genetic evidence for widespread population size expansion in North American boreal birds prior to the Last Glacial Maximum

Pleistocene climate cycles are well documented to have shaped contemporary species distributions and genetic diversity. Northward range expansions in response to deglaciation following the Last Glacial Maximum (LGM; approximately 21 000 years ago) are surmised to have led to population size expansions in terrestrial taxa and changes in seasonal migratory behaviour. Recent findings, however, suggest that some northern temperate populations may have been more stable than expected through the LGM. We modelled the demographic history of 19 co-distributed boreal-breeding North American bird species from full mitochondrial gene sets and species-specific molecular rates. We used these demographic reconstructions to test how species with different migratory strategies were affected by glacial cycles. Our results suggest that effective population sizes increased in response to Pleistocene deglaciation earlier than the LGM, whereas genetic diversity was maintained throughout the LGM despite shifts in geographical range. We conclude that glacial cycles prior to the LGM have most strongly shaped contemporary genetic diversity in these species. We did not find a relationship between historic population dynamics and migratory strategy, contributing to growing evidence that major switches in migratory strategy during the LGM are unnecessary to explain contemporary migratory patterns.

Genomic diversity and differentiation between island and mainland populations of white-tailed eagles (Haliaeetus albicilla)

Divergence in the face of high dispersal capabilities is a documented but poorly understood phenomenon. The white-tailed eagle (Haliaeetus albicilla) has a large geographic dispersal capability and should theoretically be able to maintain genetic homogeneity across its dispersal range. However, following analysis of the genomic variation of white-tailed eagles, from both historical and contemporary samples, clear signatures of ancient biogeographic substructure across Europe and the North-East Atlantic is observed. The greatest genomic differentiation was observed between island (Greenland and Iceland) and mainland (Denmark, Norway and Estonia) populations. The two island populations share a common ancestry from a single mainland population, distinct from the other sampled mainland populations, and despite the potential for high connectivity between Iceland and Greenland they are well separated from each other and are characterized by inbreeding and little variation. Temporal differences also highlight a pattern of regional populations persisting despite the potential for admixture. All sampled populations generally showed a decline in effective population size over time, which may have been shaped by four historical events: (1) Isolation of refugia during the last glacial period 110-115,000 years ago, (2) population divergence following the colonization of the deglaciated areas ~10,000 years ago, (3) human population expansion, which led to the settlement in Iceland ~1100 years ago, and (4) human persecution and exposure to toxic pollutants during the last two centuries.

Assembly, annotation, and comparative analysis of Ipomoea chloroplast genomes provide insights into the parasitic characteristics of Cuscuta species

In the Convolvulaceae family, around 1650 species belonging to 60 genera are widely distributed globally, mainly in the tropical and subtropical regions of America and Asia. Although a series of chloroplast genomes in Convolvulaceae were reported and investigated, the evolutionary and genetic relationships among the chloroplast genomes of the Convolvulaceae family have not been extensively elucidated till now. In this study, we first reported the complete chloroplast genome sequence of Ipomoea pes-caprae, a widely distributed coastal plant with medical values. The chloroplast genome of I. pes-caprae is 161667 bp in length, and the GC content is 37.56%. The chloroplastic DNA molecule of I. pes-caprae is a circular structure composed of LSC (large-single-copy), SSC (small-single-copy), and IR (inverted repeat) regions, with the size of the three regions being 88210 bp, 12117 bp, and 30670 bp, respectively. The chloroplast genome of I. pes-caprae contains 141 genes, and 35 SSRs are identified in the chloroplast genome. Our research results provide important genomic information for the molecular phylogeny of I. pes-caprae. The Phylogenetic analysis of 28 Convolvulaceae chloroplast genomes showed that the relationship of I. pes-caprae with I. involucrata or I. obscura was much closer than that with other Convolvulaccae species. Further comparative analyses between the Ipomoea species and Cuscuta species revealed the mechanism underlying the formation of parasitic characteristics of Cuscuta species from the perspective of the chloroplast genome.

Genetic variation of endangered Jankowski’s Bunting (Emberiza jankowskii): High connectivity and a moderate history of demographic decline

Introduction

Continued discovery of “mismatch” patterns between population size and genetic diversity, involving wild species such as insects, amphibians, birds, mammals, and others, has raised issues about how population history, especially recent dynamics under human disturbance, affects currently standing genetic variation. Previous studies have revealed high genetic diversity in endangered Jankowski’s Bunting. However, it is unclear how the demographic history and recent habitat changes shape the genetic variation of Jankowski’s Bunting.

Methods

To explore the formation and maintenance of high genetic diversity in endangered Jankowski’s Bunting, we used a mitochondrial control region (partial mtDNA CR) and 15 nuclear microsatellite markers to explore the recent demographic history of Jankowski’s Bunting, and we compared the historical and contemporary gene flows between populations to reveal the impact of habitat change on population connectivity. Specifically, we aimed to test the following hypotheses: (1) Jankowski’s Bunting has a large historical Ne and a moderate demographic history; and (2) recent habitat change might have no significant impact on the species’ population connectivity.

Results

The results suggested that large historical effective population size, as well as severe but slow population decline, may partially explain the high observable genetic diversity. Comparison of historical (over the past 4Ne generations) and contemporary (1–3 generations) gene flow indicated that the connectivity between five local populations was only marginally affected by landscape changes.

Discussion

Our results suggest that high population connectivity and a moderate history of demographic decline are powerful explanations for the rich genetic variation in Jankowski’s Bunting. Although there is no evidence that the genetic health of Jankowski’s Bunting is threatened, the time-lag effects on the genetic response to recent environmental changes is a reminder to be cautious about the current genetic characteristics of this species. Where possible, factors influencing genetic variation should be integrated into a systematic framework for conducting robust population health assessments. Given the small contemporary population size, inbreeding, and ecological specialization, we recommend that habitat protection be maintained to maximize the genetic diversity and population connectivity of Jankowski’s Bunting.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Jack of all trades: Genome assembly of Wild Jack and comparative genomics of Artocarpus

Artocarpus (Moraceae), known as breadfruits for their diverse nutritious fruits, is prized for its high-quality timber, medicinal value, and economic importance. Breadfruits are native to Southeast Asia but have been introduced to other continents. The most commonly cultivated species are Artocarpus heterophyllus (Jackfruit) and Artocarpus altilis (Breadfruit). With numerous smaller but nutritionally comparable fruits on a larger tree, Artocarpus hirsutus, also called "Wild Jack" or "Ayani", is an elusive forest species endemic to Indian Western Ghats. In this study, we sequenced and assembled the whole genome of Artocarpus hirsutus sampled from the sacred groves of Coorg, India. To decipher demographic and evolutionary history, we compared our Wild Jack genome with previously published Jackfruit and Breadfruit genomes. Demographic history reconstruction indicates a stronger effect of habitat rather than phylogeny on the population histories of these plants. Repetitive genomic regions, especially LTR Copia, strongly affected the demographic trajectory of A. heterophyllus. Upon further investigation, we found a recent lineage-specific accumulation of LTR Copia in A. heterophyllus, which had a major contribution to its larger genome size. Several genes from starch, sucrose metabolism, and plant hormone signal transduction pathways, in Artocarpus species had signatures of selection and gene family evolution. Our comparative genomic framework provides important insights by incorporating endemic species such as the Wild Jack.

Generation of raptor diversity in Europe: linking speciation with climate changes and the ability to migrate

Europe holds a rich community of diurnal birds of prey, and the highest proportion of transcontinental migratory raptorial species of any landmass. This study will test the hypotheses that the high diversification of the raptor assemblage in Europe is a recent event, that closely related species sharing the same trophic niches can only coexist in sympatry during the breeding period, when food availability is higher, and finally that migration is a function of size, with the smaller species in every trophic group moving further. A consensus molecular phylogeny for the 38 regular breeding species of raptors in Europe was obtained from BirdTree (www.birdtree.org). For the same species, a trophic niche cluster dendrogram was constructed. Size and migratory strategy were introduced in the resulting phylogeny, where trophic groups were also identified. Multispecific trophic groups tended to be composed of reciprocal sister species of different sizes, while monospecific groups (n = 3) were composed of highly specialized species. Many speciation events took place recently, during the glacial cycles of the Quaternary, and size divergence among competing species may be due to character displacement. Nowadays, the smaller species in every trophic group migrate to sub-Saharan Africa. This investigation illustrates how the rich assemblage of diurnal birds of prey in Europe, more diverse and more migratory than, for instance, the North American assemblage at equivalent latitudes, has emerged recently due to the multiplication of look-alike species with similar trophic ecologies, possibly in climate refugia during cold periods.

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.

Genomic novelty within a "great speciator" revealed by a high-quality reference genome of the collared kingfisher (Todiramphus chloris collaris)

Islands are natural laboratories for studying patterns and processes of evolution. Research on island endemic birds has revealed elevated speciation rates and rapid phenotypic evolution in several groups (e.g. white-eyes, Darwin's finches). However, understanding the evolutionary processes behind these patterns requires an understanding of how genotypes map to novel phenotypes. To date, there are few high-quality reference genomes for species found on islands. Here, we sequence the genome of one of Ernst Mayr's "great speciators," the collared kingfisher (Todiramphus chloris collaris). Utilizing high molecular weight DNA and linked-read sequencing technology, we assembled a draft high-quality genome with highly contiguous scaffolds (scaffold N50 = 19 Mb). Based on universal single-copy orthologs, we estimated a gene space completeness of 96.6% for the draft genome assembly. The population demographic history analyses reveal a distinct pattern of contraction and expansion in population size throughout the Pleistocene. Comparative genomic analysis of gene family evolution revealed that species-specific and rapidly expanding gene families in the collared kingfisher (relative to other Coraciiformes) are mainly involved in the ErbB signaling pathway and focal adhesion. Todiramphus kingfishers are a species-rich group that has become a focus of speciation research. This draft genome will be a platform for future taxonomic, phylogeographic, and speciation research in the group. For example, target genes will enable testing of changes in sensory structures associated with changes in vision and taste genes across kingfishers.

Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum

Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effects of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet, few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.

Genomic basis of insularity and ecological divergence in barn owls (Tyto alba) of the Canary Islands

Islands, and the particular organisms that populate them, have long fascinated biologists. Due to their isolation, islands offer unique opportunities to study the effect of neutral and adaptive mechanisms in determining genomic and phenotypical divergence. In the Canary Islands, an archipelago rich in endemics, the barn owl (Tyto alba), present in all the islands, is thought to have diverged into a subspecies (T. a. gracilirostris) on the eastern ones, Fuerteventura and Lanzarote. Taking advantage of 40 whole-genomes and modern population genomics tools, we provide the first look at the origin and genetic makeup of barn owls of this archipelago. We show that the Canaries hold diverse, long-standing and monophyletic populations with a neat distinction of gene pools from the different islands. Using a new method, less sensitive to structure than classical F_ST, to detect regions involved in local adaptation to insular environments, we identified a haplotype-like region likely under selection in all Canaries individuals and genes in this region suggest morphological adaptations to insularity. In the eastern islands, where the subspecies is present, genomic traces of selection pinpoint signs of adapted body proportions and blood pressure, consistent with the smaller size of this population living in a hot arid climate. In turn, genomic regions under selection in the western barn owls from Tenerife showed an enrichment in genes linked to hypoxia, a potential response to inhabiting a small island with a marked altitudinal gradient. Our results illustrate the interplay of neutral and adaptive forces in shaping divergence and early onset speciation.

Potential millennial-scale avian declines by humans in southern China

Mounting observational records demonstrate human-caused faunal decline in recent decades, while accumulating archaeological evidence suggests an early biodiversity impact of human activities during the Holocene. A fundamental question arises concerning whether modern wildlife population declines began during early human disturbance. Here, we performed a population genomic analysis of six common forest birds in East Asia to address this question. For five of them, demographic history inference based on 25-33 genomes of each species revealed dramatic population declines by 4- to 48-fold over millennia (e.g. 2000-5000 thousand years ago). Nevertheless, summary statistics detected nonsignificant correlations between these population size trajectories and Holocene temperature variations, and ecological niche models explicitly predicted extensive range persistence during the Holocene, implying limited demographic consequence of Holocene climate change. Further analyses suggest high negative correlations between the reconstructed population declines and human disturbance intensities and indicate a potential driver of human activities. These findings provide a deep-time and large-scale insight into the recently recognized avifaunal decline and support an early origin hypothesis of human effects on biodiversity. Overall, our study sheds light on the current biodiversity crisis in the context of long-term human-environment interactions and offers a multi-evidential framework for quantitatively assessing the ecological consequences of human disturbance.

Rainfall and sea level drove the expansion of seasonally flooded habitats and associated bird populations across Amazonia

Spatial arrangement of distinct Amazonian environments through time and its effect on specialized biota remain poorly known, fueling long-lasting debates about drivers of biotic diversification. We address the late Quaternary sediment deposition that assembled the world's largest seasonally flooded ecosystems. Genome sequencing was used to reconstruct the demographic history of bird species specialized in either early successional vegetation or mature floodplain forests. Sediment deposition that built seasonally flooded habitats accelerated throughout the Holocene (last 11,700 years) under sea level highstand and intensification of the South American Monsoon, at the same time as global increases in atmospheric methane concentration. Bird populations adapted to seasonally flooded habitats expanded due to enlargement of Amazonian river floodplains and archipelagos. Our findings suggest that the diversification of the biota specialized in seasonally flooded habitats is coupled to sedimentary budget changes of large rivers, which rely on combined effects of sea level and rainfall variations.

This study found that millennial periods of higher rainfall combined with rising sea level enhanced sediment accumulation in Amazonian rivers valleys. This fuelled synchronous expansion of vegetation adapted to seasonally flooded substrates and its specialized bird populations, showing how global climate changes can affect specific Amazonian species.

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.

14,000 year old DNA reveals the evolutionary dynamics and adaptations of South American vultures.

De Novo Whole-Genome Sequencing and Assembly of the Yellow-Throated Bunting (Emberiza elegans) Provides Insights into Its Evolutionary Adaptation

Simple Summary

We report the genomic sequence of Emberiza elegans for understanding the evolutionary mechanisms of environmental adaptation and for studying a more effective genetic monitoring of this species. The E. elegans assembly was approximately 1.14 Gb, with a scaffold N50 of 28.94 Mb. About 15,868 protein-coding genes were predicted, and 16.62% of the genome was identified as having repetitive elements. Our genomic evolution analyses found considerable numbers of adaptive genes that may help the yellow-throated bunting cope with migratory behavior and environmental stressors of diseases. These results provide us with new insights into genomic evolution and adaptation, thus providing a valuable resource for further studies of population genetic diversity and genome evolution in this species.

Abstract

Yellow-throated bunting is a small migratory songbird unique to the Palearctic region. However, the genetic studies of this species remain limited, with no nuclear genomic sequence reported to date. In this study, the genomic DNA from the bird was sequenced in long reads using Nanopore sequencing technology. Combining short-read sequencing, the genome was well-assembled and annotated. The final length of the assembly is approximately 1.14 Gb, with a scaffold N50 of 28.94 Mb. About 15,868 protein-coding genes were predicted, and 16.62% of the genome was identified as having repetitive elements. Comparative genomic analysis showed numerous expanded gene families and positively selected genes significantly enriched in those KEGG pathways that are associated with migratory behavior adaptation and immune response. Here, this newly generated de novo genome of the yellow-throated bunting using long reads provide the research community with a valuable resource for further studies of population genetic diversity and genome evolution in this species.

Impact of host demography and evolutionary history on endosymbiont molecular evolution: A test in carpenter ants (genus Camponotus) and their Blochmannia endosymbionts

Obligate endosymbioses are tight associations between symbionts and the hosts they live inside. Hosts and their associated obligate endosymbionts generally exhibit codiversification, which has been documented in taxonomically diverse insect lineages. Host demography (e.g., effective population sizes) may impact the demography of endosymbionts, which may lead to an association between host demography and the patterns and processes of endosymbiont molecular evolution. Here, we used whole-genome sequencing data for carpenter ants (Genus Camponotus; subgenera Camponotus and Tanaemyrmex) and their Blochmannia endosymbionts as our study system to address whether Camponotus demography shapes Blochmannia molecular evolution. Using whole-genome phylogenomics, we confirmed previous work identifying codiversification between carpenter ants and their Blochmannia endosymbionts. We found that Blochmannia genes have evolved at a pace ~30× faster than that of their hosts' molecular evolution and that these rates are positively associated with host rates of molecular evolution. Using multiple tests for selection in Blochmannia genes, we found signatures of positive selection and shifts in selection strength across the phylogeny. Host demography was associated with Blochmannia shifts toward increased selection strengths, but not associated with Blochmannia selection relaxation, positive selection, genetic drift rates, or genome size evolution. Mixed support for relationships between host effective population sizes and Blochmannia molecular evolution suggests weak or uncoupled relationships between host demography and Blochmannia population genomic processes. Finally, we found that Blochmannia genome size evolution was associated with genome-wide estimates of genetic drift and number of genes with relaxed selection pressures.

Birds adapted to cold conditions show greater changes in range size related to past climatic oscillations than temperate birds

Investigation of ecological responses of species to past climate oscillations provides crucial information to understand the effects of global warming. In this work, we investigated how past climate changes affected the distribution of six bird species with different climatic requirements and migratory behaviours in the Western Palearctic and in Africa. Species Distribution Models and Marine Isotopic Stage (MIS) 2 fossil occurrences of selected species were employed to evaluate the relation between changes in range size and species climatic tolerances. The Last Glacial Maximum (LGM) range predictions, generally well supported by the MIS 2 fossil occurrences, suggest that cold-dwelling species considerably expanded their distribution in the LGM, experiencing more pronounced net changes in range size compared to temperate species. Overall, the thermal niche proves to be a key ecological trait for explaining the impact of climate change in species distributions. Thermal niche is linked to range size variations due to climatic oscillations, with cold-adapted species currently suffering a more striking range reduction compared to temperate species. This work also supports the persistence of Afro-Palearctic migrations during the LGM due to the presence of climatically suitable wintering areas in Africa even during glacial maxima.

The First Genome of the Balearic Shearwater (Puffinus mauretanicus) Provides a Valuable Resource for Conservation Genomics and Sheds Light on Adaptation to a Pelagic lifestyle

The Balearic shearwater (Puffinus mauretanicus) is the most threatened seabird in Europe and a member of the most speciose group of pelagic seabirds, the order Procellariiformes, which exhibit extreme adaptations to a pelagic lifestyle. The fossil record suggests that human colonisation of the Balearic Islands resulted in a sharp decrease of the Balearic shearwater population size. Currently, populations of the species continue to be decimated mainly due to predation by introduced mammals and bycatch in longline fisheries, with some studies predicting its extinction by 2070. Here, using a combination of short and long reads, we generate the first high-quality reference genome for the Balearic shearwater, with a completeness amongst the highest across available avian species. We used this reference genome to study critical aspects relevant to the conservation status of the species and to gain insights into the adaptation to a pelagic lifestyle of the order Procellariiformes. We detected relatively high levels of genome-wide heterozygosity in the Balearic shearwater despite its reduced population size. However, the reconstruction of its historical demography uncovered an abrupt population decline potentially linked to a reduction of the neritic zone during the Penultimate Glacial Period (∼194-135 ka). Comparative genomics analyses uncover a set of candidate genes that may have played an important role into the adaptation to a pelagic lifestyle of Procellariiformes, including those for the enhancement of fishing capabilities, night vision, and the development of natriuresis. The reference genome obtained will be the crucial in the future development of genetic tools in conservation efforts for this Critically Endangered species.