A CONTACT ZONE WITH NONCOINCIDENT CLINES FOR SEX-SPECIFIC MARKERS IN THE FIELD VOLE (MICROTUS AGRESTIS)

A field vole (Microtus agrestis) population characterized by unique mitochondrial DNA (mtDNA) and Y chromosome markers occurs in southwest Sweden. A contact zone between this Lund (Lu) population and a standard (St) population was examined with two sex-specific genetic markers. The field vole mtDNA and Y chromosome clines display a remarkable lack of coincidence, rarely observed in contact zones. The cline width for both markers is about 50 km, but the two clines are displaced from each other: the mtDNA cline is found in the central part of the study area, whereas the cline for the Y chromosome is located in the eastern part. Thus, the absolute width of the Lu-St zone spans about 95–110 km. As a result of the cline shift, all male hybrids carry the Lu-Y chromosome and St-type mtDNA. The other possible male hybrid class is lacking. The distinct noncoincidence of the mtDNA and Y chromosome clines is most likely explained by selection against male hybrids with the St-Y chromosome and Lu-mtDNA. It is possible that incompatibilities between the maternal genome of Lu-type animals and the paternal genome of St-type individuals exist. However, alternative explanations based on neutrality or selective advantage cannot be totally dismissed.     

Natural hybridization between extremely divergent chromosomal races of the common shrew (Sorex araneus, Soricidae, Soricomorpha): hybrid zone in Siberia

Chromosomal races of the common shrew differ in sets of metacentric chromosomes and on contact may produce hybrids with extraordinarily complex configurations at meiosis I that are associated with reduced fertility. There is an expectation that these may be some of the most extreme tension zones available for study and therefore are of interest as potential sites for reproductive isolation. Here, we analyse one of these zones, between the Novosibirsk race (characterized by metacentrics go, hn, ik, jl, mp and qr) and the Tomsk race (metacentrics gk, hi, jl and mn and acrocentrics o, p, q and r), which form hybrids with a chain-of-nine (CIX) and a chain-of-three (CIII) configuration at meiosis I. At the Novosibirsk-Tomsk hybrid zone, the CIX chromosomes form clines of 8.53 km standardized width on average, whereas the cline for the CIII chromosomes was 52.83 km wide. The difference in these cline widths fits with the difference in meiotic errors expected with the CIX and CIII configuration, and we produce estimates of selection against hybrids with these types of configurations, which we relate to dispersal and age of the hybrid zone. The hybrid zone is located at the isocline at 200 m altitude above sea level; this relationship between the races and altitude is suggested at both coarse and fine scales. This indicates adaptive differences between the races that may in turn have been promoted by the chromosome differences. Thus, the extreme chromosomal divergence between the Novosibirsk and Tomsk may be associated with genic differentiation, but it is still striking that, despite the large chromosomal differences, reproductive isolation between the Novosibirsk and Tomsk races has not occurred.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

A comparison of nuclear and mitochondrial cline shapes in a hybrid zone in the Sceloporus grammicus complex (Squamata; Phrynosomatidae)

The F5 and FM2 chromosome races of the Sceloporus grammicus complex form a hybrid zone in the Mexican state of Hidalgo. Previous studies of this zone have assessed genetic structure by averaging estimates of shape and width across three diagnostic chromosome markers. This approach is likely to mask subtle differences in cline shape among loci (e.g. selected vs. neutral), and obscure any displacement of cline centres (if present). Here we use maximum likelihood methods to construct the best fitting individual clines for three chromosomal markers, and also add two new markers; the mitochondrial DNA (mtDNA) locus, and the nuclear ribosomal DNA (rDNA) repeat. For each locus, hybrid zone models were fitted by cline shape and width, and the position and number of segments describing the centre of the zone. Pairwise comparisons between all clines revealed concordance between chromosomes 2 and 6, but significant discordance in cline structure among all other paired combinations. The concordance of chromosomes 2 and 6 suggests that these clines are maintained by genome-wide forces. The discordance of the chromosome 1 cline suggests an influence of asymmetric introgression, while the mtDNA cline is probably influenced by selection and drift. The rDNA locus reveals a pattern best explained by either extreme asymmetric introgression or gene conversion. The structure of zone indicates that genome-wide processes and locus specific selective forces as well as drift, are operating to different degrees on different loci. The locus-by-locus approach used here permits a finer discrimination among possible mechanisms responsible for the maintenance of the individual clines.     

A hybrid zone between the Kirillov and Petchora chromosomal races of the common shrew (Sorex araneus L., 1758) in northeastern European Russia: a preliminary description

