Spreading introgression in the wake of a moving contact zone

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Assortative mating and epistatic mating-trait architecture induce complex movement of the crow hybrid zone

Hybrid zones provide a window into the evolutionary processes governing species divergence. Yet, the contribution of mate choice to the temporal and spatial stability of hybrid zones remains poorly explored. Here, we investigate the effects of assortative mating on hybrid-zone dynamics by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and gray-coated hooded crows. In the best-fit model, narrow clines of the two mating-trait loci were maintained by a moderate degree of assortative mating inducing pre- and postzygotic isolation via positive frequency-dependent selection. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favoring gray-coated phenotypes ∼ 1 200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of the genetic mating-trait architecture for spatio-temporal hybrid-zone dynamics.

Geographic patterns in morphometric and genetic variation for coyote populations with emphasis on southeastern coyotes

Prior to 1900, coyotes (Canis latrans) were restricted to the western and central regions of North America, but by the early 2000s, coyotes became ubiquitous throughout the eastern United States. Information regarding morphological and genetic structure of coyote populations in the southeastern United States is limited, and where data exist, they are rarely compared to those from other regions of North America. We assessed geographic patterns in morphology and genetics of coyotes with special consideration of coyotes in the southeastern United States. Mean body mass of coyote populations increased along a west-to-east gradient, with southeastern coyotes being intermediate to western and northeastern coyotes. Similarly, principal component analysis of body mass and linear body measurements suggested that southeastern coyotes were intermediate to western and northeastern coyotes in body size but exhibited shorter tails and ears from other populations. Genetic analyses indicated that southeastern coyotes represented a distinct genetic cluster that differentiated strongly from western and northeastern coyotes. We postulate that southeastern coyotes experienced lower immigration from western populations than did northeastern coyotes, and over time, genetically diverged from both western and northeastern populations. Coyotes colonizing eastern North America experienced different selective pressures than did stable populations in the core range, and we offer that the larger body size of eastern coyotes reflects an adaptation that improved dispersal capabilities of individuals in the expanding range.

Genetic analyses reveal cryptic introgression in secretive marsh bird populations

Hybridization is common in bird populations but can be challenging for management, especially if one of the two parent species is of greater conservation concern than the other. King rails (Rallus elegans) and clapper rails (R. crepitans) are two marsh bird species with similar morphologies, behaviors, and overlapping distributions. The two species are found along a salinity gradient with the king rail in freshwater marshes and the clapper in estuarine marshes. However, this separation is not absolute; they are occasionally sympatric, and there are reports of interbreeding. In Virginia, USA, both king and clapper rails are identified by the state as Species of Greater Conservation Need, although clappers are thought to be more abundant and king rails have a higher priority ranking. We used a mitochondrial DNA marker and 13 diagnostic nuclear single nucleotide polymorphisms (SNPs) to identify species, classify the degree of introgression, and explore the evolutionary history of introgression in two putative clapper rail focal populations along a salinity gradient in coastal Virginia. Genetic analyses revealed cryptic introgression with site-specific rates of admixture. We identified a pattern of introgression where clapper rail alleles predominate in brackish marshes. These results suggest clapper rails may be displacing king rails in Virginia coastal waterways, most likely as a result of ecological selection. As introgression can result in various outcomes from outbreeding depression to local adaptation, continued monitoring of these populations would allow further exploration of hybrid fitness and inform conservation management.

Sharp acoustic boundaries across an altitudinal avian hybrid zone despite asymmetric introgression

Birdsong is a sexually selected trait that could play an important evolutionary role when related taxa come into secondary contact. Many songbird species, however, learn their songs through copying one or more tutors, which complicates the evolutionary outcome of such contact. Two subspecies of a presumed vocal learner, the grey-breasted wood-wren (Henicorhina leucophrys), replace each other altitudinally across the western slope of the Ecuadorian Andes. These subspecies are morphologically very similar, but show striking differences in their song. We examined variation in acoustic traits and genetic composition across the altitudinal range covered by both subspecies and between two allopatric populations. The acoustic boundary between the subspecies was found to be highly abrupt across a narrow elevational range with virtually no evidence of song convergence. Mixed singing and use of hetero-subspecific song occurred in the contact zone and was biased towards the use of leucophrys song types. Hetero-subspecific song copying by hilaris and not by leucophrys reflected a previously found asymmetric pattern of response to song playback. Using amplified fragment length polymorphisms (AFLP) markers, we detected hybridization in the contact zone and asymmetric introgression in parapatric populations, with more leucophrys alleles present in hilaris populations than vice versa. This pattern may be a trail of introgression due to upslope displacement of leucophrys by hilaris. Our data suggest that song learning may impact speciation and hybridization in contrasting ways at different spatial scales: although learning may speed up population divergence in songs, thereby enhancing assortative mating and reducing gene flow, it may at a local level also lead to the copying of heterospecific songs, therefore allowing some level of hybridization and introgression.

