Historical and size-dependent genetic variation in hybrid mussel populations

Integration of natural selection across the life cycle stabilizes a marine mussel hybrid zone

Hybridization among related species is now recognized as common but it remains unclear how hybrid zones persist for prolonged periods. Here, we test the hypothesis that selection in different components of the life cycle may stabilize a hybrid zone. A hybrid zone occurs in southwest England between the marine mussels Mytilus edulis and M. galloprovincialis. Previous studies have found strong directional selection against alleles from M. edulis occurs among hybrids in the adult stage. Traditional hybrid zone models argue that alleles that are selected within the hybrid zone are replaced by migration from neighboring parental population into the hybrid zone. In this system, however, migration occurs out of this hybrid zone into neighboring parental populations. This hybrid zone should therefore be unstable and dissipate, yet this zone has persisted for more than 30 years. We tested and rejected the hypothesis that differences in fecundity may select for M. edulis alleles within this hybrid zone and thus counter the selection observed against these alleles among adults. We also tested the hypothesis that selection during the larval stage may counter selection against M. edulis alleles in the adult stage. We found that selection favors M. edulis alleles during the veliger stage of larval development. The direction and strength of selection during the larval stage are sufficient to counter strong selection during the adult portion of the life cycle. This hybrid zone is stabilized by opposing forms of directional selection operating in different portions of the life cycle.

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

Genome-wide insights into introgression and its consequences for genome-wide heterozygosity in the Mytilus species complex across Europe

The three mussel species comprising the Mytilus complex are widespread across Europe and readily hybridize when they occur in sympatry, resulting in a mosaic of populations with varying genomic backgrounds. Two of these species, M. edulis and M. galloprovincialis, are extensively cultivated across Europe, with annual production exceeding 230,000 tonnes. The third species, M. trossulus, is considered commercially damaging as hybridization with this species results in weaker shells and poor meat quality. We therefore used restriction site associated DNA sequencing to generate high-resolution insights into the structure of the Mytilus complex across Europe and to resolve patterns of introgression. Inferred species distributions were concordant with the results of previous studies based on smaller numbers of genetic markers, with M. edulis and M. galloprovincialis predominating in northern and southern Europe respectively, while introgression between these species was most pronounced in northern France and the Shetland Islands. We also detected traces of M. trossulus ancestry in several northern European populations, especially around the Baltic and in northern Scotland. Finally, genome-wide heterozygosity, whether quantified at the population or individual level, was lowest in M. edulis, intermediate in M. galloprovincialis, and highest in M. trossulus, while introgression was positively associated with heterozygosity in M. edulis but negatively associated with heterozygosity in M. galloprovincialis. Our study will help to inform mussel aquaculture by providing baseline information on the genomic backgrounds of different Mytilus populations across Europe and by elucidating the effects of introgression on genome-wide heterozygosity, which is known to influence commercially important traits such as growth, viability, and fecundity in mussels.

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.

Empirical study of hybrid zone movement

Hybrid zones are 'natural laboratories' for studying the origin, maintenance and demise of species. Theory predicts that hybrid zones can move in space and time, with significant consequences for both evolutionary and conservation biology, though such movement is often perceived as rare. Here, a review of empirical studies of moving hybrid zones in animals and plants shows 23 examples with observational evidence for movement, and a further 16 where patterns of introgression in molecular markers could be interpreted as signatures of movement. The strengths and weaknesses of methods used for detecting hybrid zone movement are discussed, including long-term replicated sampling, historical surveys, museum/herbarium collections, patterns of relictual populations and introgression of genetic markers into an advancing taxon. Factors governing hybrid zone movement are assessed in the light of the empirical studies, including environmental selection, competition, asymmetric hybridization, dominance drive, hybrid fitness, human activity and climate change. Hybrid zone movement means that untested assumptions of stability in evolutionary studies on hybrid zone can lead to mistaken conclusions. Movement also means that conservation effort aimed at protecting against introgression could unwittingly favour an invading taxon. Moving hybrid zones are of wide interest as examples of evolution in action and possible indicators of environmental change. More long-term experimental studies are needed that incorporate reciprocal transplants, hybridization experiments and surveys of molecular markers and population densities on a range of scales.

