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.     

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.     

Patterns of introgression across an expanding hybrid zone: analysing historical patterns of gene flow using nonequilibrium approaches

Previous studies suggest that the hybrid zone between two taxa in the Piriqueta caroliniana complex in central Florida (south-eastern North America) has recently expanded with hybrids replacing parental genotypes across a broad region of the Florida peninsula. Here I use patterns of genetic disequilibria and levels of differentiation among populations to infer historical patterns of introgression across this broad hybrid zone. There were strong positive associations among taxon-specific alleles at the southern extent of hybridization, but disequilibria values were close to zero across the central and northern portions of the hybrid zone. Levels of among-population differentiation in the central portion of the hybrid zone were relatively low, and increased towards the northern extent of hybridization. The high levels of disequilibria are coincident with the sharp clines at the southern end of the hybrid zone, suggesting that there is a tension zone in this region that is maintained by selection against hybrids and dispersal from parental regions. The levels of disequilibria within populations and patterns of differentiation among populations are consistent with historical introgression and northward expansion of this hybrid zone, which may have slowed or ceased in recent generations.     

Ecological differentiation and habitat unsuitability maintaining a ground beetle hybrid zone

Exogenous selection via interactions between organisms and environments may influence the dynamics of hybrid zones between species in multiple ways. Two major models of a hybrid zone allowed us to hypothesize that environmental conditions influence hybrid zone dynamics in two ways. In the first model, an environmental gradient determines the mosaic distribution at the boundary between ecologically differentiated species (mosaic hybrid zone model). In the second model, a patch of unsuitable habitat traps a hybrid zone between species whose hybrids are unfit (tension zone model). To test these, we examined the environmental factors influencing the spatial structure of a hybrid zone between the ground beetles Carabus maiyasanus and C. iwawakianus using GIS‐based quantification of environmental factors and a statistical comparison of species distribution models (SDMs). We determined that both of the hypothetical processes can be important in the hybrid zone. We detected interspecific differences in the environmental factors in presence localities and their relative contribution in SDMs. SDMs were not identical between species even within contact areas, but tended to be similar within the range of each species. These results suggest an association between environments and species, and provide evidence that ecological differentiation between species plays a role in the maintenance of the hybrid zone. Contact areas were characterized by a relatively high temperature, low precipitation, and high topological wetness. Thus, the contact areas were regarded as being located in an unsuitable habitat with a drier climate, where those populations are likely to occur in patches with limited precipitation concentrated. A comparison of spatial scales suggests that exogenous selection via environmental factors may be weaker than endogenous selection via genitalic incompatibility.     

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.     

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.     

Populations composed entirely of hybrid colonies: bidirectional hybridization and polyandry in harvester ants

In eusocial Hymenoptera, haplodiploid life cycles, obligate sterile castes, and polyandry may facilitate selection for hybridization. We analyzed a broad hybrid zone between the ecologically distinct seed‐harvester ants Pogonomyrmex occidentalis (Cresson) and Pogonomyrmex maricopa (Wheeler) using mitochondrial (mt)DNA sequence data, eight morphological markers, and 14 random amplified polymorphic DNA (RAPD) markers. Average mtDNA sequence divergence among parental species was 11.34%, indicating secondary contact. RAPD markers were significantly correlated with morphological variation, confirming the interspecific hybrid origin of all morphologically putative hybrid colonies. A morphological hybrid index indicates an abundance of both F₁ hybrids and parental morphotypes within colonies. Individual character frequencies plotted against distance show coincident and concordant clines, suggesting little to no introgression. The structure of the hybrid zone is two‐fold. Within the western region, stark reversals in character frequencies coincide with overt soil differences, indicating a mosaic hybrid zone structure. The eastern region is a riparian habitat where four adjacent populations were composed entirely of hybrid colonies. These habitat associations suggest that hybrid worker genomes permit dispersal into intermediate environments that select against one or both parental species. The present study suggests that, in addition to retaining reproductive compatibility, ecologically distinct species of ants may generate hybrid colonies maintained by environmental selection. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 320–336.     

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.     

GENOMIC AND MORPHOLOGICAL ANALYSIS OF A SEMIPERMEABLE AVIAN HYBRID ZONE SUGGESTS ASYMMETRICAL INTROGRESSION OF A SEXUAL SIGNAL

