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.

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.

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.

A high-quality genome assembly and annotation of the dark-eyed junco Junco hyemalis, a recently diversified songbird.

The dark-eyed junco (Junco hyemalis) is one of the most common passerines of North America, and has served as a model organism in studies related to ecophysiology, behavior, and evolutionary biology for over a century. It is composed of at least 6 distinct, geographically structured forms of recent evolutionary origin, presenting remarkable variation in phenotypic traits, migratory behavior, and habitat. Here, we report a high-quality genome assembly and annotation of the dark-eyed junco generated using a combination of shotgun libraries and proximity ligation Chicago and Dovetail Hi-C libraries. The final assembly is ∼1.03 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome scaffolds, and with a N50/L50 of 71.3 Mb/5 scaffolds. We identified 19,026 functional genes combining gene prediction and similarity approaches, of which 15,967 were associated to GO terms. The genome assembly and the set of annotated genes yielded 95.4% and 96.2% completeness scores, respectively when compared with the BUSCO avian dataset. This new assembly for J. hyemalis provides a valuable resource for genome evolution analysis, and for identifying functional genes involved in adaptive processes and speciation.

Rapid evolutionary divergence of a songbird population following recent colonization of an urban area.

Colonization of a novel environment by a few individuals can lead to rapid evolutionary change, yet there is scarce evidence of the relative contributions of neutral and selective factors in promoting divergence during the early stages of colonization. Here we explore the role of neutral and selective forces in the divergence of a unique urban population of the dark-eyed junco (Junco hyemalis), which became established on the campus of the University of California at San Diego (UCSD) in the early 1980s. Previous studies based on microsatellite loci documented significant genetic differentiation of the urban population as well as divergence in phenotypic traits relative to nearby montane populations, yet the geographical origin of the colonization and the contributing factors remained uncertain. Our genome-wide single nucleotide polymorphism data set confirmed the marked genetic differentiation of the UCSD population, and we identified the coastal subspecies pinosus from central California as its sister group instead of the neighbouring mountain population. Demographic inference recovered a separation from pinosus as recent as 20-32 generations ago after a strong bottleneck, suggesting a role for drift in genetic differentiation. However, we also found significant associations between habitat variables and genome-wide variants linked to functional genes, some of which have been reported as potentially adaptive in birds inhabiting modified environments. These results suggest that the interplay between founder events and selection may result in rapid shifts in neutral and adaptive loci across the genome, and reveal the UCSD junco population as a case of contemporary evolutionary divergence in an anthropogenic environment.

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.

Exploring the vertebrate fauna of the Bird's Head Peninsula (Indonesia, West Papua) through DNA barcodes.

Biodiversity knowledge is widely heterogeneous across the Earth's biomes. Some areas, due to their remoteness and difficult access, present large taxonomic knowledge gaps. Mostly located in the tropics, these areas have frequently experienced a fast development of anthropogenic activities during the last decades and are therefore of high conservation concerns. The biodiversity hotspots of Southeast Asia exemplify the stakes faced by tropical countries. While the hotspots of Sundaland (Java, Sumatra, Borneo) and Wallacea (Sulawesi, Moluccas) have long attracted the attention of biologists and conservationists alike, extensive parts of the Sahul area, in particular the island of New Guinea, have been much less explored biologically. Here, we describe the results of a DNA-based inventory of aquatic and terrestrial vertebrate communities, which was the objective of a multidisciplinary expedition to the Bird's Head Peninsula (West Papua, Indonesia) conducted between 17 October and 20 November 2014. This expedition resulted in the assembly of 1005 vertebrate DNA barcodes. Based on the use of multiple species-delimitation methods (GMYC, PTP, RESL, ABGD), 264 molecular operational taxonomic units (MOTUs) were delineated, among which 75 were unidentified and an additional 48 were considered cryptic. This study suggests that the diversity of vertebrates of the Bird's Head is severely underestimated and considerations on the evolutionary origin and taxonomic knowledge of these biotas are discussed.

The island rule explains consistent patterns of body size evolution in terrestrial vertebrates.

Island faunas can be characterized by gigantism in small animals and dwarfism in large animals, but the extent to which this so-called 'island rule' provides a general explanation for evolutionary trajectories on islands remains contentious. Here we use a phylogenetic meta-analysis to assess patterns and drivers of body size evolution across a global sample of paired island-mainland populations of terrestrial vertebrates. We show that 'island rule' effects are widespread in mammals, birds and reptiles, but less evident in amphibians, which mostly tend towards gigantism. We also found that the magnitude of insular dwarfism and gigantism is mediated by climate as well as island size and isolation, with more pronounced effects in smaller, more remote islands for mammals and reptiles. We conclude that the island rule is pervasive across vertebrates, but that the implications for body size evolution are nuanced and depend on an array of context-dependent ecological pressures and environmental conditions.

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.

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.

