Island songbirds as windows into evolution in small populations

Genomes of critically endangered saola are shaped by population structure and purging

The saola is one of the most elusive large mammals, standing at the brink of extinction. We constructed a reference genome and resequenced 26 saola individuals, confirming the saola as a basal member of the Bovini. Despite its small geographic range, we found that the saola is partitioned into two populations with high genetic differentiation (FST = 0.49). We estimate that these populations diverged and started declining 5,000-20,000 years ago, possibly due to climate changes and exacerbated by increasing human activities. The saola has long tracts without genomic diversity; however, most of these tracts are not shared by the two populations. Saolas carry a high genetic load, yet their gradual decline resulted in the purging of the most deleterious genetic variation. Finally, we find that combining the two populations, e.g., in an eventual captive breeding program, would mitigate the genetic load and increase the odds of species survival.

The Relative Roles of Selection and Drift in the Chaffinch Radiation (Aves: Fringilla) Across the Atlantic Archipelagos of Macaronesia

Island populations diverge from the mainland and from each other by both natural selection and neutral forces such as founder effects and genetic drift. In this work, we aim to determine the relative roles of selection and drift in the diversification of chaffinches (Fringilla spp.) in Macaronesia. We tested the hypothesis that taxa inhabiting Macaronesian archipelagos, which exhibit significant differences in habitat and climate compared with the mainland, should experience distinct ecological pressures, leading to divergence at loci under selection related to environmental variables. To determine the role of local adaptation in the differentiation of these taxa, we performed genotype-environment association (GEA) analyses using ten environmental variables and 52,306 single nucleotide polymorphism markers obtained from genotyping-by-sequencing (GBS) in 79 chaffinches. Redundancy analysis (RDA) revealed that genomic variation is associated with environmental variables and identified candidate genes related to phenotypic traits and environmental variables. Variables associated with habitat type and precipitation, together with drift, played an important role in the genomic differentiation between chaffinches from Macaronesia and the mainland, as well as within the Canarian archipelago. Genetic drift was identified as the main factor in the differentiation of North African chaffinches from Macaronesia and mainland Europe, as well as Madeira chaffinches from those in the Canary Islands. Finally, chaffinches from the Canary Islands show an incipient diversification process at the genetic and phenotypic level driven by both selection and neutral processes. Our results suggest that both habitat-driven local adaptation and drift have played a role in the radiation of chaffinch taxa in Macaronesia.

Genomes of Galápagos Mockingbirds Reveal the Impact of Island Size and Past Demography on Inbreeding and Genetic Load in Contemporary Populations

Restricted range size brings about noteworthy genetic consequences that may affect the viability of a population and eventually its extinction. Particularly, the question if an increase in inbreeding can avert the accumulation of genetic load via purging is hotly debated in the conservation genetic field. Insular populations with limited range sizes represent an ideal setup for relating range size to these genetic factors. Leveraging a set of eight differently sized populations of Galápagos mockingbirds (Mimus), we investigated how island size shaped effective population size (Ne), inbreeding and genetic load. We assembled a genome of M. melanotis and genotyped three individuals per population by whole-genome resequencing. Demographic inference showed that the Ne of most populations remained high after the colonisation of the archipelago 1-2 Mya. Ne decline in M. parvulus happened only 10-20 Kya, whereas the critically endangered M. trifasciatus showed a longer history of reduced Ne. Despite these historical fluctuations, the current island size determines Ne in a linear fashion. In contrast, significant inbreeding coefficients, derived from runs of homozygosity, were identified only in the four smallest populations. The index of additive genetic load suggested purging in M. parvulus, where the smallest populations showed the lowest load. By contrast, M. trifasciatus carried the highest genetic load, possibly due to a recent rapid bottleneck. Overall, our study demonstrates a complex effect of demography on inbreeding and genetic load, providing implications in conservation genetics in general and in a conservation project of M. trifasciatus in particular.

Genomic insights into the biogeography and evolution of Galápagos iguanas

Galápagos iguanas are a monophyletic group endemic to the Galápagos archipelago, comprising the marine iguana Amblyrhynchus cristatus and three species of land iguanas: Conolophus subcristatus, C. pallidus and C. marthae. The biogeographic history of the land species in relation to their current distributions remains uncertain, in particular the origins of C. marthae, which is restricted to a small area of the northern part of Isabela Island. The classification of C. pallidus as a separate species has also been debated. We analyzed DNA sequences (RADseq) to reconstruct demographic histories of selected local populations of all Galápagos iguana species and estimate their divergence times within a multispecies coalescent framework. Our results indicate an early date for the colonization of Galápagos by iguanas, relative to island formation, at ca. 10 Mya, and support a recent split of C. marthae via allopatric speciation, after the emergence of Isabela Island, at ca. 0.57 Mya. We find contrasting demographic histories in C. marthae and the syntopic population of C. subcristatus, suggesting competitive interaction between these species. We also confirm that the divergence of C. pallidus from C. subcristatus is recent (0.09 Mya) and close in time to the split between populations of C. subcristatus from different islands. Our genetic data support recent census estimates indicating a relatively small current effective population size (Ne) in all the studied populations. Our findings shed light on the evolutionary history of Galápagos iguanas and emphasize the need for targeted conservation strategies.

Island size shapes genomic diversity in a great speciator (Aves: Zosterops)

Islands have long represented natural laboratories for studying many aspects of ecology and evolutionary biology, from speciation to community assembly. One aspect that has been well documented is the correlation between island size and taxonomic diversity, likely due to decreased complexity and population size on small islands. This same logic can apply to genetic diversity, which should predictably decrease with effective population size. The island size-diversity correlation has received support over the years but often focuses on single metrics of genetic diversity. Here, we use Zosterops white-eyes in the Solomon Islands to study the correlation between island size and various metrics related to genetic diversity, including runs of homozygosity and fixation of transposable elements. We find that almost all these metrics strongly correlate with island size, and in turn with each other. We infer that island size is independently correlated with these different variables, demonstrating that population size impacts genomic metrics of diversity in a variety of ways across temporal and hierarchical scales.

Contrasting Genomic Diversity and Inbreeding Levels Among Two Closely Related Falcon Species With Overlapping Geographic Distributions

Genomic resources are valuable to examine historical demographic patterns and their effects to better inform management and conservation of threatened species. We evaluated population trends and genome-wide variation in the near-threatened Orange-breasted Falcon (Falco deiroleucus) and its more common sister species, the Bat Falcon (F. rufigularis), to explore how the two species differ in genomic diversity as influenced by their contrasting long-term demographic histories. We generated and aligned whole genome resequencing data for 12 Orange-breasted Falcons and 9 Bat Falcons to an annotated Gyrfalcon (F. rusticolus) reference genome that retained approximately 22.4 million biallelic autosomal SNPs (chromosomes 1-22). Our analyses indicated much lower genomic diversity in Orange-breasted Falcons compared to Bat Falcons. All sampled Orange-breasted Falcons were significantly more inbred than the sampled Bat Falcons, with values similar to those observed in island-mainland species comparisons. The distribution of runs of homozygosity showed variation suggesting long-term low population size and the possibility of bottlenecks in Orange-breasted Falcons contrasting with consistently larger populations in Bat Falcons. Analysis of genetic load suggests that Orange-breasted Falcons are less likely to experience inbreeding depression than Bat Falcons due to reduced inbreeding load but are at elevated risk from fixation of deleterious gene variants and perhaps a reduced adaptive potential. These genomic analyses highlight differences in the historical demography of two closely related species that have influenced their current genomic diversity and should result in differing strategies for their continued conservation.

Inferring Long-Term and Short-Term Determinants of Genetic Diversity in Honey Bees: Beekeeping Impact and Conservation Strategies

Bees are vital pollinators in natural and agricultural landscapes around the globe, playing a key role in maintaining flowering plant biodiversity and ensuring food security. Among the honey bee species, the Western honey bee (Apis mellifera) is particularly significant, not only for its extensive crop pollination services but also for producing economically valuable products such as honey. Here, we analyzed whole-genome sequence data from four Apis species to explore how honey bee evolution has shaped current diversity patterns. Using Approximate Bayesian Computation, we first reconstructed the demographic history of A. mellifera in Europe, finding support for postglacial secondary contacts, therefore predating human-mediated transfers linked to modern beekeeping. However, our analysis of recent demographic changes reveals significant bottlenecks due to beekeeping practices, which have notably affected genetic diversity. Black honey bee populations from conservatories, particularly those on islands, exhibit considerable genetic loss, highlighting the need to evaluate the long-term effectiveness of current conservation strategies. Additionally, we observed a high degree of conservation in the genomic landscapes of nucleotide diversity across the four species, despite a divergence gradient spanning over 15 million years, consistent with a long-term conservation of the recombination landscapes. Taken together, our results provide the most comprehensive assessment of diversity patterns in honey bees to date and offer insights into the optimal management of resources to ensure the long-term persistence of honey bees and their invaluable pollination services.

Avian Island Radiations Shed Light on the Dynamics of Adaptive and Nonadaptive Radiation

Understanding the mechanisms underlying species formation and differentiation is a central goal of evolutionary biology and a formidable challenge. This understanding can provide valuable insights into the origins of the astonishing diversity of organisms living on our planet. Avian evolutionary radiations on islands have long fascinated biologists as they provide the ideal variation to study the ecological and evolutionary forces operating on the continuum between incipient lineages to complete speciation. In this review, we summarize the key insights gained from decades of research on adaptive and nonadaptive radiations of both extant and extinct insular bird species. We present a new comprehensive global list of potential avian radiations on oceanic islands, based on published island species checklists, taxonomic studies, and phylogenetic analyses. We demonstrate that our understanding of evolutionary processes is being greatly enhanced through the use of genomic tools. However, to advance the field, it is critical to complement this information with a solid understanding of the ecological and behavioral traits of both extinct and extant avian island species.

A comprehensive review of livestock development: insights into domestication, phylogenetics, diversity, and genomic advances

Livestock plays an essential role in sustaining human livelihoods, offering a diverse range of species integral to food security, economic stability, and cultural traditions. The domestication of livestock, which began over 10,000 years ago, has driven significant genetic changes in species such as cattle, buffaloes, sheep, goats, and pigs. Recent advancements in genomic technologies, including next-generation sequencing (NGS), genome-wide association studies (GWAS), and genomic selection, have dramatically enhanced our understanding of these genetic developments. This review brings together key research on the domestication process, phylogenetics, genetic diversity, and selection signatures within major livestock species. It emphasizes the importance of admixture studies and evolutionary forces like natural selection, genetic drift, and gene flow in shaping livestock populations. Additionally, the integration of machine learning with genomic data offers new perspectives on the functional roles of genes in adaptation and evolution. By exploring these genomic advancements, this review provides insights into genetic variation and evolutionary processes that could inform future approaches to improving livestock management and adaptation to environmental challenges, including climate change.

Evolutionary history of an island endemic, the Azorean common quail

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.

A high-quality genome assembly of the Spectacled Fulvetta (Fulvetta ruficapilla) endemic to China

The Spectacled Fulvetta (Fulvetta ruficapilla) is the type species of Fulvetta, an evolutionarily distinct group whose species show a high degree of sympatry in distribution and phenotypic convergence. To pave the way for insights into their adaptive evolution and speciation, we have assembled the first high quality reference genome for F. ruficapilla using high-fidelity (HiFi) long-read and Hi-C sequencing technologies. The resulting assembly spans a total of ~1.21 Gb with a contig N50 of 18.8 Mb and scaffold N50 of 75.9 Mb, and has a BUSCO completeness of 97.0%. The quality assessment suggests a high standard in base accuracy, continuity, and completeness of the assembly, comparable or close to that of Vertebrate Genomes Project. On this basis, we have annotated 23,774 protein-coding genes, of which 18,832 are functionally identified. The availability of this high-quality genome provides a solid foundation for the future studies of evolution and local adaptation in birds.

Fitness consequences of structural variation inferred from a House Finch pangenome

Genomic structural variants (SVs) play a crucial role in adaptive evolution, yet their average fitness effects and characterization with pangenome tools are understudied in wild animal populations. We constructed a pangenome for House Finches (Haemorhous mexicanus), a model for studies of host-pathogen coevolution, using long-read sequence data on 16 individuals (32 de novo-assembled haplotypes) and one outgroup. We identified 887,118 SVs larger than 50 base pairs, mostly (60%) involving repetitive elements, with reduced SV diversity in the eastern US as a result of its introduction by humans. The distribution of fitness effects of genome-wide SVs was estimated using maximum likelihood approaches and revealed that SVs in both coding and noncoding regions were on average more deleterious than smaller indels or single nucleotide polymorphisms. The reference-free pangenome facilitated identification of a > 10-My-old, 11-megabase-long pericentric inversion on chromosome 1. We found that the genotype frequencies of the inversion, estimated from 135 birds widely sampled temporally and geographically, increased steadily over the 25 y since House Finches were first exposed to the bacterial pathogen Mycoplasma gallisepticum and showed signatures of balancing selection, capturing genes related to immunity and telomerase activity. We also observed shorter telomeres in populations with a greater number of years exposure to Mycoplasma. Our study illustrates the utility of long-read sequencing and pangenome methods for understanding wild animal populations, estimating fitness effects of genome-wide SVs, and advancing our understanding of adaptive evolution through structural variation.

Repeated evolution on oceanic islands: comparative genomics reveals species-specific processes in birds

Understanding the interplay between genetic drift, natural selection, gene flow, and demographic history in driving phenotypic and genomic differentiation of insular populations can help us gain insight into the speciation process. Comparing patterns across different insular taxa subjected to similar selective pressures upon colonizing oceanic islands provides the opportunity to study repeated evolution and identify shared patterns in their genomic landscapes of differentiation. We selected four species of passerine birds (Common Chaffinch Fringilla coelebs/canariensis, Red-billed Chough Pyrrhocorax pyrrhocorax, House Finch  Haemorhous mexicanus and Dark-eyed/island Junco Junco hyemalis/insularis) that have both mainland and insular populations. Changes in body size between island and mainland populations were consistent with the island rule. For each species, we sequenced whole genomes from mainland and insular individuals to infer their demographic history, characterize their genomic differentiation, and identify the factors shaping them. We estimated the relative (Fst) and absolute (dxy) differentiation, nucleotide diversity (π), Tajima's D, gene density and recombination rate. We also searched for selective sweeps and chromosomal inversions along the genome. All species shared a marked reduction in effective population size (Ne) upon island colonization. We found diverse patterns of differentiated genomic regions relative to the genome average in all four species, suggesting the role of selection in island-mainland differentiation, yet the lack of congruence in the location of these regions indicates that each species evolved differently in insular environments. Our results suggest that the genomic mechanisms involved in the divergence upon island colonization-such as chromosomal inversions, and historical factors like recurrent selection-differ in each species, despite the highly conserved structure of avian genomes and the similar selective factors involved. These differences are likely influenced by factors such as genetic drift, the polygenic nature of fitness traits and the action of case-specific selective pressures.

Purifying Selection Influences the Comparison of Heterozygosities between Populations

