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.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Geographically extensive hybridization between the forest trees American butternut and Japanese walnut

We investigate the question of naturally occurring interspecific hybrids between two forest trees: the native North American butternut (Juglans cinerea L.) and the introduced Japanese walnut (Juglans ailantifolia Carrière). Using nuclear and chloroplast DNA markers, we provide evidence for 29 F₁ and 22 advanced generation hybrids in seven locations across the eastern and southern range of the native species. Two locations show extensive admixture (95% J. ailantifolia and hybrids) while other locations show limited admixture. Hybridization appears to be asymmetrical with 90.9 per cent of hybrids having J. ailantifolia as the maternal parent. This is, to our knowledge, the first genetic data supporting natural hybridization between these species. The long-term outcome of introgression could include loss of native diversity, but could also include transfer of useful traits from the introduced species.     

Population genetic structure of Picea engelmannii, P. glauca and their previously unrecognized hybrids in the central Rocky Mountains

Areas of geographic overlap between potentially hybridizing species provide the opportunity to study interspecific gene flow and reproductive barriers. Here we identified hybrids between Picea engelmannii and P. glauca by their genetic composition at 17 microsatellite markers, and determined the broad-scale geographic distribution of hybrids in the central Rocky Mountains of North America, a geographic region where hybrids and isolation between species had not previously been studied. Parameter estimates from admixture models revealed considerable variation in ancestry within and among collection sites, suggesting that within this area of geographic overlap, the interaction of the two species varies extensively. The results document a previously unrecognized patchy distribution of hybrids between P. engelmannii and P. glauca, including locations where hybrids were not known or expected to exist. Further, the ancestry of many hybrids was consistent with multiple generations of hybridization, with probable directional backcrossing to P. engelmannii, suggesting a relatively porous species boundary. The identification and characterization of hybridization between these spruce in this region raises the question of what factors maintain barriers to gene flow in these long-lived forest trees. The current research lays the groundwork for future study of the ecological and evolutionary contexts of their hybridization, as well as of differential introgression and permeability of species boundaries.     

Hybridization and divergence in multi‐species oak (Quercus) communities

Oaks (Quercus: Fagaceae) commonly interbreed yet retain their morphological, genetic and ecological distinctiveness. Post‐zygotic isolation mechanisms, such as ecologically dependent selection on adaptive loci, may therefore limit introgression. To test this hypothesis, we quantified hybridization and genetic divergence across the contact zone of four red oaks (Quercus section Lobatae) in the Great Lakes region of North America using a suite of 259 amplified fragment length polymorphisms and 27 genic and genomic microsatellite markers. First, we identified hybrids using genetic structure analysis and confirmed the reliability of our assignments via simulations. Then, we identified candidate loci for species maintenance with three complementary tests for selection and obtained partial gene sequences linked to an outlier locus and three other loci. We detected evidence of recent hybridization among all species and considerable gene flow between Q. ellipsoidalis and Q. velutina. Overall, c. 20% of Q. velutina had recent ancestry from Q. ellipsoidalis, whereas nearly 30% of Q. ellipsoidalis had a Q. velutina ancestor. Most loci were negligibly to weakly differentiated among species, but two gene‐linked microsatellites deviated significantly from neutral expectations in multiple, complementary outlier tests. Both outlier loci were located in the same 15‐cM bin on an existing Q. robur linkage map, a region under divergent selection in other oak species. Adaptive loci in this highly differentiated genomic region may contribute to ecological divergence among species and limit introgression.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced‐complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over‐represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini ‐ the largest radiation of mimetic butterflies.     

Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double‐digest restriction‐site‐associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well‐developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California.     

Patterns of Admixture and Population Structure in Native Populations of Northwest North America

