High Calcification Costs Limit Mussel Growth at Low Salinity

In coastal temperate regions such as the Baltic Sea, calcifying bivalves dominate benthic communities playing a vital ecological role in maintaining biodiversity and nutrient recycling. At low salinities, bivalves exhibit reduced growth and calcification rates which is thought to result from physiological constraints associated with osmotic stress. Calcification demands a considerable amount of energy in calcifying molluscs and estuarine habitats provide sub-optimal conditions for calcification due to low concentrations of calcification substrates and large variations in carbonate chemistry. Therefore, we hypothesize that slow growth rates in estuarine bivalves result from increased costs of calcification, rather than costs associated with osmotic stress. To investigate this, we estimated the cost of calcification for the first time in benthic bivalve life stages and the relative energy allocation to calcification in three Mytilus populations along the Baltic salinity gradient. Our results indicate that calcification rates are significantly reduced only in 6 psu populations compared to 11 and 16 psu populations, coinciding with ca. 2–3-fold higher calcification costs at low salinity and temperature. This suggests that reduced growth of Baltic Mytilus at low salinities results from increased calcification costs rather than osmotic stress related costs. We also reveal that shell growth (both calcification and shell organic production) demands 31–60% of available assimilated energy from food, which is significantly higher than previous estimates. Energetically expensive calcification represents a major constraint on growth of mytilids in the estuarine and coastal seas where warming, acidification and desalination are predicted over the next century

Salinity driven selection and local adaptation in Baltic Sea Mytilid mussels

Baltic blue mussels can colonise and dominate habitats with far lower salinity (<10 psu) than other Mytilus congeners. Pervasive gene flow was observed between Western Baltic Mytilus edulis living at high salinity conditions and Eastern Baltic M. trossulus living at lower salinites, with highest admixture proportions within a genetic transition zone located at intermediate salinities (Darss Sill area). Yet, we do not understand the impacts of low salinity on larval performance, and how salinity may act as an early selective pressure during passive larval drift across salinity gradients. This study tested whether larvae originating from two different populations along the natural salinity cline in the Baltic Sea have highest fitness at their native salinities. Our results suggest that Eastern Baltic M. trossulus (Usedom, 7 psu) and Western Baltic M. edulis (Kiel, 16 psu) larvae display better performance (fitness components: growth, mortality, settlement success) when reared at their respective native salinities. This suggests that these populations are adapted to their local environment. Additionally, species diagnostic markers were used for genetic analyses of transition zone (Ahrenshoop, 11 psu) mussel larvae exposed to low salinity. This revealed that low salinity selection resulted in a shift towards allele frequencies more typical for Eastern Baltic M. trossulus. Thus, salinity acts as a selective pressure during the pre-settlement phase and can shape the genetic composition of Baltic mussel populations driving local adaptation to low salinity. Future climate change driven desalination, therefore, has the potential to shift the Baltic Sea hybrid gradient westward with consequences for benthic ecosystem structure

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

Abstract 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

Genome-wide association and genomic prediction of growth traits in the European flat oyster (Ostrea edulis)

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes