Proteomic signatures of the oyster metabolic response to herpesvirus OsHV-1 μVar infection

Harsh intertidal environment enhances metabolism and immunity in oyster (Crassostrea gigas) spat

The Pacific oyster Crassostrea gigas is established in the marine intertidal zone, experiencing rapid and highly dynamic environmental changes throughout the tidal cycle. Depending on the bathymetry, oysters face oxygen deprivation, lack of nutrients, and high changes in temperature during alternation of the cycles of emersion/immersion. Here we showed that intertidal oysters at a bathymetry level of 3 and 5 m delayed by ten days the onset of mortality associated with Pacific Oyster Mortality Syndrome (POMS) as compared to subtidal oysters. Intertidal oysters presented a lower growth but similar energetic reserves to subtidal oysters but induced proteomic changes indicative of a boost in metabolism, inflammation, and innate immunity that may have improved their resistance during infection with the Ostreid herpes virus. Our work highlights that intertidal harsh environmental conditions modify host-pathogen interaction and improve oyster health. This study opens new perspectives on oyster farming for mitigation strategies based on tidal height.

Pooled resequencing of larvae and adults reveals genomic variations associated with Ostreid herpesvirus 1 resistance in the Pacific oyster Crassostrea gigas

The Ostreid herpesvirus 1 (OsHV-1) is a lethal pathogen of the Pacific oyster (Crassostrea gigas), an important aquaculture species. To understand the genetic architecture of the defense against the pathogen, we studied genomic variations associated with herpesvirus-caused mortalities by pooled whole-genome resequencing of before and after-mortality larval samples as well as dead and surviving adults from a viral challenge. Analysis of the resequencing data identified 5,271 SNPs and 1,883 genomic regions covering 3,111 genes in larvae, and 18,692 SNPs and 28,314 regions covering 4,863 genes in adults that were significantly associated with herpesvirus-caused mortalities. Only 1,653 of the implicated genes were shared by larvae and adults, suggesting that the antiviral response or resistance in larvae and adults involves different sets of genes or differentiated members of expanded gene families. Combined analyses with previous transcriptomic data from challenge experiments revealed that transcription of many mortality-associated genes was also significantly upregulated by herpesvirus infection confirming their importance in antiviral response. Key immune response genes especially those encoding antiviral receptors such as TLRs and RLRs displayed strong association between variation in regulatory region and herpesvirus-caused mortality, suggesting they may confer resistance through transcriptional modulation. These results point to previously undescribed genetic mechanisms for disease resistance at different developmental stages and provide candidate polymorphisms and genes that are valuable for understanding antiviral immune responses and breeding for herpesvirus resistance.

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.

Shotgun analysis to identify differences in protein expression between granulocytes and hyalinocytes of the European flat oyster Ostrea edulis

Recovery of wild populations of the European flat oyster Ostrea edulis is important for ecosystem health and conservation of this species, because native oyster populations have dramatically declined or disappeared in most European waters. Diseases have contributed to oyster decline and are important constrains for oyster recovery. Understanding oyster immune system should contribute to design effective strategies to fight oyster diseases. Haemocytes play a pivotal role in mollusc immune responses protecting from infection. Two main types of haemocytes, granulocytes and hyalinocytes, are distinguished in O. edulis. A study aiming to explore differential functions between both haemocyte types and, thus, to enrich the knowledge of Ostrea edulis immune system, was performed by comparing the proteome of the two haemolymph cell types, using a shotgun approach through liquid chromatography (LC) coupled to mass spectrometry (MS). Cells from oyster haemolymph were differentially separated by Percoll density gradient centrifugation. Shotgun LC-MS/MS performance allowed the identification of 145 proteins in hyalinocytes and 138 in the proteome of granulocytes. After a comparative analysis, 55 proteins with main roles in defence were identified, from which 28 were representative of granulocytes and 27 of hyalinocytes, plus 11 proteins shared by both cell types. Different proteins involved in signal transduction, apoptosis, oxidative response, processes related with the cytoskeleton and structure, recognition and wound healing were identified as representatives of each haemocyte type. Important signalling pathways in the immune response such as MAPK, Ras and NF-κβ seemed to be more relevant for granulocytes, while the Wnt signalling pathway, particularly relevant for wound healing, more relevant in hyalinocytes. The differences in proteins involved in recognition and in cytoskeleton and structure suggest differential specialisation in processes of phagocytosis and internalisation of pathogens between haemocyte types. Apoptosis seemed more active in granulocytes. The differences in proteins involved in oxidative response also suggest different redox processes in each cell type.

Proteome Response of Meretrix Bivalves Hepatopancreas Exposed to Paralytic Shellfish Toxins Producing Dinoflagellate Gymnodinium catenatum

Paralytic shellfish toxins (PSTs) contamination of seafood has become a growing global problem. However, the molecular response of bivalves, some of the most popular seafoods, to PSP toxins has seldom been reported and the underlying molecular mechanisms of the interactions between Meretrix meretrix bivalves and PSTs-producing dinoflagellates are scarcely known. This study compared the protein expression profiles between PSP toxin-contaminated and non-PSP toxin contaminated M. meretrix, determined proteome responses and identified potential biomarkers based on feeding experiments. Results showed that the content of total PSP toxins in contaminated bivalves was 40.63 ± 4.08 μg saxitoxin (STX) equivalents per gram, with 95.3% in hepatopancreas, followed by gill (1.82%) and foot (1.79%). According to two-dimensional gel electrophoresis (2-DE), 15 differentially expressed proteins (at least 2-fold difference) between the hepatopancreas of bivalves with and without PSP toxins were detected. Eight of them were successfully identified by MALDI-TOF MS. These were catalase, protein ultraspiracle homolog, G2 and S phase-expression protein, paramyosin, Mn-superoxide dismutase, response regulator receiver domain-containing protein, sarcoplasmic calcium-binding protein and major facilitator superfamily transporters. The differences in the expression levels of the last three proteins involving in cell signaling, structure and membrane transport were 4.2, 5.3 and 4.9-fold, respectively. These proteins could be further developed as potential biomarkers. The other two up-regulated proteins, Mn-superoxide dismutase and catalase, were involved in cell defence mechanisms against oxidative stress, suggesting PSP toxin acts as xenobiotics and poses oxidative stress in bivalves. This study gives insights into the response of bivalves to PSP toxin-producing dinoflagellate at the proteomic level and the potential of using 2-DE to develop specific protein markers in bivalves.