Heterozygosity is a fundamental measure routinely used to compare between populations to infer the level of genetic variation and their relative effective population sizes. However, such comparison is highly influenced by the magnitude of selection pressure on the genomic regions used. Using over 2 million Single Nucleotide Variants (SNVs) from chimpanzee and mouse populations, this study shows that the heterozygosities estimated using neutrally evolving sites of large populations were two times higher than those of small populations. However, this difference was only ~1.6 times for the heterozygosities estimated using nonsynonymous sites. This suggests an excess in the nonsynonymous heterozygosities due to the segregation of deleterious variants in small populations. This excess in the nonsynonymous heterozygosities of the small populations was estimated to be 23-31%. Further analysis revealed that the magnitude of the excess is modulated by effective population size (Ne) and selection intensity (s). Using chimpanzee populations, this investigation found that the excess in nonsynonymous diversity in the small population was little (6%) when the difference between the Ne values of large and small populations was small (2.4 times). Conversely, this was high (23%) when the difference in Ne was large (5.9 times). Analysis using mouse populations showed that the excess in the nonsynonymous diversity of highly constrained genes of the small population was much higher (38%) than that observed for the genes under relaxed selective constraints (21%). Similar results were observed when the expression levels of genes were used as a proxy for selection intensity. These results emphasize the use of neutral regions, less constrained genes, or lowly expressed genes when comparing the heterozygosities between populations.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

The effect of the demographic history on the evolution of senescence: A potential new test of the mutation accumulation theory

The different evolutionary theories of senescence predict different directions for the correlation between the population size and the intensity of senescence. Using simulations, I highlighted how the effect of the population size on the intensity of senescence could be reinforced by the time since populations have been large or small. I devised a mutation-selection model in which the effect of the mutations was age-specific. Several small populations diverged from a same large population at different points in time. At the end of the simulation, the correlation between the time since the populations had been small and the rate of senescence was positive under the mutation accumulation theory and negative under the antagonistic pleiotropy theory. The phenomenon was strong enough to reverse the usually negative relationship between the intensity of senescence and the generation time. These mutually-exclusive predictions could help broaden the taxonomic support for the mutation accumulation theory of senescence, currently mostly supported in humans and lab invertebrates. I briefly mention a few potential applications in real-life systems.

The burden of anthropogenic changes and mutation load in a critically endangered harrier from the Reunion biodiversity hotspot, Circus maillardi

Anthropogenic impact is causing the decline of a large proportion of species worldwide and reduces their genetic diversity. Island species typically have smaller ranges than continental species. As a consequence, island species are particularly liable to undergo population bottlenecks, giving rise to conservation challenges such as inbreeding and unmasking of deleterious genetic load. Such challenges call for more detailed assessments of the genetic make-up of threatened island populations. The Mascarene islands (Indian Ocean) present many prime examples, being unusual in having been pristine until first human arrival ~400 years ago, following which anthropogenic pressure was unusually intense. A threatened harrier (Circus maillardi) endemic to the westernmost island of the archipelago is a good example of the challenges faced by species that have declined to small population size following intense anthropogenic pressure. In this study, we use an extensive set of population genomic tools to quantify variation at near-neutral and coding loci, in order to test the historical impact of human activity on this species, and evaluate the species' (mal)adaptive potential. We observed low but significant genetic differentiation between populations on the West and North-East sides of the island, echoing observations in other endemic species. Inbreeding was significant, with a substantial fraction of samples being first or second-degree relatives. Historical effective population sizes have declined from ~3000 to 300 individuals in the past 1000 years, with a more recent drop ~100 years ago consistent with human activity. Based on our simulations and comparisons with a close relative (Circus melanoleucos), this demographic history may have allowed purging of the most deleterious variants but is unlikely to have allowed the purging of mildly deleterious variants. Our study shows how using relatively affordable methods can reveal the massive impact that human activity may have on the genetic diversity and adaptive potential of island populations, and calls for urgent action to closely monitor the reproductive success of such endemic populations, in association with genetic studies.

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.

Evidence that genetic drift not adaptation drives fast-Z and large-Z effects in Ficedula flycatchers

The sex chromosomes have been hypothesized to play a key role in driving adaptation and speciation across many taxa. The reason for this is thought to be the hemizygosity of the heteromorphic part of sex chromosomes in the heterogametic sex, which exposes recessive mutations to natural and sexual selection. The exposure of recessive beneficial mutations increases their rate of fixation on the sex chromosomes, which results in a faster rate of evolution. In addition, genetic incompatibilities between sex-linked loci are exposed faster in the genomic background of hybrids of divergent lineages, which makes sex chromosomes contribute disproportionately to reproductive isolation. However, in birds, which show a Z/W sex determination system, the role of adaptation versus genetic drift as the driving force of the faster differentiation of the Z chromosome (fast-Z effect) and the disproportionate role of the Z chromosome in reproductive isolation (large-Z effect) are still debated. Here, we address this debate in the bird genus Ficedula flycatchers based on population-level whole-genome sequencing data of six species. Our analysis provides evidence for both faster lineage sorting and reduced gene flow on the Z chromosome than the autosomes. However, these patterns appear to be driven primarily by the increased role of genetic drift on the Z chromosome, rather than an increased rate of adaptive evolution. Genomic scans of selective sweeps and fixed differences in fact suggest a reduced action of positive selection on the Z chromosome.

Genomic signatures of convergent shifts to plunge-diving behavior in birds

