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

Unlocking the secret life of blue mussels: Exploring connectivity in the Skagerrak through biophysical modeling and population genomics

Knowledge of functional dispersal barriers in the marine environment can be used to inform a wide variety of management actions, such as marine spatial planning, restoration efforts, fisheries regulations, and invasive species management. Locations and causes of dispersal barriers can be studied through various methods, including movement tracking, biophysical modeling, demographic models, and genetics. Combining methods illustrating potential dispersal, such as biophysical modeling, with realized dispersal through, e.g., genetic connectivity estimates, provides particularly useful information for teasing apart potential causes of observed barriers. In this study, we focus on blue mussels (Mytilus edulis) in the Skagerrak-a marginal sea connected to the North Sea in Northern Europe-and combine biophysical models of larval dispersal with genomic data to infer locations and causes of dispersal barriers in the area. Results from both methods agree; patterns of ocean currents are a major structuring factor in the area. We find a complex pattern of source-sink dynamics with several dispersal barriers and show that some areas can be isolated despite an overall high dispersal capability. Finally, we translate our finding into management advice that can be used to sustainably manage this ecologically and economically important species in the future.

Coastal Mussel (Mytilus spp.) Soft Tissues as Bioindicators of Methylmercury: Exploring the Relationship Between Condition Index and Methylmercury Concentrations

We investigated total mercury (THg) and methylmercury (MeHg) concentrations in coastal mussels (Mytilus spp.) sampled from the Minas Basin, Bay of Fundy and evaluated the relationship with condition index (CI). THg concentrations were low in sediment (mean THg = 5.15 ± 2.11 ng/g dw; n = 6) and soft tissues (mean THg = 62.3 ± 13.7 ng/g; mean MeHg = 13.2 ± 6.3 ng/g; n = 57). The THg in tissues had no significant relationship with CI (Rs= -0.205, p = 0.126). MeHg in tissues were significantly and negatively correlated with condition index (Rs = -0.361, p = 0.006) indicating that healthier mussels (higher CI) have lower mercury content possibly due to elimination strategies or growth dilution.

Linking Acrosome Size and Genetic Divergence in an Inter-Oceanic Mussel from the Pacific and Atlantic Coasts: A Case of Incipient Speciation?

In recent years, advances in analyses of the sperm morphology and genetics of Perumytilus purpuratus have allowed to two evolutionary scenarios for this mussel to be suggested: (1) the scenario of cryptic species and (2) the scenario of incipient or in progress speciation. For a better understanding of the evolutionary history of P. purpuratus, we performed extensive sampling along a latitudinal gradient of ca. 7180 km of coastline-from the Southern Pacific Ocean to the Atlantic Ocean-and we delved deeper into the sperm morphology of P. purpuratus, exploring its association with the phylogeny and population genetics to determine whether the variability in sperm traits between the northern and southern regions was a signal of cryptic or incipient species. Overall, our results showed that sperm sizes were strongly correlated with the genetic structure in males of P. purpuratus. We identified at 37° S on the Pacific coast a coincident break of both sperm size and genetic disruption that can be explained by historical events and postglacial recolonization as causal phenomena for the observed divergences. Furthermore, evidence of genetic admixture between lineages was found at 38° S, suggesting the presence of an introgressive hybridization zone and incomplete reproductive isolation in an in fraganti or incipient speciation process.

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

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

Development of a novel molecular tool to study molecular ecology of Ornithomya (Hippoboscidae) avian louse flies

Louse flies (Diptera: Hippoboscidae) are obligatory hematophagous ectoparasites of birds and mammals. These widely distributed parasitic flies may have a significant impact on wild and farm animals by feeding on their blood and transmitting bloodborne pathogens. However, despite their ecological importance, louse flies are clearly underrepresented in host-parasite research and implementation of genetic approaches in this group is generally hampered by lacking molecular tools. In addition, louse flies that parasitize long-distance migrants can travel long distances with their avian hosts, facilitating the large-scale spread of pathogens across landscapes and geographic regions. Given the wide diversity of louse flies that parasitize a variety of avian hosts, their direct negative impact on host survival, and their high potential to transmit bloodborne pathogens even along avian migration routes, it is surprising that our knowledge of louse fly ecology is rather modest and incomplete. Here, we aimed to develop a novel molecular tool for polyxenous avian louse flies from the genus Ornithomya, which are among the most common and widely distributed representatives of Hippoboscidae family, to improve research of their genetic population structure and molecular ecology. Using the Illumina Mi-seq sequencing, we conducted a genome-wide scan in Ornithomya avicularia to identify putative microsatellite markers. A panel of 26 markers was selected to develop amplification protocols and assess polymorphism in the Central European population of O. avicularia, as well as to test for cross-amplification in a congeneric species (O. chloropus). A genome-scan in O. avicularia identified over 12 thousand putative microsatellite markers. Among 26 markers selected for a population-wide screening; one did not amplify successfully and three were monomorphic. 22 markers were polymorphic with at least two alleles detected. Two markers showed presence of null alleles. A cross-amplification of microsatellite markers in O. chloropus revealed allelic polymorphism at 14 loci, with the mean allelic richness of 3.78 alleles per locus (range: 2-8). Our genome-wide scan in O. avicularia provides a novel and powerful tool for molecular research in Ornithomya louse flies. Our panel of polymorphic microsatellite loci should allow genotyping of louse flies from geographically distinct populations and from a wide spectrum of avian hosts, enhancing population genetic and phylogeographic research in Ornithomya.

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Where have all the beads gone? Fate of microplastics in a closed exposure system and their effects on clearance rates in Mytilus spp

The fate of microplastic particles (MP) in exposure experiments is mostly unclear. We measured the recovery of polystyrene (PS) microbeads, which were applied in various concentrations from 0.07 to 47.47 beads/ml, from the different compartments of an experimental system with mussels (Mytilus spp.). At the end of the experiment, we detected a significant loss of MP indicating that the mussels were exposed to less particles than intended. If such a discrepancy remains un-recognized by the experimenter, observed effects are related to an inaccurate particle concentration. Additionally, we observed reduced clearance rates of the mussels in the presence of MP and the effect size increased with increasing particle concentration. This effect was more pronounced in mussels that had recently spawned than in mussels that still had mature gonads. This is a hint that effects of MP may depend on the reproductive status of an organism.

Evidence of multiple genome duplication events in Mytilus evolution

Background

Molluscs remain one significantly under-represented taxa amongst available genomic resources, despite being the second-largest animal phylum and the recent advances in genomes sequencing technologies and genome assembly techniques. With the present work, we want to contribute to the growing efforts by filling this gap, presenting a new high-quality reference genome for Mytilus edulis and investigating the evolutionary history within the Mytilidae family, in relation to other species in the class Bivalvia.

Results

Here we present, for the first time, the discovery of multiple whole genome duplication events in the Mytilidae family and, more generally, in the class Bivalvia. In addition, the calculation of evolution rates for three species of the Mytilinae subfamily sheds new light onto the taxa evolution and highlights key orthologs of interest for the study of Mytilus species divergences.

Conclusions

The reference genome presented here will enable the correct identification of molecular markers for evolutionary, population genetics, and conservation studies. Mytilidae have the capability to become a model shellfish for climate change adaptation using genome-enabled systems biology and multi-disciplinary studies of interactions between abiotic stressors, pathogen attacks, and aquaculture practises.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08575-9.

Plastic and natural inorganic microparticles do not differ in their effects on adult mussels (Mytilidae) from different geographic regions

Microplastics are ubiquitous in the marine environment and studies on their effects on benthic filter feeders at least partly revealed a negative influence. However, it is still unclear whether the effects of microplastics differ from those of natural suspended microparticles, which constitute a common stressor in many coastal environments. We present a series of experiments that compared the effects of six-week exposures of marine mussels to two types of natural particles (red clay and diatom shells) to two types of plastic particles (Polymethyl Methacrylate and Polyvinyl Chloride). Mussels of the family Mytilidae from temperate regions (Japan, Chile, Tasmania) through subtropical (Israel) to tropical environments (Cabo Verde) were exposed to concentrations of 1.5 mg/L, 15 mg/L and 150 mg/L of the respective microparticles. At the end of this period, we found significant effects of suspended particles on respiration rate, byssus production and condition index of the animals. There was no significant effect on clearance rate and survival. Surprisingly, we observed only small differences between the effects of the different types of particles, which suggests that the mussels were generally equally robust towards exposure to variable concentrations of suspended solids regardless of whether they were natural or plastic. We conclude, that microplastics and suspended solids elicit similar effects on the tested response variables, and that both types of microparticles mainly cause acute responses rather than more persistent carry-over effects.

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.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Nonconsumptive Predator Effects on Prey Demography: Recent Advances Using Intertidal Invertebrates

