The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

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.

Deep ecomorphological and genetic divergence in Steller's Jays (Cyanocitta stelleri, Aves: Corvidae)

The relationship between ecology and morphology is a cornerstone of evolutionary biology, and quantifying variation across environments can shed light on processes that give rise to biodiversity. Three morphotypes of the Steller's Jay (Cyanocitta stelleri) occupy different ecoregions in western North America, which vary in climate and landcover. These morphotypes (Coastal, Interior, Rocky Mountain) differ in size, plumage coloration, and head pattern. We sampled 1080 Steller's Jays from 68 populations (plus 11 outgroups) to address three main questions using data on morphology, plumage, genetics (mtDNA, microsatellites), and ecological niches: (1) How do phenotypic and genetic traits vary within and among populations, morphotypes, and ecoregions? (2) How do population‐level differences in Steller's Jays compare with other sister species pairs of North American birds? (3) What can we infer about the population history of Steller's Jays in relation to past climates, paleoecology, and niche evolution? We found substantial morphological, genetic, and ecological differentiation among morphotypes. The greatest genetic divergence separated Coastal and Interior morphotypes from the Rocky Mountain morphotype, which was associated with warmer, drier, and more open habitats. Microsatellites revealed additional structure between Coastal and Interior groups. The deep mtDNA split between Coastal/Interior and Rocky Mountain lineages of Steller's Jay (ND2 ~ 7.8%) is older than most North American avian sister species and dates to approximately 4.3 mya. Interior and Rocky Mountain morphotypes contact across a narrow zone with steep clines in traits and reduced gene flow. The distribution of the three morphotypes coincides with divergent varieties of ponderosa pine and Douglas fir. Species distribution models support multiple glacial refugia for Steller's Jays. Our integrative dataset combined with extensive geographic sampling provides compelling evidence for recognizing at least two species of Steller's Jay.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250 km² island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA) similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3,000 single nucleotide polymorphisms (SNPs) in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing (RADseq). Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis (RDA) while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Genomic variation, population history and within-archipelago adaptation between island bird populations

Oceanic island archipelagos provide excellent models to understand evolutionary processes. Colonization events and gene flow can interact with selection to shape genetic variation at different spatial scales. Landscape-scale variation in biotic and abiotic factors may drive fine-scale selection within islands, while long-term evolutionary processes may drive divergence between distantly related populations. Here, we examine patterns of population history and selection between recently diverged populations of the Berthelot's pipit (Anthus berthelotii), a passerine endemic to three North Atlantic archipelagos. First, we use demographic trees and f3 statistics to show that genome-wide divergence across the species range is largely shaped by colonization and bottlenecks, with evidence of very weak gene flow between populations. Then, using a genome scan approach, we identify signatures of divergent selection within archipelagos at single nucleotide polymorphisms (SNPs) in genes potentially associated with craniofacial development and DNA repair. We did not detect within-archipelago selection at the same SNPs as were detected previously at broader spatial scales between archipelagos, but did identify signatures of selection at loci associated with similar biological functions. These findings suggest that similar ecological factors may repeatedly drive selection between recently separated populations, as well as at broad spatial scales across varied landscapes.

Differential divergence in autosomes and sex chromosomes is associated with intra-island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.

Within-island diversification in a passerine bird

The presence of congeneric taxa on the same island suggests the possibility of in situ divergence, but can also result from multiple colonizations of previously diverged lineages. Here, using genome-wide data from a large population sample, we test the hypothesis that intra-island divergence explains the occurrence of four geographical forms meeting at hybrid zones in the Reunion grey white-eye (Zosterops borbonicus), a species complex endemic to the small volcanic island of Reunion. Using population genomic and phylogenetic analyses, we reconstructed the population history of the different forms. We confirmed the monophyly of the complex and found that one of the lowland forms is paraphyletic and basal relative to others, a pattern highly consistent with in situ divergence. Our results suggest initial colonization of the island through the lowlands, followed by expansion into the highlands, which led to the evolution of a distinct geographical form, genetically and ecologically different from the lowland ones. Lowland forms seem to have experienced periods of geographical isolation, but they diverged from one another by sexual selection rather than niche change. Overall, low dispersal capabilities in this island bird combined with both geographical and ecological opportunities seem to explain how divergence occurred at such a small spatial scale.

The Multispecies Coalescent Over-Splits Species in the Case of Geographically Widespread Taxa

Many recent species delimitation studies rely exclusively on limited analyses of genetic data analyzed under the multispecies coalescent (MSC) model, and results from these studies often are regarded as conclusive support for taxonomic changes. However, most MSC-based species delimitation methods have well-known and often unmet assumptions. Uncritical application of these genetic-based approaches (without due consideration of sampling design, the effects of a priori group designations, isolation by distance, cytoplasmic-nuclear mismatch, and population structure) can lead to over-splitting of species. Here, we argue that in many common biological scenarios, researchers must be particularly cautious regarding these limitations, especially in cases of well-studied, geographically variable, and parapatrically distributed species complexes. We consider these points with respect to a historically controversial species group, the American milksnakes (Lampropeltis triangulum complex), using genetic data from a recent analysis (Ruane et al. 2014). We show that over-reliance on the program Bayesian Phylogenetics and Phylogeography, without adequate consideration of its assumptions and of sampling limitations, resulted in over-splitting of species in this study. Several of the hypothesized species of milksnakes instead appear to represent arbitrary slices of continuous geographic clines. We conclude that the best available evidence supports three, rather than seven, species within this complex. More generally, we recommend that coalescent-based species delimitation studies incorporate thorough analyses of geographic variation and carefully examine putative contact zones among delimited species before making taxonomic changes.

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P ST) which largely exceeds neutral genetic differentiation (F ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.

Gene flow does not prevent personality and morphological differentiation between two blue tit populations

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait-specific Q_st and F_st . Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Q_st -F_st comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Q_st -F_st comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a 'brown' ancestral population, the dominant 'grey' allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.