Evolution of 2-DE protein patterns in a mussel hybrid zone

Protein expression patterns were compared in mussels from a hybrid zone between the two species Mytilus edulis and M. galloprovincialis using 2-DE. Significant differences in expression pattern were observed between species and between mussels within species. Hybrid mussels had more variable protein expression patterns than mussels of each species. This could be due to segregation at expression modifier loci in the hybrids. It is proposed that unusual hybrid expression patterns might contribute to postzygotic isolation, if such patterns reduce fitness. The use of proteomic data for testing evolutionary models, as in many transcriptomics studies, is explored with results consistent with the expectations of neutral theory.

Disruption of doubly uniparental inheritance of mitochondrial DNA in hybrid mussels (Mytilus edulis x M. galloprovincialis)

Blue mussels of the genus Mytilus have an unusual mode of mitochondrial DNA inheritance termed doubly uniparental inheritance (DUI). Females are homoplasmic for the F mitotype which is inherited maternally, whereas males are heteroplasmic for this and the paternally inherited M mitotype. In areas where species distributions overlap a varying degree of hybridization occurs; yet genetic differences between allopatric populations are maintained. Observations from natural populations and previous laboratory experiments suggest that DUI may be disrupted by hybridization, giving rise to heteroplasmic females and homoplasmic males. We carried out controlled laboratory crosses between Mytilus edulis and M. galloprovincialis to produce pure species and hybrid larvae of known parentage. DNA markers were used to follow the fate of the F and M mitotypes through larval development. Disruption of the mechanism which determines whether the M mitotype is retained or eliminated occurred in an estimated 38% of M. edulis x M. galloprovincialis hybrid larvae, a level double that previously observed in adult mussels from a natural M. edulis x M. galloprovincialis hybrid population. Furthermore, reciprocal hybrid crosses exhibited contrasting types of DUI disruption. The results indicate that disruption of DUI in hybrid mussels may be associated with increased mortality and hence could be a factor in the maintenance of genetic integrity for each species.

Habitat preference and the marine-speciation paradox

Marine organisms challenge the classical theories of local adaptation and speciation because their planktonic larvae have the potential to maintain high gene flow. The marine-speciation paradox is illustrated by contact zones between incipient species that are so large that allopatric divergence seems unlikely. For this reason any mechanism preventing sympatric larvae of two incipient species from coexisting in the same habitats can be a powerful promoter of speciation. The contact zone between two hybridizing taxa of mussel, Mytilus edulis and M. galloprovincialis, in Europe provides an excellent example. Although the zone itself extends over thousands of kilometres, the opportunities for interbreeding are considerably reduced by the small-scale mosaic structure of the zone, where local patches of each taxon alternate at scales of kilometres or less, in response to locally variable ecological factors. Habitat choice by settling larvae would be a less costly mechanism than post-settlement selection to maintain such a mosaic structure. Unfortunately the role of selective settlement has remained hypothetical because larvae could not be scored by classical genetic markers. PCR markers allowed us to study larvae and settlement in ecologically contrasting sites within the zone. We show that only a subset of the genotypes present in the plankton settle in some sites, and that the adults on these sites show the same genetic bias. Genetically based variation in pre-settlement processes therefore accounts for the ecological segregation observed, though it is not the only factor involved in limiting successful interbreeding. The present dataset also supports previous reports of partial spawning asynchrony.

Patterns of larval dispersal and their effect on the maintenance of a blue mussel hybrid zone in southwestern England

