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.

How structural variants shape avian phenotypes: Lessons from model systems.

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Corticosterone exposure is associated with long-term changes in DNA methylation, physiology and breeding decisions in a wild bird.

When facing challenges, vertebrates activate a hormonal stress response that can dramatically alter behaviour and physiology. Although this response can be costly, conceptual models suggest that it can also recalibrate the stress response system, priming more effective responses to future challenges. Little is known about whether this process occurs in wild animals, particularly in adulthood, and if so, how information about prior experience with stressors is encoded. One potential mechanism is hormonally mediated changes in DNA methylation. We simulated the spikes in corticosterone that accompany a stress response using non-invasive dosing in tree swallows (Tachycineta bicolor) and monitored the phenotypic effects 1 year later. In a subset of individuals, we characterized DNA methylation using reduced representation bisulfite sequencing shortly after treatment and a year later. The year after treatment, experimental females had stronger negative feedback and initiated breeding earlier-traits that are associated with stress resilience and reproductive performance in our population-and higher baseline corticosterone. We also found that natural variation in corticosterone predicted patterns of DNA methylation. Finally, corticosterone treatment influenced methylation on short (1-2 weeks) and long (1 year) time scales; however, these changes did not have clear links to functional regulation of the stress response. Taken together, our results are consistent with corticosterone-induced priming of future stress resilience and support DNA methylation as a potential mechanism, but more work is needed to demonstrate functional consequences. Uncovering the mechanisms linking experience with the response to future challenges has implications for understanding the drivers of stress resilience.

A Deep-Learning Approach for Inference of Selective Sweeps from the Ancestral Recombination Graph.

Detecting signals of selection from genomic data is a central problem in population genetics. Coupling the rich information in the ancestral recombination graph (ARG) with a powerful and scalable deep-learning framework, we developed a novel method to detect and quantify positive selection: Selection Inference using the Ancestral recombination graph (SIA). Built on a Long Short-Term Memory (LSTM) architecture, a particular type of a Recurrent Neural Network (RNN), SIA can be trained to explicitly infer a full range of selection coefficients, as well as the allele frequency trajectory and time of selection onset. We benchmarked SIA extensively on simulations under a European human demographic model, and found that it performs as well or better as some of the best available methods, including state-of-the-art machine-learning and ARG-based methods. In addition, we used SIA to estimate selection coefficients at several loci associated with human phenotypes of interest. SIA detected novel signals of selection particular to the European (CEU) population at the MC1R and ABCC11 loci. In addition, it recapitulated signals of selection at the LCT locus and several pigmentation-related genes. Finally, we reanalyzed polymorphism data of a collection of recently radiated southern capuchino seedeater taxa in the genus Sporophila to quantify the strength of selection and improved the power of our previous methods to detect partial soft sweeps. Overall, SIA uses deep learning to leverage the ARG and thereby provides new insight into how selective sweeps shape genomic diversity.

Genomic islands of differentiation in a rapid avian radiation have been driven by recent selective sweeps.

Numerous studies of emerging species have identified genomic "islands" of elevated differentiation against a background of relative homogeneity. The causes of these islands remain unclear, however, with some signs pointing toward "speciation genes" that locally restrict gene flow and others suggesting selective sweeps that have occurred within nascent species after speciation. Here, we examine this question through the lens of genome sequence data for five species of southern capuchino seedeaters, finch-like birds from South America that have undergone a species radiation during the last ∼50,000 generations. By applying newly developed statistical methods for ancestral recombination graph inference and machine-learning methods for the prediction of selective sweeps, we show that previously identified islands of differentiation in these birds appear to be generally associated with relatively recent, species-specific selective sweeps, most of which are predicted to be soft sweeps acting on standing genetic variation. Many of these sweeps coincide with genes associated with melanin-based variation in plumage, suggesting a prominent role for sexual selection. At the same time, a few loci also exhibit indications of possible selection against gene flow. These observations shed light on the complex manner in which natural selection shapes genome sequences during speciation.

