GENETIC ANALYSIS OF A HYBRID ZONE BETWEEN THE FIRE-BELLIED TOADS,BOMBINA BOMBINA AND B. VARIEGATA,NEAR CRACOW IN SOUTHERN POLAND

The fire-bellied toads Bombina bombina and B. variegata differ extensively in biochemistry, morphology, and behavior. We use a survey of five diagnostic enzyme loci across the hybrid zone near Cracow in Southern Poland to estimate the dispersal rate, selection pressures, and numbers of loci which maintain this zone. The enzyme clines coincide closely with each other and with morphological and mitochondrial DNA clines. Although the zone lies on a broad transition between environments suitable for bombina and variegata, the close concordance of diverse characters, together with increased aberrations and mortality in hybrids, suggest that the zone is maintained largely by selection against hybrids. There are strong “linkage disequilibria” between each pair of (unlinked) enzyme loci (R? = 0.129 [2-unit support limits: 0.119–0.139]). These are probably caused by gene flow into the zone, and they give an estimate of dispersal (σ = 890 [790–940] m gen^?½). The clines are sharply stepped, with most of the change occurring within 6.15 (5.45–6.45) km, but with long tails of introgression on either side. This implies that the effective selection pressure on each enzyme marker (due largely to disequilibrium with other loci) is s* = 0.17 (0.159–0.181) at the center but that the selection acting directly on the enzyme loci is weak or zero (s_e < 0.0038). The stepped pattern implies a barrier to gene flow of 220 (48–415) km. This would substantially delay neutral introgression but would have little effect on advantageous alleles; the two taxa need not evolve independently. Strong selection is needed to maintain such a barrier: hybrid populations must have their mean fitness reduced by a factor of 0.65 (0.60–0.77). This selection must be spread over a large number of loci to account for the concordant patterns and the observed cline widths (N = 300 [80–2,000]).     

Variation in complex mating signals in an “island” hybrid zone between Stenobothrus grasshopper species

Two grasshopper species Stenobothrus rubicundus and S. clavatus were previously shown to meet in a narrow hybrid zone on Mount Tomaros in northern Greece. The species are remarkable for their complex courtship songs accompanied by conspicuous movements of antennae and wings. We analyzed variations in forewing morphology, antenna shape, and courtship song across the hybrid zone using a geographic information system, and we documented three contact zones on Mount Tomaros. All male traits and female wings show abrupt transitions across the contact zones, suggesting that these traits are driven by selection rather than by drift. Male clines in antennae are displaced toward S. clavatus, whereas all clines in wings are displaced toward S. rubicundus. We explain cline discordance as depending on sexual selection via female choice. The high covariance between wings and antennae found in the centers of all contact zones results from high levels of linkage disequilibria among the underlying loci, which in turn more likely results from assortative mating than from selection against hybrids. The covariance is found to be higher in clavatus‐like than rubicundus‐like populations, which implies asymmetric assortative mating in parental‐like sites of the hybrid zone and a movement of the hybrid zone in favor of S. clavatus.     

Hybridization between the fire-bellied toads Bombina bombina and Bombina variegata in the karst regions of Slovakia and Hungary: morphological and allozyme evidence

Geographic variability and genetic interactions in the contact zone between the fire-bellied toads, Bombina bombina and B. variegata, were studied using analysis of morphological and genetic variation in sixteen samples from the Slovak Karst and Aggtelek Karst regions. Genotype frequencies at four marker loci (Ldh-1, Mdh-1, Adk, Hem) demonstrate the existence of a hybrid zone with highly variable population structures. While some samples appear to represent panmictic hybrid populations, other samples are very heterogeneous. Pure individuals of both species occurred together with hybrids at one site. Habitat segregation among genetically differentiated demes probably causes this heterogeneity of population structures. Increased frequencies of the allele Ldh-1^M, which is present in low proportion in B. bombina south of the contact zone, were found in some hybrid populations.     

The mating call of hybrids of the fire-bellied toad and yellow-bellied toad (Bombina bombina (L.), Bombina v. variegata (L.), discoglossidae,anura)

Summary Male hybrids obtained by crossing Bombina bombina with Bombina v. variegata displayed calling behavior when treated with hormones at the age of about 12 months. The mating calls of the hybrids are intermediate with respect to rate, duration, inter-call interval, and fundamental frequency between the calls of the two purebred species.     

Eight new microsatellite loci for the critically endangered fire‐bellied toad Bombina bombina and their cross‐species applicability among anurans

We describe eight new microsatellite loci for the critically endangered fire‐bellied toad, Bombina bombina. Seven of them are polymorphic with two to seven alleles per locus, an expected heterozygosity between 0.41 and 0.8, and an observed heterozygosity between 0.27 and 0.7. The yield of new loci was relatively low, presumably due to mildly repetitive sequence motifs in microsatellite flanking regions. As typical for anurans, cross‐species amplification was limited (here, to congeners Bombina orientalis and Bombina variegata). Combining these new loci with those already available provides a reasonable number of loci for population studies and pedigree analysis in Bombina.     

VARIATION IN MATING CALL ACROSS THE HYBRID ZONE BETWEEN THE FIRE-BELLIED TOADS BOMBINA BOMBINA AND B. VARIEGATA

Three components of mating call (pulse duration, cycle length, and fundamental frequency) were measured and six diagnostic enzyme loci scored across the hybrid zone between the toads Bombina bombina and B. variegata. All three call components differ significantly, but only cycle length is diagnostic. The clines in call coincide with those for enzymes, and have similar widths. This suggests that there is no strong selection on any of these characters. There are significant correlations between electrophoretic markers and call components, but these are no stronger than would be expected if the electrophoretic loci and the genes causing mating call were neutral. The selection differential on the call is no greater than 6% of the difference in mean cycle length between the two taxa. There is a substantial increase in the variance of cycle length in the center of the zone, suggesting that a small number of loci are involved (≈ three). Recombination between these loci will hinder the evolution of reinforcement and may partly be responsible for the lack of premating isolation between B. bombina and B. variegata.     

Nine new tetranucleotide microsatellite markers for the fire‐bellied toad (Bombina bombina)

We describe nine new polymorphic tetranucleotide microsatellite loci isolated from the fire‐bellied toad (Bombina bombina). The relative yield of new loci was higher than described in previous studies in amphibians: out of 12 loci initially evaluated, nine were polymorphic and amplifying reliably. Number of alleles ranged from four to 10 and observed heterozygosities from 0.47 to 0.91. Seven loci were polymorphic also in Bombina variegata and five in Bombina orientalis. Enrichment protocols yielding long flanking regions potentially overcome difficulties (i.e, low yield of reliable loci relative to number of clones screened) which have been reported in microsatellite development in anurans.     

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.     

Inheritance and linkage analysis of five enzyme loci in interspecific hybrids of toadlets,genus Bombina

Progeny produced from Bombina bombina, B. variegata, and field-collected interspecific hybrids have been analyzed for the inheritance of five enzyme loci, which are fixed for alternate alleles in the parental species. Lactate dehydrogenase (Ldh-1), NAD-dependent malate dehydrogenase (Mdh-1), creatine kinase (Ck), adenylate kinase (Ak), and glucosephosphate isomerase (Gpi) are all inherited in a Mendelian manner as codominant alleles at nuclear loci. Both parental alleles are equally functional in artificial F₁ hybrids (female B. bombina×male B. variegata) at each of the loci studied. No linkage between any pair of loci was observed. Discovery of this inherited biochemical variation combined with a technique for assaying individual genotypes without killing the animals makes feasible studies of hybrid population structure heretofore impossible.This work was partly supported by the Polish Academy of Sciences within the project MR-II/6.     

HYBRID DYSFUNCTION IN FIRE-BELLIED TOADS (BOMBINA)

Reproductive isolation between two taxa may be due to endogenous selection, which is generated by incompatibilities between the respective genomes, to exogenous selection, which is generated by differential adaptations to alternative environments, or to both. The continuing debate over the relative importance of either mode of selection has highlighted the need for unambiguous data on the fitness of hybrid genotypes. The hybrid zone between the fire-bellied toad (Bombina bombina) and the yellow-bellied toad (B. variegata) in central Europe involves adaptation to different environments, but evidence of hybrid dysfunction is equivocal. In this study, we followed the development under laboratory conditions of naturally laid eggs collected from a transect across the Bombina hybrid zone in Croatia. Fitness was significantly reduced in hybrid populations: Egg batches from the center of the hybrid zone showed significantly higher embryonic and larval mortality and higher frequencies of morphological abnormalities relative to either parental type. Overall mortality from day of egg collection to three weeks after hatching reached 20% in central hybrid populations, compared to 2% in pure populations. There was no significant difference in fitness between two parental types. Within hybrid populations, there was considerable variation in fitness, with some genotypes showing no evidence of reduced viability. We discuss the implications of these findings for our understanding of barriers to gene flow between species.     

