Cline Centres Research Articles

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages.

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

Are there hybrid zones in Fagus sylvatica L. sensu lato?

The complex of beech taxa in western Eurasia (Fagus sylvatica L. s.l.) represents one of the most widespread broad-leaved trees in the temperate forests of the northern hemisphere. Two subspecies are generally recognized, F. sylvatica ssp. sylvatica in Europe and F. sylvatica ssp. orientalis in southwestern Asia, but several strongly differentiated regional populations were identified within this complex. To reveal geographical and genomic patterns of hybridization and introgression, genetic variation at 6789 single-nucleotide polymorphisms (SNP) resulting from ddRAD sequencing was studied in 288 specimens representing 207 populations distributed across the whole range. Two clusters resulting from a Bayesian analysis of population Structure well coincided with the subspecies, while increasing the number of clusters allowed to distinguish divergent lineages within subspecies. Two zones of potential hybridization between such lineages were analyzed: western Balkans at the contact of both subspecies and the Alpine area at the contact of populations coming from different glacial refugia. Geographical distribution of hybrid indices in the Balkan hybrid zone indicates asymmetric gene flow preferentially directed from F. orientalis to F. sylvatica. Eight SNPs showed a significant shift of genomic cline centers, seven of them toward excess F. orientalis ancestry, suggesting that introgression irregularly affects different regions of the beech genome. On the other hand, no genomic cline outliers were found in the Alpine hybrid zone. Implications for assisted gene flow are shortly discussed.

Hybrid zone or hybrid lineage: a genomic reevaluation of Sibley's classic species conundrum in Pipilo towhees.

Hybrid zones can be studied by modeling clines of trait variation (e.g., morphology, genetics) over a linear transect. Yet, hybrid zones can also be spatially complex, can shift over time, and can even lead to the formation of hybrid lineages with the right combination of dispersal and vicariance. We reassessed Sibley's (1950) gradient between Collared Towhee (Pipilo ocai) and Spotted Towhee (Pipilo maculatus) in Central Mexico to test whether it conformed to a typical tension-zone cline model. By comparing historical and modern data, we found that cline centers for genetic and phenotypic traits have not shifted over the course of 70 years. This equilibrium suggests that secondary contact between these species, which originally diverged over 2 million years ago, likely dates to the Pleistocene. Given the amount of mtDNA divergence, parental ends of the cline have very low autosomal nuclear differentiation (FST = 0.12). Dramatic and coincident cline shifts in mtDNA and throat color suggest the possibility of sexual selection as a factor in differential introgression, while a contrasting cline shift in green back color hints at a role for natural selection. Supporting the idea of a continuum between clinal variation and hybrid lineage formation, the towhee gradient can be analyzed as one population under isolation-by-distance, as a two-population cline, and as three lineages experiencing divergence with gene flow. In the middle of the gradient, a hybrid lineage has become partly isolated, likely due to forested habitat shrinking and fragmenting as it moved upslope after the last glacial maximum and a stark environmental transition. This towhee system offers a window into the potential outcomes of hybridization across a dynamic landscape including the creation of novel genomic and phenotypic combinations and incipient hybrid lineages.

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

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

Locus-specific introgression in young hybrid swarms: Drift may dominate selection.

Closely related species that have previously inhabited geographically separated ranges are hybridizing at an increasing rate due to human disruptions. These human-mediated hybrid zones can be used to study reproductive isolation between species at secondary contact, including examining locus-specific rates of introgression. Introgression is expected to be heterogenous across the genome, reflecting variation in selection. Those loci that introgress especially slowly are good candidates for being involved in reproductive isolation, while those loci that introgress quickly may be involved in adaptive introgression. In the context of conservation, policy makers are especially concerned about introduced alleles moving quickly into the background of a native or endemic species, as these alleles could replace the native alleles in the population, leading to extinction via hybridization. We applied genomic cline analyses to 44,997 SNPs to identify loci introgressing more or less when compared to the genome wide expectation in a human-mediated hybridizing population of red deer and sika in Kintyre Scotland. We found 11.4% of SNPs had cline centres that were significantly different from the genome wide expectation, and 17.6% of all SNPs had excess rates of introgression. Based on simulations, we believe that many of these markers have diverged from the genome-wide average due to drift, rather than because of selection, and we suggest that these simulations can be useful as a null distribution for future studies of genomic clines. Future work on red deer and sika could determine the policy implications of allelic-replacement due to drift rather than selection, and could use replicate, geographically distinct hybrid zones to narrow down those loci that are responding to selection.

