Cadmium Highlights Common and Specific Responses of Two Freshwater Sentinel Species, Dreissena polymorpha and Dreissena rostriformis bugensis.

Zebra mussel (ZM), Dreissena polymorpha, commonly used as a sentinel species in freshwater biomonitoring, is now in competition for habitat with quagga mussel (QM), Dreissena rostriformis bugensis. This raises the question of the quagga mussel's use in environmental survey. To better characterise QM response to stress compared with ZM, both species were exposed to cadmium (100 µg·L-1), a classic pollutant, for 7 days under controlled conditions. The gill proteomes were analysed using two-dimensional electrophoresis coupled with mass spectrometry. For ZM, 81 out of 88 proteoforms of variable abundance were identified using mass spectrometry, and for QM, 105 out of 134. Interestingly, the proteomic response amplitude varied drastically, with 5.6% of proteoforms of variable abundance (DAPs) in ZM versus 9.4% in QM. QM also exhibited greater cadmium accumulation. Only 12 common DAPs were observed. Several short proteoforms were detected, suggesting proteolysis. Functional analysis is consistent with the pleiotropic effects of the toxic metal ion cadmium, with alterations in sulphur and glutathione metabolisms, cellular calcium signalling, cytoskeletal dynamics, energy production, chaperone activation, and membrane events with numerous proteins involved in trafficking and endocytosis/exocytosis processes. Beyond common responses, the sister species display distinct reactions, with cellular response to stress being the main category involved in ZM as opposed to calcium and cytoskeleton alterations in QM. Moreover, QM exhibited greater evidence of proteolysis and cell death. Overall, these results suggest that QM has a weaker stress response capacity than ZM.

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica).

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria.

The increasing concentration of CO2 in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO2 and reared under acidified (pH ~ 7.5, pCO2 ~ 1200 ppm) or control (pH ~ 7.9, pCO2 ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

Two parallel chromosome-level reference genomes to support restoration and aquaculture of European flat oyster Ostrea edulis.

This volume of Evolutionary Applications sees the publication of two genomes for the European native flat oyster Ostrea edulis, a species of significant evolutionary, ecological and commercial value. Each is a highly contiguous chromosome-level assembly from individuals of different genetic backgrounds, which have been benchmarked against one another. This situation has resulted from the serendipitous discovery that two independent research groups were both deep into the process of building, annotating and investigating separately produced assemblies. Due to constraints with funder requirements and the need to recognize early career researchers for their work, alongside the technical challenge of integrating assemblies from two very different genomes, there was limited capacity to merge the sequences into one publication at the stage of discovery. This issue is likely to become very common over the next few years until the technologies for working with multiple genomes at once, for example, graph genomes, become commonplace in nonmodel species. Consequently, both of our teams have decided to collaborate rather than compete, recognizing the benefit to copublishing two separate genome resources for the research community, each with distinct scientific investigations, and working collaboratively to benchmark the assemblies.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture.

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.).

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks.

BACKGROUND: The hard clam Mercenaria mercenaria is a major marine resource along the Atlantic coasts of North America and has been introduced to other continents for resource restoration or aquaculture activities. Significant mortality events have been reported in the species throughout its native range as a result of diseases (microbial infections, leukemia) and acute environmental stress. In this context, the characterization of the hard clam genome can provide highly needed resources to enable basic (e.g., oncogenesis and cancer transmission, adaptation biology) and applied (clam stock enhancement, genomic selection) sciences. RESULTS: Using a combination of long and short-read sequencing technologies, a 1.86 Gb chromosome-level assembly of the clam genome was generated. The assembly was scaffolded into 19 chromosomes, with an N50 of 83 Mb. Genome annotation yielded 34,728 predicted protein-coding genes, markedly more than the few other members of the Venerida sequenced so far, with coding regions representing only 2% of the assembly. Indeed, more than half of the genome is composed of repeated elements, including transposable elements. Major chromosome rearrangements were detected between this assembly and another recent assembly derived from a genetically segregated clam stock. Comparative analysis of the clam genome allowed the identification of a marked diversification in immune-related proteins, particularly extensive tandem duplications and expansions in tumor necrosis factors (TNFs) and C1q domain-containing proteins, some of which were previously shown to play a role in clam interactions with infectious microbes. The study also generated a comparative repertoire highlighting the diversity and, in some instances, the specificity of LTR-retrotransposons elements, particularly Steamer elements in bivalves. CONCLUSIONS: The diversity of immune molecules in M. mercenaria may allow this species to cope with varying and complex microbial and environmental landscapes. The repertoire of transposable elements identified in this study, particularly Steamer elements, should be a prime target for the investigation of cancer cell development and transmission among bivalve mollusks.

Does the environmental history of mussels have an effect on the physiological response to additional stress under experimental conditions?

Expected effects on marine biota of the ongoing elevation of water temperature and high latitudes is of major concern when considering the reliability of coastal ecosystem production. To compare the capacity of coastal organisms to cope with a temperature increase depending on their environmental history, responses of adult blue mussels (Mytilus spp.) taken from two sites differentially exposed to chemical pollution were investigated during an experimental exposure to a thermal stress. Immune parameters were notably altered by extreme warming and transcriptional changes for a broad selection of genes were associated to the temperature increase following a two-step response pattern. Site-specific responses suggested an influence of environmental history and support the possibility of a genetic basis in the physiological response. However no meaningful difference was detected between the response of hybrids and M galloprovincialis. This study brings new information about the capacity of mussels to cope with the ongoing elevation of water temperature in these coastal ecosystems.

Identification of variants associated with hard clam, Mercenaria mercenaria, resistance to Quahog Parasite Unknown disease.

Severe losses in aquacultured and wild hard clam (Mercenaria mercenaria) stocks have been previously reported in the northeastern United States due to a protistan parasite called QPX (Quahog Parasite Unknown). Previous work demonstrated that clam resistance to QPX is under genetic control. This study identifies single nucleotide polymorphism (SNP) associated with clam survivorship from two geographically segregated populations, both deployed in an enzootic site. The analysis contrasted samples collected before and after undergoing QPX-related mortalities and relied on a robust draft clam genome assembly. ~200 genes displayed significant variant enrichment at each sampling point in both populations, including 18 genes shared between both populations. Markers from both populations were identified in genes related to apoptosis pathways, protein-protein interaction, receptors, and signaling. This research begins to identify genetic markers associated with clam resistance to QPX disease, leading the way for the development of resistant clam stocks through marker-assisted selection.

Biological rhythms in the deep-sea hydrothermal mussel Bathymodiolus azoricus.

Biological rhythms are a fundamental property of life. The deep ocean covers 66% of our planet surface and is one of the largest biomes. The deep sea has long been considered as an arrhythmic environment because sunlight is totally absent below 1,000 m depth. In the present study, we have sequenced the temporal transcriptomes of a deep-sea species, the ecosystem-structuring vent mussel Bathymodiolus azoricus. We reveal that tidal cycles predominate in the transcriptome and physiology of mussels fixed directly at hydrothermal vents at 1,688 m depth at the Mid-Atlantic Ridge, whereas daily cycles prevail in mussels sampled after laboratory acclimation. We identify B. azoricus canonical circadian clock genes, and show that oscillations observed in deep-sea mussels could be either a direct response to environmental stimulus, or be driven endogenously by one or more biological clocks. This work generates in situ insights into temporal organisation in a deep-sea organism.

A review of the potential genes implicated in follicular atresia in teleost fish.

In recent years, numerous studies conducted on teleost fish have highlighted the contribution of transcriptomic studies in elucidating the physiological mechanisms underlying the molecular events of oogenesis and follicular atresia, enabling the identification of potential genes and molecular networks that participate in both the reproductive cycle and the process of follicular atresia. Atresia can affect the reproductive potential of females by reducing the healthy eggs that a female can spawn in both aquaculture and wild populations. The substantial diversity of reproductive strategies exhibited by teleost fish has contributed to the difficulty in identifying common genes between species, but a set of core genes has emerged as potential markers for atresia in relation to apoptosis/autophagy, lipid metabolism, oxidative metabolism and other physiological processes similar to those identified in other vertebrates, even mammals. We review the current status of the genes that have been identified in ovaries with atretic oocytes. Our primary goal is to review the current status regarding gene expression during gonadal development and follicular atresia. This information will enable us to understand the factors and expression patterns involved in the follicular atresia of teleost fish.

Global host molecular perturbations upon in situ loss of bacterial endosymbionts in the deep-sea mussel Bathymodiolus azoricus assessed using proteomics and transcriptomics.

BACKGROUND: Colonization of deep-sea hydrothermal vents by most invertebrates was made efficient through their adaptation to a symbiotic lifestyle with chemosynthetic bacteria, the primary producers in these ecosystems. Anatomical adaptations such as the establishment of specialized cells or organs have been evidenced in numerous deep-sea invertebrates. However, very few studies detailed global inter-dependencies between host and symbionts in these ecosystems. In this study, we proposed to describe, using a proteo-transcriptomic approach, the effects of symbionts loss on the deep-sea mussel Bathymodiolus azoricus' molecular biology. We induced an in situ depletion of symbionts and compared the proteo-transcriptome of the gills of mussels in three conditions: symbiotic mussels (natural population), symbiont-depleted mussels and aposymbiotic mussels. RESULTS: Global proteomic and transcriptomic results evidenced a global disruption of host machinery in aposymbiotic organisms. We observed that the total number of proteins identified decreased from 1118 in symbiotic mussels to 790 in partially depleted mussels and 761 in aposymbiotic mussels. Using microarrays we identified 4300 transcripts differentially expressed between symbiont-depleted and symbiotic mussels. Among these transcripts, 799 were found differentially expressed in aposymbiotic mussels and almost twice as many in symbiont-depleted mussels as compared to symbiotic mussels. Regarding apoptotic and immune system processes - known to be largely involved in symbiotic interactions - an overall up-regulation of associated proteins and transcripts was observed in symbiont-depleted mussels. CONCLUSION: Overall, our study showed a global impairment of host machinery and an activation of both the immune and apoptotic system following symbiont-depletion. One of the main assumptions is the involvement of symbiotic bacteria in the inhibition and regulation of immune and apoptotic systems. As such, symbiotic bacteria may increase their lifespan in gill cells while managing the defense of the holobiont against putative pathogens.