Effect of Different Species of Prorocentrum Genus on the Japanese Oyster Crassostrea gigas Proteomic Profile

This paper assesses the effects of exposure to toxic concentrations (1200 to 6000 cells/mL) of the dinoflagellates Prorocentrum lima, Prorocentrum minimum, and Prorocentrum rhathymum and several concentrations of aqueous and organic extracts obtained from the same species (0 to 20 parts per thousand) on the Crassostrea gigas (5-7 mm) proteomic profile. Through comparative proteomic map analyses, several protein spots were detected with different expression levels, of which eight were selected to be identified by liquid chromatography-mass spectrometry (LC-MS/MS) analyses. The proteomic response suggests that, after 72 h of exposure to whole cells, the biological functions of C. gigas affected proteins in the immune system, stress response, contractile systems and cytoskeletal activities. The exposure to organic and aqueous extracts mainly showed effects on protein expressions in muscle contraction and cytoskeleton morphology. These results enrich the knowledge on early bivalve developmental stages. Therefore, they may be considered a solid base for new bioassays and/or generation of specific analytical tools that allow for some of the main effects of algal proliferation phenomena on bivalve mollusk development to be monitored, characterized and elucidated.

In situ characterisation of pathogen dynamics during a Pacific oyster mortality syndrome episode

Significant mortality of Crassostrea gigas juveniles is observed systematically every year worldwide. Pacific Oyster Mortality Syndrome (POMS) is caused by Ostreid Herpesvirus 1 (OsHV-1) infection leading to immune suppression, followed by bacteraemia caused by a consortium of opportunistic bacteria. Using an in-situ approach and pelagic chambers, our aim in this study was to identify pathogen dynamics in oyster flesh and in the water column during the course of a mortality episode in the Mediterranean Thau lagoon (France). OsHV-1 concentrations in oyster flesh increased before the first clinical symptoms of the disease appeared, reached maximum concentrations during the moribund phase and the mortality peak. The structure of the bacterial community associated with oyster flesh changed in favour of bacterial genera previously associated with oyster mortality including Vibrio, Arcobacter, Psychrobium, and Psychrilyobacter. During the oyster mortality episode, releases of OsHV-1 and opportunistic bacteria were observed, in succession, in the water surrounding the oyster lanterns. These releases may favour the spread of disease within oyster farms and potentially impact other marine species, thereby reducing marine biodiversity in shellfish farming areas.

Comparative Proteomics of Ostreid Herpesvirus 1 and Pacific Oyster Interactions With Two Families Exhibiting Contrasted Susceptibility to Viral Infection

Massive mortality outbreaks affecting Pacific oysters (Crassostrea gigas) spat/juveniles are often associated with the detection of a herpesvirus called ostreid herpesvirus type 1 (OsHV-1). In this work, experimental infection trials of C. gigas spat with OsHV-1 were conducted using two contrasted Pacific oyster families for their susceptibility to viral infection. Live oysters were sampled at 12, 26, and 144 h post infection (hpi) to analyze host-pathogen interactions using comparative proteomics. Shotgun proteomics allowed the detection of seven viral proteins in infected oysters, some of them with potential immunomodulatoy functions. Viral proteins were mainly detected in susceptible oysters sampled at 26 hpi, which correlates with the mortality and viral load observed in this oyster family. Concerning the Pacific oyster proteome, more than 3,000 proteins were identified and contrasted proteomic responses were observed between infected A- and P-oysters, sampled at different post-injection times. Gene ontology (GO) and KEGG pathway enrichment analysis performed on significantly modulated proteins uncover the main immune processes (such as RNA interference, interferon-like pathway, antioxidant defense) which contribute to the defense and resistance of Pacific oysters to viral infection. In the more susceptible Pacific oysters, results suggest that OsHV-1 manipulate the molecular machinery of host immune response, in particular the autophagy system. This immunomodulation may lead to weakening and consecutively triggering death of Pacific oysters. The identification of several highly modulated and defense-related Pacific oyster proteins from the most resistant oysters supports the crucial role played by the innate immune system against OsHV-1 and the viral infection. Our results confirm the implication of proteins involved in an interferon-like pathway for efficient antiviral defenses and suggest that proteins involved in RNA interference process prevent viral replication in C. gigas. Overall, this study shows the interest of multi-omic approaches applied on groups of animals with differing sensitivities and provides novel insight into the interaction between Pacific oyster and OsHV-1 with key proteins involved in viral infection resistance.

Autophagy Dually Induced by AMP Surplus and Oxidative Stress Enhances Hemocyte Survival and Bactericidal Capacity via AMPK Pathway in Crassostrea hongkongensis

