ASYMMETRIC INTROGRESSION OF MITOCHONDRIAL DNA AMONG EUROPEAN POPULATIONS OF BLUE MUSSELS ( MYTILUS SPP.)

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

First description of a widespread Mytilus trossulus-derived bivalve transmissible cancer lineage in M. trossulus itself

Two lineages of bivalve transmissible neoplasia (BTN), BTN1 and BTN2, are known in blue mussels Mytilus. Both lineages derive from the Pacific mussel M. trossulus and are identified primarily by their unique genotypes of the nuclear gene EF1α. BTN1 is found in populations of M. trossulus from the Northeast Pacific, while BTN2 has been detected in populations of other Mytilus species worldwide but not in M. trossulus itself. Here we examined M. trossulus from the Sea of Japan (Northwest Pacific) for the presence of BTN. Using hemocytology and flow cytometry of the hemolymph, we confirmed the presence of disseminated neoplasia in our specimens. Cancerous mussels possessed the BTN2 EF1α genotype and two mitochondrial haplotypes with different recombinant control regions, similar to that of common BTN2 lineages. This is the first report of BTN2 in its original host species M. trossulus. A comparison of all available BTN and M. trossulus COI sequences suggests a common and recent origin of BTN2 diversity in populations of M. trossulus outside the Northeast Pacific, possibly in the Northwest Pacific.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Mimicry and mitonuclear discordance in nudibranchs: New insights from exon capture phylogenomics

Phylogenetic inference and species delimitation can be challenging in taxonomic groups that have recently radiated and where introgression produces conflicting gene trees, especially when species delimitation has traditionally relied on mitochondrial data and color pattern. Chromodoris, a genus of colorful and toxic nudibranch in the Indo-Pacific, has been shown to have extraordinary cryptic diversity and mimicry, and has recently radiated, ultimately complicating species delimitation. In these cases, additional genome-wide data can help improve phylogenetic resolution and provide important insights about evolutionary history. Here, we employ a transcriptome-based exon capture approach to resolve Chromodoris phylogeny with data from 2,925 exons and 1,630 genes, derived from 15 nudibranch transcriptomes. We show that some previously identified mimics instead show mitonuclear discordance, likely deriving from introgression or mitochondrial capture, but we confirm one "pure" mimic in Western Australia. Sister-species relationships and species-level entities were recovered with high support in both concatenated maximum likelihood (ML) and summary coalescent phylogenies, but the ML topologies were highly variable while the coalescent topologies were consistent across datasets. Our work also demonstrates the broad phylogenetic utility of 149 genes that were previously identified from eupulmonate gastropods. This study is one of the first to (a) demonstrate the efficacy of exon capture for recovering relationships among recently radiated invertebrate taxa, (b) employ genome-wide nuclear markers to test mimicry hypotheses in nudibranchs and (c) provide evidence for introgression and mitochondrial capture in nudibranchs.

Cycles of trans-Arctic dispersal and vicariance, and diversification of the amphi-boreal marine fauna

The amphi-boreal faunal element comprises closely related species and conspecific populations with vicarious distributions in the North Atlantic and North Pacific basins. It originated from an initial trans-Arctic dispersal in the Pliocene after the first opening of the Bering Strait, and subsequent inter-oceanic vicariance through the Pleistocene when the passage through the Arctic was severed by glaciations and low sea levels. Opportunities for further trans-Arctic dispersal have risen at times, however, and molecular data now expose more complex patterns of inter-oceanic affinities and dispersal histories. For a general view on the trans-Arctic dynamics and of the roles of potential dispersal-vicariance cycles in generating systematic diversity, we produced new phylogeographic data sets for amphi-boreal taxa in 21 genera of invertebrates and vertebrates, and combined them with similar published data sets of mitochondrial coding gene variation, adding up to 89 inter-oceanic comparisons involving molluscs, crustaceans, echinoderms, polychaetes, fishes and mammals. Only 39% of the cases correspond to a simple history of Pliocene divergence; in most taxonomical groups, the range of divergence estimates implies connections through the entire Pliocene-Pleistocene-Holocene time frame. Repeated inter-oceanic exchange was inferred for 23 taxa, and the latest connection was usually post-glacial. Such repeated invasions have usually led to secondary contacts and occasionally to widespread hybridization between the different invasion waves. Late- or post-glacial exchange was inferred in 46% of the taxa, stressing the importance of the relatively recent invasions to the current diversity in the North Atlantic. Individual taxa also showed complex idiosyncratic patterns and histories, and several instances of cryptic speciation were recognized. In contrast to a simple inter-oceanic vicariance scenario underlying amphi-boreal speciation, the data expose complex patterns of reinvasion and reticulation that complicate the interpretation of taxon boundaries in the region.

