Genetic divergence and a hybrid zone between Baltic and North Sea Mytilus populations (Mytilidae: Mollusca)

Populations of the common mussel ( Mytilus edulis ) from the North Sea area (Skakerrak-Kattegat) and those from the Baltic Sea are almost diagnostically differentiated at five out of 22 studied allozyme loci; at a further seven loci, alleles predominant or common in one area are nearly absent in the other. Genetic distance was estimated at 0.28; this is similar to the distances of these populations to the Mediterranean mussel M. galloprovincialis. The three mussel types obviously represent equal evolutionary divergence from one another, and should also be taxonomically equally separated; a semispecies rank within a more comprehensive M. edulis complex or superspecies is suggested. The age of the Baltic mussel type ( 'M. trossulus' ), as an independent evolutionary lineage, is probably far greater than that of the post-glacial Baltic Sea.
Allele frequencies change gradually and in parallel when entering from the Kattegat through the Sound into the Baltic. Only a slight Wahlund effect at the strongly diverged Gpi and Pgm loci was found in intermediate populations, indicating that extensive hybridization of the two taxa takes place in the area. However, strong interlocus genotypic associations suggest that selection against hybrids is intense in later generations; the c. 100 km wide hybrid zone is narrow relative to the dispersal distance. The genotypic structure of the Lap locus does not conform with those of the other loci studied in the hybrid zone; it cannot be viewed merely as a neutral marker of the process of hybridization.     

Interspecies transfer of female mitochondrial DNA is coupled with role-reversals and departure from neutrality in the mussel Mytilus trossulus.

Mussels of the genus Mytilus have distinct and highly diverged male and female mitochondrial DNA (mtDNA) genomes with separate routes of inheritance. Previous studies of European populations of Mytilus trossulus demonstrated that 33% of males are heteroplasmic for a second mtDNA genome of increased length and that hybridization with Mytilus edulis does not block mtDNA introgression, in contrast to reports for American populations. Here, we demonstrate that the female mtDNA type of M. edulis has replaced the resident female mtDNA type of European M. trossulus. This is supported by COIII sequence data indicating that the female mtDNA of European M. trossulus is very similar to that of M. edulis and that in phylogenetic trees, the mtDNAs of these two species cluster together but separately from American M. trossulus sequences, the latter not being disturbed by introgressive hybridization. We also provide evidence that the mtDNA genome of increased length found in heteroplasmic males of European M. trossulus derives from a recent partition of an introgressed M. edulis female type into the male route of transmission. Neutrality tests reveal that European populations of M. trossulus display an excess of replacement polymorphism within the female mtDNA type with respect to conspecific American populations, as well as a significant excess of rare variants, of a similar magnitude to those previously reported for the invading European M. edulis mtDNA. Results are consistent with a nearly neutral model of molecular evolution and suggest that selection acting on European M. trossulus mtDNA is largely independent of the nuclear genetic background.     

Doubly uniparental inheritance is associated with high polymorphism for rearranged and recombinant control region haplotypes in Baltic Mytilus trossulus

Many bivalve species, including mussels of the genus Mytilus, are unusual in having two mtDNA genomes, one inherited maternally (the F genome) and the other inherited paternally (the M genome). The sequence differences between the genomes are usually great, indicating ancient divergence predating speciation events. However, in Mytilus trossulus from the Baltic, both genomes are similar to the F genome from the closely related M. edulis. This study analyzed the mtDNA control region structure in male and female Baltic M. trossulus mussels. We show that a great diversity of structural rearrangements is present in both sexes. Sperm samples are dominated by recombinant haplotypes with M. edulis M-like control region segments, some having large duplications. By contrast, the rearranged haplotypes that dominate in eggs lack segments from this M genome. The rearrangements can be explained by a combination of tandem duplication, deletion, and intermolecular recombination. An evolutionary pathway leading to the recombinant haplotypes is suggested. The data are also considered in relation to the hypothesis that the M. edulis M-like control region sequence is necessary to confer the paternal role on genomes that are otherwise F-like.     

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.     

Study of two species of mussels, Mytilus trossulus and Mytilus galloprovincialis (Bivalvia, Mytilidae) and their hybrids, using PCR markers, in Peter the Great Bay of the Sea of Japan]

Bivalve mollusks of the genus Mytilus (M. trossulus and M. galloprovincialis) occurring in Peter the Great Bay of the Sea of Japan were first studied in Russia. A region of nonrepetitive sequences of the gene encoding the polyphenolic adhesive protein bissus was used as a species-specific genetic marker. After amplification using specific primers, a 126-bp fragment was found to amplify in all representatives of M. galloprovincialis collected from driftwood in the gulf Posset (the southwestern part of Peter the Great Bay). M. trossulus specimens from the same region were shown to have a 168-bp fragment. In Vostok Gulf (the eastern part of Peter the Great Bay), both artificially grown mussels and those from natural habitats contained a 168-bp fragment or two fragments (126- and 168-bp) that corresponded to a hybrid form between the above species. The possibility of using this genetic marker to identify closely related Mytilus strains and their hybrids in similar habitats, near the Primorye coast in particular, was demonstrated. The presence of approximately 9% of hybrid specimens confirms that a zone of hybridization between M. trossulus and M. galloprovincialis may exist in this region.     

Comparative analysis of different satellite DNAs in four Mytilus species.

We report the characterization of three satellite DNAs in four species of mussel: Mytilus edulis, Mytilus galloprovincialis, Mytilus trossulus, and Mytilus californianus. The monomers of the Apa I satellite DNAs were 173, 161, and 166 bp long. These satellite monomers were used to construct phylogenetic trees to infer relationships among these species. The topologies obtained clearly indicate that M. californianus is the most divergent species with respect to the other three. Furthermore, localization of satellite DNAs on metaphase chromosomes was performed using fluorescent in situ hybridization (FISH). Fluorescent signals revealed a different organization and distribution of these three satellite DNAs.     

Species-Specific Segregation of Gender-Associated Mitochondrial DNA Types in an Area Where Two Mussel Species (Mytilus Edulis and M. Trossulus) Hybridize

In each of the mussel species Mytilus edulis and M. trossulus there exist two types of mtDNA, the F type transmitted through females and the M type transmitted through males. Because the two species produce fertile hybrids in nature, F and M types of one may introgress into the other. We present the results from a survey of a population in which extensive hybridization occurs between these two species. Among specimens classified as ``pure'''' M. edulis or ``pure'''' M. trossulus on the basis of allozyme analysis, we observed no animal that carried the F or the M mitotype of the other species. In most animals of mixed nuclear background, an individual''s mtDNA came from the species that contributed the majority of the individual''s nuclear genes. Most importantly, the two mtDNA types in post-F(1) male hybrids were of the same species origin. We interpret this to mean that there are intrinsic barriers to the exchange of mtDNA between these two species. Because such barriers were not noted in other hybridizing species pairs (many being even less interfertile than M. edulis and M. trossulus), their presence in Mytilus could be another feature of the unusual mtDNA system in this genus.     

Comparison of gamete compatibility between two blue mussel species in sympatry and in allopatry

Recent demonstrations of positive selection on genes controlling gamete compatibility have resulted in a proliferation of hypotheses concerning the sources of selection. We tested a prediction of one prominent hypothesis, selection to avoid hybridization (i.e., reinforcement), by comparing heterospecific gamete compatibility in two Mytilus edulis populations: one population in Cobscook Bay, Maine, in which the close congener, M. trossulus, is abundant (a region of sympatry), and one population in Kittery, Maine, in which M. trossulus is absent (a region of allopatry). Three diagnostic nuclear DNA markers were used to identify mussels to species and to estimate the frequency of both species and their hybrids in the two populations. Controlled crosses were then conducted by combining eggs of M. edulis females with a range of M.edulis and M. trossulus sperm concentrations. Results were not consistent with the reinforcement hypothesis. M. edulis females collected from the region of sympatry were no more incompatible with M. trossulus males than were M. edulis females collected from the region of allopatry. A trend in the opposite direction, toward greater compatibility in sympatry, suggests that introgression of M. trossulus genes that control egg compatibility, such as those encoding receptors for sperm, may influence evolution of gametic isolation in hybridizing populations.     

Complete sequences of mitochondrial genomes from the Baltic mussel Mytilus trossulus

Marine mussels Mytilus possess two mitochondrial (mt) genomes, which undergo doubly uniparental inheritance (DUI). Female (F) and male (M) genomes are usually highly diverged at the sequence level. Both genomes contain the same set of metazoan genes (for 12 proteins, 2 rRNAs and 23 tRNAs), both lack the atp8 gene and have two tRNAs for methionine. However, recently recombination between those variants has been reported. Both original F and M mt genomes of M. trossulus were replaced by M. edulis mtDNA in the Baltic populations. Highly diverged M genome occurs rarely in the Baltic mussels. Full sequences of the M genome identified in males (sperm) and F genome in females (eggs) were obtained. Both genomes were diverged by 24% in nucleotide sequence, but had similar nucleotide composition and codon usage bias. Constant domain (CD) of the control region (CR), the tRNA and rRNA genes were the most conserved. The most diverged was the variable domain 1 (VD1) of the control region. The F genome was longer than M by 147 bp. and the main difference was localised in the VD1 region. No recombination was observed in whole mtDNA of both studied variants. Nuclear mitochondrial pseudogenes (numts) have not been found by hybridisation with probes complementary to several fragments of the Baltic M. trossulus mtDNA.     

The Effects of Natural Hybridization on the Regulation of Doubly Uniparental Mtdna Inheritance in Blue Mussels (Mytilus Spp.)

Blue mussels in the Mytilus edulis species complex have a doubly uniparental mode of mtDNA inheritance with separate maternal and paternal mtDNA lineages. Female mussels inherit their mtDNA solely from their mother, while males inherit mtDNA from both parents. In the male gonad the paternal mtDNA is preferentially replicated so that only paternal mtDNA is transmitted from fathers to sons. Hybridization is common among differentiated blue mussel taxa; whenever it involves M. trossulus, doubly uniparental mtDNA inheritance is disrupted. We have found high frequencies of males without and females with paternal mtDNA among hybrid mussels produced by interspecific matings between M. galloprovincialis and M. trossulus. In contrast, hybridization between M. galloprovincialis and M. edulis does not affect doubly uniparental inheritance, indicating a difference in the divergence of the mechanisms regulating mtDNA inheritance among the three blue mussel taxa. Our data indicate a high frequency of disrupted mtDNA transmission in F(1) hybrids and suggest that two separate mechanisms, one regulating the transmission of paternal mtDNA to males and another inhibiting the establishment of paternal mtDNA in females, act to regulate doubly uniparental inheritance. We propose a model for the regulation of doubly uniparental inheritance that is consistent with these observations.     

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.     

Evidence of selective mortality in favour of the Mytilus galloprovincialis Lmk phenotype in British mussel populations

Samples of mussels ( Mytilus ) were collected from 17 localities within hybrid zones of Mytilus edulis and Mytilus galloprovincialis in south-west and north-east England. The study of two polymorphic allozyme loci ( esterase-D and octopine dehydrogenase ), which are partially diagnostic for the two forms of mussel, reveal the existence of widespread length-dependent allele frequency variation. Larger mussels tend to have a higher frequency of alleles characteristically at high frequency in Mytilus galloprovincialis. Also at a given shell length galloprovincialis alleles have a higher frequency higher up the shore. Computer simulation is used to demonstrate that length-dependent variation may be generated not only by differential mortality but also by differential growth and in models including or excluding immigration. Evidence supports the hypothesis that selective mortality acting in favour of the galloprovincialis phenotype within hybrid populations in Britain is balanced by immigration of the more abundant Mytilus edulis.     

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°‐13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.     

Comparative Genomics of Marine Mussels (Mytilus spp.) Gender Associated mtDNA: Rapidly Evolving atp8

The unusual mode of mitochondrial DNA inheritance, with two separate: maternal (F) and paternal (M) lineages, gives unique opportunities to study the evolution of the mitochondrial genome. This system was first discovered in the marine mussels Mytilus. The three related species: Mytilus edulis, Mytilus galloprovincialis and Mytilus trossulus form a complex in which the divergence of M and F lineages pre-dates the speciation. The complete mitochondrial genomes of both lineages were known for all species except Pacific M. trossulus. Here we report, for the first time, the complete sequences of both mitochondrial genomes of Pacific M. trossulus, filling the gap. While the reported M and F genomes are highly diverged (26%), they have similar organisation. The only difference is the translocation of one tRNA gene into the long, mosaic control region of the F genome. Consistent presence of an ORF which most likely represents the atp8 gene was confirmed in both genomes. The predicted protein has characteristics expected of the functional atp8 even though the M and F versions are markedly different in length. Comparative analysis involving all three species led to the conclusion that the cause of a faster evolution of atp8 and Mytilus mtDNA in general is most likely the Compensation-Draft Feedback process coupled with relatively relaxed selection in the M lineage. Thus, we postulate that the adaptive changes may have played a role in the emergence of highly diverged, barely recognizable atp8 in Mytilus mussels.     

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.     

Seasonal variability of free amino acids in two marine bivalves, Macoma balthica and Mytilus spp., in relation to environmental and physiological factors

