Strong reproductive isolation and narrow genomic tracts of differentiation among three woodpecker species in secondary contact

Bayesian hybrid index and genomic cline estimation with the R package gghybrid

Admixture, the interbreeding of individuals from differentiated source populations, is now known to be a widespread phenomenon. Genomic studies of natural hybridisation can help to answer many questions on the impacts of admixture on adaptive evolution, reproductive isolation, and speciation. When a large variety of admixture proportions between two source populations exist, both geographic and genomic cline analysis are suitable methods for inferring biased, restricted or excessive gene flow at individual loci into the foreign genomic background, providing evidence for reproductive isolation, selection across an environmental transition, balancing selection, and adaptive introgression. Genomic cline analysis replaces geographic location with genome-wide hybrid index and is therefore useable in circumstances that violate geographic cline assumptions. Here, I introduce gghybrid, an R package for simple and flexible Bayesian estimation of Buerkle's hybrid index and Fitzpatrick's logit-logistic genomic clines using bi-allelic data, suitable for both small and large datasets. gghybrid allows any ploidy and uses Structure input file format. It has separate functions for hybrid index and cline estimation, treating each individual and locus respectively as an independent analysis, making it highly parallelisable. Admixture proportions from other software can alternatively be used in cline analysis, alongside parental allele frequencies. Parameters can be fixed and samples pooled for statistical model comparison with AIC or waic. Here, I describe the functions, pipeline, and statistical properties of gghybrid. Simulations reveal that model comparison with waic is preferred, and use of Bayesian posterior distributions and p values to select candidate non-null loci is problematic and should be avoided.

Habitat and climate influence hybridization among three genetically distinct Canada jay (Perisoreus canadensis) morphotypes in an avian hybrid zone complex

Examining the frequency and distribution of hybrids across contact zones provide insights into the factors mediating hybridization. In this study, we examined the effect of habitat and climate on hybridization patterns for three phenotypically, genetically, and ecologically distinct groups of the Canada jay (Perisoreus canadensis) in a secondary contact zone in western North America. Additionally, we tested whether the frequency of hybridization involving the three groups (referred to as Boreal, Pacific and Rocky Mountain morphotypes) is similar across the hybrid zones or whether some pairs have hybridized more frequently than others. We reanalyzed microsatellite, mtDNA and plumage data, and new microsatellite and plumage data for 526 individuals to identify putative genetic and phenotypic hybrids. The genetically and phenotypically distinct groups are associated with different habitats and occupy distinct climate niches across the contact zone. Most putative genetic hybrids (86%) had Rocky Mountain ancestry. Hybrids were observed most commonly in intermediate climate niches and in habitats where Engelmann spruce (Picea engelmannii) overlaps broadly with boreal and subalpine tree species. Our finding that hybrids occupy intermediate climate niches relative to parental morphotypes matches patterns for other plant and animal species found in this region. This study demonstrates how habitat and climate influence hybridization patterns in areas of secondary contact and adds to the growing body of research on tri-species hybrid zones.

Population genomics of an emergent tri-species hybrid zone

Isolating barriers that drive speciation are commonly studied in the context of two-species hybrid zones. There is, however, evidence that more complex introgressive relationships are common in nature. Here, we use field observations and genomic analysis, including the sequencing and assembly of a novel reference genome, to study an emergent hybrid zone involving two colliding hybrid zones of three woodpecker species: red-breasted, red-naped, and yellow-bellied sapsuckers (Sphyrapicus ruber, S. nuchalis, and S. varius). Surveys of the area surrounding Prince George, British Columbia, Canada, show that all three species are sympatric, and Genotyping-by-Sequencing identifies hybrids from each species pair and birds with ancestry from all three species. Observations of phenotypes and genotypes of mated pairs provide evidence for assortative mating, though there is some heterospecific pairing. Hybridization is more extensive in this tri-species hybrid zone than in two di-species hybrid zones. However, there is no evidence of a hybrid swarm and admixture is constrained to contact zones, so we classify this region as a tension zone and invoke selection against hybrids as a likely mechanism maintaining species boundaries. Analysis of sapsucker age classes does not show disadvantages in hybrid survival to adulthood, so we speculate the selection upholding the tension zone may involve hybrid fecundity. Gene flow among all sapsuckers in di-species hybrid zones suggests introgression probably occurred before the formation of this tri-species hybrid zone, and might result from bridge hybridization, vagrancies, or other three-species interactions.

Weak coupling among barrier loci and waves of neutral and adaptive introgression across an expanding hybrid zone

Hybridization can serve as an evolutionary stimulus, but we have little understanding of introgression at early stages of hybrid zone formation. We analyze reproductive isolation and introgression between a range‐limited and a widespread species. Reproductive barriers are estimated based on differences in flowering time, ecogeographic distributions, and seed set from crosses. We find an asymmetrical mating barrier due to cytonuclear incompatibility that is consistent with observed clusters of coincident and concordant tension zone clines (barrier loci) for mtDNA haplotypes and nuclear SNPs. These groups of concordant clines are spread across the hybrid zone, resulting in weak coupling among barrier loci and extensive introgression. Neutral clines had nearly equal introgression into both species’ ranges, whereas putative cases of adaptive introgression had exceptionally wide clines with centers shifted toward one species. Analyses of cline shape indicate that secondary contact was initiated within the last 800 generations with the per‐generation dispersal between 200 and 400 m, and provide some of the first estimates of the strength of selection required to account for observed levels of adaptive introgression. The weak species boundary between these species appears to be in early stages of dissolution, and ultimately will precipitate genetic swamping of the range‐limited species.

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell's wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.

Cryptic genetic diversity and cytonuclear discordance characterize contact among Canada jay (Perisoreus canadensis) morphotypes in western North America

Three distinct Canada jay (Perisoreus canadensis) morphotypes with easily recognizable plumage traits come into contact in western North America. Recent work demonstrated high genetic structure across the species’ range; however, patterns of genetic variation in these contact zones remain unknown. We categorized 605 individuals into one of three morphotypes (Pacific, Rocky Mountain, and Boreal) based on plumage, and genotyped individuals at the mtDNA control region and 12 microsatellite loci to assess the extent of hybridization between morphotypes. Our data showed cryptic genetic diversity and high cytonuclear discordance among morphotypes within contact zones, which is likely the result of recent and historical admixture. The distributions of the Boreal and Pacific morphotypes each showed a strong association with a single, distinct genetic group, whereas the Rocky Mountain morphotype exhibited higher genetic diversity and was associated with multiple genotypes. Our analyses show the importance of considering both plumage and genetic traits when examining contact zones between closely related taxa. Finally the data presented in this study reaffirm that the Pacific morphotype is distinct from the Boreal and Rocky Mountain morphotypes based on genetic, phenotypic and ecological data, indicating that the Pacific morphotype should be re-elevated to a full species.

Recent introgression between Taiga Bean Goose and Tundra Bean Goose results in a largely homogeneous landscape of genetic differentiation

Several studies have uncovered a highly heterogeneous landscape of genetic differentiation across the genomes of closely related species. Specifically, genetic differentiation is often concentrated in particular genomic regions (“islands of differentiation”) that might contain barrier loci contributing to reproductive isolation, whereas the rest of the genome is homogenized by introgression. Alternatively, linked selection can produce differentiation islands in allopatry without introgression. We explored the influence of introgression on the landscape of genetic differentiation in two hybridizing goose taxa: the Taiga Bean Goose (Anser fabalis) and the Tundra Bean Goose (A. serrirostris). We re-sequenced the whole genomes of 18 individuals (9 of each taxon) and, using a combination of population genomic summary statistics and demographic modeling, we reconstructed the evolutionary history of these birds. Next, we quantified the impact of introgression on the build-up and maintenance of genetic differentiation. We found evidence for a scenario of allopatric divergence (about 2.5 million years ago) followed by recent secondary contact (about 60,000 years ago). Subsequent introgression events led to high levels of gene flow, mainly from the Tundra Bean Goose into the Taiga Bean Goose. This scenario resulted in a largely undifferentiated genomic landscape (genome-wide F_ST = 0.033) with a few notable differentiation peaks that were scattered across chromosomes. The summary statistics indicated that some peaks might contain barrier loci while others arose in allopatry through linked selection. Finally, based on the low genetic differentiation, considerable morphological variation and incomplete reproductive isolation, we argue that the Taiga and the Tundra Bean Goose should be treated as subspecies.

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.

Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early-generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re-established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context-dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.

