Transcriptome analysis reveals strong and complex antiviral response in a mollusc

Mechanisms of HAHV-1 Interaction with Hemocytes in Haliotis diversicolor supertexta: An In Vitro Study

Haliotid herpesvirus 1 (HAHV-1) causes significant damage to the abalone aquaculture industry. Knowledge of HAHV-1 invasion and host defense mechanisms is limited due to the lack of stable molluscan cell lines. The present study established an in vitro infection model of HAHV-1 using the primary suspension cultures of hemocytes from Haliotis diversicolor supertexta and Haliotis discus hannai. The cytopathic effects of HAHV-1 on adherent-cultured hemocytes of both species were also investigated. The HAHV-1 DNA loads were firstly monitored by means of quantitative PCR during the development of viral infection, and subsequently the mechanism of interaction between HAHV-1 and hemocytes was explored by means of a transcriptome analysis. H. diversicolor supertexta hemocytes exhibited a high degree of susceptibility to HAHV-1, with viral loads reaching a peak of 4.0 × 10⁷ copies/ng DNA. In contrast, no significant replication was observed in H. discus hannai hemocytes. Transcriptome analysis revealed that HAHV-1 evades the host immune response in the early stages of infection, and hijacks the host's energy and redox metabolism to promote its replication at the late stages. Consequently, this study provides a valuable reference point for the investigation of virus-host interaction between HAHV-1 and abalone in vitro.

Transcriptomic and functional analysis of the antiviral response of mussels after a poly I:C stimulation

The study of mussels (Mytilus galloprovincialis) has grown in importance in recent years due to their high economic value and resistance to pathogens. Because of the biological characteristics revealed by mussel genome sequencing, this species is a valuable research model. The high genomic variability and diversity, particularly in immune genes, may be responsible for their resistance to pathogens found in seawater and continuously filtered and internalized by them. These facts, combined with the lack of proven mussel susceptibility to viruses in comparison to other bivalves such as oysters, result in a lack of studies on mussel antiviral response. We used RNA-seq to examine the genomic response of mussel hemocytes after they were exposed to poly I:C, simulating immune cell contact with viral dsRNA. Apoptosis and the molecular axis IRFs/STING-IFI44/IRGC1 were identified as the two main pathways in charge of the response but we also found a modulation of lncRNAs. Finally, in order to obtain new information about the response of mussels to putative natural challenges, we used VHSV virus (Viral Hemorrhagic Septicemia Virus) to run some functional analysis and confirm poly I:C's activity as an immunomodulator in a VHSV waterborne stimulation. Both, poly I:C as well as an injury stimulus (filtered sea water injection) accelerated the viral clearance by hemocytes and altered the expression of several immune genes, including IL-17, IRF1 and viperin.

Localization and Tissue Tropism of Ostreid Herpesvirus 1 in Blood Clam Anadara broughtonii

OsHV-1 caused detrimental infections in a variety of bivalve species of major importance to aquaculture worldwide. Since 2012, there has been a notable increase in the frequency of mass mortality events of the blood clam associated with OsHV-1 infection. The pathological characteristics, tissue and cellular tropisms of OsHV-1 in A. broughtonii remain unknown. In this study, we sought to investigate the distribution of OsHV-1 in five different organs (mantle, hepatopancreas, gill, foot, and adductor muscle) of A. broughtonii by quantitative PCR, histopathology and in situ hybridization (ISH), to obtain insight into the progression of the viral infection. Our results indicated a continuous increase in viral loads with the progression of OsHV-1 infection, reaching a peak at 48 h or 72 h post-infection according to different tissues. Tissue damage and necrosis, as well as colocalized OsHV-1 ISH signals, were observed primarily in the connective tissues of various organs and gills. Additionally, minor tissue damage accompanied by relatively weak ISH signals was detected in the foot and adductor muscle, which were filled with muscle tissue. The predominant cell types labeled by ISH signals were infiltrated hemocytes, fibroblastic-like cells, and flat cells in the gill filaments. These results collectively illustrated the progressive alterations in pathological confusion and OsHV-1 distribution in A. broughtonii, which represent most of the possible responses of cells and tissues to the virus.

MicroRNA transcriptome analysis reveals the immune regulatory mechanism of Crassostrea hongkongesis against Vibrio harveyi infection

MicroRNAs (miRNAs) are small non-coding RNA molecules that modulate target-genes expression and play crucial roles in post-transcriptional regulation and immune system regulation. The Hong Kong oyster (Crassostrea hongkongesis), as the main marine aquaculture shellfish in the South China Sea, not only has high economic and ecological value, but also is an ideal model for conducting research on pathogen host interaction. Vibrio harveyi, a Gram negative luminescent marine bacterium, is widely distributed in coastal water environments and can cause large-scale death of C. hongkongesis. However, little in formation is available on the immune regulatory mechanisms of C. hongkongesis infected with V. harveyi. Therefore, we performed microRNA transcriptome analysis for elucidating the immunoregulation mechanism of C. hongkongesis infected with V. harveyi. The results show that a total of 308468208 clean reads and 288371159 clean tags were obtained. 222 differentially expressed miRNAs were identified. A total of 388 target genes that were differentially expressed and negatively correlated with miRNA expression were predicted by 222 DEmiRs. GO enrichment analysis of 388 DETGs showed that they were mainly enriched in the immune-related term of membrane-bounded vesicle, endocytic vesicle lumen, antigen processing and presentation of exogenous peptide antigen via MHC class I, antigen processing and presentation of peptide antigen via MHC class I, and other immune-related term. KEGG enrichment analysis showed that DETGs were mainly enriched in the Complement and coagulation cascades, Herpes simplex virus 1 infection, Bacterial invasion of epithelial cells, Antigen processing and presentation and NOD-like receptor signaling pathway. The 16 key DEmiRs and their target genes form a regulatory network for seven immune-related pathways. These results suggest that V. harveyi infection induces a complex miRNA response with wide-ranging effects on immune gene expression in the C. hongkongesis. This study explored the immune response of C. hongkongesis to V. harveyi infection at the level of miRNAs, which provides new ideas for the healthy culture and selective breeding of C. hongkongesis.

Transcriptomic Response of the Atlantic Surfclam (Spisula solidissima) to Acute Heat Stress

There is clear evidence that the oceans are warming due to anthropogenic climate change, and the northeastern coast of USA contains some of the fastest warming areas. This warming is projected to continue with serious biological and social ramifications for fisheries and aquaculture. One species particularly vulnerable to warming is the Atlantic surfclam (Spisula solidissima). The surfclam is a critically important species, linking marine food webs and supporting a productive, lucrative, and sustainable fishery. The surfclam is also emerging as an attractive candidate for aquaculture diversification, but the warming of shallow coastal farms threatens the expansion of surfclam aquaculture. Little is known about the adaptive potential of surfclams to cope with ocean warming. In this study, the surfclam transcriptome under heat stress was examined. Two groups of surfclams were subjected to heat stress to assess how artificial selection may alter gene expression. One group of clams had been selected for greater heat tolerance (HS) and the other was composed of random control clams (RC). After a 6-h exposure to 16 or 29 °C, gill transcriptome expression profiles of the four temperature/group combinations were determined by RNA sequencing and compared. When surfclams experienced heat stress, they exhibited upregulation of heat shock proteins (HSPs), inhibitors of apoptosis (IAPs), and other stress-response related genes. RC clams differentially expressed 1.7 times more genes than HS clams, yet HS clams had a stronger response of key stress response genes, including HSPs, IAPs, and genes involved with mitigating oxidative stress. The findings imply that the HS clams have a more effective response to heat stress after undergoing the initial selection event due to genetic differences created by the selection, epigenetic memory of the first heat shock, or both. This work provides insights into how surfclams adapt to heat stress and should inform future breeding programs that attempt to breed surfclam for greater heat tolerance, and ultimately bring greater resiliency to shellfish farms.

Transcriptome analysis of digestive diverticula of Hong Kong oyster (Crassostrea hongkongesis) infected with Vibrio harveyi

The Hong Kong oyster (Crassostrea hongkongesis), as the main marine aquaculture shellfish in the South China Sea, not only has high economic and ecological value, but also is an ideal model for conducting research on pathogen host interaction. However, diseases caused by Vibrio pose a serious impediment to the culture of C. hongkongesis. In this study, we performed transcriptome analysis of digestive diverticula of C. hongkongesis infected with V. harveyi. A total of 977, 689, 912 high quality reads and 955, 208, 562 valid reads were obtained. At 12, 24, 48 and 72 h post-infection, 1402, 2168, 2727 and 1398 differentially expressed genes (DEGs) were captured, respectively. GO enrichment analysis showed that DEGs were significantly enriched in cellular processes, catalytic activity, cell part and other terms. KEGG enrichment analysis revealed that these DEGs were mainly closely related to Necroptosis, RIG-I-like receptor signaling pathway, NF-kappa B signaling pathway, Toll-like receptor signaling pathway and other pathways are related. The results of WGCNA analysis indicated that THBS1, CA10, Trpm2, THAP12, PTPRT, HSPA12A, and ADAM10 were the hub genes in the gene co-expression network. This study will provide new ideas at the transcriptome level for the immune regulatory mechanisms and adaptability of the C. hongkongesis to V. infection, as well as for achieving selective breeding for Vibrio resistance in the C. hongkongesis.

The Resistance to Lethal Challenge with Ostreid herpesvirus-1 of Pacific Oysters (Crassostrea gigas) Previously Exposed to This Virus

Pacific oyster (Crassostrea gigas) aquaculture has been economically impacted in many countries by Pacific oyster mortality syndrome (POMS), a disease initiated by Ostreid herpesvirus 1. The objectives of this study were to determine whether naturally exposed, adult C. gigas could act as reservoirs for OsHV-1 and explain the recurrent seasonal outbreaks of POMS and to test whether or not they were resistant to OsHV-1. In a laboratory infection experiment using thermal shock, OsHV-1 replication was not reactivated within the tissues of such oysters and the virus was not transmitted to naïve cohabitating spat. The adult oysters were resistant to intramuscular injection with a lethal dose of OsHV-1 and had 118 times lower risk of mortality than naïve oysters. Considered together with the results of other studies in C. gigas, natural exposure or laboratory exposure to OsHV-1 may result in immunity during subsequent exposure events, either in the natural environment or the laboratory. While adult C. gigas can carry OsHV-1 infection for lengthy periods, reactivation of viral replication leading to mortality and transmission of the virus to naïve oysters may require specific conditions that were not present in the current experiment. Further investigation is required to evaluate the mechanisms responsible for resistance to disease in oysters previously exposed to OsHV-1, whether immunity can be exploited commercially to prevent POMS outbreaks and to determine the source of the virus for recurrent seasonal outbreaks.

cGLRs are a diverse family of pattern recognition receptors in innate immunity

Cyclic GMP-AMP synthase (cGAS) is an enzyme in human cells that controls an immune response to cytosolic DNA. Upon binding DNA, cGAS synthesizes a nucleotide signal 2'3'-cGAMP that activates STING-dependent downstream immunity. Here, we discover that cGAS-like receptors (cGLRs) constitute a major family of pattern recognition receptors in innate immunity. Building on recent analysis in Drosophila, we identify >3,000 cGLRs present in nearly all metazoan phyla. A forward biochemical screening of 150 animal cGLRs reveals a conserved mechanism of signaling including response to dsDNA and dsRNA ligands and synthesis of isomers of the nucleotide signals cGAMP, c-UMP-AMP, and c-di-AMP. Combining structural biology and in vivo analysis in coral and oyster animals, we explain how synthesis of distinct nucleotide signals enables cells to control discrete cGLR-STING signaling pathways. Our results reveal cGLRs as a widespread family of pattern recognition receptors and establish molecular rules that govern nucleotide signaling in animal immunity.

Mussel antiviral transcriptome response and elimination of viral haemorrhagic septicaemia virus (VHSV)

As filter-feeding bivalves, mussels have been traditionally studied as possible vectors of different bacterial or viral pathogens. The absence of a known viral pathogen in these bivalves makes it particularly interesting to study the interaction of the mussel innate immune system with a virus of interest. In the present work, mussels were challenged with viral haemorrhagic septicaemia virus (VHSV), which is a pathogen in several fish species. The viral load was eliminated after 24 h and mussels evidenced antiviral activity towards VHSV, demonstrating that the virus was recognized and eliminated by the immune system of the host and confirming that mussels are not VHSV vectors in the marine environment. The transcriptome activating the antiviral response was studied, revealing the involvement of cytoplasmic viral sensors with the subsequent activation of the JAK-STAT pathway and several downstream antiviral effectors. The inflammatory response was inhibited with the profound downregulation of MyD88, shifting the immune balance towards antiviral functions. High modulation of retrotransposon activity was observed, revealing a mechanism that facilitates the antiviral response and that had not been previously observed in these species. The expression of several inhibitors of apoptosis and apoptosis-promoting genes was modulated, although clear inhibition of apoptosis in bivalves after severe viral infection and subsequent disease was not observed in this study. Finally, the modulated expression of several long noncoding RNAs that were correlated with genes involved in the immune response was detected.

Development and Evaluation of High-Density SNP Arrays for the Eastern Oyster Crassostrea virginica

