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

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.

Estimating blue mussel (Mytilus edulis) connectivity and settlement capacity in mid-latitude fjord regions

The mussel industry faces challenges such as low and inconsistent levels of larvae settlement and poor-quality spat, leading to variable production. However, mussel farming remains a vital sustainable and environmentally responsible method for producing protein, fostering ecological responsibility in the aquaculture sector. We investigate the population connectivity and larval dispersion of blue mussels (Mytilus edulis) in Scottish waters, as a case study, using a multidisciplinary approach that combined genetic data and particle modelling. This research allows us to develop a thorough understanding of blue mussel population dynamics in mid-latitude fjord regions, to infer gene-flow patterns, and to estimate population divergence. Our findings reveal a primary south-to-north particle transport direction and the presence of five genetic clusters. We discover a significant and continuous genetic material exchange among populations within the study area, with our biophysical model's outcomes aligning with our genetic observations. Additionally, our model reveals a robust connection between the southwest coast and the rest of the west coast. This study will guide the preservation of mussel farming regions, ensuring sustainable populations that contribute to marine ecosystem health and resilience.

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.

Spawning acts as a metabolic stressor enhanced by hypoxia and independent of sex in a broadcast marine spawner

Broadcast spawners, like the blue mussel Mytilus edulis, experience substantial energy expenditure during spawning due to extensive gamete release that can divert energy from other functions. This energetic cost might be intensified by environmental stressors, including hypoxia that suppress aerobic metabolism. However, the energy implications of spawning in marine broadcast spawners have not been well studied. We examined the effects of short-term hypoxia (7 days) and spawning on mitochondrial activity, reactive oxygen species (ROS) production, and cellular energy allocation (ratio of tissue energy reserves to energy demand) in somatic tissues of M. edulis. Under normoxic conditions, post-spawning (72 h) recovery correlated with increased phosphorylation (OXPHOS) rate in mitochondria from the digestive gland, while hypoxia inhibited this response. Regardless of oxygen levels, mitochondrial ROS production decreased after spawning, indicating M. edulis' ability to prevent oxidative stress. Spawning led to reduced energy reserves in somatic tissues (the gills and the digestive gland), highlighting significant energy cost of spawning primarily fueled by lipid and protein breakdown. Additionally, cellular energy allocation dropped 3 h post-spawning, indicating a shift in energy demand and supply. Normoxic conditions allowed recovery in 72 h, but hypoxia hindered recuperation. These findings underscore spawning's bioenergetic challenge for broadcast spawners like M. edulis, potentially elevating post-spawning mortality risk, especially in hypoxic coastal habitats.

Filtration and respiration responses of mussels (Mytilus edulis) to trematode parasite infections (Renicola roscovita) and transient heat exposure

The mussel Mytilus edulis, a host to various trematode species, experiences performance decrements due to these infections. Yet, the impact magnitude and potential interactions with environmental stressors remain largely unexplored. This study scrutinizes the effect of Renicola roscovita infections on mussel filtration and respiration. We first assessed performance in both uninfected and lab-infected mussels at a mild temperature (16 °C), following an acute heat ramp to 30.5 °C and subsequent cooling. The experiment revealed neither a significant direct impact of the infection on the mussels' performance, nor any significant interplay between the infection and temperature variations. To account for possible infection effects obscured by low sample sizes or mussel size disparities, we conducted a reassessment at 16 °C using both small and large mussels. Infection notably hampered filtration in large mussels, with a marginal impact on smaller ones. A positive correlation was found between infection intensity and mussel filtration capacity, though the infection had no discernible impact on respiration. Our consistent finding of an 11-12 % infection effect size across all experiments indicates a slight reduction in mussel filtration due to trematode infections. While the exacerbating effect of transient heat stress on the infection's impact on filtration was not statistically significant, future investigations should explore potential interactions with prolonged heat stress. Our findings underscore the nuanced ways in which parasitic infections can influence marine bivalve physiology, emphasizing the need for more comprehensive studies that incorporate environmental stressors, such as heat stress, to fully elucidate the impact of parasitism on marine ecosystem health and resilience.

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.

Gone with sunscreens: Responses of blue mussels (Mytilus edulis) to a wide concentration range of a UV filter ensulizole