Trans-Atlantic Distribution and Introgression as Inferred from Single Nucleotide Polymorphism: Mussels Mytilus and Environmental Factors

Large-scale climate changes influence the geographic distribution of biodiversity. Many taxa have been reported to extend or reduce their geographic range, move poleward or displace other species. However, for closely related species that can hybridize in the natural environment, displacement is not the only effect of changes of environmental variables. Another option is subtler, hidden expansion, which can be found using genetic methods only. The marine blue mussels Mytilus are known to change their geographic distribution despite being sessile animals. In addition to natural dissemination at larval phase—enhanced by intentional or accidental introductions and rafting—they can spread through hybridization and introgression with local congeners, which can create mixed populations sustaining in environmental conditions that are marginal for pure taxa. The Mytilus species have a wide distribution in coastal regions of the Northern and Southern Hemisphere. In this study, we investigated the inter-regional genetic differentiation of the Mytilus species complex at 53 locations in the North Atlantic and adjacent Arctic waters and linked this genetic variability to key local environmental drivers. Of seventy-nine candidate single nucleotide polymorphisms (SNPs), all samples were successfully genotyped with a subset of 54 SNPs. There was a clear interregional separation of Mytilus species. However, all three Mytilus species hybridized in the contact area and created hybrid zones with mixed populations. Boosted regression trees (BRT) models showed that inter-regional variability was important in many allele models but did not prevail over variability in local environmental factors. Local environmental variables described over 40% of variability in about 30% of the allele frequencies of Mytilus spp. For the 30% of alleles, variability in their frequencies was only weakly coupled with local environmental conditions. For most studied alleles the linkages between environmental drivers and the genetic variability of Mytilus spp. were random in respect to “coding” and “non-coding” regions. An analysis of the subset of data involving functional genes only showed that two SNPs at Hsp70 and ATPase genes correlated with environmental variables. Total predictive ability of the highest performing models (r² between 0.550 and 0.801) were for alleles that discriminated most effectively M. trossulus from M. edulis and M. galloprovincialis, whereas the best performing allele model (BM101A) did the best at discriminating M. galloprovincialis from M. edulis and M. trossulus. Among the local environmental variables, salinity, water temperature, ice cover and chlorophyll a concentration were by far the greatest predictors, but their predictive performance varied among different allele models. In most cases changes in the allele frequencies along these environmental gradients were abrupt and occurred at a very narrow range of environmental variables. In general, regions of change in allele frequencies for M. trossulus occurred at 8–11 psu, 0–10 °C, 60%–70% of ice cover and 0–2 mg m⁻³ of chlorophyll a, M. edulis at 8–11 and 30–35 psu, 10–14 °C and 60%–70% of ice cover and for M. galloprovincialis at 30–35 psu, 14–20 °C.

Twin introductions by independent invader mussel lineages are both associated with recent admixture with a native congener in Australia

Introduced species can impose profound impacts on the evolution of receiving communities with which they interact. If native and introduced taxa remain reproductively semi-isolated, human-mediated secondary contact may promote genetic exchange across newly created hybrid zones, potentially impacting native genetic diversity and invasive species spread. Here, we investigate the contributions of recent divergence histories and ongoing (post-introduction) gene flow between the invasive marine mussel, Mytilus galloprovincialis, and a morphologically indistinguishable and taxonomically contentious native Australian taxon, Mytilus planulatus. Using transcriptome-wide markers, we demonstrate that two contemporary M. galloprovincialis introductions into south-eastern Australia originate from genetically divergent lineages from its native range in the Mediterranean Sea and Atlantic Europe, where both introductions have led to repeated instances of admixture between introduced and endemic populations. Through increased genome-wide resolution of species relationships, combined with demographic modelling, we validate that mussels sampled in Tasmania are representative of the endemic Australian taxon (M. planulatus), but share strong genetic affinities to M. galloprovincialis. Demographic inferences indicate late-Pleistocene divergence times and historical gene flow between the Tasmanian endemic lineage and northern M. galloprovincialis, suggesting that native and introduced taxa have experienced a period of historical isolation of at least 100,000 years. Our results demonstrate that many genomic loci and sufficient sampling of closely related lineages in both sympatric (e.g. Australian populations) and allopatric (e.g. northern hemisphere Mytilus taxa) ranges are necessary to accurately (a) interpret patterns of intraspecific differentiation and to (b) distinguish contemporary invasive introgression from signatures left by recent divergence histories in high dispersal marine species. More broadly, our study fills a significant gap in systematic knowledge of native Australian biodiversity and sheds light on the intrinsic challenges for invasive species research when native and introduced species boundaries are not well defined.