EVOLUTIONARY IMPLICATIONS OF DIVERGENT CLINES IN AN AVIAN (MANACUS: AVES) HYBRID ZONE

Abstract A previous study of the hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus) documented the unidirectional introgression of vitellinus male secondary sexual traits across the zone. Here, we examine the hybrid zone in greater genetic and morphological detail. Statistical comparisons of clines are performed using maximum‐likelihood and nonparametric bootstrap methods. Our results demonstrate that an array of six molecular and two morphometric markers agree in cline position and width. Clines for male collar and belly color are similar in width to the first eight clines, but are shifted in position by at least five cline widths. The result is that birds in intervening populations are genetically and morphometrically very like parental candei, but males have the plumage color of parental vitellinus. Neither neutral diffusion nor nonlinearity of color scales appear to be viable explanations for the large cline shifts. Genetic dominance of vitellinus plumage traits is another potential explanation that will require breeding experiments to test. Sexual selection remains a plausible explanation for the observed introgression of vitellinus color traits in these highly dimorphic, polygynous, lek‐mating birds. Two other clines, including a nondiagnostic isozyme locus, are similar in position to the main cluster of clines, but are broader in width. Thus, introgression at some loci is greater than that detected with diagnostic markers. Assuming that narrow clines are maintained by selection, variation in cline width indicates that selection is not uniform throughout the genome and that diagnostic markers are under more intense selective pressure. The traditional focus on diagnostic markers in studies of hybrid zones may therefore lead to underestimates of average introgression. This effect may be more pronounced in organisms with low levels of genetic divergence between hybridizing taxa.     

Inheritance of allozyme loci inBombina: One linkage group established

Two hybridizing European species of fire-bellied toads,Bombina bombina andB. variegata, have alternate electromorphs fixed at a number of allozyme loci. Segregation of alleles at seven allozyme loci (Ldh-1, Mdh-1, Ak, Ck, Gpi, Np, andEst-) was studied in a backross progeny of an F₁ interspecific hybrid male and aB. bombina female. Mendelian inheritance of allozyme forms at all seven loci was ascertained. Except for two loci,Gpi andEst-, which were found to be tightly linked (1 cM apart), other loci showed independent segregation.     

THE GENETIC STRUCTURE OF THE HYBRID ZONE BETWEEN THE FIRE-BELLIED TOADS BOMBINA BOMBINA AND B. VARIEGATA: COMPARISONS BETWEEN TRANSECTS AND BETWEEN LOCI

We compare the pattern of morphological and electrophoretic variation in the hybrid zone between Bombina bombina and B. variegata across two transects: one near Cracow and one 200 km away, near Przemy?l in southeastern Poland. Morphological variation across the Przemy?l transect had been surveyed more than 50 years ago; though we found a significant shift at one site, there is no evidence for gross movement over this period. Morphological and electrophoretic changes coincide, and the average shape of the clines is the same across both transects. At the center, most of the change in frequency of six diagnostic allozymes occurs within w = 6.05 km (2-unit support limits 5.56–6.54 km). These steep gradients are generated not by selection on the allozymes themselves, but by associations with other loci: though these markers are unlinked, they are in strong linkage disequilibrium with each other [R = D/√pquv = 0.22 (0.15–0.29) at the center]. Disequilibria are broken up as alleles diffuse away from the zone and flow into the new genetic background. The net barrier to the flow of genes from bombina into variegata, which is generated by these disequilibria, is B = 51 (22–81) km. The fitness of hybrids must be substantially reduced to produce such a barrier [W?_H/W?_P = 0.58 (0.54–0.68)], and this selection must be spread over many loci [N = 55 (26–88)]. Alleles introgress significantly less far than would be expected from the age of the zone and the estimated dispersal rate [σ = 0.99 (0.82–1.14) km gen.^1/2]: this implies selection of s_e = 0.37 (0.15–0.58)% on the enzymes themselves. There is weak but significant linkage disequilibrium well away from the center of the zone; this, together with the presence of parental and F₁ genotypes, suggests some long-range migration. However, such migration is not likely to cause significant introgression.     

HABITAT PREFERENCE IN THE BOMBINA HYBRID ZONE IN CROATIA