Understanding the genetic basis of convergence at broad phylogenetic scales remains a key challenge in biology. Kingfishers (Aves: Alcedinidae) are a cosmopolitan avian radiation with diverse colors, diets, and feeding behaviors-including the archetypal plunge-dive into water. Given the sensory and locomotor challenges associated with air-water transitions, kingfishers offer a powerful opportunity to explore the effects of convergent behaviors on the evolution of genomes and phenotypes, as well as direct comparisons between continental and island lineages. Here, we use whole-genome sequencing of 30 diverse kingfisher species to identify the genomic signatures associated with convergent feeding behaviors. We show that species with smaller ranges (i.e., on islands) have experienced stronger demographic fluctuations than those on continents, and that these differences have influenced the dynamics of molecular evolution. Comparative genomic analyses reveal positive selection and genomic convergence in brain and dietary genes in plunge-divers. These findings enhance our understanding of the connections between genotype and phenotype in a diverse avian radiation.

Impacts of Population Size and Domestication Process on Genetic Diversity and Genetic Load in Genus Ovis

In theoretical biology, a prevailing hypothesis posits a profound interconnection between effective population size (Ne), genetic diversity, inbreeding, and genetic load. The domestication and improvement processes are believed to be pivotal in diminishing genetic diversity while elevating levels of inbreeding and increasing genetic load. In this study, we performed a whole genome analysis to quantity genetic diversity, inbreeding, and genetic load across seven wild Ovis species and five domesticated sheep breeds. Our research demonstrates that the genetic load and diversity of species in the genus Ovis have no discernible impact on recent Ne, and three species within the subgenus Pachyceros tend to carry a higher genetic load and lower genetic diversity patterns. The results coincide with these species' dramatic decline in population sizes within the subgenus Pachyceros ~80-250 thousand years ago. European mouflon presented with the lowest Ne, lower genetic diversity, and higher individual inbreeding coefficient but a lower genetic load (missense and LoF). This suggests that the small Ne of European mouflon could reduce harmful mutations compared to other species within the genus Ovis. We showed lower genetic diversity in domesticated sheep than in Asiatic mouflon, but counterintuitive patterns of genetic load, i.e., lower weak genetic load (missense mutation) and no significant difference in strong genetic load (LoF mutation) between domestic sheep and Asiatic mouflon. These findings reveal that the "cost of domestication" during domestication and improvement processes reduced genetic diversity and purified weak genetic load more efficiently than wild species.

Adaptation to the High-Arctic island environment despite long-term reduced genetic variation in Svalbard reindeer

Typically much smaller in number than their mainland counterparts, island populations are ideal systems to investigate genetic threats to small populations. The Svalbard reindeer (Rangifer tarandus platyrhynchus) is an endemic subspecies that colonized the Svalbard archipelago ca. 6,000-8,000 years ago and now shows numerous physiological and morphological adaptations to its arctic habitat. Here, we report a de-novo chromosome-level assembly for Svalbard reindeer and analyze 133 reindeer genomes spanning Svalbard and most of the species' Holarctic range, to examine the genomic consequences of long-term isolation and small population size in this insular subspecies. Empirical data, demographic reconstructions, and forward simulations show that long-term isolation and high inbreeding levels may have facilitated the reduction of highly deleterious-and to a lesser extent, moderately deleterious-variation. Our study indicates that long-term reduced genetic diversity did not preclude local adaptation to the High Arctic, suggesting that even severely bottlenecked populations can retain evolutionary potential.

Adaptive phenotypic and genomic divergence in the common chaffinch (Fringilla coelebs) following niche expansion within a small oceanic island

According to models of ecological speciation, adaptation to adjacent, contrasting habitat types can lead to population divergence given strong enough environment-driven selection to counteract the homogenizing effect of gene flow. We tested this hypothesis in the common chaffinch (Fringilla coelebs) on the small island of La Palma, Canary Islands, where it occupies two markedly different habitats. Isotopic (δ13 C, δ15 N) analysis of feathers indicated that birds in the two habitats differed in ecosystem and/or diet, and analysis of phenotypic traits revealed significant differences in morphology and plumage colouration that are consistent with ecomorphological and ecogeographical predictions respectively. A genome-wide survey of single-nucleotide polymorphism revealed marked neutral structure that was consistent with geography and isolation by distance, suggesting low dispersal. In contrast, loci putatively under selection identified through genome-wide association and genotype-environment association analyses, revealed amarked adaptive divergence between birds in both habitats. Loci associated with phenotypic and environmental differences among habitats were distributed across the genome, as expected for polygenic traits involved in local adaptation. Our results suggest a strong role for habitat-driven local adaptation in population divergence in the chaffinches of La Palma, a process that appears to be facilitated by a strong reduction in effective dispersal distances despite the birds' high dispersal capacity.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F _ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N _e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Using conservation genetics to prioritise management options for an endangered songbird