Concerted variation in melanogenesis genes underlies emergent patterning of plumage in capuchino seedeaters.

Coloration traits are central to animal communication; they often govern mate choice, promote reproductive isolation and catalyse speciation. Specific genetic changes can cause variation in coloration, yet far less is known about how overall coloration patterns-which involve combinations of multiple colour patches across the body-can arise and are genomically controlled. We performed genome-wide association analyses to link genomic changes to variation in melanin (eumelanin and pheomelanin) concentration in feathers from different body parts in the capuchino seedeaters, an avian radiation with diverse colour patterns despite remarkably low genetic differentiation across species. Cross-species colour variation in each plumage patch is associated with unique combinations of variants at a few genomic regions, which include mostly non-coding (presumably regulatory) areas close to known pigmentation genes. Genotype-phenotype associations can vary depending on patch colour and are stronger for eumelanin pigmentation, suggesting eumelanin production is tightly regulated. Although some genes are involved in colour variation in multiple patches, in some cases, the SNPs associated with colour changes in different patches segregate spatially. These results suggest that coloration patterning in capuchinos is generated by the modular combination of variants that regulate multiple melanogenesis genes, a mechanism that may have promoted this rapid radiation.

Genomic phylogeography of the White-crowned Manakin Pseudopipra pipra (Aves: Pipridae) illuminates a continental-scale radiation out of the Andes.

The complex landscape history of the Neotropics has generated opportunities for population isolation and diversification that place this region among the most species-rich in the world. Detailed phylogeographic studies are required to uncover the biogeographic histories of Neotropical taxa, to identify evolutionary correlates of diversity, and to reveal patterns of genetic connectivity, disjunction, and potential differentiation among lineages from different areas of endemism. The White-crowned Manakin (Pseudopipra pipra) is a small suboscine passerine bird that is broadly distributed through the subtropical rainforests of Central America, the lower montane cloud forests of the Andes from Colombia to central Peru, the lowlands of Amazonia and the Guianas, and the Atlantic forest of southeast Brazil. Pseudopipra is currently recognized as a single, polytypic biological species. We studied the effect of the Neotropical landscape on genetic and phenotypic differentiation within this species using genomic data derived from double digest restriction site associated DNA sequencing (ddRAD), and mitochondrial DNA. Most of the genetic breakpoints we identify among populations coincide with physical barriers to gene flow previously associated with avian areas of endemism. The phylogenetic relationships among these populations imply a novel pattern of Andean origination for this group, with subsequent diversification into the Amazonian lowlands. Our analysis of genomic admixture and gene flow reveals a complex history of introgression between some western Amazonian populations. These reticulate processes confound our application of standard concatenated and coalescent phylogenetic methods and raise the question of whether a lineage in the western Napo area of endemism should be considered a hybrid species. Lastly, analysis of variation in vocal and plumage phenotypes in the context of our phylogeny supports the hypothesis that Pseudopipra is a species-complex composed of at least 8, and perhaps up to 17 distinct species which have arisen in the last ~2.5 Ma.

Variable Signatures of Selection Despite Conserved Recombination Landscapes Early in Speciation.

