Narrow hybrid zones in spite of very low population differentiation in neutral markers in an island bird species complex

Patterns of phenotypic and genic frequencies across hybrid zones provide insight into the origin and evolution of reproductive isolation. The Reunion grey white‐eye, Zosterops borbonicus, exhibits parapatrically distributed plumage colour forms across the lowlands of the small volcanic island of Reunion (Mascarene archipelago). These forms meet and hybridize in regions that are natural barriers to dispersal (rivers, lava fields). Here, we investigated the relationship among patterns of differentiation at neutral genetic (microsatellite) markers, phenotypic traits (morphology and plumage colour) and niche characteristics across three independent hybrid zones. Patterns of phenotypic divergence revealed that these hybrid zones are among the narrowest ever documented in birds. However, the levels of phenotypic divergence stand in stark contrast to the lack of clear population neutral genetic structure between forms. The position of the hybrid zones coincides with different natural physical barriers, yet is not associated with steep changes in vegetation and related climatic variables, and major habitat transitions are shifted from these locations by at least 18 km. This suggests that the hybrid zones are stabilized over natural dispersal barriers, independently of environmental boundaries, and are not associated with niche divergence. A striking feature of these hybrid zones is the very low levels of genetic differentiation in neutral markers between forms, suggesting that phenotypic divergence has a narrow genetic basis and may reflect recent divergence at a few linked genes under strong selection, with a possible role for assortative mating in keeping these forms apart.     

Environmental factors, microgeographic patterns of endemism and hybrid zones in Monobella grassei (Insecta: Collembola: Neanuridae)

A complex pattern of narrow parapatric distributions with hybrid zones is described among subspecies of the Collembola species Monobella grassei in the eastern Pyrenees. Though hybrid zones are usually narrow, they vary from simple contact zones without recorded hybrids to relatively wide hybrid swarms, suggesting that large differences exist locally in the level of genetic divergence between subspecies. The course and structure of contact zones are analysed on morphological grounds, in relation to environmental (rainfall, vegetation, substrate) and topographical (altitude, hydrography) parameters. The results show that, unlike most cases documented so far, the courses of the contact zones in M. grassei are not connected to major ecological boundaries nor to ecological ecotones, on either the regional or the habitat scale. It is suggested that this is due to the great ecological similarity between the different parapatric subspecies of M. grassei , and their broad ecological requirements. There is no conclusive evidence that the distribution observed is linked to recent palaeoclimatic events such as glaciation. Endogenous selection could therefore be prevalent in this model to explain parapatry and the hybrid zone pattern.     

Hybrid zones between chromosome races of the common shrew in West Siberia

In West Siberia, the whole species range of the common shrew (Sorex araneus L.) is shared by two parapatric chromosome races — Novosibirsk and Tomsk. These races form a hybrid zone with each other. In addition, on the western margin of the range there is a hybrid zone between the Novosibirsk race and the Uralian race Serov, and in the east, the Tomsk race forms a hybrid zone with the East Siberian race Strelka hybrid. The structures of the three hybrid zones are very different and depend on the karyotypic state of the races in contact. A comparative analysis of the hybrid zones between the chromosome races of the common shrew in West Siberia is discussed together with the role of these zones in maintaining variability in natural populations.     

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

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

Species distribution models contribute to determine the effect of climate and interspecific interactions in moving hybrid zones

Climate is a major factor delimiting species’ distributions. However, biotic interactions may also be prominent in shaping geographical ranges, especially for parapatric species forming hybrid zones. Determining the relative effect of each factor and their interaction of the contact zone location has been difficult due to the lack of broad scale environmental data. Recent developments in species distribution modelling (SDM) now allow disentangling the relative contributions of climate and species’ interactions in hybrid zones and their responses to future climate change. We investigated the moving hybrid zone between the breeding ranges of two parapatric passerines in Europe. We conducted SDMs representing the climatic conditions during the breeding season. Our results show a large mismatch between the realized and potential distributions of the two species, suggesting that interspecific interactions, not climate, account for the present location of the contact zone. The SDM scenarios show that the southerly distributed species, Hippolais polyglotta, might lose large parts of its southern distribution under climate change, but a similar gain of novel habitat along the hybrid zone seems unlikely, because interactions with the other species (H. icterina) constrain its range expansion. Thus, whenever biotic interactions limit range expansion, species may become ‘trapped’ if range loss due to climate change is faster than the movement of the contact zone. An increasing number of moving hybrid zones are being reported, but the proximate causes of movement often remain unclear. In a global context of climate change, we call for more interest in their interactions with climate change.     