This paper demonstrates the effect of habitat heterogeneity and a habitat preference on the genetic structure of a hybrid zone between the toads Bombina bombina and B. variegata (Anura: Discoglossidae); 1613 toads from 85 sites across a transect near Pe??enica, Croatia, were scored for five unlinked diagnostic allozyme markers. These were found to be largely concordant. Aside from minor systematic deviations, there was little variance in allele frequency among loci within sites. Yet the allele frequencies did not follow a smooth cline, but formed a mosaic in the center, such that neighboring sites could differ markedly in their enzyme score. A detailed ecological survey revealed a correlation between this pattern and habitat. In keeping with the typical breeding sites of the parental taxa, B. bombina-like hybrids were found more often in ponds, whereas B. variegata-like hybrids were more common in puddles. In addition, there was significant heterozygote deficit (F_IS) and strong linkage disequilibrium (R), both of which were stronger on the B. bombina side of the transect, and stronger in puddles than ponds. Mark-recapture data showed: (1) that the animals disperse beyond the scale of the habitat pattern; (2) frequent turn-over of individuals within sites; and (3) nonrandom movement between two sites of different habitat type. We conclude that an active habitat preference must contribute to the observed association between marker alleles and habitat. As a consequence, there is incomplete mixing of the two gene pools, which could explain the high level of F_IS and R. The asymmetry in these parameters may reflect asymmetry in the preference or in the distribution of habitats across the zone. We discuss the implications of habitat preference for the dynamics of hybrid zones.     

Variation in complex courtship traits across a hybrid zone between grasshopper species of the Chorthippus albomarginatus group

Closely related grasshopper species of the Chorthippus albomarginatus group are notable for their extremely complex courtship songs, accompanied by a visual display. Two species of this group, Ch. albomarginatus and Ch. oschei, were previously shown to hybridize in a wide mosaic hybrid zone in Ukraine and Moldova. In this paper, variation in five courtship song characters, one character of visual display and the number of stridulatory pegs were analysed across the hybrid zone to estimate selection against hybrids and strength of assortative mating. Comparison of cline width and position across the hybrid zone showed concordant and coincident clines in four traits, such as three song characters and one morphological character, and discordant and non‐coincident clines in two other song characters and the character of visual display. Concordance of clines in different characters suggests an equal strength of selection acting on underlying loci. Increase of variance and covariance between phenotypic traits at the cline centre could more likely result from assortative mating than from selection against hybrids. Most pairwise cases showed the highest covariance for the oschei‐like, than for the albomarginatus‐like hybrid populations. This indicates that introgression of the oschei genes into the albomarginatus genome is stronger than vice versa, and may be evidence of the movement of the hybrid zone in favour of Ch. albomarginatus. Analysis of associations between phenotype and local vegetation showed that mosaic structure of the hybrid zone is explained to a great extent by habitat–phenotype associations. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 275–291.     

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,000 km 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 100 km 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.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

Phylogeographic analysis reveals northerly refugia for the riverine amphibian Triturus dobrogicus (Caudata: Salamandridae)

We investigated the recent evolutionary history of the Danube crested newt, Triturus dobrogicus through reconstructions of: (1) the number and position of refugia at the last glacial maximum, (2) the role of major central European rivers in pattern of post‐glacial dispersal, and (3) the present‐day distribution pattern. We analysed sequences of mitochondrial DNA (ND2, 1065 bp) and six microsatellite loci in 363 T. dobrogicus individuals from 58 populations covering the range of the species. Our analyses suggested that T. dobrogicus survived the last glacial maximum in two separate refugia positioned in northwestern Pannonia and in Southern Pannonia from where its range expanded along the Danube and Tisza Rivers. Our findings also confirmed that rivers played an important role in shaping the evolutionary history of amphibian species in Central Europe. We compared the T. dobrogicus range with another lowland amphibian, the fire‐bellied toad Bombina bombina, using species distribution modelling. In line with these models, the isolated mountains inside Pannonia are occupied not by B. bombina, but by its mountain‐dwelling sister‐species B. variegata. However, in contrast to the model, crested newts (including T. dobrogicus) are absent from these mountains. We attribute this biogeographical discrepancy to the positioning of the species’ refugia at the last glacial maximum.     

Steep,coincident, and concordant clines in mitochondrial and nuclear‐encoded genes in a hybrid zone between subspecies of Atlantic killifish,Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus‐specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non‐neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non‐neutral nuclear clines, SNPs in two nuclear‐encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin‐binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non‐unimodal, suggesting selection against advanced‐generation hybrids, possibly due to effects on processes involved in oxidative metabolism.