Suppression substractive hybridisation (SSH) and real time PCR reveal differential gene expression in the Pacific cupped oyster, Crassostrea gigas, challenged with Ostreid herpesvirus 1

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log₁₀ (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log₁₀ (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log₁₀ (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log₁₀ (Foxl2) and Log₁₀ (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

Voltage-dependent anion channel 2 (VDAC2) facilitates the accumulation of rice stripe virus in the vector Laodelphax striatellus

Rice stripe virus (RSV) causes enormous losses in rice production and is transmitted by the small brown planthopper, Laodelphax striatellus, in a persistent-propagative manner. RSV accumulation within the gut lumen of the vector is indispensable for the successful transmission to rice and insects. In this study, we obtained a 1464 bp full-length cDNA of a voltage-dependent anion channel 2 from L. striatellus (LsVDAC2), which encodes a 283 amino acid protein. RSV infection increased the expression of LsVDAC2 in the midguts and ovaries of L. striatellus by 260% and 228%, respectively. Silencing of LsVDAC2 resulted in a 88% reduction of RSV loads at 24 h after RNAi, indicating that LsVDAC2 facilitates RSV accumulation in the vector. Yeast two-hybrid and GST pulldown assays demonstrated that LsVDAC2 interacted with RSV RNA-dependent RNA polymerase, RdRp. Furthermore, experiments in vivo and in vitro showed that LsVDAC2 induced the apoptotic response in RSV-infected insects and tissues. Silencing of LsVDAC2 via RNAi significantly reduced the expression of genes for apoptosis-related caspases 1a and 1c by 62% and 78%, respectively, in RSV-infected vectors. Whether LsVDAC2-induced RSV accumulation is related to RSV RdRp and LsVDAC2-induced cell apoptosis deserves further investigation.

Evaluation of tangential flow filtration coupled to long-read sequencing for ostreid herpesvirus type 1 genome assembly

Whole-genome sequencing is widely used to better understand the transmission dynamics, the evolution and the emergence of new variants of viral pathogens. This can bring crucial information to stakeholders for disease management. Unfortunately, aquatic virus genomes are usually difficult to characterize because most of these viruses cannot be easily propagated in vitro. Developing methodologies for routine genome sequencing of aquatic viruses is timely given the ongoing threat of disease emergence. This is particularly true for pathogenic viruses infecting species of commercial interest that are widely exchanged between production basins or countries. For example, the ostreid herpesvirus type 1 (OsHV-1) is a Herpesvirus widely associated with mass mortality events of juvenile Pacific oyster Crassostrea gigas. Genomes of Herpesviruses are large and complex with long direct and inverted terminal repeats. In addition, OsHV-1 is unculturable. It therefore accumulates several features that make its genome sequencing and assembly challenging. To overcome these difficulties, we developed a tangential flow filtration (TFF) method to enrich OsHV-1 infective particles from infected host tissues. This virus purification allowed us to extract high molecular weight and high-quality viral DNA that was subjected to Illumina short-read and Nanopore long-read sequencing. Dedicated bioinformatic pipelines were developed to assemble complete OsHV-1 genomes with reads from both sequencing technologies. Nanopore sequencing allowed characterization of new structural variations and major viral isomers while having 99,98 % of nucleotide identity with the Illumina assembled genome. Our study shows that TFF-based purification method, coupled with Nanopore sequencing, is a promising approach to enable in field sequencing of unculturable aquatic DNA virus.

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.

A nuclear receptor heterodimer, CgPPAR2-CgRXR, acts as a regulator of carotenoid metabolism in Crassostrea gigas

Nuclear receptors (NRs) are mostly ligand-activated transcription factors in animals and play essential roles in metabolism and homeostasis. The NR heterodimer composed of PPAR/RXR (peroxisome proliferator-activated receptor/retinoid X receptor) is considered a key regulator of lipid metabolism in vertebrate. However, in molluscs, how this heterodimer is involved in carotenoid metabolism remains unclear. To elucidate how this heterodimer regulates carotenoid metabolism, we identified a PPAR gene in C. gigas, designated as CgPPAR2 (LOC105323212), and functionally characterized it using two-hybrid and reporter systems. CgPPAR2 is a direct orthologue of vertebrate PPARs and the second PPAR gene identified in C. gigas genome in addition to CgPPAR1 (LOC105317849). The results demonstrated that CgPPAR2 protein can form heterodimer with C. gigas RXR (CgRXR), and then regulate carotenoid metabolism by controlling carotenoid cleavage oxygenases with different carotenoid cleavage efficiencies. This regulation can be affected by retinoid ligands, i.e., carotenoid derivatives, validating a negative feedback regulation mechanism of carotenoid cleavage for retinoid production. Besides, organotins may disrupt this regulatory process through the mediation of CgPPAR2/CgRXR heterodimer. This is the first report of PPAR/RXR heterodimer regulating carotenoid metabolism in mollusks, contributing to a better understanding of the evolution and conservation of this nuclear receptor heterodimer.

Transcriptome analysis based on dietary beta-carotene supplement reveals genes potentially involved in carotenoid metabolism in Crassostrea gigas

Carotenoids are essential micronutrients for animals, and they can only be obtained from the diet for mollusk as well as other animals. In the body, carotenoids undergo processes including absorption, transport, deposition, and metabolic conversion; however, knowledge of the involved genes is still limited. To elucidate the molecular mechanisms of carotenoid processing and identify the related genes in Pacific oyster (Crassostrea gigas), we performed a comparative transcriptome analysis using digestive gland tissues of oysters on a beta-carotene supplemented diet or a normal diet. A total of 718 differentially expressed genes were obtained, including 505 upregulated and 213 downregulated genes in the beta-carotene supplemented group. Function Annotation and enrichment analyses revealed enrichment in genes possibly involved in carotenoid transport and storage (e.g., LOC105342035), carotenoid cleavage (e.g., LOC105341121), retinoid homeostasis (e.g., LOC105339597) and PPAR signaling pathway (e.g., LOC105323212). Notably, down-regulation of mRNA expressions of two apolipoprotein genes (LOC105342035 and LOC105342186) by RNA interference significantly decreased the carotenoid level in the digestive gland, supporting their role in carotenoid transport and storage. Based on these differentially expressed genes, we propose that there may be a negative feedback mechanism regulated by nuclear receptor transcription factors controlling carotenoid oxygenases. Our findings provide useful hints for elucidating the molecular basis of carotenoid metabolism and functions of carotenoid-related genes in the oyster.

Monitoring Autophagy at Cellular and Molecular Level in Crassostrea gigas During an Experimental Ostreid Herpesvirus 1 (OsHV-1) Infection

Mortality outbreaks of young Pacific oysters, Crassostrea gigas, have seriously affected the oyster-farming economy in several countries around the world. Although the causes of these mortality outbreaks appear complex, a viral agent has been identified as the main factor: a herpesvirus called ostreid herpesvirus 1 (OsHV-1). Autophagy is an important degradation pathway involved in the response to several pathologies including viral diseases. In C. gigas, recent studies indicate that this pathway is conserved and functional in at least haemocytes and the mantle. Furthermore, an experimental infection in combination with compounds known to inhibit or induce autophagy in mammals revealed that autophagy is involved in the response to OsHV-1 infection. In light of these results, the aim of this study was to determine the role of autophagy in the response of the Pacific oyster to infection by virus OsHV-1. For this purpose, an experimental infection in combination with a modulator of autophagy was performed on Pacific oysters known to have intermediate susceptibility to OsHV-1 infection. In haemolymph and the mantle, the autophagy response was monitored by flow cytometry, western blotting, and real-time PCR. At the same time, viral infection was evaluated by quantifying viral DNA and RNA amounts by real-time PCR. Although the results showed activation of autophagy in haemolymph and the mantle 14 hours post infection (after viral replication was initiated), they were also indicative of different regulatory mechanisms of autophagy in the two tissues, thus supporting an important function of autophagy in the response to virus OsHV-1.

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.

Identification and characterization of an apoptosis-inducing factor 1 involved in apoptosis and immune defense of oyster, Crassostrea gigas

The apoptosis-inducing factor (AIF) is a phylogenetically old protein with classic function of inducing caspase-independent apoptosis, which extensively present in all primary kingdoms. In the present study, an AIF homologue (designated as CgAIF1) was identified from oyster Crassostrea gigas. The open reading frame of CgAIF1 cDNA was of 1836 bp encoding a peptide of 611 amino acid residues. There are a Pyr_redox_2 domain and an AIF_C domain in the predicted CgAIF1 protein. The deduced amino acid sequence of CgAIF1 shared 35.44%-79.22% similarity with AIF1s from other species. In the phylogenetic tree, CgAIF1 firstly clustered with mollusc AIF1s, and then with insect AIF1s, displaying separation from vertebrate AIF1s. The mRNA transcripts of CgAIF1 were constitutively distributed in all the tested oyster tissues, with the highest level in gills (12.98-fold of that in haemocytes, p < 0.05). After LPS and Poly (I:C) stimulation, the mRNA transcripts of CgAIF1 in gills were significantly increased at 6 h and 24 h (5.79-fold, p < 0.001, and 21.96-fold compared to the control group, p < 0.05), respectively. In immunocytochemical assay, the CgAIF1 positive signals were mainly distributed in the cytoplasm of haemocytes, while after Poly (I:C) stimulation, the increased CgAIF1 positive signals were observed in the nucleus. Moreover, in the HEK293T cells transfected with pcDNA3.1-CgAIF1 recombinant plasmid, green signal of CgAIF1 were observed in both the cytoplasm and nucleus. The cell mortality rate, cell shrinking and the phosphatidylserine (PS) ectropion (Annexin V⁺/PI^- cells and Annexin V⁺/PI⁺ cells) of CgAIF1 transfected HEK293T cells were significantly increased, compared to the groups with or without pcDNA3.1 transfection. These results collectively suggested that CgAIF1 was a conserved AIF1 member in oysters, and participated in immune response by inducing cell apoptosis.

Viral Decoys: The Only Two Herpesviruses Infecting Invertebrates Evolved Different Transcriptional Strategies to Deflect Post-Transcriptional Editing

The highly versatile group of Herpesviruses cause disease in a wide range of hosts. In invertebrates, only two herpesviruses are known: the malacoherpesviruses HaHV-1 and OsHV-1 infecting gastropods and bivalves, respectively. To understand viral transcript architecture and diversity we first reconstructed full-length viral genomes of HaHV-1 infecting Haliotis diversicolor supertexta and OsHV-1 infecting Scapharca broughtonii by DNA-seq. We then used RNA-seq over the time-course of experimental infections to establish viral transcriptional dynamics, followed by PacBio long-read sequencing of full-length transcripts to untangle viral transcript architectures at two selected time points. Despite similarities in genome structure, in the number of genes and in the diverse transcriptomic architectures, we measured a ten-fold higher transcript variability in HaHV-1, with more extended antisense gene transcription. Transcriptional dynamics also appeared different, both in timing and expression trends. Both viruses were heavily affected by post-transcriptional modifications performed by ADAR1 affecting sense-antisense gene pairs forming dsRNAs. However, OsHV-1 concentrated these modifications in a few genomic hotspots, whereas HaHV-1 diluted ADAR1 impact by elongated and polycistronic transcripts distributed over its whole genome. These transcriptional strategies might thus provide alternative potential roles for sense-antisense transcription in viral transcriptomes to evade the host's immune response in different virus-host combinations.

Establishment of primary cell culture of Ruditapes decussatus haemocytes for metal toxicity assessment

In ecotoxicology, in vitro testing on cell cultures represents an ideal alternative to in vivo strategies for emerging contaminants. These tests have limited use particularly with marine invertebrates like the clams Ruditapes decussatus. In the present study, a primary culture of R. decussatus haemocytes was realized for the first time in order to determine the effect of metals (copper, zinc, and cobalt) on haemocyte parameters like viability and phagocytosis. Results showed that (i) among the studied medium, the modified Leibovitz (L-15) is the best for R. decussatus haemocytes primary culture. (ii) The primary culture system used here represents a suitable in vitro model for assessing cytotoxic responses, (iii) a decrease of cell viability and phagocytosis after 24 h exposure to 100 μg mL^-1 CoSO₄ and an increase of phagocytosis after 24 h exposure to 50 μg mL^-1CuSO₄.

First detection of OsHV-1 in the cephalopod Octopus vulgaris. Is the octopus a dead-end for OsHV-1?

The ostreid herpes virus (OsHV-1), associated with massive mortalities in the bivalve Crassostrea gigas, was detected for the first time in the cephalopod Octopus vulgaris. Wild adult animals from a natural breeding area in Spain showed an overall prevalence of detection of 87.5% between 2010 and 2015 suggesting an environmental source of viral material uptake. Overall positive PCR detections were significantly higher in adult animals (p = 0.031) compared to newly hatched paralarvae (62%). Prevalence in embryos reached 65%. Sequencing of positive amplicons revealed a match with the variant OsHV-1 µVar showing the genomic features that distinguish this variant in the ORF4. Gill tissues from adult animals were also processed for in situ hybridization and revealed positive labelling. Experimental exposure trials in octopus paralarvae were carried out by cohabitation with virus injected oysters and by immersion in viral suspension observing a significant decrease in paralarval survival in both experiments. An increase in the number of OsHV-1 positive animals was detected in dead paralarvae after cohabitation with virus injected oysters. No signs of viral replication were observed based on lack of viral gene expression or visualization of viral structures by transmission electron microscopy. The octopus response against OsHV-1 was evaluated by gene expression of previously reported transcripts involved in immune response in C. gigas suggesting that immune defences in octopus are also activated after exposure to OsHV-1.

Freshwater mussels (Unionidae) brought into captivity exhibit up-regulation of genes involved in stress and energy metabolism

Approximately two thirds of freshwater mussel species in the United States and Canada are imperiled, and populations are declining rapidly. Translocation and captive management are commonly used to mitigate losses of freshwater mussel biodiversity, but these conservation tools may result in decreased growth and increased mortality. This study uses RNA-Seq to determine how translocation into captivity affects gene expression in Amblema plicata. Mussels were collected from the Muskingum River in Ohio, USA and brought into a captive holding facility. RNA was extracted from gill tissue 11 months post translocation from mussels in captivity and the Muskingum River on the same day. RNA was sequenced on an Illumina HiSeq 2500, and differential expression analysis was performed on de novo assembled transcripts. More than 1200 transcripts were up-regulated in captive mussels, and 246 were assigned functional annotations. Many up-regulated transcripts were involved in energy metabolism and the stress response, such as heat shock proteins and antioxidants. More than 500 transcripts were down-regulated in captive mussels, and 41 were assigned functional annotations. We observed an over-representation of down-regulated transcripts associated with immune response. Our work suggests that A. plicata experienced moderate levels of stress and altered energy metabolism and immune response for at least 11 months post translocation into captivity.

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.

Functional characterization and gene expression profile of perforin-2 in starry flounder (Platichthys stellatus)

The membrane attack complex/perforin (MACPF) superfamily consists of multifunctional proteins that form pores on the membrane surface of microorganisms to induce their death and have various immune-related functions. PFN2 is a perforin-like protein with an MACPF domain, and humans with deficient PFN2 levels have increased susceptibility to bacterial infection, which can lead to fatal consequences for some patients. Therefore, in this study, we confirmed the antimicrobial function of PFN2 in starry flounder (Platichthys stellatus). The molecular properties were confirmed based on the verified amino acid sequence of PsPFN2. In addition, the expression characteristics of tissue-specific and pathogen-specific PsPFN2 mRNA were also confirmed. The recombinant protein was produced using Escherichia coli, and the antimicrobial activity was then confirmed. The coding sequence of PFN2 (PsPFN2) in P. stellatus consists of 710 residues. The MACPF domain was conserved throughout evolution, as shown by multiple sequence alignment and phylogenetic analysis. PsPFN2 mRNA is abundantly distributed in immune-related organs such as the spleen and gills of healthy starry flounder, and significant expression changes were confirmed after artificial infection by bacteria or viruses. We cloned the MACPF domain region of PFN2 to produce a recombinant protein (rPFN2) and confirmed its antibacterial effect against a wide range of bacterial species and the parasite (Miamiensis avidus).

The neonicotinoid insecticide imidacloprid, but not salinity, impacts the immune system of Sydney rock oyster, Saccostrea glomerata

The broad utilisation of neonicotinoids, particularly imidacloprid (IMI), in agriculture has led to unplanned contamination of aquatic systems around the world. The sublethal effects of individual pesticides on the immune system of oysters, as well as their combined effects with other environmental stressors that fluctuate in estuarine environments, such as salinity, are yet to be investigated in ecotoxicology. We investigated the acute (4 d) toxicity of IMI in two salinity regimes on the immune parameters of Sydney rock oysters (SRO), including total hemocyte counts (THC), differential hemocyte counts (DHC), phagocytosis and hemocyte aggregation (HA), hemolymph protein expression and enzyme (catalase (CAT), glutathione S-transferase (GST) and acetylcholinesterase (AChE)) activities. Environmentally relevant concentrations of IMI were found to cause an increase in THC, induce GST activity, reduce HA, and inhibit AChE activity. However, DHC, CAT activity and phagocytosis were not significantly impacted at any test concentration at either salinity. IMI concentrations ≥0.01 mg/L significantly altered the expression of 28 proteins in the hemolymph of SRO, including an increase in the relative expression of extracellular superoxide dismutase, severin, ATP synthase subunit beta, as well as stress response proteins (heat shock proteins, serine/threonine-protein kinase DCLK3 and peroxiredoxin-1), and a decrease/absence of collagen alpha-4 (VI) and alpha-6 (VI) chain, metalloendopeptidase, L-ascorbate oxidase, transporter, CEP209_CC5 domain-containing protein and actin. This study indicates that the immune system of SRO can be impacted at environmentally relevant concentrations of IMI, but reduced salinity does not appear to influence the toxicity of this insecticide.

Ancient but Not Forgotten: New Insights Into MPEG1, a Macrophage Perforin-Like Immune Effector

Macrophage-expressed gene 1 [MPEG1/Perforin-2 (PRF2)] is an ancient metazoan protein belonging to the Membrane Attack Complex/Perforin (MACPF) branch of the MACPF/Cholesterol Dependent Cytolysin (CDC) superfamily of pore-forming proteins (PFPs). MACPF/CDC proteins are a large and extremely diverse superfamily that forms large transmembrane aqueous channels in target membranes. In humans, MACPFs have known roles in immunity and development. Like perforin (PRF) and the membrane attack complex (MAC), MPEG1 is also postulated to perform a role in immunity. Indeed, bioinformatic studies suggest that gene duplications of MPEG1 likely gave rise to PRF and MAC components. Studies reveal partial or complete loss of MPEG1 causes an increased susceptibility to microbial infection in both cells and animals. To this end, MPEG1 expression is upregulated in response to proinflammatory signals such as tumor necrosis factor α (TNFα) and lipopolysaccharides (LPS). Furthermore, germline mutations in MPEG1 have been identified in connection with recurrent pulmonary mycobacterial infections in humans. Structural studies on MPEG1 revealed that it can form oligomeric pre-pores and pores. Strikingly, the unusual domain arrangement within the MPEG1 architecture suggests a novel mechanism of pore formation that may have evolved to guard against unwanted lysis of the host cell. Collectively, the available data suggest that MPEG1 likely functions as an intracellular pore-forming immune effector. Herein, we review the current understanding of MPEG1 evolution, regulation, and function. Furthermore, recent structural studies of MPEG1 are discussed, including the proposed mechanisms of action for MPEG1 bactericidal activity. Lastly limitations, outstanding questions, and implications of MPEG1 models are explored in the context of the broader literature and in light of newly available structural data.

Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in Crassostrea gigas infected by Ostreid herpesvirus 1

Background

Since 2008, the aquaculture production of Crassostrea gigas was heavily affected by mass mortalities associated to Ostreid herpesvirus 1 (OsHV-1) microvariants worldwide. Transcriptomic studies revealed the major antiviral pathways of the oyster immune response while other findings suggested that also small non-coding RNAs (sncRNA) such as microRNAs might act as key regulators of the oyster response against OsHV-1. To explore the explicit connection between small non-coding and protein-coding transcripts, we performed paired whole transcriptome analysis of sncRNA and messenger RNA (mRNA) in six oysters selected for different intensities of OsHV-1 infection.

Results

The mRNA profiles of the naturally infected oysters were mostly governed by the transcriptional activity of OsHV-1, with several differentially expressed genes mapping to the interferon, toll, apoptosis, and pro-PO pathways. In contrast, miRNA profiles suggested more complex regulatory mechanisms, with 15 differentially expressed miRNAs (DE-miRNA) pointing to a possible modulation of the host response during OsHV-1 infection. We predicted 68 interactions between DE-miRNAs and oyster 3′-UTRs, but only few of them involved antiviral genes. The sncRNA reads assigned to OsHV-1 rather resembled mRNA degradation products, suggesting the absence of genuine viral miRNAs.

Conclusions

We provided data describing the miRNAome during OsHV-1 infection in C. gigas. This information can be used to understand the role of miRNAs in healthy and diseased oysters, to identify new targets for functional studies and, eventually to disentangle cause and effect relationships during viral infections in marine mollusks.

A novel Adiponectin receptor (AdipoR) involved in regulating cytokines production and apoptosis of haemocytes in oyster Crassostrea gigas

Adiponectin receptors (AdipoRs) comprise a seven-transmembrane domain-containing protein family, which specifically recognize adiponectin (APN) and play critical roles in the immunological and physiological processes in vertebrates. In the present study, a novel AdipoR is identified from oyster Crassostrea gigas (designated as CgAdipoR). The full-length cDNA of CgAdipoR is of 1209 bp encoding a polypeptide of 343 amino acids. There is an N-terminal domain, a Hly III domain, and a C-terminal domain in CgAdipoR. After the transfection of CgAdipoR, the level of intracellular Ca²⁺ into HEK293T cells increases significantly (1.36-fold, p < 0.05) after APN incubation. The mRNA transcripts of CgAdipoR are widely distributed in all the tested tissues, with the highest expression level in haemocytes (3.20-fold of that in hepatopancreas, p < 0.05). After lipopolysaccharide (LPS), Vibrio splendidus and polyinosinic-polycytidylic acid (poly (I:C)) stimulations, the mRNA expression of CgAdipoR in haemocytes is significantly up-regulated and reached the highest level at 24 h (15.07-fold, p < 0.01), 6 h (4.39-fold, p < 0.01) and 24 h (5.62-fold, p < 0.01) compared to control group, respectively. After CgAdipoR is interfered by specific CgAdipoR-dsRNA, the expression level of interleukins (CgIL17-1, CgIL17-2, CgIL17-3 and CgIL17-5) in haemocytes decreases significantly (p < 0.01) at 24 h post LPS stimulation, while the expression level of CgTNF-1 increases significantly (1.68-fold, p < 0.01), compared to that in the dsEGFP group. In CgAdipoR dsRNA-injected oysters, the mRNA expressions of anti-apoptotic B-cell lymphoma-2 (Bcl-2) in haemocytes significantly decreases at 24 h after LPS challenge, which is (0.58-fold, p < 0.05) of that in dsEGFP-injected oysters, while the apoptotic rate of haemocytes is significantly up-regulated (1.93-fold of that in dsEGFP group, p < 0.05). These results collectively suggest that CgAdipoR plays an important role in the immune response of oysters by regulating the expressions of inflammatory cytokines and haemocyte apoptosis.

Identification of candidate microRNAs from Ostreid herpesvirus-1 (OsHV-1) and their potential role in the infection of Pacific oysters (Crassostrea gigas)

Oyster production is an economic activity of great interest worldwide. Recently, oysters have been suffering significant mortalities from OsHV-1infection, which has resulted in substantial economic loses in several countries around the world. Understanding viral pathogenicity mechanisms is of central importance for the establishment of disease control measures. Thus, the present work aimed to identify and characterize miRNAs from OsHV-1 as well as to predict their target transcripts in the virus and the host. OsHV-1 genome was used for the in silico discovery of pre-miRNAs. Subsequently, viral and host target transcripts of the OsHV-1 miRNAs were predicted according to the base pairing interaction between mature miRNAs and mRNA 3' untranslated regions (UTRs). Six unique pre-miRNAs were found in different regions of the viral genome, ranging in length from 85 to 172 nucleotides. A complex network of self-regulation of viral gene expression mediated by the miRNAs was identified. These sequences also seem to have a broad ability to regulate the expression of host immune-related genes, especially those associated with pathogen recognition. Our results suggest that OsHV-1 encodes miRNAs with important functions in the infection process, inducing self-regulation of viral transcripts, as well as affecting the regulation of Pacific oyster transcripts related to immunity. Understanding the molecular basis of host-pathogen interactions can help mitigate the recurrent events of oyster mass mortalities by OsHV-1 observed worldwide.

The Complicated Evolutionary Diversification of the Mpeg-1/Perforin-2 Family in Cnidarians

The invertebrate innate immune system is surprisingly complex, yet our knowledge is limited to a few select model systems. One understudied group is the phylum Cnidaria (corals, sea anemones, etc.). Cnidarians are the sister group to Bilateria and by studying their innate immunity repertoire, a better understanding of the ancestral state can be gained. Corals in particular have evolved a highly diverse innate immune system that can uncover evolutionarily basal functions of conserved genes and proteins. One rudimentary function of the innate immune system is defense against harmful bacteria using pore forming proteins. Macrophage expressed gene 1/Perforin-2 protein (Mpeg-1/P2) is a particularly important pore forming molecule as demonstrated by previous studies in humans and mice, and limited studies in non-bilaterians. However, in cnidarians, little is known about Mpeg-1/P2. In this perspective article, we will summarize the current state of knowledge of Mpeg-1/P2 in invertebrates, analyze identified Mpeg-1/P2 homologs in cnidarians, and demonstrate the evolutionary diversity of this gene family using phylogenetic analysis. We will also show that Mpeg-1 is upregulated in one species of stony coral in response to lipopolysaccharides and downregulated in another species of stony coral in response to white band disease. This data presents evidence that Mpeg-1/P2 is conserved in cnidarians and we hypothesize that it plays an important role in cnidarian innate immunity. We propose that future research focus on the function of Mpeg-1/P2 family in cnidarians to identify its primary role in innate immunity and beyond.

The involvement of a regucalcin in suppressing hemocyte apoptosis in Pacific oyster Crassostrea gigas

Regucalcin (RGN), also known as senescence marker protein-30 (SMP30), plays a vital role in the regulation of Ca²⁺ homeostasis. In the present study, a regucalcin (designated as CgRGN) was identified from Pacific oyster Crassostrea gigas. The complete cDNA sequence of CgRGN was of 1059 bp, containing an open reading frame of 933 bp which encoded a protein of 310 amino acids. The deduced amino acid sequence of CgRGN shared similarity with other RGNs from the genome of C. gigas as well as other species. The mRNA transcripts of CgRGN were universally detected in all tested tissues, with higher level in hepatopancreas, labial palp, and gills. The relative expression level of CgRGN in hemocytes was significantly up-regulated (p < 0.05) at 3, 12, 72, and 96 h after the stimulation of lipopolysaccharide (LPS). After CgRGN expression was interfered by specific CgRGN-dsRNA, the hemocytes apoptosis rate increased dramatically at 12 h post LPS stimulation (1.56 fold, p < 0.01), compared to the control group. The caspase-3 activity in hemocytes and NO concentration in hemolymph increased significantly (p < 0.05) in dsCgRGN injection oysters. These results collectively indicated that CgRGN could suppress LPS-induced apoptosis and be involved in the immune response of oysters.

Developments in marine invertebrate primary culture reveal novel cell morphologies in the model bivalve Crassostrea gigas

Cell culture provides useful model systems used in a wide range of biological applications, but its utility in marine invertebrates is limited due to the lack of immortalised cell lines. Primary cell and tissue cultures are typically used but remain poorly characterised for oysters, which can cause issues with experimental consistency and reproducibility. Improvements to methods of repeatable isolation, culture, and characterisation of oyster cells and tissues are required to help address these issues. In the current study, systematic improvements have been developed to facilitate the culture of primary cells from adult Pacific oyster tissues and identify novel cell morphologies that have not been reported previously. Cultures analysed by light microscopy, qPCR, and live cell imaging demonstrated maintenance of live, metabolically active Pacific oyster cells for several weeks post-explant. Interestingly, whole hearts dissected from adult oysters were found to continue contracting rhythmically up to 8 weeks after being transferred to a tissue culture system. Mantle tissue explants were also actively moving in the culture system. These improvements in primary cell culture of bivalves may be beneficial for research in ecotoxicology, virology, immunology, and genetic resistance to disease.

An activating transcription factor 6 beta (ATF6β) regulates apoptosis of hemocyte during immune response in Crassostrea gigas

The homeostasis of immune cells during immune response is vital for hosts to defend against invaders. Activating transcription factor 6 (ATF6) is an important transcription factor in the unfolded protein response (UPR) to maintaining cellular homeostasis. In the present study, one ATF6 homologue was identified from Pacific oyster Crassostrea gigas (designated as CgATF6β). The full length cDNA of CgATF6β was of 2645 bp with a 1596 bp open reading frame (ORF) encoding a polypeptide of 531 amino acids. The deduced amino acid sequence of CgATF6β was predicted to contain a transmembrane region, a conserved basic leucine zipper (bZIP) domain, a site 1 protease cleavage site, a site 2 protease cleavage site, and a Golgi localization signal. CgATF6β mRNA was constitutively expressed in hemocytes, gill, mantle, gonad, hepatopancreas and labial palp, with a slightly higher expression level in muscle (2.45-fold of that in gill, p < 0.05). After oysters were challenged with Vibrio splendidus, the mRNA expression levels of CgATF6β in hemocytes were significantly up-regulated at 3 h (2.68-fold of that in seawater group, p < 0.01) and peaked at 12 h (3.14-fold of that in seawater group, p < 0.01). The endogenic CgATF6β protein was mainly located in the cytoplasm of oyster hemocytes, and it was significantly transported into the nuclei of hemocytes at 1.5 h after the challenge with V. splendidus. After an injection with CgATF6β dsRNA, the mRNA expression of CgATF6β was knocked down to 0.26-fold of that in dsGFP group (p < 0.01). In CgATF6β dsRNA-injected oysters, the mRNA expressions of glucose-regulated protein 78 (GRP78), calnexin (CNX) and anti-apoptotic B-cell lymphoma-2 (Bcl-2) in hemocytes were significantly decreased at 12 h after V. splendidus challenge, which were 0.65-fold (p < 0.01), 0.54-fold (p < 0.01) and 0.17-fold (p < 0.01) of that in dsGFP-injected oysters, while the apoptotic rate of hemocytes was significantly up-regulated (1.97-fold of that in dsGFP group, p < 0.05). Collectively, these results suggested that CgATF6β was involved in apoptosis inhibition of oyster hemocytes upon V. splendidus challenge by regulating the expression of CgGRP78, CgCNX and CgBcl-2.

Assessing the effects of neonicotinoid insecticide on the bivalve mollusc Mytilus galloprovincialis

In the present work, the marine invertebrate Mytilus galloprovincialis was used as model organism to evaluate the toxic effects of the neonicotinoid Calypso 480 SC (CAL) following 20 days of exposure to sub-lethal concentrations of 7.77 mg L^-1 (0.1% 96 h-LC50) and 77.70 mg L^-1 (1% 96 h-LC50), and a recovery period of 10 days in uncontaminated seawater. Results revealed that exposure to both concentrations of CAL increased significantly mortality rate in the cells of haemolymph and digestive gland, while digestive gland cells were no longer able to regulate cell volume. Exposure significantly reduced haemolymph parameters (Cl^-, Na⁺), affected the enzymatic activities of superoxide dismutase of digestive gland and catalase of gill, and caused also histopathological alterations in digestive gland and gills. Main histological damages detected in mussels were lipofuscin accumulation, focal points of necrosis, mucous overproduction and infiltrative inflammations. Interestingly, alterations persisted after the recovery period in CAL-free water, especially for haemocyte parameters (K⁺, Na⁺, Ca²⁺, lactate dehydrogenase, glucose). A slight recovery of histological conditions was detected. These findings suggested that sub-chronic exposure to the neonicotinoid insecticide caused significant alterations in both cell and tissue parameters of M. galloprovincialis. Considering the ecologically and commercially important role of mussels in coastal waters, a potential risk posed by neonicotinoids to this essential aquatic resource can be highlighted.

Host-microbiota interactions shed light on mortality events in the striped venus clam Chamelea gallina

Mass mortalities due to disease outbreaks have recently affected a number of major taxa in marine ecosystems. Climate- and pollution-induced stress may compromise host immune defenses, increasing the risk of opportunistic diseases. Despite growing evidence that mass mortality events affecting marine species worldwide are strongly influenced by the interplay of numerous environmental factors, the reductionist approaches most frequently used to investigate these factors hindered the interpretation of these multifactorial pathologies. In this study, we propose a broader approach based on the combination of RNA-sequencing and 16S microbiota analyses to decipher the factors underlying mass mortality in the striped venus clam, Chamelea gallina, along the Adriatic coast. On one hand, gene expression profiling and functional analyses of microbial communities showed the over-expression of several genes and molecular pathways involved in xenobiotic metabolism, suggesting potential chemical contamination in mortality sites. On the other hand, the down-regulation of several genes involved in immune and stress response, and the over-representation of opportunistic pathogens such as Vibrio and Photobacterium spp. indicates that these microbial species may take advantage of compromised host immune pathways and defense mechanisms that are potentially affected by chemical exposure, resulting in periodic mortality events. We propose the application of our approach to interpret and anticipate the risks inherent in the combined effects of pollutants and microbes on marine animals in today's rapidly changing environment.

Exogenous Ca2+ mitigates the toxic effects of TiO2 nanoparticles on phagocytosis, cell viability, and apoptosis in haemocytes of a marine bivalve mollusk, Tegillarca granosa

Phagocytosis suppression induced by nanoparticles (NPs) exposure is increasingly reported in marine species. However, the mechanisms underlying this impact remain poorly understood. In order to improve our present understanding of the immunotoxicity of NPs, acute (96 h) TiO₂ NP exposure and rescue trials via exogenous supply of Ca²⁺ were performed in the blood clam, Tegillarca granosa. The results show that the phagocytosis rate, cell viability, and intracellular Ca²⁺ concentration of haemocytes were significantly suppressed, whereas the intracellular ROS concentration of haemocytes significantly increased upon nTiO₂ exposure. Exposure to nTiO₂ also led to the significant downregulation of Caspase-3, Caspase-6, apoptosis regulator Bcl-2, Bcl-2-associated X, calmodulin kinase II, and calmodulin kinase kinase II. Furthermore, the toxic impacts of nTiO₂ were partially mitigated by the addition of exogenous Ca²⁺, as indicated by the recovery tendency in almost all the measured parameters. The present study indicates that Ca²⁺ signaling could be one of the key pathways through which nTiO₂ attacks phagocytosis.

Immune Tolerance in Mytilus galloprovincialis Hemocytes After Repeated Contact With Vibrio splendidus

Mediterranean mussels (Mytilus galloprovincialis) are sessile filter feeders that live in close contact with numerous marine microorganisms. As is the case in all invertebrates, mussels lack an adaptive immune system, but they respond to pathogens, injuries or environmental stress in a very efficient manner. However, it is not known if they are able to modify their immune response when they reencounter the same pathogen. In this work, we studied the transcriptomic response of mussel hemocytes before and after two consecutive sublethal challenges with Vibrio splendidus. The first exposure significantly regulated genes related to inflammation, migration and response to bacteria. However, after the second exposure, the differentially expressed genes were related to the control and inhibition of ROS production and the resolution of the inflammatory response. Our results also show that the second injection with V. splendidus led to changes at the transcriptional (control of the expression of pro-inflammatory transcripts), cellular (shift in the hemocyte population distribution), and functional levels (inhibition of ROS production). These results suggest that a modified immune response after the second challenge allowed the mussels to tolerate rather than fight the infection, which minimized tissue damage.

RNA Interference by Ingested dsRNA-Expressing Bacteria to Study Shell Biosynthesis and Pigmentation in Crassostrea gigas

RNA interference (RNAi) is an important molecular tool for analysis of gene function in vivo. Although the Pacific oyster Crassostrea gigas is an economically important species with fully sequenced genome, very few mechanistic studies have been carried out due to the lack of molecular techniques to alter gene expression without inducing stress. In this present study, we used unicellular alga Platymonas subcordiformis and Nitzschia closterium f. minutissima as a vector to feed oysters with Escherichia coli strain HT115 engineered to express double-stranded RNAs (dsRNAs) targeting specific genes involved in shell pigmentation. A C. gigas strain with black shell was used to target tyrosinase or peroxidase gene expression by RNAi using the above-mentioned approach. The results showed that feeding oyster with dsRNA of tyrosinase could knock down the expression of corresponding tyrosinase and hinder the developed shell growth. Feeding oyster with dsRNA of peroxidase could knock down the expression of the corresponding peroxidase and result in reduced black pigmentation in the newly developed shell. This non-invasive RNAi study demonstrated that tyrosinase played a vital role in the assembly and maturation of shell matrices and peroxidase was essential for black pigmentation in the shell. Moreover, the RNA interference by ingested dsRNA-expressing bacteria is a relatively simple and effective method for knockdown of a gene expression in adult oysters, thus further advances the use of C. gigas as model organism in functional genomic studies.

Comparative Transcriptome Analysis Reveals Molecular Basis Underlying Fast Growth of the Selectively Bred Pacific Oyster, Crassostrea gigas

Fast growth is one of the most desired traits for all food animals, which affects the profitability of animal production. The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with the largest annual production. Growth of the Pacific oyster has been greatly improved by artificial selection breeding, but molecular mechanisms underlying growth remains poorly understood, which limited the molecular integrative breeding of fast growth with other superior traits. In this study, comparative transcriptome analyses between the fast-growing selectively bred Pacific oyster and unselected wild Pacific oysters were conducted by RNA-Seq. A total of 1,303 protein-coding genes differentially expressed between fast-growing oysters and wild controls were identified, of which 888 genes were expressed at higher levels in the fast-growing oysters. Functional analysis of the differentially expressed genes (DEGs) indicated that genes involved in microtubule motor activity and biosynthesis of nucleotides and proteins are potentially important for growth in the oyster. Positive selection analysis of genes at the transcriptome level showed that a significant number of ribosomal protein genes had undergone positive selection during the artificial selection breeding process. These results also indicated the importance of protein biosynthesis and metabolism for the growth of oysters. The alternative splicing (AS) of genes was also compared between the two groups of oysters. A total of 3,230 differential alternative splicing events (DAS) were identified, involved in 1,818 genes. These DAS genes were associated with specific functional pathways related to growth, such as “long-term potentiation,” “salivary secretion,” and “phosphatidylinositol signaling system.” The findings of this study will be valuable resources for future investigation to unravel molecular mechanisms underlying growth regulation in the oyster and other marine invertebrates and to provide solid support for breeding application to integrate fast growth with other superior traits in the Pacific oyster.

Exploring First Interactions Between Ostreid Herpesvirus 1 (OsHV-1) and Its Host, Crassostrea gigas: Effects of Specific Antiviral Antibodies and Dextran Sulfate

Viral entry mechanisms of herpesviruses constitute a highly complex process which implicates several viral glycoproteins and different receptors on the host cell surfaces. This initial infection stage was currently undescribed for Ostreid herpes virus 1 (OsHV-1), a herpesvirus infecting bivalves including the Pacific oyster, Crassostrea gigas. To identify OsHV-1 glyproteins implicated in the attachment of the virus to oyster cells, three viral putative membrane proteins, encoded by ORF 25, 41, and 72, were selected and polyclonal antibodies against these targets were used to explore first interactions between the virus and host cells. In addition, effects of dextran sulfate, a negative charged sulfated polysaccharide, were investigated on OsHV-1 infection. Effects of antiviral antibodies and dextran sulfate were evaluated by combining viral DNA and RNA detection in spat (in vivo trials) and in oyster hemolymph (in vitro trials). Results showed that viral protein encoded by ORF 25 appeared to be involved in interaction between OsHV-1 and host cells even if other proteins are likely implicated, such as proteins encoded by ORF 72 and ORF 41. Dextran sulfate at 30 μg/mL significantly reduced the spat mortality rate in the experimental conditions. Taken together, these results contribute to better understanding the pathogenesis of the viral infection, especially during the first stage of OsHV-1 infection, and open the way toward new approaches to control OsHV-1 infection in confined facilities.

Dual Transcriptomic Analysis Reveals a Delayed Antiviral Response of Haliotis diversicolor supertexta against Haliotid Herpesvirus-1

Haliotid herpesvirus-1 (HaHV-1) is the first identified gastropod herpesvirus, causing a highly lethal neurologic disease of abalone species. The genome of HaHV-1 has been sequenced, but the functions of the putative genes and their roles during infection are still poorly understood. In the present study, transcriptomic profiles of Haliotis diversicolor supertexta at 0, 24 and 60 h post injection (hpi) with HaHV-1 were characterized through high-throughput RNA sequencing. A total of 448 M raw reads were obtained and assembled into 2.08 × 10⁵ unigenes with a mean length of 1486 bp and an N50 of 2455 bp. Although we detected increased HaHV-1 DNA loads and active viral expression at 24 hpi, this evidence was not linked to significant changes of host transcriptomic profiles between 0 and 24 hpi, whereas a rich immune-related gene set was over-expressed at 60 hpi. These results indicate that, at least at the beginning of HaHV-1 infection, the virus can replicate with no activation of the host immune response. We propose that HaHV-1 may evolve more effective strategies to modulate the host immune response and hide during replication, so that it could evade the immune surveillance at the early stage of infection.

The immune system of the freshwater zebra mussel, Dreissena polymorpha, decrypted by proteogenomics of hemocytes and plasma compartments

The immune system of bivalves is of great interest since it reflects the health status of these organisms during stressful conditions. While immune molecular responses are well documented for marine bivalves, few information is available for continental bivalves such as the zebra mussel, Dreissena polymorpha. A proteogenomic approach was conducted on both hemocytes and plasma to identified immune proteins of this non-model species. Combining transcriptomic sequences with mass spectrometry data acquired on proteins is a relevant strategy since 3020 proteins were identified, representing the largest protein inventory for this sentinel organism. Functional annotation and gene ontology (GO) analysis performed on the identified proteins described the main molecular players of hemocytes and plasma in immunity. GO analysis highlights the complementary immune functions of these two compartments in the management of micro-organisms. Functional annotation revealed new mechanisms in the immune defence of the zebra mussel. Proteins rarely observed in the hemolymph of bivalves were pinpointed such as natterin-like and thaumatin-like proteins. Furthermore, the high abundance of complement-related proteins observed in plasma suggested a strong implication of the complement system in the immune defence of D. polymorpha. This work brings a better understanding of the molecular mechanisms involved in zebra mussel immunity. SIGNIFICANCE: Although the molecular mechanisms of marine bivalves are widely investigated, little information is known for continental bivalves. Moreover, few proteomic studies described the complementarity of both hemolymphatic compartments (cellular and plasmatic) in the immune defence of invertebrates. The recent proteogenomics concept made it possible to discover proteins in non-model organisms. Here, we propose a proteogenomic strategy with the zebra mussel, a key sentinel species for biomonitoring of freshwater, to identify and describe the molecular actors involved in the immune system in both hemocytes and plasma compartments. More widely, this study provided new insight into bivalve immunity.

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.

RNA-seq of HaHV-1-infected abalones reveals a common transcriptional signature of Malacoherpesviruses

Haliotid herpesvirus-1 (HaHV-1) is the viral agent causative of abalone viral ganglioneuritis, a disease that has severely affected gastropod aquaculture. Although limited, the sequence similarity between HaHV-1 and Ostreid herpesvirus-1 supported the assignment of both viruses to Malacoherpesviridae, a Herpesvirales family distantly related with other viruses. In this study, we reported the first transcriptional data of HaHV-1, obtained from an experimental infection of Haliotis diversicolor supertexta. We also sequenced the genome draft of the Chinese HaHV-1 variant isolated in 2003 (HaHV-1-CN2003) by PacBio technology. Analysis of 13 million reads obtained from 3 RNA samples at 60 hours post injection (hpi) allowed the prediction of 51 new ORFs for a total of 117 viral genes and the identification of 207 variations from the reference genome, consisting in 135 Single Nucleotide Polymorphisms (SNPs) and 72 Insertions or Deletions (InDels). The pairing of genomic and transcriptomic data supported the identification of 60 additional SNPs, representing viral transcriptional variability and preferentially grouped in hotspots. The expression analysis of HaHV-1 ORFs revealed one putative secreted protein, two putative capsid proteins and a possible viral capsid protease as the most expressed genes and demonstrated highly synchronized viral expression patterns of the 3 infected animals at 60 hpi. Quantitative reverse transcription data of 37 viral genes supported the burst of viral transcription at 30 and 60 hpi during the 72 hours of the infection experiment, and allowed the distinction between early and late viral genes.

Dual transcriptomic analysis of Ostreid herpesvirus 1 infected Scapharca broughtonii with an emphasis on viral anti-apoptosis activities and host oxidative bursts

The ark shell, Scapharca (Anadara) broughtonii, is an economically important marine shellfish species in Northwestern Pacific. Mass mortalities of ark shell adults related to Ostreid herpesvirus-1 (OsHV-1) infection have occurred frequently since 2012. However, due to the lack of transcriptomic resource of ark shells, the molecular mechanisms underpinning the virus-host interaction remains largely undetermined. In the present study, we resolved the dual transcriptome changes of OsHV-1 infected ark shell with Illumina sequencing. A total of 44 M sequence reads were generated, of which 67,119 reads were mapped to the OsHV-1 genome. De novo assembly of host reads resulted in 276,997 unigenes. 74,529 (26.90%), 47,653 (17.20%) and 19, 611 (7.07%) unigenes were annotated into GO, KOG and KEGG database, respectively. According to RSEM expression values, we identified 2998 differentially expressed genes (DEGs) between control and challenged groups, which included 2065 up-regulated unigenes and 933 down-regulated unigenes. Further analysis of functional pathways indicated that OsHV-1 could inhibit host cell apoptosis mainly by the up-regulation of inhibitor of apoptosis protein (IAP), and thus facilitating its successful replication. While host hemoglobins could induce oxidative burst by suppressing its peroxidase activity, and thus defense against OsHV-1 infection. Although we reported a narrow expression of the OsHV-1 genome compared to Crassostrea gigas infection, we highlighted several common viral genes highly expressed in the two hosts, suggesting an important functional role. This study offers insights into the pathogenesis mechanisms of OsHV-1 infection in bivalve mollusks of the Arcidae family.

Waterborne Cd2+ weakens the immune responses of blood clam through impacting Ca2+ signaling and Ca2+ related apoptosis pathways

Exposure to heavy metals such as Cadmium (Cd) may exert detrimental impacts on the immune responses of marine bivalve species. However, the immunotoxicity of Cd on blood clams remains unknown to date. Furthermore, though Cd²⁺ is known to compete with calcium (Ca²⁺) ions for their binding sites in cells and inhibit Ca²⁺ influx, whether Cd²⁺ weakens the immune responses of marine bivalves through inducing intracellular Ca²⁺ disorders still remains unclear. Therefore, the immunotoxicity of Cd²⁺ at different waterborne Ca²⁺ concentrations on blood clam, Tegillarca granosa, were investigated in the present study. Results obtained demonstrated that the total number, phagocytic activity, and red granulocytes ratio of the haemocytes were all significantly reduced after 10 days exposure of individuals to 25 μg/L Cd²⁺. However, when the waterborne Ca²⁺ concentrations were elevated by 10% and 20% (approximately 370 and 410 mg/L, respectively), mitigation effects on the immune responses of individuals were detected. In addition, though the expressions of genes from the Ca²⁺ signaling and Ca²⁺-related apoptosis pathways were significantly altered by Cd²⁺ exposure, the expression patterns of these genes were similar to that of the control when the waterborne Ca²⁺ concentrations were elevated, suggesting a relieving effect of waterborne Ca²⁺ on Cd²⁺ induced toxicity to haemocytes. The results obtained in the present study revealed that waterborne Cd²⁺ may hamper the immune responses of T. granosa through influencing Ca²⁺ signaling and Ca²⁺-related apoptosis pathways, which can be partially mitigated by elevating the waterborne Ca²⁺ concentrations.

Selection of reference genes for quantitative real-time PCR in Octopus minor (Cephalopoda: Octopoda) under acute ammonia stress

High concentration of ammonia is a common issue in the aquaculture industry and often causes detrimental effects to aquatic products. Exploring expression regulation of genes involved in acute ammonia stress can help to understand the molecular mechanisms of ammonia response. Quantitative real-time PCR (qRT-PCR) with proper reference genes is an effective way to normalize the expression of target genes. To identify the most suitable reference genes for Octopus minor (Mollusca: Cephalopoda: Octopoda) under acute ammonia stress and the normal culturing, nine candidate genes were selected for the validation: OD, RPS18b, RPL29, RPS5, EF-1α, RPL6, AA4, ACT, TUB. The results showed that the stability of candidate genes varied considerably in gill, digestive gland, brain and hemolymph. Thus, the reference genes were determined separately in different tissues. RPS5 and RPL6 showed relatively high stability in gill, while RPL6, TUB, RPL29 and OD were four suitable reference genes in digestive gland. EF-1α, TUB and RPL6 were the best combination in brain and EF-1α and RPS18b were the most appropriate reference genes in hemolymph.

A conserved interferon regulation factor 1 (IRF-1) from Pacific oyster Crassostrea gigas functioned as an activator of IFN pathway

Interferon regulatory factors (IRFs), a family of transcription factors with a novel helix-turn-helix DNA-binding motif, play important roles in regulating the expression of interferons (IFNs) and IFN-stimulated genes. In the present study, an interferon regulation factor 1 was identified from oyster Crassostrea gigas (designated CgIRF-1), and its immune function was characterized to understand the regulatory mechanism of interferon system against viral infection in invertebrates. The open reading frame (ORF) of CgIRF-1 was 990 bp, encoding a polypeptide of 329 amino acids with a typical IRF domain (also known as DNA-binding domain). The mRNA transcripts of CgIRF-1 were detected in all the tested tissues with the highest expression level in hemocyte. CgIRF-1 protein was distributed in both nucleus and cytoplasm of the oyster hemocyte. The mRNA expression of CgIRF-1 in hemocytes was significantly up-regulated at 48 h after poly (I:C) stimulation (p < 0.05). The recombinant CgIRF-1 (rCgIRF-1) could interact with classically IFN-stimulated response elements (ISRE) in vitro. The relative luciferase activity of interferon-like protein promotor reporter gene (pGL-CgIFNLP promotor) was significantly (p < 0.05) enhanced in HEK293T cell after transfection of CgIRF-1. These results indicated that CgIRF-1 could bind ISRE and regulate the expression of CgIFNLP as a transcriptional regulatory factor, and participated in the antiviral immune response of oysters.

Antiviral Defense and Innate Immune Memory in the Oyster

The Pacific oyster, Crassostrea gigas, is becoming a valuable model for investigating antiviral defense in the Lophotrochozoa superphylum. In the past five years, improvements to laboratory-based experimental infection protocols using Ostreid herpesvirus I (OsHV-1) from naturally infected C. gigas combined with next-generation sequencing techniques has revealed that oysters have a complex antiviral response involving the activation of all major innate immune pathways. Experimental evidence indicates C. gigas utilizes an interferon-like response to limit OsHV-1 replication and spread. Oysters injected with a viral mimic (polyI:C) develop resistance to OsHV-1. Improved survival following polyI:C injection was found later in life (within-generational immune priming) and in the next generation (multi-generational immune priming). These studies indicate that the oyster's antiviral defense system exhibits a form of innate immune-memory. An important priority is to identify the molecular mechanisms responsible for this phenomenon. This knowledge will motivate the development of practical and cost-effective treatments for improving oyster health in aquaculture.

The oyster immunity

Oysters, the common name for a number of different bivalve molluscs, are the worldwide aquaculture species and also play vital roles in the function of ecosystem. As invertebrate, oysters have evolved an integrated, highly complex innate immune system to recognize and eliminate various invaders via an array of orchestrated immune reactions, such as immune recognition, signal transduction, synthesis of antimicrobial peptides, as well as encapsulation and phagocytosis of the circulating haemocytes. The hematopoietic tissue, hematopoiesis, and the circulating haemocytes have been preliminary characterized, and the detailed annotation of the Pacific oyster Crassostrea gigas genome has revealed massive expansion and functional divergence of innate immune genes in this animal. Moreover, immune priming and maternal immune transfer are reported in oysters, suggesting the adaptability of invertebrate immunity. Apoptosis and autophagy are proved to be important immune mechanisms in oysters. This review will summarize the research progresses of immune system and the immunomodulation mechanisms of the primitive catecholaminergic, cholinergic, neuropeptides, GABAergic and nitric oxidase system, which possibly make oysters ideal model for studying the origin and evolution of immune system and the neuroendocrine-immune regulatory network in lower invertebrates.

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.

Transcriptional profiling of long non-coding RNAs in mantle of Crassostrea gigas and their association with shell pigmentation

Long non-coding RNAs (lncRNAs) play crucial roles in diverse biological processes and have drawn extensive attention in the past few years. However, lncRNAs remain poorly understood about expression and roles in Crassostrea gigas, a potential model organism for marine molluscan studies. Here, we systematically identified lncRNAs in the mantles of C. gigas from four full-sib families characterized by white, black, golden, and partially pigmented shell. Using poly(A)-independent and strand-specific RNA-seq, a total of 441,205,852 clean reads and 12,243 lncRNA transcripts were obtained. LncRNA transcripts were relatively short with few exons and low levels of expression in comparison to protein coding mRNA transcripts. A total of 427 lncRNAs and 349 mRNAs were identified to differentially express among six pairwise groups, mainly involving in biomineralization and pigmentation through functional enrichment. Furthermore, a total of 6 mRNAs and their cis-acting lncRNAs were predicted to involve in synthesis of melanin, carotenoid, tetrapyrrole, or ommochrome. Of them, chorion peroxidase and its cis-acting lincRNA TCONS_00951105 are implicated in playing an essential role in the melanin synthetic pathway. Our studies provided the first systematic characterization of lncRNAs catalog expressed in oyster mantle, which may facilitate understanding the molecular regulation of shell colour diversity and provide new insights into future selective breeding of C. gigas for aquaculture.

Effects of phenanthrene on early development of the Pacific oyster Crassostrea gigas (Thunberg, 1789)

Phenanthnere (PHE) is a polycyclic aromatic hydrocarbon continuously discarded in the marine environment and bioavailable to many aquatic species. Although studies about PHE toxicity have been documented for adult oysters, the effects on early developmental stages are poorly characterized in bivalves. In this study, the effects of PHE (0.02 and 2.0μg.L^-1) were evaluated on the embryogenesis and larval development of Crassostrea gigas. Toxicity bioassays, growth and deformities assessment, analysis of shell calcium abundance and transcript levels of genes related to xenobiotic biotransformation (CYP2AU2, CYP30C1), immune system (Cg-Tal) and tissue growth and shell formation (Ferritin, Insulin-like, Cg-Try, Calmodulin and Nacrein) were assayed in D-shape larvae after 24h of PHE exposure. At the highest concentration (2.0μg.L^-1), PHE decreased the frequency of normal development (19.7±2.9%) and shell size (53.5±2.8mm). Developmental deformities were mostly related to abnormal mantle and shell formation. Lower calcium levels in oyster shells exposed to PHE 2.0μg.L^-1 were observed, suggesting effects on shell structure. At this same PHE concentration, CYP30C1, Cg-Tal, Cg-Tyr, Calmodulin were upregulated and CYP2AU2, Ferritin, Nacrein, and Insulin-Like were downregulated compared to control larvae. At the lowest PHE concentration (0.02μg.L^-1), it was observed a minor decrease in normal larval development (89,6±6%) and the remaining parameters were not affected. This is the first study to provide evidences that exposure to PHE can affect early oyster development at the molecular and morphological levels, possibly threatening this bivalve species.

Recent findings on phenoloxidases in bivalves

The production of melanin is a complex process involving biochemical cascades, such as the pro-phenoloxidase (proPO) system, and enzymes, such as phenoloxidases (POs). Different studies have shown a strong correlation between the decrease in PO activities and the occurrence of diseases in bivalve invertebrates, leading to mortalities in the host. Results of these studies suggest that POs could play a fundamental role in defense mechanisms in bivalves. This article reviews the fundamental knowledge on the proPO system in bivalves and the methods used to assess PO activities. Finally, this is the first report on the major findings of laboratory and field studies that indicate that a type of PO in bivalves, the laccase enzyme, is inducible and involved in the 1) immune 2) antioxidant and 3) detoxification roles in bivalves, and might be an ecological potential biomarker of environmental stress.

Haemocytes from Crassostrea gigas and OsHV-1: A promising in vitro system to study host/virus interactions

Since 2008, mass mortality outbreaks associated with the detection of particular variants of OsHV-1 have been reported in Crassostrea gigas spat and juveniles in several countries. Recent studies have reported information on viral replication during experimental infection. Viral DNA and RNA were also detected in the haemolymph and haemocytes suggesting that the virus could circulate through the circulatory system. However, it is unknown if the virus is free in the haemolymph, passively associated at the surface of haemocytes, or able to infect and replicate inside these cells inducing (or not) virion production. In the present study, we collected haemocytes from the haemolymphatic sinus of the adductor muscle of healthy C. gigas spat and exposed them in vitro to a viral suspension. Results showed that viral RNAs were detectable one hour after contact and the number of virus transcripts increased over time in association with an increase of viral DNA detection. These results suggested that the virus is able to initiate replication rapidly inside haemocytes maintained in vitro. These in vitro trials were also used to carry out a dual transcriptomic study. We analyzed concomitantly the expression of some host immune genes and 15 viral genes. Results showed an up regulation of oyster genes currently studied during OsHV-1 infection. Additionally, transmission electron microscopy examination was carried out and did not allow the detection of viral particles. Moreover, All the results suggested that the in vitro model using haemocytes can be valuable for providing new perspective on virus-oyster interactions.