GENETIC, MORPHOLOGICAL, AND ECOLOGICAL CHARACTERIZATION OF A HYBRID ZONE THAT SPANS A MIGRATORY DIVIDE

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

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.

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.

Multi-isotope (δ2H, δ13C, δ15N) feather profiles and morphometrics inform patterns of migratory connectivity in three species of North American swallows

Aerial insectivorous birds have suffered steep population declines in North America over the last 60 years. A lack of information on migratory connectivity between breeding and non-breeding grounds for these species limits our ability to interpret factors affecting their population-specific trends. We determined likely Latin American non-breeding regions of Bank (Riparia riparia), Barn (Hirundo rustica) and Cliff (Petrochelidon pyrrhonota) swallow from populations across their breeding ranges. We used predicted feather hydrogen (δ²H_f) and carbon (δ¹³C_f) isoscapes for winter-grown feathers to indicate areas of highest probability of moult origin and incorporated these results into a cluster analysis to determine likely broad non-breeding regions. We also assessed variation in wing length among populations to determine the potential for this metric to differentiate population moult origins. We then investigated patterns of multi-isotopic (δ²H_f, δ¹³C_f, δ¹⁵N_f) and wing-length niche occupancy by quantifying niche size and overlap among populations under the assumption that broad niches were consistent with low within-species migratory connectivity and narrow and non-overlapping niches with higher connectivity. Multivariate assignment identified different non-breeding regions and potential clusters of moult origin generally corresponding to Central America and northern South America, eastern and south-central South America, and the western and southern part of that continent, with variation within and among populations and species. Separate niche space indicated different wintering habitat or areas used by species or populations whereas niche overlap indicated only potential spatial similarity. Wing length varied significantly among populations by species, being longer in the west and north for Bank and Cliff Swallow and longer in eastern Canadian Barn Swallow populations. Barn Swallow occupied consistently larger isotopic and wing length niche space than the other species. Comparisons among populations across species showed variable isotopic and wing-length niche overlap generally being greater within breeding regions and lower between western and eastern breeding populations supporting a general North American continental divide for all species with generally low migratory connectivity for all species. We present a novel approach to assessing connectivity using inexpensive and broad isotopic approaches that provides the basis for hypothesis testing using more spatially explicit expensive techniques.

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.

Genetic and ecological drivers of molt in a migratory bird

The ability of animals to sync the timing and location of molting (the replacement of hair, skin, exoskeletons or feathers) with peaks in resource availability has important implications for their ecology and evolution. In migratory birds, the timing and location of pre-migratory feather molting, a period when feathers are shed and replaced with newer, more aerodynamic feathers, can vary within and between species. While hypotheses to explain the evolution of intraspecific variation in the timing and location of molt have been proposed, little is known about the genetic basis of this trait or the specific environmental drivers that may result in natural selection for distinct molting phenotypes. Here we take advantage of intraspecific variation in the timing and location of molt in the iconic songbird, the Painted Bunting (Passerina ciris) to investigate the genetic and ecological drivers of distinct molting phenotypes. Specifically, we use genome-wide genetic sequencing in combination with stable isotope analysis to determine population genetic structure and molting phenotype across thirteen breeding sites. We then use genome-wide association analysis (GWAS) to identify a suite of genes associated with molting and pair this with gene-environment association analysis (GEA) to investigate potential environmental drivers of genetic variation in this trait. Associations between genetic variation in molt-linked genes and the environment are further tested via targeted SNP genotyping in 25 additional breeding populations across the range. Together, our integrative analysis suggests that molting is in part regulated by genes linked to feather development and structure (GLI2 and CSPG4) and that genetic variation in these genes is associated with seasonal variation in precipitation and aridity. Overall, this work provides important insights into the genetic basis and potential selective forces behind phenotypic variation in what is arguably one of the most important fitness-linked traits in a migratory bird.

As the Goose Flies: Migration Routes and Timing Influence Patterns of Genetic Diversity in a Circumpolar Migratory Herbivore

