The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks.

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

Expansion and diversity of caspases in Mytilus coruscus contribute to larval metamorphosis and environmental adaptation.

BACKGROUND: Apoptosis is involved (directly and indirectly) in several physiological processes including tissue remodeling during the development, the turnover of immune cells, and a defense against harmful stimuli. The disordered apoptotic process participates in the pathogenesis of various diseases, such as neoplasms, and chronic inflammatory or systemic autoimmune diseases, which are associated with its inadequate regulation. Caspases are vital components of the apoptotic pathway that are involved in developmental and immune processes. However, genome-wide identification and functional analysis of caspase have not been conducted in Mytilus coruscus, which is an economically important bivalve. RESULTS: Here, 47 caspase genes were identified from the genomes of M. coruscus, and the expansion of caspase-2/9 and caspase-3/6/7 genes were observed. Tandem duplication acts as an essential driver of gene expansion. The expanded caspase genes were highly diverse in terms of sequence, domain structure, and spatiotemporal expression profiles, suggesting their functional differentiation. The high expression of the expanded caspase genes at the pediveliger larvae stage and the result of apoptosis location in the velum suggest that the apoptosis mediated by them plays a critical role in the metamorphosis of M. coruscus larvae. In gill, caspase genes respond differently to the challenge of different strains, and most caspase-2/9 and caspase-3/6/7 genes were induced by copper stress, whereas caspase-8/10 genes were suppressed. Additionally, most caspase genes were upregulated in the mantle under ocean acidification which could weaken the biomineralization capacity of the mantle tissue. CONCLUSIONS: These results provide a comprehensive overview of the evolution and function of the caspase family and enhanced the understanding of the biological function of caspases in M. coruscus larval development and response to biotic and abiotic challenges.

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

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

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.

Transgressive gene expression and expression plasticity under thermal stress in a stable hybrid zone.

Interspecific hybridization can lead to myriad outcomes, including transgressive phenotypes in which the hybrids are more fit than either parent species. Such hybrids may display important traits in the context of climate change, able to respond to novel environmental conditions not previously experienced by the parent populations. While this has been evaluated in an agricultural context, the role of transgressive hybrids under changing conditions in the wild remains largely unexplored; this is especially true regarding transgressive gene expression. Using the blue mussel species complex (genus Mytilus) as a model system, we investigated the effects of hybridization on temperature induced gene expression plasticity by comparing expression profiles in parental species and their hybrids following a 2-week thermal challenge. Hybrid expression plasticity was most often like one parent or the other (50%). However, a large fraction of genes (26%) showed transgressive expression plasticity (i.e. the change in gene expression was either greater or lesser than that of both parent species), while only 2% were intermediately plastic in hybrids. Despite their close phylogenetic relationship, there was limited overlap in the differentially expressed genes responding to temperature, indicating interspecific differences in the responses to high temperature in which responses from hybrids are distinct from both parent species. We also identified differentially expressed long non-coding RNAs (lncRNAs), which we suggest may contribute to species-specific differences in thermal tolerance. Our findings provide important insight into the impact of hybridization on gene expression under warming. We propose transgressive hybrids may play an important role in population persistence under future warming conditions.

Transcriptome analysis reveals tissue-specific responses of Mytilus unguiculatus to Vibrio alginolyticus infection.