Karyotypes of the Petchora and Kirillov chromosomal races of the common shrew differ by six Robertsonian metacentrics with monobrachial homology, such that interracial F₁ hybrids produce a ring-of-six configuration at meiosis I and are expected to suffer infertility. Mapping of 52 karyotyped individuals by using a unique global positioning system (GPS) revealed that the Kirillov-Petchora hybrid zone is positioned close to the river Mezen, which separated these races, and so may limit the migration of shrews across the contact zone. Although the population density of shrews was found to be markedly different with respect to habitats, the zone runs through a mosaic of habitats that are similar for both the Petchora and the Kirillov sides. This is one of the narrowest chromosomal hybrid zones among those studied in Sorex araneus with a standard cline width of about 1 km. The center of the cline is located on a bank occupied by the Petchora race at a distance of 0.4 km away from a riverine barrier. Interestingly, both the Kirillov race and hybrid individuals were found on a small island in the middle of a river fully flooded each spring. The frequencies of karyotypic variants allow us to consider the zone as an example of a bimodal zone. New Robertsonian and de novo whole-arm reciprocal translocations (WART) chromosomal variants found in the zone could be regarded as evidences of current evolutionary process in chromosomal hybrid zones.     

Phenotypic Variation across Chromosomal Hybrid Zones of the Common Shrew (Sorex araneus) Indicates Reduced Gene Flow

Sorex araneus, the Common shrew, is a species with more than 70 karyotypic races, many of which form parapatric hybrid zones, making it a model for studying chromosomal speciation. Hybrids between races have reduced fitness, but microsatellite markers have demonstrated considerable gene flow between them, calling into question whether the chromosomal barriers actually do contribute to genetic divergence. We studied phenotypic clines across two hybrid zones with especially complex heterozygotes. Hybrids between the Novosibirsk and Tomsk races produce chains of nine and three chromosomes at meiosis, and hybrids between the Moscow and Seliger races produce chains of eleven. Our goal was to determine whether phenotypes show evidence of reduced gene flow at hybrid zones. We used maximum likelihood to fit tanh cline models to geometric shape data and found that phenotypic clines in skulls and mandibles across these zones had similar centers and widths as chromosomal clines. The amount of phenotypic differentiation across the zones is greater than expected if it were dissipating due to unrestricted gene flow given the amount of time since contact, but it is less than expected to have accumulated from drift during allopatric separation in glacial refugia. Only if heritability is very low, N_e very high, and the time spent in allopatry very short, will the differences we observe be large enough to match the expectation of drift. Our results therefore suggest that phenotypic differentiation has been lost through gene flow since post-glacial secondary contact, but not as quickly as would be expected if there was free gene flow across the hybrid zones. The chromosomal tension zones are confirmed to be partial barriers that prevent differentiated races from becoming phenotypically homogenous.     

Narrow hybrid zone between Moscow and Western Dvina chromosomal races and specific features of population isolation in common shrew Sorex araneus (mammalia)

The contact zone between Moscow and Western Dvina chromosomal races of common shrew Sorex araneus L. at the south of the Valdai Hights was traced over a distance of 20 km. Within this, close to parapatric, contact zone of chromosomal races the width of sympatry zone was about 500 m (the narrowest among currently known hybrid zones), and the proportion of hybrids was 24.3%. It was shown that in bimodal hybrid zones between chromosomal races of common shrew the width of sympatry zones varied from 0.5 to 13 km. This width does not correlate with the cytogenetic features of the hybrids, and seems to be determined by competitive relations between the races. The hybrid proportion is determined by the type of hybrid heterozygosity, and decreased in the race sympatry zone from 33-40 to 21.5-25.2%. The decrease of the hybrid proportion can be associated with the abnormal fertility of either the first generation, or the backcross hybrids.     

Narrow hybrid zone between Moscow and Western Dvina chromosomal races and specific features of population isolation in common shrew <Emphasis Type="Italic">Sorex araneus</Emphasis> (Mammalia)

The contact zone between Moscow and Western Dvina chromosomal races of common shrew Sorex araneus L. at the south of the Valdai Hights was traced over a distance of 20 km. Within this, close to parapatric, contact zone of chromosomal races the width of sympatry zone was about 500 m (the narrowest among currently known hybrid zones), and the proportion of hybrids was 24.3%. It was shown that in bimodal hybrid zones between chromosomal races of common shrew the width of sympatry zones varied from 0.5 to 13 km. This width does not correlate with the cytogenetic features of the hybrids, and seems to be determined by competitive relations between the races. The hybrid proportion is determined by the type of hybrid heterozygosity, and decreased in the race sympatry zone from 33–40 to 21.5–25.2%. The decrease of the hybrid proportion can be associated with the abnormal fertility of either the first generation, or the backcross hybrids.     