Identifying contact zone hotspots of passerine birds in the Palaearctic region

Birds are known to be a group rich in pairs of closely related species that have parapatric or allopatric distributions with relatively narrow contact zones. Here we analyse the geographical distribution of these contact zones for parapatric species pairs of passerine birds of the Palaearctic region. Their contact zones are located mainly in southwestern, northern and central-southern Asia, and in northwestern Africa, with a hotspot in the Middle East. A mid-domain effect null model, where contact zone hotspots are a neutral correlate of continental geometry, had a low explanatory power of 3.8%; the observed distribution of contact zones was not sufficiently predicted. Hypotheses involving range contractions and secondary contact in areas of high topographic and habitat diversity may offer more convincing explanations and offer promising perspectives for future studies.     

Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time

Abstract Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds ( B. bombina ) or puddles ( B. variegata ), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.     

Geographic variation in the structure of oak hybrid zones provides insights into the dynamics of speciation

Studying geographic variation in the rate of hybridization between closely related species could provide a useful window on the evolution of reproductive isolation. Reinforcement theory predicts greater prezygotic isolation in areas of prolonged contact between recently diverged species than in areas of recent contact, which implies that old contact zones would be dominated by parental phenotypes with few hybrids (bimodal hybrid zones), whereas recent contact zones would be characterized by hybrid swarms (unimodal hybrid zones). Here, we investigate how the hybrid zones of two closely related Chinese oaks, Quercus mongolica and Q. liaotungensis, are structured geographically using both nuclear and chloroplast markers. We found that populations of Q. liaotungensis located around the Changbai Mountains in Northeast China, an inferred glacial refugium, were introgressed by genes from Q. mongolica, suggesting historical contact between the two species in this region. However, these introgressed populations form sharp bimodal hybrid zones with Q. mongolica. In contrast, populations of Q. liaotungensis located in North China, which show no sign of ancient introgression with Q. mongolica, form unimodal hybrid zones with Q. mongolica. These results are consistent with the hypothesis that selection against hybrids has had sufficient time to reinforce the reproductive barriers between Q. liaotungensis and Q. mongolica in Northeast China but not in North China.     

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.     

When invaders meet locally adapted types: rapid moulding of hybrid zones between sculpins (Cottus, Pisces) in the Rhine system

Hybrid zones are commonly studied to dissect the processes that drive divergence among lineages, which have incomplete barriers of reproduction. Most hybrid zones have existed for an extended time making inferences on the initial mode of formation difficult. It is a priori unclear how fast a nascent hybrid zone would form as a response to endogenous and exogenous factors. We have studied several hybrid zones between two lineages of sculpins (Cottus spp.), which emerged due to a recent range expansion of one of the lineages along the river Rhine in the early 1990s. Applying a dense sampling across two contact areas and using a highly informative set of 45 microsatellite markers we found pronounced genetic structure. Steep genetic clines suggest that strong selective forces have shaped the respective hybrid zones from the beginning. We find that the zones are coupled to ecological transitions from small streams to larger rivers. The width of these zones is much smaller than estimates of annual individual dispersal distances, as estimated outside of the hybrid zones. The pattern is apparently not strongly affected by pre- or postzygotic reproductive isolation because numerous backcross hybrids occur within the zones. This suggests that strong natural selection acts against immigrant genotypes. The study exemplifies how local adaptation can play a key role in preventing admixture in dependence of the ecological context.     

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental‐like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.     

Hybridization during altitudinal range shifts: nuclear introgression leads to extensive cyto‐nuclear discordance in the fire salamander

Ecological models predict that, in the face of climate change, taxa occupying steep altitudinal gradients will shift their distributions, leading to the contraction or extinction of the high‐elevation (cold‐adapted) taxa. However, hybridization between ecomorphologically divergent taxa commonly occurs in nature and may lead to alternative evolutionary outcomes, such as genetic merger or gene flow at specific genes. We evaluate this hypothesis by studying patterns of divergence and gene flow across three replicate contact zones between high‐ and low‐elevation ecomorphs of the fire salamander (Salamandra salamandra) that have experienced altitudinal range shifts over the current postglacial period. Strong population structure with high genetic divergence in mitochondrial DNA suggests that vicariant evolution has occurred over several glacial–interglacial cycles and that it has led to cryptic differentiation within ecomorphs. In current parapatric boundaries, we do not find evidence for local extinction and replacement upon postglacial expansion. Instead, parapatric taxa recurrently show discordance between mitochondrial and nuclear markers, suggesting nuclear‐mediated gene flow across contact zones. Isolation with migration models support this hypothesis by showing significant gene flow across all five parapatric boundaries. Together, our results suggest that, while some genomic regions, such as the mitochondria, may follow morphologic species traits and retreat to isolated mountain tops, other genomic regions, such as nuclear markers, may flow across parapatric boundaries, sometimes leading to a complete genetic merger. We show that despite high ecologic and morphologic divergence over prolonged periods of time, hybridization allows for evolutionary outcomes alternative to extinction and replacement of taxa in response to climate change.     

Distinctiveness in the face of gene flow: hybridization between specialist and generalist gartersnakes

Patterns of divergence and polymorphism across hybrid zones can provide important clues as to their origin and maintenance. Unimodal hybrid zones or hybrid swarms are composed predominantly of recombinant individuals whose genomes are patchworks of alleles derived from each parental lineage. In contrast, bimodal hybrid zones contain few identifiable hybrids; most individuals fall within distinct genetic clusters. Distinguishing between hybrid swarms and bimodal hybrid zones can be important for taxonomic and conservation decisions regarding the status and value of hybrid populations. In addition, the causes of bimodality are important in understanding the generation and maintenance of biological diversity. For example, are distinct clusters mostly reproductively isolated and co‐adapted gene complexes, or can distinctiveness be maintained by a few ‘genomic islands’ despite rampant gene flow across much of the genome? Here we focus on three patterns of distinctiveness in the face of gene flow between gartersnake taxa in the Great Lakes region of North America. Bimodality, the persistence of distinct clusters of genotypes, requires strong barriers to gene flow and supports recognition of distinct specialist (Thamnophis butleri) and generalist (Thamnophis radix) taxa. Concordance of DNA‐based clusters with morphometrics supports the hypothesis that trophic morphology is a key component of divergence. Finally, disparity in the level of differentiation across molecular markers (amplified fragment length polymorphisms) indicates that distinctiveness is maintained by strong selection on a few traits despite high gene flow currently or in the recent past.     

Some genetic consequences of ice ages, and their role, in divergence and speciation

The genetic effects of pleistocene ice ages are approached by deduction from paleoenvironmental information, by induction from the genetic structure of populations and species, and by their combination to infer likely consequences. (1) Recent palaeoclimatic information indicate rapid global reversals and changes in ranges of species which would involve elimination with spreading from the edge. Leading edge colonization during a rapid expansion would be leptokurtic and lead to homozygosity and spatial assortment of genomes. In Europe and North America, ice age contractions were into southern refugia, which would promote genome reorganization. (2) The present day genetic structure of species shows frequent geographic subdivision, with parapatric genomes, hybrid zones and suture zones. A survey of recent DNA phylogeographic information supports and extends earlier work. (3) The grasshopper Chorthippus parallelus is used to illustrate such data and processes. Its range in Europe is divided on DNA sequences into five parapatric races, with southern genomes showing greater haplotype diversity - probably due to southern mountain blocks acting as refugia and northern expansion reducing diversity. (4) Comparison with other recent studies shows a concordance of such phylogeographic data over pleistocene time scales. (5) The role that ice age range changes may have played in changing adaptations is explored, including the limits of range, rapid change in new invasions and refugial differentiation in a variety of organisms. (6) The effects of these events in causing divergence and speciation are explored using Chorthippus as a paradigm. Repeated contraction and expansion would accumulate genome differences and adaptations, protected from mixing by hybrid zones, and such a composite mode of speciation could apply to many organisms.     

The maintenance of hybrid zones across a disturbance gradient

The parapatric distribution of genetically divergent lineages in hybrid zones can be maintained by ecological differences (dispersal-independent 'ecotonal' hybrid zones), by frequency- and density-dependent interference when they intermingle and mate (dispersal-dependent 'tension' hybrid zones), or by both processes acting together. One potentially important ecological factor that has received little theoretical attention is gradients in habitat disturbance. Such gradients may be particularly important in contact zones in which the interacting lineages differ in their sexual system (e.g., self-fertile versus obligately outcrossing) because self-fertility promotes the colonization of open patches. Here we use a spatially explicit metapopulation model to examine the dynamics of a dispersal-dependent ecotonal hybrid zone across a gradient in the rate of habitat disturbance, where competing lineages differ in their sexual system. We found that self-fertility promoted the maintenance of one lineage over its outcrossing counterpart at high extinction rates, predominantly because self-fertility confers reproductive assurance. Additionally, greater seed and pollen production promoted a lineage's persistence by reducing the seed fertility of its counterpart through hybridization. Our results draw attention to the joint effects of ecological and endogenous selection in regulating the location of hybrid zones. Our study also casts new light on the maintenance of the parapatric distribution of incompatible lineages of Spanish populations of the plant Mercurialis annua. In particular, we expect the rate of movement of a contact zone in eastern Spain to increase as it moves further south, contrary to earlier predictions.     

A hybrid zone dominated by fertile F1s: maintenance of species barriers in Rhododendron

Isolating barriers between interbreeding sympatric or parapatric interfertile species are maintained by processes that occur within their hybrid zones. Although the effects of intrinsic selection on hybrid fitness are well known, less is understood about extrinsic fitness variation. At Tiryal Dag, northeast Turkey, Rhododendron x sochadzeae (R. ponticum x caucasicum) forms large populations in which neither segregation nor backcrossing occur, in habitats intermediate between those of its parents. Using single-copy species-specific random amplified polymorphic DNA and inter simple sequence repeat markers, it was determined that most or all R. x sochadzeae plants are F1s, and that there are many separate genets present. Hand pollination and germination experiments showed that R. x sochadzeae plants can produce viable seed of F2s or backcrosses in either direction. Furthermore, adult backcrosses have been observed in habitats atypical for R. x sochadzeae. From this, all non-F1 hybrid derivatives appear to be eliminated in the hybrid zone at Tiryal Dag as a result of postgermination selection. This absence of post-F1 hybrid derivatives apparently prevents introgression. This type of hybrid population is here termed an F1-dominated hybrid zone (F1DZ), and also occurs in Encelia. The observed dominance of F1s within a narrow habitat range is best explained by habitat-mediated superiority of F1s over all other genotype classes. Therefore, habitat-mediated selection against the second hybrid generation might be preventing interspecific gene flow in R. x sochadzeae. F1DZ formation is postulated to require the formation of F1s in quantity, habitat-mediated superiority in F1s, and highly specific habitat conditions.     

Homoploid hybridization of plants in the Hengduan mountains region

The Hengduan Mountains Region (HMR) is a major global biodiversity hotspot. Complex tectonic and historical climatic conditions created opportunities for natural interspecific hybridization. Likewise, anthropogenic disturbance potentially raises the frequency of hybridization. Among species studies to date, the frequency of homoploid hybridization appears in the HMR. Of nine taxa in which natural hybridization has been detected, three groups are involved in homoploid hybrid speciation, and species pairs from the remaining six genera suggest that continuous gene flow occurs in hybrid zones. Reproductive isolation may greatly affect the dynamic and architecture of hybrid zones in the HMR. Asymmetrical hybridization and introgression can primarily be attributed to both prezygotic and postzygotic barriers. The frequent observation of such asymmetry may imply that reproductive barrier contributes to maintaining species boundaries in the alpine region. Ecological isolations with environmental disturbance may promote breeding barriers between parental species and hybrids. Hybrid zones may be an important phase for homoploid hybrid speciation. Hybrid zones potentially provided abundant genetic resources for the diversification of the HMR flora. The ecological and molecular mechanisms of control and mediation for natural hybridization will help biologists to understand the formation of biodiversity in the HMR. More researches from ecological and molecular aspects were required in future studies.     