Comparison between single and multi-locus approaches for specimen identification in Mytilus mussels

Mytilus mussels have been the object of much research given their sentinel role in coastal ecosystems and significant value as an aquaculture resource appreciated for both, its flavour and nutritional content. Some of the most-studied Mytilus species are M. edulis, M. galloprovincialis, M. chilensis and M. trossulus. As species identification based on morphological characteristics of Mytilus specimens is difficult, molecular markers are often used. Single-locus markers can give conflicting results when used independently; not all markers differentiate among all species, and the markers target genomic regions with different evolutionary histories. We evaluated the concordance between the PCR-RFLP markers most commonly-used for species identification in mussels within the Mytilus genus (Me15-16, ITS, mac-1, 16S rRNA and COI) when used alone (mono-locus approach) or together (multi-locus approach). In this study, multi-locus strategy outperformed the mono-locus methods, clearly identifying all four species and also showed similar specimen identification performance than a 49 SNPs panel. We hope that these findings will contribute to a better understanding of DNA marker-based analysis of Mytilus taxa. These results support the use of a multi-locus approach when studying this important marine resource, including research on food quality and safety, sustainable production and conservation.

Genetic diversity and connectivity within Mytilus spp. in the subarctic and Arctic

Climate changes in the Arctic are predicted to alter distributions of marine species. However, such changes are difficult to quantify because information on present species distribution and the genetic variation within species is lacking or poorly examined. Blue mussels, Mytilus spp., are ecosystem engineers in the coastal zone globally. To improve knowledge of distribution and genetic structure of the Mytilus edulis complex in the Arctic, we analyzed 81 SNPs in 534 Mytilus spp. individuals sampled at 13 sites to provide baseline data for distribution and genetic variation of Mytilus mussels in the European Arctic. Mytilus edulis was the most abundant species found with a clear genetic split between populations in Greenland and the Eastern Atlantic. Surprisingly, analyses revealed the presence of Mytilus trossulus in high Arctic NW Greenland (77°N) and Mytilus galloprovincialis or their hybrids in SW Greenland, Svalbard, and the Pechora Sea. Furthermore, a high degree of hybridization and introgression between species was observed. Our study highlights the importance of distinguishing between congener species, which can display local adaptation and suggests that information on dispersal routes and barriers is essential for accurate predictions of regional susceptibility to range expansions or invasions of boreal species in the Arctic.

A First Insight into the Genome of the Filter-Feeder Mussel Mytilus galloprovincialis

Mussels belong to the phylum Mollusca, one of the largest and most diverse taxa in the animal kingdom. Despite their importance in aquaculture and in biology in general, genomic resources from mussels are still scarce. To broaden and increase the genomic knowledge in this family, we carried out a whole-genome sequencing study of the cosmopolitan Mediterranean mussel (Mytilus galloprovincialis). We sequenced its genome (32X depth of coverage) on the Illumina platform using three pair-end libraries with different insert sizes. The large number of contigs obtained pointed out a highly complex genome of 1.6 Gb where repeated elements seem to be widespread (~30% of the genome), a feature that is also shared with other marine molluscs. Notwithstanding the limitations of our genome sequencing, we were able to reconstruct two mitochondrial genomes and predict 10,891 putative genes. A comparative analysis with other molluscs revealed a gene enrichment of gene ontology categories related to multixenobiotic resistance, glutamate biosynthetic process, and the maintenance of ciliary structures.

2-DE Mapping of the Blue Mussel Gill Proteome: The Usual Suspects Revisited

The Blue Mussel (Mytilus edulis, L. 1758) is an ecologically important and commercially relevant bivalve. Because of its ability to bioconcentrate xenobiotics, it is also a widespread sentinel species for environmental pollution, which has been used in ecotoxicological studies for biomarker assessment. Consequently, numerous proteomics studies have been carried out in various research contexts using mussels of the genus Mytilus, which intended to improve our understanding of complex physiological processes related to reproduction, adaptation to physical stressors or shell formation and for biomarker discovery. Differential-display 2-DE proteomics relies on an extensive knowledge of the proteome with as many proteoforms identified as possible. To this end, extensive characterization of proteins was performed in order to increase our knowledge of the Mytilus gill proteome. On average, 700 spots were detected on 2-DE gels by colloidal blue staining, of which 122 different, non-redundant proteins comprising 203 proteoforms could be identified by tandem mass spectrometry. These proteins could be attributed to four major categories: (i) “metabolism”, including antioxidant defence and degradation of xenobiotics; (ii) “genetic information processing”, comprising transcription and translation as well as folding, sorting, repair and degradation; (iii) “cellular processes”, such as cell motility, transport and catabolism; (iv) “environmental information processing”, including signal transduction and signalling molecules and interaction. The role of cytoskeleton proteins, energetic metabolism, chaperones/stress proteins, protein trafficking and the proteasome are discussed in the light of the exigencies of the intertidal environment, leading to an enhanced stress response, as well as the structural and physiological particularities of the bivalve gill tissue.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

Glacial history of the European marine mussels Mytilus, inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.

Karyotypic diversification in Mytilus mussels (Bivalvia: Mytilidae) inferred from chromosomal mapping of rRNA and histone gene clusters

Background

Mussels of the genus Mytilus present morphologically similar karyotypes that are presumably conserved. The absence of chromosome painting probes in bivalves makes difficult verifying this hypothesis. In this context, we comparatively mapped ribosomal RNA and histone gene families on the chromosomes of Mytilus edulis, M. galloprovincialis, M. trossulus and M. californianus by fluorescent in situ hybridization (FISH).

Results

Major rRNA, core and linker histone gene clusters mapped to different chromosome pairs in the four taxa. In contrast, minor rRNA gene clusters showed a different behavior. In all Mytilus two of the 5S rDNA clusters mapped to the same chromosome pair and one of them showed overlapping signals with those corresponding to one of the histone H1 gene clusters. The overlapping signals on mitotic chromosomes became a pattern of alternate 5S rRNA and linker histone gene signals on extended chromatin fibers. Additionally, M. trossulus showed minor and major rDNA clusters on the same chromosome pair.

Conclusion

The results obtained suggest that at least some of the chromosomes bearing these sequences are orthologous and that chromosomal mapping of rRNA and histone gene clusters could be a good tool to help deciphering some of the many unsolved questions in the systematic classification of Mytilidae.

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.

Physiological energetics and biogeographic range limits of three congeneric mussel species

Closely related species with different physiological tolerances and distributions make ideal systems for documenting range shifts in response to a changing climate. Mytilus edulis, M. trossulus, and M. galloprovincialis are sibling species of marine mussels with distinct biogeographical ranges that are correlated with sea surface temperatures. We determined the scope for growth of these three species at a range of temperatures to determine if energetics could predict their distributions. Scope for growth (SFG) represents energy available for growth and/or reproduction above that necessary for maintenance requirements. The SFG of M. galloprovincialis, the species known to inhabit the warmest habitats, was shifted towards warmer temperatures compared to the other two species, remaining positive until nearly 30 °C. M. edulis, a cold-temperate species, maintained a positive SFG up to 23 °C. M. trossulus, a boreal species, generally was not able to maintain a positive SFG above 17 °C. The warm end of each species' range correlated strongly with the point at which that species' SFG became negative in summer and fall. Energetics at cold temperatures did not predict the cold end of the species' ranges, as there was no clear SFG advantage to explain the dominance of M. trossulus in cold habitats. As sea surface temperatures continue to warm with climate change, the energetics of these three species provide a basis for developing mechanistic models predicting future distribution and productivity changes in mussel populations.

