Differential gene flow of mitochondrial and nuclear DNA markers among chromosomal races of Australian morabine grasshoppers (Vandiemenella, viatica species group)

Recent theoretical developments have led to a renewed interest in the potential role of chromosomal rearrangements in speciation. Australian morabine grasshoppers (genus Vandiemenella, viatica species group) provide an excellent study system to test this potential role of chromosomal rearrangements because they show extensive chromosomal variation and formed the basis of a classic chromosomal speciation model. There are three chromosomal races, viatica19, viatica17, and P24(XY), on Kangaroo Island, South Australia, forming five parapatric populations with four putative contact zones among them. We investigate the extent to which chromosomal variation among these populations may be associated with barriers to gene flow. Population genetic and phylogeographical analyses using 15 variable allozyme loci and the elongation factor-1alpha (EF-1alpha) gene indicate that the three races represent genetically distinct taxa. In contrast, analyses of the mitochondrial cytochrome c oxidase subunit I (COI) gene show the presence of three distinctive and geographically localized groups that do not correspond with the distribution of the chromosomal races. These discordant population genetic patterns are likely to result from introgressive hybridization between the chromosomal races and range expansions/contractions. Overall, these results suggest that reduction of nuclear gene flow may be associated with chromosomal variation, or underlying genetic variation linked with chromosomal variation, whereas mitochondrial gene flow appears to be independent of this variation in these morabine grasshoppers. The identification of an intact contact zone between P24(XY) and viatica17 offers considerable potential for further investigation of molecular mechanisms that maintain distinct nuclear genomes among the chromosomal races.     

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.     

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.     

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.     

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, morphological, and ecological characterization of a hybrid zone that spans a migratory divide

This study characterizes a hybrid zone that spans a migratory divide between subspecies of the Swainson's thrush (Catharus ustulatus), a long distance migratory songbird, in the Coast Mountains of British Columbia. To assess the potential for a barrier to gene flow between the subspecies, I: (1) analyzed the shape and width of genetic and morphological clines relative to estimates of dispersal distance, (2) assessed the ratio of parental to hybrid genotypes across the hybrid zone, (3) estimated population density across the hybrid zone, and (4) compared the spatial relationship between the hybrid zone and an existing environmental gradient. The results indicate that the hybrid zone is characterized by mostly concordant character clines that are narrow relative to dispersal, the absence of a hybrid swarm, and low population density at the center of the zone. This hybrid zone and additional regions of contact between these subspecies are found on the border between coastal and interior climatic regions throughout the Pacific Northwest. An identified shift in the location, but not the width, of the mtDNA cline relative to the nuclear clines is consistent with asymmetrical hybridization. Neutral diffusion of populations following secondary contact and hybrid superiority within an ecotone are insufficient explanations for the observed patterns. The hypothesis that best fits the data is that the Swainson's thrush hybrid zone is a tension zone maintained by dispersal and ecologically mediated barriers to gene flow.     

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.     

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.     

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.     

MITOCHONDRIAL DNA AND PROTEIN DIFFERENTIATION BETWEEN HYBRIDIZING CYTOTYPES OF THE WHITE-FOOTED MOUSE,PEROMYSCUS LEUCOPUS

Restriction-enzyme analysis of mitochondrial DNA and protein electrophoresis were used to document patterns of gene flow across a hybrid zone between chromosomal races of Peromyscus leucopus. Chromosomal markers (three inversions) are such that individuals can be classified as potential F₁'s, backcrosses, or parental types. Allozymic characterization of the hybrid zone is congruent with the chromosomal data (Stangl, 1986) and indicates an assymetrical distribution of markers, with the northeastern markers being distributed at a higher frequency into southwestern populations. Restriction patterns of mtDNA suggest that the two cytotypes may have had different evolutionary histories, and the distribution of haplotypes is concordant with other genetic markers used to identify the hybrid zone. Concordant changes in chromosomes, allozymes, and mtDNA suggest that the most viable hypothesis for the origin of the zone is secondary contact. A unique aspect of this study is that the same individuals were used for protein electrophoresis, mtDNA analysis, and chromosomal analysis. Thus, patterns of genetic variation can be interpreted free of any historical bias associated with samples collected at different times.     

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.     

Direct selection on allozymes is not required to explain heterogeneity among marker loci across a Mytilus hybrid zone

