Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Hybrid zones”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Clustering of contact zones, hybrid zones, and phylogeographic breaks in North America.

A recent test for the existence of suture zones in North America, based on hybrid zones studied since 1970, found support for only two of the 13 suture zones identified by Remington in 1968 (Swenson and Howard 2004). One limitation of that recent study was the relatively small number of hybrid zones available for mapping. In this study, we search for evidence of clustering of contact zones between closely related taxa using data not only from hybrid zones but from species range maps of trees, birds, and mammals and from the position of phylogeographic breaks within species. Digital geographic range maps and a geographic information system approach allowed for accurate and rapid mapping of distributional data. Areas of contact between closely related species and phylogeographic breaks within species clustered into areas characterized by common physiographic features or predicted by previously hypothesized glacial refugia. The results underscore the general importance of geographic barriers to dispersal (mountain chains) and climate change (periods of cooling alternating with periods of warming, which lead to the contraction and expansion of species ranges) in species evolution.

Geography↗

Patterns of Genomic Divergence and Introgression in Two Primulina Hybrid Zones.

Hybrid zones have long been promoted as natural laboratories for understanding the mechanisms of speciation. Multiple or replicated hybrid zones are particularly informative, as they allow for assessing the consistency of genomic divergence and introgression across different environmental contexts and demographic histories, thereby improving our understanding of the factors that drive or hinder speciation on a broader scale. Here, using whole-genome resequencing data, we compare the patterns of genomic divergence and introgression in two Primulina hybrid zones. We found that genomic divergence in both hybrid zones is largely shaped by neutral processes, with only a few genomic regions showing signatures of balancing or lineage-specific selection. Genomic cline analyses identified numerous SNPs that showed significantly steeper clines and biased centres than the genome-wide expectation in both hybrid zones, consistent with the existence of reproductive barriers. Within regions of restricted gene flow, we identified 21 genes shared between the two hybrid zones. Annotation of gene function revealed that several genes are involved in reproductive processes. In addition, many zone-specific outlier loci were linked to genes associated with pollen and flower development, suggesting that these barriers may contribute to reproductive isolation under localised ecological conditions. Overall, these findings suggest that while certain reproductive barriers remain consistent across independent hybrid zones, others may be contingent on local environmental contexts. Our results demonstrate that both general and zone-specific mechanisms contribute to reproductive isolation in Primulina, providing empirical evidence that some genomic barriers recur across independent hybrid zones while others arise through localised adaptation.

Lamiales↗

Estimates of selection and gene flow from measures of cline width and linkage disequilibrium in heliconius hybrid zones.

Hybrid zones can yield estimates of natural selection and gene flow. The width of a cline in gene frequency is approximately proportional to gene flow (sigma) divided by the square root of per-locus selection (square root of s). Gene flow also causes gametic correlations (linkage disequilibria) between genes that differ across hybrid zones. Correlations are stronger when the hybrid zone is narrow, and rise to a maximum roughly equal to s. Thus cline width and gametic correlations combine to give estimates of gene flow and selection. These indirect measures of sigma and s are especially useful because they can be made from collections, and require no field experiments. The method was applied to hybrid zones between color pattern races in a pair of Peruvian Heliconius butterfly species. The species are Müllerian mimics of one another, and both show the same changes in warning color pattern across their respective hybrid zones. The expectations of cline width and gametic correlation were generated using simulations of clines stabilized by strong frequency-dependent selection. In the hybrid zone in Heliconius erato, clines at three major color pattern loci were between 8.5 and 10.2 km wide, and the pairwise gametic correlations peaked at R approximately 0.35. These measures suggest that s approximately 0.23 per locus, and that sigma approximately 2.6 km. In erato, the shapes of the clines agreed with that expected on the basis of dominance. Heliconius melpomene has a nearly coincident hybrid zone. In this species, cline widths at four major color pattern loci varied between 11.7 and 13.4 km. Pairwise gametic correlations peaked near R approximately 1.00 for tightly linked genes, and at R approximately 0.40 for unlinked genes, giving s approximately 0.25 per locus and sigma approximately 3.7 km. In melpomene, cline shapes did not perfectly fit theoretical shapes based on dominance; this deviation might be explained by long-distance migration and/or strong epistasis. Compared with erato, sample sizes in melpomene are lower and the genetics of its color patterns are less well understood. In spite of these problems, selection and gene flow are clearly of the same order of magnitude in the two species. The relatively high per locus selection coefficients agree with "major gene" theories for the evolution of Müllerian mimicry, but the genetic architecture of the color patterns does not. These results show that the genetics and evolution of mimicry are still only sketchily understood.