Crassostrea hongkongensis (Hong Kong oyster) is an ecologically and economically valuable shellfish endemic to South/Southeast Asia. Due to ocean acidification and warming waters, they have become increasingly vulnerable to invading microbes including Vibrio parahaemolyticus, a significant foodborne human pathogen. In recent years, outbreaks of V. parahaemolyticus have emerged as a perennial phenomenon in parts of the world, necessitating to better understand the biology of host-pathogen interactions in this under-examined marine invertebrate. Although an immunologically relevant autophagy apparatus has been identified in Crassostrea gigas, an evolutionarily close mollusk cousin, the precise mechanistic details of C. hongkongensis autophagy during V. parahaemolyticus infection are still wanting. Here, we compellingly demonstrated that in vivo V. parahaemolyticus challenge robustly triggered autophagic signaling in C. hongkongensis hemocytes peaking at 6 h post-infection, which subsequently promoted bacterial clearance and dampened premature apoptosis. Simultaneously, a large surplus of adenosine monophosphate (AMP) and elevations in reactive oxygen species (ROS, specifically mitochondrial O₂^– and cellular H₂O₂) formation were observed post-infection. Extrinsically applied AMP and ROS could synergistically induce AMP-activated protein kinase (AMPK) phosphorylation to stimulate downstream autophagic events. V. parahaemolyticus infection-induced autophagy was pharmacologically shown to be AMPK-dependent in vivo. Overall, our results establish autophagy as a crucial arm of host defense against Vibrio infections in mollusks, and provide new insights into the underappreciated roles of ROS and AMP as co-regulators of autophagy.

OsHV-1 infection leads to mollusc tissue lesion and iron redistribution, revealing a strategy of iron limitation against pathogen

The mass mortality of molluscs caused by OsHV-1 infection has frequently occurred worldwide in recent years. Meanwhile the interaction between OsHV-1 and its host is largely unknown. Innate immunity mainly makes up the mollusc defense system, due to the lack of adaptive immunity in invertebrates. The iron limitation strategy is an indispensable facet of innate immunity across vertebrate and invertebrate species. In this study, an iron limitation strategy was interestingly found to contribute to mollusc innate immune responses against OsHV-1 infection. Firstly, ark clams, Scapharca broughtonii, were experimentally infected with OsHV-1, and serious hyperaemia in hepatopancreases and the erosion of gills were observed post OsHV-1 infection according to a histology assay. Meanwhile, based on quantification and Prussian blue staining, the process of iron efflux from ark clams was described post OsHV-1 infection. Secondly, ferritin, as an important iron storage protein, was characterized in ark clams and showed significant iron binding activity. According to the results of an immunohistochemistry assay, ferritin was supposed to be responsible for the iron translocation in ark clams post OsHV-1 infection. Its expression level was significantly fluctuant in response to OsHV-1 infection. Finally, oxidative stress was assessed by the analyses of H2O2 content, total antioxidant capacity and MDA level post OsHV-1 infection. Supplementary iron was found to promote ROS generation and death of hemocytes in vivo. These results highlighted that microenvironment changes in the essential nutrient iron should be an important aspect of the pathogenesis of OsHV-1 disease.

High temperature induces transcriptomic changes in Crassostrea gigas that hinders progress of Ostreid herpesvirus (OsHV-1) and promotes survival

Among all the environmental factors, seawater temperature plays a decisive role in triggering marine diseases. Like fever in vertebrates, high seawater temperature could modulate the host response to the pathogens in ectothermic animals. In France, massive mortality of Pacific oysters Crassostrea gigas caused by the ostreid herpesvirus 1 (OsHV-1) is markedly reduced when temperatures exceed 24°C in the field. In the present study we assess how high temperature influences the host response to the pathogen by comparing transcriptomes (RNA-sequencing) during the course of experimental infection at 21°C (reference) and 29°C. We show that high temperature induced host physiological processes that are unfavorable to the viral infection. Temperature influenced the expression of transcripts related to the immune process and increased the transcription of genes related to apoptotic process, synaptic signaling, and protein processes at 29°C. Concomitantly, the expression of genes associated to catabolism, metabolites transport, macromolecules synthesis and cell growth remained low since the first stage of infection at 29°C. Moreover, viral entry into the host might have been limited at 29°C by changes in extracellular matrix composition and protein abundance. Overall, these results provide new insights into how environmental factors modulate the host-pathogen interactions.

Investigating the natural resistance of blackfoot p-a%%KERN_ERR%%ua Haliotis iris to abalone viral ganglioneuritis using whole transcriptome analysis

The natural resistance of New Zealand blackfoot p-a%%%%%%%%%%%%%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%ua Haliotis iris to infection by haliotid herpesvirus-1 (HaHV-1) and to the disease abalone viral ganglioneuritis was investigated in experimentally challenged p-aua using high throughput RNA-sequencing. HaHV-1-challenged p-aua up-regulated broad classes of genes that contained chitin-binding peritrophin-A domains, which seem to play diverse roles in the p-aua immune response. The p-aua also up-regulated vascular adhesion protein-1 (VAP-1), an important adhesion molecule for lymphocytes, and chitotriosidase-1 (CHIT-1), an immunologically important gene in mammalian immune systems. Moreover, several blood coagulation pathways were dysregulated in the p-aua, possibly indicating viral modulation. We also saw several indications that neurological tissues were specifically affected by HaHV-1, including the dysregulation of beta-1,4-N-acetylgalactosaminyltransferase (B4GALNT), GM2 ganglioside, neuroligin-4 and the Notch signalling pathway. This research may support the development of molecular therapeutics useful to control and/or manage viral outbreaks in abalone culture.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

[The oyster Crassostrea gigas, a new model against cancer]

