Brain size diversity in adaptation and speciation of subterranean mole rats

We have tested brain size diversity and encephalization in the actively speciating subterranean mole rats of the Spalax ehrenbergi superspecies in Israel. Our sample involved 171 individuals comprising 12 populations and 4 chromosomal species (2n = 52, 54, 58 and 60) distributed parapatrically from the northern Mediterranean region southward (2n = 52, 54→+58→60) into increasingly more arid and unpredictable climatic regimes, approaching the Negev Desert. Our results indicate that relative brain size and encephalization are highest in 2n = 60 as compared with 2n = 52, 54 and 58. We hypothesize that this pattern is adaptive and molded by natural selection. Brain evolution and higher encephalization in the S. ehrenbergi complex appears to be associated with increasing ecological stresses of aridity and climatic unpredictability.     

Sexual selection and natural selection in body size differentiation of subterranean mole rats

We exlored indirectly, the operation of sexual selection in subterranean mole rats of the Spalax ebrenbergi superspecies in Israel comprising four chromosomal species, 2n = 52, 54, 58 and 60. We reanalzed two previously available data sets of 1. body size differentiation (Nevo et al. 1986a) and 2. the intensity of “Total Aggression” in mole rats (Nevo et al. 1986b). We correlated the mean size difference between the two sexes, in each of the 12 populations of the chromosomal species, with the mean level of agression, and with climatic factors, both displaying significant correlations. The results indicated that for 2n = 52, 54, 58 and 60, the population averae difference in body weight between the sexes decreases southward as follows: 37.7g (30.8 % of females body weight), 39.3g(29.0%) 26.3g(22.8%) and 20.3g (19.3%), respectively. We interpret the higher body size diherential ketween the sexes in the north as due to sexual selection.     

Total body water and adaptive water turnover rate in four chromosomal species of subterranean mole rats of the Spalax ehrenbergi superspecies in Israel

Total body water (TBW) and water turnover rate (WTR) were measured in 24 subterranean mole rats comprising four populations, each belonging to a different chromosomal species of the Spalax ehrenbergi superspecies in Israel. The four species range in different climates: humid-cool (2n = 52); semiarid-cool (2n = 54); humid-warm (2n = 58) and arid-warm (2n = 60). TBW, as a percentage of body weight, measured by tritiated water (HTO), was 72.4%±4–7 in 2n = 52, significantly ( P < 005) higher than the similar estimates 61.7%± 7.2, and 59.4%± 5–3, for 2n = 60 and 58, respectively. A comparison of HTO space, as a percentage of TBW, closely approximated TBW, ranging from 97% to 108%. WTR was high, 218.1 and 230.9 ml/kg^0-75/day in the mesic populations of 2n = 58 and 52, respectively. By contrast, WTR estimates were significantly lower , ( P < 0.001), 150.2 and 148.9 ml/kg^0-75/day in the xeric populations of 2n = 54 and 60, respectively. The biological half-life time, T_1/2, was similar and faster, 32.7 and 27.9 hours in the mesic populations of 2n = 52 and 58, as compared with slower, 47.9 and 40.8 hours in the xeric populations of 2n = 54 and 60, respectively. Urine osmolarity in the most xeric northern Negev steppic population of 2n = 60 (737 ± 45 mmol/kg) was significantly ( P < 0.001) higher than in the other species.
We conclude that adaptive radiation in the Spalax ehrenbergi superspecies involves speciation in semiarid (2n = 54) and arid (2n = 60) climates by physiological adaptations of kidney water conservation, along with multiple morphological, physiological and behavioural syndromes of climatic adaptations to increasing aridity (Nevo, 1986).     

Karyotype and genetic evolution in speciation of subterranean mole rats of the genus Spalax in Turkey

Karyotype (2fl) and allozymc diversity at 37 gene loci were determined in 69 subterranean mole rats in Turkey belonging to the two superspecies: the ancestor Spalax leucodon ( n = 55; 20 populations) and the descendant S. ehrenbergi ( n = 14; lour populations). We identified remarkable variation of diploid chromosome numbers in the S. leucodon superspecies: 2 n = 38. 40. 50, 54, 60 and 62; and in the S. ehrenbergi superspecies: 2 n = 52, 56 and 58. Genetic diversity indices were low on average in both S. leucodon and S. ehrenbergi superspecies: Allele diversity, A = 1.081 and 1.074; polymorphism, P-5%= 0.077 and 0.068; heterozygosity, H = 0.038 and 0.027; and gene diversity, H, = 0.038 and 0.034, respectively. H ranged from 0 in mesic or semimesic regions to 0.088 in arid Anatolia. We consider the populations with different diploid chromosome numbers, 2 n , as good biological species. Karyotypic diversity may mark extensive ecological speciation. Nei's genetic distances, D (average 0.174, range 0.002 0.422) and ecogeographical criteria suggest that almost each population may represent a different biological species, but critical future testing is necessary to support this claim. Karyotypes and allozymes are nonrandomly distributed across Turkey, displaying remarkable correlations with climatic and biolir factors. Both In and H are significantly correlated with aridity stress (2 n /rainfall, r , =–0.74; P < 0.001), and in our region also with climatic unpredictability. These results support the niche-width genetic variation hypothesis in space and time. Climatic selection in Turkey appears to be a major architect of karyotype and genetic (allozymc) diversity and divergence in mole rat evolution, in both speciation and adaptation.     

DNA sequence variation in the mitochondrial control region of subterranean mole rats,Spalax ehrenbergi superspecies,in Israel

The complete mitochondrial control region was sequenced for 60 individuals representing different populations for each of the four species of the subterranean mole rat Spalax ehrenbergi superspecies in Israel: Spalax galili (2n = 52), S. golani (2n = 54), S. carmeli (2n = 58), and S. judaei (2n = 60). The control region of all species and populations is very similar both in length (979 to 983 bp) and in base composition. As in agreement with previous surveys on mitochondrial control regions on mammals, the mole rat control region can be divided into a central domain and two flanking domains, ETAS (extended termination associated sequences) and CSB (conserved sequence blocks). Along with the common conserved blocks found in these domains (ETAS1, ETAS2, CSB1, CSB2, and CSB3), we have also detected in all individuals an ETAS1-like and a CSB1-like element, both in the ETAS domain. The most conserved region was the central domain, followed by the CSB and ETAS domains, showing important differences in the four species analyzed. Phylogenetic analysis supported the existence of two clades. One clade contained individuals belonging to Spalax galili (2n = 52) and S. golani (2n = 54), separated in two different branches depending on the species. The other clade contained individuals belonging to S. carmeli (2n = 58) and S. judaei (2n = 60) mixed together, suggesting a more recent event of speciation. Within species we have observed a southward trend of increasing variability. These results have been explained as a consequence of the adaptation of the species to ecological factors such as aridity and temperature stresses.     

Speciation and adaptive radiation of subterranean mole rats, Spalax ehrenbergi superspecies, in Jordan

The major initial mechanism of speciation in subterranean blind mole rats, Spalacidae, is chromosomal, primarily through Robertsonian rearrangements. Here we highlight another scenario of chromosomal rearrangement leading to ecological speciation and adaptive radiation apparently initiated by pericentric inversions and genie divergence to different ecologies in mole rats in Jordan. We analysed karyotype, allozyme, size and ecological diversity across the range of mole rats in Jordan from mesic Irbid in the north to xeric Wadi Musa (Petra region) in the south, a transect of 250 km. We examined mole rats for chromosome ( N =71), size ( N =76), and allozyme ( N =67) diversities, encoded by 32 loci, in 12 populations of the Spalax ehrenbetgi superspecies in Jordan. By a combination of chromosome morphology, genetic distance, body size and ecogeography, we identified four new putative biological species. All species (except two animals in Madaba) share 2 n = 60 but vary in chromosome morphology, caused by pericentric inversions and/or centromeric shifts. The 'north Moav' species is karyotypically polymorphic for 2 n (2 n = 60; including locally also two animals with 2 n = 62). The distribution of the four species is associated with ecogeographical different domains and climatic diversity. Genetic diversity indices were low, but like chromosome arms (NFa) were positively correlated with aridity stress. Discriminant analysis correctly classified 91% of the individuals into the four species utilizing combinatorially chromosome, allozyme and size diversities. It is hypothesized that mole rat evolution underground is intimately associated with climatic diversity stress above ground.     

Microsatellite diversity in populations of blind subterranean mole rats (Spalax ehrenbergi superspecies) in Israel: speciation and adaptation

The genetic diversity, structure, and relationships of 12 populations of blind subterranean mole rats belonging to four species of the Spalax ehrenbergi superspecies were examined. Diversity levels and divergence of 19 nuclear microsatellite markers were analysed. Results indicate that the 12 populations fall into two major clusters representing the species pairs Spalax galili (2 n  = 52) and S. golani (2 n  = 54) vs. S. carmeli (2 n  = 58) and S. judaei (2 n  = 60). There is very little gene flow between the species pairs, except for one population (Afiq) of S. carmeli (2 n  = 58), which is located near the hybrid zone between S. golani (2 n  = 54) and S. carmeli (2 n  = 58). Higher values of the mean for alleles per locus, gene diversity, and number of unique alleles were observed for the southern species pair, S. carmeli and S. judaei, than for the northern species pair, S. galili and S. golani . Microsatellite diversity increases southward toward the Negev Desert and is positively correlated with aridity stress. Climatic factors of water availability coupled with temperature account for 76.7% of microsatellite variance. Natural selection appears to adaptively determine microsatellite evolution in Spalax regionally across Israel.  © 2004 The Linnean Society of London , Biological Journal of the Linnean Society , 2004, 83 , 229–241.     

Mhc Class II DNA Polymorphisms within and between Chromosomal Species of the Spalax Ehrenbergi Superspecies in Israel

Restriction fragment length polymorphisms (RFLPs) of two major histocompatibility class II genes (P alpha 1 and Q beta) were studied in 13 populations of four chromosomal species (2n = 52, 54, 58 and 60) of the mole rat, Spalax ehrenbergi superspecies in Israel. A substantial frequency of allelic fragments was found in both genes for all populations, including a desert isolate. In the P alpha 1 gene, one allelic fragment is a result of a deletion mutation which is diagnostic of the 2n = 52 chromosomal species. All other ten allelic variants are the result of point mutations. All mutations are located in a short region flanking the 3' end of the gene. Based on Mhc polymorphisms we confirm earlier evidence that gene flow does not occur between the older chromosomal species (2n = 52, 54, 58), and that reproductive isolation decreases, progressively from the oldest to the youngest species (2n = 60).     

MtDNA polymorphisms: evolutionary significance in adaptation and speciation of subterranean mole rats

Patterns of mtDNA diversity in subterranean mole rats of the Spalax ehrenbergi superspecies (2n = 52, 54, 58 and 60) were previously studied in the main ranges of the four chromosomal species, and specifically in the 2n = 60 species and its peripheral steppe semi-isolates and desert isolates. In the present study we correlated mtDNA diversity indices, nucleon diversity, h , and nucleotide divergence, π, with physical (climatic), biotic (parasites) and biological (genetical, morphological, physiological and behavioural) factors, showing that mtDNA diversity is structured ecogeographically and biologically. The following significant correlations of mtDNA diversity were indicated with: (i) climatic heterogeneity and unpredictability; (ii) levels of ecto- and endoparasites; and (iii) biological diversities, primarily with physiological diversity associated with the energy budget. Small steppe semi-isolates and desert isolates harbour high levels of mtDNA haplotype diversity, some novel, which may be a prerequisite for future speciation events. We conclude that the ecogeographical and biological correlates, as well as the maintenance of mtDNA polymorphisms in small isolated populations, strongly suggest that mtDNA diversity is not neutral. Diversifying natural selection appears to be an important differentiating factor of mtDNA diversity in the twin evolutionary processes of adaptive radiation and active speciation. We suggest critical experiments to substantiate our conclusions and highlight the contribution of mtDNA diversity to fitness, i.e. to the biological function of mtDNA diversity in the evolutionary process.     

Habitat selection in evolving mole rats

Summary The actively speciating four chromosomal species of fossorial mole rats of the Spalax ehrenbergi complex in Israel (2n=52, 58, 54 and 60) which inhabit an increasingly arid environment in this order were tested to determine their habitat preference. The testing apparatus simulated four climatic regimes based on temperature and humidity combinations corresponding to the climatic origins of the four chromosomal species: coolhumid, cool-dry, warm-humid and warm-dry, respectively. The tests involved 175 adults comprising all four chromosome species and representing 10 populations. Out of the 139 analyzed animals 88% selected the warm cages and only 12% selected the cool cages. The four karyotype progressively preferred the warm-dry cage in the following order: 53, 59, 60 and 72% for 2n=58, 52, 54 and 60 respectively, largely in accord with their increasingly arid climatic origins. Even larger differences were found in populations within karyotypes in accord with the local climatic variation within a karyotype range.The results of our analysis indicate that the chromosomal species and populations select their climatic habitat in accord with the climatic conditions of their geographic localities. The humidity index appears to be the prime differentiator of habitat selection and may have been a substantial ecological factor in species differentiation and distribution of the Spalax ehrenbergi complex in Israel.     

Mole rat hemoglobin: primary structure and evolutionary aspects in a second karyotype of Spalax ehrenbergi, Rodentia, (2n = 52)

The hemoglobins of the four karyotypes of Spalax ehrenbergi (2n = 52, 54, 58, 60) did not show any differences in their electrophoretic pattern and in high performance liquid chromatography. The complete amino-acid sequence of mole rat hemoglobin (Spalax ehrenbergi), chromosome species 2n = 52, is presented. It was elucidated by automatic Edman degradation of the chains, the tryptic peptides, and the C-terminal peptide obtained by acid hydrolysis of the Asp-Pro bond in beta-chains. The alpha- and beta-chains are identical with those of the chromosome species 2n = 60. A comparison of the hemoglobins of mole rat, mouse, and other rodents shows homology but no indication of adaptation to subterranean life. In all probability alpha 11(A9)Arg and alpha 120(H3)Gly, unique in mole rat among all mammalian hemoglobins, are not involved in high oxygen affinity. The construction of a phylogenetic tree by the maximum parsimony method, based on hemoglobin sequences, made it possible to show that Rodentia originated as a monophyletic clade, and to find the phylogenetic relationship of Spalacidae to other Rodentia (Mus, Rattus, Ondatra, Mesocricetus, Citellus, and Cavia). Among all rodents the slowest rate of nucleotide replacements occurred in the lineage to Spalax (20%) and the fastest in the lineage to Cavia (59%).     

Ribosomal DNA non-transcribed spacer polymorphism in subterranean mole rats: genetic differentiation,environmental correlates and phylogenetic relationships

Summary An analysis is presented of genetic differentiation in the non-transcribed spacers of ribosomal DNA (NTS rDNA). Diversity, environmental correlates and the phylogenetic relationships are examined within and between species of the actively speciating subterranean mole rat, superspeciesSpalax ehrenbergi (2n=52, 54, 58, 60) in Israel. This analysis is based on a previous study of the geographic distribution of restriction fragment length polymorphisms of NTS rDNA. Here we present results indicating that NTS rDNA diversity exists mostly (66%) within populations, while 20% is between populations within species, and 14% between species. Multivariate discriminant analysis succeeded in separating 10 of the 13 populations (77%) into their correct chromosomal species, on the basis of the combination of three NTS rDNA repetypes. The phylogenetic relationships suggest that the complex involves two pairs of closely related species (2n=52–54 and 2n=58–60). NTS rDNA diversity, as well as the decrease southward in frequency of repetype C, are correlated with climatic factors of humidity and temperature. These data are discussed in terms of the evolutionary forces of migration and selection which may cause NTS rDNA differentiation. Climatic selection appears to be the major differentiating factor of NTS rDNA.     

Hypoxic survival differs between two mole rat species (Spalax ehrenbergi) of humid and arid habitats.

1. Two chromosomal species, 2n = 52 and 2n = 60 of the mole rat superspecies (Spalax ehrenbergi), occupy humid (2n = 52) and arid (2n = 60) habitats in Israel. 2. Gas conductivity of the soil of the 2n = 52 mole rat is lower than that of the 2n = 60 mole rat, and the 2n = 52 mole rat is better adapted to hypoxia. 3. The hypothesis that the 2n = 52 mole rat can survive to a lower pO2 than the 2n = 60 mole rat was tested. 4. Terminal pO2 (Torr) of four females 2n = 52 was lower, 18.0 +/- 2.9 (SD), than the terminal pO2 of five females 2n = 60, 28.2 +/- 5.1 (SD). 5. The hypoxic survival of the 2n = 52 mole rat as compared to that of the 2n = 60 mole rat correlates with other physiological traits: breathing and heart frequencies, blood hemoglobin and tissue gas tensions.     

Mitochondrial DNA polymorphisms in subterranean mole-rats of the Spalax ehrenbergi superspecies in Israel, and its peripheral isolates

Patterns of mitochondrial DNA (mtDNA) variation were examined in 133 mole-rats constituting all four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of the Spalax ehrenbergi superspecies in Israel, as well as the peripheral isolates of 2n = 60. In the main range of the complex, a total of 28 mtDNA haplotypes were found in 64 mole-rats, with most haplotypes being unique to either a single chromosomal species or population. mtDNA divergence increased from low to high diploid number in a north-to-south direction in Israel. Overall levels of mtDNA diversity were unexpectedly the highest in the 2n = 60, the youngest species of the complex. The mtDNA haplotypes can be separated into two major groups, 2n = 52-54 and 2n = 58-60, and a phylogenetic analysis for each group revealed evidence of a few haplotypes not sorted by diploid number. The overall patterns of mtDNA divergence seen within and among the four chromosomal species are consistent with the parapatric mode of speciation as suggested from previous studies of allozyme and DNA hybridization. In a separate data set the patterns of mtDNA variation were examined across the main geographic range and across peripheral semi-isolates and isolates of the 2n = 60 chromosomal species. Fifteen haplotypes were found in 69 mole-rats. High levels of mtDNA diversity characterized the main range, semi-isolated, and even some desert isolated populations. The peripheral isolates contain much mtDNA diversity, including novel haplotypes.      

Adaptive variation of pelage color within and between species of the subterranean mole rat (Spalax ehrenbergi) in Israel

Summary Color notations of dorsal pelage were analyzed in 451 adult subterranean mole rats, comprising 19 populations and 4 chromosomal species (2n=52, 54, 58 and 60) of the Spalax ehrenbergi superspecies in Israel. In addition, the color notations of soil samples from the collection sites were measured. In an attempt to evaluate the degree of correlation between pelage color of mole rats and the local soil color, each color was compared both macro- and microgeographically. The macrogeographic comparisons were among populations of the 4 chromosomal species, and the microgeographic comparisons were among pairs of geographically neighboring subpopulations which occur under similar climatic conditions but differ in soil type and color. The results indicate that 1. no differences in pelage color notation exist between the sexes; 2. pelage of mole rats is generally characterized by an over-all grayish color, but its shades vary regionally and locally in accordance with varying soil color; 3. a positive correlation was found between pelage and local soil colors in the macrogeographic analysis and this correlation was verified in a microgeographic comparison; 4. specimens from populations of the northern chromosomal species (2n=52, 54) mainly inhabit terra rossa and basalt soils, respectively, with reddish brown tones, and they tend to have a more reddish shade, whereas animals from the central species (2n=58), which inhabits mainly dark soils including alluvial clays, terra rossa and brown rendzina and the souther species 2n=60, which occurs mainly in light soil types such as pale rendzina, sandy loams and loess, tend to have a more yellowish shade. The general appearance of specimens from 2n=58 was dark gray while that of 2n=60 specimens was light gray; 5. mole rats living in xeric environments (particularly 2n=60, but also 2n=54) are lighter than those which live in mesic environments (2n=52, 58). The results support the conclusion that pelage color of strictly subterranean mole rats is subject to selection pressures of overground predation of disharmonious types with their background soil color. The possibility that the pelage color variation and patterns also contributes to better thermoregulation while mole rats are above ground is discussed.     

Chromosomal polymorphisms in subterranean mole rats: origins and evolutionary significance

Two distinct chromosomal polymorphisms, among others, characterize the short (p) and long (qh) arms of chromosome number 1 of the Spalax ehrenbergi complex in Israel. We have studied the geographic distributions of these two polymorphisms in 60 animals belonging to four chromosomal species (2 n = 52, 54, 58 and 60). These comprised 15 populations: 12 from continuously distributed populations, two from semi-isolates, and one isolate in the northern Negev Desert.
Our results indicate that: (i) the two polymorphisms are widespread not only within populations and species but most strikingly between different bone marrow cells within an individual; this may reflect a diversity generating mechanism; (ii) the two polymorphisms of p and qh, primarily of p, are correlated with climatic factors of water availability and temperature, as well as geographical distances; these correlations are significantly above those expected by chance, (iii) The semi-isolates and the isolate display high levels of polymorphism in both p and qh, but particularly in qh, despite the small size population effects.
We suggest that at least the polymorphism in p is involved in the adaptive radiation of mole rats into diverse climatic regimes. Furthermore, we speculate that the forces generating chromosomal polymorphisms relevant to speciation (i.e. Robertsonian mutations), coupled with the evolutionary forces operating in small peripherally isolated populations, may be appropriate for the origin of new chromosomal species through peripatric speciation.     

Morphometrics of speciating mole rats: Adaptive differentiation in ecological speciation

The morphometrics of subterranean mole rats, Spalax ehrenbergi superspecies in Israel, were studied. A total of 42 skull and body variables of 327 adults from 44 populations across the ranges of the 4 chromosomal species (2n = 52, 54, 58 and 60) were measured. The results showed: a. in general, significantly higher morphological values in males than in females; b. significant interspecific quantitative diffferences in 15 skull variables of males, and in 3 of females; but no qualitative diagnosis of any of the species; c. an intercorrelation between most skull variables at different orgamsmal levels, and a uniform and complete correlation with weight; d. a statistical discrimination in both metric and nonmetric multivariate analyses between most species pairs; e. an indication of high morphological similarity based on the relatively small multivariate Mahalanobis distances of dissimilarity between species; f. a possible explanation for a significant part of the variance in most skull and body variables resulting from a combination of temperature variables and water availability; g. a declination in size in skull and body variables between northern and southern species; h. the order of speciation events to be 2n = (54, 52)→58→60. Based on these results we conclude: 1. no morphological breaks are associated with the speciation of the S. ehrenbergi complex, and morphological differentiation is quantitative rather than qualitative between the species: 2. morphological diversity is significantly explained by climatic selection, and is therefore adaptive and explicable on even very low selective pressures over evolutionary time: 3. the morphospecies seem unable to reliably reflect the biological species in the S. ehrenbergi complex, and morphological differentiation, at least in this superspecies, evolves in a gradual rather than a punctual way: this would conform with Neo-Darwinian evolutionary theory.     

Variations in ribosomal DNA and mitochondrial DNA among chromosomal species of subterranean mole rats

Restriction site variations in nuclear ribosomal DNA (rDNA) spacers and mitochondrial DNA (mtDNA) were examined in several populations of mole rats with variable numbers of chromosomes, which represented the two superspecies Spalax leucodon (2n = 38, 54, or 62) and Spalax ehrenbergi (2n = 52, 54, 58, or 60). Sequence divergence of rDNA spacers between the members of the superspecies was approximately 8%, while the variation within each superspecies was 4% on average. The intrasuperspecies differentiation of rDNA spacers was generally associated with changes in the diploid number of chromosomes. However, substantial divergence (approximately 1.5%) was also detected among populations with the same diploid number of chromosomes. The sequence divergence of mtDNA among haplotypes of S. ehrenbergi was 10% or higher and among haplotypes of S. leucodon it was approximately 12%. By contrast, the range of sequence divergence between superspecies was 7.4%-12%. The large divergence of mtDNAs within each superspecies of mole rats may be explained by the preservation of ancient mtDNA polymorphisms.      

Ribosomal DNA (rDNA) spacer polymorphism in mole rats

Genetic differentiation of ribosomal DNA (rDNA) nontranscribed-spacer (NTS) polymorphism was analyzed in 50 individuals from 13 populations among the four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of subterranean mole rats of the Spalax ehrenbergi complex in Israel. Southern blot analysis with a mouse rDNA probe and two restriction enzymes, EcoRI and BamHI, revealed various sizes of major restriction fragments. The assumption that this variation is due to length polymorphism of NTS DNA was supported by the construction of restriction-site maps. On the basis of the EcoRI and BamHi fragment lengths, we could characterize the major types of NTS rDNA repeating units in each individual. Each member of the central population in the four chromosomal groups of mole rats has a characteristic combination of the NTS types, suggesting that the karyotype groups were genetically diverged. Some near-hybrid-zone populations reflect similarities with the rDNA spectra of a neighbor chromosomal group. This might have resulted from gene flow across the hybrid zones.