Migration schedules and the timing of other annual events (e.g., pair formation and molt) can affect the distribution of genetic diversity as much as where these events occur. The greater white-fronted goose (Anser albifrons) is a circumpolar goose species, exhibiting temporal and spatial variation of events among populations during the annual cycle. Previous range-wide genetic assessments of the nuclear genome based on eight microsatellite loci suggest a single, largely panmictic population despite up to five subspecies currently recognized based on phenotypic differences. We used double digest restriction-site associated DNA (ddRAD-seq) and mitochondrial DNA (mtDNA) sequence data to re-evaluate estimates of spatial genomic structure and to characterize how past and present processes have shaped the patterns of genetic diversity and connectivity across the Arctic and subarctic. We uncovered previously undetected inter-population differentiation with genetic clusters corresponding to sampling locales associated with current management groups. We further observed subtle genetic clustering within each management unit that can be at least partially explained by the timing and directionality of migration events along with other behaviors during the annual cycle. The Tule Goose (A. a. elgasi) and Greenland subspecies (A. a. flavirostris) showed the highest level of divergence among all sampling locales investigated. The recovery of previously undetected broad and fine-scale spatial structure suggests that the strong cultural transmission of migratory behavior restricts gene flow across portions of the species’ range. Our data further highlight the importance of re-evaluating previous assessments conducted based on a small number of highly variable genetic markers in phenotypically diverse species.

A migratory divide spanning two continents is associated with genomic and ecological divergence

Migratory divides are contact zones between breeding populations with divergent migratory strategies during the nonbreeding season. These locations provide an opportunity to evaluate the role of seasonal migration in the maintenance of reproductive isolation, particularly the relationship between population structure and features associated with distinct migratory strategies. We combine light-level geolocators, genomic sequencing, and stable isotopes to investigate the timing of migration and migratory routes of individuals breeding on either side of a migratory divide coinciding with genomic differentiation across a hybrid zone between barn swallow (Hirundo rustica) subspecies in China. Individuals west of the hybrid zone, with H. r. rustica ancestry, had comparatively enriched stable-carbon and hydrogen isotope values and overwintered in eastern Africa, whereas birds east of the hybrid zone, with H. r. gutturalis ancestry, had depleted isotope values and migrated to southern India. The two subspecies took divergent migratory routes around the high-altitude Karakoram Range and arrived on the breeding grounds over 3 weeks apart. These results indicate that assortative mating by timing of arrival and/or selection against hybrids with intermediate migratory traits may maintain reproductive isolation between the subspecies, and that inhospitable geographic features may have contributed to the diversification of Asian avifauna by influencing migratory patterns.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

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.

Autumn migration direction of juvenile willow warblers (Phylloscopus t. trochilus and P. t. acredula) and their hybrids assessed by qPCR SNP genotyping

BACKGROUNDS

Geographic regions, where two closely related taxa with different migration routes come into contact, are known as migratory divides. Hybrids originating from migratory divides are hypothesized to migrate intermediately relative to the parental populations. Few studies have tested this hypothesis in wild birds, and only in hybrids that have completed the migration back to the breeding grounds. Here, we make use of the well-established migration routes of willow warblers (Phylloscopus trochilus), for which the subspecies trochilus and acredula have migration-associated genetic markers on chromosomes 1 and 5. The genetic approach enabled us to analyze the geographic distribution of juveniles during their first autumn migration, predicting that hybrids should be more frequent in the central flyway over Italy than along the typical SW routes of trochilus and SE routes of acredula.

METHODS

Blood and feather samples were collected from wintering birds in Africa (n = 69), and from juveniles during autumn migration in Portugal (n = 33), Italy (n = 38) and Bulgaria (n = 32). Genotyping was carried out by qPCR SNP assays, on one SNP each on chromosome 1 (SNP 65) and chromosome 5 (SNP 285). Both these SNPs have alternative alleles that are highly fixed (> 97%) in each of the subspecies.

RESULTS

The observed combined genotypes of the two SNPs were associated with the known migration routes and wintering distributions of trochilus and acredula, respectively. We found hybrids (HH) among the juveniles in Italy (5/38) and in Portugal (2/33). The proportion of hybrids in Italy was significantly higher than expected from a background rate of hybrid genotypes (1.5%) in allopatric populations of the subspecies.

CONCLUSIONS

Our genetic approach to assign individuals to subspecies and hybrids allowed us to investigate migration direction in juvenile birds on their first migration, which should better reflect the innate migratory direction than studies restricted to successful migrants. The excess of hybrids in Italy, suggests that they employ an intermediate route relative to the parental populations. Our qPCR SNP genotyping method is efficient for processing large sample sizes, and will therefore be useful in migration research of species with known population genetic structure.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