The eastern oyster Crassostrea virginica is a major aquaculture species for the USA. The sustainable development of eastern oyster aquaculture depends upon the continued improvement of cultured stocks through advanced breeding technologies. The Eastern Oyster Breeding Consortium (EOBC) was formed to advance the genetics and breeding of the eastern oyster. To facilitate efficient genotyping needed for genomic studies and selection, the consortium developed two single-nucleotide polymorphism (SNP) arrays for the eastern oyster: one screening array with 566K SNPs and one breeders' array with 66K SNPs. The 566K screening array was developed based on whole-genome resequencing data from 292 oysters from Atlantic and Gulf of Mexico populations; it contains 566,262 SNPs including 47K from protein-coding genes with a marker conversion rate of 48.34%. The 66K array was developed using best-performing SNPs from the screening array, which contained 65,893 oyster SNPs including 22,984 genic markers with a calling rate of 99.34%, a concordance rate of 99.81%, and a much-improved marker conversion rate of 92.04%. Null alleles attributable to large indels were found in 13.1% of the SNPs, suggesting that copy number variation is pervasive. Both arrays provided easy identification and separation of selected stocks from wild progenitor populations. The arrays contain 31 mitochondrial SNPs that allowed unambiguous identification of Gulf mitochondrial genotypes in some Atlantic populations. The arrays also contain 756 probes from 13 oyster and human pathogens for possible detection. Our results show that marker conversion rate is low in high polymorphism species and that the two-step process of array development can greatly improve array performance. The two arrays will advance genomic research and accelerate genetic improvement of the eastern oyster by delineating genetic architecture of production traits and enabling genomic selection. The arrays also may be used to monitor pedigree and inbreeding, identify selected stocks and their introgression into wild populations, and assess the success of oyster restoration.

Understanding Crassostrea virginica tolerance of Perkinsus marinus through global gene expression analysis

Disease tolerance, a host's ability to limit damage from a given parasite burden, is quantified by the relationship between pathogen load and host survival or reproduction. Dermo disease, caused by the protozoan parasite P. marinus, negatively impacts survival in both wild and cultured eastern oyster (C. virginica) populations. Resistance to P. marinus has been the focus of previous studies, but tolerance also has important consequences for disease management in cultured and wild populations. In this study we measured dermo tolerance and evaluated global expression patterns of two sensitive and two tolerant eastern oyster families experimentally challenged with distinct doses of P. marinus (0, 106, 107, and 108 parasite spores per gram wet weight, n = 3-5 individuals per family per dose). Weighted Gene Correlation Network Analysis (WGCNA) identified several modules correlated with increasing parasite dose/infection intensity, as well as phenotype. Modules positively correlated with dose included transcripts and enriched GO terms related to hemocyte activation and cell cycle activity. Additionally, these modules included G-protein coupled receptor, toll-like receptor, and tumor necrosis factor pathways, which are important for immune effector molecule and apoptosis activation. Increased metabolic activity was also positively correlated with treatment. The module negatively correlated with infection intensity was enriched with GO terms associated with normal cellular activity and growth, indicating a trade-off with increased immune response. The module positively correlated with the tolerant phenotype was enriched for transcripts associated with "programmed cell death" and contained a large number of tripartite motif-containing proteins. Differential expression analysis was also performed on the 108 dosed group using the most sensitive family as the comparison reference. Results were consistent with the network analysis, but signals for "programmed cell death" and serine protease inhibitors were stronger in one tolerant family than the other, suggesting that there are multiple avenues for disease tolerance. These results provide new insight for defining dermo response traits and have important implications for applying selective breeding for disease management.

Inactivated ostreid herpesvirus-1 induces an innate immune response in the Pacific oyster, Crassostrea gigas, hemocytes

Infectious diseases are a major constraint to the expansion of shellfish production worldwide. Pacific oyster mortality syndrome (POMS), a polymicrobial disease triggered by the Ostreid herpesvirus-1 (OsHV-1), has devastated the global Pacific oyster (Crassostrea gigas) aquaculture industry. Recent ground-breaking research revealed that C. gigas possess an immune memory, capable of adaption, which improves the immune response upon a second exposure to a pathogen. This paradigm shift opens the door for developing 'vaccines' to improve shellfish survival during disease outbreaks. In the present study, we developed an in-vitro assay using hemocytes - the main effectors of the C. gigas immune system - collected from juvenile oysters susceptible to OsHV-1. The potency of multiple antigen preparations (e.g., chemically and physically inactivated OsHV-1, viral DNA, and protein extracts) to stimulate an immune response in hemocytes was evaluated using flow cytometry and droplet digital PCR to measure immune-related subcellular functions and gene expression, respectively. The immune response to the different antigens was benchmarked against that of hemocytes treated with Poly (I:C). We identified 10 antigen preparations capable of inducing immune stimulation in hemocytes (ROS production and positively expressed immune- related genes) after 1 h of exposure, without causing cytotoxicity. These findings are significant, as they evidence the potential for priming the innate immunity of oysters using viral antigens, which may enable cost-effective therapeutic treatment to mitigate OsHV-1/POMS. Further testing of these antigen preparations using an in-vivo infection model is essential to validate promising candidate pseudo-vaccines.

RNA-Seq comparative study reveals molecular effectors linked to the resistance of Pinna nobilis to Haplosporidium pinnae parasite

With the intensification of maritime traffic, recently emerged infectious diseases have become major drivers in the decline and extinction of species. Since 2016, mass mortality events have decimated the endemic Mediterranean Sea bivalve Pinna nobilis, affecting ca. 100% of individuals. These events have largely been driven by Haplosporidium pinnae's infection, an invasive species which was likely introduced by shipping. While monitoring wild populations of P. nobilis, we observed individuals that survived such a mass mortality event during the summer of 2018 (France). We considered these individuals resistant, as they did not show any symptoms of the disease, while the rest of the population in the area was devastated. Furthermore, the parasite was not detected when we conducted a PCR amplification of a species-specific fragment of the small subunit ribosomal DNA. In parallel, the transcriptomic analysis showed evidence of some parasite RNA indicating that the resistant individuals had been exposed to the parasite without proliferating. To understand the underlying mechanisms of resistance in these individuals, we compared their gene expression with that of susceptible individuals. We performed de novo transcriptome assembly and annotated the expressed genes. A comparison of the transcriptomes in resistant and susceptible individuals highlighted a gene expression signature of the resistant phenotype. We found significant differential expressions of genes involved in immunity and cell architecture. This data provides the first insights into how individuals escape the pathogenicity associated with infection.

Ecological significance of G protein-coupled receptors in the Pacific oyster (Crassostrea gigas): Pervasive gene duplication and distinct transcriptional response to marine environmental stresses

Marine ecosystems with ocean warming and industry pollution threaten the survival and adaptation of organisms. G protein-coupled receptors (GPCRs) play critical roles in various physiological and toxicological processes in vertebrates and invertebrates. The Pacific oyster (Crassostrea gigas) was widely used to study the adaptation of marine molluscs to coastal environments. In this work, we identified a total of 586 GPCRs in C. gigas genome. The C. gigas GPCRs were divided into five classes (including class A, B, C, E and F) with different degrees of expansion. Meta-analysis of multiple RNA-seq datasets revealed that transcriptional expression patterns of GPCRs in C. gigas were distinct in response to high temperature, salinity, air exposure, heavy metal, ostreid herpes virus 1 (OsHV-1) and Vibrio challenge. This work for the first time characterized the GPCR gene family and provided insights into the potential roles of GPCRs in adaptation of marine molluscs to stressful coastal environment.

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

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

The first identified invertebrate LGP2-like homolog gene in the Pacific oyster Crassostrea gigas

The LGP2 (Laboratory of Genetics and Physiology 2) protein is a member of the retinoic acid-inducible gene I (RIG-I)-like receptor (RLRs) family, which is a class of antiviral pattern recognition receptors located in the cytoplasm. However, few studies have investigated the function of LGP2 in invertebrates. In this study, the complete coding sequence of the LGP2 gene of the Pacific oyster, Crassostrea gigas, was obtained and named CgLGP2-like. Sequence analysis revealed that CgLGP2-like encodes 803 amino acids, and the encoded protein contains a DEXDc, HELICc, and C-terminal regulatory domains. Multiple sequence alignment demonstrated that the sequences of these key protein functional domains were relatively conserved. Phylogenetic analysis revealed that CgLGP2-like was a new member of the animal LGP2 family. Quantitative real-time PCR results showed that CgLGP2-like mRNA was expressed in all tested oyster tissues, with the highest expression observed in the labial palpus and digestive glands. CgLGP2-like expression in gill tissues was significantly induced after the poly(I:C) challenge. Furthermore, multiple IRF and NF-κB binding sites were identified in the CgLGP2-like promoter region, which may be one of the reasons why CgLGP2-like responds to poly(I:C) stimulation. Finally, the results of dual-luciferase reporter gene assays revealed that overexpression of CgLGP2-like may have a regulatory effect on the human IFN, AP-1, and oyster CgIL-17 genes in HEK293T cells. Overall, our results preliminarily elucidate the immune functions of invertebrate LGP2 protein and provide valuable information for the development of comparative immunology.

Gut Microbiota and Transcriptomics Reveal the Effect of Human Norovirus Bioaccumulation on Oysters (Crassostrea gigas)

Human norovirus (HuNoV) is a major foodborne pathogen that causes acute viral gastroenteritis, and oysters are one of the main carriers of HuNoV transmission. While progress has been made toward understanding the pattern of oyster-bioaccumulated HuNoV, the response of oysters to HuNoV bioaccumulation, including changes in gene expression and gut microbiota, is unclear. In this study, histo-blood group antigen (HBGA)-like molecule expression and gene regulation features and the HuNoV-microbiome interactions of oysters during HuNoV bioaccumulation were characterized. With the prolongation of bioaccumulation time, the HuNoV content and expression of type A HBGA-like molecules in oysters increased and stabilized. HuNoV also altered the expression of immunity- and glycosphingolipid biosynthesis-related genes. Prolonged bioaccumulation of HuNoV can reduce the abundance and change the composition of the oyster gut microbiota. In particular, with the extension of bioaccumulation time, the abundance of Blautia, Agathobacter, Faecalibacterium, Terrisporobacter, Bifidobacterium, Lactobacillus, and Ruminococcus decreased, while the abundance of Vibrio and Alphaproteobacteria increased. This study provides potential candidates for identifying functional genes involved in the bioaccumulation of HuNoV in oysters. More importantly, it provides the first description of the changes in gut microbiota during HuNoV bioaccumulation in oysters. IMPORTANCE The role of the oyster gut microbiota in HuNoV bioaccumulation is poorly understood. This study revealed, for the first time, the changes in gut microbiota and gene expression of oysters with HuNoV bioaccumulation. This study enriches the understanding of the impact of HuNoV bioaccumulation on oysters and provides a new direction for the study of the molecular mechanism of HuNoV bioaccumulation in oysters.

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.

Effects of environmental concentrations of the fragrance amyl salicylate on the mediterranean mussel Mytilus galloprovincialis

Amyl salicylate (AS) is a fragrance massively used as a personal care product and following the discharged in wastewaters may end up in the aquatic environment representing a potential threat for the ecosystem and living organisms. AS was recently detected in water of the Venice Lagoon, a vulnerable area continuously subjected to the income of anthropogenic chemicals. The lagoon is a relevant area for mollusc farming, including the Mediterranean mussels (Mytilus galloprovincialis) having an important economic and ecological role. Despite high levels of AS occurred in water of the Lagoon of Venice, no studies investigated the possible consequences of AS exposures on species inhabiting this ecosystem to date. For the first time, we applied a multidisciplinary approach to investigate the potential effects of the fragrance AS on Mediterranean mussels. To reach such a goal, bioaccumulation, cellular, biochemical, and molecular analyses (RNA-seq and microbiota characterization) were measured in mussels treated for 7 and 14 days with different AS Venice lagoon environmental levels (0.1 and 0.5 μg L^-1). Despite chemical investigations suggested low AS bioaccumulation capability, cellular and molecular analyses highlighted the disruption of several key cellular processes after the prolonged exposures to the high AS concentration. Among them, potential immunotoxicity and changes in transcriptional regulation of pathways involved in energy metabolism, stress response, apoptosis and cell death regulations have been observed. Conversely, exposure to the low AS concentration demonstrated weak transcriptional changes and transient increased representation of opportunistic pathogens, as Arcobacter genus and Vibrio aestuarianus. Summarizing, this study provides the first overview on the effects of AS on one of the most widely farmed mollusk species.

Integrated application of transcriptomics and metabolomics provides insights into condition index difference mechanisms in the Pacific oyster (Crassostrea gigas)

The condition index (CI) is an economically important tool for assessing the quality of oysters, such as the Pacific oyster Crassostrea gigas. However, little is known about the mechanisms that underlie differences in CI between different C. gigas populations. In this study, we integrated transcriptomic and metabolomic profiling to investigate the mechanisms that underlie the differences between high- and low-CI groups in one- and two-year-old populations of C. gigas. The results indicate that differences in CI were associated with the regulation of growth-related genes, the FoxO signaling pathway, and the complex regulation of carbohydrate, lipid, amino acid, and energy metabolism. Moreover, the mechanisms underlying these differences differed between the populations. This study is the first to elucidate the molecular and chemical mechanisms associated with CI, and the results will be helpful for breeding higher quality oysters.

Synergistic Interaction of Low Salinity Stress With Vibrio Infection Causes Mass Mortalities in the Oyster by Inducing Host Microflora Imbalance and Immune Dysregulation

A sudden drop in salinity following extreme precipitation events usually causes mass mortality of oysters exposed to pathogens in ocean environment. While how low salinity stress interacts with pathogens to cause mass mortality remains obscure. In this study, we performed an experiment by low salinity stress and pathogen infection with Vibrio alginolyticus to investigate their synergistic effect on the mortality of the Pacific oyster toward understanding of the interaction among environment, host, and pathogen. We showed that low salinity stress did not significantly affect proliferation and virulence of V. alginolyticus, but significantly altered microbial composition and immune response of infected oysters. Microbial community profiling by 16S rRNA amplicon sequencing revealed disrupted homeostasis of digestive bacterial microbiota with the abundance of several pathogenic bacteria being increased, which may affect the pathogenesis in infected oysters. Transcriptome profiling of infected oysters revealed that a large number of genes associated with apoptosis and inflammation were significantly upregulated under low salinity, suggesting that low salinity stress may have triggered immune dysregulation in infected oysters. Our results suggest that host-pathogen interactions are strongly affected by low salinity stress, which is of great significance for assessing future environmental risk of pathogenic diseases, decoding the interaction among environment, host genetics and commensal microbes, and disease surveillance 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.

