Characterizing patterns of introgressive hybridization between two species of Tyrannus following concurrent range expansion

Hybridization between two species typically occurs when allopatric or ecologically dissimilar species expand into areas of secondary contact or habitat transitions. However, as species' ranges shift rapidly in response to environmental change, the potential for novel types of ephemeral hybrid zones exists. Here, we document and describe the occurrence, prevalence and symmetry of a previously undocumented hybrid zone involving two sympatric, ecologically similar sister species that have been expanding their ranges eastward in the central USA: Scissor‐tailed Flycatchers Tyrannus forficatus and Western Kingbirds Tyrannus verticalis. We identified cases of hybridization and introgression using analyses of eight microsatellite DNA loci and a single mitochondrial gene. We also evaluated short‐term reproductive consequences of hybridization for both species by surveying for both species and potential hybrids at the periphery of their ranges in northeastern Arkansas and western Tennessee, USA. Genetic data revealed bi‐directional backcrossing at the periphery of the species’ ranges, including a cryptic hybrid. We also analysed DNA of putative ‘pure’ individuals from other parts of their ranges and detected two cryptic admixed individuals, suggesting backcrossed individuals from the periphery may be dispersing to breed or that hybridization events have occurred in the core. Finally, our results suggest that there are no short‐term reproductive consequences of hybridization for the two species. In total, hybrid zones that occur at the edges of expanding, sympatric ranges may be ephemeral; we suggest they play an important role in introgression and may have long‐standing impacts for sympatric sister species. Exploring the extent of hybrid zones such as this for other range‐expanding taxa will elucidate whether this type of hybrid zone is unique or a common occurrence.     

Hybridization between two sympatrically breeding species of fur seal at Iles Crozet revealed by genetic analysis

Introgressive hybridization occurs when closely related taxa overlap in distribution and is often associated with historically isolated populations coming into contact as a result of anthropogenic disturbance. There is evolutionary and conservation interest in detecting hybridization to determine its implications on future species composition, especially for threatened and recovering taxa such as subantarctic (Arctocephalus tropicalis) and Antarctic (A. gazella) fur seals, which were driven to the brink of extinction by human exploitation. Hybridization between these species has been reported at two locations and they breed sympatrically at a third site, Iles Crozet. While hybrid individuals have previously been identified based on phenotype, individuals can be difficult to classify based on these characteristics alone. Genotypic hybrid identification has been successful in several species, including fur seals. In this study we conducted an assignment test using microsatellite data to identify hybrids and to measure the frequency of hybridization at Iles Crozet. Samples were collected from 372 individuals and screened with 6 polymorphic microsatellite markers. MtDNA genotypes were also determined for individuals identified as hybrids or backcrosses based on microsatellite genotype. Phenotype, microsatellite and mtDNA genotype were then compared in order to identify hybrids. The results indicate that 1% of the population have hybrid genotypes and at a minimum, 2.4% of the population are backcrossed to parental species. We found that the two species are genetically distinct from one another and given the low rate of hybridization it is unlikely that they will fuse. These results suggest that there is a mechanism for species recognition that acts as a barrier to hybridization. It therefore seems unlikely that fur seals are threatened by significant introgression. Further investigation of fur seal mating systems would provide valuable insight into the mechanisms governing hybridization and species recognition in mate choice.     

Preliminary analysis of a hybrid zone between two subspecies of Zootermopsis nevadensis

For species largely allopatric in distribution, zones of contact provide an opportunity for hybridization, testing grounds for species boundaries, and may result in the formation of a new species. Thus, hybrid zones have the potential to provide important insights into speciation. In this study, we performed a preliminary analysis of a hybrid zone between two subspecies of the dampwood termite, Zootermopsis nevadensis (Z. n. nuttingi and Z. n. nevadensis) near Bartle, CA, using 12 microsatellite loci and a mitochondrial gene. Fifty-seven colonies collected in 36 locations were analyzed. The analysis of genetic markers revealed a large hybrid zone approximately 104 km in width. Although stepped clines best explained the data, we are unable to rule out the existence of a mosaic hybrid zone. We inferred weak selection (s* < 0.05%) against hybrids, but the data also suggested the existence of a barrier to gene flow from Z. n. nevadensis to Z. n. nuttingi, but not in the other direction. Given the large zone of contact, extensive sampling is needed to obtain a more complete characterization of the hybrid zone. However, the results of this study suggest that despite the accumulation of phenotypic and genetic differences and time since divergence (~2 million years), Z. n. nuttingi and Z. n. nevadensis are capable of extensive hybridization.     

Molecular patterns of introgression in a classic hybrid zone between the Australian tree frogs,Litoria ewingii and L. paraewingi: evidence of a tension zone

Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization.     

Window of opportunity: an episode of recruitment in a Banksia hybrid zone demonstrates continuing hybridization and phenotypic plasticity

Background and Aims

In perennial plants (especially post-fire resprouters), extant populations may reflect recruitment events in the distant past. This is true of hybrid zones formed by two Banksia species of swamps and woodlands in south-eastern Australia, Banksia robur and B. oblongifolia. Both resprout after fire but recruitment is dependent on periodic fires. Although plants of intermediate morphology have also been identified as hybrids using allozyme markers, the extent of ongoing hybridization is unknown. This study investigates whether both microsatellite markers and morphological measurements can be used to distinguish between the two species and their hybrids. A recent recruitment event and microsatellite markers allow the frequency of ongoing hybridization to be estimated, and also the effects of environmental variation on the morphology of plants and seedlings to be tested.

Methods

Variation at seven microsatellite loci was scored and seven leaf characteristics within putatively pure stands and mixed stands of both species were measured, revealing that the two species were genetically and morphologically distinct and that mixed stands also contained genetically and sometimes morphologically distinct hybrids. An opportunity created by wildfires was used to analyse the genetics and morphometrics of adults and seedlings from two hybrid zones.

Key Results

Approximately 9 % of adults and 21 % of seedlings were identified as genetic hybrids in both hybrid zones. Within these sites, the genotype of mature plants correlated well with morphology, except for some hybrid plants that had parental morphology. However, seedling morphology was highly variable and insufficient to describe the composition of the hybrid zone in this cohort. Greater phenotypic plasticity was evident among seedlings growing within the hybrid zones than seedlings growing in pots.

Conclusions

The hybrid zones are complex and the range of genotypes detected in seedlings reveals both continuing hybridization and introgression.     

Low levels of hybridization between sympatric Arctic char (Salvelinus alpinus) and Dolly Varden char (Salvelinus malma) highlights their genetic distinctiveness and ecological segregation

Understanding the extent of interspecific hybridization and how ecological segregation may influence hybridization requires comprehensively sampling different habitats over a range of life history stages. Arctic char (Salvelinus alpinus) and Dolly Varden (S. malma) are recently diverged salmonid fishes that come into contact in several areas of the North Pacific where they occasionally hybridize. To better quantify the degree of hybridization and ecological segregation between these taxa, we sampled over 700 fish from multiple lake (littoral and profundal) and stream sites in two large, interconnected southwestern Alaskan lakes. Individuals were genotyped at 12 microsatellite markers, and genetic admixture (Q) values generated through Bayesian‐based clustering revealed hybridization levels generally lower than reported in a previous study (<0.6% to 5% of samples classified as late‐generation hybrids). Dolly Varden and Arctic char tended to make different use of stream habitats with the latter apparently abandoning streams for lake habitats after 2–3 years of age. Our results support the distinct biological species status of Dolly Varden and Arctic char and suggest that ecological segregation may be an important factor limiting opportunities for hybridization and/or the ecological performance of hybrid char.     

Hybridization between <Emphasis Type="Italic">Schoenoplectus</Emphasis> sedges across Chesapeake Bay marshes

Hybridization is an evolutionary mechanism capable of enhancing adaptive potential, especially among species in fragmented or disturbed ecosystems like coastal marshes. In this study, we evaluated whether hybridization might influence adaptive responses in coastal marshes that are susceptible to the effects of global environmental change. To do so, we examined the extent and nature of hybridization between Schoenoplectus americanus and S. pungens, two ecologically dominant sedges in low-lying marshes across Chesapeake Bay (USA). Observed patterns of variation at genetically based morphological traits, cpDNA and nuclear microsatellite markers confirm that introgressive hybridization occurs between the two species. Comparisons of microsatellite and cpDNA profiles found that hybridization is reciprocal, although a disproportionate number of hybrids exhibit genomic asymmetries favoring S. americanus. AIC model selection consistently identified latitude as the strongest explanatory variable for the distribution of parental species, although discriminant analysis indicated that distributions also correspond to variation in environmental conditions. Discriminant analysis further indicated that ecological correlates of hybrid and S. americanus genotypes are similar, but not uniformly so. These findings indicate that the boundary between S. americanus and S. pungens is porous, and that hybridization could influence responses of one or both species to changing environmental regimes.     