Organic UV filters have emerged as a new threat to marine organisms, but ecotoxicological studies have so far focused on only a few substances despite the chemical diversity of these synthetic sunscreen agents. Here we examined the responses of blue mussels Mytilus edulis to ensulizole, a non-lipophilic UV filter commonly found in the Baltic Sea. Mussels were exposed for three weeks to five ensulizole concentrations of 10, 10², 10³, 10⁴, and 10⁵ ng/L. Stress on stress response was evaluated by subjecting mussels to air exposure. A battery of biomarkers related to detoxification and antioxidant defense, oxidative stress damage, energy reserves and metabolism, autophagy, apoptosis, inflammation, and DNA damage was measured in the gills and the digestive gland. In general, ensulizole affected the antioxidant response, energy storage, and cell death-related processes in mussel tissues. Mussels exposed to low, environmentally relevant concentrations of ensulizole had a shorter air survival time than the control. Ensulizole often showed the non-monotonic concentration-response curves, suggesting the complex effects of this UV filter at molecular, biochemical, and organismal levels.

Biomarker-based assessment of sublethal toxicity of organic UV filters (ensulizole and octocrylene) in a sentinel marine bivalve Mytilus edulis

The global occurrence of organic UV filters in the marine environment is of increasing ecotoxicological concern. Here we assessed the toxicity of UV filters ensulizole and octocrylene in the blue mussels Mytilus edulis exposed to 10 or 100 μg l^-1 of octocrylene and ensulizole for two weeks. An integrated battery of biochemical and molecular biomarkers related to xenobiotics metabolism and cellular toxicity (including oxidative stress, DNA damage, apoptosis, autophagy and inflammation) was used to assess the toxicity of these UV filters in the mussels. Octocrylene (but not ensulizole) accumulated in the mussel tissues during the waterborne exposures. Both studied UV filters induced sublethal toxic effects in M. edulis at the investigated concentrations. These effects involved induction of oxidative stress, genotoxicity (indicated by upregulation of DNA damage sensing and repair markers), upregulation of apoptosis and inflammation, and dysregulation of the xenobiotic biotransformation system. Octocrylene induced cellular stress in a concentration-dependent manner, whereas ensulizole appeared to be more toxic at the lower (10 μg l^-1) studied concentration than at 100 μg l^-1. The different concentration-dependence of sublethal effects and distinct toxicological profiles of ensulizole and octocrylene show that the environmental toxicity is not directly related to lipophilicity and bioaccumulation potential of these UV filters and demonstrate the importance of using bioassays for toxicity assessment of emerging pollutants in coastal marine ecosystems.

Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas

In the Himalayas, a number of secondary contact zones have been described for vicariant vertebrate taxa. However, analyses of genetic divergence and admixture are missing for most of these examples. In this study, we provide a population genetic analysis for the coal tit (Periparus ater) hybrid zone in Nepal. Intermediate phenotypes between the distinctive western "spot-winged tit" (P. a. melanolophus) and Eastern Himalayan coal tits (P. a. aemodius) occur across a narrow range of <100 km in western Nepal. As a peculiarity, another distinctive cinnamon-bellied form is known from a single population so far. Genetic admixture of western and eastern mitochondrial lineages was restricted to the narrow zone of phenotypically intermediate populations. The cline width was estimated 46 km only with a center close to the population of the cinnamon-bellied phenotype. In contrast, allelic introgression of microsatellite loci was asymmetrical from eastern P. a. aemodius into far western populations of phenotypic P. a. melanolophus but not vice versa. Accordingly, the microsatellite cline was about 3.7 times wider than the mitochondrial one.

Salinity-dependent effects of ZnO nanoparticles on bioenergetics and intermediate metabolite homeostasis in a euryhaline marine bivalve, Mytilus edulis