The expanded inhibitor of apoptosis gene family in oysters possesses novel domain architectures and may play diverse roles in apoptosis following immune challenge

Background

Apoptosis plays important roles in a variety of functions, including immunity and response to environmental stress. The Inhibitor of Apoptosis (IAP) gene family of apoptosis regulators is expanded in molluscs, including eastern, Crassostrea virginica, and Pacific, Crassostrea gigas, oysters. The functional importance of IAP expansion in apoptosis and immunity in oysters remains unknown.

Results

Phylogenetic analysis of IAP genes in 10 molluscs identified lineage specific gene expansion in bivalve species. Greater IAP gene family expansion was observed in C. virginica than C. gigas (69 vs. 40), resulting mainly from tandem duplications. Functional domain analysis of oyster IAP proteins revealed 3 novel Baculoviral IAP Repeat (BIR) domain types and 14 domain architecture types across gene clusters, 4 of which are not present in model organisms. Phylogenetic analysis of bivalve IAPs suggests a complex history of domain loss and gain. Most IAP genes in oysters (76% of C. virginica and 82% of C. gigas), representing all domain architecture types, were expressed in response to immune challenge (Ostreid Herpesvirus OsHV-1, bacterial probionts Phaeobacter inhibens and Bacillus pumilus, several Vibrio spp., pathogenic Aliiroseovarius crassostreae, and protozoan parasite Perkinsus marinus). Patterns of IAP and apoptosis-related differential gene expression differed between the two oyster species, where C. virginica, in general, differentially expressed a unique set of IAP genes in each challenge, while C. gigas differentially expressed an overlapping set of IAP genes across challenges. Apoptosis gene expression patterns clustered mainly by resistance/susceptibility of the oyster host to immune challenge. Weighted Gene Correlation Network Analysis (WGCNA) revealed unique combinations of transcripts for 1 to 12 IAP domain architecture types, including novel types, were significantly co-expressed in response to immune challenge with transcripts in apoptosis-related pathways.

Conclusions

Unprecedented diversity characterized by novel BIR domains and protein domain architectures was observed in oyster IAPs. Complex patterns of gene expression of novel and conserved IAPs in response to a variety of ecologically-relevant immune challenges, combined with evidence of direct co-expression of IAP genes with apoptosis-related transcripts, suggests IAP expansion facilitates complex and nuanced regulation of apoptosis and other immune responses in oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08233-6.

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.

Withering syndrome induced gene expression changes and a de-novo transcriptome for the Pinto abalone, Haliotis kamtschatkana

In the abalone and Candidatus Xenohaliotis californiensis (Ca. Xc) system, the Ca. Xc bacterium infects abalone digestive tissues and leads to extreme starvation and a characteristic "withering" of the gastropod foot. First identified in black abalone in California after an El Niño event, withering syndrome (WS) has caused large declines in wild black and captive white abalone on the northeastern Pacific coast, but disease resistance levels are species-, and possibly population-specific. This study compared gene expression patterns in the digestive gland of Ca. Xc-exposed and unexposed (control) Pinto abalone (Haliotis kamtschatkana), a particularly susceptible species. Lab-induced Ca. Xc infections were followed over 7 months and RNAseq was used to identify differential gene expression. Exposed Pinto abalone showed distinct changes in expression of 68 genes at 3 and 7 months post-infection relative to those in control animals. Upregulation of an orexin-like receptor (which is involved in feeding signaling) and a zinc peptidase-like region (many amino peptidases are zinc peptidases) in animals infected for 7 months indicates that animals with Ca. Xc infection may be starving and upregulating processes associated with feeding and digestion. Other groups of differentially expressed genes (DEGs) were upregulated or downregulated across control and exposed individuals over the 7-month experiment, including DEG groups that likely correspond to early disease state and to general stress response of being held in captivity. No patterns emerged in genes known to be involved in molluscan immune response, despite this being an expectation during a 7-month infection; digestion-related genes and unannotated DEGs were identified as targets for future research on potential immune response to WS in abalone.

The Role of Anti-Viral Effector Molecules in Mollusc Hemolymph

Molluscs are major contributors to the international and Australian aquaculture industries, however, their immune systems remain poorly understood due to limited access to draft genomes and evidence of divergences from model organisms. As invertebrates, molluscs lack adaptive immune systems or 'memory', and rely solely on innate immunity for antimicrobial defence. Hemolymph, the circulatory fluid of invertebrates, contains hemocytes which secrete effector molecules with immune regulatory functions. Interactions between mollusc effector molecules and bacterial and fungal pathogens have been well documented, however, there is limited knowledge of their roles against viruses, which cause high mortality and significant production losses in these species. Of the major effector molecules, only the direct acting protein dicer-2 and the antimicrobial peptides (AMPs) hemocyanin and myticin-C have shown antiviral activity. A better understanding of these effector molecules may allow for the manipulation of mollusc proteomes to enhance antiviral and overall antimicrobial defence to prevent future outbreaks and minimize economic outbreaks. Moreover, effector molecule research may yield the description and production of novel antimicrobial treatments for a broad host range of animal species.

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

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.

Marine mollusc extracts-Potential source of SARS-CoV-2 antivirals

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) is a novel human coronavirus and the causative agent of coronavirus disease 2019 (Covid‐19). There is an urgent need for effective antivirals to treat current Covid‐19 cases and protect those unable to be vaccinated against SARS‐CoV‐2. Marine molluscs live in an environment containing high virus densities (>10⁷ virus particles per ml), and there are an estimated 100,000 species in the phylum Mollusca, demonstrating the success of their innate immune system. Mollusc‐derived antivirals are yet to be used clinically despite the activity of many extracts, including against human viruses, being demonstrated in vitro. Hemolymph of the Pacific oyster (Crassostrea gigas) has in vitro antiviral activity against herpes simplex virus and human adenovirus, while antiviral action against SARS‐CoV‐2 has been proposed by in silico studies. Such evidence suggests that molluscs, and in particular C. gigas hemolymph, may represent a source of antivirals for human coronaviruses.

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.

Genomic analysis and functional characterization of immune genes from the RIG-I- and MAVS-mediated antiviral signaling pathway in lamprey

Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) are well-known viral RNA sensors in the cytoplasm. RIG-I-mediated antiviral signals are activated by interacting with the adapter protein mitochondrial antiviral signaling (MAVS), which triggers interferon (IFN) responses via a signaling cascade. Although the complete RIG-I receptor signaling pathway has been traced back to teleosts, definitive evidence of its presence in lampreys is lacking. Here, we identified 13 pivotal molecules in the RIG-I signaling pathway in lamprey, and demonstrated that the original RIG-I/MAVS signaling pathway was activated and mediated the expression of unique immunity factors such as RRP4, to inhibit viral proliferation after viral infection in vivo and in vitro. This study confirmed the conservation of the RIG-I pathway, and the uniqueness of the RRP4 effector molecule in lamprey, and further clarified the evolutionary process of the RIG-I antiviral signaling pathway, providing evidence on the origins of innate antiviral immunity in vertebrates.

Functional characterization of the cnidarian antiviral immune response reveals ancestral complexity

Animals evolved a broad repertoire of innate immune sensors and downstream effector cascades for defense against RNA viruses. Yet, this system varies greatly among different bilaterian animals, masking its ancestral state. In this study, we aimed to characterize the antiviral immune response of the cnidarian Nematostella vectensis and decipher the function of the retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) known to detect viral double-stranded RNA (dsRNA) in bilaterians but activate different antiviral pathways in vertebrates and nematodes. We show that polyinosinic:polycytidylic acid (poly(I:C)), a mimic of long viral dsRNA and a primary ligand for the vertebrate RLR melanoma differentiation-associated protein 5 (MDA5), triggers a complex antiviral immune response bearing features distinctive for both vertebrate and invertebrate systems. Importantly, a well-characterized agonist of the vertebrate RIG-I receptor does not induce a significant transcriptomic response that bears signature of the antiviral immune response, which experimentally supports the results of a phylogenetic analysis indicating clustering of the two N. vectensis RLR paralogs (NveRLRa and NveRLRb) with MDA5. Furthermore, the results of affinity assays reveal that NveRLRb binds poly(I:C) and long dsRNA and its knockdown impairs the expression of putative downstream effector genes including RNA interference components. Our study provides for the first time the functional evidence for the conserved role of RLRs in initiating immune response to dsRNA that originated before the cnidarian–bilaterian split and lay a strong foundation for future research on the evolution of the immune responses to RNA viruses.

Type II Grass Carp Reovirus Infects Leukocytes but Not Erythrocytes and Thrombocytes in Grass Carp (Ctenopharyngodon idella)

Grass carp reovirus (GCRV) causes serious losses to the grass carp industry. At present, infectious tissues of GCRV have been studied, but target cells remain unclear. In this study, peripheral blood cells were isolated, cultured, and infected with GCRV. Using quantitative real-time polymerase chain reaction (qRT-PCR), Western Blot, indirect immunofluorescence, flow cytometry, and transmission electron microscopy observation, a model of GCRV infected blood cells in vitro was established. The experimental results showed GCRV could be detectable in leukocytes only, while erythrocytes and thrombocytes could not. The virus particles in leukocytes are wrapped by empty membrane vesicles that resemble phagocytic vesicles. The empty membrane vesicles of leukocytes are different from virus inclusion bodies in C. idella kidney (CIK) cells. Meanwhile, the expression levels of IFN1, IL-1β, Mx2, TNFα were significantly up-regulated in leukocytes, indicating that GCRV could cause the production of the related immune responses. Therefore, GCRV can infect leukocytes in vitro, but not infect erythrocytes and thrombocytes. Leukocytes are target cells in blood cells of GCRV infections. This study lays a theoretical foundation for the study of the GCRV infection mechanism and anti-GCRV immunity.

The primitive interferon-like system and its antiviral function in molluscs

The phylum mollusca is a very important group in the animal kingdom for the large number and diversified species. Recently, interest in molluscan immunity has increased due to their phylogenetic position and importance in worldwide aquaculture and aquatic environment. As the main aquaculture animal, most molluscs live in the water environment and they have to cope with many pathogen challenges, in which virus is one of the primary causes for the mass mortality. In vertebrates, interferon (IFN) system is generally recognized as the first line of defence against viral infection, while the antiviral mechanisms in molluscs remain to be clearly illuminated. Recently, some IFN-like proteins and IFN-related components have been characterized from molluscs, such as pattern recognition receptors (PRRs), interferon regulatory factors (IRFs), IFN-like receptors, JAK/STAT and IFN-stimulated genes (ISGs), which reinforce the existence of IFN-like system in molluscs. This system can be activated by virus or poly (I:C) challenges and further regulate the antiviral response of haemocytes in molluscs. This review summarizes the research progresses of IFN-like system in molluscs with the emphases on the uniformity and heterogeneity of IFN-like system of molluscs compared to that of other animals, which will be helpful for elucidating the antiviral modulation in molluscs and understanding the origin and evolution of IFN system.

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.

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.

The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia

BACKGROUND

Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa.

RESULTS

Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors.

CONCLUSIONS

Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.

Single-molecule long-read (SMRT) transcriptome sequencing of Mercenaria mercenaria reveals a powerful anti-apoptotic system critical for air exposure endurance

Mercenaria mercenaria is an economically important clam species and exhibits an outstanding resistance to multiple environmental stressors. However, our understanding of their stress adaptability is limited due to a lack of genomic information, such as transcriptome resources. In this study, single-molecule long-read (SMRT) mRNA sequencing was performed to obtain the full-length gill transcriptome reference sequences of M. mercenaria under air exposure stress. In all, 14.5 G subreads were obtained and assembled into 64,603 unigenes, among which 50,613 were successfully annotated. Additionally, 56,295 SSRs, 1457 transcription factors, and 5924 lncRNAs were identified in M. mercenaria transcriptome. Furthermore, numerous apoptosis-related transcripts were identified according to Swiss-Prot annotation and their numbers were counted. We also found that most apoptosis-related transcripts exhibited typical domains of a certain protein family through conserved domain prediction. Additionally, eight typical sequences related to apoptosis pathway were detected by RT-PCR, with the aim to show the sequential variation of gene expression levels under air exposure. These results implied that the complicated apoptosis system, especially the powerful anti-apoptotic system was critical for M. mercenaria to endure air exposure.

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.

Cellular and molecular complementary immune stress markers for the model species Dreissena polymorpha

This study aimed to combine cellular and molecular analyses for better detail the effects of various stresses on a sentinel species of freshwater invertebrate. For this purpose, the hemocytes of the zebra mussel, Dreissena polymorpha, were exposed to different stresses at two different intensities, high or low: chemical (cadmium and ionomycin), physical (ultraviolet B), or biological ones (Cryptosporidium parvum and Toxoplasma gondii). After exposure, flow cytometry and droplet digital PCR analyses were performed on the same pools of hemocytes. Several responses related to necrosis, apoptosis, phagocytosis, production of nitric oxide and expression level of several genes related to the antioxidant, detoxification and immune systems were evaluated. Results showed that hemocyte integrity was compromised by both chemical and physical stress, and cellular markers of phagocytosis reacted to ionomycin and protozoa. While cadmium induced oxidative stress and necrosis, ionomycin tends to modulate the immune response of hemocytes. Although both biological stresses led to a similar immune response, C. parvum oocysts induced more effects than T. gondii, notably through the expression of effector caspases gene and an increase in hemocyte necrosis. This suggests different management of the two protozoa by the cell. This work provides new knowledge of biomarkers in the zebra mussel, at both cellular and molecular levels, and contributes to elucidate the mechanisms of action of different kinds of stress in this species.

The noncoding small RNA SsrA is released by Vibrio fischeri and modulates critical host responses