Identification and gene expression of multiple peptidoglycan recognition proteins (PGRPs) in the deep-sea mussel Bathymodiolus azoricus, involvement in symbiosis?

The relationship between the deep-sea mussel Bathymodiolus azoricus and its thiotrophic (SOX) and methanotrophic (MOX) symbionts has been ecologically and functionally well studied. Endosymbiosis is common in deep-sea hydrothermal vent fauna, yet little is known about the molecular mechanisms underlying the regulation of interactions between host and symbionts. In this study we focused on a group of pattern recognition receptors (PRR), called PGRPs that are able to recognize the peptidoglycan of bacterial cell wall. We first characterised the different PGRPs isoforms in B. azoricus gills and identified five paralogs. Among them two displayed a signal peptide. Then, specific probes designed for each paralog were used to perform real-time PCR quantification in gills of individuals showing various bacterial content as a result of in situ experimental procedures. Overall we found a decrease of PGRPs expression when symbionts amount decreases, suggesting an implication of PGRPs in the regulation of symbionts in B. azoricus gills. We therefore hypothesize that secreted proteins could act as cooperation signals to induce colonisation of symbiotic tissue while non-secreted proteins may regulate the density of endosymbionts within the gill tissue.

Differential Gene Expression in Five Isolates of the Clam Pathogen, Quahog Parasite Unknown (QPX).

Quahog parasite unknown (QPX) is a thraustochytrid protist that infects the hard clam, Mercenaria mercenaria, causing significant economic losses along the northeastern coast of North America. Previous investigations noted differences in growth dynamics and virulence in QPX cells from different geographic locations. In order to probe the molecular determinants for these variations, we investigated the transcriptomic profiles of five geographically distinct QPX isolates using custom 15k 60-mer oligonucleotide arrays. A total of 1,263 transcripts were differentially expressed (DE) among the five QPX isolates. The hierarchical clustering of gene expression profiles showed that the QPX isolates from Raritan Bay (RB, NY) and from Provincetown Harbor (MA) were more similar to each other and diverged from QPX isolates from Peconic Bay (PB, NY) and Old Plantation Creek (VA), which had more similar gene expression profiles. The most prominent difference was based on 78 transcripts coding for heat shock proteins DE between the five QPX isolates. The study generated contrasting transcriptomic profiles for QPX isolated from northern (MA) and deeper (RB, NY) locations as compared to southern (VA) and shallower (PB, NY) areas, suggesting the adaptation of the parasite to local environmental, in particular temperature, conditions.

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages.

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.

Multiple I-Type Lysozymes in the Hydrothermal Vent Mussel Bathymodiolus azoricus and Their Role in Symbiotic Plasticity.

The aim of this study was first to identify lysozymes paralogs in the deep sea mussel Bathymodiolus azoricus then to measure their relative expression or activity in different tissue or conditions. B. azoricus is a bivalve that lives close to hydrothermal chimney in the Mid-Atlantic Ridge (MAR). They harbour in specialized gill cells two types of endosymbiont (gram-bacteria): sulphide oxidizing bacteria (SOX) and methanotrophic bacteria (MOX). This association is thought to be ruled by specific mechanism or actors of regulation to deal with the presence of symbiont but these mechanisms are still poorly understood. Here, we focused on the implication of lysozyme, a bactericidal enzyme, in this endosymbiosis. The relative expression of Ba-lysozymes paralogs and the global anti-microbial activity, were measured in natural population (Lucky Strike--1700 m, Mid-Atlantic Ridge), and in in situ experimental conditions. B. azoricus individuals were moved away from the hydrothermal fluid to induce a loss of symbiont. Then after 6 days some mussels were brought back to the mussel bed to induce a re-acquisition of symbiotic bacteria. Results show the presence of 6 paralogs in B. azoricus. In absence of symbionts, 3 paralogs are up-regulated while others are not differentially expressed. Moreover the global activity of lysozyme is increasing with the loss of symbiont. All together these results suggest that lysozyme may play a crucial role in symbiont regulation.

