A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata).

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7-10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Tadpoles of hybridising fire-bellied toads (B. bombina and B. variegata) differ in their susceptibility to predation.

The role of adaptive divergence in the formation of new species has been the subject of much recent debate. The most direct evidence comes from traits that can be shown to have diverged under natural selection and that now contribute to reproductive isolation. Here, we investigate differential adaptation of two fire-bellied toads (Anura, Bombinatoridae) to two types of aquatic habitat. Bombina bombina and B. variegata are two anciently diverged taxa that now reproduce in predator-rich ponds and ephemeral aquatic sites, respectively. Nevertheless, they hybridise extensively wherever their distribution ranges adjoin. We show in laboratory experiments that, as expected, B. variegata tadpoles are at relatively greater risk of predation from dragonfly larvae, even when they display a predator-induced phenotype. These tadpoles spent relatively more time swimming and so prompted more attacks from the visually hunting predators. We argue in the discussion that genomic regions linked to high activity in B. variegata should be barred from introgression into the B. bombina gene pool and thus contribute to gene flow barriers that keep the two taxa from merging into one.

Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time.

Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds (B. bombina) or puddles (B. variegata), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata): Inference from transcriptome-wide coalescence analyses.

Ancient origins, profound ecological divergence, and extensive hybridization make the fire-bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small-scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes.

Local dynamics and dispersal in a structured population of the whirligig beetle Deneutus assimilis.

The study illustrates the ecological determinants and evolutionary consequences of dispersal in the pond-living water beetle Dineutus assimilis (Coleoptera: Gyrinidae). Over 2 years, local populatiopn dynamics were studied in 51 ponds within a 60-km2 study area. In most of the 31 occupied ponds, and even in large populations, abundances changed dramatically from one year to the next. Nine extinction and nine colonisation events were observed. These temporal patterns show no sign of spatial autocorrelation. Such a habitat distribution should favour high dispersal rates. Indeed, D. assimilis was found to be a very effective coloniser of newly available sites (mean propagule size: 23). A mark-recapture study showed that most dispersal occurred after diapause and over distances ranging from 100 m to at least 20 km. Yet despite frequent movement, the local variability in environmental conditions maintiins a large variance in average reproductive success per pond. Furthermore, immigration rates vary widely within a season. The apparent lack of correlation between these two sources of variation should greatly strengthen the role of drift in this system. A companion paper (Nürnberger and Harrison 1995) documents a non-random distribution of mitochondrial haplotypes due to recent population bottlenecks.

SPATIAL POPULATION STRUCTURE IN THE WHIRLIGIG BEETLE DINEUTUS ASSIMILIS: EVOLUTIONARY INFERENCES BASED ON MITOCHONDRIAL DNA AND FIELD DATA.

The spatial population structure of the pond-living water beetle Dineutus assimilis (Coleoptera: Gyrinidae) was investigated through a field study of population dynamics and dispersal, with a concurrent assessment of the spatial distribution of mitochondrial DNA (mtDNA) restriction-fragment-length polymorphism (RFLP). A comprehensive 2-yr survey within a 60-km2 study area revealed pronounced fluctuations in local abundances, including extinctions and colonizations. The recapture of marked individuals showed that dispersal among ponds is frequent in both males and females and connects populations on a large geographic scale (maximum observed flight distance: 20 km). The population structure of D. assimilis is thus characterized by both pronounced genetic drift and frequent gene flow. Together, these two forces generate a pattern of very local and transient genetic differentiation. Mitochondrial DNA samples collected within a few kilometers indicate highly significant spatial structure, if newly founded demes or those that experienced recent bottlenecks are included. These results based on four demes within the study area were placed into a regional context by further samples collected at distances of 100 km and 200 km. Fst estimates computed on increasing spatial scales were variable but showed no increasing trend. Thus, gene flow exerts a strong homogenizing force over a wide geographic range but is counteracted locally by genetic drift. These findings highlight the need to supplement estimates of Fst with additional data to arrive at valid interpretations of the genetic information. More generally, this study raises questions about how to capture the relevant features of dynamic, subdivided populations to understand their evolutionary dynamics.

NATURAL SELECTION ON QUANTITATIVE TRAITS IN THE BOMBINA HYBRID ZONE.

Observations on the means, variances, and covariances of quantitative traits across hybrid zones can give information similar to that from Mendelian markers. In addition, they can identify particular traits through which the cline is maintained. We describe a survey of six traits across the hybrid zone between Bombina bombina and Bombina variegata (Amphibia: Discoglossidae) near Pescenica in Croatia. We obtained laboratory measuments of the belly pattern, skin thickness, mating call, skeletal form, egg size, and the developmental time of tadpoles. Although offspring from hybrid populations showed no evidence of reduced viability, a third of the F1 families failed completely, irrespective of the direction of the cross. All traits differed significantly between the taxa. Clines in belly pattern, skin thickness, mating call, and skeletal form were closely concordant with clines in four diagnostic enzyme loci. However, the cline in developmental time was displaced towards bombina, and the cline in egg size was displaced towards variegata. This discordance could be because the traits are not inherited additively or because they are subject to different selection pressures. We favor the latter explanation in the case of developmental time. We show that moderate selection acting directly on a trait suffices to shift its position; rather stronger selection is needed to change its width appreciably. Within hybrid populations, there are significant associations among quantitative traits, and between traits and enzymes. Phenotypic variances also increase in hybrid populations. These observations can be explained by linkage disequilibria among the underlying loci. However, the average magnitude of the covariance between traits is about half that expected from the linkage disequilibria between enzyme loci. The discrepancy is not readily explained by nonadditive gene action. This puzzle is now unresolved and calls for further investigation.

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 Pescenica, 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 (FIS ) 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 FIS 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.