Comparative Metagenomics of Palearctic and Neotropical Avian Cloacal Viromes Reveal Geographic Bias in Virus Discovery.

Our understanding about viruses carried by wild animals is still scarce. The viral diversity of wildlife may be best described with discovery-driven approaches to the study of viral diversity that broaden research efforts towards non-canonical hosts and remote geographic regions. Birds have been key organisms in the transmission of viruses causing important diseases, and wild birds are threatened by viral spillovers associated with human activities. However, our knowledge of the avian virome may be biased towards poultry and highly pathogenic diseases. We describe and compare the fecal virome of two passerine-dominated bird assemblages sampled in a remote Neotropical rainforest in French Guiana (Nouragues Natural Reserve) and a Mediterranean forest in central Spain (La Herrería). We used metagenomic data to quantify the degree of functional and genetic novelty of viruses recovered by examining if the similarity of the contigs we obtained to reference sequences differed between both locations. In general, contigs from Nouragues were significantly less similar to viruses in databases than contigs from La Herrería using Blastn but not for Blastx, suggesting that pristine regions harbor a yet unknown viral diversity with genetically more singular viruses than more studied areas. Additionally, we describe putative novel viruses of the families Picornaviridae, Reoviridae and Hepeviridae. These results highlight the importance of wild animals and remote regions as sources of novel viruses that substantially broaden the current knowledge of the global diversity of viruses.

Acceptance and knowledge of evolutionary theory among third-year university students in Spain.

The theory of evolution is one of the greatest scientific achievements in the intellectual history of humankind, yet it is still contentious within certain social groups. Despite being as robust and evidence-based as any other notable scientific theory, some people show a strong reluctance to accept it. In this study, we used the Measure of Acceptance of the Theory of Evolution (MATE) and Knowledge of Evolution Exam (KEE) questionnaires with university students from four academic degree programs (Chemistry, English, History, and Biology) of ten universities from Spain to measure, respectively, acceptance and knowledge of evolutionary theory among third-year undergraduate students (nMATE = 978; nKEE = 981). Results show that acceptance of evolution is relatively high (87.2%), whereas knowledge of the theory is moderate (5.4 out of 10) although there are differences across degrees (Biology>Chemistry>History>English), and even among various universities (ranging from 4.71 to 5.81). Statistical analysis reveals that knowledge of evolutionary theory among Biology students is partially explained by the relative weight of evolutionary themes within the curriculum, suggesting that an increase in the number of hours dedicated to this topic could have a direct influence on students' knowledge of it. We also found that religion may have a significant-although relatively small-negative influence on evolutionary theory acceptance. The moderate knowledge of evolution in our undergraduate students, together with the potential problem of acceptance in certain groups, suggests the need for a revision of the evolutionary concepts in the teaching curricula of our students since primary school.

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post-glacial radiation of the songbird genus Junco.

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome-wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark-eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome-wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.

Phylogenetic relationships and systematics of a subclade of Mesoamerican emerald hummingbirds (Aves: Trochilidae: Trochilini).

Phylogenetic relationships among major hummingbird clades are relatively well resolved, yet due to the lack of morphological synapomorphies and relative phenotypic homogeneity, the systematics of several hummingbird groups remain unresolved. Here, we present the results of a multilocus study of a clade of emerald hummingbirds composed of Cynanthus, Chlorostilbon, Riccordia and Phaeoptila (sensu Stiles et al. 2017). We include taxa not analyzed in previous studies (C. lawrencei, C. auriceps and C. forficatus, from Tres Marías Islands, Western Mexico, and Cozumel Island, respectively), and this allows us to develop a new hypothesis for the phylogenetic relationships within this group. We found that this clade originated in Mesoamerica about 12 million years ago, and comprises four geographically congruent clades: (a) the most basal clade, Phaeoptila sordida, of the Balsas River basin, Mexico; (b) Riccordia bicolor, R. maugaeus, R. ricordii and R. swainsonii of the West Indies; (c) Chlorostilbon assimilis, C. aureoventris, C. melanorhynchus, C. mellisugus, C. poortmani and C. pucherani, all of South America, except C. assimilis of Central America; and (d) a Mesoamerican clade with two Cynanthus subclades: a widespread Mexican clade that includes species formerly belonging to the genus Chlorostilbon: Cynanthus auriceps, C. canivetii, C. doubledayi and C. forficatus; and a clade restricted to the west coast of Mexico and the Tres Marías Islands and composed of C. latirostris and C. lawrencei. Our results help clarify the systematics of this group of emeralds, reconstruct its true evolutionary history, and advance understanding of phenotypic evolution in hummingbirds.

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.

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.