Unequal differentiation between two types of loci (allozyme and DNA markers) across a Mytilus hybrid zone has recently been claimed as evidence for direct selection on some allozyme loci. We provide here a counter-example: a noncoding DNA locus that exhibits as much differentiation as the incriminated allozymes do. The levels of genetic differentiation varied widely among both allozymes and noncoding DNA markers and no clear difference emerged between the two types of markers. This suggests that the strong interlocus variance in genetic differentiation has been confounded with a discrepancy between marker types as a result of an insufficient and unbalanced locus sampling. Heterogeneity in differentiation among neutral loci can be created by stochastic variance during the allopatric divergence preceding a secondary contact. In hybrid zones, a further source of variance is differential introgression among chromosomal regions after the secondary contact owing to the local influence of selected genes on more or less distant markers. However, the degree of differentiation alone gives no way to distinguish indirect pseudo-selection (a regular and ubiquitous feature of hybrid zones) from direct selection. More generally, we suggest that comparative neutrality tests based on discrepancies among marker types have to be applied with caution when the presence of semi-permeable genetic barriers to gene exchange is suspected.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

Differential barrier strength and allele frequencies in hybrid zones maintained by sex-biased hybrid incompatibilities

In animals, hybrid sterility and inviability between closely related species often affect only the heterogametic sex (XY). This widespread phenomenon, known as Haldane's rule, is an early speciation event found across broad taxa, but the role of heterogametic hybrid incompatibilities, as opposed to homogametic ones, as a barrier in a speciation process remains obscure. It has been hypothesized that heterogametic incompatibility may be a more efficient mechanism in driving speciation. The population dynamics after (rather than before) the occurrence of sex-biased incompatibilities may account for Haldane's rule. In this study, a recursion model of hybrid zones was developed to investigate the differences between heterogametic and homogametic incompatibilities. The selection strengths and selection patterns of sex chromosome-linked, two-locus Bateson-Dobzhansky-Muller (BDM) genetic incompatibilities were examined. It is noted that a sex-biased hybrid incompatibility in a hybrid zone confers asymmetric and uneven impedance to gene flow. The clines of different loci in such a hybrid zone displayed diverse differentiation in their width, steepness and asymmetry. Alleles involved in the incompatibility face much stronger resistance to cross a hybrid zone. Different sex-biased BDM incompatibilities also affect the flow of neutral alleles differently. Compared to a homogametic one, heterogametic incompatibility is a weaker but more asymmetric barrier. These unique patterns of gene flow may explain uneven divergence among different genomic regions during speciation between some closely related species.     

Analysis of molecular differentiation in a hybrid zone between chromosomally distinct races of the common shrew Sorex araneus (Insectivora: Soricidae) suggests their common ancestry

During postglacial colonization, populations that diverged in different refugia produced a patchwork of genomes, often delimited with sharp hybrid zones. The outcome of hybridization following the secondary contact of two genetically distinct populations is hard to predict. In this context, the present study investigated the genetic structure of the hybrid zone between the Drnholec and Białowieża chromosome races of the common shrew ( Sorex araneus ) in Poland using biparentally inherited (seven autosomal microsatellites) and uniparentally inherited (Y-linked microsatellite and mtDNA) molecular markers. On the basis of diagnostic chromosomes, the Drnholec and Białowieża races were classified to different karyotypic groups, which were believed to have independent glacial histories. It was found that genetic differentiation between the Drnholec and Białowieża races was weak and nonsignificant with respect to all molecular markers. However, these results are in contrast with the chromosomal structure of this hybrid zone. The very sharp frequency clines of the diagnostic chromosomes strongly suggest that gene flow between the Drnholec and Białowieża races was reduced. Nonsignificant correlations between genetic differentiation and both the presence of an environmental barrier and geographical distance reveal that only differences in karyotypes might be a reason for limited gene exchange between the races. It is assumed that a lack of molecular differences between the Drnholec and Białowieża races results from a shared ancestral variation.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 79–90.     

Genetic analysis of a contact zone between two lineages of the ocellated lizard (Lacerta lepida Daudin 1802) in south‐eastern Iberia reveal a steep and narrow hybrid zone