Differential introgression and the maintenance of species boundaries in an advanced generation avian hybrid zone

Background

Evolutionary processes, including selection and differential fitness, shape the introgression of genetic material across a hybrid zone, resulting in the exchange of some genes but not others. Differential introgression of molecular or phenotypic markers can thus provide insight into factors contributing to reproductive isolation. We characterized patterns of genetic variation across a hybrid zone between two tidal marsh birds, Saltmarsh (Ammodramus caudacutus) and Nelson’s (A. nelsoni) sparrows (n = 286), and compared patterns of introgression among multiple genetic markers and phenotypic traits.

Results

Geographic and genomic cline analyses revealed variable patterns of introgression among marker types. Most markers exhibited gradual clines and indicated that introgression exceeds the spatial extent of the previously documented hybrid zone. We found steeper clines, indicating strong selection for loci associated with traits related to tidal marsh adaptations, including for a marker linked to a gene region associated with metabolic functions, including an osmotic regulatory pathway, as well as for a marker related to melanin-based pigmentation, supporting an adaptive role of darker plumage (salt marsh melanism) in tidal marshes. Narrow clines at mitochondrial and sex-linked markers also offer support for Haldane’s rule. We detected patterns of asymmetrical introgression toward A. caudacutus, which may be driven by differences in mating strategy or differences in population density between the two species.

Conclusions

Our findings offer insight into the dynamics of a hybrid zone traversing a unique environmental gradient and provide evidence for a role of ecological divergence in the maintenance of pure species boundaries despite ongoing gene flow.

Electronic supplementary material

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

Evolutionary consequences of climate-induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short- and long-term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short-term processes acting on reproductive barriers between species and (ii) long-term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate-induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.

Climate-induced range shifts and possible hybridisation consequences in insects

Many ectotherms have altered their geographic ranges in response to rising global temperatures. Current range shifts will likely increase the sympatry and hybridisation between recently diverged species. Here we predict future sympatric distributions and risk of hybridisation in seven Mediterranean ischnurid damselfly species (I. elegans, I. fountaineae, I. genei, I. graellsii, I. pumilio, I. saharensis and I. senegalensis). We used a maximum entropy modelling technique to predict future potential distribution under four different Global Circulation Models and a realistic emissions scenario of climate change. We carried out a comprehensive data compilation of reproductive isolation (habitat, temporal, sexual, mechanical and gametic) between the seven studied species. Combining the potential distribution and data of reproductive isolation at different instances (habitat, temporal, sexual, mechanical and gametic), we infer the risk of hybridisation in these insects. Our findings showed that all but I. graellsii will decrease in distributional extent and all species except I. senegalensis are predicted to have northern range shifts. Models of potential distribution predicted an increase of the likely overlapping ranges for 12 species combinations, out of a total of 42 combinations, 10 of which currently overlap. Moreover, the lack of complete reproductive isolation and the patterns of hybridisation detected between closely related ischnurids, could lead to local extinctions of native species if the hybrids or the introgressed colonising species become more successful.

Species distribution models contribute to determine the effect of climate and interspecific interactions in moving hybrid zones

