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.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

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.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

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.     

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

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

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

Limited mitochondrial DNA introgression in a Bombina hybrid zone

We analysed the variation at mitochondrial DNA (mtDNA) and six allozyme loci in two transects across a hybrid zone between the fire-bellied toads Bombina bombina and Bombina variegata in southern Poland. The mtDNA cline was narrower than allozyme clines in one transect (Przemyśl, β = 0.435, w  = 3.42 km) and shifted to the B .  bombina side in both (near Kraków by 1.32 km, α = 0.437, and near Przemyśl by 0.97 km, α = 0.319). Cytonuclear associations were weak. Narrowed mtDNA clines could be a by-product of female demography and lowered effective population size of mtDNA. Alternatively, restricted mtDNA introgression is a likely result of negative epistasis in recombinants augmented by environment-dependent selection because divergence of Bombina mitochondrial and nuclear genomes is large. The shift of the mtDNA cline, contrary to expectations from initial mating preferences and fecundity differences between the species, suggests that on the B .  bombina side of the zone hybrid females with B .  variegata mtDNA have a higher chance of leaving progeny.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 295–306.     

When vicars meet: a narrow contact zone between morphologically cryptic phylogeographic lineages of the rainforest skink, Carlia rubrigularis

Phylogeographic analyses of the fauna of the Australian wet tropics rainforest have provided strong evidence for long-term isolation of populations among allopatric refugia, yet typically there is no corresponding divergence in morphology. This system provides an opportunity to examine the consequences of geographic isolation, independent of morphological divergence, and thus to assess the broader significance of historical subdivisions revealed through mitochondrial DNA phylogeography. We have located and characterized a zone of secondary contact between two long isolated (mtDNA divergence > 15%) lineages of the skink Carlia rubrigularis using one mitochondrial and eight nuclear (two intron, six microsatellite) markers. This revealed a remarkably narrow (width < 3 km) hybrid zone with substantial linkage disequilibrium and strong deficits of heterozygotes at two of three nuclear loci with diagnostic alleles. Cline centers were coincident across loci. Using a novel form of likelihood analysis, we were unable to distinguish between sigmoidal and stepped cline shapes except at one nuclear locus for which the latter was inferred. Given estimated dispersal rates of 90-133 m x gen(-1/2) and assuming equilibrium, the observed cline widths suggest effective selection against heterozygotes of at least 22-49% and possibly as high as 70%. These observations reveal substantial postmating isolation, although the absence of consistent deviations from Hardy-Weinberg equilibrium at diagnostic loci suggests that there is little accompanying premating isolation. The tight geographic correspondence between transitions in mtDNA and those for nuclear genes and corresponding evidence for selection against hybrids indicates that these morphologically cryptic phylogroups could be considered as incipient species. Nonetheless, we caution against the use of mtDNA phylogeography as a sole criterion for defining species boundaries.     

Hybridization between multiple fence lizard lineages in an ecotone: locally discordant variation in mitochondrial DNA, chromosomes, and morphology

We investigated a hybrid zone between two major lineages of fence lizards (Sceloporus cowlesi and Sceloporus tristichus) in the Sceloporus undulatus species complex in eastern Arizona. This zone occurs in an ecotone between Great Basin Grassland and Conifer Woodland habitats. We analysed spatial variation in mtDNA (N=401; 969 bp), chromosomes (N=217), and morphology (N=312; 11 characters) to characterize the hybrid zone and assess species limits. A fine-scale population level phylogenetic analysis refined the boundaries between these species and indicated that four nonsister mtDNA clades (three belonging to S. tristichus and one to S. cowlesi) are sympatric at the centre of the zone. Estimates of cytonuclear disequilibria in the population closest to the centre of the hybrid zone suggest that the S. tristichus clades are randomly mating, but that the S. cowlesi haplotype has a significant nonrandom association with nuclear alleles. Maximum-likelihood cline-fitting analyses suggest that the karyotype, morphology, and dorsal colour pattern clines are all coincident, but the mtDNA cline is skewed significantly to the south. A temporal comparison of cline centres utilizing karyotype data collected in the early 1970s and in 2002 suggests that the cline may have shifted by approximately 1.5 km to the north over a 30-year period. The recent northward expansion of juniper trees into the Little Colorado River Basin resulting from intense cattle overgrazing provides a plausible mechanism for a shifting hybrid zone and the introgression of the mtDNA haplotypes, which appear to be selectively neutral. It is clear that complex interactions are operating simultaneously in this contact zone, including the formation of hybrids between populations within S. tristichus having diagnostic mtDNA, morphology, karyotypes, and dorsal colour patterns, and secondary contact between these and a distantly related yet morphologically cryptic mtDNA lineage (S. cowlesi).     

Divergent character clines across a recent secondary contact zone in a Hispaniolan lizard

Studies of genetic contact zones provide valuable information regarding the processes of population divergence, adaptation and speciation. In this paper, I examine transitions in morphology, mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) haplotypes across a recent secondary contact zone in a Hispaniolan lizard Ameiva chrysolaema . Maximum likelihood cline fitting analyses suggest non-coincidence of cline centers and that the mtDNA cline is significantly displaced to the west of the remaining clines. nDNA and morphological clines are coincident and tend to be associated with the prevailing environmental gradient. The lack of cytonuclear disequilibrium near the center of the contact zone and the non-coincidence of character clines suggest that this zone does not conform to a tension zone model of hybridization; thus, gene flow across the zone does not seem to be impeded. The extremely narrow width of the dorsal scale size cline and the close association of this cline with the steepness of the environmental (precipitation) gradient suggest that this character may be under environmental selection. Taken together, this contact zone appears to be structured by a combination of mtDNA introgression, possibly associated with eastward movement of the zone, and environmental selection on some characters.     

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

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

Genetic architecture and genomic patterns of gene flow between hybridizing species of Picea

Hybrid zones provide an opportunity to study the effects of selection and gene flow in natural settings. We employed nuclear microsatellites (single sequence repeat (SSR)) and candidate gene single-nucleotide polymorphism markers (SNPs) to characterize the genetic architecture and patterns of interspecific gene flow in the Picea glauca × P. engelmannii hybrid zone across a broad latitudinal (40–60 degrees) and elevational (350–3500 m) range in western North America. Our results revealed a wide and complex hybrid zone with broad ancestry levels and low interspecific heterozygosity, shaped by asymmetric advanced-generation introgression, and low reproductive barriers between parental species. The clinal variation based on geographic variables, lack of concordance in clines among loci and the width of the hybrid zone points towards the maintenance of species integrity through environmental selection. Congruency between geographic and genomic clines suggests that loci with narrow clines are under strong selection, favoring either one parental species (directional selection) or their hybrids (overdominance) as a result of strong associations with climatic variables such as precipitation as snow and mean annual temperature. Cline movement due to past demographic events (evidenced by allelic richness and heterozygosity shifts from the average cline center) may explain the asymmetry in introgression and predominance of P. engelmannii found in this study. These results provide insights into the genetic architecture and fine-scale patterns of admixture, and identify loci that may be involved in reproductive barriers between the species.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

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.     

Cline coupling and uncoupling in a stickleback hybrid zone

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

A comparison of five hybrid zones of the weta Hemideina thoracica (Orthoptera: Anostostomatidae): degree of cytogenetic differentiation fails to predict zone width

Abstract Tension zones are maintained by the interaction between selection against hybrids and dispersal of individuals. Investigating multiple hybrid zones within a single species provides the opportunity to examine differences in zone structure on a background of differences in extrinsic factors (e.g., age of the zone, ecology) or intrinsic factors (e.g., chromosomes). The New Zealand tree weta Hemideina thoracica comprises at least eight distinct chromosomal races with diploid numbers ranging from 2n = 11 (XO) to 2n = 23 (XO). Five independent hybrid zones were located that involve races differing from one another by a variety of chromosomal rearrangements. The predicted negative correlation between extent of karyotypic differentiation (measured in terms of both percent of genome and number of rearrangements) and zone width was not found. Conversely, the widest zones were those characterized by two chromosome rearrangements involving up to 35% of the genome. The narrowest zone occurred where the two races differ by a single chromosome rearrangement involving approximately 2% of the genome. The five estimates of chromosomal cline width ranged from 0.5 km to 47 km. A comparative investigation of cline width for both chromosomal and mitochondrial markers revealed a complex pattern of zone characteristics. Three of the five zones in this study showed cline concordance for the nuclear and cytoplasmic markers, and at two of the zones the clines were also coincident. Zones with the widest chromosomal clines had the widest mitochondrial DNA clines. It appears that, even within a single species, the extent of karyotypic differentiation between pairs of races is not a good predictor of the level of disadvantage suffered by hybrids.     

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7-10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3-7 m) nor the wider nuclear cline (RAG-2; 34-42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2-3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.     

Evidence for strong assortative mating,limited gene flow,and strong differentiation across the blue‐footed/Peruvian booby hybrid zone in northern Peru

Hybrid zones represent natural laboratories in which the processes of divergence and genetic isolation can be examined. The generation and maintenance of a hybrid zone requires mispairing and successful reproduction between organisms that differ in one or more heritable traits. Understanding the dynamics of hybridization between two species requires an understanding of the extent to which they have diverged genetically, the frequency of misparing and hybrid production, and the extent of introgression. Three hundred and twenty one blue‐footed Sula nebouxii and Peruvian S. variegata boobies from the eastern tropical Pacific Ocean were analyzed using 19 putatively neutral genetic markers to evaluate interspecific differentiation, to classify morphological hybrids using Bayesian assignments, and to characterize hybridization using cline theory and Bayesian assignments. The species were well differentiated at mitochondrial and nuclear microsatellites, the hybrid zone was bimodal (contained a high frequency of each parental species but a low frequency of hybrids), and morphologically intermediate individuals were most likely F1 hybrids resulting from mating between female Peruvian boobies and male blue‐footed boobies. Clines in allele frequency could be constrained to share a common geographic centre but could not be constrained to share a common width. Peruvian and blue‐footed boobies hybridize infrequently, potentially due to strong premating reproductive isolation; however, backcrossing appears to facilitate introgression from blue‐footed to Peruvian boobies in this hybrid system.