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.

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.

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.

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.

Candidate genes for migration do not distinguish migratory and non-migratory birds.

Migratory traits in birds have been shown to have a strong heritable component and several candidate genes have been suggested to control these migratory traits. To investigate if the genetic makeup of one or a set of these candidate genes can be used to identify a general pattern between migratory and non-migratory birds, we extracted genomic sequence data for 25 hypothesised candidate genes for migration from 70 available genomes across all orders of Aves and characterised sequence divergence between migratory and non-migratory phenotypes. When examining each gene separately across all species, we did not identify any genetic variants in candidate genes that distinguished migrants from non-migrants; any resulting pattern was driven by the phylogenetic signal. This was true for each gene analysed independently, but also for concatenated sequence alignments of all candidate genes combined. We also attempted to distinguish between migrant and non-migrants using structural features at four candidate genes that have previously been reported to show associated with migratory behaviour but did not pick up a signal for migratory phenotype here either. Finally, a screen for dN/dS ratio across all focal candidate genes to probe for putative features of selection did not uncover a pattern, though this might not be expected given the broad phylogenetic scale used here. Our study demonstrates the potential of public genomic data to test for general patterns of migratory gene candidates in a cross-species comparative context, and raise questions on the applicability of candidate gene approaches in a macro-evolutionary context to understand the genetic architecture of migratory behaviour.

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species.

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within-group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within-population selection-induced reductions in variation. An important implication of this 'sweep-before-differentiation' model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.

Phenotypic divergence during speciation is inversely associated with differences in seasonal migration.

Differences in seasonal migration might promote reproductive isolation and differentiation by causing populations in migratory divides to arrive on the breeding grounds at different times and/or produce hybrids that take inferior migratory routes. We examined this question by quantifying divergence in song, colour, and morphology between sister pairs of North American migratory birds. We predicted that apparent rates of phenotypic differentiation would differ between pairs that do and do not form migratory divides. Consistent with this prediction, results from mixed effects models and Ornstein-Uhlenbeck models of evolution showed different rates of divergence between these groups; surprisingly, differentiation was greater among non-divide pairs. We interpret this finding as a result of variable rates of population blending and fusion between partially diverged forms. Ancient pairs of populations that subsequently fused are now observed as a single form, whereas those that did not fuse are observable as pairs and included in our study. We propose that fusion of two populations is more likely to occur when they have similar migratory routes and little other phenotypic differentiation that would cause reproductive isolation. By contrast, pairs with migratory divides are more likely to remain reproductively isolated, even when differing little in other phenotypic traits. These findings suggest that migratory differences may be one among several isolating barriers that prevent divergent populations from fusing and thereby increase the likelihood that they will continue differentiating as distinct species.

Genomic analysis of a migratory divide reveals candidate genes for migration and implicates selective sweeps in generating islands of differentiation.

Differential gene flow, reductions in diversity following linked selection and/or features of the genome can structure patterns of genomic differentiation during the process of speciation. Possible sources of reproductive isolation are well studied between coastal and inland subspecies groups of Swainson's thrushes, with differences in seasonal migratory behaviour likely playing a key role in reducing hybrid fitness. We assembled and annotated a draft reference genome for this species and generated whole-genome shotgun sequence data for populations adjacent to the hybrid zone between these groups. We documented substantial genomewide heterogeneity in relative estimates of genetic differentiation between the groups. Within population diversity was lower in areas of high relative differentiation, supporting a role for selective sweeps in generating this pattern. Absolute genetic differentiation was reduced in these areas, further suggesting that recurrent selective sweeps in the ancestral population and/or between divergent populations following secondary contact likely occurred. Relative genetic differentiation was also higher near centromeres and on the Z chromosome, suggesting that features of the genome also contribute to genomewide heterogeneity. Genes linked to migratory traits were concentrated in islands of differentiation, supporting previous suggestions that seasonal migration is under divergent selection between Swainson's thrushes. Differences in migratory behaviour likely play a central role in the speciation of many taxa; we developed the infrastructure here to permit future investigations into the role several candidate genes play in reducing gene flow between not only Swainson's thrushes but other species as well.

Hybrid songbirds employ intermediate routes in a migratory divide.

Migratory divides are contact zones between populations that use different routes to navigate around unsuitable areas on seasonal migration. Hybrids in divides have been predicted to employ intermediate and potentially inferior routes. We provide the first direct test of this hypothesis, using light-level geolocators to track birds breeding in a hybrid zone between Swainson's thrushes in western Canada. Compared to parental forms, hybrids exhibited increased variability in their migratory routes, with some using intermediate routes that crossed arid and mountainous regions, and some using the same routes as one parental group on fall migration and the other on spring migration. Hybrids also tended to use geographically intermediate wintering sites. Analysis of genetic variation across the hybrid zone suggests moderately strong selection against hybrids. These results indicate that seasonal migratory behaviour might be a source of selection against hybrids, supporting a possible role for migration in 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 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.

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.

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.

Dramatic intraspecific differences in migratory routes, stopover sites and wintering areas, revealed using light-level geolocators.

Migratory divides are contact zones between breeding populations that use divergent migratory routes and have been described in a variety of species. These divides are of major importance to evolution, ecology and conservation but have been identified using limited band recovery data and/or indirect methods. Data from band recoveries and mitochondrial haplotypes suggested that inland and coastal Swainson's thrushes (Catharus ustulatus) form a migratory divide in western North America. We attached light-level geolocators to birds at the edges of this contact zone to provide, to our knowledge, the first direct test of a putative divide using data from individual birds over the entire annual cycle. Coastal thrushes migrated along the west coast to Mexico, Guatemala and Honduras. Some of these birds used multiple wintering sites. Inland thrushes migrated across the Rocky Mountains, through central North America to Columbia and Venezuela. These birds migrated longer distances than coastal birds and performed a loop migration, navigating over the Gulf of Mexico in autumn and around this barrier in spring. These findings support the suggestion that divergent migratory behaviour could contribute to reproductive isolation between migrants, advance our understanding of their non-breeding ecology, and are integral to development of detailed conservation strategies for this group.

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.