Pleistocene glacial cycles and physical barriers influence phylogeographic structure in Black-capped Chickadees (Poecile atricapillus), a widespread North American passerine

The nonmigratory Black-capped Chickadee (Poecile atricapillus (Linnaeus, 1766)) has a continent-wide distribution extending across large parts of North America. To investigate the phylogeographic structure and verify possible refugia during the last glacial maximum, we sequenced a 678 bp region of the mitochondrial control region from 633 Black-capped Chickadees at 35 sites across North America and performed paleoecological distribution modeling. Two genetically distinct groups were found using multiple analyses: one in Newfoundland (Canada) and a widespread continental group, with additional substructure evident in western continental populations. While gene flow is low throughout the range, it is especially low in peripheral populations. The Newfoundland population has remained isolated from continental populations for at least 65 000 years and contains a number of fixed nucleotide differences. Within the continental populations, Black-capped Chickadees are subdivided into Pacific Coast, Alaska (USA), southeast Rockies, and main-northeast groups consistent with late Pleistocene vicariance events. Evidence of secondary contact was identified between Pacific and main-northeast populations in northwest British Columbia (Canada) and between southeast Rockies and main-northeast groups in Montana (USA). Paleoecological distribution modeling predicted suitable habitat in Alaska, off the coast of Newfoundland, and several locations across the southern United States during the last glacial maximum, whereas suitable habitat during the last interglacial was more similar to the contemporary distribution.

A long winter for the Red Queen: rethinking the evolution of seasonal migration

This paper advances an hypothesis that the primary adaptive driver of seasonal migration is maintenance of site fidelity to familiar breeding locations. We argue that seasonal migration is therefore principally an adaptation for geographic persistence when confronted with seasonality - analogous to hibernation, freeze tolerance, or other organismal adaptations to cyclically fluctuating environments. These ideas stand in contrast to traditional views that bird migration evolved as an adaptive dispersal strategy for exploiting new breeding areas and avoiding competitors. Our synthesis is supported by a large body of research on avian breeding biology that demonstrates the reproductive benefits of breeding-site fidelity. Conceptualizing migration as an adaptation for persistence places new emphasis on understanding the evolutionary trade-offs between migratory behaviour and other adaptations to fluctuating environments both within and across species. Seasonality-induced departures from breeding areas, coupled with the reproductive benefits of maintaining breeding-site fidelity, also provide a mechanism for explaining the evolution of migration that is agnostic to the geographic origin of migratory lineages (i.e. temperate or tropical). Thus, our framework reconciles much of the conflict in previous research on the historical biogeography of migratory species. Although migratory behaviour and geographic range change fluidly and rapidly in many populations, we argue that the loss of plasticity for migration via canalization is an overlooked aspect of the evolutionary dynamics of migration and helps explain the idiosyncratic distributions and migratory routes of long-distance migrants. Our synthesis, which revolves around the insight that migratory organisms travel long distances simply to stay in the same place, provides a necessary evolutionary context for understanding historical biogeographic patterns in migratory lineages as well as the ecological dynamics of migratory connectivity between breeding and non-breeding locations.

Genomic islands of divergence or opportunities for introgression?

In animals, introgression between species is often perceived as the breakdown of reproductive isolating mechanisms, but gene flow between incipient species can also represent a source for potentially beneficial alleles. Recently, genome-wide datasets have revealed clusters of differentiated loci ('genomic islands of divergence') that are thought to play a role in reproductive isolation and therefore have reduced gene flow. We use simulations to further examine the evolutionary forces that shape and maintain genomic islands of divergence between two subspecies of the migratory songbird, Swainson's thrush (Catharus ustulatus), which have come into secondary contact since the last glacial maximum. We find that, contrary to expectation, gene flow is high within islands and is highly asymmetric. In addition, patterns of nucleotide diversity at highly differentiated loci suggest selection was more frequent in a single ecotype. We propose a mechanism whereby beneficial alleles spread via selective sweeps following a post-glacial demographic expansion in one subspecies and move preferentially across the hybrid zone. We find no evidence that genomic islands are the result of divergent selection or reproductive isolation, rather our results suggest that differentiated loci both within and outside islands could provide opportunities for adaptive introgression across porous species boundaries.