Animals↗

Morphological and molecular evidence for hybridization and introgression in a willow (Salix) hybrid zone.

Hybrid zones provide biologists with the opportunity to examine genetic and ecological interactions between differentiated populations. Accurate identification of hybrid genealogies is considered a necessary prerequisite to understanding observed patterns of hybridization-related phenomena. We analysed molecular and morphological data from individuals in a hybrid zone between two species of willows (Salix sericea Marshall and S. eriocephala Michaux) and report the use of randomly amplified polymorphic DNA (RAPD), chloroplast DNA (cpDNA), and ribosomal DNA (rDNA) markers, as well as vegetative morphology and foliar chemistry data to identify individuals in terms of hybrid genealogy and to infer the direction and extent of backcrossing and introgression within the hybrid zone. A novel version of a maximum likelihood estimate approach (developed for this study) was used to calculate hybrid index scores from RAPD marker data; this method produced results similar to those obtained using traditional arithmetic methods. Distribution of rDNA, cpDNA, and chemistry data were examined within the graphical context of RAPD-based hybrid index score histograms and principal component analyses (PCA) on RAPD and morphology data. Seven of the 21 plants classified as S. eriocephala in the field were possible introgressants. Another plant presented an unequivocal example of backcrossed S. sericea chemistry and RAPD markers. Inter- and intraspecific chloroplast diversity found within the hybrid zone suggests both historic introgression (perhaps in a glacial refugium), and contemporary hybridization. Patterns of inheritance and expression within the hybrid zone suggest that morphological characters are often not expressed in a simple additive fashion, and problems associated with both morphological and molecular data are considered.

Crosses, Genetic↗

Clines maintained by overdominant selection in hybrid zones.

Hybrid zone models often consider environment-independent selection to operate against all hybrids. However, empirical studies suggest that hybrids may be as fit or fitter than the hybridizing parental taxa in some environments. In this study we develop a novel mathematical model to explore the effects of one form of hybrid superiority on the genetic structure of hybrid zones. Our primary goals were to investigate the allele frequency clines at a locus experiencing overdominant selection and at a linked neutral or underdominant locus. Our results indicate that overdominant selection results in flat equilibrium allele frequency clines throughout the hybrid zone and an excess of heterozygosity relative to Hardy-Weinberg equilibrium. However, the genetic clines at linked neutral or underdominant loci tend not to reflect this overdominance even when the loci are tightly linked. Overall, we conclude that overdominance is unlikely to be detected in genetic surveys unless many loci are assayed.

Alleles↗

Ecophysiology of first and second generation hybrids in a natural plant hybrid zone.

Hybrids between related species vary widely in relative fitness, and that fitness can depend upon the environment. We investigated aspects of physiology that might influence fitness patterns in a plant hybrid zone. Seeds of Ipomopsis aggregata, I. tenuituba, F1 hybrids, F2 hybrids, and offspring of crosses between natural hybrids were planted into the relatively mesic site of origin for I. aggregata and the drier site for natural hybrids. We measured rates of photosynthesis (Amax), transpiration (E), instantaneous (A/E) and long-term (delta13C) indices of water use efficiency (WUE), and leaf nitrogen and carbon. We also examined correlations of these traits with plant size. Photosynthetic rate and A/E were higher in vegetative than flowering plants. WUE varied between sites and years, but differences among genotypic classes were spatially and temporally consistent. Instantaneous WUE was higher for F1 hybrids than for the average of the parental species, thereby showing heterosis. There was no evidence of hybrid breakdown, as WUE was no different in the F2 than the average across the F1 and parental species. Nor did WUE depend on cross direction in producing F1 progeny. Carbon isotope discrimination revealed higher long-term water use efficiency in I. tenuituba than I. aggregata. Leaf nitrogen was higher in I. tenuituba than I. aggregata, and higher in offspring of natural hybrids than in the F2. Results indicate heterosis for water use efficiency, with no hybrid breakdown. Heterosis in WUE may help to explain the relatively high survival of both reciprocal F1 hybrids in dry sites within the natural hybrid zone.

Analysis of Variance↗

Cytonuclear disequilibria in a spatially structured hybrid zone.