Hybrid zones are geographic regions where differentiated taxa meet and potentially exchange genes. Increasingly, genomic analyses have demonstrated that many hybrid zones are semipermeable boundaries across which introgression is highly variable. In some cases, certain alleles penetrate across the hybrid zone in only one direction, recombining into the alternate genome. We investigated this phenomenon using genomic (genotyping‐by‐sequencing) and morphological (plumage reflectance spectrophotometry) analyses of the hybrid zone between two subspecies of the red‐backed fairy‐wren (Malurus melanocephalus) that differ conspicuously in a sexual signal, male back plumage color. Geographic cline analyses revealed a highly variable pattern of differential introgression, with many narrow coincident clines combined with several significantly wider clines, suggesting that the hybrid zone is a semipermeable tension zone. The plumage cline was shifted significantly into the genomic background of the orange subspecies, consistent with sexual selection driving asymmetrical introgression of red plumage alleles across the hybrid zone. This interpretation is supported by previous experimental work demonstrating an extra‐pair mating advantage for red males, but the role of genetic dominance in driving this pattern remains unclear. This study highlights the potential for sexual selection to erode taxonomic boundaries and promote gene flow, particularly at an intermediate stage of divergence.     

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

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

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.     

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long‐distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present‐day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.     

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological,genetic and climatic data to examine stability

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south‐eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90‐km transect through the hybrid zone and examined variation in 26 microsatellites, the D‐loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection‐gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.     

A dispersal-dependent zone of introgressive hybridization between weakfish, Cynoscion regalis, and sand seatrout, C. arenarius, (Sciaenidae) in the Florida Atlantic

Five diagnostic codominant nuclear DNA markers and a diagnostic mitochondrial DNA marker were used to survey weakfish (Cynoscion regalis) and sand seatrout (C. arenarius), with particular focus on heretofore uncharacterized juvenile populations along the Florida (FL) Atlantic coast. Geographic and reproductive ranges of weakfish and sand seatrout were shown to overlap on the Atlantic coast along north and central FL. An active bidirectional zone of introgressive hybridization exists between these taxa, centered in the St Johns River, FL. Strong patterns of Hardy-Weinberg, linkage, and cytonuclear disequilibrium and a bimodal hybrid index distribution were observed for juvenile cohorts in the zone center, coupled with narrow (～240 km) concordant clines. Parental forms had disparate habitat preferences; hybrid forms occurred predominantly in intermediate habitats. All genetic data were consistent with the hypothesis that the C. arenarius-C. regalis hybrid zone is maintained by a dynamic equilibrium between continued interspecific gene flow and one or more opposing forces. Cytonuclear analyses indicated that parental forms mate assortatively in the zone but that mate recognition was imperfect. Ethological mating dynamics are likely stabilized by some form of endogenous or exogenous postfertilization selection against hybrids such that parental taxa will likely continue to evolve independently.     

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.     

Concordance of genetic and phenotypic characters across a sapsucker hybrid zone

Hybridization has presented a challenge for taxonomists and conservation biologists, since hybridizing forms could be stable evolutionary entities or ephemeral forms that are blending together. However, hybrid zones also provide a unique opportunity for evolutionary biologists who study the interaction between gene flow and reproductive isolation in speciation. Three forms of woodpeckers (sapsuckers; genus Sphyrapicus) in North America that are mostly geographically separated but hybridize with each other where they come into contact present a remarkable system for the study of hybridization. We provide the first comprehensive analysis of phenotypic and genetic variation across a hybrid zone between two of these forms, the red‐breasted Sphyrapicus ruber and yellow‐bellied S. varius sapsuckers. The objective was to infer whether selection maintains the differences between forms. Our analysis of eight morphometric and 20 plumage traits, and two molecular markers showed clear differences between the forms and roughly concordant clinal variation across a narrow hybrid zone. Thirty percent of sampled birds in the hybrid zone had mixed west/east genotypes at the genetic markers examined. The center of the genetic cline was located 20 km west of the crest of the Rocky Mountains. The width of the zone was 122 km, narrower than would be expected under neutral blending given reasonable estimates of the age of the zone and individual dispersal distances. Heterozygote deficit and cytonuclear disequilibrium at the centre of the hybrid zone suggested nonrandom mating or limited hybridization. Given these patterns and lack of evidence for habitat segregation we conclude that this hybrid zone is maintained by selection, most likely in the form of hybrid inferiority. This study provides an illustrative example of extensive hybridization between stable entities, providing additional evidence against the historical practice of treating hybridizing forms as members of the same species.     

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.     

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.