Engineered nanoparticles including ZnO nanoparticles (nZnO) are important emerging pollutants in aquatic ecosystems creating potential risks to coastal ecosystems and associated biota. The toxicity of nanoparticles and its interaction with the important environmental stressors (such as salinity variation) are not well understood in coastal organisms and require further investigation. Here, we examined the interactive effects of 100 μg l^-1 nZnO or dissolved Zn (as a positive control for Zn²⁺ release) and salinity (normal 15, low 5, and fluctuating 5-15) on bioenergetics and intermediate metabolite homeostasis of a keystone marine bivalve, the blue mussel Mytilus edulis from the Baltic Sea. nZnO exposures did not lead to strong disturbances in energy or intermediate metabolite homeostasis regardless of the salinity regime. Dissolved Zn exposures suppressed the mitochondrial ATP synthesis capacity and coupling as well as anaerobic metabolism and modified the free amino acid profiles in the mussels indicating that dissolved Zn is metabolically more damaging than nZnO. The environmental salinity regime strongly affected metabolic homeostasis and altered physiological and biochemical responses to nZnO or dissolved Zn in the mussels. Exposure to low (5) or fluctuating (5-15) salinity affected the physiological condition, energy metabolism and homeostasis, as well as amino acid metabolism in M. edulis. Generally, fluctuating salinity (5-15) appeared bioenergetically less stressful than constantly hypoosmotic stress (salinity 5) in M. edulis indicating that even short (24 h) periods of recovery might be sufficient to restore the metabolic homeostasis in this euryhaline species. Notably, the biological effects of nZnO and dissolved Zn became progressively less detectable as the salinity stress increased. These findings demonstrate that habitat salinity must be considered in the biomarker-based assessment of the toxic effects of nanopollutants on coastal organisms.

Transcriptomic analysis of shell repair and biomineralization in the blue mussel, Mytilus edulis

BACKGROUND

Biomineralization by molluscs involves regulated deposition of calcium carbonate crystals within a protein framework to produce complex biocomposite structures. Effective biomineralization is a key trait for aquaculture, and animal resilience under future climate change. While many enzymes and structural proteins have been identified from the shell and in mantle tissue, understanding biomieralization is impeded by a lack of fundamental knowledge of the genes and pathways involved. In adult bivalves, shells are secreted by the mantle tissue during growth, maintenance and repair, with the repair process, in particular, amenable to experimental dissection at the transcriptomic level in individual animals.

RESULTS

Gene expression dynamics were explored in the adult blue mussel, Mytilus edulis, during experimentally induced shell repair, using the two valves of each animal as a matched treatment-control pair. Gene expression was assessed using high-resolution RNA-Seq against a de novo assembled database of functionally annotated transcripts. A large number of differentially expressed transcripts were identified in the repair process. Analysis focused on genes encoding proteins and domains identified in shell biology, using a new database of proteins and domains previously implicated in biomineralization in mussels and other molluscs. The genes implicated in repair included many otherwise novel transcripts that encoded proteins with domains found in other shell matrix proteins, as well as genes previously associated with primary shell formation in larvae. Genes with roles in intracellular signalling and maintenance of membrane resting potential were among the loci implicated in the repair process. While haemocytes have been proposed to be actively involved in repair, no evidence was found for this in the M. edulis data.

CONCLUSIONS

The shell repair experimental model and a newly developed shell protein domain database efficiently identified transcripts involved in M. edulis shell production. In particular, the matched pair analysis allowed factoring out of much of the inherent high level of variability between individual mussels. This snapshot of the damage repair process identified a large number of genes putatively involved in biomineralization from initial signalling, through calcium mobilization to shell construction, providing many novel transcripts for future in-depth functional analyses.

Exposure to dissolved TNT causes multilevel biological effects in Baltic mussels (Mytilus spp.)

Baltic mussels (Mytilus spp.) were exposed to the explosive trinitrotoluene (TNT) for 96 h (0.31-10.0 mg/L) and 21 d (0.31-2.5 mg/L). Bioaccumulation of TNT and its degradation products (2- and 4-ADNT) as well as biological effects ranging from the gene and cellular levels to behaviour were investigated. Although no mortality occurred in the concentration range tested, uptake and metabolism of TNT and responses in antioxidant enzymes and histochemical biomarkers were observed already at the lowest concentrations. The characteristic shell closure behaviour of bivalves at trigger concentrations led to complex exposure patterns and non-linear responses to the exposure concentrations. Conclusively, exposure to TNT exerts biomarker reponses in mussels already at 0.31 mg/L while effects are recorded also after a prolonged exposure although no mortality occurs. Finally, more attention should be paid on shell closure of bivalves in exposure studies since it plays a marked role in definining toxicity threshold levels.

