Genetic, Ecological and Morphological Distinctness of the Blue Mussels Mytilus trossulus Gould and M. edulis L. in the White Sea

Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae)

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

In vitro proliferation of Mytilus edulis male germ cell progenitors

Our understanding of basic cellular processes has mostly been provided by mammalian cell culture, and by some non-mammalian vertebrate and few invertebrate cell culture models. Developing reliable culture conditions for non-model organisms is essential to allow investigation of more unusual cellular processes. Here, we investigate how cells isolated from different tissues of the marine mussel Mytilus edulis thrive and survive in vitro in the hope of establishing a suitable laboratory model for the investigation of cellular mechanisms specific to these bivalve mollusks. We found that cells dissociated from mantle tissue attached to the culture vessels and proliferated well in vitro, whereas cells isolated from gills, although remaining viable, did not maintain divisions over three to four weeks in culture. We used antibodies against the germ-line marker DEAD-box helicase 4 (DDX4), also known as VASA, and the epithelial cell marker cytokeratin to distinguish different cell types in culture. DDX4-positive cells were predominant in 25-day-old cultures from male mantles. Cells from other tissues remained in low numbers and did not seem to change in composition over time. Overall, the culture conditions described here allow an efficient selection of male germ cells that could be used to study specific cellular mechanisms in vitro.

SNP genotyping revealed a hybrid zone between Mytilus chilensis and M. platensis in southern South America (the Strait of Magellan, Isla Grande de Tierra del Fuego and the Falkland Islands)

Hybrid zones among mussel species have been extensively studied in the northern hemisphere. In South America, it has only recently become possible to study the natural hybrid zones, due to the clarification of the taxonomy of native mussels of the Mytilus genus. Analysing 54 SNP markers, we show the genetic species composition and admixture in the hybrid zone between M. chilensis and M. platensis in the southern end of South America. Bayesian, non-Bayesian clustering and re-assignment algorithms showed that the natural hybrid zone between M. chilensis and M. platensis in the Strait of Magellan, Isla Grande de Tierra del Fuego and the Falkland Islands shows clinal architecture. The hybrid zone can be divided into three different areas: the first one is on the Atlantic coast where only pure M. platensis and hybrid were found. In the second one, inside the Strait of Magellan, pure individuals of both species and mussels with variable degrees of hybridisation coexist. In the last area at the Strait in front of Punta Arenas City, fjords on the Isla Grande de Tierra del Fuego, and at the Beagle Channel, only M. chilensis and a low number of hybrids were found. According to the proportion of hybrids, bays with protected conditions away from strong currents would give better conditions for hybridisation. We do not find evidence of any other mussel species such as M. edulis, M. galloprovincialis, M. planulatus or M. trossulus in the zone.

Species identification based on a semi-diagnostic marker: Evaluation of a simple conchological test for distinguishing blue mussels Mytilus edulis L. and M. trossulus Gould

Cryptic and hybridizing species may lack diagnostic taxonomic characters leaving researchers with semi-diagnostic ones. Identification based on such characters is probabilistic, the probability of correct identification depending on the species composition in a mixed population. Here we test the possibilities of applying a semi-diagnostic conchological character for distinguishing two cryptic species of blue mussels, Mytilus edulis and M. trossulus. These ecologically, stratigraphically and economically important molluscs co-occur and hybridize in many areas of the North Atlantic and the neighboring Arctic. Any cues for distinguishing them in sympatry without genotyping would save much research effort. Recently these species have been shown to statistically differ in the White Sea, where a simple character of the shell was used to distinguish two mussel morphotypes. In this paper, we analyzed the associations between morphotypes and species-specific genotypes based on an abundant material from the waters of the Kola Peninsula (White Sea, Barents Sea) and a more limited material from Norway, the Baltic Sea, Scotland and the Gulf of Maine. The performance of the "morphotype test" for species identification was formally evaluated using approaches from evidence-based medicine. Interspecific differences in the morphotype frequencies were ubiquitous and unidirectional, but their scale varied geographically (from 75% in the White Sea to 15% in the Baltic Sea). In addition, salinity-related variation of this character within M. edulis was revealed in the Arctic Barents Sea. For every studied region, we established relationships between the proportions of the morphotypes in the populations as well as between the proportions of the morphotypes in samples and the probabilities of mussels of different morphotypes being M. trossulus and M. edulis. We provide recommendations for the application of the morphotype test to mussels from unstudied contact zones and note that they may apply equally well to other taxa identified by semi-diagnostic traits.

