Genetic structure of hybrid zones between Pinus pumila and P. parviflora var. pentaphylla (Pinaceae) revealed by molecular hybrid index analysis

Pinus species have three differently inherited genomes: paternal chloroplast, maternal mitochondrial, and biparental nuclear. Our previous study on the hybrid zones between alpine Pinus pumila and montane to subalpine P. parviflora var. pentaphylla demonstrated contrasting patterns of introgression of two cytoplasmic genomes, i.e., the paternal cpDNA flowed from P. parviflora var. pentaphylla to P. pumila, and the maternal mtDNA flowed in the reverse direction. In the present study, we developed codominant nuclear DNA markers diagnostic or mostly diagnostic for each parental species by single-strand conformation polymorphism (SSCP) of polymerase chain reaction (PCR) products, using expressed sequence tag (EST) primers of Pinus taeda. To describe the introgressive patterns of the nuclear genes, the molecular hybrid index (MHI) showing the overall proportion of alleles inferred to be derived from P. pumila was determined for each plant collected in hybrid zones on Mt. Asahidake and Mt. Higashiazuma, Japan. At Mt. Asahidake, the MHI values changed clinally according to the altitudes at which the plants were collected. However, at Mt. Higashiazuma, there was a gap in the MHI values between the plants above and below the Abies and Tsuga forest zone (alt. 1800-1900 m). This suggested that the zone plays a role in creating an effective barrier to gene flow in the hybrid zone.     

Reproductive isolation and hybrid pollen disadvantage in Ipomopsis

One cause of reproductive isolation is gamete competition, in which conspecific pollen has an advantage over heterospecific pollen in siring seeds, thereby decreasing the formation of F1 hybrids. Analogous pollen interactions between hybrid pollen and conspecific pollen can contribute to post-zygotic isolation. The herbaceous plants Ipomopsis aggregata and I. tenuituba frequently hybridize in nature. Hand-pollination of I. aggregata with pollen from F1 or F2 hybrids produced as many seeds as hand-pollination with conspecific pollen, suggesting equal pollen viability. However, when mixed pollen loads with 50% conspecific pollen and 50% hybrid pollen were applied to I. aggregata stigmas, fewer than half of the seeds had hybrid sires. Such pollen mixtures are frequently received if plants of the two species and F1 and F2 hybrids are intermixed, suggesting that this advantage of conspecific over hybrid pollen reduces backcrossing and contributes to reproductive isolation.     

Spatial distribution of cpDNA and mtDNA haplotypes in a hybrid zone betweenPinus pumila andP. parviflora var.pentaphylla (Pinaceae)

Pinus species exhibit paternal chloroplast inheritance and maternal mitochondrial inheritance. This independent inheritance of two cytoplasmic genomes provides an exceptional environment for discriminating female (seeds) and male (pollen) components of gene flow across hybridizing species. We obtained mitochondrial genetic markers diagnostic toP. parviflora var.pentaphylla andP. pumila by PCR amplification of the intron ofnad1 on mtDNA, and examined the spatial-distribution pattern of the mtDNA haplotypes in a hybrid zone betweenP. parviflora var.pentaphylla andP. pumila in the Tanigawa Mountains of Japan. These data, in conjunction with previous information on cpDNA haplotypes and needle morphology, revealed contrastive patterns of introgression of two cytoplasmic genomes. CpDNA introgression has occurred uni-directionally fromP. parviflora var.pentaphylla toP. pumila. Conversely, mtDNA introgression has occurred in the opposite direction, fromP. pumila toP. parviflora var.pentaphylla. Levels of introgression are roughly equivalent for cpDNA and mtDNA. The contrastive spatial distribution pattern of cpDNA and mtDNA haplotypes could be caused by differential movement of seeds and pollen for interspecific genetic exchange.     

Secondary contact between Lycaeides idas and L. melissa in the Rocky Mountains: extensive admixture and a patchy hybrid zone