Genetic data can be highly informative for answering questions relevant to practical conservation efforts, but remain one of the most neglected aspects of species recovery plans. Framing genetic questions with reference to practical and tractable conservation objectives can help bypass this limitation of the application of genetics in conservation. Using a single-nucleotide polymorphism dataset from reduced-representation sequencing (DArTSeq), we conducted a genetic assessment of remnant populations of the endangered forty-spotted pardalote (Pardalotus quadragintus), a songbird endemic to Tasmania, Australia. Our objectives were to inform strategies for the conservation of genetic diversity in the species and estimate effective population sizes and patterns of inter-population movement to identify management units relevant to population conservation and habitat restoration. We show population genetic structure and identify two small populations on mainland Tasmania as 'satellites' of larger Bruny Island populations connected by migration. Our data identify management units for conservation objectives relating to genetic diversity and habitat restoration. Although our results do not indicate the immediate need to genetically manage populations, the small effective population sizes we estimated for some populations indicate that they are vulnerable to genetic drift, highlighting the urgent need to implement habitat restoration to increase population size and to conduct genetic monitoring. We discuss how our genetic assessment can be used to inform management interventions for the forty-spotted pardalote and show that by assessing contemporary genetic aspects, valuable information for conservation planning and decision-making can be produced to guide actions that account for genetic diversity and increase chances of recovery in species of conservation concern.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

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.

Small islands and large biogeographic barriers have driven contrasting speciation patterns in Indo-Pacific sunbirds (Aves: Nectariniidae)

Birds of the Indo-Pacific have provided biologists with many foundationalinsights. This study presents evidence for strong phylogeographic structure in two sunbird species from the heart of this region, the olive-backed sunbird, Cinnyris jugularis, and the black sunbird, Leptocoma aspasia. We assessed population divergence using morphological, plumage, bioacoustic and molecular data (mitochondrial ND2/ND3). Our findings indicate that the olive-backed sunbird should be recognized as multiple species, because birds from Sulawesi and the Sahul Shelf are closely related to each other, but widely separated from those in other regions. In addition, we provide evidence for an endemic species on the Wakatobi Islands, an archipelago of deep-sea islands off south-east Sulawesi. That a small bird could exhibit a range all the way from Sulawesi to Australia, while diverging on a small archipelago within this range, illustrates the complex interplay between dispersal and speciation. Our black sunbird genetic data also suggest unrecognized population structure, despite relatively weak plumage divergence. Black sunbirds in Sulawesi are likely to be a separate species from those in New Guinea, with a mean genetic distance of 9.1%. Current taxonomy suggests these sunbird species transcend classic biogeographic barriers, but our results suggest that these barriers are not easily bypassed.

Orthologous microsatellites, transposable elements, and DNA deletions correlate with generation time and body mass in neoavian birds

The rate of mutation accumulation in germline cells can be affected by cell replication and/or DNA damage, which are further related to life history traits such as generation time and body mass. Leveraging the existing datasets of 233 neoavian bird species, here, we investigated whether generation time and body mass contribute to the interspecific variation of orthologous microsatellite length, transposable element (TE) length, and deletion length and how these genomic attributes affect genome sizes. In nonpasserines, we found that generation time is correlated to both orthologous microsatellite length and TE length, and body mass is negatively correlated to DNA deletions. These patterns are less pronounced in passerines. In all species, we found that DNA deletions relate to genome size similarly as TE length, suggesting a role of body mass dynamics in genome evolution. Our results indicate that generation time and body mass shape the evolution of genomic attributes in neoavian birds.

Population genetics reveals divergent lineages and ongoing hybridization in a declining migratory fish species complex

Deciphering the effects of historical and recent demographic processes responsible for the spatial patterns of genetic diversity and structure is a key objective in evolutionary and conservation biology. Using population genetic analyses, we investigated the demographic history, the contemporary genetic diversity and structure, and the occurrence of hybridization and introgression of two species of anadromous fish with contrasting life history strategies and which have undergone recent demographic declines, the allis shad (Alosa alosa) and the twaite shad (Alosa fallax). We genotyped 706 individuals from 20 rivers and 5 sites at sea in Southern Europe at thirteen microsatellite markers. Genetic structure between populations was lower for the nearly semelparous species A. alosa, which disperses greater distances compared to the iteroparous species, A. fallax. Individuals caught at sea were assigned at the river level for A. fallax and at the region level for A. alosa. Using an approximate Bayesian computation framework, we inferred that the most likely long term historical divergence scenario between both species and lineages involved historical separation followed by secondary contact accompanied by strong population size decline. Accordingly, we found evidence for contemporary hybridization and bidirectional introgression due to gene flow between both species and lineages. Moreover, our results support the existence of at least one distinct species in the Mediterrannean sea: A. agone in Golfe du Lion area, and another divergent lineage in Corsica. Overall, our results shed light on the interplay between historical and recent demographic processes and life history strategies in shaping population genetic diversity and structure of closely related species. The recent demographic decline of these species' populations and their hybridization should be carefully considered while implementing conservation programs.

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.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Age-specific survivorship and fecundity shape genetic diversity in marine fishes