Cytonuclear genetic structure of a hybrid zone in lizards of the Sceloporus grammicus complex (Sauria, Phrynosomatidae)

Lizards of the Sceloporus grammicus complex are comprised of multiple chromosome races that form several zones of parapatric hybridization in central Mexico. We scored diagnostic mitochondrial DNA (mtDNA) haplotypes and autosomal chromosome markers in a sample of 342 lizards from one well-defined zone between 2n = 34 and 2n = 46 races. A two-part analysis was performed on this data set in an attempt to infer the predominant evolutionary forces shaping the cytonuclear structure of this zone. The complications posed by its spatial structure were addressed by analysing a hierarchical series of smaller subsamples chosen to approximate single mating units. Two critical conclusions were drawn from this first-stage analysis. First and foremost, the three chromosomes have largely concordant cytonuclear disequilibrium patterns within each subsample with adequate numbers of individuals for detecting nonrandom cytonuclear associations. This suggests that the cytonuclear structure of this zone is predominantly a result of deterministic genome-wide forces rather than genetic drift or deterministic forces specific to individual chromosomes or loci. Second, the fit of a series of migration models to the data shows that the cytonuclear structure of the subsamples is well accounted for by continued gene flow from the two parental races alone, with random mating with respect to cytonuclear genotype and no other evolutionary forces. These results motivate several further empirical and theoretical investigations to refine our understanding of the relative roles of migration and other potentially important forces such as natural selection and genetic drift, in this and other hybrid zones.     

INCIPIENT SPECIATION DESPITE LITTLE ASSORTATIVE MATING: THE YELLOW‐RUMPED WARBLER HYBRID ZONE

Hybrid zones between recently diverged taxa are natural laboratories for speciation research, allowing us to determine whether there is reproductive isolation between divergent forms and the causes of that isolation. We present a study of a classic avian hybrid zone in North America between two subspecies of the yellow‐rumped warbler (Dendroica coronata). Although previous work has shown very little differentiation in mitochondrial DNA across this hybrid zone, we identified two nuclear loci (one sex‐linked and one autosomal) that show fixed differences across the hybrid zone, in a close concordance with patterns of plumage variation. Temporal stability and limited width of the hybrid zone, along with substantial linkage disequilibrium between these two diagnostic markers in the center of the zone, indicate that there is moderate reproductive isolation between these populations, with an estimated strength of selection maintaining the zone of 18%. Pairing data indicate that assortative mating is either very weak or absent, suggesting that this reproductive isolation is largely due to postmating barriers. Thus, despite extensive hybridization the two forms are distinct evolutionary groups carrying genes for divergent adaptive peaks, and this situation appears relatively stable.     

Structure of a mosaic hybrid zone between the field crickets Gryllus firmus and G. pennsylvanicus

Hybrid zones provide insight into the nature of species boundaries and the evolution of barriers to gene exchange. Characterizing multiple regions within hybrid zones is essential for understanding both their history and current dynamics. Here, we describe a previously uncharacterized region of a well‐studied hybrid zone between two species of field crickets, Gryllus pennsylvanicus and G. firmus. We use a combination of mitochondrial DNA sequencing, morphological data, and modeling of environmental variables to identify the ecological factors structuring the hybrid zone and define patterns of hybridization and introgression. We find an association between species distribution and natural habitat; Gryllus pennsylvanicus occupies natural habitat along forest edges and natural clearings, whereas G. firmus occupies more disturbed areas in agricultural and suburban environments. Hybridization and introgression occur across patch boundaries; there is evidence of substantial admixture both in morphological characters and mtDNA, over a broad geographic area. Nonetheless, the distribution of morphological types is bimodal. Given that F₁ hybrids are viable and fertile in the lab, this suggests that strong pre‐zygotic barriers are operating in this portion of the hybrid zone.