Predators influence prey demography through consumption, but the mere presence of predators may trigger behavioural changes in prey that, if persistent or intense, may also influence prey demography. A tractable system to study such nonconsumptive effects (NCEs) of predators involves intertidal invertebrates. This mini review summarises recent research using barnacles and mussels as prey and dogwhelks as predators. The field manipulation of dogwhelk density revealed that pelagic barnacle larvae avoid benthic settlement near dogwhelks, which limits barnacle recruitment, a relevant outcome because recruitment is the only source of population replenishment for barnacles, as they are sessile. This avoidance behaviour is likely triggered by waterborne dogwhelk cues and may have evolved to limit future predation risk. Increasing densities of barnacle recruits and adults can prevent such NCEs from occurring, seemingly because benthic barnacles attract conspecific larvae through chemical cues. Barnacle recruit density increased with the abundance of coastal phytoplankton (food for barnacle larvae and recruits), so barnacle food supply seems to indirectly limit dogwhelk NCEs. By inhibiting barnacle feeding, dogwhelk cues also limited barnacle growth and reproductive output. Wave action weakens dogwhelk NCEs likely through hydrodynamic influences. Dogwhelk cues also limit mussel recruitment, as mussel larvae also exhibit predator avoidance behaviour. The NCEs on recruitment are weaker for mussels than for barnacles, possibly because mussel larvae can detach themselves after initial settlement, an ability that barnacle larvae lack. Overall, these field experiments provide evidence of predator NCEs on prey demography for coastal marine systems.

Natural distribution of pure and hybrid Mytilus sp. along the south Mediterranean and North-east Atlantic coasts and sensitivity of D-larvae stages to temperature increases and metal pollution

The distribution of the Mediterranean mussel Mytilus galloprovincialis extends more and more northwards in the Atlantic. Crossings are frequently observed with the blue mussel Mytilus edulis along the French and English coasts. The aim of this study is firstlyto identify the co-presence of M. galloprovincialis, M. edulis, and their hybrids in different sites of the Mediterranean and Atlantic coasts, and to provide insights for the thermal tolerance and toxicant susceptibility of Mytilus edulis, Mytilus galloprovincialis and their hybrids. Mussels were collected from the shore at 20 sampling sitesin Europe and Tunisia and identified using Me 15/16 primers targeting the adhesive protein gene sequence. Samples were screened for the presence of Mytilus edulis, Mytilus galloprovincialis, and hybrids alleles using PCR. To get more information on hybrids sensitivities to temperature and metals, freshly fertilized eggs of the two species and their hybrids were reared at four temperatures 18, 20, 22, and 24 °C and exposed to concentrations of Cu, Ag, and a mixture of both metals. Arrests of development and malformations were recorded after 48 h of exposure. The genotypic identification of the two species on 20 sites of the Mediterranean and Atlantic coasts carried out during this study confirms the presence of pure and hybrid species of mussel. Our results highlighted that hybrid larvae from a female of M. galloprovincialis are significantly more tolerant to temperature increases than pure larvae of M. galloprovincialis and pure and hybrid larvae of M. edulis. No significant interspecies-differences of sensitivity were noted for metal exposure alone. However, a co-exposure of larvae to both metal and high temperature highlighted the higher tolerance of hybrid larvae from a female of M. galloprovincialis to both stresses. The overall results could allow the prediction of the future evolution of mussel populations facing environmental changes.

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.

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

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

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.

Molecular mechanisms underlying the effects of temperature increase on Mytilus sp. and their hybrids at early larval stages

