Survival and divergence in a small group: the extraordinary genomic history of the endangered Apennine brown bear stragglers 2 AbstractAbout 100 km east of Rome, in the Central Apennine mountains, a critically endangered population of approximately fifty brown bears live in complete isolation. Mating outside this population is prevented by several hundred kilometers of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. First, we found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In Central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show non-random peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles a) can be an important driver of divergence in isolation, b) can be tolerated when balancing selection prevents random loss of variation at important genes and c) is followed by or results directly in favorable behavioral changes. SignificanceA small and relict population of brown bears lives in complete isolation in the Italian Apennine mountains, providing a unique opportunity to study the impact of drift and selection on the genomes of a large endangered mammal and to reconstruct the phenotypic consequences and the conservation implications of such evolutionary processes. The Apennine bear is highly inbred and harbors very low genomic variation. Several deleterious mutations have been accumulated by drift. We found evidence that this is a consequence of habitat fragmentation in the Neolithic, when human expansion and land clearance shrank its habitat, and that retention of variation at immune system and olfactory receptor genes, as well as changes in diet and behavior, prevented the extinction of the Apennine bear.

Among birds, white-eyes (genus Zosterops) have diversified so extensively that Jared Diamond and Ernst Mayr referred to them as the “great speciator.” The Zosterops lineage exhibits some of the fastest rates of species diversification among vertebrates, and its members are the most prolific passerine island colonizers. We present a high-quality genome assembly for the silvereye (Zosterops lateralis), a white-eye species consisting of several subspecies distributed across multiple islands. We investigate the genetic basis of rapid diversification in white-eyes by conducting genomic analyses at varying taxonomic levels. First, we compare the silvereye genome with those of birds from different families and searched for genomic features that may be unique to Zosterops. Second, we compare the genomes of different species of white-eyes from Lifou island (South Pacific), using whole genome resequencing and restriction site associated DNA. Third, we contrast the genomes of two subspecies of silvereye that differ in plumage color. In accordance with theory, we show that white-eyes have high rates of substitutions, gene duplication, and positive selection relative to other birds. Below genus level, we find that genomic differentiation accumulates rapidly and reveals contrasting demographic histories between sympatric species on Lifou, indicative of past interspecific interactions. Finally, we highlight genes possibly involved in color polymorphism between the subspecies of silvereye. By providing the first whole-genome sequence resources for white-eyes and by conducting analyses at different taxonomic levels, we provide genomic evidence underpinning this extraordinary bird radiation.

In the post-genomic era, much of phylogenetic analyses still relies on mitochondrial DNA, either alone or in combination with few nuclear genes. Although this approach often makes it possible to construct well-supported trees, it is limited because mtDNA describes the history of a single locus, and nuclear phylogenies based on a few loci may be biased, leading to inaccurate tree topologies and biased estimations of species divergence time. In this study, we perform a phylogenomic analysis of the Daphniidae family (Crustacea: Branchiopoda: Anomopoda) including some of the most frequently studied model organisms (Daphnia magna and D. pulex) whose phylogenetic relationships have been based primarily on an assessment of a few mtDNA genes. Using high-throughput sequencing, we were able to assemble 38 whole mitochondrial genomes and draft nuclear genomes for 18 species, including at least one species for each known genus of the family Daphniidae. Here we present phylogenies based on 636 nuclear single-copy genes shared among all sampled taxa and based on whole mtDNA genomes. The phylogenies we obtained were highly supported and showed some discrepancies between nuclear and mtDNA based trees at deeper nodes. We also identified a new candidate sister lineage of Daphnia magna. Our time-calibrated genomic trees, which we constructed using both fossil records and substitution rates, yielded very different estimates of branching event times compared to those based on mtDNA. By providing multi-locus, fossil-calibrated trees of the Daphniidae, our study contributes to an improved phylogenetic framework for ecological and evolutionary studies that use water fleas as a model system.

Spatial variation in pathogen-mediated selection is predicted to influence the evolutionary trajectory of host populations and lead to spatial variation in their immunogenetic composition. However, to date few studies have been able to directly link small-scale spatial variation in infection risk to host immune gene evolution in natural, nonhuman populations. Here, we use a natural rodent-Borrelia system to test for associations between landscape-level spatial variation in Borrelia infection risk along replicated elevational gradients in the Swiss Alps and Toll-like receptor 2 (TLR2) evolution, a candidate gene for Borrelia resistance, across bank vole (Myodes glareolus) populations. We found that Borrelia infection risk (i.e., the product of Borrelia prevalence in questing ticks and the average tick load of voles at a sampling site) was spatially variable and significantly negatively associated with elevation. Across sampling sites, Borrelia prevalence in bank voles was significantly positively associated with Borrelia infection risk along the elevational clines. We observed a significant association between naturally occurring TLR2 polymorphisms in hosts and their Borrelia infection status. The TLR2 variant associated with a reduced likelihood of Borrelia infection was most common in rodent populations at lower elevations that face a high Borrelia infection risk, and its frequency changed in accordance with the change in Borrelia infection risk along the elevational clines. These results suggest that small-scale spatial variation in parasite-mediated selection affects the immunogenetic composition of natural host populations, providing a striking example that the microbial environment shapes the evolution of the host's immune system in the wild.