From hybridization to introgression between two closely related sympatric ant species

Interspecific hybridization is becoming more frequent worldwide due to increasing global changes and translocations of organisms. For individual organisms, the most significant negative consequences are sterility or inviability of hybrid offspring. However, hybridization sometimes leads to fertile offspring, promoting introgression from one species into another. In such situations, hybridization can play a key role in evolution and speciation. Combining hypervariable DNA (microsatellites) and mitochondrial DNA markers with the use of several modeling methods allow an efficient detection of hybridization processes. The present study therefore investigates hybridization between two ant species, Tetramorium immigrans and T. caespitum, using multiple methods, and systematically comparing results with simulated data to ensure accurate identification of hybrids. Introgression was revealed both by backcross detection based on 14 nuclear microsatellite loci and by mitochondrial‐nuclear discordance based on comparison with mitochondrial DNA (cytochrome c oxidase subunit I). Results were spatially consistent, with hybrids located at latitudes where parental species are sympatric. The causes and consequences of hybridization and introgression between T. caespitum and T. immigrans remain to be further investigated, especially because T. immigrans could be an invasive species in France.     

High gene flow despite opposite chirality in hybrid zones between enantiomorphic door snails

We studied differentiation and geneflow patterns between enantiomorphic door‐snail species in two hybrid zones in the Bucegi Mountains (Romania) to investigate the effects of intrinsic barriers (complications in copulation) and extrinsic selection by environmental factors. A mitochondrial gene tree confirmed the historical separation of the examined populations into the dextral Alopia livida and the sinistral Alopia straminicollis in accordance with the morphological classification, but also indicated gene flow between the species. By contrast, a network based on amplified fragment length polymorphisms (AFLP) markers revealed local groups of populations as units independent of their species affiliation. Admixture analyses based on AFLP data showed that the genomes of most individuals in the hybrid zones are composed of parts of the genomes of both parental taxa. The introgression patterns of a notable fraction of the examined markers deviated from neutral introgression. However, the patterns of most non‐neutral markers were not concordant between the two hybrid zones. There was also no concordance between non‐neutral markers in the two genomic clines and markers that were correlated with environmental variables or markers that were correlated with the proportion of dextral individuals in the populations. Neither extrinsic selection by environmental factors nor intrinsic barriers resulting from positive frequency‐dependent selection of the prevailing coiling direction were sufficient to maintain the distinctness of A. straminicollis and A. livida. Despite being historically separated units, we conclude that these taxa now merge where they come into contact.     

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

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise F_ST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.     

Unexpected hybridization patterns in Near Eastern terrapins (Mauremys caspica,M. rivulata) indicate ancient gene flow across the Fertile Crescent

Recent studies indicate that hybridization in animals occurs more frequently than previously thought and that it may play an important evolutionary role. Chelonians are capable of extensive hybridization, raising the question how chelonian species evolve and maintain genetic integrity despite hybridization. Here, we use two sister species with parapatric distribution, Mauremys caspica and M. rivulata, as our model. These taxa are estimated to have diverged some 5.3–7.0 million years ago. Using rangewide sampling and 13 unlinked polymorphic microsatellite markers, five nuclear loci and one mitochondrial gene, we show that hybridization is rare along the contact zone of the two species in Turkey. However, we discovered an unexpected hybrid swarm in the southern Levant that has been hitherto identified with M. rivulata. This hybrid swarm is separated from the inland species M. caspica by a 700‐km‐wide distribution gap corresponding to the Syrian Desert. Ecological palaeomodelling suggests that during more humid climatic episodes in the Last Glacial Maximum and mid‐Holocene, the current contact zone extended into the southern Levant, facilitating the establishment of the now isolated hybrid swarm. Our results support that there is not necessarily a general hybridization pattern in a given species couple and that the extent of gene flow may differ considerably in different parts of the distribution range. Moreover, our results highlight that studies on hybridization should not focus only on extant contact and hybrid zones, but should use rangewide sampling to detect signals of ancient hybridization that might otherwise be missed.     

Detection of hybridization between two loach species (<Emphasis Type="Italic">Paramisgurnus dabryanus</Emphasis> and <Emphasis Type="Italic">Misgurnus anguillicaudatus</Emphasis>) in wild populations

