Mapping seasonal migration in a songbird hybrid zone -- heritability, genetic correlations, and genomic patterns linked to speciation.

Seasonal migration is a widespread behavior relevant for adaptation and speciation, yet knowledge of its genetic basis is limited. We leveraged advances in tracking and sequencing technologies to bridge this gap in a well-characterized hybrid zone between songbirds that differ in migratory behavior. Migration requires the coordinated action of many traits, including orientation, timing, and wing morphology. We used genetic mapping to show these traits are highly heritable and genetically correlated, explaining how migration has evolved so rapidly in the past and suggesting future responses to climate change may be possible. Many of these traits mapped to the same genomic regions and small structural variants indicating the same, or tightly linked, genes underlie them. Analyses integrating transcriptomic data indicate cholinergic receptors could control multiple traits. Furthermore, analyses integrating genomic differentiation further suggested genes underlying migratory traits help maintain reproductive isolation in this hybrid zone.

Reduced hybrid survival in a migratory divide between songbirds.

Migratory divides, hybrid zones between populations that use different seasonal migration routes, are hypothesised to contribute to speciation. Specifically, relative to parental species, hybrids at divides are predicted to exhibit (1) intermediate migratory behaviour and (2) reduced fitness as a result. We provide the first direct test of the second prediction here with one of the largest existing avian tracking datasets, leveraging a divide between Swainson's thrushes where the first prediction is supported. Using detection rates as a proxy for survival, our results supported the migratory divide hypothesis with lower survival rates for hybrids than parental forms. This finding was juvenile-specific (vs. adults), suggesting selection against hybrids is stronger earlier in life. Reduced hybrid survival was not explained by selection against intermediate phenotypes or negative interactions among phenotypes. Additional work connecting specific features of migration is needed, but these patterns provide strong support for migration as an ecological driver of speciation.

Night-time neuronal activation of Cluster N in a North American songbird.

Night-migrating songbirds utilize the Earth's magnetic field to help navigate to and from their breeding sites each year. A region of the avian forebrain called Cluster N has been shown to be activated during night migratory behavior and it has been implicated in processing geomagnetic information. Previous studies with night-migratory European songbirds have shown that neuronal activity at Cluster N is higher at night than during the day. Comparable work in North American migrants has only been performed in one species of swallows, so extension of examination for Cluster N in other migratory birds is needed. In addition, it is unclear if Cluster N activation is lateralized and the full extent of its boundaries in the forebrain have yet to be described. We used sensory-driven gene expression based on ZENK and the Swainson's thrush, a night-migratory North American songbird, to fill these knowledge gaps. We found elevated levels of gene expression in night- vs. day-active thrushes and no evidence for lateralization in this region. We further examined the anatomical extent of neural activation in the forebrain using 3D reconstruction topology. Our findings demonstrate that Swainson's thrushes possess an extensive bilateral night-activated Cluster N region in the forebrain similar to other European avian species, suggesting that Cluster N is highly conserved in nocturnal migrants.

Gene regulation and speciation in a migratory divide between songbirds.

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration - an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns - no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation.

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

Structural genomic variation and migratory behavior in a wild songbird.

Structural variants (SVs) are a major source of genetic variation; and descriptions in natural populations and connections with phenotypic traits are beginning to accumulate in the literature. We integrated advances in genomic sequencing and animal tracking to begin filling this knowledge gap in the Eurasian blackcap. Specifically, we (a) characterized the genome-wide distribution, frequency, and overall fitness effects of SVs using haplotype-resolved assemblies for 79 birds, and (b) used these SVs to study the genetics of seasonal migration. We detected >15 K SVs. Many SVs overlapped repetitive regions and exhibited evidence of purifying selection suggesting they have overall deleterious effects on fitness. We used estimates of genomic differentiation to identify SVs exhibiting evidence of selection in blackcaps with different migratory strategies. Insertions and deletions dominated the SVs we identified and were associated with genes that are either directly (e.g., regulatory motifs that maintain circadian rhythms) or indirectly (e.g., through immune response) related to migration. We also broke migration down into individual traits (direction, distance, and timing) using existing tracking data and tested if genetic variation at the SVs we identified could account for phenotypic variation at these traits. This was only the case for 1 trait-direction-and 1 specific SV (a deletion on chromosome 27) accounted for much of this variation. Our results highlight the evolutionary importance of SVs in natural populations and provide insight into the genetic basis of seasonal migration.

The Ecology of Hybrid Incompatibilities.

Ecologically mediated selection against hybrids, caused by hybrid phenotypes fitting poorly into available niches, is typically viewed as distinct from selection caused by epistatic Dobzhansky-Muller hybrid incompatibilities. Here, we show how selection against transgressive phenotypes in hybrids manifests as incompatibility. After outlining our logic, we summarize current approaches for studying ecology-based selection on hybrids. We then quantitatively review QTL-mapping studies and find traits differing between parent taxa are typically polygenic. Next, we describe how verbal models of selection on hybrids translate to phenotypic and genetic fitness landscapes, highlighting emerging approaches for detecting polygenic incompatibilities. Finally, in a synthesis of published data, we report that trait transgression-and thus possibly extrinsic hybrid incompatibility in hybrids-escalates with the phenotypic divergence between parents. We discuss conceptual implications and conclude that studying the ecological basis of hybrid incompatibility will facilitate new discoveries about mechanisms of speciation.

Learning and memory in hybrid migratory songbirds: cognition as a reproductive isolating barrier across seasons.

Hybrid zones can be used to identify traits that maintain reproductive isolation and contribute to speciation. Cognitive traits may serve as post-mating reproductive isolating barriers, reducing the fitness of hybrids if, for example, misexpression occurs in hybrids and disrupts important neurological mechanisms. We tested this hypothesis in a hybrid zone between two subspecies of Swainson's thrushes (Catharus ustulatus) using two cognitive tests-an associative learning spatial test and neophobia test. We included comparisons across the sexes and seasons (spring migration and winter), testing if hybrid females performed worse than males (as per Haldane's rule) and if birds (regardless of ancestry or sex) performed better during migration, when they are building navigational maps and encountering new environments. We documented reduced cognitive abilities in hybrids, but this result was limited to males and winter. Hybrid females did not perform worse than males in either season. Although season was a significant predictor of performance, contrary to our prediction, all birds learned faster during the winter. The hypothesis that cognitive traits could serve as post-mating isolating barriers is relatively new; this is one of the first tests in a natural hybrid zone and non-food-caching species. We also provide one of the first comparisons of cognitive abilities between seasons. Future neurostructural and neurophysiological work should be used to examine mechanisms underlying our behavioral observations.

Genomic architecture of migration timing in a long-distance migratory songbird.

The impact of climate change on spring phenology poses risks to migratory birds, as migration timing is controlled predominantly by endogenous mechanisms. Despite recent advances in our understanding of the underlying genetic basis of migration timing, the ways that migration timing phenotypes in wild individuals may map to specific genomic regions requires further investigation. We examined the genetic architecture of migration timing in a long-distance migratory songbird (purple martin, Progne subis subis) by integrating genomic data with an extensive dataset of direct migratory tracks. A moderate to large amount of variance in spring migration arrival timing was explained by genomics (proportion of phenotypic variation explained by genomics = 0.74; polygenic score R2 = 0.24). On chromosome 1, a region that was differentiated between migration timing phenotypes contained genes that could facilitate nocturnal flights and act as epigenetic modifiers. Overall, these results advance our understanding of the genomic underpinnings of migration timing.

The genetics of bird migration.

Bird migration is one of the most amazing biological phenomena exhibited by organisms today, with birds as small as three grams travelling thousands of kilometers twice each year. Most people are familiar with this behaviour; they likely recognize the seasonal movements of birds each year and may have even witnessed specific migratory events (e.g., geese flying in their characteristic V-formations). What many people may not know is that bird migration has a strong genetic basis, especially in groups like songbirds who migrate alone and at night, preventing juveniles from following more experienced adult birds on their first migration. Successful migration also requires the coordinated function of many traits, not only behavioural traits like timing and orientation but also morphological and physiological traits to permit these long distance treks. Together, these traits comprise what has been called the 'migratory syndrome'. In this primer, we provide a broad overview of research on the genetics of migratory traits, from early experimental work to future functional assays and their importance for bird conservation. We focus primarily on songbirds as most work on the genetics of migration has been conducted with this group.

Population-specific association of Clock gene polymorphism with annual cycle timing in stonechats.

Timing is essential for survival and reproduction of organisms across the tree of life. The core circadian clock gene Clk is involved in the regulation of annual timing events and shows highly conserved sequence homology across vertebrates except for one variable region of poly-glutamine repeats. Clk genotype varies in some species with latitude, seasonal timing and migration. However, findings are inconsistent, difficult to disentangle from environmental responses, and biased towards high latitudes. Here we combine field data with a common-garden experiment to study associations of Clk polymorphism with latitude, migration and annual-cycle timing within the stonechat species complex across its trans-equatorial distribution range. Our dataset includes 950 records from 717 individuals from nine populations with diverse migratory strategies. Gene diversity was lowest in resident African and Canary Island populations and increased with latitude, independently of migration distance. Repeat length and annual-cycle timing was linked in a population-specific way. Specifically, equatorial African stonechats showed delayed timing with longer repeat length for all annual-cycle stages. Our data suggest that at low latitudes with nearly constant photoperiod, Clk genotype might orchestrate a range of consistent, individual chronotypes. In contrast, the influence of Clk on annual-cycle timing at higher latitudes might be mediated by its interactions with genes involved in (circadian) photoperiodic pathways.

Estimating hybridization rates in the wild: Easier said than done?

Hybridization has important effects on the evolutionary trajectories of natural populations but estimates of this process in the wild and at the individual-level are lacking. Justyn et al. attempted to fill this gap using the citizen science database eBird but there are limitations to this approach. Here, we outline and directly test these limitations using literature searches, case studies, and a comparison between eBird and Birds of North America (BNA), a database that documents hybridization using the scientific literature. We use a hybrid zone between Lazuli and Indigo buntings to highlight the importance of considering geographic range when estimating rates of hybridization and two literature searches to show the importance of considering cryptic hybrids (those that cannot be identified using phenotypic traits) when quantifying these rates. We also use BNA and a case study of hybrid White-faced and Glossy Ibises to show that citizen scientists are underreporting hybrids compared with experts. Justyn et al. highlighted an important gap in the literature, but their results likely represent the lower limit of hybridization between birds and a more nuanced interpretation of their results (e.g., considering extrinsic postzygotic selection) is needed.

Individual variability and versatility in an eco-evolutionary model of avian migration.

Seasonal migration is a complex and variable behaviour with the potential to promote reproductive isolation. In Eurasian blackcaps (Sylvia atricapilla), a migratory divide in central Europe separating populations with southwest (SW) and southeast (SE) autumn routes may facilitate isolation, and individuals using new wintering areas in Britain show divergence from Mediterranean winterers. We tracked 100 blackcaps in the wild to characterize these strategies. Blackcaps to the west and east of the divide used predominantly SW and SE directions, respectively, but close to the contact zone many individuals took intermediate (S) routes. At 14.0° E, we documented a sharp transition from SW to SE migratory directions across only 27 (10-86) km, implying a strong selection gradient across the divide. Blackcaps wintering in Britain took northwesterly migration routes from continental European breeding grounds. They originated from a surprisingly extensive area, spanning 2000 km of the breeding range. British winterers bred in sympatry with SW-bound migrants but arrived 9.8 days earlier on the breeding grounds, suggesting some potential for assortative mating by timing. Overall, our data reveal complex variation in songbird migration and suggest that selection can maintain variation in migration direction across short distances while enabling the spread of a novel strategy across a wide range.

The evolutionary history and genomics of European blackcap migration.

Seasonal migration is a taxonomically widespread behaviour that integrates across many traits. The European blackcap exhibits enormous variation in migration and is renowned for research on its evolution and genetic basis. We assembled a reference genome for blackcaps and obtained whole genome resequencing data from individuals across its breeding range. Analyses of population structure and demography suggested divergence began ~30,000 ya, with evidence for one admixture event between migrant and resident continent birds ~5000 ya. The propensity to migrate, orientation and distance of migration all map to a small number of genomic regions that do not overlap with results from other species, suggesting that there are multiple ways to generate variation in migration. Strongly associated single nucleotide polymorphisms (SNPs) were located in regulatory regions of candidate genes that may serve as major regulators of the migratory syndrome. Evidence for selection on shared variation was documented, providing a mechanism by which rapid changes may evolve.

Every year as the seasons change, thousands of animals migrate huge distances in search of food or better climates. As far as migrations go, there might be none so impressive as the trans-oceanic flights made by small migrating songbirds. These birds can weigh as little as three grams and travel up to 15,000 kilometres. Most migrate alone and at night and yet still manage to return to the same location each year. Several strands of research suggest there could be a genetic basis to their migratory behaviour, but exactly which genes control this phenomenon remains poorly understood. One small songbird that has been studied for decades is the European blackcap. This species exhibits a real variety of migration patterns. Some blackcaps travel rather short distances, others much further, and some populations do not migrate at all. Populations that share the same breeding grounds in the summer may migrate in different directions in the autumn. These features make it a good species to study the genetic variation between populations that migrate in different directions and over different distances. However, only in recent years has advancing technology made it possible to comprehensively study an animal's entire genome, leaving no gene unturned. Now, Delmore et al. have used high-throughput sequencing technologies to trace the evolutionary history of migration in European blackcap and started by assembling a reference genome for the species. Then, the genomes of 110 blackcaps from several populations that take different annual migrations were compared to the reference. This revealed that the populations began to diverge some 30,000 years ago and that there was some apparent gene mixing between groups of migrating and resident blackcaps around 5,000 years ago. The analysis showed only a small set of genes code for their differences in migration. Additionally, while the candidate genes were shown to be common among blackcaps, the genes identified did not match those reported from studies of other migrating songbirds. Finally, Delmore et al. also noted that the differences between the populations tend to be in the parts of the genome that control whether a given gene is switched on or off, which could explain how new migratory behaviours can rapidly evolve. This study is one of the most comprehensive genomic analysis of migration to date. It is important work as songbirds, like other animals, are responding to increasing pressures of environmental and climate change. In time, the findings could be used to support conservation efforts whereby genetic analyses could determine if certain populations possess enough variation to respond to coming changes in their habitats.

Shared Patterns of Genome-Wide Differentiation Are More Strongly Predicted by Geography Than by Ecology.

Closely related populations often display similar patterns of genomic differentiation, yet it remains an open question which ecological and evolutionary forces generate these patterns. The leading hypothesis is that this similarity in divergence is driven by parallel natural selection. However, several recent studies have suggested that these patterns may instead be a product of the depletion of genetic variation that occurs as result of background selection (i.e., linked negative selection). To date, there have been few direct tests of these competing hypotheses. To determine the relative contributions of background selection and parallel selection to patterns of repeated differentiation, we examined 24 independently derived populations of freshwater stickleback occupying a variety of niches and estimated genomic patterns of differentiation in each relative to their common marine ancestor. Patterns of genetic differentiation were strongly correlated across pairs of freshwater populations adapting to the same ecological niche, supporting a role for parallel natural selection. In contrast to other recent work, our study comparing populations adapting to the same niche produced no evidence signifying that similar patterns of genomic differentiation are generated by background selection. We also found that overall patterns of genetic differentiation were considerably more similar for populations found in closer geographic proximity. In fact, the effect of geography on the repeatability of differentiation was greater than that of parallel selection. Our results suggest that shared selective landscapes and ancestral variation are the key drivers of repeated patterns of differentiation in systems that have recently colonized novel environments.

Avian Population Genomics Taking Off: Latest Findings and Future Prospects.

Birds are one of the most recognizable and diverse groups of organisms on earth. This group has played an important role in many fields, including the development of methods in behavioral ecology and evolutionary theory. The use of population genomics took off following the advent of high-throughput sequencing in various taxa. Several features of avian genomes make them particularly amenable for work in this field, including their nucleated red blood cells permitting easy DNA extraction and small, compact genomes. We review the latest findings in the population genomics of birds here, emphasizing questions related to behavior, ecology, evolution, and conservation. Additionally, we include insights in trait mapping and the ability to obtain accurate estimates of important summary statistics for conservation (e.g., genetic diversity and inbreeding). We highlight roadblocks that will need to be overcome in order to advance work on the population genomics of birds and prospects for future work. Roadblocks include the assembly of more contiguous reference genomes using long-reads and optical mapping. Prospects include the integration of population genomics with additional fields (e.g., landscape genetics, phylogeography, and genomic mapping) along with studies beyond genetic variants (e.g., epigenetics).

Cross-decades stability of an avian hybrid zone.

Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long-term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S. occidentalis) in the Washington Cascades was previously inferred to be moving from northern S. townsendi southwards towards S. occidentalis, based on plumage and behavioural patterns as well as a 2000-km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2-3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987-94 and 2015-16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a 'tension zone' model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2-3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species.

Comparative analysis examining patterns of genomic differentiation across multiple episodes of population divergence in birds.

Heterogeneous patterns of genomic differentiation are commonly documented between closely related populations and there is considerable interest in identifying factors that contribute to their formation. These factors could include genomic features (e.g., areas of low recombination) that promote processes like linked selection (positive or purifying selection that affects linked neutral sites) at specific genomic regions. Examinations of repeatable patterns of differentiation across population pairs can provide insight into the role of these factors. Birds are well suited for this work, as genome structure is conserved across this group. Accordingly, we reestimated relative (FST ) and absolute (dXY ) differentiation between eight sister pairs of birds that span a broad taxonomic range using a common pipeline. Across pairs, there were modest but significant correlations in window-based estimates of differentiation (up to 3% of variation explained for FST and 26% for dXY ), supporting a role for processes at conserved genomic features in generating heterogeneous patterns of differentiation; processes specific to each episode of population divergence likely explain the remaining variation. The role genomic features play was reinforced by linear models identifying several genomic variables (e.g., gene densities) as significant predictors of FST and dXY repeatability. FST repeatability was higher among pairs that were further along the speciation continuum (i.e., more reproductively isolated) providing further insight into how genomic differentiation changes with population divergence; early stages of speciation may be dominated by positive selection that is different between pairs but becomes integrated with processes acting according to shared genomic features as speciation proceeds.

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.

Gene flow and selection interact to promote adaptive divergence in regions of low recombination.

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (FST and dXY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.