The regulatory noncoding small RNAs (sRNAs) of bacteria are key elements influencing gene expression; however, there has been little evidence that beneficial bacteria use these molecules to communicate with their animal hosts. We report here that the bacterial sRNA SsrA plays an essential role in the light-organ symbiosis between Vibrio fischeri and the squid Euprymna scolopes. The symbionts load SsrA into outer membrane vesicles, which are transported specifically into the epithelial cells surrounding the symbiont population in the light organ. Although an SsrA-deletion mutant (ΔssrA) colonized the host to a normal level after 24 h, it produced only 2/10 the luminescence per bacterium, and its persistence began to decline by 48 h. The host's response to colonization by the ΔssrA strain was also abnormal: the epithelial cells underwent premature swelling, and host robustness was reduced. Most notably, when colonized by the ΔssrA strain, the light organ differentially up-regulated 10 genes, including several encoding heightened immune-function or antimicrobial activities. This study reveals the potential for a bacterial symbiont's sRNAs not only to control its own activities but also to trigger critical responses promoting homeostasis in its host. In the absence of this communication, there are dramatic fitness consequences for both partners.

Contrasting Immunomodulatory Effects of Probiotic and Pathogenic Bacteria on Eastern Oyster, Crassostrea Virginica, Larvae

Several Vibrio spp. cause acute and severe mortality events in hatcheries where larvae of bivalve mollusks are reared, potentially leading to subsequent shortage of bivalve seed for the grow-out industry. In particular, strains of Vibrio coralliilyticus have been identified as a major cause of disease in Pacific, Crassostrea gigas, and eastern, C. virginica, oyster hatcheries in the United States of America. Probiotic bacteria are an inexpensive, practical, and natural method of disease control. Previous research shows that pretreatment of larval oysters with probiotic bacteria Bacillus pumilus RI06-95 (RI) and Phaeobacter inhibens S4 (S4) significantly decreases mortality caused by experimental challenge with the bacterial pathogen V. coralliilyticus RE22 (RE22). This study aims to characterize the immune response of 6-10-day-old eastern oyster larvae to experimental challenge with pathogen V. coralliilyticus RE22 and probionts RI and S4. Treatments included (a) pathogen and probiont exposure at a concentration of 5 × 10⁴ CFU per mL (~2500 bacterial cells per larva) for a duration of 6 h, (b) probiont exposure at the same concentration for a duration of 24 h, and (c) probiont RI daily treatment of larvae in the hatchery for 4, 11, and 15 days. Differential gene expression analysis compared pathogen or probiotic-treated transcriptomes to unexposed controls. Probiotic and pathogen treatment led to upregulation of transcripts coding for several immune pattern recognition receptors (PRRs) involved in environmental sensing and detection of microbes in oyster larvae. Larval oyster responses to pathogen RE22 suggested suppression of expression of genes in immune signaling pathways (myd88, tak1, nkap), failure in upregulation of immune effector genes, high metabolic demand, and oxidative stress that potentially contributed to mortality. On the other hand, the transcriptomic response to probiotic bacteria RI and S4 suggested activation of immune signaling pathways and expression of immune effectors (e.g., Cv-spi2, mucins and perforin-2). These key features of the host immune response to probiotic bacteria were shared despite the length of probiotic exposure, probiotic species, and the type of environment in which exposures were conducted. This study suggests that pre-exposure of eastern oyster larvae to probiotics for 6-24 h prior to pathogenic challenge leads to a robust and effective immune response that may contribute to protecting larvae from subsequent challenge with V. coralliilyticus RE22. This research provides new insights into host-microbe interactions in larval oysters that could be applied in the management of vibriosis in bivalve hatcheries.

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.

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.

Immune Control of Herpesvirus Infection in Molluscs

Molluscan herpesviruses that are capable of infecting economically important species of abalone and oysters have caused significant losses in production due to the high mortality rate of infected animals. Current methods in preventing and controlling herpesviruses in the aquacultural industry are based around biosecurity measures which are impractical and do not contain the virus as farms source their water from oceans. Due to the lack of an adaptive immune system in molluscs, vaccine related therapies are not a viable option; therefore, a novel preventative strategy known as immune priming was recently explored. Immune priming has been shown to provide direct protection in oysters from Ostreid herpesvirus-1, as well as to their progeny through trans-generational immune priming. The mechanisms of these processes are not completely understood, however advancements in the characterisation of the oyster immune response has assisted in formulating potential hypotheses. Limited literature has explored the immune response of abalone infected with Haliotid herpesvirus as well as the potential for immune priming in these species, therefore, more research is required in this area to determine whether this is a practical solution for control of molluscan herpesviruses in an aquaculture setting.

CgMyD88s Serves as an Innate Immune System Plug During Ostreid Herpesvirus 1 Infection in the Pacific Oyster (Crassostrea gigas)

Ostreid herpesvirus-1 microvariant (OsHV-1 μVar) is considered a major infectious microbe that can reduce the survival of natural or cultured oysters in summer. Because they lack an adaptive immune system, oysters are dependent on their innate immune systems to fight pathogens. The duplication and functional divergence of innate immune genes in the oyster have been studied, but the contribution of molecular mechanisms underlying innate immunity remains to be defined. Here, we identified the interacting proteins associated with Crassostrea gigas Toll-like receptors (CgTLR) using a yeast two-hybrid (Y2H) screening system. A total of eight proteins were identified that could interact with CgTLR. Three of these appeared at least four times in the screening and were related to MyD88. Two genes encoding these MyD88-like proteins, CgMyD88-1 and CgMyD88-2, possessed typical death and TIR domains. The third gene encoding an MyD88-like protein possessed only a TIR domain, and we named it CgMyD88s. CgMyD88s interacted only with CgTLR, but not CgMyD88-1 or CgMyD88-2. Both CgMyD88-1 and CgMyD88-2 mRNAs were upregulated after OsHV-1 μVar infection, whereas the expression of CgMyD88s decreased. When overexpressed in HEK293T cells, CgMyD88-1 and CgMyD88-2 activated an NF-κB reporter, whereas CgMyD88s impaired activation induced by CgMyD88-1 or CgMyD88-2. Intriguingly, the silencing of CgMyD88s using double-stranded RNA (dsRNA)-mediated RNA interference increased the expression of CgMyD88-1 and CgMyD88-2. Taken together, our results revealed that CgMyD88-1, CgMyD88-2, and CgMyD88s may all participate in the TLR-mediated innate immune pathway and that CgMyD88s served as a plug to avoid oysters from excessive inflammatory response during OsHV-1 μVar infections.

A Sustained Immune Response Supports Long-Term Antiviral Immune Priming in the Pacific Oyster, Crassostrea gigas