The Warburg effect is one of the hallmarks of cancer cells in humans. It is a true metabolic reprogramming to aerobic glycolysis, allowing cancer cells to meet their particular energy needs for growth, proliferation, and resistance to apoptosis, depending on the microenvironment they encounter within the tumor. We have recently discovered that the Crassostrea gigas oyster can naturally reprogram its metabolism to the Warburg effect. Thus, the oyster becomes a new invertebrate model useful for cancer research. Due to its lifestyle, the oyster C. gigas has special abilities to adapt its metabolism to the extreme changes in the environment in which it is located. The oyster C. gigas is therefore a model of interest to study how the environment can control the Warburg effect under conditions that could not be explored in vertebrate model species.

Dual Analysis of Virus-Host Interactions: The Case of Ostreid herpesvirus 1 and the Cupped Oyster Crassostrea gigas

Dual analyses of the interactions between Ostreid herpesvirus 1 (OsHV-1) and the bivalve Crassostrea gigas during infection can unveil events critical to the onset and progression of this viral disease and can provide novel strategies for mitigating and preventing oyster mortality. Among the currently used “omics” technologies, dual transcriptomics (dual RNA-seq) coupled with the analysis of viral DNA in the host tissues has greatly advanced the knowledge of genes and pathways mostly contributing to host defense responses, expression profiles of annotated and unknown OsHV-1 open reading frames (ORFs), and viral genome variability. In addition to dual RNA-seq, proteomics and metabolomics analyses have the potential to add complementary information, needed to understand how a malacoherpesvirus can redirect and exploit the vital processes of its host. This review explores our current knowledge of “omics” technologies in the study of host-pathogen interactions and highlights relevant applications of these fields of expertise to the complex case of C gigas infections by OsHV-1, which currently threaten the mollusk production sector worldwide.

Transgenerational plasticity and antiviral immunity in the Pacific oyster (Crassostrea gigas) against Ostreid herpesvirus 1 (OsHV-1)

The oyster's immune system is capable of adapting upon exposure to a pathogen-associated molecular pattern (PAMP) to have an enhanced secondary response against the same type of pathogen. This has been demonstrated using poly(I:C) to elicit an antiviral response in the Pacific oyster (Crassostrea gigas) against Ostreid herpesvirus (OsHV-1). Improved survival following exposure to poly(I:C) has been found in later life stages (within-generational immune priming) and in the next generation (transgenerational immune priming). The mechanism that the oyster uses to transfer immunity to the next generation is unknown. Here we show that oyster larvae have higher survival to OsHV-1 when their mothers, but not their fathers, are exposed to poly(I:C) prior to spawning. RNA-seq provided no evidence to suggest that parental exposure to poly(I:C) reconfigures antiviral gene expression in unchallenged larvae. We conclude that the improved survival of larvae might occur via maternal provisioning of antiviral compounds in the eggs.

Oncogenic virus-induced aerobic glycolysis and tumorigenesis

Enhanced glycolysis under normoxic conditions is known as aerobic glycolysis or the Warburg effect and is a hallmark of many tumors. Viral infection may also induce aerobic glycolysis as it is required for replication and survival. Tumor viruses inducing aerobic glycolysis and lactate production during latent infection suggest a potential role of virus-induced glycolysis in tumorigenesis. Virus or virus-encoded proteins regulate glucose uptake and lactate export, increase the activity of glycolytic enzymes, and modulate glucose metabolic signals. Accumulating evidence suggests that virus-induced glycolysis may facilitate cell growth, transformation, migration, and invasion, but its significance in tumorigenesis remains unclear. We summarize the effects of oncogenic viruses on the metabolic shift to aerobic glycolysis and discuss the possible association of this metabolic reprogramming with tumor development and progression.

Metabolic and immunological responses of male and female new Zealand Greenshell™ mussels (Perna canaliculus) infected with Vibrio sp

Massive mortalities due to pathogens are routinely reported in bivalve cultivation that have significant economic consequences for the global aquaculture industry. However, host-pathogen interactions and infection mechanisms that mediate these interactions are poorly understood. In addition, gender-specific immunological responses have been reported for some species, but the reasons for such differences have not been elucidated. In this study, we used a GC/MS-based metabolomics platform and flow cytometry approach to characterize metabolic and immunological responses in haemolymph of male and female mussels (Perna canaliculus) experimentally infected with Vibrio sp. Sex-based differences in immunological responses were identified, with male mussels displaying higher mortality, oxidative stress and apoptosis after pathogen exposure. However, central metabolic processes appeared to be similar between sexes at 24 h post injection with Vibrio sp. DO1. Significant alterations in relative levels of 37 metabolites were detected between infected and uninfected mussels. These metabolites are involved in major perturbations on the host's innate immune system. In addition, there were alterations of seven metabolites in profiles of mussels sampled on the second day and mussels that survived six days after exposure. These metabolites include itaconic acid, isoleucine, phenylalanine, creatinine, malonic acid, glutaric acid and hydroxyproline. Among these, itaconic acid has the potential to be an important biomarker for Vibrio sp. DO1 infection. These findings provide new insights on the mechanistic relationship between a bivalve host and a pathogenic bacterium and highlight the need to consider host sex as a biological variable in future immunological studies.

Proteomic profile of Ostrea edulis haemolymph in response to bonamiosis and identification of candidate proteins as resistance markers