A Novel and Divergent Gyrovirus with Unusual Genomic Features Detected in Wild Passerine Birds from a Remote Rainforest in French Guiana.

Sequence-independent amplification techniques have become important tools for virus discovery, metagenomics, and exploration of viral diversity at the global scale, especially in remote areas. Here, we describe the detection and genetic characterization of a novel gyrovirus, named GyV11, present in cloacal, oral, and blood samples from neotropical wild birds in French Guiana. The molecular epidemiology revealed the presence of GyV11 only in passerine birds from three different species at a low prevalence (0.73%). This is the first characterization and prevalence study of a gyrovirus carried out in resident wild bird populations in a remote region, and provides evidence of the fecal-oral route transmission and local circulation of the virus. The molecular phylogeny of gyroviruses reveals the existence of two distinct gyrovirus lineages in which GyV11 is phylogenetically distinct from previously reported gyroviruses. Furthermore, GyV11 is placed basal in the gyrovirus phylogeny, likely owing to its ancestral origin and marked divergence. This study also provides important insights into the ecology, epidemiology, and genomic features of gyroviruses in a remote neotropical rainforest. The pathogenesis of this virus in avian species or whether GyV11 can infect humans and/or chickens needs to be further investigated.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird.

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

A comparison of genomic islands of differentiation across three young avian species pairs.

Detailed evaluations of genomic variation between sister species often reveal distinct chromosomal regions of high relative differentiation (i.e., "islands of differentiation" in FST ), but there is much debate regarding the causes of this pattern. We briefly review the prominent models of genomic islands of differentiation and compare patterns of genomic differentiation in three closely related pairs of New World warblers with the goal of evaluating support for the four models. Each pair (MacGillivray's/mourning warblers; Townsend's/black-throated green warblers; and Audubon's/myrtle warblers) consists of forms that were likely separated in western and eastern North American refugia during cycles of Pleistocene glaciations and have now come into contact in western Canada, where each forms a narrow hybrid zone. We show strong differences between pairs in their patterns of genomic heterogeneity in FST , suggesting differing selective forces and/or differing genomic responses to similar selective forces among the three pairs. Across most of the genome, levels of within-group nucleotide diversity (πWithin ) are almost as large as levels of between-group nucleotide distance (πBetween ) within each pair, suggesting recent common ancestry and/or gene flow. In two pairs, a pattern of the FST peaks having low πBetween suggests that selective sweeps spread between geographically differentiated groups, followed by local differentiation. This "sweep-before-differentiation" model is consistent with signatures of gene flow within the yellow-rumped warbler species complex. These findings add to our growing understanding of speciation as a complex process that can involve phases of adaptive introgression among partially differentiated populations.

Differential gene regulation underlies variation in melanic plumage coloration in the dark-eyed junco (Junco hyemalis).

Colour plays a prominent role in species recognition; therefore, understanding the proximate basis of pigmentation can provide insight into reproductive isolation and speciation. Colour differences between taxa may be the result of regulatory differences or be caused by mutations in coding regions of the expressed genes. To investigate these two alternatives, we studied the pigment composition and the genetic basis of coloration in two divergent dark-eyed junco (Junco hyemalis) subspecies, the slate-coloured and Oregon juncos, which have evolved marked differences in plumage coloration since the Last Glacial Maximum. We used HPLC and light microscopy to investigate pigment composition and deposition in feathers from four body areas. We then used RNA-seq to compare the relative roles of differential gene expression in developing feathers and sequence divergence in transcribed loci under common-garden conditions. Junco feathers differed in eumelanin and pheomelanin content and distribution. Within subspecies, in lighter feathers melanin synthesis genes were downregulated (including PMEL, TYR, TYRP1, OCA2 and MLANA), and ASIP was upregulated. Feathers from different body regions also showed differential expression of HOX and WNT genes. Feathers from the same body regions that differed in colour between the two subspecies showed differential expression of ASIP and three other genes (MFSD12, KCNJ13 and HAND2) associated with pigmentation in other taxa. Sequence variation in the expressed genes was not related to colour differences. Our findings support the hypothesis that differential regulation of a few genes can account for marked differences in coloration, a mechanism that may facilitate the rapid phenotypic diversification of juncos.