Over the last decade, innate immune priming has been evidenced in many invertebrate phyla. If mechanistic models have been proposed, molecular studies aiming to substantiate these models have remained scarce. We reveal here the transcriptional signature associated with immune priming in the oyster Crassostrea gigas Oysters were fully protected against Ostreid herpesvirus 1 (OsHV-1), a major oyster pathogen, after priming with poly(I·C), which mimics viral double-stranded RNA. Global analysis through RNA sequencing of oyster and viral genes after immune priming and viral infection revealed that poly(I·C) induces a strong antiviral response that impairs OsHV-1 replication. Protection is based on a sustained upregulation of immune genes, notably genes involved in the interferon pathway and apoptosis, which control subsequent viral infection. This persistent antiviral alert state remains active over 4 months and supports antiviral protection in the long term. This acquired resistance mechanism reinforces the molecular foundations of the sustained response model of immune priming. It further opens the way to applications (pseudovaccination) to cope with a recurrent disease that causes dramatic economic losses in the shellfish farming industry worldwide.IMPORTANCE In the last decade, important discoveries have shown that resistance to reinfection can be achieved without a functional adaptive immune system, introducing the concept of innate immune memory in invertebrates. However, this field has been constrained by the limited number of molecular mechanisms evidenced to support these phenomena. Taking advantage of an invertebrate species, the Pacific oyster (Crassostrea gigas), in which we evidenced one of the longest and most effective periods of protection against viral infection observed in an invertebrate, we provide the first comprehensive transcriptomic analysis of antiviral innate immune priming. We show that priming with poly(I·C) induced a massive upregulation of immune-related genes, which control subsequent viral infection, and it was maintained for over 4 months after priming. This acquired resistant mechanism reinforces the molecular foundations of the sustained response model of immune priming. It opens the way to pseudovaccination to prevent the recurrent diseases that currently afflict economically or ecologically important invertebrates.

Variation in global transcriptomic response to Perkinsus marinus infection among eastern oyster families highlights potential mechanisms of disease resistance

Dermo disease, caused by the protozoan parasite Perkinsus marinus, negatively impacts wild and cultured Eastern oyster populations, yet our knowledge of the mechanistic bases for parasite pathogenicity and the Eastern oyster's response to it is limited. To better understand host responses to the parasite and identify molecular mechanisms underlying disease-resistance phenotypes, we experimentally challenged two families exhibiting divergent Dermo-resistance phenotypes with the parasite, generated global expression profiles using RNAseq and identified differentially expressed transcripts between control and challenged oysters from each family at multiple time points post-parasite injection. The susceptible and resistant families exhibited strikingly different transcriptomic responses to the parasite over a 28-day time period. The resistant family exhibited a strong, focused, early response to P. marinus infection, where many significantly upregulated transcripts were associated with the biological processes "regulation of proteolysis" and "oxidation-reduction process." P. marinus virulence factors are mainly comprised of proteases that facilitate parasite invasion and weaken host humoral defenses, thus host upregulation of transcripts associated with negative regulation of proteolysis is consistent with a Dermo-resistant phenotype. In contrast, the susceptible family mounted a very weak, disorganized, initial response to the parasite. Few transcripts were differentially expressed between control and injected oysters, and no functional enrichment was detected among them. At the final 28 d time point 2450 differentially expressed transcripts were identified and were associated with either "G-protein coupled receptor activity" (upregulated) or "microtubule-based process" (downregulated). A handful of protease inhibitors were differentially expressed between control and injected susceptible oysters, but this function was not enriched in the susceptible data set. The differential expression patterns observed in this study provide valuable insight into the functional basis of Dermo resistance and suggest that the timing of expression is just as important as the transcripts being expressed.

Massive expansion and diversity of nicotinic acetylcholine receptors in lophotrochozoans

Background

Nicotinic acetylcholine receptors (nAChRs) are among the oldest and most conserved transmembrane receptors involved in signal transduction. Despite the prevalence and significance of cholinergic signaling, the diversity and evolution of nAChRs are not fully understood.

Result

By comparative genomic analysis, we found massive expansions of nAChR genes in molluscs and some other lophotrochozoans. The expansion is particularly pronounced in stationary bivalve molluscs with simple nervous systems, with the number of nAChR genes ranging from 99 to 217 in five bivalves, compared with 10 to 29 in five ecdysozoans and vertebrates. The expanded molluscan nAChR genes tend to be intronless and in tandem arrays due to retroposition followed by tandem duplication. Phylogenetic analysis revealed diverse nAChR families in the common ancestor of bilaterians, which subsequently experienced lineage-specific expansions or contractions. The expanded molluscan nAChR genes are highly diverse in sequence, domain structure, temporal and spatial expression profiles, implying diversified functions. Some molluscan nAChR genes are expressed in early development before the development of the nervous system, while others are involved in immune and stress responses.

Conclusion

The massive expansion and diversification of nAChR genes in bivalve molluscs may be a compensation for reduced nervous systems as part of adaptation to stationary life under dynamic environments, while in vertebrates a subset of specialized nAChRs are retained to work with advanced nervous systems. The unprecedented diversity identified in molluscs broadens our view on the evolution and function of nAChRs that are critical to animal physiology and human health.

Oyster Versatile IKKα/βs Are Involved in Toll-Like Receptor and RIG-I-Like Receptor Signaling for Innate Immune Response

IκB kinases (IKKs) play critical roles in innate immunity through signal-induced activation of the key transcription factors nuclear factor-κB (NF-κB) and interferon regulatory factors (IRFs). However, studies of invertebrate IKK functions remain scarce. In this study, we performed phylogenetic analysis of IKKs and IKK-related kinases encoded in the Pacific oyster genome. We then cloned and characterized the oyster IKKα/β-2 gene. We found that oyster IKKα/β-2, a homolog of human IKKα/IKKβ, responded to challenge with lipopolysaccharide (LPS), peptidoglycan (PGN), and polyinosinic-polycytidylic acid [poly(I:C)]. As a versatile immune molecule, IKKα/β-2 activated the promoters of NF-κB, TNFα, and IFNβ, as well as IFN-stimulated response element (ISRE)-containing promoters, initiating an antibacterial or antiviral immune state in mammalian cells. Importantly, together with the cloned oyster IKKα/β-1, we investigated the signal transduction pathways mediated by these two IKKα/β proteins. Our results showed that IKKα/β-1 and IKKα/β-2 could interact with the oyster TNF receptor-associated factor 6 (TRAF6) and that IKKα/β-2 could also bind to the oyster myeloid differentiation factor 88 (MyD88) protein directly, suggesting that oyster IKKα/βs participate in both RIG-I-like receptor (RLR) and Toll-like receptor (TLR) signaling for the reception of upstream immune signals. The fact that IKKα/β-1 and IKKα/β-2 formed homodimers by interacting with themselves and heterodimers by interacting with each other, along with the fact that both oyster IKKα/β proteins interacted with NEMO protein, indicates that oyster IKKα/βs and the scaffold protein NEMO form an IKK complex, which may be a key step in phosphorylating IκB proteins and activating NF-κB. Moreover, we found that oyster IKKα/βs could interact with IRF8, and this may be related to the IKK-mediated activation of ISRE promotors and their involvement in the oyster "interferon (IFN)-like" antiviral pathway. Moreover, the expression of oyster IKKα/β-1 and IKKα/β-2 may induce the phosphorylation of IκB proteins to activate NF-κB. These results reveal the immune function of oyster IKKα/β-2 and establish the existence of mollusk TLR and RLR signaling mediated by IKKα/β proteins for the first time. Our findings should be helpful in deciphering the immune mechanisms of invertebrates and understanding the development of the vertebrate innate immunity network.