Reticulate evolution and marine organisms: the final frontier?

The role that reticulate evolution (i.e., via lateral transfer, viral recombination and/or introgressive hybridization) has played in the origin and adaptation of individual taxa and even entire clades continues to be tested for all domains of life. Though falsified for some groups, the hypothesis of divergence in the face of gene flow is becoming accepted as a major facilitator of evolutionary change for many microorganisms, plants and animals. Yet, the effect of reticulate evolutionary change in certain assemblages has been doubted, either due to an actual dearth of genetic exchange among the lineages belonging to these clades or because of a lack of appropriate data to test alternative hypotheses. Marine organisms represent such an assemblage. In the past half-century, some evolutionary biologists interested in the origin and trajectory of marine organisms, particularly animals, have posited that horizontal transfer, introgression and hybrid speciation have been rare. In this review, we provide examples of such genetic exchange that have come to light largely as a result of analyses of molecular markers. Comparisons among these markers and between these loci and morphological characters have provided numerous examples of marine microorganisms, plants and animals that possess the signature of mosaic genomes.

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.

Molecular population genetics of the male and female mitochondrial DNA molecules of the California sea mussel, Mytilus californianus

The presence of two gender-associated mitochondrial genomes in marine mussels provides a unique opportunity to investigate the dynamics of mtDNA evolution without complications inherent in interspecific comparisons. Here, we assess the relative importance of selection, mutation, and differential constraint in shaping the patterns of polymorphism within and divergence between the male (M) and female (F) mitochondrial genomes of the California sea mussel, Mytilus californianus. Partial sequences were obtained from homologous regions of four genes (nad2, cox1, atp6, and nad5) totaling 2307 bp in length. The M and F mtDNA molecules of M. californianus exhibited extensive levels of nucleotide polymorphism and were more highly diverged than observed in other mytilids (overall Tamura-Nei distances >40%). Consistent with previous studies, the M molecule had significantly higher levels of silent and replacement polymorphism relative to F. Both genomes possessed large numbers of singleton and low-frequency mutations that gave rise to significantly negative Tajima's D values. Mutation-rate scalars estimated for silent and replacement mutations were elevated in the M genome but were not sufficient to account for its higher level of polymorphism. McDonald-Kreitman tests were highly significant at all loci due to excess numbers of fixed replacement mutations between molecules. Strong purifying selection was evident in both genomes in keeping the majority of replacement mutations at low population frequencies but appeared to be slightly relaxed in M. Our results suggest that a reduction in selective constraint acting on the M genome remains the best explanation for its greater levels of polymorphism and faster rate of evolution.

Genetic Divergence Detected by ISSR Markers and Characterization of Microsatellite Regions in Mytilus Mussels

The wide distribution of microsatellites in mussels of the Mytilus edulis complex (Mytilidae) enables the analysis of inter-simple-sequence repeat (ISSR) markers. The aim of this investigation was to assess genetic differentiation in six sampling localities distributed along the European Atlantic coast to expose the potential of these markers in genetic studies requiring the detection of low polymorphism and as a source of sequences for developing microsatellite markers. We detected low genetic structuring within each member of the Mytilus edulis complex. Nei and Li distances and AMOVA clustered the individuals studied into two groups. On the basis of these results two sampling localities coming from the M. edulis x M. galloprovincialis mosaic hybrid zone in Western Europe were assigned to one species. On the other hand, mussels of a sampling locality in the Baltic Sea were not significantly different from a pure M. edulis locality supporting an extensive introgression of M. edulis in these individuals. Finally, 148 microsatellites were found in the sequences of 51 ISSR markers, and two polymorphic microsatellite markers were developed.

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.

Nonhomologous recombination between the large unassigned region of the male and female mitochondrial genomes in the mussel, Mytilus trossulus