The blue mussels Mytilus edulis and M. galloprovincialis hybridize in southwestern England. Within this hybrid zone environmentally based directional selection favors individuals with alleles specific to M. galloprovincialis. What forces are countering this directional selection and allowing for the maintenance of a stable hybrid population are unknown. We used both the genetics of recently settled larvae and a fine-scale model of the physical oceanography of the region to determine the patterns of larval dispersal throughout the hybrid zone and the bordering parental populations. Evidence from both the model and the genetics suggests that the hybrid zone lies between two barriers to dispersal. Start Point separates the M. edulis population from the hybrid zone and allows minimal dispersal from the hybrid zone into the M. edulis population, but none in the other direction. Likewise, the M. galloprovincialis populations along the northern coast of Cornwall regularly receive immigrating larvae from the hybrid zone, but larvae from the M. galloprovincialis population do not enter the hybrid zone. However, larvae settling at hybrid zone sites have high frequencies of alleles specific to M. edulis, suggesting that reproductive barriers, selection in the larval stage, or gene flow from an undetermined source is effectively balancing the directional selection observed in the adults.

Introgression patterns in the mosaic hybrid zone between Mytilus edulis and M. galloprovincialis

Hybrid zones are fascinating systems to investigate the structure of genetic barriers. Marine hybrid zones deserve more investigation because of the generally high dispersion potential of planktonic larvae which allows migration on scales unrivalled by terrestrial species. Here we analyse the genetic structure of the mosaic hybrid zone between the marine mussels Mytilus edulis and M. galloprovincialis, using three length-polymorphic PCR loci as neutral and diagnostic markers on 32 samples along the Atlantic coast of Europe. Instead of a single genetic gradient from M. galloprovincialis on the Iberian Peninsula to M. edulis populations in the North Sea, three successive transitions were observed in France. From South to North, the frequency of alleles typical of M. galloprovincialis first decreases in the southern Bay of Biscay, remains low in Charente, then increases in South Brittany, remains high in most of Brittany, and finally decreases again in South Normandy. The two enclosed patches observed in the midst of the mosaic hybrid zone in Charente and Brittany, although predominantly M. edulis-like and M. galloprovincialis-like, respectively, are genetically original in two respects. First, considering only the various alleles typical of one species, the patches show differentiated frequencies compared to the reference external populations. Second, each patch is partly introgressed by alleles of the other species. When introgression is taken into account, linkage disequilibria appear close to their maximum possible values, indicating a strong genetic barrier within all transition zones. Some pre- or postzygotic isolation mechanisms (habitat specialization, spawning asynchrony, assortative fertilization and hybrid depression) have been documented in previous studies, although their relative importance remains to be evaluated. We also provided evidence for a recent migratory 'short-cut' connecting M. edulis-like populations of the Charente patch to an external M. edulis population in Normandy and thought to reflect artificial transfer of spat for aquaculture.

Assortative fertilization and selection at larval stage in the mussels Mytilus edulis and M. galloprovincialis

Assortative mating (prezygotic isolation) and reduced hybrid fitness (postzygotic isolation) are typically invoked to explain the stability of hybrid zones. In the tension zone model, these factors work in opposition to migration, which promotes genetic homogeneity. Many marine animals migrate over long distances through a planktonic larval stage. Therefore, strong reproductive isolation is needed to maintain stable marine hybrid zones. However, surprisingly little is known about mating preferences and hybrid fitness in marine organisms. Smooth-shelled mussels (Mytilus spp.) form a well-known species complex, with hybridization over extensive areas such as the contact zone of M. edulis and M. galloprovincialis around European Atlantic coasts. This paper reports direct experimental evidence of assortative fertilization, hybrid larval inviability, and early heterosis for growth rate in M. edulis and M. galloprovincialis. Four crosses between pure M. edulis and M. galloprovincialis were analyzed with a new polymerase-chain-reaction-based diagnostic marker. Gamete competition between taxa was allowed in two out of the four crosses. Genotype frequencies observed at an early stage (36 h after fertilization) unambiguously revealed assortative fertilization when gamete competition was allowed. A significant reduction in hybrid viability was subsequently observed during the larval stage. At the same stage an antagonistic effect, heterosis, was observed on growth rate. However, even if heterosis is observed in the F1, it is expected to vanish in subsequent hybrid generations. Although specialization for different habitats and asynchronous spawning have been mentioned as factors contributing to the maintenance of the blue mussel hybrid zone in Europe, we argue that assortative fertilization and reduced hybrid fitness are important factors that also contribute to the stabilization of this zone. These results emphasize that multiple factors may act concomitantly in a barrier to gene flow, especially in complex life cycles. Furthermore, they show that assortative mating through gamete preference, as already demonstrated for sea urchins, may play a role in speciation processes taking place in the sea.