Deciphering shell proteome within different Baltic populations of mytilid mussels illustrates important local variability and potential consequences in the context of changing marine conditions

Molluscs defend themselves against predation and environmental stressors through the possession of mineralized shells. Mussels are widely used to predict the effects of abiotic factors such as salinity and pH on marine calcifiers in the context of changing ocean conditions. Shell matrix proteins are part of the molecular control regulating the biomineralization processes underpinning shell production. Under changing environmental conditions, differential expression of these proteins leads to the phenotypic plasticity of shells seen in many mollusc species. Low salinity decreases the availability of calcium and inorganic carbon in seawater and consequently energetic constraints often lead to thin, small and fragile shells in Mytilid mussels inhabiting Baltic Sea. To understand how the modulation of shell matrix proteins alters biomineralization, we compared the shell proteomes of mussels living under full marine conditions in the North Sea to those living in the low saline Baltic Sea. Modulation of proteins comprising the Mytilus biomineralization tool kit is observed. These data showed a relative increase in chitin related proteins, decrease in SD-rich, GA-rich shell matrix proteins indicating that altered protein scaffolding and mineral nucleation lead to impaired shell microstructures influencing shell resistance in Baltic Mytilid mussels. Interestingly, proteins with immunity domains in the shell matrix are also found to be modulated. Shell traits such as periostracum thickness, organic content and fracture resistance qualitatively correlates with the modulation of SMPs in Mytilid mussels providing key insights into control of biomineralization at molecular level in the context of changing marine conditions.

Effects of first intermediate host density, host size and salinity on trematode infections in mussels of the south-western Baltic Sea

Trematode prevalence and abundance in hosts are known to be affected by biotic drivers as well as by abiotic drivers. In this study, we used the unique salinity gradient found in the south-western Baltic Sea to: (i) investigate patterns of trematode infections in the first intermediate host, the periwinkle Littorina littorea and in the downstream host, the mussel Mytilus edulis, along a regional salinity gradient (from 13 to 22) and (ii) evaluate the effects of first intermediate host (periwinkle) density, host size and salinity on trematode infections in mussels. Two species dominated the trematode community, Renicola roscovita and Himasthla elongata. Salinity, mussel size and density of infected periwinkles were significantly correlated with R. roscovita, and salinity and density correlated with H. elongata abundance. These results suggest that salinity, first intermediate host density and host size play an important role in determining infection levels in mussels, with salinity being the main major driver. Under expected global change scenarios, the predicted freshening of the Baltic Sea might lead to reduced trematode transmission, which may be further enhanced by a potential decrease in periwinkle density and mussel size.

Complementarity and discriminatory power of genotype and otolith shape in describing the fine-scale population structure of an exploited fish, the common sole of the Eastern English Channel

Marine organisms show population structure at a relatively fine spatial scale, even in open habitats. The tools commonly used to assess subtle patterns of connectivity have diverse levels of resolution and can complement each other to inform on population structure. We assessed and compared the discriminatory power of genetic markers and otolith shape to reveal the population structure on evolutionary and ecological time scales of the common sole (Solea solea), living in the Eastern English Channel (EEC) stock off France and the UK. First, we genotyped fish with Single Nucleotide Polymorphisms to assess population structure at an evolutionary scale. Then, we tested for spatial segregation of the subunits using otolith shape as an integrative tracer of life history. Finally, a supervised machine learning framework was applied to genotypes and otolith phenotypes to probabilistically assign adults to subunits and assess the discriminatory power of each approach. Low but significant genetic differentiation was found among subunits. Moreover, otolith shape appeared to vary spatially, suggesting spatial population structure at fine spatial scale. However, results of the supervised discriminant analyses failed to discriminate among subunits, especially for otolith shape. We suggest that the degree of population segregation may not be strong enough to allow for robust fish assignments. Finally, this study revealed a weak yet existing metapopulation structure of common sole at the fine spatial scale of the EEC based on genotypes and otolith shape, with one subunit being more isolated. Our study argues for the use of complementary tracers to investigate marine population structure.

The explosive trinitrotoluene (TNT) induces gene expression of carbonyl reductase in the blue mussel (Mytilus spp.): a new promising biomarker for sea dumped war relicts?

