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

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.

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.

Beyond morphs: Inter-individual colour variation despite strong genetic determinism of colour morphs in a wild bird

Categorizing individuals into discrete forms in colour polymorphic species can overlook more subtle patterns in coloration that can be of functional significance. Thus, quantifying inter-individual variation in these species at both within- and between-morph levels is critical to understand the evolution of colour polymorphisms. Here we present analyses of inter-individual colour variation in the Reunion grey white-eye (Zosterops borbonicus), a colour polymorphic wild bird endemic to the island of Reunion in which all highland populations contain two sympatric colour morphs, with birds showing predominantly grey or brown plumage, respectively. We first quantified colour variation across multiple body areas by using a continuous plumage colour score to assess variation in brown-grey coloration as well as smaller scale variation in light patches. To examine the possible causes of among-individual variation, we tested if colour variation in plumage component elements could be explained by genotypes at two markers near a major-effect locus previously related to back coloration in this species, and by other factors such as age, sex and body condition. Overall, grey-brown coloration was largely determined by genetic factors and was best described by three distinct clusters that were associated to genotypic classes (homozygotes and heterozygote), with no effect of age or sex, whereas variation in smaller light patches was primarily related to age and sex. Our results highlight the importance of characterizing subtle plumage variation beyond morph categories that are readily observable since multiple patterns of colour variation may be driven by different mechanisms, have different functions and will likely respond in different ways to selection.

Selective sweeps on different pigmentation genes mediate convergent evolution of island melanism in two incipient bird species

Insular organisms often evolve predictable phenotypes, like flightlessness, extreme body sizes, or increased melanin deposition. The evolutionary forces and molecular targets mediating these patterns remain mostly unknown. Here we study the Chestnut-bellied Monarch (Monarcha castaneiventris) from the Solomon Islands, a complex of closely related subspecies in the early stages of speciation. On the large island of Makira M. c. megarhynchus has a chestnut belly, whereas on the small satellite islands of Ugi, and Santa Ana and Santa Catalina (SA/SC) M. c. ugiensis is entirely iridescent blue-black (i.e., melanic). Melanism has likely evolved twice, as the Ugi and SA/SC populations were established independently. To investigate the genetic basis of melanism on each island we generated whole genome sequence data from all three populations. Non-synonymous mutations at the MC1R pigmentation gene are associated with melanism on SA/SC, while ASIP, an antagonistic ligand of MC1R, is associated with melanism on Ugi. Both genes show evidence of selective sweeps in traditional summary statistics and statistics derived from the ancestral recombination graph (ARG). Using the ARG in combination with machine learning, we inferred selection strength, timing of onset and allele frequency trajectories. MC1R shows evidence of a recent, strong, soft selective sweep. The region including ASIP shows more complex signatures; however, we find evidence for sweeps in mutations near ASIP, which are comparatively older than those on MC1R and have been under relatively strong selection. Overall, our study shows convergent melanism results from selective sweeps at independent molecular targets, evolving in taxa where coloration likely mediates reproductive isolation with the neighboring chestnut-bellied subspecies.

Chestnut-bellied Monarchs (Monarcha castaneiventris ugiensis) from two archipelagos in the Solomon Islands have evolved entirely black plumage from a chestnut ancestor (Monarcha castaneiventris megarhynchus), a phenomenon known as island melanism. We obtain and analyze whole genome sequences using traditional summary statistics and new methods that combine inference of the ancestral recombination graph with machine learning. We find multiple lines of evidence for independent selective sweeps on the MC1R and ASIP genes, a receptor/ligand pair which regulates the production of melanin. Melanism on each archipelago is mediated by mutations in one of these two genes. Mutations in and around MC1R underwent a recent soft sweep experiencing strong selection on the islands of Santa Ana and Santa Catalina, whereas selection was also strong but comparatively older for ASIP on the island of Ugi. We show how melanism originated under positive selection on independent molecular targets, evolving convergently in taxa where coloration mediates reproductive isolation.

bric à brac controls sex pheromone choice by male European corn borer moths

The sex pheromone system of ~160,000 moth species acts as a powerful form of assortative mating whereby females attract conspecific males with a species-specific blend of volatile compounds. Understanding how female pheromone production and male preference coevolve to produce this diversity requires knowledge of the genes underlying change in both traits. In the European corn borer moth, pheromone blend variation is controlled by two alleles of an autosomal fatty-acyl reductase gene expressed in the female pheromone gland (pgFAR). Here we show that asymmetric male preference is controlled by cis-acting variation in a sex-linked transcription factor expressed in the developing male antenna, bric à brac (bab). A genome-wide association study of preference using pheromone-trapped males implicates variation in the 293 kb bab intron 1, rather than the coding sequence. Linkage disequilibrium between bab intron 1 and pgFAR further validates bab as the preference locus, and demonstrates that the two genes interact to contribute to assortative mating. Thus, lack of physical linkage is not a constraint for coevolutionary divergence of female pheromone production and male behavioral response genes, in contrast to what is often predicted by evolutionary theory.

