Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry

Evolutionarily labile dispersal behavior and discontinuous habitats enhance population differentiation in island versus continentally distributed swallows

The causes of population divergence in vagile groups remain a paradox in evolutionary biology: dispersive species should be able to colonize new areas, a prerequisite for allopatric speciation, but dispersal also facilitates gene flow, which erodes population differentiation. Strong dispersal ability has been suggested to enhance divergence in patchy habitats and inhibit divergence in continuous landscapes, but empirical support for this hypothesis is lacking. Here we compared patterns of population divergence in a dispersive clade of swallows distributed across both patchy and continuous habitats. The Pacific Swallow (Hirundo tahitica) has an insular distribution throughout Southeast Asia and the Pacific, while its sister species, the Welcome Swallow (H. neoxena), has a continental distribution in Australia. We used whole-genome data to demonstrate strong genetic structure and limited introgression among insular populations, but not among continental populations. Demographic models show that historic changes in habitat connectivity have contributed to population structure within the clade. Swallows appear to exhibit evolutionarily labile dispersal behavior in which they reduce dispersal propensity after island colonization despite retaining strong flight ability. Our data support the hypothesis that fragmented habitats enhance population differentiation in vagile groups, and suggest that labile dispersal behavior is a key mechanism underlying this pattern.

Quantitative Analyses of Coupling in Hybrid Zones

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

Population genetic structure and hybrid zone analyses for species delimitation in the Japanese toad (Bufo japonicus)

Hybridization following secondary contact may produce different outcomes depending on the extent to which genetic diversity and reproductive barriers have accumulated during isolation. The Japanese toad, Bufo japonicus, is distributed on the main islands of Japan. In the present study, we applied multiplexed inter-simple sequence repeat genotyping by sequencing to achieve the fine-scale resolution of the genetic cluster in B. j. japonicus and B. j. formosus. We also elucidated hybridization patterns and gene flow degrees across contact zones between the clusters identified. Using SNP data, we found four genetic clusters in B. j. japonicus and B. j. formosus and three contact zones of the cluster pairs among these four clusters. The two oldest diverged lineages, B. j. japonicus and B. j. formosus, formed a narrow contact zone consistent with species distinctiveness. Therefore, we recommend that these two subspecies be elevated to the species level. In contrast, the less diverged pairs of two clusters in B. j. japonicus and B. j. formosus, respectively, admixed over a hundred kilometers, suggesting that they have not yet developed strong reproductive isolation and need to be treated as conspecifics. These results will contribute to resolving taxonomic confusion in Japanese toads.

Conspicuous coloration of toxin-resistant predators implicates additional trophic interactions in a predator-prey arms race

Antagonistic coevolution between natural enemies can produce highly exaggerated traits, such as prey toxins and predator resistance. This reciprocal process of adaptation and counter-adaptation may also open doors to other evolutionary novelties not directly involved in the phenotypic interface of coevolution. We tested the hypothesis that predator-prey coevolution coincided with the evolution of conspicuous coloration on resistant predators that retain prey toxins. In western North America, common garter snakes (Thamnophis sirtalis) have evolved extreme resistance to tetrodotoxin (TTX) in the coevolutionary arms race with their deadly prey, Pacific newts (Taricha spp.). TTX-resistant snakes can retain large amounts of ingested TTX, which could serve as a deterrent against the snakes' own predators if TTX toxicity and resistance are coupled with a conspicuous warning signal. We evaluated whether arms race escalation covaries with bright red coloration in snake populations across the geographic mosaic of coevolution. Snake colour variation departs from the neutral expectations of population genetic structure and covaries with escalating clines of newt TTX and snake resistance at two coevolutionary hotspots. In the Pacific Northwest, bright red coloration fits an expected pattern of an aposematic warning to avian predators: TTX-resistant snakes that consume highly toxic newts also have relatively large, reddish-orange dorsal blotches. Snake coloration also seems to have evolved with the arms race in California, but overall patterns are less intuitively consistent with aposematism. These results suggest that interactions with additional trophic levels can generate novel traits as a cascading consequence of arms race coevolution across the geographic mosaic.

A chromosome-level reference genome and pangenome for barn swallow population genomics

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