A hybrid zone model was developed that uses a stepping-stone model of population structure along with both nuclear and cytoplasmic genotypes to evaluate the effect that migration has on random mating organisms when no selection is present. Numerical simulations indicate that a number of different allele frequency and cytonuclear disequilibrium patterns can be found across the hybrid zone when it has reached equilibrium, and these results are discussed in the context of relative migration rates of the two sexes and the two source populations. The importance of numbers of subpopulations sampled, census time, and the persistence of pure species individuals into the hybrid zone are also discussed. Although this model does not consider selection or assortative mating in the hybrid zone, it should provide a useful baseline for evaluating joint cytonuclear genetic data from a structurally complex hybrid zone.

Animals↗

Spatial structure and habitat variation in a grasshopper hybrid zone.

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

Animals↗

Reproductive division of labor between hybrid and nonhybrid offspring in a fire ant hybrid zone.

Interspecific hybridization can often impose a substantial fitness cost due to reduced hybrid viability or fecundity. In social insects, however, such costs disproportionately impact reproductive offspring, whereas hybrids who become sterile workers can be functional, and even beneficial, colony members. Genomic imprinting of the paternal genome in reproductive, but not worker female offspring has been proposed as a mechanism to avoid genomic incompatibilities in hybrid queens in a hybrid zone between two fire ant species, Solenopsis geminata and S. xyloni. A study of allozyme variation demonstrated differences between the worker caste displaying a hybrid phenotype, and the winged queen caste displaying only the mother's phenotype. In this study, we investigate whether these differences are caused by genomic imprinting or genetic differences between castes by comparing variability of proteins to that of microsatellite markers. Workers and winged queens differed genetically at both classes of marker, indicating that allozyme differences were caused by underlying genetic differences between castes rather than differences in gene expression due to imprinting. Workers were F1 S. geminata x S. xyloni hybrids, whereas nearly all winged queens were of pure S. xyloni ancestry. Thus, S. xyloni within the hybrid zone appears to have evolved social hybridogenesis, in which the loss of worker potential in pure-species offspring necessitates hybridization for worker production, but prevents hybrids from being represented in the reproductive caste.

Animals↗

AdmixLD: fast genome-scale inference of ancestry disequilibrium in hybrid zones.

SUMMARY: Hybrid zones represent powerful natural systems for studying reproductive isolation and speciation. One key genomic signature of genetic incompatibilities and epistatic interactions is linkage disequilibrium (LD)-non-random associations between loci from different lineage backgrounds, generated by selection against maladaptive allele combinations. However, admixture alone induces strong genome-wide LD in hybrid populations, obscuring selection-driven signals. Here, we present AdmixLD, a fast, scalable C++ tool for genome-wide LD scanning in hybrid zones that estimates LD using partial correlation to control for individual hybrid index. By removing admixture-driven covariance, AdmixLD enhances detection of locus-specific associations and enables genome-scale identification of candidate barrier loci and interacting genomic regions. AVAILABILITY: The software and its code source are available at https://github.com/yzfranci/AdmixLD, and scripts for the data analysis are available at https://github.com/yzfranci/AdmixLDAnalysis.

Linkage Disequilibrium↗

Morphological variation within a macaque hybrid zone.

A hybrid zone exists between Macaca tonkeana and Macaca hecki (Primates: Cercopithecidae), centered along the Tawaeli-Toboli road in the narrow isthmus that connects North and Central Sulawesi, Indonesia. The current study demonstrates morphological substructure from north to south across the hybrid zone. Macaques to the northwest of the Tawaeli-Toboli road more closely resemble M. hecki, and macaques to the southeast of the road resemble M. tonkeana. While morphology shifted for both males and females over a distance of 1,500-2,000 m, adult males were significantly more M. tonkeana-like across the morphological gradient. This suggests that in the study area, males of M. tonkeana-like morphology are dispersing into hybrid groups at the expense of M. hecki-like males. A permutation analysis of diagnostic character states indicated that associations existed among several morphological traits. This could be due to the operation of one or several nonexclusive evolutionary processes, including recent secondary contact, pleiotropic effects, physical linkage of loci, natural selection against hybrids, the influx of parental types, or assortative mating. Continued environmental perturbation associated with the Tawaeli-Toboli road is likely to be a significant factor in the future of the M. tonkeana/M. hecki hybrid interaction.

Animals↗

Patterns of male sterility in a grasshopper hybrid zone imply accumulation of hybrid incompatibilities without selection.

