Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis)

Population Genetic Divergence among Worldwide Gene Pools of the Mediterranean Mussel Mytilus galloprovincialis

The Mediterranean mussel Mytilus galloprovincialis is distributed in both hemispheres either natively or introduced. The updated population genetic distribution of this species provides a useful knowledge against which future distribution shifts could be assessed. This study, performed with seven microsatellite markers and three reference species (M. edulis, M. chilensis and M. trossulus), aimed to determine the scenario of genetic divergence between 15 samples of M. galloprovincialis from 10 localities in Europe, Africa, Asia, Australia, North America and South America. In agreement with previous data, M. trossulus was the most divergent taxon of the genus, but M. chilensis appeared as an intermediate taxon between M. edulis and M. galloprovincialis, though closer to this latter. M. galloprovincialis from the Atlantic Northeast appears as the most likely source of worldwide exotic settlements instead of the previously thought Mediterranean population. The successful worldwide establishment of M. galloprovincialis suggests it is a flexible evolutionary species (FES), i.e., a species or population whose genetic background allows it to rapidly adapt to changing environments. This natural endowed plastic adaptation makes it a candidate resilient species amidst the ongoing climatic change.

Limited population genetic variation but pronounced seascape genetic structuring in populations of the Mediterranean mussel (Mytilus galloprovincialis) from the eastern Adriatic Sea

Population genetic analysis of variation at five neutral microsatellite loci for Mediterranean mussels (Mytilus galloprovincialis) from 18 sites along the eastern Adriatic Sea revealed little or no spatial variation. In contrast, seascape genetics analysis revealed a pronounced locus-specific gradient in allelic and genotypic frequencies across the study region. At a sixth locus, MGE7, the frequencies of two alleles, MGE7²⁴³ and MGE7²⁴⁹, were strongly associated, negatively and positively, respectively, with a single environmental variable - minimum salinity (minSAL). The frequency of the MGE7^243/243 homozygous genotype was strongly negatively associated with minSAL, whereas the frequencies of the MGE7^246/249 and the MGE7^249/249 genotypes were strongly positively correlated with minSAL. Interpretation of these pronounced gradients is confounded by the fact that minSAL and another environmental variable, maximum sea surface temperature (maxSST), are highly correlated (R = -.911) and are therefore not necessarily acting independently. BLAST searches of the MGE7 locus against M. galloprovincialis whole genome shotgun sequence returned an alignment with contig mg10_S01094 (accession UYJE01010330.1) and 7 predicted M. galloprovincialis proteins VDI82194.1 - VDI82200.1. Conserved domain searches revealed a similar structure to the transcriptional regulator Msx2-interacting protein. The BLASTp search also returned significant alignments to Msx2-interacting proteins in Mytilus coruscus, Crassostrea virginica, and Haliotis rubra. The existence of the MGE7 gradient highlights the role that environmental variation may play in retarding gene flow among wild M. galloprovincialis populations, and also how the success of collection of young mussels (spat) from one site and their transfer to another site (the farm) may be influenced by a single factor such as minSAL or maxSST on a localized scale.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

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.

Convergent Evolution and Structural Adaptation to the Deep Ocean in the Protein-Folding Chaperonin CCTα

The deep ocean is the largest biome on Earth and yet it is among the least studied environments of our planet. Life at great depths requires several specific adaptations; however, their molecular mechanisms remain understudied. We examined patterns of positive selection in 416 genes from four brittle star (Ophiuroidea) families displaying replicated events of deep-sea colonization (288 individuals from 216 species). We found consistent signatures of molecular convergence in functions related to protein biogenesis, including protein folding and translation. Five genes were recurrently positively selected, including chaperonin-containing TCP-1 subunit α (CCTα), which is essential for protein folding. Molecular convergence was detected at the functional and gene levels but not at the amino-acid level. Pressure-adapted proteins are expected to display higher stability to counteract the effects of denaturation. We thus examined in silico local protein stability of CCTα across the ophiuroid tree of life (967 individuals from 725 species) in a phylogenetically corrected context and found that deep-sea-adapted proteins display higher stability within and next to the substrate-binding region, which was confirmed by in silico global protein stability analyses. This suggests that CCTα displays not only structural but also functional adaptations to deep-water conditions. The CCT complex is involved in the folding of ∼10% of newly synthesized proteins and has previously been categorized as a "cold-shock" protein in numerous eukaryotes. We thus propose that adaptation mechanisms to cold and deep-sea environments may be linked and highlight that efficient protein biogenesis, including protein folding and translation, is a key metabolic deep-sea adaptation.