High heterogeneity in genomic differentiation between phenotypically divergent songbirds: a test of mitonuclear co-introgression

Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This discordance is best explained by adaptive mtDNA introgression-a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand these discordant differentiation patterns and characterize nuclear differentiation in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further investigated mitonuclear gene co-introgression between yellowhammers and pine buntings and found support for this process in the direction of pine buntings into yellowhammers. Genomic signals indicative of co-introgression were common in mitonuclear genes coding for subunits of the mitoribosome and electron transport chain complexes. Such introgression of mitochondrial DNA and mitonuclear genes provides a possible explanation for the patterns of high genomic heterogeneity in genomic differentiation seen among some species groups.

Ecological outcomes of hybridization vary extensively in Catostomus fishes

Hybridization outcomes vary geographically and can depend on the environment. Hybridization can also reshape biotic interactions, leading to ecological shifts. If hybrids function differently ecologically in ways that enhance or reduce fitness, and those ecological roles vary geographically, ecological factors might explain variation in hybridization outcomes. However, relatively few studies have focused on ecological traits of hybrids. We compared the feeding ecology of Catostomus fish species and hybrids by using stable isotopes (δ13 C and δ15 N) as a proxy for diet and habitat use, and compared two native species, an introduced species, and three interspecific hybrid crosses. We included hybrids and parental species from seven rivers where hybridization outcomes vary. Relative isotopic niches of native species varied geographically, but native species did not fully overlap in isotopic space in any river sampled, suggesting little overlap of resource use between historically sympatric species. The introduced species overlapped with one or both native species in every river, suggesting similar resource use and potential competition. Hybrids occupied intermediate, matching, or more transgressive isotopic niches, and varied within and among rivers. Ecological outcomes of hybridization varied across locations, implying that hybridization might have unpredictable, idiosyncratic ecological effects.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica)

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (< 20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

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.

Population genomics of free-ranging Great Plains white-tailed and mule deer reflects a long history of interspecific hybridization

Hybridization is a natural process at species-range boundaries that may variably promote the speciation process or break down species barriers but minimally will influence management outcomes of distinct populations. White-tailed deer (Odocoileus virginianus) and mule deer (Odocoileus hemionus) have broad and overlapping distributions in North America and a recognized capacity for interspecific hybridization. In response to contemporary environmental change to any of one or multiple still-unknown factors, mule deer range is contracting westward accompanied by a westward expansion of white-tailed deer, leading to increasing interactions, opportunities for gene flow, and associated conservation implications. To quantify genetic diversity, phylogenomic structure, and dynamics of hybridization in sympatric populations of white-tailed and mule deer, we used mitochondrial cytochrome b data coupled with SNP loci discovered with double-digest restriction site-associated DNA sequencing. We recovered 25,018 SNPs across 92 deer samples from both species, collected from two regions of western Kansas. Eight individuals with unambiguous external morphology representing both species were of hybrid origin (8.7%), and represented the product of multi-generational backcrossing. Mitochondrial data showed both ancient and recent directional discordance with morphological species assignments, reflecting a legacy of mule deer males mating with white-tailed deer females. Mule deer had lower genetic diversity than white-tailed deer, and both mitochondrial and nuclear data suggest contemporary mule deer effective population decline. Landscape genetic analyses show relative isolation between the two study regions for white-tailed deer, but greater connectivity among mule deer, with predominant movement from north to south. Collectively, our results suggest a long history of gene flow between these species in the Great Plains and hint at evolutionary processes that purge incompatible functional genomic elements as a result of hybridization. Surviving hybrids evidently may be reproductive, but with unknown consequences for the future integrity of these species, population trajectories, or relative susceptibility to emerging pathogens.

A de novo genome assembly and annotation of the southern flying squirrel (Glaucomys volans)

Northern (Glaucomys sabrinus) and southern (Glaucomys volans) flying squirrels are widespread species distributed across North America. Northern flying squirrels are common inhabitants of the boreal forest, also occurring in coniferous forest remnants farther south, whereas the southern flying squirrel range is centered in eastern temperate woodlands. These two flying squirrel species exhibit a hybrid zone across a latitudinal gradient in an area of recent secondary contact. Glaucomys hybrid offspring are viable and can successfully backcross with either parental species, however, the fitness implications of such events are currently unknown. Some populations of G. sabrinus are endangered, and thus, interspecific hybridization is a key conservation concern in flying squirrels. To provide a resource for future studies to evaluate hybridization and possible introgression, we sequenced and assembled a de novo long-read genome from a G. volans individual sampled in southern Ontario, Canada, while four short-read genomes (two G. sabrinus and two G. volans, all from Ontario) were resequenced on Illumina platforms. The final genome assembly consisted of approximately 2.40 Gb with a scaffold N50 of 455.26 Kb. Benchmarking Universal Single-Copy Orthologs reconstructed 3,742 (91.2%) complete mammalian genes and genome annotation using RNA-Seq identified the locations of 19,124 protein-coding genes. The four short-read individuals were aligned to our reference genome to investigate the demographic history of the two species. A principal component analysis clearly separated resequenced individuals, while inferring population size history using the Pairwise Sequentially Markovian Coalescent model noted an approximate species split 1 million years ago, and a single, possibly recently introgressed individual.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Sex-linked genetic diversity and differentiation in a globally distributed avian species complex

Sex chromosomes often bear distinct patterns of genetic variation due to unique patterns of inheritance and demography. The processes of mutation, recombination, genetic drift and selection also influence rates of evolution on sex chromosomes differently than autosomes. Measuring such differences provides information about how these processes shape genomic variation and their roles in the origin of species. To test hypotheses and predictions about patterns of autosomal and sex-linked genomic diversity and differentiation, we measured population genetic statistics within and between populations and subspecies of the barn swallow (Hirundo rustica) and performed explicit comparisons between autosomal and Z-linked genomic regions. We first tested for evidence of low Z-linked genetic diversity and high Z-linked population differentiation relative to autosomes, then for evidence that the Z chromosome bears greater ancestry information due to faster lineage sorting. Finally, we investigated geographical clines across hybrid zones for evidence that the Z chromosome is resistant to introgression due to selection against hybrids. We found evidence that the barn swallow mating system, demographic history and linked selection each contribute to low Z-linked diversity and high Z-linked differentiation. While incomplete lineage sorting is rampant across the genome, our results indicate faster sorting of ancestral polymorphism on the Z. Finally, hybrid zone analyses indicate barriers to introgression on the Z chromosome, suggesting that sex-linked traits are important in reproductive isolation, especially in migratory divide regions. Our study highlights how selection, gene flow and demography shape sex-linked genetic diversity and underlines the relevance of the Z chromosome in speciation.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Population genetic structure underlying the geographic variation in beetle structural colour with multiple transition zones

We studied the population genetic structure underlying the geographic variation in the structural colour of the geotrupid dung beetle, Phelotrupes auratus, which exhibits metallic body colours of different reflectance wavelengths perceived as red, green and indigo. These forms occur parapatrically in an area of Japan. The colour variation was not related to variation in climatic factors. Using single nucleotide polymorphisms (SNPs) from restriction-site-associated DNA sequences, we discriminated five groups of populations (west/red, south/green, south/indigo, south/red and east/red) by a combination of genetic clusters (west, south and east) and three colour forms. There were three transition zones for the colour forms: two between the red and green forms were hybrid zones with steep genetic clines, which implies the existence of barriers to gene flow between regions with different colours. The remaining transition zone between the green and indigo forms lacked genetic differentiation, despite the evident colour changes, which implies regionally specific selection on the different colours. In a genomewide association study, we identified four SNPs that were associated with the red/green or indigo colour and were not linked with one another, which implies that the coloration was controlled by multiple loci, each affecting the expression of a different colour range. These loci may have controlled the transitions between different combinations of colours. Our study demonstrates that geographic colour variation within a species can be maintained by nonuniform interactions among barriers to gene flow, locally specific selection on different colours, and the effects of different colour loci.

Genetic diversity and structure of the hedgehogs Erinaceus europaeus and Erinaceus roumanicus: evidence for ongoing hybridization in Eastern Europe

Secondary contact zones between related species are key to understanding speciation mechanisms. The Central European sympatry zone of West European (Erinaceus europaeus) and northern white-breasted (Erinaceus roumanicus) hedgehogs is well studied, whereas data on the Eastern European sympatry zone are scarce. We examined the genetic variation in Russian populations using the mitochondrial Cytb gene, TTR intron 1 and 11 microsatellites to assess genetic variability and distribution patterns. In contrast to the Central European sympatry zone, we found evidence of ongoing hybridization between the two species in the sympatry zone of European Russia, where the proportion of individuals with mixed ancestry was c. 20%. Our data indicate bi-directional mtDNA introgression, but with a higher frequency of E. europaeus haplotypes in hybrids. The proportion of pure specimens with introgressed mitotypes is higher in E. roumanicus than in E. europaeus. Nuclear data showed the prevalence of the genetic contribution from E. roumanicus in admixed individuals. Demographic analyses indicated recent population growth in E. europaeus and little change in E. roumanicus, suggesting that E. europaeus colonized East Europe later than E. roumanicus.

CYP2J19 mediates carotenoid colour introgression across a natural avian hybrid zone

It has long been of interest to identify the phenotypic traits that mediate reproductive isolation between related species, and more recently, the genes that underpin them. Much work has focused on identifying genes associated with animal colour, with the candidate gene CYP2J19 identified in laboratory studies as the ketolase converting yellow dietary carotenoids to red ketocarotenoids in birds with red pigments. However, evidence that CYP2J19 explains variation between red and yellow feather coloration in wild populations of birds is lacking. Hybrid zones provide the opportunity to identify genes associated with specific traits. Here we investigate genomic regions associated with colour in red-fronted and yellow-fronted tinkerbirds across a hybrid zone in southern Africa. We sampled 85 individuals, measuring spectral reflectance of forecrown feathers and scoring colours from photographs, while testing for carotenoid presence with Raman spectroscopy. We performed a genome-wide association study to identify associations with carotenoid-based coloration, using double-digest RAD sequencing aligned to a short-read whole genome of a Pogoniulus tinkerbird. Admixture mapping using 104,933 single nucleotide polymorphisms (SNPs) identified a region of chromosome 8 that includes CYP2J19 as the only locus with more than two SNPs significantly associated with both crown hue and crown score, while Raman spectra provided evidence of ketocarotenoids in red feathers. Asymmetric backcrossing in the hybrid zone suggests that yellow-fronted females mate more often with red-fronted males than vice versa. Female red-fronted tinkerbirds mating assortatively with red-crowned males is consistent with the hypothesis that converted carotenoids are an honest signal of quality.

Migration, pathogens and the avian microbiome: A comparative study in sympatric migrants and residents

Animals generally benefit from their gastrointestinal microbiome, but the factors that influence the composition and dynamics of their microbiota remain poorly understood. Studies of nonmodel host species can illuminate how microbiota and their hosts interact in natural environments. We investigated the role of migratory behaviour in shaping the gut microbiota of free-ranging barn swallows (Hirundo rustica) by studying co-occurring migrant and resident subspecies sampled during the autumn migration at a migratory bottleneck. We found that within-host microbial richness (α-diversity) was similar between migrant and resident microbial communities. In contrast, we found that microbial communities (β-diversity) were significantly different between groups regarding both microbes present and their relative abundances. Compositional differences were found for 36 bacterial genera, with 27 exhibiting greater abundance in migrants and nine exhibiting greater abundance in residents. There was heightened abundance of Mycoplasma spp. and Corynebacterium spp. in migrants, a pattern shared by other studies of migratory species. Screens for key regional pathogens revealed that neither residents nor migrants carried avian influenza viruses and Newcastle disease virus, suggesting that the status of these diseases did not underlie observed differences in microbiome composition. Furthermore, the prevalence and abundance of Salmonella spp., as determined from microbiome data and cultural assays, were both low and similar across the groups. Overall, our results indicate that microbial composition differs between migratory and resident barn swallows, even when they are conspecific and sympatrically occurring. Differences in host origins (breeding sites) may result in microbial community divergence, and varied behaviours throughout the annual cycle (e.g., migration) could further differentiate compositional structure as it relates to functional needs.

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental-like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.

The geographic mosaic of arms race coevolution is closely matched to prey population structure

Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of "trait remixing" in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

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.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.

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.

Biochemical and structural insights into the cytochrome P450 reductase from Candida tropicalis

Cytochrome P450 reductases (CPRs) are diflavin oxidoreductases that supply electrons to type II cytochrome P450 monooxygenases (CYPs). In addition, it can also reduce other proteins and molecules, including cytochrome c, ferricyanide, and different drugs. Although various CPRs have been functionally and structurally characterized, the overall mechanism and its interaction with different redox acceptors remain elusive. One of the main problems regarding electron transfer between CPRs and CYPs is the so-called "uncoupling", whereby NAD(P)H derived electrons are lost due to the reduced intermediates' (FAD and FMN of CPR) interaction with molecular oxygen. Additionally, the decay of the iron-oxygen complex of the CYP can also contribute to loss of reducing equivalents during an unproductive reaction cycle. This phenomenon generates reactive oxygen species (ROS), leading to an inefficient reaction. Here, we present the study of the CPR from Candida tropicalis (CtCPR) lacking the hydrophobic N-terminal part (Δ2-22). The enzyme supports the reduction of cytochrome c and ferricyanide, with an estimated 30% uncoupling during the reactions with cytochrome c. The ROS produced was not influenced by different physicochemical conditions (ionic strength, pH, temperature). The X-ray structures of the enzyme were solved with and without its cofactor, NADPH. Both CtCPR structures exhibited the closed conformation. Comparison with the different solved structures revealed an intricate ionic network responsible for the regulation of the open/closed movement of CtCPR.

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by-product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced-representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus-specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by-product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.

SMRT long reads and Direct Label and Stain optical maps allow the generation of a high-quality genome assembly for the European barn swallow (Hirundo rustica rustica)

Background

The barn swallow (Hirundo rustica) is a migratory bird that has been the focus of a large number of ecological, behavioral, and genetic studies. To facilitate further population genetics and genomic studies, we present a reference genome assembly for the European subspecies (H. r. rustica).

Findings

As part of the Genome10K effort on generating high-quality vertebrate genomes (Vertebrate Genomes Project), we have assembled a highly contiguous genome assembly using single molecule real-time (SMRT) DNA sequencing and several Bionano optical map technologies. We compared and integrated optical maps derived from both the Nick, Label, Repair, and Stain technology and from the Direct Label and Stain (DLS) technology. As proposed by Bionano, DLS more than doubled the scaffold N50 with respect to the nickase. The dual enzyme hybrid scaffold led to a further marginal increase in scaffold N50 and an overall increase of confidence in the scaffolds. After removal of haplotigs, the final assembly is approximately 1.21 Gbp in size, with a scaffold N50 value of more than 25.95 Mbp.

Conclusions

This high-quality genome assembly represents a valuable resource for future studies of population genetics and genomics in the barn swallow and for studies concerning the evolution of avian genomes. It also represents one of the very first genomes assembled by combining SMRT long-read sequencing with the new Bionano DLS technology for scaffolding. The quality of this assembly demonstrates the potential of this methodology to substantially increase the contiguity of genome assemblies.

Demographic inference in barn swallows using whole-genome data shows signal for bottleneck and subspecies differentiation during the Holocene

Accounting for historical demographic features is vital for many types of evolutionary inferences, including the estimation of divergence times between closely related populations. In barn swallow, Hirundo rustica, inferring historical population sizes and subspecies divergence times can shed light on the recent co-evolution of this species with humans. Pairwise sequentially Markovian coalescent uncovered population growth beginning on the order of one million years ago-which may reflect the radiation of the broader Hirundo genus-and a more recent population decline. Additionally, we used approximate Bayesian computation to evaluate hypotheses about recent timescale barn swallow demography, including population growth due to human commensalism, and a potential founder event associated with the onset of nesting on human structures. We found signal for a bottleneck event approximately 7,700 years ago, near the time that humans began building substantial structures, although there was considerable uncertainty associated with this estimate. Subspecies differentiation and subsequent growth occurred after the bottleneck in the best-supported model, an order of magnitude more recently than previous estimates in this system. We also compared results obtained from whole-genome sequencing versus reduced representation sequencing, finding many similar results despite substantial allelic dropout in the reduced representation data, which may have affected estimates of some parameters. This study presents the first genetic evidence of a potential barn swallow founder effect and subspecies divergence coinciding with the Holocene, which is an important step in analysing the biogeographical history of a well-known human commensal species.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.