European flat oyster Ostrea edulis populations have suffered extensive mortalities caused by bonamiosis. The protozoan parasite Bonamia ostreae is largely responsible for this disease in Europe, while its congener B. exitiosa has been detected more recently in various European countries. Both of these intracellular parasites are able to survive and proliferate within haemocytes, the main cellular effectors of the immune system in molluscs. Two-dimensional electrophoresis was used to compare the haemolymph protein profile between Bonamia spp.-infected and non-infected oysters within 3 different stocks, a Galician stock of oysters selected for resistance against bonamiosis, a non-selected Galician stock and a selected Irish stock. Thirty-four proteins with a presumably relevant role in the oyster-Bonamia spp. interaction were identified; they were involved in major metabolic pathways, such as energy production, respiratory chain, oxidative stress, signal transduction, transcription, translation, protein degradation and cell defence. Furthermore, the haemolymph proteomic profiles of the non-infected oysters of the 2 Galician stocks were compared. As a result, 7 proteins representative of the non-infected Galician oysters selected for resistance against bonamiosis were identified; these 7 proteins could be considered as candidate markers of resistance to bonamiosis, which should be further assessed.

Metabolism of the Pacific oyster, Crassostrea gigas, is influenced by salinity and modulates survival to the Ostreid herpesvirus OsHV-1

The Pacific oyster, Crassostrea gigas, is an osmoconforming bivalve exposed to wide salinity fluctuations. The physiological mechanisms used by oysters to cope with salinity stress are energy demanding and may impair other processes, such as defense against pathogens. This oyster species has been experiencing recurrent mortality events caused by the Ostreid herpesvirus 1 (OsHV-1). The objectives of this study were to investigate the effect of salinity (10, 15, 25 and 35‰) on energetic reserves, key enzyme activities and membrane fatty acids, and to identify the metabolic risk factors related to OsHV-1-induced mortality of oysters. Acclimation to low salinity led to increased water content, protein level, and energetic reserves (carbohydrates and triglycerides) of oysters. The latter was consistent with lower activity of hexokinase, the first enzyme involved in glycolysis, up-regulation of AMP-activated protein kinase, a major regulator of cellular energy metabolism, and lower activity of catalase, an antioxidant enzyme involved in management of reactive oxygen species. Acclimation to salinity also involved a major remodeling of membrane fatty acids. Particularly, 20:4n-6 decreased linearly with decreasing salinity, likely reflecting its mobilization for prostaglandin synthesis in oysters. The survival of oysters exposed to OsHV-1 varied from 43% to 96% according to salinity ( Fuhrmann et al., 2016). Risk analyses showed that activity of superoxide dismutase and levels of proteins, carbohydrates, and triglycerides were associated with a reduced risk of death. Therefore, animals with a higher antioxidant activity and a better physiological condition seemed less susceptible to OsHV-1.

The Voltage-Dependent Anion Channel (VDAC) of Pacific Oysters Crassostrea gigas Is Upaccumulated During Infection by the Ostreid Herpesvirus-1 (OsHV-1): an Indicator of the Warburg Effect

Voltage-dependent anion channel (VDAC) is a key mitochondrial protein. VDAC drives cellular energy metabolism by controlling the influx and efflux of metabolites and ions through the mitochondrial membrane, playing a role in its permeabilization. This protein exerts a pivotal role during the white spot syndrome virus (WSSV) infection in shrimp, through its involvement in a particular metabolism that plays in favor of the virus, the Warburg effect. The Warburg effect corresponds to an atypical metabolic shift toward an aerobic glycolysis that provides energy for rapid cell division and resistance to apoptosis. In the Pacific oyster Crassostrea gigas, the Warburg effect occurs during infection by Ostreid herpesvirus (OsHV-1). At present, the role of VDAC in the Warburg effect, OsHV-1 infection and apoptosis is unknown. Here, we developed a specific antibody directed against C. gigas VDAC. This tool allowed us to quantify the tissue-specific expression of VDAC, to detect VDAC oligomers, and to follow the amount of VDAC in oysters deployed in the field. We showed that oysters sensitive to a mortality event in the field presented an accumulation of VDAC. Finally, we propose to use VDAC quantification as a tool to measure the oyster susceptibility to OsHV-1 depending on its environment.

Oyster RNA-seq Data Support the Development of Malacoherpesviridae Genomics

The family of double-stranded DNA (dsDNA) Malacoherpesviridae includes viruses able to infect marine mollusks and detrimental for worldwide aquaculture production. Due to fast-occurring mortality and a lack of permissive cell lines, the available data on the few known Malacoherpesviridae provide only partial support for the study of molecular virus features, life cycle, and evolutionary history. Following thorough data mining of bivalve and gastropod RNA-seq experiments, we used more than five million Malacoherpesviridae reads to improve the annotation of viral genomes and to characterize viral InDels, nucleotide stretches, and SNPs. Both genome and protein domain analyses confirmed the evolutionary diversification and gene uniqueness of known Malacoherpesviridae. However, the presence of Malacoherpesviridae-like sequences integrated within genomes of phylogenetically distant invertebrates indicates broad diffusion of these viruses and indicates the need for confirmatory investigations. The manifest co-occurrence of OsHV-1 genotype variants in single RNA-seq samples of Crassostrea gigas provide further support for the Malacoherpesviridae diversification. In addition to simple sequence motifs inter-punctuating viral ORFs, recombination-inducing sequences were found to be enriched in the OsHV-1 and AbHV1-AUS genomes. Finally, the highly correlated expression of most viral ORFs in multiple oyster samples is consistent with the burst of viral proteins during the lytic phase.

Differential expression of novel metabolic and immunological biomarkers in oysters challenged with a virulent strain of OsHV-1