The initial contact of European populations with indigenous populations of the Americas produced diverse admixture processes across North, Central, and South America. Recent studies have examined the genetic structure of indigenous populations of Latin America and the Caribbean and their admixed descendants, reporting on the genomic impact of the history of admixture with colonizing populations of European and African ancestry. However, relatively little genomic research has been conducted on admixture in indigenous North American populations. In this study, we analyze genomic data at 475,109 single-nucleotide polymorphisms sampled in indigenous peoples of the Pacific Northwest in British Columbia and Southeast Alaska, populations with a well-documented history of contact with European and Asian traders, fishermen, and contract laborers. We find that the indigenous populations of the Pacific Northwest have higher gene diversity than Latin American indigenous populations. Among the Pacific Northwest populations, interior groups provide more evidence for East Asian admixture, whereas coastal groups have higher levels of European admixture. In contrast with many Latin American indigenous populations, the variance of admixture is high in each of the Pacific Northwest indigenous populations, as expected for recent and ongoing admixture processes. The results reveal some similarities but notable differences between admixture patterns in the Pacific Northwest and those in Latin America, contributing to a more detailed understanding of the genomic consequences of European colonization events throughout the Americas.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Isolation and characterization of 10 polymorphic microsatellites in saltcedars (Tamarix chinensis and Tamarix ramosissima)

Tamarix ramosissima and Tamarix chinensis are invasive weed species in western North America. Previous studies based on single locus DNA sequence data revealed some information about the invasion process, but multilocus markers can provide additional information about levels of introgression and genotype origins. We have developed primers that amplify 10 polymorphic microsatellite loci from T. ramosissima; these primer pairs also successfully amplify polymorphic microsatellites from the closely related T. chinensis, a species that forms hybrids with T. ramosissima in the western USA.     

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves,Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early‐breeding house sparrows and hybrids and rural late‐breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z‐chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land‐use changes in a mosaic landscape.     

Wolves in the Great Lakes region: a phylogeographic puzzle

Empirical studies demonstrate that natural hybridization in animals is more common than thought so far ( Mallet 2005 ), particularly among species that originated recently through cycles of population contraction–expansion arising from climate changes over the last glacial period, the Pleistocene. In addition, the post‐glacial global growth of human populations has fostered anthropogenic hybridization events, mediated by habitat changes, the persecution of large predators and the introduction of alien species ( Allendorf et al. 2001 ). The Canis lineage shows cases of both natural and anthropogenic hybridization, exacerbating the controversy about the number of species that should be formally validated in the taxonomic lists, the evolutionary role of genetic introgression and the ways to manage hybrids with invading wild or domesticated populations. The study by Wheeldon et al. (2010) , published in this issue of Molecular Ecology, adds a new piece to the intricate puzzle of evolution and taxonomy of Canis in North America. They show that sympatric wolves (C. lupus) and coyotes (C. latrans) are not (extensively) hybridizing in the western North American Great Lakes region (GLR). Widespread hybridization between coyotes and a genetically distinct, but closely related, wolf‐like population (the eastern wolf) occurred in the northeastern regions of North America. In Wheeldon et al.’s (2010) opinion, these data should prove definitely that two different species of wolf (the western gray wolf C. lupus and the eastern wolf C. lycaon) and their hybrids are distributed across the GLR.     

Genetic analysis of admixture and patterns of introgression in foundation cottonwood trees (Salicaceae) in southwestern Colorado,USA

