Mosaic Hybrid Zone Research Articles

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.

Heterogeneous genome divergence, differential introgression, and the origin and structure of hybrid zones.

Hybrid zones have been promoted as windows on the evolutionary process and as laboratories for studying divergence and speciation. Patterns of divergence between hybridizing species can now be characterized on a genomewide scale, and recent genome scans have focused on the presence of 'islands' of divergence. Patterns of heterogeneous genomic divergence may reflect differential introgression following secondary contact and provide insights into which genome regions contribute to local adaptation, hybrid unfitness and positive assortative mating. However, heterogeneous genome divergence can also arise in the absence of any gene flow, as a result of variation in selection and recombination across the genome. We suggest that to understand hybrid zone origins and dynamics, it is essential to distinguish between genome regions that are divergent between pure parental populations and regions that show restricted introgression where these populations interact in hybrid zones. The latter, more so than the former, reveal the likely genetic architecture of reproductive isolation. Mosaic hybrid zones, because of their complex structure and multiple contacts, are particularly good subjects for distinguishing primary intergradation from secondary contact. Comparisons among independent hybrid zones or transects that involve the 'same' species pair can also help to distinguish between divergence with gene flow and secondary contact. However, data from replicate hybrid zones or replicate transects do not reveal consistent patterns; in a few cases, patterns of introgression are similar across independent transects, but for many taxa, there is distinct lack of concordance, presumably due to variation in environmental context and/or variation in the genetics of the interacting populations.

Genotype-environment associations support a mosaic hybrid zone between two tidal marsh birds.

Local environmental features can shape hybrid zone dynamics when hybrids are bounded by ecotones or when patchily distributed habitat types lead to a corresponding mosaic of genotypes. We investigated the role of marsh-level characteristics in shaping a hybrid zone between two recently diverged avian taxa - Saltmarsh (Ammodramus caudacutus) and Nelson's (A.nelsoni) sparrows. These species occupy different niches where allopatric, with caudacutus restricted to coastal marshes and nelsoni found in a broader array of wetland and grassland habitats and co-occur in tidal marshes in sympatry. We determined the influence of habitat types on the distribution of pure and hybrid sparrows and assessed the degree of overlap in the ecological niche of each taxon. To do this, we sampled and genotyped 305 sparrows from 34 marshes across the hybrid zone and from adjacent regions. We used linear regression to test for associations between marsh characteristics and the distribution of pure and admixed sparrows. We found a positive correlation between genotype and environmental variables with a patchy distribution of genotypes and habitats across the hybrid zone. Ecological niche models suggest that the hybrid niche was more similar to that of A.nelsoni and habitat suitability was influenced strongly by distance from coastline. Our results support a mosaic model of hybrid zone maintenance, suggesting a role for local environmental features in shaping the distribution and frequency of pure species and hybrids across space.

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.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M.galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Unexpected ancestry of Populus seedlings from a hybrid zone implies a large role for postzygotic selection in the maintenance of species.

In the context of potential interspecific gene flow, the integrity of species will be maintained by reproductive barriers that reduce genetic exchange, including traits associated with prezygotic isolation or poor performance of hybrids. Hybrid zones can be used to study the importance of different reproductive barriers, particularly when both parental species and hybrids occur in close spatial proximity. We investigated the importance of barriers to gene flow that act early vs. late in the life cycle of European Populus by quantifying the prevalence of homospecific and hybrid matings within a mosaic hybrid zone. We obtained genotypic data for 11 976 loci from progeny and their maternal parents and constructed a Bayesian model to estimate individual admixture proportions and hybrid classes for sampled trees and for the unsampled pollen parent. Matings that included one or two hybrid parents were common, resulting in admixture proportions of progeny that spanned the whole range of potential ancestries between the two parental species. This result contrasts strongly with the distribution of admixture proportions in adult trees, where intermediate hybrids and each of the parental species are separated into three discrete ancestry clusters. The existence of the full range of hybrids in seedlings is consistent with weak reproductive isolation early in the life cycle of Populus. Instead, a considerable amount of selection must take place between the seedling stage and maturity to remove many hybrid seedlings. Our results highlight that high hybridization rates and appreciable hybrid fitness do not necessarily conflict with the maintenance of species integrity.