Early lifestages of the Pacific oyster (Crassostrea gigas) are highly susceptible to infection by OsHV-1 μVar, but little information exists regarding metabolic or pathophysiological responses of larval hosts. Using a metabolomics approach, we identified a range of metabolic and immunological responses in oyster larvae exposed to OsHV-1 μVar; some of which have not previously been reported in molluscs. Multivariate analyses of entire metabolite profiles were able to separate infected from non-infected larvae. Correlation analysis revealed the presence of major perturbations in the underlying biochemical networks and secondary pathway analysis of functionally-related metabolites identified a number of prospective pathways differentially regulated in virus-exposed larvae. These results provide new insights into the pathogenic mechanisms of OsHV-1 infection in oyster larvae, which may be applied to develop disease mitigation strategies and/or as new phenotypic information for selective breeding programmes aiming to enhance viral resistance.

High throughput sequencing of RNA transcriptomes in Ruditapes philippinarum identifies genes involved in osmotic stress response

Ruditapes philippinarum, is an economically important marine bivalve species. The ability to cope with low salinity stress is quite important for the survival of aquatic species under natural conditions. In this study, the transcriptional response of the Manila clam to low salinity stress was characterized using RNA sequencing. The transcriptomes of a low salinity-treatment group (FRp1, FRp2), which survived under low salinity stress, and control group (SRp1, SRp2), which was not subjected to low salinity stress, were sequenced with the Illumina HiSeq platform. A total of 196,578 unigenes were generated. GO and KEGG analyses revealed that signal transduction, immune response, cellular component organization or biogenesis, and energy production processes were the most highly enriched pathways among the genes that were differentially expressed under low salinity stress. All these pathways could be assigned to the following biological functions in the low salinity tolerant Manila clam: signal response to low salinity stress, antioxidant response, intracellular free amino acid transport and metabolism, energy production and conversion, cell signaling pathways, and regulation of ionic content and cell volume. In summary, this is the first study using high-throughput sequencing to identify gene targets that could explain osmotic regulation mechanisms under salinity stress in R. philippinarum.

Viruses infecting marine molluscs

Although a wide range of viruses have been reported in marine molluscs, most of these reports rely on ultrastructural examination and few of these viruses have been fully characterized. The lack of marine mollusc cell lines restricts virus isolation capacities and subsequent characterization works. Our current knowledge is mostly restricted to viruses affecting farmed species such as oysters Crassostrea gigas, abalone Haliotis diversicolor supertexta or the scallop Chlamys farreri. Molecular approaches which are needed to identify virus affiliation have been carried out for a small number of viruses, most of them belonging to the Herpesviridae and birnaviridae families. These last years, the use of New Generation Sequencing approach has allowed increasing the number of sequenced viral genomes and has improved our capacity to investigate the diversity of viruses infecting marine molluscs. This new information has in turn allowed designing more efficient diagnostic tools. Moreover, the development of experimental infection protocols has answered some questions regarding the pathogenesis of these viruses and their interactions with their hosts. Control and management of viral diseases in molluscs mostly involve active surveillance, implementation of effective bio security measures and development of breeding programs. However factors triggering pathogen development and the life cycle and status of the viruses outside their mollusc hosts still need further investigations.

Proteomic analysis of the adaptative response of Mucor spp. to cheese environment

In the cheese industry context, Mucor species exhibit an ambivalent behavior as some species are essential "technological" organisms of some cheeses while others can be spoiling agents. Previously, we observed that cheese "technological" species exhibited higher optimal growth rates on cheese related matrices than on synthetic media. This growth pattern combined with morphological differences raise the question of their adaptation to cheese. In this study, using a comparative proteomic approach, we described the metabolic pathways of three Mucor strains considered as "technological" or "contaminant" in the cheese environment (M. lanceolatus UBOCC-A-109153, M. racemosus UBOCC-A-109155, M. circinelloides CBS 277-49) as well as a non-cheese related strain (M. endophyticus CBS 385-95). Overall, 15.8 to 19.0% of the proteomes showed a fold change ≥1.6 in Potato Dextrose Agar (PDA) versus Cheese Agar (CA), a cheese mimicking-medium. The 289 differentially expressed proteins identified by LC MS-MS analysis were mostly assigned to energy and amino-acid metabolisms in PDA whereas a higher diversity of biological processes was observed for cheese related strains in CA. Surprisingly, the vast majority (72.9%) of the over-accumulated proteins were different according to the considered medium and strain. These results strongly suggest that the observed better adaptative response of "technological" strains to cheese environment is mediated by species-specific proteins.

BIOLOGICAL SIGNIFICANCE

The Mucor genus consists of a multitude of poorly known species. In the food context, few species are known for their positive role in the production of various food products, including cheese, while others are spoiling agents. The present study focused on the analysis of morphological and proteome differences of various Mucor spp. representative strains known as either positively (hereafter referred as "technological") or negatively (hereafter referred as "contaminant") associated with cheese or non-related to cheese (endophyte) on two different media, a synthetic medium and a cheese-mimicking medium. The main goal was to assess if adaptative traits of "technological" strains to the cheese environment could be identified. This work was based on observations we did in a recently published physiological study (Morin-Sardin et al., 2016). One of the important innovative aspects lies in the use for the first time of an extensive 2-DE approach to compare proteome variations for 4 strains on two different media. Results obtained offered an insight in the metabolic mechanisms associated with growth on a given medium and showed that adaptation to cheese environment is probably supported by species-specific proteins. The obtained data represent an essential step point for more targeted studies at the genomic and transcriptomic levels.

Identification of a regulation network in response to cadmium toxicity using blood clam Tegillarca granosa as model