Selection and hybridization shaped the rapid spread of African honey bee ancestry in the Americas.

Recent biological invasions offer 'natural' laboratories to understand the genetics and ecology of adaptation, hybridization, and range limits. One of the most impressive and well-documented biological invasions of the 20th century began in 1957 when Apis mellifera scutellata honey bees swarmed out of managed experimental colonies in Brazil. This newly-imported subspecies, native to southern and eastern Africa, both hybridized with and out-competed previously-introduced European honey bee subspecies. Populations of scutellata-European hybrid honey bees rapidly expanded and spread across much of the Americas in less than 50 years. We use broad geographic sampling and whole genome sequencing of over 300 bees to map the distribution of scutellata ancestry where the northern and southern invasions have presently stalled, forming replicated hybrid zones with European bee populations in California and Argentina. California is much farther from Brazil, yet these hybrid zones occur at very similar latitudes, consistent with the invasion having reached a climate barrier. At these range limits, we observe genome-wide clines for scutellata ancestry, and parallel clines for wing length that span hundreds of kilometers, supporting a smooth transition from climates favoring scutellata-European hybrid bees to climates where they cannot survive winter. We find no large effect loci maintaining exceptionally steep ancestry transitions. Instead, we find most individual loci have concordant ancestry clines across South America, with a build-up of somewhat steeper clines in regions of the genome with low recombination rates, consistent with many loci of small effect contributing to climate-associated fitness trade-offs. Additionally, we find no substantial reductions in genetic diversity associated with rapid expansions nor complete dropout of scutellata ancestry at any individual loci on either continent, which suggests that the competitive fitness advantage of scutellata ancestry at lower latitudes has a polygenic basis and that scutellata-European hybrid bees maintained large population sizes during their invasion. To test for parallel selection across continents, we develop a null model that accounts for drift in ancestry frequencies during the rapid expansion. We identify several peaks within a larger genomic region where selection has pushed scutellata ancestry to high frequency hundreds of kilometers past the present cline centers in both North and South America and that may underlie high-fitness traits driving the invasion.

Postzygotic isolation drives genomic speciation between highly cryptic Hypocnemis antbirds from Amazonia.

How species evolve reproductive isolation in the species-rich Amazon basin is poorly understood in vertebrates. Here, we sequenced a reference genome and used a genome-wide sample of SNPs to analyze a hybrid zone between two highly cryptic species of Hypocnemis warbling-antbirds-the Rondonia warbling-antbird (H. ochrogyna) and Spix's warbling-antbird (H. striata)-in a headwater region of southern Amazonia. We found that both species commonly hybridize, producing F1 s and a variety of backcrosses with each species but we detected only one F2 -like hybrid. Patterns of heterozygosity, hybrid index, and interchromosomal linkage disequilibrium in hybrid populations closely match expectations under strong postzygotic isolation. Hybrid zone width (15.4 km) was much narrower than expected (211 km) indicating strong selection against hybrids. A remarkably high degree of concordance in cline centers and widths across loci, and a lack of reduced interspecific Fst between populations close to versus far from the contact zone, suggest that genetic incompatibilities have rendered most of the genome immune to introgression. These results support intrinsic postzygotic isolation as a driver of speciation in a moderately young cryptic species pair from the Amazon and suggest that species richness of the Amazon may be grossly underestimated.

Social and economic factors responsible for environmental performance: A global analysis.

Several factors influence the environmental performance simultaneously but ethnic fractionalization, political freedom, financial development and institutional quality have a substantial impact to explain the environmental performance across economies. This study focuses to explore that how environmental performance is affected by these economic, political and social indicators by using the annual data of 163 developed and developing countries covering the time period of 1996-2016.The data is collected from World Development Indicators, World Governance Indicators, Freedom House and Cline Centre. The stationarity of variables is analyzed through LLC, IPS and ADF Fisher Chi-square test. Before applying panel ARDL approach to find out the long run relationship among variables, order of integration is determined through Pedroni’s cointegration test. The findings of study highlight that ethnic diversity; institutional quality and political freedom play a significant role to decrease CO2emissions while energy consumption, GDP growth and financial development are increasing the environmental degradation. Ethnic diversity is a source of creative and innovative approaches about problem solving of environmental degradation. Political freedom allows people to participate in decision making that posits much compliance with environmental agreements. The foreign direct investment is attracted by good quality institutions which cause to advent of more environment friendly technology along with attractions for further innovations that may helpful to reduce CO2 emissions. Contrary, GDP growth, financial development and energy consumption enhance the industrialization and urbanization leading to increase the level of CO2 emissions. It is suggested for policy makers that cohesion among different ethnic groups; improving institutional quality; providing political freedom to people and inclusive financial sector will acknowledge the less pollutant environment.

Cross-decades stability of an avian hybrid zone.

Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long-term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S.occidentalis) in the Washington Cascades was previously inferred to be moving from northern S.townsendi southwards towards S.occidentalis, based on plumage and behavioural patterns as well as a 2000-km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2-3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987-94 and 2015-16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a 'tension zone' model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2-3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species.

Tracing the footprints of a moving hybrid zone under a demographic history of speciation with gene flow.

A lack of optimal gene combinations, as well as low levels of genetic diversity, is often associated with the formation of species range margins. Conservation efforts rely on predictive modelling using abiotic variables and assessments of genetic diversity to determine target species and populations for controlled breeding, germplasm conservation and assisted migration. Biotic factors such as interspecific competition and hybridization, however, are largely ignored, despite their prevalence across diverse taxa and their role as key evolutionary forces. Hybridization between species with well‐developed barriers to reproductive isolation often results in the production of offspring with lower fitness. Generation of novel allelic combinations through hybridization, however, can also generate positive fitness consequences. Despite this possibility, hybridization‐mediated introgression is often considered a threat to biodiversity as it can blur species boundaries. The contribution of hybridization towards increasing genetic diversity of populations at range margins has only recently gathered attention in conservation studies. We assessed the extent to which hybridization contributes towards range dynamics by tracking spatio‐temporal changes in the central location of a hybrid zone between two recently diverged species of pines: Pinus strobiformis and P. flexilis. By comparing geographic cline centre estimates for global admixture coefficient with morphological traits associated with reproductive output, we demonstrate a northward shift in the hybrid zone. Using a combination of spatially explicit, individual‐based simulations and linkage disequilibrium variance partitioning, we note a significant contribution of adaptive introgression towards this northward movement, despite the potential for differences in regional population size to aid hybrid zone movement. Overall, our study demonstrates that hybridization between recently diverged species can increase genetic diversity and generate novel allelic combinations. These novel combinations may allow range margin populations to track favourable climatic conditions or facilitate adaptive evolution to ongoing and future climate change.

Testing an hypothesis of hybrid zone movement for toads in France.

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B.spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B.spinosus to B.bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B.spinosus and B.bufo gene pools at the southern (B.spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Sex chromosome inversions enforce reproductive isolation across an avian hybrid zone.

Across hybrid zones, the sex chromosomes are often more strongly differentiated than the autosomes. This is regularly attributed to the greater frequency of reproductive incompatibilities accumulating on sex chromosomes and their exposure in the heterogametic sex. Working within an avian hybrid zone, we explore the possibility that chromosome inversions differentially accumulate on the Z chromosome compared to the autosomes and thereby contribute to Z chromosome differentiation. We analyse the northern Australian hybrid zone between two subspecies of the long-tailed finch (Poephila acuticauda), first described based on differences in bill colour, using reduced-representation genomic sequencing for 293 individuals over a 1,530-km transect. Autosomal differentiation between subspecies is minimal. In contrast, 75% of the Z chromosome is highly differentiated and shows a steep genomic cline, which is displaced 350km to the west of the cline in bill colour. Differentiation is associated with two or more putative chromosomal inversions, each predominating in one subspecies. If inversions reduce recombination between hybrid incompatibilities, they are selectively favoured and should therefore accumulate in hybrid zones. We argue that this predisposes inversions to differentially accumulate on the Z chromosome. One genomic region affecting bill colour is on the Z, but the main candidates are on chromosome 8. This and the displacement of the bill colour and Z chromosome cline centres suggest that bill colour has not strongly contributed to inversion accumulation. Based on cline width, however, the Z chromosome and bill colour both contribute to reproductive isolation established between this pair of subspecies.

Dobzhansky-Muller incompatibilities, dominance drive, and sex-chromosome introgression at secondary contact zones: A simulation study.

Dobzhansky-Muller (DM) incompatibilities involving sex chromosomes have been proposed to account for Haldane's rule (lowered fitness among hybrid offspring of the heterogametic sex) as well as Darwin's corollary (asymmetric fitness costs with respect to the direction of the cross). We performed simulation studies of a hybrid zone to investigate the effects of different types of DM incompatibilities on cline widths and positions of sex-linked markers. From our simulations, X-Y incompatibilities generate steep clines for both X-linked and Y-linked markers; random effects may produce strong noise in cline center positions when migration is high relative to fitness costs, but X- and Y-centers always coincide strictly. X-autosome and Y-autosome incompatibilities also generate steep clines, but systematic shifts in cline centers occur when migration is high relative to selection, as a result of a dominance drive linked to Darwin's corollary. Interestingly, sex-linked genes always show farther introgression than the associated autosomal genes. We discuss ways of disentangling the potentially confounding effects of sex biases in migration, we compare our results to those of a few documented contact zones, and we stress the need to study independent replicates of the same contact zone.

Morphologically cryptic Amazonian bird species pairs exhibit strong postzygotic reproductive isolation.

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

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies.

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome-wide divergence. Variation in only one trait-head plumage patterning-was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome-wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.

Evidence for concerted movement of nuclear and mitochondrial clines in a lizard hybrid zone.

Moving hybrid zones provide compelling examples of evolution in action, yet long-term studies that test the assumptions of hybrid zone stability are rare. Using replicated transect samples collected over a 10-year interval from 2002 to 2012, we find evidence for concerted movement of genetic clines in a plateau fence lizard hybrid zone (Sceloporus tristichus) in Arizona. Cline-fitting analyses of SNP and mtDNA data both provide evidence that the hybrid zone shifted northward by approximately 2km during the 10-year interval. For each sampling period, the mtDNA cline centre is displaced from the SNP cline centre and maintaining an introgression distance of approximately 3km. The northward expansion of juniper trees into the Little Colorado River Basin in the early 1900s provides a plausible mechanism for hybrid zone formation and movement, and a broadscale quantification of recent land cover change provides support for increased woody species encroachment at the southern end of the hybrid zone. However, population processes can also contribute to hybrid zone movement, and the current stability of the ecotone habitats in the centre of the hybrid zone suggests that movement could decelerate in the future.

Evidence for a bimodal distribution of hybrid indices in a hybrid zone with high admixture.

The genetic structure of a hybrid zone can provide insights into the relative roles of the various factors that maintain the zone. Here, we use a multilocus approach to characterize a hybrid zone between two subspecies of killifish (Fundulus heteroclitus, Walbaum 1792) found along the Atlantic coast of North America. We first analysed clinal variation along the Atlantic coast using a single-nucleotide polymorphism in the mitochondrial DNA (mtDNA) displacement loop (D-loop) and a panel of nine nuclear microsatellite markers. A model constraining all clines to the same width and centre was not significantly different from a model in which the clines were allowed to vary independently. Locus-by-locus analysis indicated that the majority of nuclear clines shared the same centre as the mtDNA cline, and the widths of these clines were also narrower than that predicted by a neutral model, suggesting that selection is operating to maintain the hybrid zone. However, two of the nuclear clines had widths greater than the neutral prediction and had centres that were displaced relative to the mtDNA cline centre. We also found that a marsh located near the centre of the mtDNA cline demonstrated a bimodal distribution of nuclear hybrid index values, suggesting a deficit of first-generation hybrids and backcrossed genotypes. Thus, selection against hybrid genotypes may be playing a role in maintaining this hybrid zone and the associated steep nuclear and mtDNA clines.

Differences in gene flow in a twofold secondary contact zone of pond turtles in southern Italy (Testudines: Emydidae:Emys orbicularis galloitalica,E. o. hellenica,E. trinacris)

Using virtually range‐wide sampling for three pond turtle taxa (Emys orbicularis galloitalica,E. o. hellenica,E. trinacris), we analyse gene flow across their southern Italian contact zone. Based on population genetic analyses of 15 highly polymorphic microsatellite loci and a mitochondrial marker, we show that the general genetic pattern matches well with the current taxon delimitation. Yet, single individuals with conflicting genetic identity suggest translocation of turtles by humans. In addition, we identify in south‐western France and the vicinity of Rome populations being heavily impacted by introduced turtles. Cline analyses reveal that the major genetic break betweenE. o. galloitalicaandE. o. hellenicacorresponds well with the currently accepted intergradation zone in southern Italy. However, introgression is largely unidirectional fromE. o. galloitalicaintoE. o. hellenica. In the distribution range of the latter subspecies, genetic footprints ofE. o. galloitalicaare evident along most of the Italian east coast. Our results corroborate thatE. o. galloitalicawas introduced long ago in Corsica and Sardinia and naturalized there. Gene flow betweenE. orbicularisandE. trinacrisis negligible, with the Strait of Messina matching well with the narrow cline centre between the two species. This contrasts with other Mediterranean freshwater turtle species with extensive transoceanic gene flow. Compared to the two subspecies ofE. orbicularis, the SicilianE. trinacrisshows an unexpectedly strong population structuring, a finding also of some relevance for conservation. The differences between the two taxon pairsE. orbicularis/E. trinacrisandE. o. galloitalica/E. o. hellenicasupport their current taxonomic classification and make them attractive objects for follow‐up studies to elucidate the underlying mechanisms of speciation by comparing their properties.

Genetic and phenotypic characterization of a hybrid zone between polyandrous Northern and Wattled Jacanas in Western Panama.

BackgroundHybridization provides a unique perspective into the ecological, genetic and behavioral context of speciation. Hybridization is common in birds, but has not yet been reported among bird species with a simultaneously polyandrous mating system; a mating system where a single female defends a harem of males who provide nearly all parental care. Unlike simple polyandry, polyandrous mating is extremely rare in birds, with only 1% of bird species employing this mating system. Although it is classically held that females are “choosy” in avian hybrid systems, nearly-exclusive male parental care raises the possibility that female selection against heterospecific matings might be reduced compared to birds with other mating systems.ResultsWe describe a narrow hybrid zone in southwestern Panama between two polyandrous freshwater waders: Northern Jacana, Jacana spinosa and Wattled Jacana, J. jacana. We document coincident cline centers for three phenotypic traits, mtDNA, and one of two autosomal introns. Cline widths for these six markers varied from seven to 142 km, with mtDNA being the narrowest, and five of the six markers having widths less than 100 km. Cline tails were asymmetrical, with greater introgression of J. jacana traits extending westward into the range of J. spinosa. Likewise, within the hybrid zone, the average hybrid index of phenotypic hybrids was significantly biased towards J. spinosa. Species distribution models indicate that the hybrid zone is located at the edge of a roughly 100 km wide overlap where habitat is predicted to be suitable for both species, with more westerly areas suitable only for spinosa and eastward habitats suitable only for J. jacana.ConclusionThe two species of New World jacanas maintain a narrow, and persistent hybrid zone in western Panama. The hybrid zone may be maintained by the behavioral dominance of J. spinosa counterbalanced by unsuitable habitat for J. spinosa east of the contact zone. Although the two parental species are relatively young, mitochondrial cline width was extremely narrow. This result suggests strong selection against maternally-inherited markers, which may indicate either mitonuclear incompatibilities and/or female choice against heterospecific matings typical of avian hybrid systems, despite jacana sex role reversal.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0227-7) contains supplementary material, which is available to authorized users.