Millions of tons of all kind of munitions, including mines, bombs and torpedoes have been dumped after World War II in the marine environment and do now pose a new threat to the seas worldwide. Beside the acute risk of unwanted detonation, there is a chronic risk of contamination, because the metal vessels corrode and the toxic and carcinogenic explosives (trinitrotoluene (TNT) and metabolites) leak into the environment. While the mechanism of toxicity and carcinogenicity of TNT and its derivatives occurs through its capability of inducing oxidative stress in the target biota, we had the idea if TNT can induce the gene expression of carbonyl reductase in blue mussels. Carbonyl reductases are members of the short-chain dehydrogenase/reductase (SDR) superfamily. They metabolize xenobiotics bearing carbonyl functions, but also endogenous signal molecules such as steroid hormones, prostaglandins, biogenic amines, as well as sugar and lipid peroxidation derived reactive carbonyls, the latter providing a defence mechanism against oxidative stress and reactive oxygen species (ROS). Here, we identified and cloned the gene coding for carbonyl reductase from the blue mussel Mytilus spp. by a bioinformatics approach. In both laboratory and field studies, we could show that TNT induces a strong and concentration-dependent induction of gene expression of carbonyl reductase in the blue mussel. Carbonyl reductase may thus serve as a biomarker for TNT exposure on a molecular level which is useful to detect TNT contaminations in the environment and to perform a risk assessment both for the ecosphere and the human seafood consumer.

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.

Ocean winter warming induced starvation of predator and prey

Ocean warming impacts the fitness of marine ectothermic species, leading to poleward range shifts, re-shuffling of communities, and changes in ecosystem services. While the detrimental effects of summer heat waves have been widely studied, little is known about the impacts of winter warming on marine species in temperate regions. Many species benefit from low winter temperature-induced reductions in metabolism, as these permit conservation of energy reserves that are needed to support reproduction in spring. Here, we used a unique outdoor mesocosm system to expose a coastal predator-prey system, the sea star Asterias and the blue mussel Mytilus, to different winter warming scenarios under near-natural conditions. We found that the body condition of mussels decreased in a linear fashion with increasing temperature. Sea star growth also decreased with increasing temperature, which was a function of unaltered predation rates and decreased mussel body condition. Asterias relative digestive gland mass strongly declined over the studied temperature interval (ca twofold). This could have severe implications for reproductive capacity in the following spring, as digestive glands provide reserve compounds to maturing gonads. Thus, both predator and prey suffered from a mismatch of energy acquisition versus consumption in warmer winter scenarios, with pronounced consequences for food web energy transfer in future oceans.

Interactive effects of salinity variation and exposure to ZnO nanoparticles on the innate immune system of a sentinel marine bivalve, Mytilus edulis

ZnO nanoparticles (nZnO) are released into the coastal environment from multiple sources, yet their toxicity to marine organisms is not well understood. We investigated the interactive effects of salinity (normal 15, low 5, and fluctuating 5-15) and nZnO (100 μg l-1) on innate immunity of the blue mussels Mytilus edulis from a brackish area of the Baltic Sea. Exposure to ionic Zn (100 μg l-1) was used to test whether the toxic effects of nZnO can be attributed to the potential release of Zn2+. Functional parameters and the expression of key immune-related genes were investigated in the mussels exposed to nZnO or ionic Zn under different salinity regimes for 21 days. nZnO exposures elevated hemocyte mortality, suppressed adhesion, stimulated phagocytosis, and led to an apparent increase in lysosomal volume. At salinity 15, nZnO suppressed the mRNA expression of the Toll-like receptors TLRb and c, C-lectin, and the complement system component C3q indicating impaired ability for pathogen recognition. In contrast, the mRNA levels of an antimicrobial peptide defensin increased during nZnO exposure at salinity 15. At fluctuating salinity (5-15), nZnO exposure increased expression of multiple immune-related genes in hemocytes including the complement system components C1 and C3q, and the Toll-like receptors TLRa, b and c. Low salinity (5) had strong immunosuppressive effects on the functional and molecular immune traits of M. edulis that overshadowed the effects of nZnO. The salinity-dependent modulation of immune response to nZnO cannot be attributed to the differences in the aggregation or solubility of nZnO, and likely reflects the interaction of the toxic effects of nanoparticles and physiological effects of the osmotic stress. These findings have implications for the environmental risk assessment of nanomaterials and the development of the context-specific biomarker baselines for coastal pollution monitoring.

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.

