Local adaptation and species segregation in two mussel (Mytilus edulis X M. trossulus) hybrid zones

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and M. galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; M. edulis and M. trossulus are sympatric and hybridize on both western and eastern coasts of the Atlantic Ocean. We review the dynamics of hybridization in these two hybrid zones and evaluate the role of local adaptation for maintaining species boundaries. In Scandinavia, hybridization and gene introgression is so extensive that no individuals with pure M. trossulus genotypes have been found. However, M. trossulus alleles are maintained at high frequencies in the extremely low salinity Baltic Sea for some allozyme genes. A synthesis of reciprocal transplantation experiments between different salinity regimes shows that unlinked Gpi and Pgm alleles change frequency following transplantation, such that post-transplantation allelic composition resembles native populations found in the same salinity. These experiments provide strong evidence for salinity adaptation at Gpi and Pgm (or genes linked to them). In the Canadian Maritimes, pure M. edulis and M. trossulus individuals are abundant, and limited data suggest that M. edulis predominates in low salinity and sheltered conditions, whereas M. trossulus are more abundant on the wave-exposed open coasts. We suggest that these conflicting patterns of species segregation are, in part, caused by local adaptation of Scandinavian M. trossulus to the extremely low salinity Baltic Sea environment.

Using a Bayesian method to assign individuals to karyotypic taxa in shrew hybrid zones.

Individuals sampled in hybrid zones are usually analysed according to their sampling locality, morphology, behaviour or karyotype. But the increasing availability of genetic information more and more favours its use for individual sorting purposes and numerous assignment methods based on the genetic composition of individuals have been developed. The shrews of the Sorex araneus group offer good opportunities to test the genetic assignment on individuals identified by their karyotype. Here we explored the potential and efficiency of a Bayesian assignment method combined or not with a reference dataset to study admixture and individual assignment in the difficult context of two hybrid zones between karyotypic species of the Sorex araneus group. As a whole, we assigned more than 80% of the individuals to their respective karyotypic categories (i.e. 'pure' species or hybrids). This assignment level is comparable to what was obtained for the same species away from hybrid zones. Additionally, we showed that the assignment result for several individuals was strongly affected by the inclusion or not of a reference dataset. This highlights the importance of such comparisons when analysing hybrid zones. Finally, differences between the admixture levels detected in both hybrid zones support the hypothesis of an impact of chromosomal rearrangements on gene flow.

Chromosomal rearrangements do not seem to affect the gene flow in hybrid zones between karyotypic races of the common shrew (Sorex araneus).

Chromosomal rearrangements are proposed to promote genetic differentiation between chromosomally differentiated taxa and therefore promote speciation. Due to their remarkable karyotypic polymorphism, the shrews of the Sorex araneus group were used to investigate the impact of chromosomal rearrangements on gene flow. Five intraspecific chromosomal hybrid zones characterized by different levels of karyotypic complexity were studied using 16 microsatellites markers. We observed low levels of genetic differentiation even in the hybrid zones with the highest karyotypic complexity. No evidence of restricted gene flow between differently rearranged chromosomes was observed. Contrary to what was observed at the interspecific level, the effect of chromosomal rearrangements on gene flow was undetectable within the S. araneus species.

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones.

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.

Alps, genes, and chromosomes: their role in the formation of species in the Sorex araneus group (Mammalia, Insectivora), as inferred from two hybrid zones.

During the Pleistocene glaciations, the Alps were an efficient barrier to gene flow between isolated populations, often leading to allopatric speciation. Afterwards, the Alps strongly influenced the post-glacial recolonization of Europe and represent a major suture zone between differentiated populations. Two hybrid zones in the Swiss and French Alps between genetically and chromosomally well-differentiated species-the Valais shrew, Sorex antinorii, and the common shrew, S. araneus-were studied karyotypically and by analyzing the distribution of seven microsatellite loci. In the center of the Haslital hybrid zone the two species coexist over a distance of 900 m. Hybrid karyotypes, among them the most complex known in Sorex, are rare. F-statistics based on microsatellite data revealed a strong heterozygote deficit only in the center of the zone, due to the sympatric distribution of the two species with little hybridization between them. Structuring within the species (both F(IS) and F(ST)) was low. An hierarchical analysis showed a high level of interspecific differentiation. Results were compared with those previously reported in another hybrid zone located at Les Houches in the French Alps. Genetic structuring within and between species was comparable in both hybrid zones, although chromosomal incompatibilities are more important in Haslital, where a linkage block of the race-specific chromosomes should additionally impede gene flow. Evidence for a more restricted gene flow in Haslital comes from the genetically intermediate hybrid karyotypes, whereas in Les Houches, hybrid karyotypes are genetically identical to individuals of the pure karyotypic races. Genic and chromosomal introgression was observed in Les Houches, but not in Haslital. The possible influence of a river, separating the two species at Les Houches, on gene flow is discussed.

Estimating gene flow across hybrid zones: how reliable are microsatellites?

Due to their high polymorphism, microsatellites have become one of the most valued genetic markers in population biology. We review the first two published studies on hybrid zones of the common shrew based on microsatellites. Both reveal surprisingly high interracial gene flow. It can be shown that these are overestimates. Indeed, in classical population genetics models as F-statistics, mutation is neglected. This constitutes an acceptable assumption as long as migration is higher than mutation. However, in hybrid zones where genetic exchanges can be rare, neglecting mutation will lead to strong overestimates of migration when working with microsatellites which display mutation rates up to 10(-3). As there seems to be no straightforward way to correct for this bias, interracial gene flow estimates based on microsatellites should be taken with caution. This problem should however not conceal the enormous potential of microsatellites to unravel the genetics of hybrid zones.

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.