Inferences from a rapidly moving hybrid zone

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e.. all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the 'wake" of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline movement is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter pure fatima genotypes. This movement can be explained if random dispersal of 7-19 km x gen(1/2) is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.     

A MODEL OF A HYBRID ZONE BETWEEN TWO CHROMOSOMAL RACES OF THE COMMON SHREW (SOREX ARANEUS)

The common shrew (Sorex araneus) is subdivided into several karyotypic races in Britain. Two of these races meet near Oxford o form the “Oxford-Hermitage” hybrid zone. We present a model which describes this system a; a “tension zone,” i.e., a set of clines maintained by a balance between dispersal and selection against chromosomal heterozygotes. The Oxford and Hermitage races differ by Robertsonian fusions with monobrachial homology (kq, no versus ko), and so Fl hybrids between them would have low fertility. However, the acrocentric karyotype is found at high frequency within the hybrid zone, so that complex Robertsonian heterozygotes (kg no/q ko n) are replaced by more fertile combinations, such as (kg no/k q n o). This suggests that the hybrid zone has been modified so as to increase hybrid fitness. Mathematical analysis and simulation show that, if selection against complex heterozygotes is sufficiently strong relative to selection against simple heterozygotes, acrocentrics increase, and displace the clines for kg and no from the cline for ko. Superimposed on this separation is a tendency for the hybrid zone to move in favor of the Oxford (kg no) race. We compare the model with estimates of linkage disequilibrium and cline shape made from field data.     

A hybrid zone comprising staggered chromosomal clines in the house mouse (Mus musculus domesticus).

In the vicinity of John o'Groats (Caithness, Scotland) there is a small karyotypic race of the house mouse (Mus musculus domesticus) characterized by a diploid number of 32 chromosomes, including the metacentrics 4.10, 9.12, 6.13 and 11.14. This race forms a hybrid zone with the standard British race (fully acrocentric chromosomes, 2n = 40). Although hybrid zones normally consist of several (or many) narrow character clines at the same position, this zone is unusual in that the chromosomal clines do not coincide. The cline for arm combination 11.14 is staggered relative to the 6.13 cline and both are separate from the clines for 4.10 and 9.12 (which may or may not coincide). A variety of explanations for the structure of the hybrid zone are discussed. It is possible that this may be a case of 'zonal raciation'.     

THE GENETIC STRUCTURE OF A HYBRID ZONE BETWEEN TWO CHROMOSOME RACES OF THE SCELOPORUS GRAMMICUS COMPLEX (SAURIA,PHRYNOSOMATIDAE) IN CENTRAL MEXICO

The F5 (2n = 34) and FM2 (2n = 44–46) chromosome races of the Sceloporus grammicus complex form a parapatric hybrid zone in the Mexican state of Hidalgo, characterized by steep concordant clines among three diagnostic chromosome markers across a straight-line distance of about 2 km. Here, we show that this zone is actually structured into local patches in which hybridization extends over an extremely irregular front. The distribution of hybrid-index (HI) scores across the transect reveals some hybridization at almost all localities mapped in a central 7 km × 3 km area. Pooling the central samples produces both a strong heterozygote deficit for all diagnostic markers and strong linkage disequilibria between all pairwise combinations of these (unlinked) markers. Moreover, a highly significant association exists between the habitat on which each individual was caught and its karyotype (F5 chromosomes are more likely to be found on oak). Analysis of genotype frequencies over a range of spatial scales shows that there is no significant heterozygote deficit or habitat association within local areas of less than about 200 m; however, there is significant linkage disequilibrium over the smallest scales (R = D (pquv)^1/2 = 0.29, support limits, 0.18–0.36) over 100 m. These patterns suggest that lizards mate and choose habitats randomly within local patches. This conclusion is supported by mark-recapture estimates of dispersal (≈ 80 m in a generation) and by inference of matings from embryo and maternal karyotypes. Closer examination of the two-dimensional pattern reveals a convoluted cline for all three markers, with a width of 830 m (support limits 770 m–930 m). This cline width, combined with the strength of local linkage disequilibrium, implies a dispersal rate of σ = 160 m in a generation and an effective selection pressure of 30% on each chromosome marker. The proportion of inviable embryos is greater in females from the center of the hybrid zone; this is caused by effects associated with both karyotype and location. The hybrid zone is likely to be maintained by selection against chromosomal heterozygotes, by other kinds of selection against hybrids, and by selection adapting the chromosome races to different habitats. The structure of the contact may be caused by both random drift and by selection in relation to habitat.     