The present work aims to investigate the effects of water temperature increase on Mytilus galloprovincilis and Mytilus edulis pure larvae (PG, PE) and their hybrids (HFG, HFE). D-larvae were maintained at 18 °C or exposed to a higher temperature of 22 °C for 48 h. Initially, Embryotoxicity test was evaluated. Second, a transcriptomic analysis using a recently developed microarray platform was applied to determine the main biological processes involved in early life stages responses to temperature increase. Finally, an immunofluorescence investigation was performed to bridge the gap between transcriptomic regulation and the real changes at cellular/tissue levels. Embryotoxicity test revealed a higher sensitivity of M. edulis (PE) D-larvae as well as hybrids from females M. edulis (HFE) to temperature increase, with the highest rate of larval malformations. Transcriptomic results indicated a lack of an adequate heat shock protein (Hsp) response in PE and HFE larvae (the high expression was observed in PG larvae); the differential expression of gene involved in translation, energy metabolism and oxidative stress response may contribute to explain the observed complex alterations in the studied conditions. As revealed by immunohistochemistry, cytoskeleton proteins changes associated with a drastic decrease of Histidine-Rich Glycoprotein (HRG) may elucidate the larval abnormalities in shell development observed for PE and HFE larvae. Overall, the results indicate that each type of pure larva (PG and PE) and their respective female hybrid (HFG and HFE) react similarly to the temperature increase. Our data should be carefully considered in view of the water temperature increase in marine ecosystems and especially for the mussel's species in confluence zones.

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.

Biomineralization plasticity and environmental heterogeneity predict geographical resilience patterns of foundation species to future change

Although geographical patterns of species' sensitivity to environmental changes are defined by interacting multiple stressors, little is known about compensatory processes shaping regional differences in organismal vulnerability. Here, we examine large-scale spatial variations in biomineralization under heterogeneous environmental gradients of temperature, salinity and food availability across a 30° latitudinal range (3,334 km), to test whether plasticity in calcareous shell production and composition, from juveniles to large adults, mediates geographical patterns of resilience to climate change in critical foundation species, the mussels Mytilus edulis and M. trossulus. We find shell calcification decreased towards high latitude, with mussels producing thinner shells with a higher organic content in polar than temperate regions. Salinity was the best predictor of within-region differences in mussel shell deposition, mineral and organic composition. In polar, subpolar, and Baltic low-salinity environments, mussels produced thin shells with a thicker external organic layer (periostracum), and an increased proportion of calcite (prismatic layer, as opposed to aragonite) and organic matrix, providing potentially higher resistance against dissolution in more corrosive waters. Conversely, in temperate, higher salinity regimes, thicker, more calcified shells with a higher aragonite (nacreous layer) proportion were deposited, which suggests enhanced protection under increased predation pressure. Interacting effects of salinity and food availability on mussel shell composition predict the deposition of a thicker periostracum and organic-enriched prismatic layer under forecasted future environmental conditions, suggesting a capacity for increased protection of high-latitude populations from ocean acidification. These findings support biomineralization plasticity as a potentially advantageous compensatory mechanism conferring Mytilus species a protective capacity for quantitative and qualitative trade-offs in shell deposition as a response to regional alterations of abiotic and biotic conditions in future environments. Our work illustrates that compensatory mechanisms, driving plastic responses to the spatial structure of multiple stressors, can define geographical patterns of unanticipated species resilience to global environmental change.

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves, Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early-breeding house sparrows and hybrids and rural late-breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z-chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land-use changes in a mosaic landscape.

The demographic history of Atlantic salmon (Salmo salar) across its distribution range reconstructed from approximate Bayesian computations

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo-genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.

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.

New diagnostic SNP molecular markers for the Mytilus species complex

The development of diagnostic markers has been a long-standing interest of population geneticists as it allows clarification of taxonomic uncertainties. Historically, there has been much debate on the taxonomic status of species belonging to the Mytilus species complex (M. edulis, M. galloprovincialis and M. trossulus), and whether they are discrete species. We analysed reference pure specimens of M. edulis, M. galloprovincialis and M. trossulus, using Restriction site associated DNA (RAD) sequencing and identified over 6,000 SNP markers separating the three species unambiguously. We developed a panel of diagnostic SNP markers for the genotyping of Mytilus species complex as well as the identification of hybrids and interspecies introgression events in Mytilus species. We validated a panel of twelve diagnostic SNP markers which can be used for species genotyping. Being able to accurately identify species and hybrids within the Mytilus species complex is important for the selective mussel stock management, the exclusion of invasive species, basic physiology and bio-diversity studies.

The demographic history of Atlantic salmon (Salmo salar) across its distribution range reconstructed from approximate Bayesian computations

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo-genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.

Combined genotyping, microbial diversity and metabolite profiling studies on farmed Mytilus spp. from Kiel Fjord