High-throughput quantification of protein structural change reveals potential mechanisms of temperature adaptation in Mytilus mussels

Background

Temperature exerts a strong influence on protein evolution: species living in thermally distinct environments often exhibit adaptive differences in protein structure and function. However, previous research on protein temperature adaptation has focused on small numbers of proteins and on proteins adapted to extreme temperatures. Consequently, less is known about the types and quantity of evolutionary change that occurs to proteins when organisms adapt to small shifts in environmental temperature. In this study, these uncertainties were addressed by developing software that enabled comparison of structural changes associated with temperature adaptation (hydrogen bonding, salt bridge formation, and amino acid use) among large numbers of proteins from warm- and cold-adapted species of marine mussels, Mytilus galloprovincialis and Mytilus trossulus, respectively.

Results

Small differences in habitat temperature that characterize the evolutionary history of Mytilus mussels were sufficient to cause protein structural changes consistent with temperature adaptation. Hydrogen bonds and salt bridges that increase stability and protect against heat-induced denaturation were more abundant in proteins from warm-adapted M. galloprovincialis compared with proteins from cold-adapted M. trossulus. These structural changes were related to deviations in the use of polar and charged amino acids that facilitate formation of hydrogen bonds and salt bridges within proteins, respectively. Enzymes, in particular those within antioxidant and cell death pathways, were over-represented among proteins with the most hydrogen bonds and salt bridges in warm-adapted M. galloprovincialis. Unlike extremophile proteins, temperature adaptation in Mytilus proteins did not involve substantial changes in the number of hydrophobic or large volume amino acids, nor in the content of glycine or proline.

Conclusions

Small shifts in organism temperature tolerance, such as that needed to cope with climate warming, may result from structural and functional changes to a small percentage of the proteome. Proteins in which function is dependent on large conformational change, notably enzymes, may be particularly sensitive to temperature perturbation and represent foci for natural selection. Protein temperature adaptation can occur through different types and frequencies of structural change, and adaptive mechanisms used to cope with small shifts in habitat temperature appear different from mechanisms used to retain protein function at temperature extremes.

Physiological responses to heat stress in an invasive mussel Mytilus galloprovincialis depend on tidal habitat

Mussels are ecologically important organisms that can survive in subtidal and intertidal zones where they experience thermal stress. We know little about how mussels from different tidal habitats respond to thermal stress. We used the mussel Mytilus galloprovincialis from separate subtidal and intertidal populations to test whether heart rate and indicators of potential aerobic (citrate synthase activity) and anaerobic (cytosolic malate dehydrogenase activity) metabolic capacity are affected by increased temperatures while exposed to air or submerged in water. Subtidal mussels were affected by warming when submerged in water (decreased heart rate) but showed no effect in air. In contrast, intertidal mussels were affected by exposure to air (increased anaerobic capacity) but not by warming. Overall, physiological responses of mussels to thermal stress were dependent on their tidal habitat. These results highlight the importance of considering the natural habitat of mussels when assessing their responses to environmental challenges.

Cryptic diversity in smooth-shelled mussels on Southern Ocean islands: connectivity, hybridisation and a marine invasion

Background

Large numbers of endemic species inhabit subantarctic continental coasts and islands that are characterised by highly variable environmental conditions. Southern hemisphere populations of taxa that are morphologically similar to northern counterparts have traditionally been considered to be extensions of such Northern hemisphere taxa, and may not exhibit differentiation amongst geographically isolated populations in the Southern Ocean. Smooth-shelled blue mussels of the genus Mytilus that exhibit an anti-tropical distribution are a model group to study phylogeography, speciation and hybridisation in the sea, and contribute to the theory and practice of marine biosecurity.

Methods

We used a single nucleotide polymorphism (SNPs) panel that has the ability to accurately identify reference Northern and Southern hemisphere Mytilus taxa to test for evolutionary differentiation amongst native Southern Ocean island populations.

Results

Native mussels from the Falkland Islands and the Kerguelen Islands exhibited greatest affinity to native M. platensis d’Orbigny 1846 from the Atlantic coast of South America. The major Southern Ocean current flow from west to east is likely to explain the spreading of M. platensis to remote offshore islands, as adults via the process of rafting or perhaps directly as larvae. SNPs variation revealed that mussels from Tasmania were native and clearly differentiated from all other blue mussel groups in the Southern and Northern hemispheres. The native mussels M. planulatus from Tasmania and from mainland New Zealand (NZ), and tentatively M. aoteanus from the two NZ Southern Ocean offshore island groups (the Auckland Islands and Campbell Island), formed a distinct M. galloprovincialis–like Southern hemisphere group with closest affinity to Northern hemisphere M. galloprovincialis from the Mediterranean Sea. In all cases, the SNPs revealed evidence of hybridisation between two or more distinct taxa. The invasive Northern hemisphere M. galloprovincialis was identified only in Tasmania, amongst native mussels of a distinct Australian M. planulatus lineage.

Conclusion

Overall, our results reveal that Southern hemisphere island mussels have mixed genome ancestry and are native, not introduced by human activities. The preservation of distinct evolutionary lineages of Southern hemisphere species needs to be an ongoing focus of conservation efforts, given that population sizes on some of the remote offshore oceanic islands will be small and may be more easily adversely affected by invasion and subsequent hybridisation and introgression than larger populations elsewhere.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0332-y) contains supplementary material, which is available to authorized users.

Genomics and immunity of the Mediterranean mussel Mytilus galloprovincialis in a changing environment

The Mediterranean mussel (Mytilus galloprovincialis) is a marine invasive species cultured all over the world. Mussels are an appreciated resource in local aquaculture enterprises because of their robust production and resilience that translates into a reliable economic value. So far, no massive mortalities have been reported in natural or cultured populations of this species. In the last years, the knowledge about its immune system has greatly improved but there are still many questions to be answered. One of them is why mussels, with their high filtering activity, are able to be exposed to a high number of potential pathogens without getting infected and without developing an elevated inflammatory response. The sequencing of the mussel genome has revealed a very complex organization with high heterozygosity, abundance of repetitive sequences and extreme intraspecific sequence diversity among individuals, mainly in immune related genes. Among those genes, antimicrobial peptides are the most expressed gene families in mussels, highly polymorphic and with antimicrobial effect against molluscs pathogens, but also against pathogens of lower vertebrates and humans. The combination of a complex genome with the adaptation of mussel immune system to a changing environment could explain this high variability, not only in immune-related genes, but also in the functional response among individuals sampled in the same location and date.

RNA-seq coupled to proteomic analysis reveals high sperm proteome variation between two closely related marine mussel species

Speciation mechanisms in marine organisms have attracted great interest because of the apparent lack of substantial barriers to genetic exchange in marine ecosystems. Marine mussels of the Mytilus edulis species complex provide a good model to study mechanisms underlying species formation. They hybridise extensively at many localities and both pre- and postzygotic isolating mechanisms may be operating. Mussels have external fertilisation and sperm cells should show specific adaptations for survival and successful fertilisation. Sperm thus represent key targets in investigations of the molecular mechanisms underlying reproductive isolation. We undertook a deep transcriptome sequencing (RNA-seq) of mature male gonads and a 2DE/MS-based proteome analysis of sperm from Mytilus edulis and M. galloprovincialis raised in a common environment. We provide evidence of extensive expression differences between the two mussel species, and general agreement between the transcriptomic and proteomic results in the direction of expression differences between species. Differential expression is marked for mitochondrial genes and for those involved in spermatogenesis, sperm motility, sperm-egg interactions, the acrosome reaction, sperm capacitation, ATP reserves and ROS production. Proteins and their corresponding genes might thus be good targets in further genomic analysis of reproductive barriers between these closely related species. SIGNIFICANCE: Model systems for the study of fertilization include marine invertebrates with external fertilisation, such as abalones, sea urchins and mussels, because of the ease with which large quantities of gametes released into seawater can be collected after induced spawning. Unlike abalones and sea urchins, hybridisation has been reported between mussels of different Mytilus spp., which thus makes them very appealing for the study of reproductive isolation at both pre- and postzygotic levels. There is a lack of empirical proteomic studies on sperm samples comparing different Mytilus species, which could help to advance this study. A comparative analysis of sperm proteomes across different taxa may provide important insights into the fundamental molecular processes and mechanisms involved in reproductive isolation. It might also contribute to a better understanding of sperm function and of the adaptive evolution of sperm proteins in different taxa. There is now growing evidence from genomics studies that multiple protein complexes and many individual proteins might have important functions in sperm biology and the fertilisation process. From an applied perspective, the identification of sperm-specific proteins could also contribute to the improved understanding of fertility problems and as targets for fertility control.