Origin and cross-century dynamics of an avian hybrid zone

Background

Characterizations of the dynamics of hybrid zones in space and time can give insights about traits and processes important in population divergence and speciation. We characterized a hybrid zone between tanagers in the genus Ramphocelus (Aves, Thraupidae) located in southwestern Colombia. We evaluated whether this hybrid zone originated as a result of secondary contact or of primary differentiation, and described its dynamics across time using spatial analyses of molecular, morphological, and coloration data in combination with paleodistribution modeling.

Results

Models of potential historical distributions based on climatic data and genetic signatures of demographic expansion suggested that the hybrid zone likely originated following secondary contact between populations that expanded their ranges out of isolated areas in the Quaternary. Concordant patterns of variation in phenotypic characters across the hybrid zone and its narrow extent are suggestive of a tension zone, maintained by a balance between dispersal and selection against hybrids. Estimates of phenotypic cline parameters obtained using specimens collected over nearly a century revealed that, in recent decades, the zone appears to have moved to the east and to higher elevations, and may have become narrower. Genetic variation was not clearly structured along the hybrid zone, but comparisons between historical and contemporary specimens suggested that temporal changes in its genetic makeup may also have occurred.

Conclusions

Our data suggest that the hybrid zone likey resulted from secondary contact between populations. The observed changes in the hybrid zone may be a result of sexual selection, asymmetric gene flow, or environmental change.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12862-017-1096-7) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome sequences of the northern spotted owl (Strix occidentalis caurina) and the barred owl (Strix varia; Aves: Strigiformes: Strigidae) confirm the presence of a duplicated control region

We report here the successful assembly of the complete mitochondrial genomes of the northern spotted owl (Strix occidentalis caurina) and the barred owl (S. varia). We utilized sequence data from two sequencing methodologies, Illumina paired-end sequence data with insert lengths ranging from approximately 250 nucleotides (nt) to 9,600 nt and read lengths from 100–375 nt and Sanger-derived sequences. We employed multiple assemblers and alignment methods to generate the final assemblies. The circular genomes of S. o. caurina and S. varia are comprised of 19,948 nt and 18,975 nt, respectively. Both code for two rRNAs, twenty-two tRNAs, and thirteen polypeptides. They both have duplicated control region sequences with complex repeat structures. We were not able to assemble the control regions solely using Illumina paired-end sequence data. By fully spanning the control regions, Sanger-derived sequences enabled accurate and complete assembly of these mitochondrial genomes. These are the first complete mitochondrial genome sequences of owls (Aves: Strigiformes) possessing duplicated control regions. We searched the nuclear genome of S. o. caurina for copies of mitochondrial genes and found at least nine separate stretches of nuclear copies of gene sequences originating in the mitochondrial genome (Numts). The Numts ranged from 226–19,522 nt in length and included copies of all mitochondrial genes except tRNA^Pro, ND6, and tRNA^Glu. Strix occidentalis caurina and S. varia exhibited an average of 10.74% (8.68% uncorrected p-distance) divergence across the non-tRNA mitochondrial genes.

Linking the wintering and breeding grounds of warblers along the Pacific Flyway

Long‐distance migration is a behavior that is exhibited by many animal groups. The evolution of novel migration routes can play an important role in range expansions, ecological interactions, and speciation. New migration routes may evolve in response to selection in favor of reducing distance between breeding and wintering areas, or avoiding navigational barriers. Many migratory changes are likely to evolve gradually and are therefore difficult to study. Here, we attempt to connect breeding and wintering populations of myrtle warblers (Setophaga coronata coronata) to better understand the possible evolution of distinct migration routes within this species. Myrtle warblers, unlike most other warblers with breeding ranges primarily in eastern North America, have two disjunct overwintering concentrations—one in the southeastern USA and one along the Pacific Coast—and presumably distinct routes to‐and‐from these locations. We studied both myrtle and Audubon's warblers (S. c. auduboni) captured during their spring migration along the Pacific Coast, south of the narrow region where these two taxa hybridize. Using stable hydrogen isotopes and biometric data, we show that those myrtle warblers wintering along the southern Pacific Coast of North America are likely to breed at high latitudes in Alaska and the Yukon rather than in Alberta or further east. Our interpretation is that the evolution of this wintering range and migration route along the Pacific Coast may have facilitated the breeding expansion of myrtle warblers into northwestern North America. Moreover, these data suggest that there may be a migratory divide within genetically similar populations of myrtle warblers.

Migratory orientation in a narrow avian hybrid zone

BACKGROUND

Zones of contact between closely related taxa with divergent migratory routes, termed migratory divides, have been suggested as areas where hybrid offspring may have intermediate and inferior migratory routes, resulting in low fitness of hybrids and thereby promoting speciation. In the Rocky Mountains of Canada there is a narrow hybrid zone between Audubon's and myrtle warblers that is likely maintained by selection against hybrids. Band recoveries and isotopic studies indicate that this hybrid zone broadly corresponds to the location of a possible migratory divide, with Audubon's warblers migrating south-southwest and myrtle warblers migrating southeast. We tested a key prediction of the migratory divide hypothesis: that genetic background would be predictive of migratory orientation among warblers in the center of the hybrid zone.

METHODS

We recorded fall migratory orientation of wild-caught migrating warblers in the center of the hybrid zone as measured by video-based monitoring of migratory restlessness in circular orientation chambers. We then tested whether there was a relationship between migratory orientation and genetic background, as measured using a set of species-specific diagnostic genetic markers.

RESULTS

We did not detect a significant association between orientation and genetic background. There was large variation among individuals in orientation direction. Mean orientation was towards the NE, surprising for birds on fall migration, but aligned with the mountain valley in which the study took place.

CONCLUSIONS

Only one other study has directly analyzed migratory orientation among naturally-produced hybrids in a migratory divide. While the other study showed an association between genetic background and orientation, we did not observe such an association in yellow-rumped warblers. We discuss possible reasons, including the possibility of a lack of a strong migratory divide in this hybrid zone and/or methodological limitations that may have prevented accurate measurements of long-distance migratory orientation.

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.

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.

Wind patterns as a potential driver in the evolution and maintenance of a North American migratory suture zone

Suture zones are areas where range contact zones and hybrid zones of multiple taxa are clustered. Migratory divides, contact zones between divergent populations that breed adjacent to one another but use different migratory routes, are a particular case of suture zones. Although multiple hypotheses for both the formation and maintenance of migratory divides have been suggested, quantitative tests are scarce. Here, we tested whether a novel factor, prevailing winds, was sufficient to explain both the evolution and maintenance of the Cordilleran migratory divide using individual-based models. Empirical observations of eastern birds suggest a circuitous migratory route across Canada before heading south. Western breeders, however, travel south along the Pacific coast to their wintering grounds. We modeled the effect of wind on bird migratory flights by allowing them to float at elevation using spatially explicit modeled wind data. Modeled eastern birds had easterly mean trajectories, whereas western breeders showed significantly more southern trajectories. We also determined that a mean airspeed of 18.5 m s(-1) would be necessary to eliminate this difference in trajectory, a speed that is achieved by waterfowl and shorebirds, but is faster than songbird flight speeds. These results lend support for the potential importance of wind in shaping the phylogeographic history of North American songbirds.

Ecological differentiation and habitat unsuitability maintaining a ground beetle hybrid zone

Exogenous selection via interactions between organisms and environments may influence the dynamics of hybrid zones between species in multiple ways. Two major models of a hybrid zone allowed us to hypothesize that environmental conditions influence hybrid zone dynamics in two ways. In the first model, an environmental gradient determines the mosaic distribution at the boundary between ecologically differentiated species (mosaic hybrid zone model). In the second model, a patch of unsuitable habitat traps a hybrid zone between species whose hybrids are unfit (tension zone model). To test these, we examined the environmental factors influencing the spatial structure of a hybrid zone between the ground beetles Carabus maiyasanus and C. iwawakianus using GIS‐based quantification of environmental factors and a statistical comparison of species distribution models (SDMs). We determined that both of the hypothetical processes can be important in the hybrid zone. We detected interspecific differences in the environmental factors in presence localities and their relative contribution in SDMs. SDMs were not identical between species even within contact areas, but tended to be similar within the range of each species. These results suggest an association between environments and species, and provide evidence that ecological differentiation between species plays a role in the maintenance of the hybrid zone. Contact areas were characterized by a relatively high temperature, low precipitation, and high topological wetness. Thus, the contact areas were regarded as being located in an unsuitable habitat with a drier climate, where those populations are likely to occur in patches with limited precipitation concentrated. A comparison of spatial scales suggests that exogenous selection via environmental factors may be weaker than endogenous selection via genitalic incompatibility.

Repeated Reticulate Evolution in North American Papilio machaon Group Swallowtail Butterflies

Hybridization between distinct populations or species is increasingly recognized as an important process for generating biodiversity. However, the interaction between hybridization and speciation is complex, and the diverse evolutionary outcomes of hybridization are difficult to differentiate. Here we characterize potential hybridization in a species group of swallowtail butterflies using microsatellites, DNA sequences, and morphology, and assess whether adaptive introgression or homoploid hybrid speciation was the primary process leading to each putative hybrid lineage. Four geographically separated hybrid populations were identified in the Papilio machaon species group. One distinct mitochondrial DNA clade from P. machaon was fixed in three hybrid taxa (P. brevicauda, P. joanae, and P. m. kahli), while one hybrid swarm (P. zelicaon x machaon) exhibited this hybrid mtDNA clade as well as widespread parental mtDNA haplotypes from both parental species. Microsatellite markers and morphology showed variable admixture and intermediacy, ranging from signatures of prolonged differential introgression from the paternal species (P. polyxenes/P. zelicaon) to current gene flow with both parental species. Divergences of the hybrid lineages dated to early- to mid-Pleistocene, suggesting that repeated glaciations and subsequent range shifts of parental species, particularly P. machaon hudsonianus, facilitated initial hybridization. Although each lineage is distinct, P. joanae is the only taxon with sufficient evidence (ecological separation from parental species) to define it as a homoploid hybrid species. The repetition of hybridization in this group provides a valuable foundation for future research on hybridization, and these results emphasize the potential for hybridization to drive speciation in diverse ways.

Improving effectiveness of systematic conservation planning with density data

Systematic conservation planning aims to design networks of protected areas that meet conservation goals across large landscapes. The optimal design of these conservation networks is most frequently based on the modeled habitat suitability or probability of occurrence of species, despite evidence that model predictions may not be highly correlated with species density. We hypothesized that conservation networks designed using species density distributions more efficiently conserve populations of all species considered than networks designed using probability of occurrence models. To test this hypothesis, we used the Zonation conservation prioritization algorithm to evaluate conservation network designs based on probability of occurrence versus density models for 26 land bird species in the U.S. Pacific Northwest. We assessed the efficacy of each conservation network based on predicted species densities and predicted species diversity. High-density model Zonation rankings protected more individuals per species when networks protected the highest priority 10-40% of the landscape. Compared with density-based models, the occurrence-based models protected more individuals in the lowest 50% priority areas of the landscape. The 2 approaches conserved species diversity in similar ways: predicted diversity was higher in higher priority locations in both conservation networks. We conclude that both density and probability of occurrence models can be useful for setting conservation priorities but that density-based models are best suited for identifying the highest priority areas. Developing methods to aggregate species count data from unrelated monitoring efforts and making these data widely available through ecoinformatics portals such as the Avian Knowledge Network will enable species count data to be more widely incorporated into systematic conservation planning efforts.

No evidence for assortative mating within a willow warbler migratory divide

Introduction

In contact zones, genetic mixing of two taxa can be restricted by prezygotic (e.g. assortative mating) or postzygotic (lower fitness of hybrid offspring) barriers, or a combination of the two. A hybrid zone between two willow warbler subspecies (Phylloscopus trochilus trochilus, P. t. acredula) with distinctive migratory strategies occurs in central Sweden. These subspecies exhibit differences in migratory direction and distance, resulting in geographically distinct wintering areas in Africa. The subspecies may have diverged from a common refuge after the last ice age, and neutral genetic markers are homogeneous across their range. By contrast, several phenotypic traits and genetic markers of two chromosomal regions previously identified show steep clines across the divide. The evolutionary forces that maintain this migratory divide remain unknown. Here we use plumage colour, morphology, genetic markers and feather stable nitrogen-isotopes (δ¹⁵N) to assess if assortative mating between migratory phenotypes could be acting as a possible mechanism for keeping the two forms genetically separate and maintaining the migratory divide. We colour-ringed a willow warbler breeding population in the central part of the hybrid zone and observed the breeding population to assess phenotypic and genotypic traits of social pairs.