Genetic diversity varies among species due to a range of eco-evolutionary processes that are not fully understood. The neutral theory predicts that the amount of variation in the genome sequence between different individuals of the same species should increase with its effective population size (N e ). In real populations, multiple factors that modulate the variance in reproductive success among individuals causeN e to differ from the total number of individuals ( N ). Among these, age-specific mortality and fecundity rates are known to have a direct impact on theN e / N ratio. However, the extent to which vital rates account for differences in genetic diversity among species remains unknown. Here, we addressed this question by comparing genome-wide genetic diversity across 16 marine fish species with similar geographic distributions but contrasted lifespan and age-specific survivorship and fecundity curves. We sequenced the whole genome of 300 individuals to high coverage and assessed their genome-wide heterozygosity with a reference-free approach. Genetic diversity varied from 0.2% to 1.4% among species, and showed a negative correlation with adult lifespan, with a large negative effect (s l o p e = - 0.089 per additional year of lifespan) that was further increased when brooding species providing intense parental care were removed from the dataset (s l o p e = - 0.129 per additional year of lifespan). Using published vital rates for each species, we showed that theN e / N ratio resulting simply from life tables parameters can predict the observed differences in genetic diversity among species. Using simulations, we further found that the extent of reduction inN e / N with increasing adult lifespan is particularly strong under Type III survivorship curves (high juvenile and low adult mortality) and increasing fecundity with age, a typical characteristic of marine fishes. Our study highlights the importance of vital rates as key determinants of species genetic diversity levels in nature.

Response to Kratochvíl and Rovatsos

Thibault Leroy and Benoit Nabholz respond to the letter by Lukáš Kratochvíl and Michail Rovatsos that comments on the original authors' study of island songbirds.

Ratios can be misleading for detecting selection

Lukáš Kratochvíl and Michail Rovatsos respond to the recent study by Thibault Leroy and colleagues on island songbirds.

OUP accepted manuscript

The nearly neutral theory is a common framework to describe natural selection at the molecular level. This theory emphasizes the importance of slightly deleterious mutations by recognizing their ability to segregate and eventually get fixed due to genetic drift in spite of the presence of purifying selection. As genetic drift is stronger in smaller than in larger populations, a correlation between population size and molecular measures of natural selection is expected within the nearly neutral theory. However, this hypothesis was originally formulated under equilibrium conditions. As most natural populations are not in equilibrium, testing the relationship empirically may lead to confounded outcomes. Demographic nonequilibria, for instance following a change in population size, are common scenarios that are expected to push the selection–drift relationship off equilibrium. By explicitly modeling the effects of a change in population size on allele frequency trajectories in the Poisson random field framework, we obtain analytical solutions of the nonstationary allele frequency spectrum. This enables us to derive exact results of measures of natural selection and effective population size in a demographic nonequilibrium. The study of their time-dependent relationship reveals a substantial deviation from the equilibrium selection–drift balance after a change in population size. Moreover, we show that the deviation is sensitive to the combination of different measures. These results therefore constitute relevant tools for empirical studies to choose suitable measures for investigating the selection–drift relationship in natural populations. Additionally, our new modeling approach extends existing population genetics theory and can serve as foundation for methodological developments.

Life-history traits and habitat availability shape genomic diversity in birds: implications for conservation

More than 25% of species assessed by the International Union for Conservation of Nature (IUCN) are threatened with extinction. Understanding how environmental and biological processes have shaped genomic diversity may inform management practices. Using 68 extant avian species, we parsed the effects of habitat availability and life-history traits on genomic diversity over time to provide a baseline for conservation efforts. We used published whole-genome sequence data to estimate overall genomic diversity as indicated by historical long-term effective population sizes (Ne) and current genomic variability (H), then used environmental niche modelling to estimate Pleistocene habitat dynamics for each species. We found that Ne and H were positively correlated with habitat availability and related to key life-history traits (body mass and diet), suggesting the latter contribute to the overall genomic variation. We found that H decreased with increasing species extinction risk, suggesting that H may serve as a leading indicator of demographic trends related to formal IUCN conservation status in birds. Our analyses illustrate that genome-wide summary statistics estimated from sequence data reflect meaningful ecological attributes relevant to species conservation.

Adaptive divergence in bill morphology and other thermoregulatory traits is facilitated by restricted gene flow in song sparrows on the California Channel Islands

Disentangling the effects of neutral and adaptive processes in maintaining phenotypic variation across environmental gradients is challenging in natural populations. Song sparrows (Melospiza melodia) on the California Channel Islands occupy a pronounced east-west climate gradient within a small spatial scale, providing a unique opportunity to examine the interaction of genetic isolation (reduced gene flow) and the environment (selection) in driving variation. We used reduced representation genomic libraries to infer the role of neutral processes (drift and restricted gene flow) and divergent selection in driving variation in thermoregulatory traits with an emphasis on the mechanisms that maintain bill divergence among islands. Analyses of 22,029 neutral SNPs confirm distinct population structure by island with restricted gene flow and relatively large effective population sizes, suggesting bill differences are probably not a product of genetic drift. Instead, we found strong support for local adaptation using 3294 SNPs in differentiation-based and environmental association analyses coupled with genome-wide association tests. Specifically, we identified several putatively adaptive and candidate loci in or near genes involved in bill development pathways (e.g., BMP, CaM, Wnt), confirming the highly complex and polygenic architecture underlying bill morphology. Furthermore, we found divergence in genes associated with other thermoregulatory traits (i.e., feather structure, plumage colour, and physiology). Collectively, these results suggest strong divergent selection across an island archipelago results in genomic changes in a suite of traits associated with climate adaptation over small spatial scales. Future research should move beyond studying univariate traits to better understand multidimensional responses to complex environmental conditions.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

Quantifying the relationship between genetic diversity and population size suggests natural selection cannot explain Lewontin's Paradox