Recently diverged taxa often exhibit heterogeneous landscapes of genomic differentiation, characterized by regions of elevated differentiation on an otherwise homogeneous background. While divergence peaks are generally interpreted as regions responsible for reproductive isolation, they can also arise due to background selection, selective sweeps unrelated to speciation, and variation in recombination and mutation rates. To investigate the association between patterns of recombination and landscapes of genomic differentiation during the early stages of speciation, we generated fine-scale recombination maps for six southern capuchino seedeaters (Sporophila) and two subspecies of White Wagtail (Motacilla alba), two recent avian radiations in which divergent selection on pigmentation genes has likely generated peaks of differentiation. We compared these recombination maps to those of Collared (Ficedula albicollis) and Pied Flycatchers (Ficedula hypoleuca), non-sister taxa characterized by moderate genomic divergence and a heterogenous landscape of genomic differentiation shaped in part by background selection. Although recombination landscapes were conserved within all three systems, we documented a weaker negative correlation between recombination rate and genomic differentiation in the recent radiations. All divergence peaks between capuchinos, wagtails, and flycatchers were located in regions with lower-than-average recombination rates, and most divergence peaks in capuchinos and flycatchers fell in regions of exceptionally reduced recombination. Thus, co-adapted allelic combinations in these regions may have been protected early in divergence, facilitating rapid diversification. Despite largely conserved recombination landscapes, divergence peaks are specific to each focal comparison in capuchinos, suggesting that regions of elevated differentiation have not been generated by variation in recombination rate alone.

A test of the riverine barrier hypothesis in the largest subtropical river basin in the Neotropics.

The riverine barrier hypothesis proposes that large rivers represent geographical barriers to gene flow for terrestrial organisms, leading to population differentiation and ultimately allopatric speciation. Here we assess for the first time if the subtropical Paraná-Paraguay River system in the Del Plata basin, second in size among South American drainages, acts as a barrier to gene flow for birds. We analysed the degree of mitochondrial and nuclear genomic differentiation in seven species with known subspecies divided by the Paraná-Paraguay River axis. Only one species showed genetic differentiation concordant with the current river channel, but another five species have an east/west genetic split broadly coincident with the Paraná River's dynamic palaeochannel, suggesting this fluvial axis has had a past role in shaping present-day genetic structure. Moreover, dating analyses show that these splits have been asynchronous, with species responding differently to the riverine barrier. Comparisons informed by the geological history of the Paraná River and its influence on the ecological and climatic differences among ecoregions in the study area further bolster the finding that responses to this geographical barrier have been species-specific.

Genome-wide variation in DNA methylation is associated with stress resilience and plumage brightness in a wild bird.

Individuals often differ in their ability to cope with challenging environmental and social conditions. Evidence from model systems suggests that patterns of DNA methylation are associated with variation in coping ability. These associations could arise directly if methylation plays a role in controlling the physiological response to stressors by, among other things, regulating the release of glucocorticoids in response to challenges. Alternatively, the association could arise indirectly if methylation and resilience have a common cause, such as early-life conditions. In either case, methylation might act as a biomarker for coping ability. At present, however, relatively little is known about whether variation in methylation is associated with organismal performance and resilience under natural conditions. We studied genome-wide patterns of DNA methylation in free-living female tree swallows (Tachycineta bicolor) using methylated DNA immunoprecipitation (MeDIP) and a tree swallow genome that was assembled for this study. We identified areas of the genome that were differentially methylated with respect to social signal expression (breast brightness) and physiological traits (ability to terminate the glucocorticoid stress response through negative feedback). We also asked whether methylation predicted resilience to a subsequent experimentally imposed challenge. Individuals with brighter breast plumage and higher stress resilience had lower methylation at differentially methylated regions across the genome. Thus, widespread differences in methylation predicted both social signal expression and the response to future challenges under natural conditions. These results have implications for predicting individual differences in resilience, and for understanding the mechanistic basis of resilience and its environmental and social mediators.

Contrasting evolutionary histories in Neotropical birds: Divergence across an environmental barrier in South America.

Avian diversity in the Neotropics has been traditionally attributed to the effect of vicariant forces promoting speciation in allopatry. Recent studies have shown that phylogeographical patterns shared among codistributed species cannot be explained by a single vicariant event, as species responses to a common barrier depend on the biological attributes of each taxon. The open vegetation corridor (OVC) isolates Amazonia and the Andean forests from the Atlantic Forest, creating a notorious pattern of avian taxa that are disjunctly codistributed in these forests. Here, we studied and compared the evolutionary histories of Ramphotrigon megacephalum and Pipraeidea melanonota, two passerines with allopatric populations east and west of the OVC that represent different subspecies. These species differ in their biological attributes: R. megacephalum is a sedentary, forest specialist mostly confined to bamboo understorey, whereas P. melanonota is a seasonal migrant and generalist species that ranges in a variety of closed and semi-open environments. We performed genetic and genomic analyses, complemented with the study of coloration and behavioural differentiation, to assess population divergence across the OVC. We found that the evolutionary histories of both R. megacephalum and P. melanonota have been shaped by this environmental barrier. However, these species responded in different and asynchronous manners to the establishment of the OVC and to past connections between the currently isolated South American forests, which can be mostly explained by their distinct ecologies and dispersal abilities. Our results support the fact that the biological attributes of species can make their evolutionary histories idiosyncratic.

Phylogeographic variation within the Buff-browed Foliage-gleaner (Aves: Furnariidae: Syndactyla rufosuperciliata) supports an Andean-Atlantic forests connection via the Cerrado.

We evaluated whether the Andean and the Atlantic forests acted as refugia during the Quaternary, and tested biogeographic hypotheses about the regions involved in the connectivity between those biomes (through the Chaco or the Cerrado). To achieve these goals we selected the Buff-browed Foliage-gleaner Syndactyla rufosuperciliata (Aves, Furnariidae) as a study system, a taxon distributed between the Andean and Atlantic forest. We first explored the historical connectivity between regions through niche modeling. We subsequently used DNA sequences (n = 71 individuals) and genomic analyses (ddRADseq, n = 33 individuals) to evaluate population genetic structure and gene flow within this species. Finally, we performed population model selection using Approximate Bayesian Computation. Our findings indicate that the Andean and the Atlantic forests acted as refugia, and that the populations of the focal species from both regions contacted through the Cerrado region, thus suggesting that the historical dynamics of Andean and Atlantic forests are important for the evolution of forest birds in the region. The results are in agreement with studies of other organisms and may indicate a more general pattern of connectivity among biomes in the Neotropics. Finally, we recommend recognizing both the Andean and the Altantic forests lineages of S. rufosuperciliata as independent species.

Correlated patterns of genetic diversity and differentiation across an avian family.

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole-genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.

Genomic evidence of demographic fluctuations and lack of genetic structure across flyways in a long distance migrant, the European turtle dove.

BACKGROUND: Understanding how past climatic oscillations have affected organismic evolution will help predict the impact that current climate change has on living organisms. The European turtle dove, Streptopelia turtur, is a warm-temperature adapted species and a long distance migrant that uses multiple flyways to move between Europe and Africa. Despite being abundant, it is categorized as vulnerable because of a long-term demographic decline. We studied the demographic history and population genetic structure of the European turtle dove using genomic data and mitochondrial DNA sequences from individuals sampled across Europe, and performing paleoclimatic niche modelling simulations. RESULTS: Overall our data suggest that this species is panmictic across Europe, and is not genetically structured across flyways. We found the genetic signatures of demographic fluctuations, inferring an effective population size (Ne) expansion that occurred between the late Pleistocene and early Holocene, followed by a decrease in the Ne that started between the mid Holocene and the present. Our niche modelling analyses suggest that the variations in the Ne are coincident with recent changes in the availability of suitable habitat. CONCLUSIONS: We argue that the European turtle dove is prone to undergo demographic fluctuations, a trait that makes it sensitive to anthropogenic impacts, especially when its numbers are decreasing. Also, considering the lack of genetic structure, we suggest all populations across Europe are equally relevant for conservation.

Distinguishing noise from signal in patterns of genomic divergence in a highly polymorphic avian radiation.

Recently diverged taxa provide the opportunity to search for the genetic basis of the phenotypes that distinguish them. Genomic scans aim to identify loci that are diverged with respect to an otherwise weakly differentiated genetic background. These loci are candidates for being past targets of selection because they behave differently from the rest of the genome that has either not yet differentiated or that may cross species barriers through introgressive hybridization. Here we use a reduced-representation genomic approach to explore divergence among six species of southern capuchino seedeaters, a group of recently radiated sympatric passerine birds in the genus Sporophila. For the first time in these taxa, we discovered a small proportion of markers that appeared differentiated among species. However, when assessing the significance of these signatures of divergence, we found that similar patterns can also be recovered from random grouping of individuals representing different species. A detailed demographic inference indicates that genetic differences among Sporophila species could be the consequence of neutral processes, which include a very large ancestral effective population size that accentuates the effects of incomplete lineage sorting. As these neutral phenomena can generate genomic scan patterns that mimic those of markers involved in speciation and phenotypic differentiation, they highlight the need for caution when ascertaining and interpreting differentiated markers between species, especially when large numbers of markers are surveyed. Our study provides new insights into the demography of the southern capuchino radiation and proposes controls to distinguish signal from noise in similar genomic scans.

Genomic islands of speciation harbor genes underlying coloration differences in a pair of Neotropical seedeaters.

Incomplete speciation can be leveraged to associate phenotypes with genotypes, thus providing insights into the traits relevant to the reproductive isolation of diverging taxa. We investigate the genetic underpinnings of the phenotypic differences between Sporophila plumbea and Sporophila beltoni. Sporophila beltoni has only recently been described based, most notably, based on differences in bill coloration (yellow vs. black in S. plumbea). Both species are indistinguishable through mtDNA or reduced-representation genomic data, and even whole-genome sequencing revealed low genetic differentiation. Demographic reconstructions attribute this genetic homogeneity to gene flow, despite divergence in the order of millions of generations. We found a narrow hybrid zone in southern Brazil where genetically, yet not phenotypically, admixed individuals appear to be prevalent. Despite the overall low genetic differentiation, we identified 3 narrow peaks along the genome with highly differentiated SNPs. These regions harbor 6 genes, one of which is involved in pigmentation (EDN3) and is a candidate for controlling bill color. Within the outlier peaks, we found signatures of resistance to gene flow, as expected for islands of speciation. Our study shows how genes related to coloration traits are likely involved in generating prezygotic isolation and establishing species boundaries early in speciation.

Population structure and connectivity among coastal and freshwater Kelp Gull (Larus dominicanus) populations from Patagonia.

The genetic identification of evolutionary significant units and information on their connectivity can be used to design effective management and conservation plans for species of concern. Despite having high dispersal capacity, several seabird species show population structure due to both abiotic and biotic barriers to gene flow. The Kelp Gull is the most abundant species of gull in the southern hemisphere. In Argentina it reproduces in both marine and freshwater environments, with more than 100,000 breeding pairs following a metapopulation dynamic across 140 colonies in the Atlantic coast of Patagonia. However, little is known about the demography and connectivity of inland populations. We aim to provide information on the connectivity of the largest freshwater colonies (those from Nahuel Huapi Lake) with the closest Pacific and Atlantic populations to evaluate if these freshwater colonies are receiving immigrants from the larger coastal populations. We sampled three geographic regions (Nahuel Huapi Lake and the Atlantic and Pacific coasts) and employed a reduced-representation genomic approach to genotype individuals for single-nucleotide polymorphisms (SNPs). Using clustering and phylogenetic analyses we found three genetic groups, each corresponding to one of our sampled regions. Individuals from marine environments are more closely related to each other than to those from Nahuel Huapi Lake, indicating that the latter population constitutes the first freshwater Kelp Gull colony to be identified as an evolutionary significant unit in Patagonia.

Continental phylogeography of an ecologically and morphologically diverse Neotropical songbird, Zonotrichia capensis.

BACKGROUND: The Neotropics are exceptionally diverse, containing roughly one third of all extant bird species on Earth. This remarkable species richness is thought to be a consequence of processes associated with both Andean orogenesis throughout the Tertiary, and climatic fluctuations during the Quaternary. Phylogeographic studies allow insights into how such events might have influenced evolutionary trajectories of species and ultimately contribute to a better understanding of speciation. Studies on continentally distributed species are of particular interest because different populations of such taxa may show genetic signatures of events that impacted the continent-wide biota. Here we evaluate the genealogical history of one of the world's most broadly-distributed and polytypic passerines, the rufous-collared sparrow (Zonotrichia capensis). RESULTS: We obtained control region DNA sequences from 92 Zonotrichia capensis individuals sampled across the species' range (Central and South America). Six additional molecular markers, both nuclear and mitochondrial, were sequenced for a subset of individuals with divergent control region haplotypes. Median-joining network analysis, and Bayesian and maximum parsimony phylogenetic analyses all recovered three lineages: one spanning Middle America, the Dominican Republic, and north-western South America; one encompassing the Dominican Republic, Roraima (Venezuela) and La Paz (Bolivia) south to Tierra del Fuego, Argentina; and a third, including eastern Argentina and Brazil. Phylogenetic analyses suggest that the Middle American/north-western South American clade is sister to the remaining two. Bayesian and maximum likelihood coalescent simulations used to study lineage demographic history, diversification times, migration rates and population expansion together suggested that diversification of the three lineages occurred rapidly during the Pleistocene, with negligible gene flow, leaving genetic signatures of population expansions. CONCLUSIONS: The Pleistocene history of the rufous-collared sparrow involved extensive range expansion from a probable Central American origin. Its remarkable morphological and behavioral diversity probably represents recent responses to local conditions overlying deeper patterns of lineage diversity, which are themselves produced by isolation and the history of colonization of South America.

Comparative genomics of two Empidonax flycatchers reveal candidate genes for bird song production.

Whole-genome-level comparisons of sister taxa that vary in phenotype against a background of high genomic similarity can be used to identify the genomic regions that might underlie their phenotypic differences. In wild birds, this exploratory approach has detected markers associated with plumage coloration, beak and wing morphology, and complex behavioral traits like migration. Here, we use genomic comparisons of two closely related suboscine flycatchers (Empidonax difficilis and E. occidentalis) and their hybrids to search for candidate genes underlying their variation in innate vocal signals. We sequenced the genomes of 20 flycatchers that sang one of two species-specific pure song types and 14 putative hybrid individuals with intermediate song types. In the resulting genomic comparisons, we found six areas of high differentiation that may be associated with variation in nonlearned songs. These narrow regions of genomic differentiation contain a total of 67 described genes, of which three have been previously associated with forms of language impairment and dyslexia in humans and 18 are known to be differentially expressed in the song nuclei regions of the avian brain compared with adjacent parts of the avian brain. This "natural experiment" therefore may help identify loci associated with song differences that merit further study across bird lineages with both learned and innate vocalizations.

Rapid phenotypic evolution during incipient speciation in a continental avian radiation.

Adaptive radiations have helped shape how we view animal speciation, particularly classic examples such as Darwin's finches, Hawaiian fruitflies and African Great Lakes cichlids. These 'island' radiations are comparatively recent, making them particularly interesting because the mechanisms that caused diversification are still in motion. Here, we identify a new case of a recent bird radiation within a continentally distributed species group; the capuchino seedeaters comprise 11 Sporophila species originally described on the basis of differences in plumage colour and pattern in adult males. We use molecular data together with analyses of male plumage and vocalizations to understand species limits of the group. We find marked phenotypic variation despite lack of mitochondrial DNA monophyly and few differences in other putatively neutral nuclear markers. This finding is consistent with the group having undergone a recent radiation beginning in the Pleistocene, leaving genetic signatures of incomplete lineage sorting, introgressive hybridization and demographic expansions. We argue that this apparent uncoupling between neutral DNA homogeneity and phenotypic diversity is expected for a recent group within the framework of coalescent theory. Finally, we discuss how the ecology of open habitats in South America during the Pleistocene could have helped promote this unique and ongoing radiation.