Artificial interspecific hybrids between large scale loach P. dabryanus and tetraploid pond loach M. anguillicaudatus (Cobitidae, Cypriniformes) are viable. To detect the occurrence of possible natural hybridization, genetic analyses by using microsatellite markers were performed for natural populations of large scale loach and pond loach, the reciprocal laboratory hybrids, and “supposed hybrids” with ambiguous morphology. The fertility of the artificial hybrids was also tested. At one diagnostic microsatellite (Mac50), one out of 20 “supposed hybrids” was identified to be F₁ hybrid between the two loach species because it had the same genotype as that of the laboratory hybrids. The triploid hybrids between the two species were confirmed to be female-sterile. The results show that rare hybridization has occurred between diploid large scale loach and tetraploid pond loach in nature although it may have little effect in genetic introgression. This study is helpful for fish conservation and encourages further investigation on natural hybridization and introgression of loaches.     

Characterization of a contemporaneous hybrid zone between two darter species (Etheostoma bison and E. caeruleum) in the Buffalo River System

Hybrid zones have long intrigued evolutionary biologists and provide a natural laboratory to explore the evolution of reproductive isolation (speciation). Molecular characterization of hybrid zone dynamics can provide insight into the strength of reproductive isolation as well as the underlying evolutionary processes shaping gene flow. Approximately one-third of darter species naturally hybridize making this species-rich North American freshwater teleost fish clade an ideal system to investigate the extent and direction of hybridization. The objective of this study was to use diagnostic microsatellite markers to calculate genetic hybrid index scores of two syntopic, but distantly related darter species, Etheostoma bison and Etheostoma caeruleum. A combination of hybrid index scores, assignment tests, and mitochondrial haplotype profiles uncovered mixed ancestry in approximately 6 % of sampled adult individuals, supporting contemporaneous hybridization that was previously undocumented in E. bison. Moreover, hybrids were not limited to the F₁ generation, but encompassed the entire suite of hybrid categories (F₁, F₂ and backcross hybrids). The low number of hybrids assigned to each hybrid category represents a bimodal hybrid zone, suggesting reproductive isolation is strong (but incomplete) and also advocates for the ability of hybrids to produce second-generation hybrids and backcross into both parental species, mediating introgression across species boundaries. To this end, cytonuclear profiles of the sampled parental species and hybrids were consistent with bidirectional gene flow, although there was an overall trend of asymmetric hybridization between E. caeruleum females and E. bison males. The spatiotemporal variation in hybridization rates and resulting cytonuclear patterns expanded on in this study provide a comparative genetic framework on which future studies can begin to elucidate the underlying processes that not only generate a mosaic hybrid zone, but maintain the distinctness of species in the face of gene flow.     

Low levels of hybridization across two contact zones among three species of woodpeckers (Sphyrapicus sapsuckers)

Three species of closely related woodpeckers (sapsuckers; Sphyrapicus) hybridize where they come into contact, presenting a rare ‘λ‐shape’ meeting of hybrid zones. Two of the three arms of this hybrid zone are located on either side of the Interior Plateau of British Columbia, Canada bordering the foothills of the Coast Mountains and the Rocky Mountains. The third arm is located in the eastern foothills of the Rocky Mountains. The zones of hybridization present high variability of phenotypes and alleles in relatively small areas and provide an opportunity to examine levels of reproductive isolation between the taxa involved. We examined phenotypes (morphometric traits and plumage) and genotypes of 175 live birds across the two hybrid zones. We used the Genotyping By Sequencing (GBS) method to identify 180 partially diagnostic single nucleotide polymorphisms (SNPs) to generate a genetic hybrid index (GHI) for each bird. Phenotypically diverged S. ruber and S. nuchalis are genetically closely related, while S. nuchalis and S. varius have similar plumage but are well separated at the genetic markers studied. The width of both hybrid zones is narrower than expected under neutrality, and analyses of both genotypes and phenotypes indicate that hybrids are rare in the hybrid zone. Rarity of hybrids indicates assortative mating and/or some form of fitness reduction in hybrids, which might maintain the species complex despite close genetic distance and introgression. These findings further support the treatment of the three taxa as distinct species.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Hybridization between two cryptic filamentous brown seaweeds along the shore: analysing pre‐ and postzygotic barriers in populations of individuals with varying ploidy levels

We aimed to study the importance of hybridization between two cryptic species of the genus Ectocarpus, a group of filamentous algae with haploid–diploid life cycles that include the principal genetic model organism for the brown algae. In haploid–diploid species, the genetic structure of the two phases of the life cycle can be analysed separately in natural populations. Such life cycles provide a unique opportunity to estimate the frequency of hybrid genotypes in diploid sporophytes and meiotic recombinant genotypes in haploid gametophytes allowing the effects of reproductive barriers preventing fertilization or preventing meiosis to be untangle. The level of hybridization between E. siliculosus and E. crouaniorum was quantified along the European coast. Clonal cultures (568 diploid, 336 haploid) isolated from field samples were genotyped using cytoplasmic and nuclear markers to estimate the frequency of hybrid genotypes in diploids and recombinant haploids. We identified admixed individuals using microsatellite loci, classical assignment methods and a newly developed Bayesian method (XPloidAssignment), which allows the analysis of populations that exhibit variations in ploidy level. Over all populations, the level of hybridization was estimated at 8.7%. Hybrids were exclusively observed in sympatric populations. More than 98% of hybrids were diploids (40% of which showed signs of aneuploidy) with a high frequency of rare alleles. The near absence of haploid recombinant hybrids demonstrates that the reproductive barriers are mostly postzygotic and suggests that abnormal chromosome segregation during meiosis following hybridization of species with different genome sizes could be a major cause of interspecific incompatibility in this system.     

Hybridization and introgression in two ecologically dissimilar Fundulus hybrid zones

Hybridization and introgression appear more common in rapidly evolving groups, suggesting an important role in the evolutionary process. Detailed studies of how extrinsic or intrinsic forces regulate hybridization and introgression have the potential for broadening our understanding of mechanisms generating diversity. Species in the Fundulus notatus species complex have broad overlapping ranges and occur in replicated hybrid zones along predictable stream gradients. Typical hybrid zone structure has Fundulus olivaceus in headwaters, F. notatus downstream, and hybrid zones near confluences or abrupt shifts in habitat. Rarely, the typical upstream‐downstream orientation is reversed raising questions as to how hybrid zones are formed and maintained. We used next‐generation sequencing data to study hybridization and introgression in hybrid zones in neighboring drainages that differ in orientation (typical and reversed). We predicted extrinsic forces linked to stream gradients would result in noticeable differences between the two. Contrary to predictions, the data indicate the hybrid zones are remarkably similar. We used individual‐based simulations to explore the potential role of intrinsic and extrinsic forces in generating and maintaining typical and reversed hybrid zones. Simulation results were consistent with reversed hybrid zones being formed from stochastic processes combined with strong intrinsic forces and weak extrinsic forces.     

Differential introgression across newt hybrid zones: Evidence from replicated transects

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small‐scale (~2 km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.     

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.     

Genetic monitoring of two decades of hybridization between allis shad (Alosa alosa) and twaite shad (Alosa fallax)

Habitat alteration has been implicated in driving hybridization between the sympatric migratory shads Alosa alosa and Alosa fallax. Morphological and molecular evidence is consistent with hybridization across the overlapping range of these species, but the temporal extent of hybrid occurrence and genetic consequences for populations have not been explored. Using eight nuclear microsatellite loci and samples collected between 1989 and 2008 in the Solway Firth (UK), we genetically identified hybrids, studied temporal changes in their frequency, and explored changes in allele frequencies of parental populations. These molecular data confirmed the hybrid status of individuals identified using morphology (number of rakers on the outer gill arch), and enabled separation of hybrids from purebred individuals. Mitochondrial cytochrome-b sequencing revealed the presence of two haplogroups, each predominantly occurring in one species. Heterospecific haplotypes were found in 22.3 and 12.8% of A. alosa and A. fallax individuals, respectively, consistent with backcrossing and suggesting that hybrids are fertile. On average, microsatellite-identified hybrids comprised 12.7% of all samples, but when individuals with cytonuclear discordance were also considered introgressed on average 25.4% of individuals were of hybrid ancestry. Overall, allelic richness remained largely unchanged within species, but there were declines in the inbreeding coefficient (F _IS) of both species and episodes of significant temporal allelic frequency change. Hybrids sampled between 2004 and 2008 showed no evidence of lower fecundity relative to purebred individuals. Together, results suggest that hybridization between shad species in northern Europe is prevalent, and has been ongoing over at least two decades. The challenge is now to understand the extent to which observed patterns are linked to immigration from other populations, and the mechanisms that have prevented species collapse despite apparent hybrid fertility and longstanding introgression of neutral markers.