Results

Our data suggest that wintering area and genetic ancestry had an effect on male arrival time to the breeding grounds which could contribute to assortment. However, evidence for assortative mating could not be detected based on a comparison of plumage colour, morphology and δ¹⁵N between social mates.

Conclusion

This finding was strengthened by analyses of subspecies-specific genetic markers, which allowed us to identify the presence of a large proportion of potential hybrids and backcrosses at the study site. Our results supported the hypothesis that pre-mating isolation in willow warblers is weak, resulting in extensive hybridisation across the migratory divide.

A role for migration-linked genes and genomic islands in divergence of a songbird

Next-generation sequencing has made it possible to begin asking questions about the process of divergence at the level of the genome. For example, recently, there has been a debate around the role of 'genomic islands of divergence' (i.e. blocks of outlier loci) in facilitating the process of speciation-with-gene-flow. The Swainson's thrush, Catharus ustulatus, is a migratory songbird with two genetically distinct subspecies that differ in a number of traits known to be involved in reproductive isolation in birds (plumage coloration, song and migratory behaviour), despite contemporary gene flow along a secondary contact zone. Here, we use RAD-PE sequencing to test emerging hypotheses about the process of divergence at the level of the genome and identify genes and gene regions involved in differentiation in this migratory songbird. Our analyses revealed distinct genomic islands on 15 of the 23 chromosomes and an accelerated rate of divergence on the Z chromosome, one of the avian sex chromosomes. Further, an analysis of loci linked to traits known to be involved in reproductive isolation in songbirds showed that genes linked to migration are significantly more differentiated than expected by chance, but that these genes lie primarily outside the genomic islands. Overall, our analysis supports the idea that genes linked to migration play an important role in divergence in migratory songbirds, but we find no compelling evidence that the observed genomic islands are facilitating adaptive divergence in migratory behaviour.

Unravelling migratory connections: the next level

Migration is widespread among birds, and the strength of the link between the breeding and wintering grounds, migratory connectivity, influences many ecological and evolutionary processes. Despite its importance, migratory connectivity is poorly estimated for most species. Traditionally, visual observations and bird ringing have been used to monitor migration, but these methods require more effort for relatively little return. Genetic markers and stable isotope signatures have increasingly been used to study connectivity. Each approach has its distinct strengths and weaknesses, and as is often the case, a combination may yield the most insight. In this issue of Molecular Ecology, Rundel and colleagues (2013) present a novel Bayesian statistical framework in which genetics and stable isotope data can be combined to improve the assignment of individuals to different winter or breeding regions. The development of such new statistical methods combined with the increasing number and ease of access of isotopic and genetic data sets will greatly enhance our understanding of migratory connectivity. Add to this the developments of miniature devices to track movements of individuals, and the field is destined to make major progression in the decades to come.

How migratory thrushes conquered northern North America: a comparative phylogeography approach

Five species of migratory thrushes (Turdidae) occupy a transcontinental distribution across northern North America. They have largely overlapping breeding ranges, relatively similar ecological niches, and mutualistic relationships with northern woodland communities as insectivores and seed-dispersing frugivores. As an assemblage of ecologically similar species, and given other vertebrate studies, we predicted a shared pattern of genetic divergence among these species between their eastern and western populations, and also that the timing of the coalescent events might be similar and coincident with historical glacial events. To determine how these five lineages effectively established transcontinental distributions, we used mitochondrial cytochrome b sequences to assess genetic structure and lineage coalescence from populations on each side of the continent. Two general patterns occur. Hermit and Swainson’s thrushes (Catharus guttatus and C. ustulatus) have relatively deep divergences between eastern and western phylogroups, probably reflecting shared historic vicariance. The Veery (C. fuscescens), Gray-cheeked Thrush (C. minimus), and American Robin (Turdus migratorius) have relatively shallow divergences between eastern and western populations. However, coalescent and approximate Bayesian computational analyses indicated that among all species as many as five transcontinental divergence events occurred. Divergence within both Hermit and Swainson’s thrushes resembled the divergence between Gray-cheeked Thrushes and Veeries and probably occurred during a similar time period. Despite these species’ ecological similarities, the assemblage exhibits heterogeneity at the species level in how they came to occupy transcontinental northern North America but two general continental patterns at an among-species organizational level, likely related to lineage age.