Selection on hybrids of ecologically divergent ecotypes of a marine snail: the relative importance of exogenous and endogenous barriers

Unravelling the form of selection acting on hybrids of ecotypes undergoing ecological speciation is essential to understand the mechanisms behind the evolution of reproductive isolation in the face of gene flow. Shell phenotype is known to be affected by natural selection and is involved in the fitness of the marine snail Littorina saxatilis. Here, we studied the association between shell traits and fitness in hybrids in order to determine the relative role of exogenous and endogenous selection in this hybrid zone of L. saxatilis. We show that directional selection is the predominant mode of selection among hybrids. We also show its heterogeneity, affecting different shell traits, within populations at the level of the microhabitat. Therefore, endogenous selection mechanisms are most probably lacking in this hybrid zone and exogenous barriers (pre- and post-zygotic) are possibly one of the main forces behind the evolution of barriers to gene flow between these ecologically divergent ecotypes. This study shows how this barrier might represent an important type of reproductive isolation within ecological speciation, and this should be taken into account in future studies of speciation in hybrid zones. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, , -. ADDITIONAL KEYWORDS: adaptation - introgression - Littorina - mosaic hybrid zone - natural selection - post-zygotic isolation - premating isolation.

Multilocus analyses indicate a mosaic distribution of hybrid populations in ground squirrels (genusIctidomys)

DNA sequence data from mitochondrial cytochrome-b (Cytb) and Y-linked structural maintenance of chromosomes (SmcY) genes were combined with 478 nuclear loci obtained from amplified fragment length polymorphisms (AFLP) to assess the extent of hybridization and genetic spatial structure of populations in two hybridizing species of ground squirrel (Ictidomys parvidens and Ictidomys tridecemlineatus). Based on AFLP analyses of 134 individuals from 28 populations, 10 populations were identified that possessed hybrid individuals. Overall estimates of FST values revealed strong support for population structure in the Cytb data set; however, analyses of the SmcY gene and the AFLP data indicated ongoing gene flow between species. Pairwise FST comparisons of populations were not significant for the SmcY gene; although they were significant for the Cytb gene, indicating that these populations were structured and that gene flow was minimal. Therefore, gene flow between I. parvidens and I. tridecemlineatus appeared to be restricted to populations that exhibited hybridization. In addition, the fragmented nature of the geographic landscape suggested limited gene flow between populations. As a result, the distributional pattern of interspersed parental and hybrid populations were compatible with a mosaic hybrid zone model. Because ground squirrels display female philopatry and male-biased dispersal, the ecology of these species is compatible with this hypothesis.

Geographic variation in allele frequency of the gamete recognition protein M7 lysin throughout a mosaic blue mussel hybrid zone

Divergence at gamete recognition loci is hypothesized to result in speciation in broadcast spawning invertebrates. Many gamete recognition loci evolve quickly and show patterns of positive selection, yet the advantage of divergence is rarely known. M7 lysin is a sperm protein in the Mytilus edulis species complex that shows evidence of adaptive evolution. This locus is polymorphic with two distinct clades within Mytilus galloprovincialis, one of which, the D clade, shows the strongest signal of positive selection. We tested whether the geographic patterns in allele frequency were consistent with the hypothesis that positive selection on D clade alleles (G D ) was due to reinforcement. Populations of M. edulis showed little evidence of introgression of G D alleles, but there was no consistent evidence of reproductive character displacement. The lack of consistent patterns expected of reinforcement suggests that another mechanism is likely responsible for the adaptive divergence of M7 lysin.

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 F1 hybrids are viable and fertile in the lab, this suggests that strong pre-zygotic barriers are operating in this portion of the hybrid zone.