Population Genetics and Speciation of Yellow-Bellied, Red-Naped, and Red-Breasted Sapsuckers (Sphyrapicus varius, S. nuchalis, and S. ruber)

The root of understanding speciation lies in determining the forces which drive it. In many closely-related species, including Sphyrapicus varius, S. nuchalis, and S. ruber, it is assumed that speciation occurred due to isolation in multiple Pleistocene refugia. We used genetic data from 457 samples at the control region (CR), cytochrome oxidase I (COI), and chromo-helicase DNA binding protein (CHD1Z) to examine range-wide population genetic structure and differentiation amongst these 3 species across each species' breeding range. In addition, we modeled these species' ecological niches for the Holocene (~6000 years ago), Last Glacial Maximum (~22000 years ago), and Last Interglacial (~120000-140000 years ago) to determine if Pleistocene glaciations could have contributed to allopatric distributions, therefore allowing these groups to differentiate. Population genetic data show a potential Pleistocene refugium in Haida Gwaii, an east-west split among S. varius, and low genetic differentiation within each species. Our CR data show some polyphyly, while COI and CHD1Z data show differentiation among species using composite genotypes. Ecological Niche Modeling shows a large amount of niche overlap at each time period suggesting that S. varius, S. nuchalis, and S. ruber may not have been completely allopatric, and these species likely had repeated intermittent contact. Our data support the growing body of research that suggests differentiation despite gene flow.

Morphologically cryptic Amazonian bird species pairs exhibit strong postzygotic reproductive isolation

We possess limited understanding of how speciation unfolds in the most species-rich region of the planet-the Amazon basin. Hybrid zones provide valuable information on the evolution of reproductive isolation, but few studies of Amazonian vertebrate hybrid zones have rigorously examined the genome-wide underpinnings of reproductive isolation. We used genome-wide genetic datasets to show that two deeply diverged, but morphologically cryptic sister species of forest understorey birds show little evidence for prezygotic reproductive isolation, but substantial postzygotic isolation. Patterns of heterozygosity and hybrid index revealed that hybrid classes with heavily recombined genomes are rare and closely match simulations with high levels of selection against hybrids. Genomic and geographical clines exhibit a remarkable similarity across loci in cline centres, and have exceptionally narrow cline widths, suggesting that postzygotic isolation is driven by genetic incompatibilities at many loci, rather than a few loci of strong effect. We propose Amazonian understorey forest birds speciate slowly via gradual accumulation of postzygotic genetic incompatibilities, with prezygotic barriers playing a less important role. Our results suggest old, cryptic Amazonian taxa classified as subspecies could have substantial postzygotic isolation deserving species recognition and that species richness is likely to be substantially underestimated in Amazonia.

Genomic and geographic footprints of differential introgression between two divergent fish species (Solea spp.)

Investigating gene flow between closely related species and its variation across the genome is important to understand how reproductive barriers shape genome divergence before speciation is complete. An efficient way to characterize differential gene flow is to study how the genetic interactions that take place in hybrid zones selectively filter gene exchange between species, leading to heterogeneous genome divergence. In the present study, genome-wide divergence and introgression patterns were investigated between two sole species, Solea senegalensis and Solea aegyptiaca, using restriction-associated DNA sequencing (RAD-Seq) to analyze samples taken from a transect spanning the hybrid zone. An integrative approach combining geographic and genomic clines methods with an analysis of individual locus introgression accounting for the demographic history of divergence was conducted. Our results showed that the two sole species have come into secondary contact postglacially, after experiencing a prolonged period (ca. 1.1 to 1.8 Myrs) of allopatric separation. Secondary contact resulted in the formation of a tension zone characterized by strong reproductive isolation, which only allowed introgression in a limited fraction of the genome. We found multiple evidence for a preferential direction of introgression in the S. aegyptiaca genetic background, indicating a possible recent or ongoing movement of the hybrid zone. Deviant introgression signals found in the opposite direction suggested that S. senegalensis could have possibly undergone adaptive introgression that has not yet spread throughout the entire species range. Our study thus illustrates the varied outcomes of genetic interactions between divergent gene pools that recently met after a long history of divergence.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.