Contrasting patterns of genetic differentiation among Blackcaps (Sylvia atricapilla) with divergent migratory orientations in Europe

Migratory divides are thought to facilitate behavioral, ecological, and genetic divergence among populations with different migratory routes. However, it is currently contentious how much genetic divergence is needed to maintain distinct migratory behavior across migratory divides. Here we investigate patterns of neutral genetic differentiation among Blackcap (Sylvia atricapilla) populations with different migratory strategies across Europe. We compare the level of genetic divergence of populations migrating to southwestern (SW) or southeastern (SE) wintering areas with birds wintering in the British Isles following a recently established northwesterly (NW) migration route. The migratory divide between SW and SE wintering areas can be interpreted as a result of a re-colonization process after the last glaciation. Thus we predicted greater levels of genetic differentiation among the SW/SE populations. However, a lack of genetic differentiation was found between SW and SE populations, suggesting that interbreeding likely occurs among Blackcaps with different migratory orientations across a large area; therefore the SW/SE migratory divide can be seen as diffuse, broad band and is, at best, a weak isolating barrier. Conversely, weak, albeit significant genetic differentiation was evident between NW and SW migrants breeding sympatrically in southern Germany, suggesting a stronger isolating mechanism may be acting in this population. Populations located within/near the SW/SE contact zone were the least genetically divergent from NW migrants, confirming NW migrants likely originated from within the contact zone. Significant isolation-by-distance was found among eastern Blackcap populations (i.e. SE migrants), but not among western populations (i.e. NW and SW migrants), revealing different patterns of genetic divergence among Blackcap populations in Europe. We discuss possible explanations for the genetic structure of European Blackcaps and how gene flow influences the persistence of divergent migratory behaviors.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

Modelling the spread of Wolbachia in spatially heterogeneous environments

The endosymbiont Wolbachia infects a large number of insect species and is capable of rapid spread when introduced into a novel host population. The bacteria spread by manipulating their hosts' reproduction, and their dynamics are influenced by the demographic structure of the host population and patterns of contact between individuals. Reaction-diffusion models of the spatial spread of Wolbachia provide a simple analytical description of their spatial dynamics but do not account for significant details of host population dynamics. We develop a metapopulation model describing the spatial dynamics of Wolbachia in an age-structured host insect population regulated by juvenile density-dependent competition. The model produces similar dynamics to the reaction-diffusion model in the limiting case where the host's habitat quality is spatially homogeneous and Wolbachia has a small effect on host fitness. When habitat quality varies spatially, Wolbachia spread is usually much slower, and the conditions necessary for local invasion are strongly affected by immigration of insects from surrounding regions. Spread is most difficult when variation in habitat quality is spatially correlated. The results show that spatial variation in the density-dependent competition experienced by juvenile host insects can strongly affect the spread of Wolbachia infections, which is important to the use of Wolbachia to control insect vectors of human disease and other pests.

North American transmission of hemosporidian parasites in the Swainson's thrush (Catharus ustulatus), a migratory songbird

The geographic structuring of parasite communities across the range of a single host species can illuminate patterns of host-population connectivity. To determine the location of parasite transmission in a Neotropical migrant bird species, we sampled adult and hatch-year (HY) birds across the breeding and wintering range of the Swainson's thrush (SWTH), an abundant passerine with a migratory divide. We examined the phylogenetic relationships among cytochrome b lineages of the avian blood parasite genera Haemoproteus, Plasmodium, and Leucocytozoon and determined the transmission location of unique lineages. We found that Haemoproteus and Plasmodium lineages are transmitted on California breeding grounds, whereas Leucocytozoon transmission occurs on Alaskan breeding grounds. The presence of hemosporidians on wintering grounds and shared lineages between the SWTH and resident species suggests that transmission of some of these lineages occurs on both breeding and wintering grounds. We emphasize that the sampling of HY birds and local resident heterospecifics will supplement vector studies to determine the key players in hemosporidian host switching and range-expansion events.

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.