Clam, a filter-feeding lamellibranch mollusk, is capable to accumulate high levels of trace metals and has therefore become a model for investigation the mechanism of heavy metal toxification. In this study, the effects of cadmium were characterized in the gills of Tegillarca granosa during a 96-hour exposure course using integrated metabolomic and proteomic approaches. Neurotoxicity and disturbances in energy metabolism were implicated according to the metabolic responses after Cd exposure, and eventually affected the osmotic function of gill tissue. Proteomic analysis showed that oxidative stress, calcium-binding and sulfur-compound metabolism proteins were key factors responding to Cd challenge. A knowledge-based network regulation model was constructed with both metabolic and proteomic data. The model suggests that Cd stimulation mainly inhibits a core regulation network that is associated with histone function, ribosome processing and tight junctions, with the hub proteins actin, gamma 1 and Calmodulin 1. Moreover, myosin complex inhibition causes abnormal tight junctions and is linked to the irregular synthesis of amino acids. For the first time, this study provides insight into the proteomic and metabolomic changes caused by Cd in the blood clam T. granosa and suggests a potential toxicological pathway for Cd.

Reprint of 'Diseases in marine invertebrates associated with mariculture and commercial fisheries'

Diseases in marine invertebrates are increasing in both frequency and intensity around the globe. Diseases in individuals which offer some commercial value are often well documented and subsequently well studied in comparison to those wild groups offering little commercial gain. This is particularly the case with those associated with mariculture or the commercial fisheries. Specifically, these include many Holothuroidea, and numerous crustacea and mollusca species. Pathogens/parasites consisting of both prokaryotes and eukaryotes from all groups have been associated with diseases from such organisms, including bacteria, viruses, fungi and protozoa. Viral pathogens in particular, appear to be an increasingly important group and research into this group will likely highlight a larger number of diseases and pathogens being described in the near future. Interestingly, although there are countless examples of the spread of disease usually associated with transportation of specific infected hosts for development of aquaculture practices, this process appears to be continuing with no real sign of effective management and mitigation strategies being implicated. Notably, even in well developed countries such as the UK and the US, even though live animal trade may be well managed, the transport of frozen food appears to be less well so and as evidence suggests, even these to have the potential to transmit pathogens when used as a food source for example.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

Non-additive effects of ocean acidification in combination with warming on the larval proteome of the Pacific oyster, Crassostrea gigas

Increasing atmospheric carbon dioxide results in ocean acidification and warming, significantly impacting marine invertebrate larvae development. We investigated how ocean acidification in combination with warming affected D-veliger larvae of the Pacific oyster Crassostrea gigas. Larvae were reared for 40h under either control (pH8.1, 20 °C), acidified (pH7.9, 20 °C), warm (pH8.1, 22 °C) or warm acidified (pH7.9, 22 °C) conditions. Larvae in acidified conditions were significantly smaller than in the control, but warm acidified conditions mitigated negative effects on size, and increased calcification. A proteomic approach employing two-dimensional electrophoresis (2-DE) was used to quantify proteins and relate their abundance to phenotypic traits. In total 12 differentially abundant spots were identified by nano-liquid chromatography-tandem mass spectrometry. These proteins had roles in metabolism, intra- and extra-cellular matrix formations, stress response, and as molecular chaperones. Seven spots responded to reduced pH, four to increased temperature, and six to acidification and warming. Reduced abundance of proteins such as ATP synthase and GAPDH, and increased abundance of superoxide dismutase, occurred when both pH and temperature changes were imposed, suggesting altered metabolism and enhanced oxidative stress. These results identify key proteins that may be involved in the acclimation of C. gigas larvae to ocean acidification and warming.

SIGNIFICANCE

Increasing atmospheric CO2 raises sea surface temperatures and results in ocean acidification, two climatic variables known to impact marine organisms. Larvae of calcifying species may be particularly at risk to such changing environmental conditions. The Pacific oyster Crassostrea gigas is ecologically and commercially important, and understanding its ability to acclimate to climate change will help to predict how aquaculture of this species is likely to be impacted. Modest, yet realistic changes in pH and/or temperature may be more informative of how populations will respond to contemporary climate change. We showed that concurrent acidification and warming mitigates the negative effects of pH alone on size of larvae, but proteomic analysis reveals altered patterns of metabolism and an increase in oxidative stress suggesting non-additive effects of the interaction between pH and temperature on protein abundance. Thus, even small changes in climate may influence development, with potential consequences later in life.

Molecular Basis for Adaptation of Oysters to Stressful Marine Intertidal Environments

Oysters that occupy estuarine and intertidal habitats have well-developed stress tolerance mechanisms to tolerate harsh and dynamically changing environments. In this review, we summarize common pathways and genomic features in oyster that are responsive to environmental stressors such as temperature, salinity, hypoxia, air exposure, pathogens, and anthropogenic pollutions. We first introduce the key genes involved in several pathways, which constitute the molecular basis for adaptation to stress. We use genome analysis to highlight the strong cellular homeostasis system, a unique adaptive characteristic of oysters. Next, we provide a global view of features of the oyster genome that contribute to stress adaptation, including oyster-specific gene expansion, highly inducible expression, and functional divergence. Finally, we review the consequences of interactions between oysters and the environment from ecological and evolutionary perspectives by discussing mass mortality and adaptive divergence among populations and related species of the genus Crassostrea. We conclude with prospects for future study.

Diseases in marine invertebrates associated with mariculture and commercial fisheries

Diseases in marine invertebrates are increasing in both frequency and intensity around the globe. Diseases in individuals which offer some commercial value are often well documented and subsequently well studied in comparison to those wild groups offering little commercial gain. This is particularly the case with those associated with mariculture or the commercial fisheries. Specifically, these include many Holothuroidea, and numerous crustacea and mollusca species. Pathogens/parasites consisting of both prokaryotes and eukaryotes from all groups have been associated with diseases from such organisms, including bacteria, viruses, fungi and protozoa. Viral pathogens in particular, appear to be an increasingly important group and research into this group will likely highlight a larger number of diseases and pathogens being described in the near future. Interestingly, although there are countless examples of the spread of disease usually associated with transportation of specific infected hosts for development of aquaculture practices, this process appears to be continuing with no real sign of effective management and mitigation strategies being implicated. Notably, even in well developed countries such as the UK and the US, even though live animal trade may be well managed, the transport of frozen food appears to be less well so and as evidence suggests, even these to have the potential to transmit pathogens when used as a food source for example.

Immune and stress responses in oysters with insights on adaptation

Oysters are representative bivalve molluscs that are widely distributed in world oceans. As successful colonizers of estuaries and intertidal zones, oysters are remarkably resilient against harsh environmental conditions including wide fluctuations in temperature and salinity as well as prolonged air exposure. Oysters have no adaptive immunity but can thrive in microbe-rich estuaries as filter-feeders. These unique adaptations make oysters interesting models to study the evolution of host-defense systems. Recent advances in genomic studies including sequencing of the oyster genome have provided insights into oyster's immune and stress responses underlying their amazing resilience. Studies show that the oyster genomes are highly polymorphic and complex, which may be key to their resilience. The oyster genome has a large gene repertoire that is enriched for immune and stress response genes. Thousands of genes are involved in oyster's immune and stress responses, through complex interactions, with many gene families expanded showing high sequence, structural and functional diversity. The high diversity of immune receptors and effectors may provide oysters with enhanced specificity in immune recognition and response to cope with diverse pathogens in the absence of adaptive immunity. Some members of expanded immune gene families have diverged to function at different temperatures and salinities or assumed new roles in abiotic stress response. Most canonical innate immunity pathways are conserved in oysters and supported by a large number of diverse and often novel genes. The great diversity in immune and stress response genes exhibited by expanded gene families as well as high sequence and structural polymorphisms may be central to oyster's adaptation to highly stressful and widely changing environments.

Transcriptome analysis reveals strong and complex antiviral response in a mollusc

Viruses are highly abundant in the oceans, and how filter-feeding molluscs without adaptive immunity defend themselves against viruses is not well understood. We studied the response of a mollusc Crassostrea gigas to Ostreid herpesvirus 1 µVar (OsHV-1μVar) infections using transcriptome sequencing. OsHV-1μVar can replicate extremely rapidly after challenge of C. gigas as evidenced by explosive viral transcription and DNA synthesis, which peaked at 24 and 48 h post-inoculation, respectively, accompanied by heavy oyster mortalities. At 120 h post-injection, however, viral gene transcription and DNA load, and oyster mortality, were greatly reduced indicating an end of active infections and effective control of viral replication in surviving oysters. Transcriptome analysis of the host revealed strong and complex responses involving the activation of all major innate immune pathways that are equipped with expanded and often novel receptors and adaptors. Novel Toll-like receptor (TLR) and MyD88-like genes lacking essential domains were highly up-regulated in the oyster, possibly interfering with TLR signal transduction. RIG-1/MDA5 receptors for viral RNA, interferon-regulatory factors, tissue necrosis factors and interleukin-17 were highly activated and likely central to the oyster's antiviral response. Genes related to anti-apoptosis, oxidation, RNA and protein destruction were also highly up-regulated, while genes related to anti-oxidation were down-regulated. The oxidative burst induced by the up-regulation of oxidases and severe down-regulation of anti-oxidant genes may be important for the destruction of viral components, but may also exacerbate oyster mortality. This study provides unprecedented insights into antiviral response in a mollusc. The mobilization and complex regulation of expanded innate immune-gene families highlights the oyster genome's adaptation to a virus-rich marine environment.

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.

Proteomic basis of stress responses in the gills of the pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is one of the dominant sessile inhabitants of the estuarine intertidal zone, which is a physically harsh environment due to the presence of a number of stressors. Oysters have adapted to highly dynamic and stressful environments, but the molecular mechanisms underlying such stress adaptation are largely unknown. In the present study, we examined the proteomic responses in the gills of C. gigas exposed to three stressors (high temperature, low salinity, and aerial exposure) they often encounter in the field. We quantitatively compared the gill proteome profiles using iTRAQ-coupled 2-D LC-MS/MS. There were 3165 identified proteins among which 2379 proteins could be quantified. Heat shock, hyposalinity, and aerial exposure resulted in 50, 15, and 33 differentially expressed gill proteins, respectively. Venn diagram analysis revealed substantial different responses to the three stressors. Only xanthine dehydrogenase/oxidase showed a similar expression pattern across the three stress treatments, suggesting that reduction of ROS accumulation may be a conserved response to these stressors. Heat shock caused significant overexpression of molecular chaperones and production of S-adenosyl-l-methionine, indicating their crucial protective roles against protein denature. In addition, heat shock also activated immune responses, Ca(2+) binding protein expression. By contrast, hyposalinity and aerial exposure resulted in the up-regulation of 3-demethylubiquinone-9 3-methyltransferase, indicating that increase in ubiquinone synthesis may contribute to withstanding both the osmotic and desiccation stress. Strikingly, the majority of desiccation-responsive proteins, including those involved in metabolism, ion transportation, immune responses, DNA duplication, and protein synthesis, were down-regulated, indicating conservation of energy as an important strategy to cope with desiccation stress. There was a high consistency between the expression levels determined by iTRAQ and Western blotting, highlighting the high reproducibility of our proteomic approach and its great value in revealing molecular mechanisms of stress responses.