It is now widely accepted that post-zygotic reproductive isolation is the result of negative epistatic interactions between derived alleles fixed independently at different loci in diverging populations (the Dobzhansky-Muller model). What is less clear is the nature of the loci involved and whether the derived alleles increase in frequency through genetic drift, or as a result of natural or sexual selection. If incompatible alleles are fixed by selection, transient polymorphisms will be rare and clines for these alleles will be steep where divergent populations meet. If they evolve by drift, populations are expected to harbour substantial genetic variation in compatibility and alleles will introgress across hybrid zones once they recombine onto a genetic background with which they are compatible. Here we show that variation in male sterility in a naturally occurring Chorthippus parallelus grasshopper hybrid zone conforms to the neutral expectations. Asymmetrical clines for male sterility have long tails of introgression and populations distant from the zone centre show significant genetic variation for compatibility. Our data contrast with recent observations on 'speciation genes' that have diverged as a result of strong natural selection.

Animals↗

Local adaptation and species segregation in two mussel (Mytilus edulis x Mytilus trossulus) hybrid zones.

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and Mytilus galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; Mytilus edulis and Mytilus trossulus are sympatric and hybridize on both western and eastern coasts of the Atlantic Ocean. We review the dynamics of hybridization in these two hybrid zones and evaluate the role of local adaptation for maintaining species boundaries. In Scandinavia, hybridization and gene introgression is so extensive that no individuals with pure M. trossulus genotypes have been found. However, M. trossulus alleles are maintained at high frequencies in the extremely low salinity Baltic Sea for some allozyme genes. A synthesis of reciprocal transplantation experiments between different salinity regimes shows that unlinked Gpi and Pgm alleles change frequency following transplantation, such that post-transplantation allelic composition resembles native populations found in the same salinity. These experiments provide strong evidence for salinity adaptation at Gpi and Pgm (or genes linked to them). In the Canadian Maritimes, pure M. edulis and M. trossulus individuals are abundant, and limited data suggest that M. edulis predominates in low salinity and sheltered conditions, whereas M. trossulus are more abundant on the wave-exposed open coasts. We suggest that these conflicting patterns of species segregation are, in part, caused by local adaptation of Scandinavian M. trossulus to the extremely low salinity Baltic Sea environment.

Adaptation, Biological↗

The effect of epistasis on the structure of hybrid zones.

Within hybrid zones that are maintained by a balance between selection and dispersal, linkage disequilibrium is generated by the mixing of divergent populations. This linkage disequilibrium causes selection on each locus to act on all other loci, thereby steepening clines, and generating a barrier to gene flow. Diffusion models predict simple relations between the strength of linkage disequilibrium and the dispersal rate, sigma, and between the barrier to gene flow, B, and the reduction in mean fitness, W. The aim of this paper is to test the accuracy of these predictions by comparison with an exact deterministic model of unlinked loci (r = 0.5). Disruptive selection acts on the proportion of alleles from the parental populations (p,q): W = exp[-S(4pq)beta], such that the least fit genotype has fitness e-s. Where beta << 1, fitness is reduced for a wide range of intermediate genotypes; where beta >> 1, fitness is only reduced for those genotypes close to p = 0.5. Even with strong epistasis, linkage disequilibria are close to sigma 2p'ip'j/rij, where p'i, p'j are the gradients in allele frequency at loci i, j. The barrier to gene flow, which is reflected in the steepening of neutral clines, is given by [formula: see text] where r, the harmonic mean recombination rate between the neural and selected loci, is here 0.5. This is a close approximation for weak selection, but underestimates B for strong selection. The barrier is stronger for small beta, because hybrid fitness is then reduced over a wider range of p. The widths of the selected clines are harder to predict: though simple approximations are accurate for beta = 1, they become inaccurate for extreme beta because, then, fitness changes sharply with p. Estimates of gene number, made from neutral clines on the assumption that selection acts against heterozygotes, are accurate for weak selection when beta = 1; however, for strong selection, gene number is overestimated. For beta > 1, gene number is systematically overestimated and, conversely, when beta < 1, it is underestimated.

Alleles↗

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice.

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.

Animals↗

Studies of mitochondrial DNA, allozyme and morphometric variation in a house mouse hybrid zone.