Characterization of a contemporaneous hybrid zone between two darter species (Etheostoma bison and E. caeruleum) in the Buffalo River System

Hybrid zones have long intrigued evolutionary biologists and provide a natural laboratory to explore the evolution of reproductive isolation (speciation). Molecular characterization of hybrid zone dynamics can provide insight into the strength of reproductive isolation as well as the underlying evolutionary processes shaping gene flow. Approximately one-third of darter species naturally hybridize making this species-rich North American freshwater teleost fish clade an ideal system to investigate the extent and direction of hybridization. The objective of this study was to use diagnostic microsatellite markers to calculate genetic hybrid index scores of two syntopic, but distantly related darter species, Etheostoma bison and Etheostoma caeruleum. A combination of hybrid index scores, assignment tests, and mitochondrial haplotype profiles uncovered mixed ancestry in approximately 6% of sampled adult individuals, supporting contemporaneous hybridization that was previously undocumented in E. bison. Moreover, hybrids were not limited to the F1 generation, but encompassed the entire suite of hybrid categories (F1, F2 and backcross hybrids). The low number of hybrids assigned to each hybrid category represents a bimodal hybrid zone, suggesting reproductive isolation is strong (but incomplete) and also advocates for the ability of hybrids to produce second-generation hybrids and backcross into both parental species, mediating introgression across species boundaries. To this end, cytonuclear profiles of the sampled parental species and hybrids were consistent with bidirectional gene flow, although there was an overall trend of asymmetric hybridization between E. caeruleum females and E. bison males. The spatiotemporal variation in hybridization rates and resulting cytonuclear patterns expanded on in this study provide a comparative genetic framework on which future studies can begin to elucidate the underlying processes that not only generate a mosaic hybrid zone, but maintain the distinctness of species in the face of gene flow.

The geography of introgression in a patchy environment and the thorn in the side of ecological speciation

Abstract When incompletely isolated taxa coexist in a patchy environment (e.g. mosaic hybrid zones, host-race complexes), patterns of variation may differ between selected traits/genes and neutral markers. While the genetic structure of selected traits/loci tends to coincide with habitat variables (producing Genetic-Environment Association or GEA), genetic differentiation at neutral loci unlinked to any selected locus rather depends on geographic connectivity at a large scale (e.g. Isolation-By-Distance or IBD), although these loci often display GEA at a small scale. This discrepancy has been repeatedly taken as evidence for parallel primary divergence driven by local adaptation. We argue that this interpretation needs to be addressed more thoroughly by considering the alternative hypothesis that speciation was initiated in allopatry and secondary introgression has subsequently erased the signal of past differentiation at neutral loci. We present a model of neutral introgression after secondary contact in a mosaic hybrid zone, which describes how GEAs dissipate with time and how neutral variation self-organizes according to the environmental and geographic structures. We show that although neutral loci can be affected by environmental selection, they are often more affected by history and connectivity: the neutral structure retains the initial geographic separation more than it correlates with the environment during the colonization and introgression phases, and then converges to a migration-drift balance, the most frequent outcome of which is GEA at a local scale but IBD at a large scale. This is the exact pattern usually attributed to parallel ecological speciation. Introgression is heterogeneous in space and depends on the landscape structure (e.g. it is faster in small patches, which are more impacted by immigration). Furthermore, there is no directionality in the association and it is possible to observe reversed GEAs between distant regions. We argue that the history of differentiation should ideally be reconstructed with selected loci or neutral loci linked to them, not neutral ones, and review some case studies for which the hypothesis of a long co-existence of co-adapted genetic backgrounds might have been refuted too hastily.

Habitat Association and Seasonality in a Mosaic and Bimodal Hybrid Zone between Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae)

Understanding why some hybrid zones are bimodal and others unimodal can aid in identifying barriers to gene exchange following secondary contact. The hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi contains a mix of allopatric parental populations and inter-mingled bimodal and unimodal sympatric populations, and provides an ideal system to examine the roles of local selection and gene flow between populations in maintaining bimodality. However, it is first necessary to confirm, over a larger spatial scale, previously identified associations between population composition and season and habitat. Here we use cline-fitting of one morphological and one song trait along two valley transects, and intervening mountains, to confirm previously identified habitat associations (mountain versus valley) and seasonal changes in population composition. As expected from previous findings of studies on a smaller spatial scale, C. jacobsi dominated mountain habitats and mixed populations dominated valleys, and C. brunneus became more prevalent in August. Controlling for habitat and incorporating into the analysis seasonal changes in cline parameters and the standard errors of parental trait values revealed wider clines than previous studies (best estimates of 6.4 to 24.5 km in our study versus 2.8 to 4.7 km in previous studies) and increased percentage of trait variance explained (52.7% and 61.5% for transects 1 and 2 respectively, versus 17.6%). Revealing such strong and consistent patterns within a complex hybrid zone will allow more focused examination of the causes of variation in bimodality in mixed populations, in particular the roles of local selection versus habitat heterogeneity and gene flow between differentiated populations.

The coupling hypothesis: why genome scans may fail to map local adaptation genes

Genomic scans of multiple populations often reveal marker loci with greatly increased differentiation between populations. Often this differentiation coincides in space with contrasts in ecological factors, forming a genetic-environment association (GEA). GEAs imply a role for local adaptation, and so it is tempting to conclude that the strongly differentiated markers are themselves under ecologically based divergent selection, or are closely linked to loci under such selection. Here, we highlight an alternative and neglected explanation: intrinsic (i.e. environment-independent) pre- or post-zygotic genetic incompatibilities rather than local adaptation can be responsible for increased differentiation. Intrinsic genetic incompatibilities create endogenous barriers to gene flow, also known as tension zones, whose location can shift over time. However, tension zones have a tendency to become trapped by, and therefore to coincide with, exogenous barriers due to ecological selection. This coupling of endogenous and exogenous barriers can occur easily in spatially subdivided populations, even if the loci involved are unlinked. The result is that local adaptation explains where genetic breaks are positioned, but not necessarily their existence, which can be best explained by endogenous incompatibilities. More precisely, we show that (i) the coupling of endogenous and exogenous barriers can easily occur even when ecological selection is weak; (ii) when environmental heterogeneity is fine-grained, GEAs can emerge at incompatibility loci, but only locally, in places where habitats and gene pools are sufficiently intermingled to maintain linkage disequilibria between genetic incompatibilities, local-adaptation genes and neutral loci. Furthermore, the association between the locally adapted and intrinsically incompatible alleles (i.e. the sign of linkage disequilibrium between endogenous and exogenous loci) is arbitrary and can form in either direction. Reviewing results from the literature, we find that many predictions of our model are supported, including endogenous genetic barriers that coincide with environmental boundaries, local GEA in mosaic hybrid zones, and inverted or modified GEAs at distant locations. We argue that endogenous genetic barriers are often more likely than local adaptation to explain the majority of Fst-outlying loci observed in genome scan approaches - even when these are correlated to environmental variables.

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.

Morphological, Genetic, and Chromosomal Variation at a Small Spatial Scale within a Mosaic Hybrid Zone of the Grasshopper Dichroplus pratensis Bruner (Acrididae)