Mytilus unguiculatus is an important economic bivalve species with wide consumption and aquaculture value. Disease is one of the primary limiting factors in mussel aquaculture, thus understanding the response of different tissues of M. unguiculatus to pathogens is crucial for disease prevention and control. In this study, we investigated the physiological and transcriptomic responses of the gills, adductor muscle, and mantle of M. unguiculatus infected with Vibrio alginolyticus. The results showed that V. alginolyticus infection caused inflammation and tissue structure changes in the gill, adductor muscle and mantle of M. unguiculatus. Meanwhile, the activities of superoxide dismutase and catalase in the three tissues increased, while the total antioxidant capacity decreased, suggesting that M. unguiculatus have an activated defense mechanism against infection-induced oxidative stress, despite a compromised total antioxidant capacity. Transcriptomic studies reveal that infected M. unguiculatus exhibits upregulation of endocytosis, lysosome activity, cellular apoptosis, and immune-related signaling pathways, indicating that M. unguiculatus responds to pathogen invasion by upregulating efferocytosis. Compared with the gill and adductor muscle, the mantle had a higher level of mytimycin, mytilin and myticin, and the three tissues also increased the expression of mytimycin to cope with the invasion of pathogens. In addition, the analysis of genes related to taste transduction pathways and muscle contraction and relaxation found that after infection with V. alginolyticus, M. unguiculatus may reduce appetite by inhibiting taste transduction in the gill, while improving muscle contraction of the adductor muscle and keeping the shell closed, to resist further invasion of pathogens and reduce the risk of pathogen transmission in the population.

Histone derived antimicrobial peptides identified from Mytilus coruscus serum by peptidomics.

Histones and their N-terminal or C-terminal derived peptides have been studied in vertebrates and presented as potential antimicrobial agents playing important roles in the innate immune defenses. Although histones and their derived peptides had been reported as components of innate immunity in invertebrates, the knowledge about the histone derived antimicrobial peptides (HDAPs) in invertebrates are still limited. Using a peptidomic technique, a set of peptide fragments derived from the histones was identified in this study from the serum of microbes challenged Mytilus coruscus. Among the 85 identified histone-derived-peptides with high confidence, 5 HDAPs were chemically synthesized and the antimicrobial activities were verified, showing strong growth inhibition against Gram-positive bacteria, Gram-negative bacteria, and fungus. The gene expression level of the precursor histones matched by representative HDAPs were further tested using q-PCR, and the results showed a significant upregulation of the histone gene expression levels in hemocytes, gill, and mantle of the mussel after immune stress. In addition, three identified HDAPs were selected for preparation of specific antibodies, and the corresponding histones and their derived C-terminal fragments were detected by Western blotting in the blood cell and serum of immune challenged mussel, respectively, indicating the existence of HDAPs in M. coruscus. Our findings revealed the immune function of histones in Mytilus, and confirmed the existence of HDAPs in the mussel. The identified Mytilus HDAPs represent a new source of immune effector with antimicrobial function in the innate immune system, and thus provide promising candidates for the treatment of microbial infections in aquaculture and medicine.

Identification and expression analysis of two steamer-like retrotransposons in the Chilean blue mussel (Mytilus chilensis).

BACKGROUND: Disseminated neoplasia (DN) is a proliferative cell disorder of the circulatory system of bivalve mollusks. The disease is transmitted between individuals and can also be induced by external chemical agents such as bromodeoxyuridine. In Mya arenaria, we have cloned and characterized an LTR-retrotransposon named Steamer. Steamer mRNA levels and gene copy number correlates with DN and can be used as a marker of the disease. So far, the only mollusk where a retrotransposon expression relates to DN is Mya arenaria. On the other hand, it has been reported that the Chilean blue mussel Mytilus chilensis can also suffers DN. Our aim was to identify retrotransposons in Mytilus chilensis and to study their expression levels in the context of disseminated neoplasia. RESULTS: Here we show that 7.1% of individuals collected in August 2018, from two farming areas, presents morphological characteristics described in DN. Using Steamer sequence to interrogate the transcriptome of M. chilensis we found two putative retrotransposons, named Steamer-like elements (MchSLEs). MchSLEs are present in the genome of M. chilensis and MchSLE1 is indeed an LTR-retrotransposon. Neither expression, nor copy number of the reported MchSLEs correlate with DN status but both are expressed at different levels among individual animals. We also report that in cultured M. chilensis haemocytes MchSLEs1 expression can be induced by bromodeoxyuridine. CONCLUSIONS: We conclude that SLEs present in Mytilus chilensis are differentially expressed among individuals and do not correlate with disseminated neoplasia. Treatment of haemocytes with a stressor like bromodeoxyuridine induces expression of MchSLE1 suggesting that in Mytilus chilensis environmental stressors can induce activation of LTR-retrotransposon.