An unusual chromosomal hybrid zone of the house mouse, Mus musculus domesticus, exists in Upper Valtellina, Northern Italy, consisting of four Robertsonian (Rb) races and the standard (all-acrocentric, or 2n = 40) race, all hybridizing freely within 10 km2. The hybrid zone in Valtellina provides an excellent opportunity to study the role of Rb fusions in reproductive isolation and speciation. This hybrid zone has already been well studied for the distribution of Rb fusions and the fertility of hybrids, but in order to understand the dynamics of the zone, a basic understanding of the origin and genetic similarity of the chromosomal races is necessary. This paper presents the results of three different methods of measuring genetic differentiation: multivariate analysis of morphological traits and analyses of allozyme variation and mitochondrial DNA sequences. The standard race is clearly distinguishable from the three Rb races by all three methods, but the Rb races are not distinguishable from one another. This provides strong evidence for our previous suggestions that the well-established Rb races in Valtellina are closely related, and that the standard race was introduced into the valley more recently from a distant source. The fact that the Rb races are indistinguishable is also consistent with our hypothesis that a within-village speciation event involving two of the races (Hauffe & Searle, 1992) was a recent occurrence. The low level of allozyme heterozygosity among the Rb races suggests that these populations are the products of at least one bottleneck. The present article substantially extends earlier studies and provides the first detailed morphological and molecular analysis of this complex hybrid zone.

Animals↗

Bounded hybrid superiority in an avian hybrid zone: effects of mate, diet, and habitat choice.

There has been considerable debate in the study of hybrid zones as to whether hybrids may be superior to parental types within the area of contact (bounded hybrid superiority). In birds, naturally occurring hybridization is relatively common, and hybridization within this group always involves mate choice. If hybrids are superior, females choosing heterospecific mates should be expected to show higher fitness under the conditions prevalent in the hybrid zone. Hybrid superiority under these circumstances would reduce reinforcement and thereby help to maintain the hybrid zone. To examine this issue, we studied reproductive performances of hybrids and parental species of gulls (Larus occidentalis and Larus glaucescens) at two colonies within a linear hybrid zone along the west coast of the United States. This hybrid zone contains predominantly gulls of intermediate phenotype. Previous studies indicated that hybrids were superior to one or both parental types, but provided no data on possible mechanisms that underlie this hybrid superiority. Using a hybrid index designed specifically for these species, we identified to phenotype more than 300 individuals associated with nests, including both individual males and females within 73 pairs in the central portion of the hybrid zone and 74 pairs in the northern portion of the hybrid zone. There was little evidence of assortative mating, and what little there was resulted solely because of pairings within intergrades. In the central hybrid zone, females paired with hybrid males produced larger clutches and hatched and fledged more chicks compared with females paired to western gull males. This was a result of heavy predation on eggs in sand habitat, where male western gulls established territories. In contrast, many hybrid males established territories in vegetated cover that was less vulnerable to predation. In the northern part of the hybrid zone, clutch size did not differ among pair categories, however, there were differences in hatching and fledging success, with females paired to hybrid males showing better success compared to females paired to glaucous-winged gull males. Hybrids showed better hatching and fledging success in the north because hybrids are more like western gulls than glaucous-winged gulls in foraging behavior, taking a higher percentage of fish in their diet, which enhances chick growth and survival. This is believed to be the first documentation of bounded hybrid superiority that delineates the mechanisms that underlie hybrid superiority.

Animals↗

A fine-scale spatial analysis of the mosaic hybrid zone between Gryllus firmus and Gryllus pennsylvanicus.

The pattern of character variation within a hybrid zone, the hybrid zone structure, has been used to infer the processes that maintain hybrid zones. Unfortunately it is difficult to infer process from structure alone because many different processes can produce the same pattern of character variation. Mosaic hybrid zones may be maintained by exogenous selection in a heterogeneous environment and/or endogenous selection against hybrid individuals; habitat preference, premating isolating barriers and/or fertility selection can also contribute. The spatial scale at which a hybrid zone is sampled affects its apparent structure; a hybrid zone may appear clinal at one scale and mosaic at another. Here, we sample the mosaic hybrid zone between two field crickets, Gryllus firmus and G. pennsylvanicus, at a scale that spans the boundaries between individual soil-habitat patches. From our analysis, we find that at fine scales, the mosaic hybrid zone resolves into a set of steep clines across patch boundaries. Both morphological and molecular traits exhibit sharp and generally concordant clines. However, clines for mitochondrial DNA and one anonymous nuclear marker are clearly displaced as a result of current hybridization or past introgression (the "ghost of hybridization past"). Thus, scale is important for the structure of this and probably other hybrid zones. The extremely sharp, concordant clines across patch boundaries indicate that the cricket hybrid zone is undoubtedly structured by selection. However, the detailed mechanisms responsible for the maintenance of the hybrid zone--whether endogenous selection against hybrids, exogenous selection by the environment, and/or behavioral preferences for mates or habitats--remain to be elucidated. Determining these mechanisms will depend on closer inspection of the organisms themselves and their interactions, as is the case for all hybrid zones.

Animals↗