Mytilus trossulus introgression and consequences for shell traits in longline cultivated mussels

Mussels belonging to the Mytilus species complex (M. edulis, ME; M. galloprovincialis, MG; and M. trossulus, MT) often occur in sympatry, facilitating introgressive hybridization. This may be further promoted by mussel aquaculture practices, with MT introgression often resulting in commercially unfavourable traits such as low meat yield and weak shells. To investigate the relationship between genotype and shell phenotype, genetic and morphological variability was quantified across depth (1 m to 7 m) along a cultivation rope at a mussel farm on the West coast of Scotland. A single nuclear marker (Me15/16) and a novel panel of 33 MT-diagnostic single nucleotide polymorphisms were used to evaluate stock structure and the extent of MT introgression across depth. Variation in shell strength, determined as the maximum compression force for shell puncture, and shell shape using geometric morphometric analysis were evaluated in relation to cultivation depth and the genetic profiles of the mussels. Overall, ME was the dominant genotype across depth, followed by ME × MG hybrids and smaller quantities of ME × MT hybrids and pure MT individuals. In parallel, we identified multiple individuals that were either predominantly homozygous or heterozygous for MT-diagnostic alleles, likely representing pure MT and first-generation ME × MT hybrids, respectively. Both the proportion of individuals carrying MT alleles and MT allele frequency declined with depth. Furthermore, MT-introgressed individuals had significantly weaker and more elongate shells than nonintrogressed individuals. This study provides detailed insights into stock structure along a cultivation rope and suggests that practical methods to assess shell strength and shape of cultivated mussels may facilitate the rapid identification of MT, limiting the impact of this commercially damaging species.

Some Examples of the Use of Molecular Markers for Needs of Basic Biology and Modern Society

Application of molecular genetic markers appeared to be very fruitful in achieving many goals, including (i) proving the theoretic basements of general biology and (ii) assessment of worldwide biodiversity. Both are provided in the present meta-analysis and a review as the main signal. One of the basic current challenges in modern biology in the face of new demands in the 21st century is the validation of its paradigms such as the synthetic theory of evolution (STE) and biological species concept (BSC). Another of most valuable goals is the biodiversity assessment for a variety of social needs including free web-based information resources about any living being, renovation of museum collections, nature conservation that recognized as a global project, iBOL, as well as resolving global trading problems such as false labeling of species specimens used as food, drug components, entertainment, etc. The main issues of the review are focused on animals and combine four items. (1) A combination of nDNA and mtDNA markers best suits the identification of hybrids and estimation of genetic introgression. (2) The available facts on nDNA and mtDNA diversity seemingly make introgression among many taxa obvious, although it is evident, that introgression may be quite restricted or asymmetric, thus, leaving at least the "source" taxon (taxa) intact. (3) If we consider sexually reproducing species in marine and terrestrial realms introgressed, as it is still evident in many cases, then we should recognize that the BSC, in view of the complete lack of gene flow among species, is inadequate because many zoological species are not biological ones yet. However, vast modern molecular data have proven that sooner or later they definitely become biological species. (4) An investigation into the fish taxa divergence using the BOLD database shows that most gene trees are basically monophyletic and interspecies reticulations are quite rare.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

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.

Replicated anthropogenic hybridisations reveal parallel patterns of admixture in marine mussels

Human-mediated transport creates secondary contacts between genetically differentiated lineages, bringing new opportunities for gene exchange. When similar introductions occur in different places, they provide informally replicated experiments for studying hybridisation. We here examined 4,279 Mytilus mussels, sampled in Europe and genotyped with 77 ancestry-informative markers. We identified a type of introduced mussels, called "dock mussels," associated with port habitats and displaying a particular genetic signal of admixture between M. edulis and the Mediterranean lineage of M. galloprovincialis. These mussels exhibit similarities in their ancestry compositions, regardless of the local native genetic backgrounds and the distance separating colonised ports. We observed fine-scale genetic shifts at the port entrance, at scales below natural dispersal distance. Such sharp clines do not fit with migration-selection tension zone models, and instead suggest habitat choice and early-stage adaptation to the port environment, possibly coupled with connectivity barriers. Variations in the spread and admixture patterns of dock mussels seem to be influenced by the local native genetic backgrounds encountered. We next examined departures from the average admixture rate at different loci, and compared human-mediated admixture events, to naturally admixed populations and experimental crosses. When the same M. galloprovincialis background was involved, positive correlations in the departures of loci across locations were found; but when different backgrounds were involved, no or negative correlations were observed. While some observed positive correlations might be best explained by a shared history and saltatory colonisation, others are likely produced by parallel selective events. Altogether, genome-wide effect of admixture seems repeatable and more dependent on genetic background than environmental context. Our results pave the way towards further genomic analyses of admixture, and monitoring of the spread of dock mussels both at large and at fine spacial scales.

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.