Alterations of the immune transcriptome in resistant and susceptible hard clams (Mercenaria mercenaria) in response to Quahog Parasite Unknown (QPX) and temperature.

Quahog Parasite Unknown (QPX) is a fatal protistan parasite that causes severe losses in the hard clam (Mercenaria mercenaria) fisheries along the northeastern coast of the US. Field and laboratory studies of QPX disease have demonstrated a major role for water temperature and M. mercenaria genetic origin in disease development. Infections are more likely to occur at cold temperatures, with clam stocks originating from southern states being more susceptible than clams from northern origin where disease is enzootic. Even though the influence of temperature on QPX infection have been examined in susceptible and resistant M. mercenaria at physiological and cellular scales, the underlying molecular mechanisms associated with host-pathogen interactions remain largely unknown. This study was carried out to explore the molecular changes in M. mercenaria in response to temperature and QPX infection on the transcriptomic level, and also to compare molecular responses between susceptible and resistant clam stocks. A M. mercenaria oligoarray (15 K Agilent) platform was produced based on our previously generated transcriptomic data and was used to compare gene expression profiles in naive and QPX-infected susceptible (Florida stock) and resistant (Massachusetts) clams maintained at temperatures favoring disease development (13 °C) or clam healing (21 °C). In addition, transcriptomic changes reflecting focal (the site of infection, mantle) and systemic (circulating hemocytes) responses were also assessed using the oligoarray platform. Results revealed significant regulation of multiple biological pathways by temperature and QPX infection, mainly associated with immune recognition, microbial killing, protein synthesis, oxidative protection and metabolism. Alterations were widely systemic with most changes in gene expression revealed in hemocytes, highlighting the role of circulating hemocytes as the first line of defense against pathogenic stress. A large number of complement-related recognition molecules with fibrinogen or C1q domains were shown to be specially induced following QPX challenge, and the expression of these molecules was significantly higher in resistant clams as compared to susceptible ones. These highly variable immune proteins may be potent candidate molecular markers for future study of M. mercenaria resistance against QPX. Beyond the specific case of clam response to QPX, this study also provides insights into the primitive complement-like system in the hard clam.

Predicting growth and mortality of bivalve larvae using gene expression and supervised machine learning.

It is commonly known that the nature of the diet has diverse consequences on larval performance and longevity, however it is still unclear which genes have critical impacts on bivalve development and which pathways are of particular importance in their vulnerability or resistance. First we show that a diet deficient in essential fatty acid (EFA) produces higher larval mortality rates, a reduced shell growth, and lower postlarval performance, all of which are positively correlated with a decline in arachidonic and eicosapentaenoic acids levels, two EFAs known as eicosanoid precursors. Eicosanoids affect the cell inflammatory reactions and are synthesized from long-chain EFAs. Second, we show for the first time that a deficiency in eicosanoid precursors is associated with a network of 29 genes. Their differential regulation can lead to slower growth and higher mortality of Mytilus edulis larvae. Some of these genes are specific to bivalves and others are implicated at the same time in lipid metabolism and defense. Several genes are expressed only during pre-metamorphosis where they are essential for muscle or neurone development and biomineralization, but only in stress-induced larvae. Finally, we discuss how our networks of differentially expressed genes might dynamically alter the development of marine bivalves, especially under dietary influence.

The use of -omic tools in the study of disease processes in marine bivalve mollusks.

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Transcriptomic response of the hydrothermal mussel Bathymodiolus azoricus in experimental exposure to heavy metals is modulated by the Pgm genotype and symbiont content.

Hydrothermal vent mussels belonging to the genus Bathymodiolus dominate communities at hydrothermal sites of the Mid-Atlantic Ridge. The mussel Bathymodiolus azoricus harbors thiotrophic and methanotrophic symbiotic bacteria in its gills and evolves in naturally highly metal contaminated environments. In the context of investigations on metal tolerance/effect in B. azoricus, we focused our work on the short-term adaptive response (15days) of mussels to different metals exposure at a molecular level using metal concentrations chosen to mimic natural situations at three vents sites. The expression of a set of 38 genes involved in different steps of the metal uptake, detoxication and various metabolisms was analysed by qPCR. Mussels were also genotyped at 10 enzyme loci to explore the relationships among natural genetic variation and gene expression. Relation between symbiont content (both sulfur-oxidizing and methanogen bacteria) and gene expression was also analysed. Our study demonstrated the influence of metal cocktail composition and time exposure on the transcriptome regulation with a specific pattern of regulation observed for the three metal cocktail tested. We also evidenced the significant influence of some specific Pgm genotype on the global gene expression in our experimental populations and a general trend of a higher gene expression in individuals carrying a high symbiont content.