Climate is a major factor delimiting species' distributions. However, biotic interactions may also be prominent in shaping geographical ranges, especially for parapatric species forming hybrid zones. Determining the relative effect of each factor and their interaction of the contact zone location has been difficult due to the lack of broad scale environmental data. Recent developments in species distribution modelling (SDM) now allow disentangling the relative contributions of climate and species' interactions in hybrid zones and their responses to future climate change. We investigated the moving hybrid zone between the breeding ranges of two parapatric passerines in Europe. We conducted SDMs representing the climatic conditions during the breeding season. Our results show a large mismatch between the realized and potential distributions of the two species, suggesting that interspecific interactions, not climate, account for the present location of the contact zone. The SDM scenarios show that the southerly distributed species, Hippolais polyglotta, might lose large parts of its southern distribution under climate change, but a similar gain of novel habitat along the hybrid zone seems unlikely, because interactions with the other species (H. icterina) constrain its range expansion. Thus, whenever biotic interactions limit range expansion, species may become 'trapped' if range loss due to climate change is faster than the movement of the contact zone. An increasing number of moving hybrid zones are being reported, but the proximate causes of movement often remain unclear. In a global context of climate change, we call for more interest in their interactions with climate change.

Asymmetric introgression in the horticultural living fossil cycas sect. Asiorientales using a genome-wide scanning approach

The Asian cycads are mostly allopatric, distributed in small population sizes. Hybridization between allopatric species provides clues in determining the mechanism of species divergence. Horticultural introduction provides the chance of interspecific gene flow between allopatric species. Two allopatrically eastern Asian Cycas sect. Asiorientales species, C. revoluta and C. taitungensis, which are widely distributed in Ryukyus and Fujian Province and endemic to Taiwan, respectively, were planted in eastern Taiwan for horticultural reason. Higher degrees of genetic admixture in cultivated samples than wild populations in both cycad species were detected based on multilocus scans by neutral AFLP markers. Furthermore, bidirectional but asymmetric introgression by horticultural introduction of C. revoluta is evidenced by the reanalyses of species associated loci, which are assumed to be diverged after species divergence. Partial loci introgressed from native cycad to the invaders were also detected at the loci of strong species association. Consistent results tested by all neutral loci, and the species-associated loci, specify the recent introgression from the paradox of sharing of ancestral polymorphisms. Phenomenon of introgression of cultivated cycads implies niche conservation among two geographic-isolated cycads, even though the habitats of the extant wild populations of two species are distinct.

Hybridization but no evidence for backcrossing and introgression in a sympatric population of great reed warblers and clamorous reed warblers

Hybridization is observed frequently in birds, but often it is not known whether the hybrids are fertile and if backcrossing occurs. The breeding ranges of the great reed warbler (Acrocephalus arundinaceus) and the clamorous reed warbler (A. stentoreus) overlap in southern Kazakhstan and a previous study has documented hybridization in a sympatric population. In the present study, we first present a large set of novel microsatellite loci isolated and characterised in great reed warblers. Secondly, we evaluate whether hybridization in the sympatric breeding population has been followed by backcrossing and introgression.We isolated 181 unique microsatellite loci in great reed warblers. Of 41 loci evaluated, 40 amplified and 30 were polymorphic. Bayesian clustering analyses based on genotype data from 23 autosomal loci recognised two well-defined genetic clusters corresponding to the two species. Individuals clustered to a very high extent to either of these clusters (admixture proportions ≥ 0.984) with the exception of four previously suggested arundinaceus-stentoreus hybrid birds that showed mixed ancestry (admixture proportions 0.495-0.619). Analyses of simulated hybrids and backcrossed individuals showed that the sampled birds do not correspond to first-fourth-generation backcrosses, and that fifth or higher generation backcrosses to a high extent resemble 'pure' birds at this set of markers.We conclude that these novel microsatellite loci provide a useful molecular resource for Acrocephalus warblers. The time to reach reproductive isolation is believed to be very long in birds, approximately 5 Myrs, and with an estimated divergence time of 2 Myrs between these warblers, some backcrossing and introgression could have been expected. However, there was no evidence for backcrossing and introgression suggesting that hybrids are either infertile or their progeny inviable. Very low levels of introgression cannot be excluded, which still may be an important factor as a source of new genetic variation.

Introgression and rapid species turnover in sympatric damselflies

Background

Studying contemporary hybridization increases our understanding of introgression, adaptation and, ultimately, speciation. The sister species Ischnura elegans and I. graellsii (Odonata: Coenagrionidae) are ecologically, morphologically and genetically similar and hybridize. Recently, I. elegans has colonized northern Spain, creating a broad sympatric region with I. graellsii. Here, we review the distribution of both species in Iberia and evaluate the degree of introgression of I. graellsii into I. elegans using six microsatellite markers (442 individuals from 26 populations) and five mitochondrial genes in sympatric and allopatric localities. Furthermore, we quantify the effect of hybridization on the frequencies of the genetically controlled colour polymorphism in females of both species.

Results

In a principal component analysis of the microsatellite data, the first two principal components summarised almost half (41%) of the total genetic variation. The first axis revealed a clear separation of I. graellsii and I. elegans populations, while the second axis separated I. elegans populations. Admixture analyses showed extensive hybridization and introgression in I. elegans populations, consistent with I. elegans backcrosses and occasional F₁-hybrids, suggesting hybridization is on-going. More specifically, approximately 58% of the 166 Spanish I. elegans individuals were assigned to the I. elegans backcross category, whereas not a single of those individuals was assigned to the backcross with I. graellsii. The mitochondrial genes held little genetic variation, and the most common haplotype was shared by the two species.

Conclusions

The results suggest rapid species turnover in sympatric regions in favour of I. elegans, corroborating previous findings that I. graellsii suffers a mating disadvantage in sympatry with I. elegans. Examination of morph frequency dynamics indicates that hybridization is likely to have important implications for the maintenance of multiple female morphs, in particular during the initial period of hybridization.

Amplified fragment length polymorphism and sequence analyses reveal massive gene introgression from the European fungal pathogen Heterobasidion annosum into its introduced congener H. irregulare

The paucity of fungal species known to be currently hybridizing has significantly hindered our understanding of the mechanisms driving gene introgression in these eukaryotic microbes. Here, we describe an area of hybridization and gene introgression between the invasive plant pathogen Heterobasidion irregulare (introduced from North America) and the native H. annosum in Italy. A STRUCTURE analysis of amplified fragment length polymorphism data for 267 individuals identified gene introgression in 8-42% of genotypes in the invasion area, depending on site. Data indicate that introgression is mostly occurring unilaterally from the native to the invasive species and is responsible for 5-45% of genomes in admixed individuals. Sequence analysis of 11 randomly selected and unlinked loci for 30 individuals identified introgression at every locus, thus confirming interspecific gene flow involves a large number of loci. In 37 cases, we documented movement of entire alleles between the two species, but in 7 cases, we also documented the creation of new alleles through intralocus recombination. Sequence analysis did not identify enrichment of either transcriptionally different nonsynonymous alleles or of transcriptionally identical synonymous alleles. These findings may suggest introgression is occurring randomly for extant alleles without an obvious enrichment process driven by selection. However, further studies are needed to ensure selection is not at work elsewhere in the genome.

Widespread introgression does not leak into allotopy in a broad sympatric zone

Species that overlap over a large part of their range and habitat requirements are challenging for the study of speciation and hybridization. In this respect, the study of broadscale introgressive hybridization has raised recent interest. Here we studied hybridization between two closely related amphibians Lissotriton helveticus and Lissotriton vulgaris that reproduce over a wide sympatric zone. We used mitochondrial and microsatellite markers on 1272 individuals in 37 sites over Europe to detect hybrids at the individual-level and to analyse Hardy-Weinberg and linkage disequilibria at the population-level. Morphological traits showed a strong bimodal distribution. Consistently, hybrid frequency was low (1.7%). We found asymmetric introgression with five times more hybrids in L. vulgaris than in L. helveticus, a pattern probably explained by an unequal effective population size in a study part wherein L. helveticus numerically predominates. Strikingly, significant levels of introgression were detected in 73% of sites shared by both species. Our study showed that introgression is widespread but remains confined to the sites where the two species reproduce at the same time. This pattern may explain why these species remain genetically distinct over a broad sympatric zone.

A hybrid zone revisited: molecular and morphological analysis of the maintenance, movement, and evolution of a Great Plains avian (Cardinalidae: Pheucticus) hybrid zone