Hybrid zones are regions where genetically different populations meet and mate, resulting in offspring of mixed characteristics. In organisms with limited dispersal, such as melanopline grasshoppers, hybrid zones can occur at small spatial scales (i.e., <500 m). We assessed levels of morphological, chromosomal, and molecular variability in adult males of the grasshopper Dichroplus pratensis Bruner (N = 137 males, 188 females) collected at 12 sites within a mosaic hybrid zone in a heterogeneous environment in Sierra de la Ventana, Argentina. In this hybrid zone, 2 Robertsonian chromosomal races, polymorphic for different centric fusions, meet (the "Northern race" at low altitudes and the "Southern race" at higher altitudes), forming hybrids that show monobrachial homologies during meiosis. High morphometric variation in 6 traits was revealed among grasshoppers of both sexes, with male body size positively and significantly correlated with increasing altitude. Frequency of Robertsonian fusions characteristic of the Southern race increased significantly with altitude. Moreover, fusion frequencies covaried between samples. Considerable genetic variation was revealed by random amplification of polymorphic DNA markers, with heterozygosity ranging from 0.3477 to 0.3745. Insects from low-altitude and high-altitude populations showed significant genetic differentiation, as indicated by F(ST) values. The proposed model for D. pratensis, involving the generation and maintenance by chromosomal fusions, of gene complexes adaptive in different environments, could explain the observed clinal patterns within the contact zone.

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.

The Flow of Antimicrobial Peptide Genes Through a Genetic Barrier Between Mytilus edulis and M. galloprovincialis

We studied the population genetics of two antimicrobial peptide (AMP) loci, called Mytilin B and Mytilus galloprovincialis defensin 2 (MGD2), in the secondary contact mosaic hybrid zone between Mytilus edulis and M. galloprovincialis. The isolation period between the two species was estimated to be approximately 1 million years (range, 0.5 million to 2 million years) long. During this period, coevolution between microbes and the immune system has likely occurred. The secondary contact, which would date back to approximately 25,000 (0-200,000) years, recently allowed these coadaptations to be rearranged through hybridization. Distinctive polymorphisms were uncovered in coding sequences of the two AMP loci such as insertion/deletion of codons or bisubstituted codons. Very low levels of differentiation were observed between populations of the two species at both loci, while other nuclear loci often showed marked structure among the same samples. The absence of population differentiation proved to be the consequence of secondary introgression of highly divergent alleles. While only a few recombinants were observed at the Mytilin B locus, the MGD2 locus showed a high intragenic recombination rate, which increased in the exon coding for the mature peptide. In addition, standard neutrality tests revealed significant deviations from the mutation-drift equilibrium at both loci. These results suggest that either balancing or directional selection is likely to play a role in the evolution of the two AMPs and introgression would be adaptive. However, evidence accumulated at the Mytilin B locus allows neither for identification of the direction of selection nor for any conclusions on whether selection acted directly on the antimicrobial peptide itself. At the MGD2 locus, a spatial variation of polymorphism patterns along the sequence suggests that selection was direct, although the precise nature of the selection (directional vs. balancing) remains unclear. This study concurs with previous reports of an effect of slight selection on AMP genes evolution in other invertebrates, although selection does not necessarily act on the mature peptides.

Variable extent of hybridization between invasive sika (Cervus nippon) and native red deer (C. elaphus) in a small geographical area

Instances of hybridization between endemic and alien species pose a threat to species integrity but also provide us with an opportunity to study the dynamics of gene flow between two species as they first meet. Here, we used variation at 22 highly differentiated microsatellite loci and one mitochondrial DNA (mtDNA) marker in a sample of 735 individuals, to investigate the genetic consequences of an introduction of Japanese sika deer (Cervus nippon) for native red deer (C. elaphus) on the Kintyre Peninsula in Scotland. We investigated population structure, estimated null-allele frequency and assigned individual hybrid scores using a Bayesian clustering algorithm implemented in structure 2.2. The dataset clearly divided into two clusters and generally, introgression into red and sika was low. However at one site, West Loch Awe, 43% of individuals were hybrids. MtDNA introgression indicated that hybridization was occurring between red-deer hinds and sika-deer stags. We argue that the pattern of differential introgression across the study area is primarily due to the rarity of hybridization events between the two species and the limited time the two species have been in contact (< 120 years). This contrasts with the causes of classic mosaic hybrid zones (selection induced by habitat variability). Currently, it seems possible that, in time, the level of hybridization found at West Loch Awe could also be found across the whole of the peninsula.