Cottonwoods are well known as foundation riparian trees that support diverse communities and drive ecosystem processes. Although hybridization naturally occurs when the distributions of two or more cottonwood species overlap, few cottonwood hybrid zones have been genetically characterized. We use genetic and genomic analyses to characterize patterns of admixture and introgression for a newly described hybrid zone at the intersection of three species (Populus L. Salicaceae—Populus deltoides, Populus fremontii, and Populus angustifolia) in southwestern Colorado, USA. Analysis of nuclear and chloroplast microsatellite marker data detected substantial genetic variation among individuals, revealing that (1) hybridization is occurring between two, not three, species (P. deltoides and P. angustifolia); (2) gene flow is bidirectional; (3) hybrids are not abundant (admixture detected in only 34 of 270 trees), with most being early-generation F₁ hybrids; (4) cytonuclear disequilibria exists and F₁ hybrids tend to retain P. deltoides—like chloroplasts; and (5) roughly 30 % of the nuclear markers deviated from a neutral pattern of introgression, suggesting that selection may play a role in shaping the genetic structure of the hybrid zone in this region. Overall, our results show that despite strong selection maintaining species divergence, transfer of allelic variants across species boundaries can occur. Our study assesses the fine-scale genetic structure of hybridization between P. angustifolia and P. deltoides and lays the foundation for examining how geographic differences in hybrid zone dynamics arise and may influence subsequent ecological and evolutionary processes.     

Spatiotemporally consistent genomic signatures of reproductive isolation in a moving hybrid zone

Studies of hybrid zone dynamics often investigate a single sampling period and draw conclusions from that temporal snapshot. Stochasticity can, however, result in loci with spurious outlier patterns, which is exacerbated by limited temporal or geographic sampling. Comparing admixed populations from different geographic regions is one way to detect repeatedly divergent genomic regions potentially involved in reproductive isolation. Temporal comparisons also allow us to control partially for the role of stochasticity, but the power of temporal sampling has not yet been adequately explored. In North America, black‐capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadees hybridize in a contact zone extending from New Jersey to Kansas. The hybrid zone is likely maintained by strong intrinsic selection against hybrids, and it is moving north. We used a reduced representation genomic approach and temporally spaced sampling—two samples of ～80 individuals separated by a decade—to determine the pattern and consistency of selection and genomic introgression in the chickadee hybrid zone. We report consistently low introgression for highly divergent loci between P. atricapillus and P. carolinensis in this moving hybrid zone. This is strong evidence that these loci may be linked to genomic regions involved in reproductive isolation between chickadees.     

Widespread hybridization and bidirectional introgression in sympatric species of coral reef fish

Coral reefs are highly diverse ecosystems, where numerous closely related species often coexist. How new species arise and are maintained in these high geneflow environments have been long‐standing conundrums. Hybridization and patterns of introgression between sympatric species provide a unique insight into the mechanisms of speciation and the maintenance of species boundaries. In this study, we investigate the extent of hybridization between two closely related species of coral reef fish: the common coral trout (Plectropomus leopardus) and the bar‐cheek coral trout (Plectropomus maculatus). Using a complementary set of 25 microsatellite loci, we distinguish pure genotype classes from first‐ and later‐generation hybrids, identifying 124 interspecific hybrids from a collection of 2,991 coral trout sampled in inshore and mid‐shelf reefs of the southern Great Barrier Reef. Hybrids were ubiquitous among reefs, fertile and spanned multiple generations suggesting both ecological and evolutionary processes are acting to maintain species barriers. We elaborate on these finding to investigate the extent of genomic introgression and admixture from 2,271 SNP loci recovered from a ddRAD library of pure and hybrid individuals. An analysis of genomic clines on recovered loci indicates that 261 SNP loci deviate from a model of neutral introgression, of which 132 indicate a pattern of introgression consistent with selection favouring both hybrid and parental genotypes. Our findings indicate genome‐wide, bidirectional introgression between two sympatric species of coral reef fishes and provide further support to a growing body of evidence for the role of hybridization in the evolution of coral reef fishes.     

Climate-influenced boreotropical survival and rampant introgressions explain the thriving of New World grapes in the north temperate zone

The north temperate region was characterized by a warm climate and a rich thermophilic flora before the Eocene, but early diversifications of the temperate biome under global climate change and biome shift remain uncertain. Moreover, it is becoming clear that hybridization/introgression is an important driving force of speciation in plant diversity. Here, we applied analyses from biogeography and phylogenetic networks to account for both introgression and incomplete lineage sorting based on geno...