Island songbirds as windows into evolution in small populations

Due to their limited ranges and inherent isolation, island species have long been recognized as crucial systems for tackling a range of evolutionary questions, including in the early study of speciation.^1,2 Such species have been less studied in the understanding of the evolutionary forces driving DNA sequence evolution. Island species usually have lower census population sizes (N) than continental species and, supposedly, lower effective population sizes (Ne). Given that both the rates of change caused by genetic drift and by selection are dependent upon Ne, island species are theoretically expected to exhibit (1) lower genetic diversity, (2) less effective natural selection against slightly deleterious mutations,^3,4 and (3) a lower rate of adaptive evolution.^5-8 Here, we have used a large set of newly sequenced and published whole-genome sequences of Passerida species (14 insular and 11 continental) to test these predictions. We confirm that island species exhibit lower census size and Ne, supporting the hypothesis that the smaller area available on islands constrains the upper bound of Ne. In the insular species, we find lower nucleotide diversity in coding regions, higher ratios of non-synonymous to synonymous polymorphisms, and lower adaptive substitution rates. Our results provide robust evidence that the lower Ne experienced by island species has affected both the ability of natural selection to efficiently remove weakly deleterious mutations and also the adaptive potential of island species, therefore providing considerable empirical support for the nearly neutral theory. We discuss the implications for both evolutionary and conservation biology.

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.

Tissue-specific geometry and chemistry of modern and fossilized melanosomes reveal internal anatomy of extinct vertebrates

Recent reports of nonintegumentary melanosomes in fossils hint at functions for melanin beyond color production, but the biology and evolution of internal melanins are poorly understood. Our results show that internal melanosomes are widespread in diverse fossil and modern vertebrates and have tissue-specific geometries and metal chemistries. Tissue-specific chemical signatures can persist in fossils despite some diagenetic overprint, allowing the reconstruction of internal soft-tissue anatomy in fossil vertebrates, and suggest that links between melanin and metal regulation have deep evolutionary origins in vertebrates.

Recent discoveries of nonintegumentary melanosomes in extant and fossil amphibians offer potential insights into the physiological functions of melanin not directly related to color production, but the phylogenetic distribution and evolutionary history of these internal melanosomes has not been characterized systematically. Here, we present a holistic method to discriminate among melanized tissues by analyzing the anatomical distribution, morphology, and chemistry of melanosomes in various tissues in a phylogenetically broad sample of extant and fossil vertebrates. Our results show that internal melanosomes in all extant vertebrates analyzed have tissue-specific geometries and elemental signatures. Similar distinct populations of preserved melanosomes in phylogenetically diverse vertebrate fossils often map onto specific anatomical features. This approach also reveals the presence of various melanosome-rich internal tissues in fossils, providing a mechanism for the interpretation of the internal anatomy of ancient vertebrates. Collectively, these data indicate that vertebrate melanins share fundamental physiological roles in homeostasis via the scavenging and sequestering of metals and suggest that intimate links between melanin and metal metabolism in vertebrates have deep evolutionary origins.

Melanin-based coloration and host-parasite interactions under global change

The role of parasites in shaping melanin-based colour polymorphism, and the consequences of colour polymorphism for disease resistance, remain debated. Here we review recent evidence of the links between melanin-based coloration and the behavioural and immunological defences of vertebrates against their parasites. First we propose that (1) differences between colour morphs can result in variable exposure to parasites, either directly (certain colours might be more or less attractive to parasites) or indirectly (variations in behaviour and encounter probability). Once infected, we propose that (2) immune variation between differently coloured individuals might result in different abilities to cope with parasite infection. We then discuss (3) how these different abilities could translate into variable sexual and natural selection in environments varying in parasite pressure. Finally, we address (4) the potential role of parasites in the maintenance of melanin-based colour polymorphism, especially in the context of global change and multiple stressors in human-altered environments. Because global change will probably affect both coloration and the spread of parasitic diseases in the decades to come, future studies should take into account melanin-based coloration to better predict the evolutionary responses of animals to changing disease risk in human-altered environments.

Heritability of plumage colour morph variation in a wild population of promiscuous, long-lived Australian magpies