Studies of hybridization have increased our understanding of the nature of species boundaries, the process of speciation, and the effects of hybridization on the evolution of populations and species. In the present study we use genetic and morphological data to determine the outcome and consequences of secondary contact and hybridization between the butterfly species Lycaeides idas and L. melissa in the Rocky Mountains. Admixture proportions estimated from structure and geographical cline analysis indicate L. idas and L. melissa have hybridized extensively in the Rocky Mountains and that reproductive isolation was insufficient to prevent introgression for much of the genome. Geographical patterns of admixture suggest that hybridization between L. idas and L. melissa has led to the formation of a hybrid zone. The hybrid zone is relatively wide, given estimates of dispersal for Lycaeides butterflies, and does not show strong evidence of cline concordance among characters. We believe the structure of the Lycaeides hybrid zone might be best explained by the patchy distribution of Lycaeides, local extinction and colonization of habitat patches, environmental variation and weak overall selection against hybrids. We found no evidence that hybridization in the Rocky Mountains has resulted in the formation of independent hybrid species, in contrast to the outcome of hybridization between L. idas and L. melissa in the Sierra Nevada. Finally, our results suggest that differences in male morphology between L. idas and L. melissa might contribute to isolation, or perhaps even that selection has favoured the spread of L. melissa male genitalia alleles.     

Absence of conspecific pollen advantage in the dynamics of an Ipomopsis (Polemoniaceae) hybrid zone

The frequency of hybrid formation in angiosperms depends on how often heterospecific pollen is transferred to the stigma and on the success of that heterospecific pollen at fertilizing ovules. Even if heterospecific pollen is capable of effecting fertilization it may perform poorly when conspecific pollen is also available on the stigma. We applied pollen mixtures to stigmas to determine how pollen interactions affect siring success and the frequency of hybrid formation between two species of Ipomopsis (Polemoniaceae) in Colorado. Plants of both parental species and natural hybrids were pollinated with I. aggregata and I. tenuituba pollen in ratios of 100:0, 80:20, 50:50, 20:80, and 0:100 by mass. Plants were homozygous for different alleles at an isozyme marker, allowing us to distinguish the type of pollen parent for 2166 viable seeds from 273 fruits. In contrast to studies of many other hybridizing taxa, there was no evidence of an advantage to conspecific pollen, nor did composition of the stigmatic pollen load affect seed set. Instead, the frequency of seeds sired by a given species was proportional to its representation in the pollen load. In this hybrid zone, both the frequency of first-generation hybrid formation and the relative male fitness of the two parental species should be predictable from the rates of pollen transfer to stigmas.     

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.     

Ants exhibit asymmetric hybridization in a mosaic hybrid zone

Research on hybridization between species provides unparalleled insights into the pre‐ and postzygotic isolating mechanisms that drive speciation. In social organisms, colony‐level incompatibilities may provide additional reproductive barriers not present in solitary species, and hybrid zones offer an opportunity to identify these barriers. Here, we use genotyping‐by‐sequencing to sequence hundreds of markers in a hybrid zone between two socially polymorphic ant species, Formica selysi and Formica cinerea. We characterize the zone, determine the frequency of hybrid workers, infer whether hybrid queens or males are produced and investigate whether hybridization is influenced by colony social organization. We also compare cuticular hydrocarbon profiles and aggression levels between the two species. The hybrid zone exhibits a mosaic structure. The asymmetric distribution of hybrids skewed towards F. cinerea suggests a pattern of unidirectional nuclear gene flow from F. selysi into F. cinerea. The occurrence of backcrossed individuals indicates that hybrid queens and/or males are fertile, and the presence of the F. cinerea mitochondrial haplotype in 97% of hybrids shows that successful F1 hybrids will generally have F. cinerea mothers and F. selysi fathers. We found no evidence that social organization contributes to speciation, because hybrids occur in both single‐queen and multiple‐queen colonies. Strongly differentiated cuticular hydrocarbon profiles and heightened interspecific aggression further reveal that species recognition cues are both present and perceived. The discovery of fertile hybrids and asymmetrical gene flow is unusual in ants, and this hybrid zone will therefore provide an ideal system with which to investigate speciation in social insects.     

Temporal cline in a hybrid zone population between Fraxinus excelsior L. and Fraxinus angustifolia Vahl

The two closely related ash species Fraxinus excelsior L. (common ash) and Fraxinus angustifolia Vahl (narrow-leaved ash) have a broad contact zone in France where they hybridize. However, little is known about the local structure of hybrid zone populations and the isolation mechanisms. We assessed the potential effect of floral phenology on the structure of a riparian ash hybrid zone population in central France. The distribution of flowering times was unimodal and lay between the flowering periods of the two species. Using microsatellite markers, we detected isolation by time, which has possibly originated from assortative mating. Multivariate analyses indicated that morphological variation is not distributed at random with respect to flowering times. Spatial autocorrelation analyses showed that temporal and spatial patterns were tightly linked. Interestingly, despite the fact that the population shows isolation by time, neighbourhood size and historical dispersal variance (sigma = 63 m) are similar to those detected in pure stands of F. excelsior where individuals flower rather synchronously and hermaphrodites are not the most frequent sexual type. Trees flowering at intermediate dates, which comprised the majority of the population, produced on average more flowers and fruits. We detected no significant differences in floral parasite infections relative to reproductive timing, although there was a tendency for late flowering trees to suffer from more gall attack. We discuss the impact of temporal variation in fitness traits and their possible role in the maintenance of the hybrid zone.     

Relative abundance and the species-specific reinforcement of male mating preference in the Chrysochus (Coleoptera: Chrysomelidae) hybrid zone

Most studies of reinforcement have focused on the evolution of either female choice or male mating cues, following the long-held view in sexual selection theory that mating mistakes are typically more costly for females than for males. However, factors such as conspecific sperm precedence can buffer females against the cost of mating mistakes, suggesting that in some hybrid zones mating mistakes may be more costly for males than for females. Thus, the historical bias in reinforcement research may underestimate its frequency. In this study, we present evidence that reinforcement has driven the evolution of male choice in a hybrid zone between the highly promiscuous leaf beetles Chrysochus cobaltinus and C. auratus, the hybrids of which have extremely low fitness. In addition, there is evidence for male choice in these beetles and that male mating mistakes may be costly, due to reduced opportunities to mate with conspecific females. The present study combines laboratory and field methods to quantify the strength of sexual isolation, test the hypothesis of reproductive character displacement, and assess the link between relative abundance and the strength of selection against hybridization. We document that, while sexual isolation is weak, it is sufficient to produce positive assortative mating. In addition, reproductive character displacement was only detected in the relatively rare species. The strong postzygotic barriers in this system are sufficient to generate the bimodality that characterizes this hybrid zone, but the weak sexual isolation is not, calling into question whether strong prezygotic isolation is necessary for the maintenance of bimodality. Growing evidence that the cost of mating mistakes is sufficient to shape the evolution of male mate choice suggests that the reinforcement of male mate choice may prove to be a widespread occurrence.     

Phenotypic and genetic introgression across a moving woodpecker hybrid zone

In hybrid zones in which two divergent taxa come into secondary contact and interbreed, selection can maintain phenotypic diversity despite widespread genetic introgression. Red‐breasted (Sphyrapicus ruber) and red‐naped (S. nuchalis) sapsuckers meet and hybridize along a narrow contact zone that stretches from northern California to southern British Columbia. We found strong evidence for changes in the structure of this hybrid zone across time, with significant temporal shifts in allele frequencies and in the proportions of parental phenotypes across the landscape. In addition to these shifts, we found that differences in plumage predict genetic differences (R² = 0.80), suggesting that plumage is a useful proxy for assessing ancestry. We also found a significant bimodal distribution of hybrids across the contact zone, suggesting that premating barriers may be driving reproductive isolation, perhaps as a result of assortative mating based on plumage differences. However, despite evidence of selection and strong patterns of population structure between parental samples, we found only weak patterns of genetic divergence. Using museum specimens and genomic data, this study of sapsuckers provides insight into the ways in which phenotypic and genetic structure have changed over a 40‐year period, as well as insight into the mechanisms that may contribute to the maintenance of the hybrid zone over time.     

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages

Diverging semi‐isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long‐snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry—i.e., in the same geographical zone—with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon‐like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus‐specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome‐wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome‐wide island of differentiation. Since Atlantic lineages do not map to lagoon‐sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts–i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.     

Assessing reproductive isolation using a contact zone between parapatric lake‐stream stickleback ecotypes

Ecological speciation occurs when populations evolve reproductive isolation as a result of divergent natural selection. This isolation can be influenced by many potential reproductive barriers, including selection against hybrids, selection against migrants and assortative mating. How and when these barriers act and interact in nature is understood for relatively few empirical systems. We used a mark–recapture experiment in a contact zone between lake and stream three‐spined sticklebacks (Gasterosteus aculeatus, Linnaeus) to evaluate the occurrence of hybrids (allowing inferences about mating isolation), the interannual survival of hybrids (allowing inferences about selection against hybrids) and the shift in lake‐like vs. stream‐like characteristics (allowing inferences about selection against migrants). Genetic and morphological data suggest the occurrence of hybrids and no selection against hybrids in general, a result contradictory to a number of other studies of sticklebacks. However, we did find selection against more lake‐like individuals, suggesting a barrier to gene flow from the lake into the stream. Combined with previous work on this system, our results suggest that multiple (most weakly and often asymmetric) barriers must be combining to yield substantial restrictions on gene flow. This work provides evidence of a reproductive barrier in lake–stream sticklebacks and highlights the value of assessing multiple reproductive barriers in natural contexts.     

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.     

The genetics of reproductive isolation and the potential for gene exchange between Drosophila pseudoobscura and D. persimilis via backcross hybrid males

Hybrid male sterility, hybrid inviability, sexual isolation, and a hybrid male courtship dysfunction reproductively isolate Drosophila pseudoobscura and D. persimilis. Previous studies of the genetic bases of these isolating mechanisms have yielded only limited information about how much and what areas of the genome are susceptible to interspecies introgression. We have examined the genetic basis of these barriers to gene exchange in several thousand backcross hybrid male progeny of these species using 14 codominant molecular genetic markers spanning the five chromosomes of these species, focusing particularly on the autosomes. Hybrid male sterility, hybrid inviability, and the hybrid male courtship dysfunction were all associated with X-autosome interactions involving primarily the inverted regions on the left arm of the X-chromosome and the center of the second chromosome. Sexual isolation from D. pseudoobscura females was primarily associated with the left arm of the X-chromosome, although both the right arm and the center of the second chromosome also contributed to it. Sexual isolation from D. persimilis females was primarily associated with the second chromosome. The absence of isolating mechanisms being associated with many autosomal regions, including some large inverted regions that separate the strains, suggests that these phenotypes may not be caused by genes spread throughout the genome. We suggest that gene flow between these species via hybrid males may be possible at loci spread across much of the autosomes.     

Male‐biased gene flow across an avian hybrid zone: evidence from mitochondrial and microsatellite DNA

Mating pattern and gene flow were studied in the contact zone between two morphologically very similar Chiffchaff taxa (Phylloscopus collybita, P. brehmii) in SW France and northern Spain. Mating was assortative in brehmii, but not in collybita. Mixed matings were strongly asymmetric (excess of collybita male × brehmii female pairs), but did produce viable offspring in some cases. Sequence divergence of the mitochondrial cytochrome b gene was 4.6%. Haplotypes segregated significantly with phenotype (only five ‘mismatches’ among 94 individuals), demonstrating that mitochondrial gene flow was very restricted. The estimated proportion of F₁ hybrids in the reproductive population was significantly lower than expected under a closed population model, indicating strong selection against hybrids. Genetic typing of 101 individuals at four microsatellite loci also showed significant population differentiation, but nuclear gene flow was estimated to be 75 times higher than mitochondrial gene flow. This strong discrepancy is probably due to unisexual hybrid sterility (Haldane’s rule). Thus, there is a strong, but incomplete, reproductive barrier between these taxa.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

Spatio-temporal dynamics of the Allonemobius fasciatus- A. socius mosaic hybrid zone: a 14-year perspective

Long-term studies of hybrid zones can provide valuable insight into a number of questions that have long attracted the attention of evolutionists. These questions range from the stability and fate of hybrid zones to the relative fitness of hybrids. In this paper we report the results of a 14-year survey of the Allonemobius fasciatus-Allonemobius socius hybrid zone. Populations were collected intensively in 1986 and 1987 and then more sporadically through the end of the 1980s and throughout the 1990s. By documenting changes in the genetic composition of populations near and within the zone during this period of time we assessed: the strength of the reproductive isolation between the two species; the relative growth rates (which can be considered a surrogate of relative fitness) of genotype classes corresponding to hybrids and to pure species individuals; and, the power of single-year and multi-year measurements of relative growth rates to predict changes in the genetic composition of mixed populations through time. In brief, we found very large year-to-year variation in the relative growth rates of pure species and hybrid individuals. This variation may reflect the fact that both species are at the edge of their range and perhaps at the limits of their ability to deal with environmental perturbations. As a consequence of the variation, even multi-year estimates of relative growth rates often provided imprecise predictions regarding the future genotypic composition of mixed populations. Despite our limited ability to predict the dynamics of individual populations, some trends are apparent. A. socius, the southern species, has clearly increased in frequency along a transect through the Appalachian Mountains, indicating that the zone is moving north in this region. In contrast, the zone appeared to be more stable along the East Coast transect. Within mixed populations, character-index profiles are often bimodal and stable through time, indicating relatively strong reproductive isolation between the two species that is not being reinforced, nor is it breaking down.