The effects of ZnO nanostructures of different morphology on bioenergetics and stress response biomarkers of the blue mussels Mytilus edulis

Biofouling causes massive economical losses in the maritime sector creating an urgent need for effective and ecologically non-harmful antifouling materials. Zinc oxide (ZnO) nanorod coatings show promise as an antifouling material; however, the toxicity of ZnO nanorods to marine organisms is not known. We compared the toxicity of suspended ZnO nanorods (NR) with that of ZnO nanoparticles (NP) and ionic Zn²⁺ in a marine bivalve Mytilus edulis exposed for two weeks to 10 or 100 μg Zn L^-1 of ZnO NPs, NRs or Zn²⁺, or to immobilized NRs. The multi-biomarker assessment included bioenergetics markers (tissue energy reserves, activity of mitochondrial electron transport system and autophagic enzymes), expression of apoptotic and inflammatory genes, and general stress biomarkers (oxidative lesions, lysosomal membrane stability and metallothionein expression). Exposure to ZnO NPs, NRs and Zn²⁺ caused accumulation of oxidative lesions in proteins and lipids, stimulated autophagy, and led to lysosomal membrane destabilization indicating toxicity. However, these responses were not specific for the form of Zn (NPs, NR or Zn²⁺) and showed no monotonous increase with increasing Zn concentrations in the experimental exposures. No major disturbance of the energy status was found in the mussels exposed to ZnO NPs, NRs, or Zn²⁺. Exposure to ZnO NPs and NRs led to a strong induction of apoptosis- and inflammation-related genes, which was not seen in Zn²⁺ exposures. Based on the integrated biomarker response, the overall toxicity as well as the pro-apoptotic and pro-inflammatory action was stronger in ZnO NPs compared with the NRs. Given the stability of ZnO NR coatings and the relatively low toxicity of suspended ZnO NR, ZnO NR coating might be considered a promising low-toxicity material for antifouling paints.

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.

In vivo characterization of bivalve larval shells: a confocal Raman microscopy study

In vivo confocal Raman microscopy (CRM), polarized light microscopy and Fourier transform infrared spectroscopy (FTIR) were used to determine if a significant amount of amorphous calcium carbonate (ACC) exists within larval shells of Baltic mytilid mussels (Mytilus edulis-like) and whether the amount of ACC varies during larval development. No evidence for ACC was found from the onset of shell deposition at 21 h post-fertilization (hpf) until 48 hpf. Larval Mytilus shells were crystalline from 21 hpf onwards and exhibited CRM and FTIR peaks characteristic of aragonite. Prior to shell deposition at 21 hpf, no evidence for carbonates was observed through in vivo CRM. We further analysed the composition of larval shells in three other bivalve species, Mercenaria mercenaria, Crassostrea gigas and Crassostrea virginica and observed no evidence for ACC, which is in contrast to previous work on the same species. Our findings indicate that larval bivalve shells are composed of crystalline aragonite and we demonstrate that conflicting results are related to sub-optimal measurements and misinterpretation of CRM spectra. Our results demonstrate that the common perception that ACC generally occurs as a stable and abundant precursor during larval bivalve calcification needs to be critically reviewed.

Global genetic diversity, lineage distribution, and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well-known example of a worldwide-distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.

Expression of calcification-related ion transporters during blue mussel larval development

The physiological processes driving the rapid rates of calcification in larval bivalves are poorly understood. Here, we use a calcification substrate-limited approach (low dissolved inorganic carbon, C T) and mRNA sequencing to identify proteins involved in bicarbonate acquisition during shell formation. As a secondary approach, we examined expression of ion transport and shell matrix proteins (SMPs) over the course of larval development and shell formation. We reared four families of Mytilus edulis under ambient (ca. 1865 µmol/kg) and low C T (ca. 941 µmol/kg) conditions and compared expression patterns at six developmental time points. Larvae reared under low C T exhibited a developmental delay, and a small subset of contigs was differentially regulated between ambient and low C T conditions. Of particular note was the identification of one contig encoding an anion transporter (SLC26) which was strongly upregulated (2.3-2.9 fold) under low C T conditions. By analyzing gene expression profiles over the course of larval development, we are able to isolate sequences encoding ion transport and SMPs to enhance our understanding of cellular pathways underlying larval calcification processes. In particular, we observe the differential expression of contigs encoding SLC4 family members (sodium bicarbonate cotransporters, anion exchangers), calcium-transporting ATPases, sodium/calcium exchangers, and SMPs such as nacrein, tyrosinase, and transcripts related to chitin production. With a range of candidate genes, this work identifies ion transport pathways in bivalve larvae and by applying comparative genomics to investigate temporal expression patterns, provides a foundation for further studies to functionally characterize the proteins involved in larval calcification.