The role of Wolbachia bacteria in reproductive incompatibilities and hybrid zones of Diabrotica beetles and Gryllus crickets

A rickettsial bacterium in the genus Wolbachia is the cause of a unidirectional reproductive incompatibility observed between two major beetle pests of maize, the western corn rootworm, Diabrotica virgifera virgifera, and the Mexican corn rootworm, D. v. zeae. These subspecies are allopatric except for two known regions of sympatry in Texas and Mexico. We demonstrate that populations of D. v. virgifera, with the exception of two populations in southern Arizona, are infected with a strain of Wolbachia. Populations of D. v. zeae are not infected. Treatment of D. v. virgifera with tetracycline eliminated the Wolbachia and removed the reproductive incompatibility. Similar patterns of reproductive incompatibility exist among taxa of the cricket genus Gryllus. Gryllus assimilis, G. integer, G. ovisopis, G. pennsylvanicus, and G. rubens are infected with Wolbachia whereas G. firmus is usually not. Populations of G. rubens and G. ovisopis carry the same Wolbachia strain, which is distinct from that of G. integer. G. pennsylvanicus is infected with two Wolbachia strains, that found in G. rubens and one unique to G. pennsylvanicus. Moreover, a proportion of G. pennsylvanicus individuals harbors both strains. Wolbachia may have influenced speciation in some members of the genus Gryllus by affecting the degree of hybridization between species. Given that Wolbachia infections are relatively common in insects, it is likely that other insect hybrid zones may be influenced by infections with Wolbachia.

Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

Study of gene flow through a hybrid zone in the common shrew (Sorex araneus) using microsatellites

The common shrew (Sorer araneus) is subdivided into several chromosomal races. As hybrid zones between them have been characterized, this organism is of particular interest in studying the role of chromosomes in speciation. Six microsatellite loci were used to evaluate the level of gene how in the S. araneus hybrid zone between the Cordon and Valais races. Most of these loci were very polymorphic, the total number of alleles detected per locus ranging from 3 to 20. Using Mantel tests, we showed that the effect of rivers as barriers to gene flow is less important at this sampling scale. The effect of the chromosomal race is of particular importantance in diminishing gene flow.

Chromosomal rearrangements and gene flow over time in an inter-specific hybrid zone of the Sorex araneus group.

Most hybrid zones have existed for hundreds or thousands of years but have generally been observed for only a short time period. Studies extending over periods long enough to track evolutionary changes in the zones or assess the ultimate outcome of hybridization are scarce. Here, we describe the evolution over time of the level of genetic isolation between two karyotypically different species of shrews (Sorex araneus and Sorex antinorii) at a hybrid zone located in the Swiss Alps. We first evaluated hybrid zone movement by contrasting patterns of gene flow and changes in cline parameters (centre and width) using 24 microsatellite loci, between two periods separated by 10 years apart. Additionally, we tested the role of chromosomal rearrangements on gene flow by analysing microsatellite loci located on both rearranged and common chromosomes to both species. We did not detect any movement of the hybrid zone during the period analysed, suggesting that the zone is a typical tension zone. However, the gene flow was significantly lower among the rearranged than the common chromosomes for the second period, whereas the difference was only marginally significant for the first period. This further supports the role of chromosomal rearrangements on gene flow between these taxa.

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup).

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3-1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5-1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.

A hybrid zone with coincident clines for autosomal and sex-specific markers in the Sorex araneus group.

In hybrid zones, endogenous counter-selection of hybrids is usually first expressed as reduced fertility or viability in hybrids of the heterogametic sex, a mechanism known as Haldane's rule. This phenomenon often leads to a differential of gene flow between sex-linked markers. Here, we address the possibility of a differential gene flow for Y chromosome, mtDNA and autosomal markers across the hybrid zone between the genetically and chromosomally well-differentiated species Sorex antinorii and Sorex araneus race Vaud. Intermarker comparison clearly revealed coincidental centre and very abrupt clines for all three types of markers. The overall level of genetic differentiation between the two species must be strong enough to hinder asymmetric introgression. Cyto-nuclear mismatches were also observed in the centre of hybrid zone. The significantly lower number of mismatches observed in males than in females possibly results from Y chromosome-mtDNA interactions. Results are compared with those previously reported in another hybrid zone between S. antinori and S. araneus race Cordon.

Isotopic variation across the Audubon's-myrtle warbler hybrid zone.

Differences in seasonal migratory behaviours are thought to be an important component of reproductive isolation in many organisms. Stable isotopes have been used with success in estimating the location and qualities of disjunct breeding and wintering areas. However, few studies have used isotopic data to estimate the movements of hybrid offspring in species that form hybrid zones. Here, we use stable hydrogen to estimate the wintering locations and migratory patterns of two common and widespread migratory birds, Audubon's (Setophaga auduboni) and myrtle (S. coronata) warblers, as well as their hybrids. These two species form a narrow hybrid zone with extensive interbreeding in the Rocky Mountains of British Columbia and Alberta, Canada, which has been studied for over four decades. Isotopes in feathers grown on the wintering grounds or early on migration reveal three important patterns: (1) Audubon's and myrtle warblers from allopatric breeding populations winter in isotopically different environments, consistent with band recovery data and suggesting that there is a narrow migratory transition between the two species, (2) most hybrids appear to overwinter in the south-eastern USA, similar to where myrtle warblers are known to winter, and (3) some hybrid individuals, particularly those along the western edge of the hybrid zone, show Audubon's-like isotopic patterns. These data suggest there is a migratory divide between these two species, but that it is not directly coincident with the centre of the hybrid zone in the breeding range. We interpret these findings and discuss them within the context of previous research on hybrid zones, speciation and migratory divides.