Neutral theory predicts that genetic diversity increases with population size, yet observed levels of diversity across metazoans vary only two orders of magnitude while population sizes vary over several. This unexpectedly narrow range of diversity is known as Lewontin's Paradox of Variation (1974). While some have suggested selection constrains diversity, tests of this hypothesis seem to fall short. Here, I revisit Lewontin's Paradox to assess whether current models of linked selection are capable of reducing diversity to this extent. To quantify the discrepancy between pairwise diversity and census population sizes across species, I combine previously-published estimates of pairwise diversity from 172 metazoan taxa with newly derived estimates of census sizes. Using phylogenetic comparative methods, I show this relationship is significant accounting for phylogeny, but with high phylogenetic signal and evidence that some lineages experience shifts in the evolutionary rate of diversity deep in the past. Additionally, I find a negative relationship between recombination map length and census size, suggesting abundant species have less recombination and experience greater reductions in diversity due to linked selection. However, I show that even assuming strong and abundant selection, models of linked selection are unlikely to explain the observed relationship between diversity and census sizes across species.

Sequential colonization of oceanic archipelagos led to a species-level radiation in the common chaffinch complex (Aves: Fringilla coelebs)

Oceanic archipelagos are excellent systems for studying speciation, yet inference of evolutionary process requires that the colonization history of island organisms be known with accuracy. Here, we used phylogenomics and patterns of genetic diversity to infer the sequence and timing of colonization of Macaronesia by mainland common chaffinches (Fringilla coelebs), and assessed whether colonization of the different archipelagos has resulted in a species-level radiation. To reconstruct the evolutionary history of the complex we generated a molecular phylogeny based on genome-wide SNP loci obtained from genotyping-by-sequencing, we ran ancestral range biogeographic analyses, and assessed fine-scale genetic structure between and within archipelagos using admixture analysis. To test for a species-level radiation, we applied a probabilistic tree-based species delimitation method (mPTP) and an integrative taxonomy approach including phenotypic differences. Results revealed a circuitous colonization pathway in Macaronesia, from the mainland to the Azores, followed by Madeira, and finally the Canary Islands. The Azores showed surprisingly high genetic diversity, similar to that found on the mainland, and the other archipelagos showed the expected sequential loss of genetic diversity. Species delimitation methods supported the existence of several species within the complex. We conclude that the common chaffinch underwent a rapid radiation across Macaronesia that was driven by the sequential colonization of the different archipelagos, resulting in phenotypically and genetically distinct, independent evolutionary lineages. We recommend a taxonomic revision of the complex that takes into account its genetic and phenotypic diversity.

The evolution of comparative phylogeography: putting the geography (and more) into comparative population genomics

Comparative population genomics is an ascendant field using genomic comparisons between species to draw inferences about forces regulating genetic variation. Comparative phylogeography, by contrast, focuses on the shared lineage histories of species codistributed geographically and is decidedly organismal in perspective. Comparative phylogeography is approximately 35 years old, and, by some metrics, is showing signs of reduced growth. Here, we contrast the goals and methods of comparative population genomics and comparative phylogeography and argue that comparative phylogeography offers an important perspective on evolutionary history that succeeds in integrating genomics with landscape evolution in ways that complement the suprageographic perspective of comparative population genomics. Focusing primarily on terrestrial vertebrates, we review the history of comparative phylogeography, its milestones and ongoing conceptual innovations, its increasingly global focus, and its status as a bridge between landscape genomics and the process of speciation. We also argue that, as a science with a strong “sense of place,” comparative phylogeography offers abundant “place-based” educational opportunities with its focus on geography and natural history, as well as opportunities for collaboration with local communities and indigenous peoples. Although comparative phylogeography does not yet require whole-genome sequencing for many of its goals, we conclude that it nonetheless plays an important role in grounding our interpretation of genetic variation in the fundamentals of geography and Earth history.

Chromosome-Level Genome Assembly of the Common Chaffinch (Aves: Fringilla coelebs): A Valuable Resource for Evolutionary Biology

The common chaffinch, Fringilla coelebs, is one of the most common, widespread, and well-studied passerines in Europe, with a broad distribution encompassing Western Europe and parts of Asia, North Africa, and the Macaronesian archipelagos. We present a high-quality genome assembly of the common chaffinch generated using Illumina shotgun sequencing in combination with Chicago and Hi-C libraries. The final genome is a 994.87-Mb chromosome-level assembly, with 98% of the sequence data located in chromosome scaffolds and a N50 statistic of 69.73 Mb. Our genome assembly shows high completeness, with a complete BUSCO score of 93.9% using the avian data set. Around 7.8% of the genome contains interspersed repetitive elements. The structural annotation yielded 17,703 genes, 86.5% of which have a functional annotation, including 7,827 complete universal single-copy orthologs out of 8,338 genes represented in the BUSCO avian data set. This new annotated genome assembly will be a valuable resource as a reference for comparative and population genomic analyses of passerine, avian, and vertebrate evolution.

Evolution: Small populations, low recombination, big trouble?

Small populations harbour less genetic diversity and more harmful mutations. They thus adapt more slowly. A new study supports these notions and suggests that reduced recombination exacerbates these effects, highlighting the impact of genome architecture on adaptability.