Effects of a common pharmaceutical, atorvastatin, on energy metabolism and detoxification mechanisms of a marine bivalve Mytilus edulis

Biologically active compounds from pharmaceuticals cause concern due to their common occurrence in water and sediments of urbanized coasts and potential threat to marine organisms. Atorvastatin (ATO), a globally prescribed drug, is environmentally stable and bioavailable to marine organisms; however, the physiological and toxic effects of this drug on ecologically important coastal species are yet to be elucidated. We studied the effect of ATO (˜1.2 μg L^-1) on bioenergetics (including whole-organism and mitochondrial respiration, as well as tissue energy reserves and mRNA expression of genes involved in mitochondrial biogenesis and fatty acid metabolism in the gills and the digestive gland) of a keystone bivalve Mytulis edulis (the blue mussel) from the Baltic Sea. Xenobiotic detoxification systems including activity and mRNA expression of P-glycoprotein, and Phase I and II biotransformation enzymes (cytochrome P450 monooxygenase CYP1A and glutathione transferase, GST) were also assessed in the gill and digestive gland of the mussels. Exposure to ATO caused rapid uptake and biotransformation of the drug by the mussels. Standard metabolic rate of ATO-exposed mussels increased by 56% indicating higher maintenance costs, yet no changes were detected in the respiratory capacity of isolated mitochondria. ATO exposure led to ˜60% decrease in the lysosomal membrane stability of hemocytes and ˜3-fold decrease in the whole-organism P-glycoprotein-driven and diffusional efflux of xenobiotics indicating altered membrane properties. The digestive gland was a major target of ATO toxicity in the mussels. Exposure of mussels to ATO led to depletion of lipid, carbohydrate and protein pools, and suppressed transcription of key enzymes involved in mitochondrial biogenesis (peroxisome proliferator-activated receptor gamma coactivator 1-alpha PGC-1α) and fatty acid metabolism (acetyl-CoA carboxylase and CYP4Y1) in the digestive gland. No bioenergetic disturbances were observed in the gills of ATO-exposed mussels, and elevated GST activity indicated enhanced ATO detoxification in this tissue. These data demonstrate that ATO can act as a metabolic disruptor and chemosensitizer in keystone marine bivalves and warrant further investigations of statins as emerging pollutants of concern in coastal marine ecosystems.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

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.

Multi biomarker analysis of pollution effect on resident populations of blue mussels from the Baltic Sea

Anthropogenic pollution including metals, petroleum, toxins, nutrients and many others is a growing problem in the marine environment. These are important factors altering the environment and by that the fate of many local populations of marine organisms. The aim of this study was to assess the impact of selected point pollution sources on resident populations of the blue mussel (Mytilus edulis trossulus) in the Baltic Sea using multiple biomarker approach. The study used a nested sampling scheme in which sites from reference (REF) habitats are geographically paired with selected sites from sewage treatment plants (STP) and harbors (HAR). The results showed that mussels from harbors had a higher frequency of histological abnormalities in the digestive gland compared to mussels from sewage effluent affected areas and reference sites. However these mussels together with mussels from STPs had higher lipid content, body mass index (BMI) and gonado-somatic index (GSI) compared to mussels from reference sites. A marked spatial variability was found with a stronger toxicity of ambient environment affecting resident mussel populations in the Gulf of Gdańsk area, while an opposite pattern was found in Tvärminne area. Yet the blue mussels sampled in the Gulf of Gdańsk were characterized by the highest GSI and BMI values compared to Askö and Tvärminne populations. No differences in analyzed biomarker response related to species identity, measured by a species-specific genetic marker, were found indicative of strong genetic introgression in the Baltic Proper.

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.

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.