Black-headed grosbeaks (Pheucticus melanocephalus) and rose-breasted grosbeaks (Pheucticus ludovicianus) are passerine bird species known to hybridize in the Great Plains of North America. Both extrinsic (environmental) and intrinsic factors (pre- and postzygotic reproductive isolation) have been credited for the generation and maintenance of the grosbeak hybrid zone, but little is known about the genetic characteristics of this hybrid zone. To investigate the stability and extent of the grosbeak hybrid zone, we constructed clines from both molecular sequence data (mtDNA, three autosomal intron loci, and one Z-linked locus) and morphological data (morphometric analyses and hybrid index scores) to determined zone centre and width. Hybrid zone centre and width were also determined for samples collected across the zone 40 years ago from morphological data. The present and past clines were compared and provided support for stability in hybrid zone location and width, and the evolutionary implications of this are discussed. Three models of hybrid zone maintenance were investigated to consider the influence of intrinsic and extrinsic factors on this zone. Our results suggest low hybrid frequencies, a stable zone location and narrow width, and reduced hybrid fitness over the past 40 years best categorize the grosbeak hybrid zone as a tension zone.

Complex hybridization dynamics between golden-winged and blue-winged warblers (Vermivora chrysoptera and Vermivora pinus) revealed by AFLP, microsatellite, intron and mtDNA markers

Blue-winged (Vermivora pinus) and golden-winged warblers (Vermivora chrysoptera) have an extensive mosaic hybrid zone in eastern North America. Over the past century, the general trajectory has been a rapid replacement of chrysoptera by pinus in a broad, northwardly moving area of contact. Previous mtDNA-based studies on these species' hybridization dynamics have yielded variable results: asymmetric and rapid introgression from pinus into chrysoptera in some areas and bidirectional maternal gene flow in others. To further explore the hybridization genetics of this otherwise well-studied complex, we surveyed variation in three nuclear DNA marker types--microsatellites, introns, and a panel of amplified fragment length polymorphisms (AFLPs)--with the goal of generating a multilocus assay of hybrid introgression. All markers were first tested on birds from phenotypically and mitochondrially pure parental-type populations from outside the hybrid zone. Searches for private alleles and assignment test approaches found no combination of microsatellite or intron markers that could separate the parental populations, but seven AFLP characters exhibited significant frequency differences among them. We then used the AFLP markers to examine the extent and pattern of introgression in a population where pinus-phenotype individuals have recently invaded a region that previously supported only a chrysoptera-phenotype population. Despite the low frequency of phenotypic hybrids at this location, the AFLP data suggest that almost a third of the phenotypically pure chrysoptera have introgressed genotypes, indicating the presence of substantial cryptic hybridization in the history of this species. The evidence for extensive cryptic introgression, combined with the lack of differentiation at other nuclear loci, cautions against hybrid assessments based on single markers or on phenotypic traits that are likely to be determined by a small number of loci. Considered in concert, these results from four classes of molecular markers indicate that pinus and chrysoptera are surprisingly weakly differentiated and that far fewer genetically 'pure' populations of chrysoptera may exist than previously assumed, two findings with broad implications for the conservation of this rapidly declining taxon.

CONTRASTING PATTERNS OF HYBRIDIZATION IN LARGE HOUSE SPIDERS (TEGENARIA ATRICA GROUP, AGELENIDAE)

The integrity of species is not fixed and may vary geographically. Here we investigate the geographic distributions and interactions of species in the Tegenaria atrica group (Araneae: Agelenidae). Detailed mapping of T. saeva and T. gigantea in England and Wales shows them to be broadly allopatric in southern England with a tightly defined, and possibly long-standing, narrow zone of parapatry in central southern England. In the north of England (Yorkshire), by contrast, the species are broadly sympatric as a result of recent range expansions. GIS techniques are used to map the species distributions and to quantify, we believe for the first time, the intimacy of interspecific interactions. The extent and nature of hybridization in these two areas is examined through regression and multivariate analyses of morphology. We show that the relative incidence of hybridization is much greater in Yorkshire than within the parapatric zone in the south. Clear patterns of asymmetric introgression are observed in both northern and southern England, with a greater impact of T. gigantea on T. saeva than vice versa. We find no sign of morphological reproductive character displacement at the zone of parapatry that might indicate reinforcement, although we cannot exclude more subtle effects, for example via cuticular pheromones. The integrity of these two species seems to be breaking down in northern England, a process that might gain momentum as the gene pools become more similar.