DNA barcodes of Antipode marine invertebrates in Bay of Biscay and Gulf of Lion ports suggest new biofouling challenges

Marine biological invasions threaten global biodiversity nowadays. In this article, we have studied fouling communities from 10 port areas of south Bay of Biscay (Atlantic Ocean) and Gulf of Lion (Mediterranean Sea). A total of 834 individuals were genetically barcoded and corresponded to 95 different species. A total of 76 native species 8 genera and 1 family were identified, 58 from the Bay of Biscay and 23 from the Gulf of Lion. Furthermore, 19 species were identified as non-indigenous or cryptogenic (18 from the Bay of Biscay and 4 from the Gulf of Lion). We found a high proportion of Antipode non-indigenous species (NIS) that represented the 19.3% of all sampled individuals and the 54.21% of NIS specimens of this study. A framework for inference of donor regions based on a phylogenetic screening of genetic sequences was proposed as a proof of concept and tested, as well as models for the relationship between NIS introductions, maritime imports and distance to NIS native range and inferred donor areas. Consistent generalized linear models (GLM) with positive association between NIS genetic diversity and distance, not with maritime growth weight imports, strongly suggest that distant NIS could pose higher invasion risk than closer species. Selection for wider tolerance ranges during the long travel -direct or stepwise, as well as environmental similarity between donor and receiving regions, may explain these results.

Discriminating Eaters: Sea Stars Asterias rubens L. Feed Preferably on Mytilus trossulus Gould in Mixed Stocks of Mytilus trossulus and Mytilus edulis L

Sea stars Asterias rubens are important natural enemies of the blue mussel Mytilus in the North Atlantic. We asked whether these predators distinguish between the cryptic species M. edulis and M. trossulus that occur sympatrically in the White Sea. In mixed experimental stocks, the odds of being eaten by sea stars were about four times greater for M. trossulus. We also showed that A. rubens preferred smaller mussels to larger ones, irrespective of their species affinity. Our findings support earlier indirect observations showing that sea stars recognize M. trossulus as a more preferable prey than M. edulis. Dramatic differences in the vulnerability to sea star predation may explain the segregation of habitats between the two mussel species in contact zones; M. trossulus usually tends to occupy habitats where the sea star predators are scarce.

Linking genotype to phenotype in a changing ocean: inferring the genomic architecture of a blue mussel stress response with genome-wide association

A key component to understanding the evolutionary response to a changing climate is linking underlying genetic variation to phenotypic variation in stress response. Here, we use a genome-wide association approach (GWAS) to understand the genetic architecture of calcification rates under simulated climate stress. We take advantage of the genomic gradient across the blue mussel hybrid zone (Mytilus edulis and Mytilus trossulus) in the Gulf of Maine (GOM) to link genetic variation with variance in calcification rates in response to simulated climate change. Falling calcium carbonate saturation states are predicted to negatively impact many marine organisms that build calcium carbonate shells - like blue mussels. We sampled wild mussels and measured net calcification phenotypes after exposing mussels to a 'climate change' common garden, where we raised temperature by 3°C, decreased pH by 0.2 units and limited food supply by filtering out planktonic particles >5 μm, compared to ambient GOM conditions in the summer. This climate change exposure greatly increased phenotypic variation in net calcification rates compared to ambient conditions. We then used regression models to link the phenotypic variation with over 170 000 single nucleotide polymorphism loci (SNPs) generated by genotype by sequencing to identify genomic locations associated with calcification phenotype, and estimate heritability and architecture of the trait. We identified at least one of potentially 2-10 genomic regions responsible for 30% of the phenotypic variation in calcification rates that are potential targets of natural selection by climate change. Our simulations suggest a power of 13.7% with our study's average effective sample size of 118 individuals and rare alleles, but a power of >90% when effective sample size is 900.

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.