The seasonal variability of the intracellular free amino acid (FAA) concentration was studied in 5 Macoma balthica populations and 7 Mytilus spp. populations along their European distribution. Because of the well known physiological role of FAA as organic osmolytes for salinity induced cell volume regulation in marine osmoconformers, FAA variations were compared in bivalve populations that were exposed to high vs. low intraannual salinity fluctuations. In general, seasonal FAA variations were more pronounced in M. balthica than in Mytilus spp. In both bivalve taxa from different locations in the Baltic Sea, highest FAA concentrations were found in autumn and winter and low FAA concentrations were measured in summer. Seasonal patterns were less pronounced in both taxa at locations with constant salinity conditions. In contrast to Baltic Sea populations, Atlantic and Mediterranean bivalves showed high FAA concentrations in summer and low values in winter, regardless of seasonal salinity fluctuations. Significant seasonal FAA variations at locations with constant salinity conditions showed that salinity appeared not to be the main factor in determining FAA concentrations. The seasonal patterns of the main FAA pool components, i.e. alanine, glycine and taurine, are discussed in the context of seasonal variations in environmental factors (salinity, temperature) and physiological state (glycogen content, reproductive stage).     

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.     

Quantitative analysis of sperm plane circular movement in the blue mussels Mytilus edulis, M. trossulus and their hybrids

Fertilization success of free spawning organisms such as Mytilus species depends on gamete interactions. Therefore, gamete traits such as sperm movement are important for determining fertilization success in free spawning organisms. Since little is known about sperm movement pattern in Mytilus species, the purpose of this study was to investigate sperm movement pattern of blue mussel M. edulis, M. trossulus and their hybrids using computer-assisted sperm movement video analysis. Sperm of all genotypes were found to conduct circular movement in a two-dimensional plane. Furthermore, new sperm movement parameters, real time radius (R), angle change rate (θ) and the center of circular track (O(t)) were developed to verify and quantitatively describe the plane circular movement pattern using software (Image-J) that may be widely applied to sperm movement study in other organisms. Angle change rate was positively correlated to fertilization success. However, no correlation between fertilization and real time radius was detected. Although no interspecific differences were found in the radius, the F1 (first generation) hybrid sperm had a lower angle change rate than M. edulis and M. trossulus. Published studies have shown that sperm circular movement is more prevalent in aquatic broadcast spawning species than in species with mating behavior or internal fertilization. Therefore, a two-dimensional circular movement pattern in sperm may represent a trait that increases fertilization success for broadcast spawning species by either increasing gamete interaction rate at a small scale and/or avoiding swimming further away from the eggs before sperm detects the chemoattractant gradient.     

First record of the microsporidian parasite Steinhausia mytilovum in Mytilus sp. (Bivalvia: Mytilidae) from France

Steinhausia mytilovum is a globally distributed microsporidian parasite which infects the oocytes of the blue mussels Mytilus edulis and M. galloprovincialis. Despite the intensive monitoring effort made on mussel populations, the parasite has not previously been reported in France. We report herein on the occurrence of S. mytilovum in Mytilus sp. from 1 cultured and 2 natural populations on the northern coast of France, thus extending the parasite's known distribution northwards. We also report on the observation in 1989 of S. mytilovum in M. galloprovincialis from the Golfe de Fos area in the Mediterranean Sea (South of France). S. mytilovum was observed in the European hybrid zone between M. edulis and M. galloprovincialis, which therefore renders the exact taxonomic status of the infected hosts unknown. The prevalence of the parasite was low, which suggests that its effect on mussel populations was probably limited.     

Positive selection on an acrosomal sperm protein,M7 lysin,in three species of the mussel genus Mytilus

Marine invertebrate sperm proteins are particularly interesting because they are characterized by positive selection and are likely to be involved in prezyogotic isolation and, thus, speciation. Here, we present the first survey of interspecific and intraspecific variation of a bivalve sperm protein among a group of species that regularly hybridize in nature. M7 lysin is found in sperm acrosomes of mussels and dissolves the egg vitelline coat, permitting fertilization. We sequenced multiple alleles of the mature protein-coding region of M7 lysin from allopatric populations of mussels in the Mytilus edulis species group (M. edulis, M. galloprovincialis, and M. trossulus). A significant McDonald-Kreitman test showed an excess of fixed amino acid replacing substitutions between species, consistent with positive selection. In addition, Kolmogorov-Smirnov tests showed significant heterogeneity in polymorphism to divergence ratios for both synonymous variation and combined synonymous and nonsynonymous variation within M. galloprovincialis. These results indicate that there has been adaptive evolution at M7 lysin and, furthermore, show that positive selection on sperm proteins can occur even when postzygotic reproductive isolation is incomplete.