Colour polymorphisms have evolutionary significance for the generation and maintenance of species diversity. Demonstrating heritability of polymorphic traits can be challenging for wild populations of long-lived species because accurate information is required on trait expression and familial relationships. The Australian magpie Cracticus tibicen has a continent-wide distribution featuring several distinct plumage morphs, differing primarily in colour of back feathers. Black or white-backed morphs occur in eastern Australia, with intermediate morphs common in a narrow hybrid zone where the two morphs meet. This study investigated heritability of back colour phenotypes in a hybrid zone population (Seymour, Victoria) based on long-term observational data and DNA samples collected over an 18 year period (1993-2010). High extra-pair paternity (~ 36% offspring), necessitated verification of parent-offspring relationships by parentage analysis. A total of 538 birds (221 parents and 317 offspring) from 36 territories were analysed. Back colour was a continuous trait scored on a five-morph scale in the field (0-4). High and consistent estimates of back colour heritability (h²) were obtained via weighted mid-parent regression (h² = 0.94) and by animal models (h² = 0.92, C.I. 0.80-0.99). Single-parent heritability estimates indicated neither maternal nor paternal non-genetic effects (e.g., parent body condition) played a large role in determining offspring back colour, and environmental effects of territory group and cohort contributed little to trait heritability. Distinctive back colouration of the Australian magpie behaves as a quantitative trait that is likely polygenic, although mechanisms responsible for maintaining these geographically structured morphs and the hybrid zone where they meet are unknown.

Admixture mapping in a hybrid zone reveals loci associated with avian feather coloration

Identifying the genetic bases for colour patterns has provided important insights into the control and expression of pigmentation and how these characteristics influence fitness. However, much more is known about the genetic bases for traits based on melanin pigments than for traits based on another major class of pigments, carotenoids. Here, we use natural admixture in a hybrid zone between Audubon's and myrtle warblers (Setophaga coronata auduboni/S. c. coronata) to identify genomic regions associated with both types of pigmentation. Warblers are known for rapid speciation and dramatic differences in plumage. For each of five plumage coloration traits, we found highly significant associations with multiple single-nucleotide polymorphisms (SNPs) across the genome and these were clustered in discrete regions. Regions near significantly associated SNPs were enriched for genes associated with keratin filaments, fibrils that make up feathers. A carotenoid-based trait that differs between the taxa-throat colour-had more than a dozen genomic regions of association. One cluster of SNPs for this trait overlaps the Scavenger Receptor Class F Member 2 (SCARF2) gene. Other scavenger receptors are presumed to be expressed at target tissues and involved in the selective movement of carotenoids into the target cells, making SCARF2 a plausible new candidate for carotenoid processing. In addition, two melanin-based plumage traits-colours of the eye line and eye spot-show very strong associations with a single genomic region mapping to chromosome 20 in the zebra finch. These findings indicate that only a subset of the genomic regions differentiated between these two warblers are associated with the plumage differences between them and demonstrate the utility of reduced-representation genomic scans in hybrid zones.

Linking genotype to phenotype in a changing ocean: inferring the genomic architecture of a blue mussel stress response with genome-wide association

A key component to understanding the evolutionary response to a changing climate is linking underlying genetic variation to phenotypic variation in stress response. Here, we use a genome-wide association approach (GWAS) to understand the genetic architecture of calcification rates under simulated climate stress. We take advantage of the genomic gradient across the blue mussel hybrid zone (Mytilus edulis and Mytilus trossulus) in the Gulf of Maine (GOM) to link genetic variation with variance in calcification rates in response to simulated climate change. Falling calcium carbonate saturation states are predicted to negatively impact many marine organisms that build calcium carbonate shells - like blue mussels. We sampled wild mussels and measured net calcification phenotypes after exposing mussels to a 'climate change' common garden, where we raised temperature by 3°C, decreased pH by 0.2 units and limited food supply by filtering out planktonic particles >5 μm, compared to ambient GOM conditions in the summer. This climate change exposure greatly increased phenotypic variation in net calcification rates compared to ambient conditions. We then used regression models to link the phenotypic variation with over 170 000 single nucleotide polymorphism loci (SNPs) generated by genotype by sequencing to identify genomic locations associated with calcification phenotype, and estimate heritability and architecture of the trait. We identified at least one of potentially 2-10 genomic regions responsible for 30% of the phenotypic variation in calcification rates that are potential targets of natural selection by climate change. Our simulations suggest a power of 13.7% with our study's average effective sample size of 118 individuals and rare alleles, but a power of >90% when effective sample size is 900.

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.