The zone of sympatry and hybridization of Mytilus edulis and M. galloprovincialis, as described by intron length polymorphism at locus mac-1

Intron-size variation at the actin gene locus mac-1 was used to characterize mussel, Mytilus spp., populations in the approximately 2000-km wide zone of contact and hybridization ('hybrid zone') between M. edulis and M. galloprovincialis in western Europe. Twenty-five samples were collected in 1995-99 in locations within the hybrid zone and from reference populations of each species. We used correspondence analysis on the matrix of allelic frequencies to determine which alleles are characteristic of each species, and to characterize samples along the genetic gradient between M. edulis and M. galloprovincialis. In the hybrid zone, some samples exhibited mac-1 allele frequencies that were typical of M. edulis; other samples were distributed along the M. edulis/M. galloprovincialis gradient and displayed variable levels of intergradation that were not correlated with geography. Some of the latter samples exhibited significant heterozygote deficiencies. The simple admixture hypothesis (Wahlund effect) could not be rejected for two-fifths of the samples. The hybrid zone thus appeared as a mosaic of populations which are either pure M. edulis, or hybrid between M. galloprovincialis and M. edulis, or a mixture of the foregoing with M. galloprovincialis individuals. These results were consistent with published allozyme data, suggesting that they can be extended to the entire nuclear genome. M. edulis mac-1 alleles were present at moderate frequency in Atlantic M. galloprovincialis, and at significantly lower frequency in some Mediterranean samples. This pattern was homogeneous over a broad geographical range within each basin. It was not evident that introgression of M. edulis into M. galloprovincialis presently occurs south of the zone of contact. We propose that the distinctness of the Atlantic M. galloprovincialis population results from past introgression by M. edulis alleles.

Evidence for intragenic recombination within a novel genetic marker that distinguishes mussels in the Mytilus edulis species complex

We have used the polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) techniques to design two genetic markers for blue mussels in the Mytilus edulis species complex. Both of these markers target the gene encoding the mussel polyphenolic adhesive protein. The first marker, Glu-5', is highly differentiated among and can be used to identify the three blue mussel species, M. edulis, M. galloprovincialis and M. trossulus. The second marker, Glu-3', can identify M. edulis and M. galloprovincialis. Using these markers we have demonstrated that hybrid mussels from Whitsand Bay, UK carry alleles for this gene that are the products of intragenic recombination. The high frequency (10 per cent) of these recombinant alleles within the hybrid population suggests that recombination is fairly frequent within this gene or that hybridization between M. edulis and M. galloprovincialis is substantial and has been occurring over considerable evolutionary time. The two novel genetic markers, Glu-5' and Glu-3' will be invaluable in additional studies regarding the importance of hybridization among blue mussels.

Differential susceptibility to a trematode parasite among genotypes of the Mytilus edulis/galloprovincialis complex

We show that parasitism by the trematode Prosorhynchus squamatus in parental and introgressed Mytilus edulis/galloprovincialis (Bivalvia) mussels occurs in individuals with a predominantly M. edulis genome. This result suggests that the restricted specificity of P. squamatus is dependent on genetic factor(s) present in M. edulis. Because of its strong pathogenic effects (i.e. total castration and possible death), this parasite may be a source of intense selection against M. edulis genomes when they are present in a site. As a consequence, it may favour the geographic extension of the M. galloprovincialis genome. Previous studies have indicated that, in hybrid zones, recombinant genotypes are more susceptible to parasitic infections than either parental genotype. We demonstrate that this is not the case for the M. edulis/M. galloprovincialis system, and that the parental genotype alone determines susceptibility.