Doubly uniparental inheritance of mtDNA (DUI) is commonly observed in several genera of bivalves. Under DUI, female offspring inherit mtDNA from their mothers, while male offspring inherit mtDNA from both parents but preferentially transmit the paternally inherited mtDNA to their sons. Several studies have shown that the female- and male-specific mtDNA lineages in blue mussels, Mytilus spp., vary by upward of 20% at the nucleotide level. In addition to high levels of nucleotide substitution, the present study observed substantial gender-based length polymorphism in the presumptive mitochondrial control region (=large unassigned region; LUR) of North American M. trossulus. In this species, female lineage LUR haplotypes are over 2 kb larger than male lineage LUR haplotypes. Analysis of sequence data for these length variants indicates that the F LUR haplotypes of North American M. trossulus contain sequences similar to the F lineage control region in the congeners M. edulis and M. galloprovincialis. Relative to the F LUR in the latter two species, however, the F lineage LUR haplotypes in M. trossulus contain two large sequence insertions, each nearly 1 kb in size. One of these insertions has high sequence similarity to the male lineage LUR of M. trossulus. The tandem arrangement of F and M control region sequences in the F lineage LUR of M. trossulus is most likely the result of nonhomologous recombination between the male and the female mitochondrial genomes in M. trossulus, a finding that has important implications regarding the transmission and evolution of blue mussel mitochondrial genomes.

INVITED REVIEW: Local adaptation and species segregation in two mussel (Mytilus edulis × Mytilus trossulus) hybrid zones

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and Mytilus galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; Mytilus edulis and Mytilus 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.

Heteroplasmy suggests paternal co-transmission of multiple genomes and pervasive reversion of maternally into paternally transmitted genomes of mussel (Mytilus) mitochondrial DNA

Marine mussels of the genus Mytilus have two types of mitochondrial DNA with separate paternal and maternal inheritance. Females are homoplasmic for an F genome that is transmitted to all offspring, whereas males are heteroplasmic for this F genome and for a highly diverged (> 20%) M genome that is transmitted only to sons. Here we provide phylogenetic evidence based on lrRNA sequence data that most of the paternal genomes in European M. trossulus have an introgressive female M. edulis origin and are nearly indistinguishable in sequence from F types of M. trossulus. This observation is best explained by the hypothesis that introgressed F type molecules have recently invaded the paternal route and have assumed the role of M molecules, then resetting to zero the time of sequence divergence between M and F lineages. European M. trossulus shows a high prevalence of males heteroplasmic for three different mitochondrial DNA types all having the same two paternal types and the same maternal type, consistent with paternal co-transmission of multiple genomes. Co-transmission of the same genomes must apparently operate uninterruptedly for several generations in spite of the very different evolutionary origin of the specific molecules that are transmitted paternally and maternally in European M. trossulus.

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.

Disjunct distribution of highly diverged mitochondrial lineage clade and population subdivision in a marine bivalve with pelagic larval dispersal

Mitochondrial DNA sequence data for 295 individuals of the marine bivalve Macoma balthica (L.) were collected from 10 sites across the European distribution, and from Alaska. The data were used to infer population subdivision history and estimate current levels of gene flow. Inferred historical biogeography was expected to be congruent with colonization of the Atlantic Ocean from the Pacific Ocean after the opening of the Bering Strait 3.5 Ma. In addition, the last glacial maximum, about 18000 years ago, was expected to have been responsible for most of the present-day distribution of molecular variation within Europe, because the area must have been recolonized after confinement to France and the south of the British Isles during the last glacial maximum. Current gene flow was hypothesized to be high, because the larvae of M. balthica spend 2-5 weeks drifting in the water column. The geographical distribution of one highly diverged haplotype clade was found to be disjunct and was encountered exclusively in samples from the Baltic Sea and Alaska. A molecular clock calibration for marine bivalve cytochrome-c-oxidase I dates this clade as having split off from the other haplotypes 9.8-39 Ma. Multiple colonizations of the Atlantic Ocean from the Pacific by M. balthica may explain the strong differences found between Baltic Sea and other European populations of this species. The sympatric occurrence of the highly diverged mitochondrial lineages in western parts of the Baltic Sea points to secondary admixture. With the use of coalescent analysis, population divergence times for French vs. other non-Baltic European populations ('Atlantic population assemblage') were estimated at a minimum of about 110000 years ago, well before the last glacial maximum 18000 years ago. Signatures of population divergence of M. balthica that appear to have originated during the Pleistocene have thus survived the last glacial maximum. Some of the populations within the Atlantic assemblage are currently isolated, while others appear to be connected by gene flow. Apparently, populations of this species can remain highly subdivided in spite of the potential for high gene flow, implying that their population and evolutionary dynamics can be independent.

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.

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.

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.