Contamination acts as a genotype-dependent barrier to gene flow, causing genetic erosion and fine-grained population subdivision in Mussels from the Strait of Istanbul.

This study provides evidence of fine-grained genetic structuring in Mediterranean mussels (Mytilus galloprovincialis) from the Strait of Istanbul, caused by barriers to gene flow via contaminant-mediated selection. In this study, mitochondrial D-loop sequences were analyzed in mussels from 8 localities, all less than 30 kilometers apart, with differing contaminant loads. The results were: 1) Intra-population genetic differentiation (ΦST) between sites with high and low contaminant loads was high (up to 0.459), even at distances of only a few kilometers. 2) Genetic diversity was negatively correlated with the contaminant load ("genetic erosion"). 3) There was evidence of selection, based on haplotype frequencies and neutrality tests (Tajima's D), with purifying selection at the most contaminated site and balancing selection at the least contaminated. 4) Genetic distance was not correlated with geographic distance (no isolation-by-distance), but was correlated with contaminant load at each site. 5) Population dendrograms and Bayesian estimators of migration indicated that gene flow between sites was affected by contamination. For the dendrograms of the sampling sites, the clades clustered according to contaminant load more than geographic distance. Overall, these results suggest that 1) contamination may serve as a genotype-dependent dispersal barrier (i.e., contamination may not affect total number of migrants, just the relative proportions of the haplotypes in the established immigrants), leading strong population differentiation over short distances, and 2) genetic erosion may occur by a combination of selection and altered patterns of haplotype-specific gene flow. These effects may be more pronounced in the Strait of Istanbul than in other locations because of the riverine nature and strong, uni-directional current of the strait.

Transcriptomic Analysis Provides Insights into Candidate Genes and Molecular Pathways Involved in Growth of Mytilus coruscus Larvae.

Growth is a fundamental aspect of aquaculture breeding programs, pivotal for successful cultivation. Understanding the mechanisms that govern growth and development differences across various stages can significantly boost seedling production of economically valuable species, thereby enhancing aquaculture efficiency and advancing the aquaculture industry. Mytilus coruscus, a commercially vital marine bivalve, underscores this importance. To decipher the intricate molecular mechanisms dictating growth and developmental disparities in marine shellfish, we conducted transcriptome sequencing and meticulously analyzed gene expression variations and molecular pathways linked to growth traits in M. coruscus. This study delved into the molecular and gene expression variations across five larval development stages, with a specific focus on scrutinizing the differential expression patterns of growth-associated genes using RNA sequencing and quantitative real-time PCR analysis. A substantial number of genes-36,044 differentially expressed genes (DEGs)-exhibited significant differential expression between consecutive developmental stages. These DEGs were then categorized into multiple pathways (Q value < 0.05), including crucial pathways such as the spliceosome, vascular smooth muscle contraction, DNA replication, and apoptosis, among others. In addition, we identified two pivotal signaling pathways-the Hedgehog (Hh) signaling pathway and the TGF-beta (TGF-ß) signaling pathway-associated with the growth and development of M. coruscus larvae. Ten key growth-related genes were pinpointed, each playing crucial roles in molecular function and the regulation of growth traits in M. coruscus. These genes and pathways associated with growth provide deep insights into the molecular basis of physiological adaptation, metabolic processes, and growth variability in marine bivalves.

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

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

Characterization of NR1J1 Paralog Responses of Marine Mussels: Insights from Toxins and Natural Activators.

The pregnane X receptor (PXR) is a nuclear hormone receptor that plays a pivotal role in regulating gene expression in response to various ligands, particularly xenobiotics. In this context, the aim of this study was to shed light on the ligand affinity and functions of four NR1J1 paralogs identified in the marine mussel Mytilus galloprovincialis, employing a dual-luciferase reporter assay. To achieve this, the activation patterns of these paralogs in response to various toxins, including freshwater cyanotoxins (Anatoxin-a, Cylindrospermopsin, and Microcystin-LR, -RR, and -YR) and marine algal toxins (Nodularin, Saxitoxin, and Tetrodotoxin), alongside natural compounds (Saint John's Wort, Ursolic Acid, and 8-Methoxypsoralene) and microalgal extracts (Tetraselmis, Isochrysis, LEGE 95046, and LEGE 91351 extracts), were studied. The investigation revealed nuanced differences in paralog response patterns, highlighting the remarkable sensitivity of MgaNR1J1γ and MgaNR1J1δ paralogs to several toxins. In conclusion, this study sheds light on the intricate mechanisms of xenobiotic metabolism and detoxification, particularly focusing on the role of marine mussel NR1J1 in responding to a diverse array of compounds. Furthermore, comparative analysis with human PXR revealed potential species-specific adaptations in detoxification mechanisms, suggesting evolutionary implications. These findings deepen our understanding of PXR-mediated metabolism mechanisms, offering insights into environmental monitoring and evolutionary biology research.

Genetic features of bivalve transmissible neoplasia in blue mussels from the Kola Bay (Barents Sea) suggest a recent trans-Arctic migration of the cancer lineages.

Ecology and biogeography of bivalve transmissible neoplasia (BTN) are underexplored due to its recent discovery and a challenging diagnostics. Blue mussels harbour two evolutionary lineages of BTN, MtrBTN1 and MtrBTN2, both derived from Mytilus trossulus. MtrBTN1 has been found only in M. trossulus from North Pacific. MtrBTN2 parasitizes different Mytilus spp. worldwide. BTN in M. trossulus in the Atlantic sector has never been studied. We looked for BTN in mussels from the Barents Sea using flow cytometry of cells, qPCR with primers specific to cancer-associated alleles and sequencing of mtDNA and nuclear loci. Both MtrBTN1 and MtrBTN2 were present in our material, though their prevalence was low (~0.4%). All cancers parasitized M. trossulus except one, MtrBTN1, which was found in a hybrid between M. trossulus and M. edulis. The mtDNA haplotypes found in both lineages were nearly identical to those known from the Northwest Pacific but not from elsewhere. Our results suggest that these two lineages may have arrived in the Barents Sea in recent decades with the maritime transport along the Northern Sea Route. A young evolutionary age of MtrBTN1 seems to indicate that it is an emerging disease in the process of niche expansion. Comparing the new and the published sequence data on tumour suppressor p53, we proved that the prevalence of BTN in mussels can reach epizootic levels. The finding of diverse recombinants between paternally and maternally inherited mtDNAs in somatic tissues of M. trossulus was an unexpected result of our study.

Shedding light on variation in reproductive success through studies of population genetic structure in a Southeast Pacific Coast mussel.

Phylogeography often focuses on the spatial dimension of genetic diversity, rarely including the temporal dynamics occurring interannually among local populations, which can provide insight into past variations in reproductive success. Currently, there is an intense aquaculture industry of Mytilus spp. on the Southeast Pacific Coast which depends entirely on the spat released by natural populations forming a relevant and sensitive social-ecological system. Temporal and spatial spat variability from natural mussel beds could be related to interannual reproductive dynamics with variable reproductive success and recruitment, which leave genetic signatures. Temporal and spatial genetic structure was evaluated in six natural beds in the Southeast Pacific (from 39°25'S to 43°07'S) on the most abundant and widespread Mytilus lineage detected, Mytilus cf. chilensis, in 4 consecutive years. Analyses included data from >180 individuals per year, with a total of 751 (mitochondrial COI) and 747 (nuclear H1) individuals, respectively. Overall, both markers showed high haplotype diversity and low spatial and temporal genetic differentiation. Likely, the high dispersal capacity of Mytilus cf. chilensis maintains population homogeneity and prevents diversity erosion. The slight differences in genetic variance of COI were better explained by differences among sites (space), and conversely, the H1 genetic variance was better explained by interannual (temporal) comparisons, which could explain temporal variability in spat availability. This study highlights the important insights achieved with the evaluation of both temporal and spatial population genetic structures in marine species with high reproductive output, which can condition the success and sustainability of the relevant social-ecological system.

Multi-omics reveals largely distinct transcript- and protein-level responses to the environment in an intertidal mussel.

Organismal responses to stressful environments are influenced by numerous transcript- and protein-level mechanisms, and the relationships between expression changes at these levels are not always straightforward. Here, we used paired transcriptomic and proteomic datasets from two previous studies from gill of the California mussel, Mytilus californianus, to explore how simultaneous transcript and protein abundance patterns may diverge under different environmental scenarios. Field-acclimatized mussels were sampled from two disparate intertidal sites; individuals from one site were subjected to three further treatments (common garden, low-intertidal or high-intertidal outplant) that vary in temperature and feeding time. Assessing 1519 genes shared between the two datasets revealed that both transcript and protein expression patterns differentiated the treatments at a global level, despite numerous underlying discrepancies. There were far more instances of differential expression between treatments in transcript only (1451) or protein only (226) than of the two levels shifting expression concordantly (68 instances). Upregulated expression of cilium-associated transcripts (likely related to feeding) was associated with relatively benign field treatments. In the most stressful treatment, transcripts, but not proteins, for several molecular chaperones (including heat shock proteins and endoplasmic reticulum chaperones) were more abundant, consistent with a threshold model for induction of translation of constitutively available mRNAs. Overall, these results suggest that the relative importance of transcript- and protein-level regulation (translation and/or turnover) differs among cellular functions and across specific microhabitats or environmental contexts. Furthermore, the degree of concordance between transcript and protein expression can vary across benign versus acutely stressful environmental conditions.

Screening of genes encoding proteins that interact with Nrf2: Probing a cDNA library from Mytilus coruscus using a yeast two-hybrid system.

The Nuclear factor Erythroid 2-related factor 2 (Nrf2) is the most important endogenous antioxidant factor in organisms, and it has been demonstrated that it exerts extensive control over the immune response by interacting with crucial innate immunity components directly or indirectly. Although Nrf2 has been widely confirmed to be involved in stress resistance in mammals and some fish, its contribution to mollusks oxidative stress resistance has not frequently been documented. In this investigation, total RNA was taken from the digestive gland of M. coruscus, and a cDNA library was constructed and screened using the GATEWAY recombination technology. The Nrf2 cDNA sequence of M. coruscus was cloned into the pGBKT7 vector to prepare the bait plasmid. Using yeast two-hybrid system, after auxotrophic medium screening, sequencing, and bioinformatics analysis, 13 binding proteins that interacted with Nrf2 were finally identified. They were QM-like protein, 40S ribosomal protein S4 (RPS4), ribosomal protein S2 (RPS2), ribosomal protein L12 (RPL12), EF1-alpha mRNA for elongation factor 1 alpha (eEF1-alpha), ferritin, alpha-amylase, trypsin, vdg3, period clock protein, cyclophilin A isoform 1 (CYP A), serine protease CFSP2, histone variant H2A.Z (H2A.Z). For a better understanding the physiological function of Nrf2 in animals and as a potential target for future research on protein roles in Nrf2 interactions, it is crucial to clarify these protein interactions.

An innovative index to incorporate transcriptomic data into weight of evidence approaches for environmental risk assessment.

The sharp decrease in the cost of RNA-sequencing and the rapid improvement in computational analysis of eco-toxicogenomic data have brought new insights into the adverse effects of chemicals on aquatic organisms. Yet, transcriptomics is generally applied qualitatively in environmental risk assessments, hampering more effective exploitation of this evidence through multidisciplinary studies. In view of this limitation, a methodology is here presented to quantitatively elaborate transcriptional data in support to environmental risk assessment. The proposed methodology makes use of results from the application of Gene Set Enrichment Analysis to recent studies investigating the response of Mytilus galloprovincialis and Ruditapes philippinarum exposed to contaminants of emerging concern. The degree of changes in gene sets and the relevance of physiological reactions are integrated in the calculation of a hazard index. The outcome is then classified according to five hazard classes (from absent to severe), providing an evaluation of whole-transcriptome effects of chemical exposure. The application to experimental and simulated datasets proved that the method can effectively discriminate different levels of altered transcriptomic responses when compared to expert judgement (Spearman correlation coefficient of 0.96). A further application to data collected in two independent studies of Salmo trutta and Xenopus tropicalis exposed to contaminants confirmed the potential extension of the methodology to other aquatic species. This methodology can serve as a proof of concept for the integration of "genomic tools" in environmental risk assessment based on multidisciplinary investigations. To this end, the proposed transcriptomic hazard index can now be incorporated into quantitative Weight of Evidence approaches and weighed, with results from other types of analysis, to elucidate the role of chemicals in adverse ecological effects.

Outer membrane vesicles induce the mussel plantigrade settlement via regulation of c-di-GMP.

Despite the importance of outer membrane vesicles (OMVs) in benthic animal settlement, the underlying molecular mechanism remains elusive. Here, the impact of OMVs and OMVs synthesis-related tolB gene in Mytilus coruscus plantigrade settlement was tested. The OMVs were extracted from Pseudoalteromonas marina through density gradient centrifugation, and a tolB knockout strain, achieved by homologous recombination, was utilized for the investigation. Our results demonstrated that OMVs could significantly enhance M. coruscus plantigrades settlement. Deleting the tolB resulted in downregulation of c-di-GMP, accompanied by a reduction of OMV production, a decline in bacterial motility and increasing biofilm-forming ability. Enzyme treatment resulted in a 61.11% reduction in OMV-inducing activity and a 94.87% reduction in LPS content. Thus, OMVs regulate mussel settlement via LPS, and c-di-GMP is responsible for the OMV-inducing capacity. These findings provide new insights into the interactions between bacteria and mussels.

Genetic Diversity, Population Structure, and Environmental Adaptation Signatures of Chinese Coastal Hard-Shell Mussel Mytilus coruscus Revealed by Whole-Genome Sequencing.

The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (FST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

Two-Component System Response Regulator ompR Regulates Mussel Settlement through Exopolysaccharides.

The outer membrane protein (OMP) is a kind of biofilm matrix component that widely exists in Gram-negative bacteria. However, the mechanism of OMP involved in the settlement of molluscs is still unclear. In this study, the mussel Mytilus coruscus was selected as a model to explore the function of ompR, a two-component system response regulator, on Pseudoalteromonas marina biofilm-forming capacity and the mussel settlement. The motility of the ΔompR strain was increased, the biofilm-forming capacity was decreased, and the inducing activity of the ΔompR biofilms in plantigrades decreased significantly (p < 0.05). The extracellular α-polysaccharide and ß-polysaccharide of the ΔompR strain decreased by 57.27% and 62.63%, respectively. The inactivation of the ompR gene decreased the ompW gene expression and had no impact on envZ expression or c-di-GMP levels. Adding recombinant OmpW protein caused the recovery of biofilm-inducing activities, accompanied by the upregulation of exopolysaccharides. The findings deepen the understanding of the regulatory mechanism of bacterial two-component systems and the settlement of benthic animals.