The blue mussel Mytilus is a popular food source with high economical value. Species of the M. edulis complex (M. edulis, M. galloprovincialis and M. trossulus) hybridise whenever their geographic ranges overlap posing difficulties to species discrimination, which is important for blue mussel aquaculture. The aim of this study was to determine the genetic structure of farmed blue mussels in Kiel Fjord. Microbial and metabolic profile patterns were studied to investigate a possible dependency on the genotype of the bivalves. Genotyping confirmed the complex genetic structure of the Baltic Sea hybrid zone and revealed an unexpected dominance of M. trossulus alleles being in contrast to the predominance of M. edulis alleles described for wild Baltic blue mussels. Culture-dependent and -independent microbial community analyses indicated the presence of a diverse Mytilus-associated microbiota, while an LC-MS/MS-based metabolome study identified 76 major compounds dominated by pigments, alkaloids and polyketides in the whole tissue extracts. Analysis of mussel microbiota and metabolome did not indicate genotypic dependence, but demonstrated high intraspecific variability of farmed mussel individuals. We hypothesise that individual differences in microbial and metabolite patterns may be caused by high individual plasticity and might be enhanced by e.g. nutritional condition, age and gender.

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.

Single nucleotide polymorphisms in native South American Atlantic coast populations of smooth shelled mussels: hybridization with invasive European Mytilus galloprovincialis

BACKGROUND

Throughout the world, harvesting of mussels Mytilus spp. is based on the exploitation of natural populations and aquaculture. Aquaculture activities include transfers of spat and live adult mussels between various geographic locations, which may result in large-scale changes in the world distribution of Mytilus taxa. Mytilus taxa are morphologically similar and difficult to distinguish. In spite of much research on taxonomy, evolution and geographic distribution, the native Mytilus taxa of the Southern Hemisphere are poorly understood. Recently, single nucleotide polymorphisms (SNPs) have been used to clarify the taxonomic status of populations of smooth shelled mussels from the Pacific coast of South America. In this paper, we used a set of SNPs to characterize, for the first time, populations of smooth shelled mussels Mytilus from the Atlantic coast of South America.

RESULTS

Mytilus spp. samples were collected from eastern South America. Six reference samples from the Northern Hemisphere were used: Mytilus edulis from USA and Northern Ireland, Mytilus trossulus from Canada, and Mytilus galloprovincialis from Spain and Italy. Two other reference samples from the Southern Hemisphere were included: M. galloprovincialis from New Zealand and Mytilus chilensis from Chile. Fifty-five SNPs were successfully genotyped, of which 51 were polymorphic. Population genetic analyses using the STRUCTURE program revealed the clustering of eight populations from Argentina (Mytilus platensis) and the clustering of the sample from Ushuaia with M. chilensis from Chile. All individuals in the Puerto Madryn (Argentina) sample were identified as M. platensis × M. galloprovincialis F2 (88.89%) hybrids, except one that was classified as Mediterranean M. galloprovincialis. No F1 hybrids were observed.

CONCLUSIONS

We demonstrate that M. platensis (or Mytilus edulis platensis) and M. chilensis are distinct native taxa in South America, which indicates that the evolutionary histories of Mytilus taxa along the Atlantic and Pacific coasts differ. M. platensis is endangered by hybridization with M. galloprovincialis that was introduced from Europe into the Puerto Madryn area in Argentina, presumably by accidental introduction via ship traffic. We confirm the occurrence of a native M. chilensis population in southern Argentina on the coast of Patagonia.

Blue mussel shell shape plasticity and natural environments: a quantitative approach

Shape variability represents an important direct response of organisms to selective environments. Here, we use a combination of geometric morphometrics and generalised additive mixed models (GAMMs) to identify spatial patterns of natural shell shape variation in the North Atlantic and Arctic blue mussels, Mytilus edulis and M. trossulus, with environmental gradients of temperature, salinity and food availability across 3980 km of coastlines. New statistical methods and multiple study systems at various geographical scales allowed the uncoupling of the developmental and genetic contributions to shell shape and made it possible to identify general relationships between blue mussel shape variation and environment that are independent of age and species influences. We find salinity had the strongest effect on the latitudinal patterns of Mytilus shape, producing shells that were more elongated, narrower and with more parallel dorsoventral margins at lower salinities. Temperature and food supply, however, were the main drivers of mussel shape heterogeneity. Our findings revealed similar shell shape responses in Mytilus to less favourable environmental conditions across the different geographical scales analysed. Our results show how shell shape plasticity represents a powerful indicator to understand the alterations of blue mussel communities in rapidly changing environments.

Blue mussel shell shape plasticity and natural environments: a quantitative approach

Shape variability represents an important direct response of organisms to selective environments. Here, we use a combination of geometric morphometrics and generalised additive mixed models (GAMMs) to identify spatial patterns of natural shell shape variation in the North Atlantic and Arctic blue mussels, Mytilus edulis and M. trossulus, with environmental gradients of temperature, salinity and food availability across 3980 km of coastlines. New statistical methods and multiple study systems at various geographical scales allowed the uncoupling of the developmental and genetic contributions to shell shape and made it possible to identify general relationships between blue mussel shape variation and environment that are independent of age and species influences. We find salinity had the strongest effect on the latitudinal patterns of Mytilus shape, producing shells that were more elongated, narrower and with more parallel dorsoventral margins at lower salinities. Temperature and food supply, however, were the main drivers of mussel shape heterogeneity. Our findings revealed similar shell shape responses in Mytilus to less favourable environmental conditions across the different geographical scales analysed. Our results show how shell shape plasticity represents a powerful indicator to understand the alterations of blue mussel communities in rapidly changing environments.

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.

Genome-wide comparisons reveal evidence for a species complex in the black-lip pearl oyster Pinctada margaritifera (Bivalvia: Pteriidae)

Evolutionary relationships in the black-lip pearl oyster Pinctada margaritifera which is highly valued for pearl production remain poorly understood. This species possesses an 18,000 km Indo-Pacific natural distribution, and its current description includes six subspecies defined exclusively on morphological characters. To evaluate its taxonomic identity using molecular data, 14 populations in both the Indian and Pacific Oceans (n = 69), and the congeneric taxa P. maxima and P. mazatlanica (n = 29 and n = 10, respectively) were sampled. Phylogenomic reconstruction was carried out using both 8,308 genome-wide SNPs and 10,000 dominant loci (DArTseq PAVs). Reconstructions using neighbour-joining (Nei's 1972 distances), maximum likelihood and Bayesian approaches all indicate that the taxonomy of P. margaritifera is quite complex, with distinct evolutionary significant units (ESUs) identified within Tanzanian and Iranian populations. Contrastingly, phylogenies generated for Pacific Ocean oysters resolved a large monophyletic clade, suggesting little support for two current morphological subspecies classifications. Furthermore, P. mazatlanica formed a basal clade closest to French Polynesian P. margaritifera, suggesting it may be conspecific. Collectively, these findings provide evidence that P. margaritifera comprises a species complex, perhaps as a result of population fragmentation and increased divergence at range limits.

Disruption of doubly uniparental inheritance of mitochondrial DNA associated with hybridization area of European Mytilus edulis and Mytilus trossulus in Norway

Doubly uniparental inheritance of mitochondria (DUI) is best known in the blue mussel Mytilus. Under this model, two types of mitochondrial DNA exist: female type (F), transmitted from females to offspring of both genders, and male type (M), transmitted exclusively from males to sons. The mitogenomes are usually highly divergent, but an occasional replacement of a typical M genome by a particular F genome has been postulated to explain reduction of this divergence. Disruption of the DUI model has been reported in hybridization areas. Here, we present a new case of DUI disruption in a hybrid M. trossulus/M. edulis population from the North Sea (Norway). No M haplotypes derived from M. trossulus were identified in this population. Typical M haplotypes derived from M. edulis (ME) were rare. Two F-type haplogroups were found: one derived from M. edulis (FE) and the second derived from M. trossulus (FT). Many haplotypes from the FT group were recombinants, with the male CR sequence coming from the M. trossulus genome (FT1 haplogroup) in contrast to M. edulis CR as in the Baltic. FT1 haplotypes were abundant in the studied population, including homoplasmic females. However, males significantly more often carried these haplotypes; therefore, male heteroplasmy involved the original FE and recombinant FT, indicating that the FT genome undergoes masculinization. Structural similarity of FT1 CR with previously reported, masculinized Baltic haplotypes, which were derived from FE/ME recombination, provides further evidence that CR M-F recombination is a prerequisite for masculinization, also in the context of native M. trossulus mtDNA.

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.

Epipodial Tentacle Gene Expression and Predetermined Resilience to Summer Mortality in the Commercially Important Greenlip Abalone, Haliotis laevigata

"Summer mortality" is a phenomenon that occurs during warm water temperature spikes that results in the mass mortality of many ecologically and economically important mollusks such as abalone. This study aimed to determine whether the baseline gene expression of abalone before a laboratory-induced summer mortality event was associated with resilience to summer mortality. Tentacle transcriptomes of 35 greenlip abalone (Haliotis laevigata) were sequenced prior to the animals being exposed to an increase in water temperature-simulating conditions which have previously resulted in summer mortality. Abalone derived from three source locations with different environmental conditions were categorized as susceptible or resistant to summer mortality depending on whether they died or survived after the water temperature was increased. We detected two genes showing significantly higher expression in resilient abalone relative to susceptible abalone prior to the laboratory-induced summer mortality event. One of these genes was annotated through the NCBI non-redundant protein database using BLASTX to an anemone (Exaiptasia pallida) Transposon Ty3-G Gag Pol polyprotein. Distinct gene expression signatures were also found between resilient and susceptible abalone depending on the population origin, which may suggest divergence in local adaptation mechanisms for resilience. Many of these genes have been suggested to be involved in antioxidant and immune-related functions. The identification of these genes and their functional roles have enhanced our understanding of processes that may contribute to summer mortality in abalone. Our study supports the hypothesis that prestress gene expression signatures are indicative of the likelihood of summer mortality.

Insights from the Shell Proteome: Biomineralization to Adaptation

Bivalves have evolved a range of complex shell forming mechanisms that are reflected by their incredible diversity in shell mineralogy and microstructures. A suite of proteins exported to the shell matrix space plays a significant role in controlling these features, in addition to underpinning some of the physical properties of the shell itself. Although, there is a general consensus that a minimum basic protein tool kit is required for shell construction, to date, this remains undefined. In this study, the shell matrix proteins (SMPs) of four highly divergent bivalves (The Pacific oyster, Crassostrea gigas; the blue mussel, Mytilus edulis; the clam, Mya truncata, and the king scallop, Pecten maximus) were analyzed in an identical fashion using proteomics pipeline. This enabled us to identify the critical elements of a "basic tool kit" for calcification processes, which were conserved across the taxa irrespective of the shell morphology and arrangement of the crystal surfaces. In addition, protein domains controlling the crystal layers specific to aragonite and calcite were also identified. Intriguingly, a significant number of the identified SMPs contained domains related to immune functions. These were often are unique to each species implying their involvement not only in immunity, but also environmental adaptation. This suggests that the SMPs are selectively exported in a complex mix to endow the shell with both mechanical protection and biochemical defense.

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.

The oceanic concordance of phylogeography and biogeography: a case study in Notochthamalus

Dispersal and adaptation are the two primary mechanisms that set the range distributions for a population or species. As such, understanding how these mechanisms interact in marine organisms in particular - with capacity for long-range dispersal and a poor understanding of what selective environments species are responding to - can provide useful insights for the exploration of biogeographic patterns. Previously, the barnacle Notochthamalus scabrosus has revealed two evolutionarily distinct lineages with a joint distribution that suggests an association with one of the two major biogeographic boundaries (~30°S) along the coast of Chile. However, spatial and genomic sampling of this system has been limited until now. We hypothesized that given the strong oceanographic and environmental shifts associated with the other major biogeographic boundary (~42°S) for Chilean coastal invertebrates, the southern mitochondrial lineage would dominate or go to fixation in locations further to the south. We also evaluated nuclear polymorphism data from 130 single nucleotide polymorphisms to evaluate the concordance of the signal from the nuclear genome with that of the mitochondrial sample. Through the application of standard population genetic approaches along with a Lagrangian ocean connectivity model, we describe the codistribution of these lineages through a simultaneous evaluation of coastal lineage frequencies, an approximation of larval behavior, and current-driven dispersal. Our results show that this pattern could not persist without the two lineages having distinct environmental optima. We suggest that a more thorough integration of larval dynamics, explicit dispersal models, and near-shore environmental analysis can explain much of the coastal biogeography of Chile.