Measuring the diffusion of genes between diverging taxa through zones of secondary contact is an essential step to understand the extent and nature of the reproductive isolation that has been achieved. Previous studies have shown that the ocellated lizard (Lacerta lepida Daudin, 1802) has endured repeated range fragmentation associated with the climatic oscillations of the Plio‐Pleistocene that promoted diversification of many different evolutionary units within the species. However, the oldest divergence within the group is estimated to have occurred much earlier, during the Miocene, around 9 Ma and corresponds to the split between the subspecies Lacerta lepida nevadensis Buchholz (1963) and Lacerta lepida lepida Daudin (1802). Although these two evolutionary units have documented genetic and morphological differentiation, most probably accumulated during periods of allopatry, little is known about patterns of gene flow between them. In this study, we performed a population genetic analysis of a putative area of secondary contact between these two taxa, using mtDNA and microsatellite data. We assessed levels of gene flow across the contact zone to clarify to what extent gene flow may be occurring. Hybridization between the subspecies was observed by the presence of genetically introgressed individuals. However, the overall coincidence of mitochondrial and multilocus nuclear clines and generally steep clines support the idea that this contact zone is acting as a barrier to gene flow. Taken together, these results suggest that L. l. lepida and L. l. nevadensis are in independent evolutionary trajectories and should be considered as two different species.     

Cytoplasmic and nuclear markers reveal contrasting patterns of spatial genetic structure in a natural Ipomopsis hybrid zone

Spatial variation in natural selection may play an important role in determining the genetic structure of hybridizing populations. Previous studies have found that F1 hybrids between naturally hybridizing Ipomopsis aggregata and Ipomopsis tenuituba in central Colorado differ in fitness depending on both genotype and environment: hybrids had higher survival when I. aggregata was the maternal parent, except in the centre of the hybrid zone where both hybrid types had high survival. Here, we developed both maternally (cpDNA PCR-RFLP) and biparentally inherited (nuclear AFLP) species-diagnostic markers to characterize the spatial genetic structure of the natural Ipomopsis hybrid zone, and tested the prediction that the majority of natural hybrids have I. aggregata cytoplasm, except in areas near the centre of the hybrid zone. Analyses of 352 individuals from across the hybrid zone indicate that cytoplasmic gene flow is bidirectional, but contrary to expectation, most plants in the hybrid zone have I. tenuituba cytoplasm. This cytotype distribution is consistent with a hybrid zone in historical transition, with I. aggregata nuclear genes advancing into the contact zone. Further, nuclear data show a much more gradual cline than cpDNA markers that is consistent with morphological patterns across the hybrid populations. A mixture of environment- and pollinator-mediated selection may contribute to the current genetic structure of this hybrid system.     

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.     

THE EFFECT OF A PARTIAL BARRIER ON THE MOVEMENT OF A HYBRID ZONE

Computer simulations of clines (Brues, 1972; Endler, 1977) as well as theoretical arguments (Nagylaki, 1975), have shown that steps in gene frequencies will be pulled to partial barriers (areas of reduced gene flow) if they form within approximately a cline width of the partial barrier. The behavior of a hybrid zone between two chromosomal taxa (“Moreton” and “Torresian”) of the acridine grasshopper Caledia captiva in southeast Queensland has been analyzed and found to comform qualitatively with a model of altered gene flow patterns. Clines in four enzyme systems were analyzed for 1983 and 1986 along a transect across the hybrid zone. The clines have shifted towards an area of regenerating forest, while homozygote frequencies have increased at this point. This forest barrier has broken the continuity of the spatial distribution of C. captiva, forming population islands in part of the hybrid zone, and thus reducing the amount of gene flow. The distance between the barrier and the original cline is approximately of the order of a cline width, so that they would be expected to interact. Historical information suggests that the secondary contact between the “Moreton” and “Torresian” taxa occurred very recently (1844–1940), due to the intensive land-clearing activities during the European settlement.     

Preliminary analysis of a hybrid zone between two subspecies of Zootermopsis nevadensis

For species largely allopatric in distribution, zones of contact provide an opportunity for hybridization, testing grounds for species boundaries, and may result in the formation of a new species. Thus, hybrid zones have the potential to provide important insights into speciation. In this study, we performed a preliminary analysis of a hybrid zone between two subspecies of the dampwood termite, Zootermopsis nevadensis (Z. n. nuttingi and Z. n. nevadensis) near Bartle, CA, using 12 microsatellite loci and a mitochondrial gene. Fifty-seven colonies collected in 36 locations were analyzed. The analysis of genetic markers revealed a large hybrid zone approximately 104 km in width. Although stepped clines best explained the data, we are unable to rule out the existence of a mosaic hybrid zone. We inferred weak selection (s* < 0.05%) against hybrids, but the data also suggested the existence of a barrier to gene flow from Z. n. nevadensis to Z. n. nuttingi, but not in the other direction. Given the large zone of contact, extensive sampling is needed to obtain a more complete characterization of the hybrid zone. However, the results of this study suggest that despite the accumulation of phenotypic and genetic differences and time since divergence (~2 million years), Z. n. nuttingi and Z. n. nevadensis are capable of extensive hybridization.