Staggered clines in a hybrid zone between two chromosome races of the harvestman Gagrellopsis nodulifera (Arachnida: Opiliones)

We analyzed a hybrid zone between two chromosome races (2n = 16 and 2n = 22) of a Japanese harvestman, Gagrellopsis nodulifera Sato and Suzuki (Arachnida: Opiliones: Phalangiidae). The hybrid zone is located in the eastern part of Tottori Prefecture, western Honshu. The width of the zone is approximately 5 to 15 km. Three independent tandem fusions/fissions seem to be the main cause of the karyotypic differences between the parental races. Ten karyotypic variants were found in the hybrid zone. They differed by numbers of diploid chromosomes and trivalents detected in meiosis. In most of the collecting sites, karyotypic heterozygotes were less common than expected. A positive correlation was found between number of trivalents in a karyotype and its deficiency rate. In some sites, the deficit of heterozygous individuals was accompanied by an excess of the intermediate homozygotes. One of the three transects across the zone was studied in detail. We found that three types of single heterozygotes (2n = 17, 2n = 19 and 2n = 21) formed a series of successive, spatially separated peaks along the transect. Two types of intermediate homozygotes (2n = 18 and 2n = 20) were also spatially separated. The most parsimonious explanation of such a structure is the staggering of clines of three tandem (or Robertsonian) fusion/fission variants that differentiate the parental races caused by selection against multiple heterozygotes. Analysis of nondisjunction in single heterozygotes demonstrated that there was a strong interindividual variation in nondisjunction rate. The mean frequency of aneuploid MII in single heterozygotes was 0.10 +/- 0.03. Crossover exchanges in some critical regions of trivalents result in abnormal chromosomal configurations: chromosomes with unequal chromatids and dicentric chromosomes. Frequency of crossover-induced chromosomal abnormalities was low in single heterozygotes (approximately equal to 4%), and was unexpectedly high in the double heterozygotes (approximately equal to 15%). Selection against karyotypic heterozygotes is considered as a main evolutionary force responsible for the structuring of the hybrid zone. A positive association between diploid chromosome number and altitude was found. The race 2n = 16 tended to occupy lower altitudes than the 2n = 22 parental race. Differences in ecological preferences may be a result of previous adaptations to different environments in allopatry. A hypothesis concerning the origin and evolution of the hybrid zone is proposed.     

Natural hybridization between extremely divergent chromosomal races of the common shrew (Sorex araneus, Soricidae, Soricomorpha): hybrid zone in European Russia

The Moscow and Seliger chromosomal races of the common shrew differ by Robertsonian fusions and possibly whole‐arm reciprocal translocations (WARTs) such that their F₁ hybrids produce a chain‐of‐eleven configuration at meiosis I and are expected to suffer substantial infertility. Of numerous hybrid zones that have been described in the common shrew, those between the Moscow and Seliger races involve the greatest chromosomal difference. We collected 211 individuals from this zone to generate a total dataset of 298 individuals from 187 unique global positioning system (GPS) locations within the vicinity of interracial contact. We used a geographic information system (GIS) to map the location of the hybrid zone, which follows a direct route between two lakes, as would be anticipated from tension zone theory. Even within the central area of the hybrid zone, there is a much higher frequency of pure race individuals than hybrid, making this a clear example of a bimodal zone in the sense of Jiggins & Mallet (2000) . The zone runs through good habitat for common shrews, but nevertheless it is very narrow (standard cline widths: 3–4 km), as would be anticipated from low hybrid fitness. There is clear potential for an interruption to gene flow and build‐up of reproductive isolation. As found in some other hybrid zones, there is a high frequency of novel genetic variants, in this case, new chromosomal rearrangements. Here, we report a de novo Robertsonian fission and a de novo reciprocal translocation, both for the first time in the common shrew. There is an extraordinarily high frequency of de novo mutations recorded in F₁ hybrids in the zone and we discuss how chromosomal instability may be associated with such hybrids. The occurrence of a de novo Robertsonian fission is of considerable significance because it provides missing evidence that fissions are the basis of the novel acrocentric forms found and apparently selected for in certain common shrew hybrid zones.     

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.     

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.     

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.     

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.     

Molecular patterns of introgression in a classic hybrid zone between the Australian tree frogs,Litoria ewingii and L. paraewingi: evidence of a tension zone

Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization.