Analysis of clinical ostreid herpesvirus 1 (Malacoherpesviridae) specimens by sequencing amplified fragments from three virus genome areas

Propidium monoazide PCR, a method to determine OsHV-1 undamaged capsids and to estimate virus Lethal Dose 50

Ostreid herpes virus 1 (OsHV-1) has been classified within the Malacoherpesviridae family from the Herpesvirales order. OsHV-1 is the etiological agent of a contagious viral disease of Pacific oysters, C. gigas, affecting also other bivalve species. Mortality rates reported associated with the viral infection vary considerably between sites and countries and depend on the age of affected stocks. A variant called μVar has been reported since 2008 in Europe and other variants in Australia and in New Zealand last decade. These variants are considered as the main causative agents of mass mortality events affecting C. gigas. Presently there is no established cell line that allows for the detection of infectious OsHV-1. In this context, a technique of propidium monoazide (PMA) PCR was developed in order to quantify "undamaged" capsids. This methodology is of interest to explore the virus infectivity. Being able to quantify viral particles getting an undamaged capsid (not only an amount of viral DNA) in tissue homogenates prepared from infected oysters or in seawater samples can assist in the definition of a Lethal Dose (LD) 50 and gain information in the experiments conducted to reproduce the viral infection. The main objectives of the present study were (i) the development/optimization of a PMA PCR technique for OsHV-1 detection using the best quantity of PMA and verifying its effectiveness through heat treatment, (ii) the definition of the percentage of undamaged capsids in four different tissue homogenates prepared from infected Pacific oysters and (iii) the approach of a LD50 during experimental viral infection assays on the basis of a number of undamaged capsids. Although the developped PMA PCR technique was unable to determine OsHV-1 infectivity in viral supensions, it could greatly improve interpretation of virus positive results obtained by qPCR. This technique is not intended to replace the quantification of viral DNA by qPCR, but it does make it possible to give a form of biological meaning to the detection of this DNA.

A core of functional complementary bacteria infects oysters in Pacific Oyster Mortality Syndrome

BACKGROUND

The Pacific oyster Crassostrea gigas is one of the main cultivated invertebrate species worldwide. Since 2008, oyster juveniles have been confronted with a lethal syndrome known as the Pacific Oyster Mortality Syndrome (POMS). POMS is a polymicrobial disease initiated by a primary infection with the herpesvirus OsHV-1 µVar that creates an oyster immunocompromised state and evolves towards a secondary fatal bacteremia.

RESULTS

In the present article, we describe the implementation of an unprecedented combination of metabarcoding and metatranscriptomic approaches to show that the sequence of events in POMS pathogenesis is conserved across infectious environments. We also identified a core bacterial consortium which, together with OsHV-1 µVar, forms the POMS pathobiota. This bacterial consortium is characterized by high transcriptional activities and complementary metabolic functions to exploit host's resources. A significant metabolic specificity was highlighted at the bacterial genus level, suggesting low competition for nutrients between members of the core bacteria.

CONCLUSIONS

Lack of metabolic competition between the core bacteria might favor complementary colonization of host tissues and contribute to the conservation of the POMS pathobiota across distinct infectious environments.

Diversity and molecular epidemiology of Ostreid herpesvirus 1 in farmed Crassostrea gigas in Australia: Geographic clusters and implications for "microvariants" in global mortality events

Since 2010, mass mortality events known as Pacific oyster mortality syndrome (POMS) have occurred in Crassostrea gigas in Australia associated with Ostreid herpesvirus 1. The virus was thought to be an OsHV-1 µVar or "microvariant", i.e. one of the dominant variants associated with POMS in Europe, but there are few data to characterize the genotype in Australia. Consequently, the genetic identity and diversity of the virus was determined to understand the epidemiology of the disease in Australia. Samples were analysed from diseased C. gigas over five summer seasons between 2011 and 2016 in POMS-affected estuaries: Georges River in New South Wales (NSW), Hawkesbury River (NSW) and Pitt Water in Tasmania. Sequencing was attempted for six genomic regions. Numerous variants were identified among these regions (n = 100 isolates) while twelve variants were identified from concatenated nucleotide sequences (n = 61 isolates). Nucleotide diversity of the seven genotypes of C region among Australian isolates (Pi 0.99 × 10-3) was the lowest globally. All Australian isolates grouped in a cluster distinct from other OsHV-1 isolates worldwide. This is the first report that Australian outbreaks of POMS were associated with OsHV-1 distinct from OsHV-1 reference genotype, µVar and other microvariants from other countries. The findings illustrate that microvariants are not the only variants of OsHV-1 associated with mass mortality events in C. gigas. In addition, there was mutually exclusive spatial clustering of viral genomic and amino acid sequence variants between estuaries, and a possible association between genotype/amino acid sequence and the prevalence and severity of POMS, as this differed between these estuaries. The sequencing findings supported prior epidemiological evidence for environmental reservoirs of OsHV-1 for POMS outbreaks in Australia.

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.

Gene delivery available in molluscan cells by strong promoter discovered from bivalve-infectious virus

Understanding gene functions in marine invertebrates has been limited, largely due to the lack of suitable assay systems. Such a system requires investigative methods that are reproducible and can be quantitatively evaluated, such as a cell line, and a strong promoter that can drive high expression of a transgene. In this study, we established primary cell culture from a marine bivalve mollusc, Mizuhopecten yessoensis. Using scallop primary cells, we optimized electroporation conditions for transfection and carried out a luciferase-based promoter activity assay to identify strong promoter sequences that can drive expression of a gene of interest. We evaluated potential promoter sequences from genes of endogenous and exogenous origin and discovered a strong viral promoter derived from a bivalve-infectious virus, ostreid herpesvirus-1 (OsHV-1). This promoter, we termed OsHV-1 promoter, showed 24.7-fold and 16.1-fold higher activity than the cytomegalovirus immediate early (CMV IE) promoter and the endogenous EF1α promoter, the two most commonly used promoters in bivalves so far. Our GFP assays showed that the OsHV-1 promoter is active not only in scallop cells but also in HEK293 cells and zebrafish embryos. The OsHV-1 promoter practically enables functional analysis of marine molluscan genes, which can contribute to unveiling gene-regulatory networks underlying astonishing regeneration, adaptation, reproduction, and aging in marine invertebrates.

Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: A first attempt to phylogeographic inference for a marine mollusc disease

The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific oyster mortality syndrome (a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled twenty-one OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis, and ancestral state reconstruction of discrete traits approaches, we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes-Oléron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

Aquatic pathogens and biofouling: pilot study of ostreid herpesvirus 1 translocation by bivalves

Ostreid herpesvirus 1 (OsHV-1) has caused mass mortalities in Pacific oysters (Crassostrea gigas) in Europe, Australia, and New Zealand. While aquaculture-associated movements of infected Pacific oysters are a well-known cause of OsHV-1 spread once established in a region, translocation via biofouling of aquaculture equipment or vessels needs further investigation to explain the more distant spread of OsHV-1. Laboratory experiments were designed to test for transmission of OsHV-1 between infected and naïve Pacific oysters via a simulated biofouling translocation scenario. Three common biofouling species [Sydney rock oysters (Saccostrea glomerata), Mediterranean mussels (Mytilus galloprovincialis) and Pacific oysters] were tested as intermediaries using a cohabitation challenge with Pacific oysters infected by injection. Transmission occurred, albeit for one of eight replicates when Pacific oysters were the intermediary species. This demonstrated a possible pathway for pathogen spread via biofouling containing Pacific oysters while highlighting the complexity of OsHV-1 transmission. Such complexities require further investigation to inform future risk assessments and management of fouled aquaculture equipment and vessels.

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 μVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses.

The first detection of a novel OsHV-1 microvariant in San Diego, California, USA

The spread, emergence, and adaptation of pathogens causing marine disease has been problematic to fisheries and aquaculture industries for the last several decades creating the need for strategic management and biosecurity practices. The Pacific oyster (Crassostrea gigas), a highly productive species globally, has been a target of disease and mortality caused by a viral pathogen, the Ostreid herpesvirus 1 (OsHV-1) and its microvariants (OsHV-1 µvars). During routine surveillance to establish health history at a shellfish aquaculture nursery system in San Diego, California, the presence of OsHV-1 in Pacific oyster juveniles was detected. Quantification of OsHV-1 in tissues of oysters revealed OsHV-1 viral loads > 10⁶ copies/mg. We characterized and identified the OsHV-1 variant by sequencing of ORFs 4 (C2/C6) and 43 (IA1/IA2), which demonstrated that this variant is a novel OsHV-1 microvariant: OsHV-1 µvar SD. A pilot transmission study indicates that OsHV-1 µvar SD is infectious with high viral loads ~ 7.57 × 10⁶ copies/mg detected in dead individuals. The detection of OsHV-1 µvar SD in a large port mirrors previous studies conducted in Australia where aquaculture farms and feral populations near port locations may be at a higher risk of OsHV-1 emergence. Further research is needed to understand the impacts of OsHV-1 µvar SD, such as transmission studies focusing on potential vectors and characterization of virulence as compared to other OsHV-1 µvars. To increase biosecurity of the global aquaculture industry, active and passive surveillance may be necessary to reduce spread of pathogens and make appropriate management decisions.

Reduction in Virulence over Time in Ostreid herpesvirus 1 (OsHV-1) Microvariants between 2011 and 2015 in Australia

Microvariant genotypes of Ostreid herpesvirus 1 (OsHV-1) are associated with mass mortality events of Pacific oysters in many countries. The OsHV-1 microvariant (µVar) emerged in France 2008 and caused significant economic losses as it became endemic and displaced the previously dominant OsHV-1 reference genotype. Recently, considerable genotypic variation has been described for OsHV-1 microvariants, however, less is known about variation in viral phenotype. This study used an in vivo laboratory infection model to assess differences in total cumulative mortality, peak viral load, transmissibility, and dose-response for three OsHV-1 isolates obtained between 2011 and 2015 from endemic waterways in Australia. This followed field observations of apparent reductions in the severity of mass mortalities over this time. Significantly higher hazard of death and cumulative mortality were observed for an isolate obtained in 2011 compared to isolates from 2014–2015. In keeping with other studies, the hazard of death was higher in oysters challenged by injection compared to challenge by cohabitation and the mortality was higher when the initial dose was 1 × 10⁴ OsHV-1 DNA copies per oyster injection compared to 1 × 10² DNA copies. There was no difference in the quantity of OsHV-1 DNA at time of death that could be related to isolate or dose, suggesting similar pathogenetic processes in the individual oysters that succumbed to end-stage disease. While the isolates examined in this study were biased towards pathogenic types of OsHV-1, as they were collected during disease outbreaks, the variation in virulence that was observed, when combined with prior data on subclinical infections, suggests that surveillance for low virulence genotypes of OsHV-1 would be rewarding. This may lead to new approaches to disease management which utilize controlled exposure to attenuated strains of OsHV-1.

The Pacific Oyster Mortality Syndrome, a Polymicrobial and Multifactorial Disease: State of Knowledge and Future Directions

The Pacific oyster (Crassostreae gigas) has been introduced from Asia to numerous countries around the world during the 20th century. C. gigas is the main oyster species farmed worldwide and represents more than 98% of oyster production. The severity of disease outbreaks that affect C. gigas, which primarily impact juvenile oysters, has increased dramatically since 2008. The most prevalent disease, Pacific oyster mortality syndrome (POMS), has become panzootic and represents a threat to the oyster industry. Recently, major steps towards understanding POMS have been achieved through integrative molecular approaches. These studies demonstrated that infection by Ostreid herpesvirus type 1 µVar (OsHV-1 µvar) is the first critical step in the infectious process and leads to an immunocompromised state by altering hemocyte physiology. This is followed by dysbiosis of the microbiota, which leads to a secondary colonization by opportunistic bacterial pathogens, which in turn results in oyster death. Host and environmental factors (e.g. oyster genetics and age, temperature, food availability, and microbiota) have been shown to influence POMS permissiveness. However, we still do not understand the mechanisms by which these different factors control disease expression. The present review discusses current knowledge of this polymicrobial and multifactorial disease process and explores the research avenues that must be investigated to fully elucidate the complexity of POMS. These discoveries will help in decision-making and will facilitate the development of tools and applied innovations for the sustainable and integrated management of oyster aquaculture.

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.

Differential Mortality and High Viral Load in Naive Pacific Oyster Families Exposed to OsHV-1 Suggests Tolerance Rather than Resistance to Infection

Pacific oysters, Crassostrea gigas, are one of the most productive aquaculture species in the world. However, they are threatened by the spread of Ostreid herpesvirus-1 (OsHV-1) and its microvariants (collectively "µvars"), which cause mass mortalities in all life stages of Pacific oysters globally. Breeding programs have been successful in reducing mortality due to OsHV-1 variants following viral outbreaks; however, an OsHV-1-resistant oyster line does not yet exist in the United States (US), and it is unknown how OsHV-1 µvars will affect US oyster populations compared to the current variant, which is similar to the OsHV-1 reference, found in Tomales Bay, CA. The goals of this study were to investigate the resistance of C. gigas juveniles produced by the Molluscan Broodstock Program (MBP) to three variants of OsHV-1: a California reference OsHV-1, an Australian µvar, and a French µvar. This is the first study to directly compare OsHV-1 µvars to a non-µvar. The survival probability of oysters exposed to the French (FRA) or Australian (AUS) µvar was significantly lower (43% and 71%, respectively) than to the reference variant and controls (96%). No oyster family demonstrated resistance to all three OsHV-1 variants, and many surviving oysters contained high copy numbers of viral DNA (mean ~3.53 × 108). These results indicate that the introduction of OsHV-1 µvars could have substantial effects on US Pacific oyster aquaculture if truly resistant lines are not achieved, and highlight the need to consider resistance to infection in addition to survival as traits in breeding programs to reduce the risk of the spread of OsHV-1 variants.

Unraveling concordant and varying responses of oyster species to Ostreid Herpesvirus 1 variants

The Ostreid herpesvirus 1 (OsHV-1) and variants, particularly the microvariants (μVars), are virulent and economically devastating viruses impacting oysters. Since 2008 OsHV-1 μVars have emerged rapidly having particularly damaging effects on aquaculture industries in Europe, Australia and New Zealand. We conducted field trials in Tomales Bay (TB), California where a non-μVar strain of OsHV-1 is established and demonstrated differential mortality of naturally exposed seed of three stocks of Pacific oyster, Crassostrea gigas, and one stock of Kumamoto oyster, C. sikamea. Oysters exposed in the field experienced differential mortality that ranged from 64 to 99% in Pacific oysters (Tasmania>Midori = Willapa stocks), which was much higher than that of Kumamoto oysters (25%). Injection trials were done using French (FRA) and Australian (AUS) μVars with the same oyster stocks as planted in the field and, in addition, two stocks of the Eastern oyster, C. virginica. No mortality was observed in control oysters. One C. virginica stock suffered ~10% mortality when challenged with both μVars tested. Two Pacific oyster stocks suffered 75 to 90% mortality, while one C. gigas stock had relatively low mortality when challenged with the AUS μVar (~22%) and higher mortality when challenged with the French μVar (~72%). Conversely, C. sikamea suffered lower mortality when challenged with the French μVar (~22%) and higher mortality with the AUS μVar (~44%). All dead oysters had higher viral loads (~1000×) as measured by quantitative PCR relative to those that survived. However, some survivors had high levels of virus, including those from species with lower mortality. Field mortality in TB correlated with laboratory mortality of the FRA μVar (69% correlation) but not with that of the AUS μVar, which also lacked correlation with the FRA μVar. The variation in response to OsHV-1 variant challenges by oyster species and stocks demonstrates the need for empirical assessment of multiple OsHV-1 variants.

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.

Contribution of Viral Genomic Diversity to Oyster Susceptibility in the Pacific Oyster Mortality Syndrome

Juvenile Pacific oysters (Crassostrea gigas) are subjected to recurrent episodes of mass mortalities that constitute a threat for the oyster industry. This mortality syndrome named “Pacific Oyster Mortality Syndrome” (POMS) is a polymicrobial disease whose pathogenesis is initiated by a primary infection by a variant of an Ostreid herpes virus named OsHV-1 μVar. The characterization of the OsHV-1 genome during different disease outbreaks occurring in different geographic areas has revealed the existence of a genomic diversity for OsHV-1 μVar. However, the biological significance of this diversity is still poorly understood. To go further in understanding the consequences of OsHV-1 diversity on POMS, we challenged five biparental families of oysters to two different infectious environments on the French coasts (Atlantic and Mediterranean). We observed that the susceptibility to POMS can be different among families within the same environment but also for the same family between the two environments. Viral diversity analysis revealed that Atlantic and Mediterranean POMS are caused by two distinct viral populations. Moreover, we observed that different oyster families are infected by distinct viral populations within a same infectious environment. Altogether these results suggest that the co-evolutionary processes at play between OsHV-1 μVar and oyster populations have selected a viral diversity that could facilitate the infection process and the transmission in oyster populations. These new data must be taken into account in the development of novel selective breeding programs better adapted to the oyster culture environment.

Oyster hemolymph is a complex and dynamic ecosystem hosting bacteria, protists and viruses

BACKGROUND

The impact of the microbiota on host fitness has so far mainly been demonstrated for the bacterial microbiome. We know much less about host-associated protist and viral communities, largely due to technical issues. However, all microorganisms within a microbiome potentially interact with each other as well as with the host and the environment, therefore likely affecting the host health.

RESULTS

We set out to explore how environmental and host factors shape the composition and diversity of bacterial, protist and viral microbial communities in the Pacific oyster hemolymph, both in health and disease. To do so, five oyster families differing in susceptibility to the Pacific oyster mortality syndrome were reared in hatchery and transplanted into a natural environment either before or during a disease outbreak. Using metabarcoding and shotgun metagenomics, we demonstrate that hemolymph can be considered as an ecological niche hosting bacterial, protist and viral communities, each of them shaped by different factors and distinct from the corresponding communities in the surrounding seawater. Overall, we found that hemolymph microbiota is more strongly shaped by the environment than by host genetic background. Co-occurrence network analyses suggest a disruption of the microbial network after transplantation into natural environment during both non-infectious and infectious periods. Whereas we could not identify a common microbial community signature for healthy animals, OsHV-1 μVar virus dominated the hemolymph virome during the disease outbreak, without significant modifications of other microbiota components.

CONCLUSION

Our study shows that oyster hemolymph is a complex ecosystem containing diverse bacteria, protists and viruses, whose composition and dynamics are primarily determined by the environment. However, all of these are also shaped by oyster genetic backgrounds, indicating they indeed interact with the oyster host and are therefore not only of transient character. Although it seems that the three microbiome components respond independently to environmental conditions, better characterization of hemolymph-associated viruses could change this picture.

Climate Change and Bivalve Mass Mortality in Temperate Regions

One of the fastest-growing global food sectors is the bivalve aquaculture industry. Bivalves particularly oysters, mussels and clams are important sources of animal protein (Tan and Ransangan 2016a, b). Bivalve aquaculture represents 14-16% of the average per capita animal protein for 1.5 billion people and supports over 200,000 livelihoods, mostly in developing countries (FAO 2018). Most of the bivalves produced around the world (89%) are from aquaculture (FAO 2016). To date, mollusc aquaculture have accounted for 21.42% (17.14 million tonnes) of the total aquaculture production, with Asia being the largest contributor (92.27%) (FAO 2018).

IKKε-like plays an important role in the innate immune signaling of the Pacific oyster (Crassostrea gigas)

IκB-related kinase ε (IKKε) plays a crucial role in the activation of nuclear factor κB (NF-κB) by phosphorylating inhibitor of NF-κB (IκB) and in the regulation of interferon (IFN) gene expression by phosphorylating IFN regulatory factors (IRFs). In this study, we cloned an IKKε homologue cDNA (designated as CgIKKε-like) from the Pacific oyster, Crassostrea gigas. The full 2896-bp cDNA sequence comprised a 2163-bp open reading frame (ORF) encoding 720 amino acids. CgIKKε-like is ubiquitously expressed, and its mRNA levels in hemocytes after poly I:C, V. alginolyticus, or OsHV-1 μVar challenge were analyzed by real-time PCR. Compared to that in the control, CgIKKε-like mRNA expression levels were significantly increased at 3 h and peaked at 6 h after OsHV-1 μVar challenge; no obvious changes were observed in expression levels until 24 h after either V. alginolyticus or poly I:C challenge, reaching a maximum at 24 h (p < 0.01) and then rapidly decreasing. CgIKKε-like transfection into human cell lines induced NF-κB and ISRE activation, while transfection with CgIKKε-like deletion mutants abolished NF-κB and ISRE reporter gene activation. Additionally, CgIKKε-like could interact with CgTBK1 and could form homodimers strongly, which may be critical for the immune signaling transduction. Last but not least, we found that CgIKKε-like may increase CgIκBs phosphorylation and could interact with CgIRF8. Together, these results suggest that CgIKKε-like could respond to pathogenic infection, participate in the immune signal transduction and activate NF-κB and ISRE reporter genes. Thus, CgIKKε-like could play an important role in the oyster immune system.

Fine-scale temporal dynamics of herpes virus and vibrios in seawater during a polymicrobial infection in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is currently being impacted by a polymicrobial disease that involves early viral infection by ostreid herpesvirus-1 (OsHV-1) followed by a secondary bacterial infection leading to death. A widely used method of inducing infection consists of placing specific pathogen-free oysters ('recipients') in cohabitation in the laboratory with diseased oysters that were naturally infected in the field ('donors'). With this method, we evaluated the temporal dynamics of pathogen release in seawater and the cohabitation time necessary for disease transmission and expression. We showed that OsHV-1 and Vibrio spp. in the seawater peaked concomitantly during the first 48 h and decreased thereafter. We found that 1.5 h of cohabitation with donors was enough time to transmit pathogens to recipients and to induce mortality later, reflecting the highly contagious nature of the disease. Finally, mortality of recipients was associated with increasing cohabitation time with donors until reaching a plateau at 20%. This reflects the cumulative effect of exposure to pathogens. The optimal cohabitation time was 5-6 d, the mortality of recipients occurring 1-2 d earlier.

A-to-I editing of Malacoherpesviridae RNAs supports the antiviral role of ADAR1 in mollusks

Background

Adenosine deaminase enzymes of the ADAR family are conserved in metazoans. They convert adenine into inosine in dsRNAs and thus alter both structural properties and the coding potential of their substrates. Acting on exogenous dsRNAs, ADAR1 exerts a pro- or anti-viral role in vertebrates and Drosophila.

Results

We traced 4 ADAR homologs in 14 lophotrochozoan genomes and we classified them into ADAD, ADAR1 or ADAR2, based on phylogenetic and structural analyses of the enzymatic domain. Using RNA-seq and quantitative real time PCR we demonstrated the upregulation of one ADAR1 homolog in the bivalve Crassostrea gigas and in the gastropod Haliotis diversicolor supertexta during Ostreid herpesvirus-1 or Haliotid herpesvirus-1 infection. Accordingly, we demonstrated an extensive ADAR-mediated editing of viral RNAs. Single nucleotide variation (SNV) profiles obtained by pairing RNA- and DNA-seq data from the viral infected individuals resulted to be mostly compatible with ADAR-mediated A-to-I editing (up to 97%). SNVs occurred at low frequency in genomic hotspots, denoted by the overlapping of viral genes encoded on opposite DNA strands. The SNV sites and their upstream neighbor nucleotide indicated the targeting of selected adenosines. The analysis of viral sequences suggested that, under the pressure of the ADAR editing, the two Malacoherpesviridae genomes have evolved to reduce the number of deamination targets.

Conclusions

We report, for the first time, evidence of an extensive editing of Malacoherpesviridae RNAs attributable to host ADAR1 enzymes. The analysis of base neighbor preferences, structural features and expression profiles of molluscan ADAR1 supports the conservation of the enzyme function among metazoans and further suggested that ADAR1 exerts an antiviral role in mollusks.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1472-6) contains supplementary material, which is available to authorized users.

Influence of environment on the pathogenesis of Ostreid herpesvirus-1 (OsHV-1) infections in Pacific oysters (Crassostrea gigas) through differential microbiome responses

The oyster microbiome is thought to contribute to the pathogenesis of mass mortality disease in Pacific oysters, associated with OsHV-1. As filter-feeders, oysters host a microbiota that can be influenced by the estuarine environment. This may alter susceptibility to OsHV-1 infections, causing variable mortality. This study aimed at: (1) differences in the microbiome of Pacific oysters with a common origin but grown in geographically distinct estuaries; (2) evaluating changes occurring in the microbiota, especially in Vibrio, and (3) differential responses of the oyster microbiome, in response to an OsHV-1 infection. Pacific oysters sourced from a single hatchery but raised separately in Patonga Creek, Shoalhaven River and Clyde River of NSW, Australia, were used and challenged with OsHV-1. The initial microbiome composition was different in the three batches and changed further, post-injection (p < 0.05). The Patonga oysters with the highest mortality also had higher OsHV-1 and Vibrio quantities compared to the other two batches (p < 0.05). The higher initial bacterial diversity in Patonga oysters decreased in moribund oysters which was not observed in the other two batches (p < 0.05). The microbiome of survivors of OsHV-1 infection and negative control oysters of two batches, did not show any changes with the relevant pre-challenged microbiome. A strong correlation was observed between the OsHV-1 and Vibrio quantities in OsHV-1 infected oysters (r = 0.6; p < 0.001). In conclusion, the Pacific oyster microbiome differed in different batches despite a common hatchery origin. Different microbiomes responded differently with a differential outcome of OsHV-1 challenge. The higher Vibrio load in oysters with higher OsHV-1 content and higher mortality, suggests a role in Vibrio in the pathogenesis of this mortality disease. This study provided insights of the potential of different estuarine environments to shape the Pacific oyster microbiome and how different microbiomes are associated with different outcomes of OsHV-1 infection.

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.

Mass mortality in Korean bay scallop (Argopecten irradians) associated with Ostreid Herpesvirus-1 μVar

Since November 2017, mass mortalities of larvae of bay scallop (Argopecten irradians) were reported in hatcheries located at the southern area of Republic of Korea. Over 90% of larvae aged 5-10 days sank to the bottom of the tank and died. The hatcheries could not produce spat, and thus artificial seed production industry incurred huge losses. We identified Ostreid Herpesvirus-1 μVar (OsHV-1 μVar) associated with mass mortality by PCR, sequencing and transmission electron microscopy (TEM). All the samples were positive for OsHV-1 μVar with 99% sequence identity to previously reported OsHV-1 μVar sequences. Partial sequence of ORF-4 of OsHV-1 detected in this study was more closely related to sequences isolated from Europe. This is the first report to confirm the mortality caused by an OsHV-1 infection in the bay scallop.

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.

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.

Immune-suppression by OsHV-1 viral infection causes fatal bacteraemia in Pacific oysters

Infectious diseases are mostly explored using reductionist approaches despite repeated evidence showing them to be strongly influenced by numerous interacting host and environmental factors. Many diseases with a complex aetiology therefore remain misunderstood. By developing a holistic approach to tackle the complexity of interactions, we decipher the complex intra-host interactions underlying Pacific oyster mortality syndrome affecting juveniles of Crassostrea gigas, the main oyster species exploited worldwide. Using experimental infections reproducing the natural route of infection and combining thorough molecular analyses of oyster families with contrasted susceptibilities, we demonstrate that the disease is caused by multiple infection with an initial and necessary step of infection of oyster haemocytes by the Ostreid herpesvirus OsHV-1 µVar. Viral replication leads to the host entering an immune-compromised state, evolving towards subsequent bacteraemia by opportunistic bacteria. We propose the application of our integrative approach to decipher other multifactorial diseases that affect non-model species worldwide.

Pacific oyster mortality syndrome is a poorly understood cause of mortality in commercially important oyster species. Here, the authors use multiple infection experiments to show that the syndrome is caused by sequential infection by herpesvirus and opportunistic bacteria.

Identification of a newly described OsHV-1 µvar from the North Adriatic Sea (Italy)

The surveillance activities for abnormal bivalve mortality events in Italy include the diagnosis of ostreid herpesvirus type 1 (OsHV-1) in symptomatic oysters. OsHV-1-positive oysters (Crassostrea gigas) were used as a source for in vivo virus propagation and a virus-rich sample was selected to perform shotgun sequencing based on Illumina technology. Starting from this unpurified supernatant sample from gills and mantle, we generated 3.5 million reads (2×300 bp) and de novo assembled the whole genome of an Italian OsHV-1 microvariant (OsHV-1-PT). The OsHV-1-PT genome encodes 125 putative ORFs, 7 of which had not previously been predicted in other sequenced Malacoherpesviridae. Overall, OsHV-1-PT displays typical microvariant OsHV-1 genome features, while few polymorphisms (0.08 %) determine its uniqueness. As little is known about the genetic determinants of OsHV-1 virulence, comparing complete OsHV-1 genomes supports a better understanding of the virus pathogenicity and provides new insights into virus-host interactions.

Real-time quantitative isothermal detection of Ostreid herpesvirus-1 DNA in Scapharca subcrenata using recombinase polymerase amplification

Ostreid herpesvirus-1 (OsHV-1) is a well-known pathogen associated with high mortality rates in hatchery-reared larvae and juveniles of different bivalve species worldwide. Early, rapid and accurate diagnosis plays a fundamental role in disease prevention and control in aquaculture. Recombinase polymerase amplification (RPA) is a novel isothermal amplification method, which can amplify detectable amount of DNA at 37 °C-39 °C within 20 min. In the present study, two sets of specific primers and probes were designed for the real-time quantitative RPA (qRPA) detection of OsHV-1 DNA. The sensitivity and specificity of detection were evaluated by comparison with quantitative polymerase chain reaction (qPCR). The detection limit for qRPA assays was shown to be 5 copies DNA/reaction for the primer set ORF95, which was lower than the 100 copies required for the qPCR test. The optimal reaction temperature and time were 37 °C for 20 min, making this approach faster than qPCR. This is the first study to apply qPCR and qRPA methods to detect OsHV-1 in Scapharca subcrenata. The percentage of viral load sample detected by the two methods was 22% and the correlation of the two virus quantitative results was 0.8. Therefore, qRPA assays is sensitive, fast, and high-temperature independent relative to qPCR and is suitable for critical clinical diagnostics use and rapid field analysis in resource-limited settings.

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.

Mass mortalities of unknown aetiology in Pacific oysters Crassostrea gigas in Port Stephens, New South Wales, Australia

From January to June 2013 and November to January 2014, mass mortalities were reported in Pacific oysters Crassostrea gigas cultivated in Port Stephens estuary, New South Wales, Australia. In some cases, 100% mortality was reported in both triploid and diploid C. gigas, although native species of oyster cultivated in the same areas remained unaffected. Histological examination of oysters sampled from the time of mortality events revealed consistent but non-specific pathology, involving a diffuse haemocytic infiltrate in the connective tissue surrounding the digestive gland, extending into the mantle in some instances, but no other signs of any infectious aetiological agent. We conducted a structured survey in early January 2014 to compare samples of C. gigas from affected and unaffected areas by bacteriology and histopathology. Quantitative PCR excluded involvement of ostreid herpesvirus-1 (OsHV-1) in these mortality events. To determine whether a directly transmissible aetiological agent was responsible for the mortalities, naïve C. gigas sourced from an estuary where no evidence of mortality was reported were challenged with material derived from affected oysters. Significant mortality was only observed in naïve C. gigas directly inoculated with purified cultures of Vibrio spp. isolated from affected oysters, but this could not be replicated by cohabitation with naïve C. gigas. Analysis of environmental data indicated that mortality events generally coincided with periods of low salinity and high temperature. The results from this study suggest that the cause of the mortality events was multifactorial in nature and not due to any single directly transmissible aetiological agent.

Oyster RNA-seq Data Support the Development of Malacoherpesviridae Genomics

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

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

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

Complete genome sequence of Ostreid herpesvirus type 1 µVar isolated during mortality events in the Pacific oyster Crassostrea gigas in France and Ireland

Infections with Ostreid herpesvirus 1 (OsHV-1) microvariants in young Pacific oysters are associated with massive mortality events and significant economic losses. Previous studies, focusing on few regions of the genome, have revealed the genomic diversity of these genotypes with respect to the reference type. We used a NGS process to sequence the whole genome of the OsHV-1 µVar in infected individuals, collected during mortality events in France and Ireland. The final genome length of OsHV-1 µVar was approximately 205kbp, shorter than the reference genotype and the overall genome organisation resembled herpes simplex viruses. 94.4% similarity was observed with the OsHV-1 reference genotype. Large indels, including five deletions and three insertions were found to induce the loss and the addition of several ORFs, summed with codon substitutions in 64% of genes shared with the reference type. This diversity raises the question of the exact origin and evolution of OsHV-1 µVar.

Haemocytes collected from experimentally infected Pacific oysters, Crassostrea gigas: Detection of ostreid herpesvirus 1 DNA, RNA, and proteins in relation with inhibition of apoptosis

Recent transcriptomic approaches focused on anti-viral immunity in molluscs lead to the assumption that the innate immune system, such as apoptosis, plays a crucial role against ostreid herpesvirus type 1 (OsHV-1), infecting Pacific cupped oyster, Crassostrea gigas. Apoptosis constitutes a major mechanism of anti-viral response by limiting viral spread and eliminating infected cells. In this way, an OsHV-1 challenge was performed and oysters were monitored at three times post injection to investigate viral infection and host response: 2h (early after viral injection in the adductor muscle), 24h (intermediate time), and 48h (just before first oyster mortality record). Virus infection, associated with high cumulative mortality rates (74% and 100%), was demonstrated in haemocytes by combining several detection techniques such as real-time PCR, real-time RT PCR, immunofluorescence assay, and transmission electron microscopy examination. High viral DNA amounts ranged from 5.46×104 to 3.68×105 DNA copies ng-1 of total DNA, were detected in dead oysters and an increase of viral transcripts was observed from 2, 24, and 48hpi for the five targeted OsHV-1 genes encoding three putative membrane proteins (ORFs 25, 41, and 72), a putative dUTPase (ORF 75), and a putative apoptosis inhibitor (ORF 87). Apoptosis was studied at molecular and cellular levels with an early marker (phosphatidyl-serine externalisation measured by flow cytometry and epifluorescence microscopy) and a later parameter (DNA fragmentation by terminal deoxynucleotidyltransferase-mediated dUTP nick end labeling assay (TUNEL)). The down-regulation of genes encoding proteins involved in the activation of the apoptotic pathway (TNF and caspase 3) and the up-regulation of genes encoding anti-apoptotic proteins (IAP-2, and Bcl-2) suggested an important anti-apoptosis phenomenon in haemocytes from OsHV-1 infected oysters at 24 and 48hpi. Additionally, more phosphatidyl-serines were externalized and more cells with DNA fragmentation were observed in haemocytes collected from artificial seawater injected oysters than in haemocytes collected from OsHV-1 infected oysters at 24 and 48hpi, suggesting an inhibition of the apoptotic process in presence of the virus. In conclusion, this study is the first to focus on C. gigas haemocytes, cells involved in the host immune defense, during an OsHV-1 challenge in controlled conditions by combining various and original approaches to investigate apoptosis at molecular and cellular levels.

Viruses infecting marine molluscs

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

Long dsRNAs promote an anti-viral response in Pacific oyster hampering ostreid herpesvirus 1 replication

Double stranded RNA-mediated genetic interference (RNAi) is a widely used reverse genetic tool for determining the loss-of-function phenotype of a gene. Here, the possible induction of an immune response by long dsRNA was tested in a marine bivalve, i.e. Crassostrea gigas, as well as the specific role of the subunit 2 of the nuclear factor κB inhibitor (IκB2). This gene is a candidate of particular interest for functional investigations in the context of massive mortality oyster events as Cg-IκB2 mRNA levels exhibited significant variation depending on the amount of ostreid herpesvirus 1 (OsHV-1) DNA detected. In the present study, dsRNAs targeting Cg-IκB2 and Green Fluorescence Protein genes were injected in vivo into oysters before being challenged by OsHV-1. Survival appeared close to 100% in both dsRNA injected conditions associated with a low detection of viral DNA and a low expression of a panel of 39 OsHV-1 genes as compared to infected control. Long dsRNA molecules, both Cg-IκB2- and GFP-dsRNA, may have induced an anti-viral state controlling the OsHV-1 replication and precluding the understanding of the Cg-IκB2 specific role. Immune-related genes including Cg-IκB1, Cg-Rel1, Cg-IFI44, Cg-PKR, and Cg-IAP appeared activated in dsRNA-injected condition potentially hampering viral replication and thus conferring a better resistance to OsHV-1 infection. We revealed that long dsRNA-mediated genetic interference triggered an anti-viral state in the oyster, emphasizing the need of new reverse genetics tools for assessing immune gene function and avoiding off-target effects in bivalves.

In situ hybridization and histopathological observations during ostreid herpesvirus-1-associated mortalities in Pacific oysters Crassostrea gigas

In a previous longitudinal study conducted during a mortality investigation associated with ostreid herpesvirus-1 (OsHV-1) microvariant in New Zealand Pacific oysters in 2010-2011, temporality of OsHV-1 nucleic acid detection by real-time PCR assay and onset of Pacific oyster mortality was observed. The present study further elucidated the role of OsHV-1 using an in situ hybridization (ISH) assay on sections of Pacific oysters collected from the same longitudinal study. Hybridization of the labelled probe with the target region of the OsHV-1 genome in infected cells was detected colorimetrically using nitro blue tetrazolium (NBT). OsHV-1 presence and distribution in spat indicated by the ISH signal was then compared with the existence of pathological changes in oyster tissues. Dark blue to purplish black NBT cell labelling was seen predominantly in the stroma of the mantle and gills at Day 5 post introduction to the farm. The distribution and location of ISH signals indicated the extent of OsHV-1-infected cells in multiple tissues. Histopathological abnormalities were mostly non-specific; however, a progressive pattern of increasingly widespread haemocytosis coincided with the appearance of OsHV-1-infected cells in spat collected at different time-points. The visualisation of an increasing number of OsHV-1-positive cells in spat prior to a marked increase in mortality indicated the strong likelihood of an on-going and active viral infection in some oysters. Further studies are recommended to elucidate OsHV-1 pathogenesis in Pacific oysters in association with other potentially causal variables, such as elevated temperature and interaction with Vibrio spp. bacteria.

Occurrence of OsHV-1 in Crassostrea gigas Cultured in Ireland during an Exceptionally Warm Summer. Selection of Less Susceptible Oysters

The occurrence of OsHV-1, a herpes virus causing mass mortality in the Pacific oyster Crassostrea gigas was investigated with the aim to select individuals with different susceptibility to the infection. Naïve spat transferred to infected areas and juveniles currently being grown at those sites were analyzed using molecular and histology approaches. The survey period distinguishes itself by very warm temperatures reaching up to 3.5°C above the average. The virus was not detected in the virus free area although a spread of the disease could be expected due to high temperatures. Overall mortality, prevalence of infection and viral load was higher in spat confirming the higher susceptibility in early life stages. OsHV-1 and oyster mortality were detected in naïve spat after 15 days of cohabitation with infected animals. Although, infection was associated with mortality in spat, the high seawater temperatures could also be the direct cause of mortality at the warmest site. One stock of juveniles suffered an event of abnormal mortality that was significantly associated with OsHV-1 infection. Those animals were infected with a previously undescribed microvariant whereas the other stocks were infected with OsHV-1 μVar. Cell lesions due to the infection were observed by histology and true infections were corroborated by in situ hybridization. Survivors from the natural outbreak were exposed to OsHV-1 μVar by intramuscular injection and were compared to naïve animals. The survival rate in previously exposed animals was significantly higher than in naïve oysters. Results derived from this study allowed the selection of animals that might possess interesting characteristics for future analysis on OsHV-1 resistance.

TANK-binding kinase-1 broadly affects oyster immune response to bacteria and viruses

As a benthic filter feeder of estuaries, the immune system of oysters provides one of the best models for studying the genetic and molecular basis of the innate immune pathway in marine invertebrates and examining the influence of environmental factors on the immune system. Here, the molecular function of molluscan TANK-binding kinase-1 (TBK1) (which we named CgTBK1) was studied in the Pacific oyster, Crassostrea gigas. Compared with known TBK1 proteins in other model organisms, CgTBK1 contains a conserved S-TKc domain and a coiled coil domain at the N- and C-terminals but lacks an important ubiquitin domain. Quantitative real-time PCR analysis revealed that the expression level of CgTBK1 was ubiquitous in all selected tissues, with highest expression in the gills. CgTBK1 expression was significantly upregulated in response to infections with Vibrio alginolyticus, ostreid herpesvirus 1 (OsHV-1 reference strain and μvar), and polyinosinic:polycytidylic acid sodium salt, suggesting its broad function in immune response. Subcellular localization showed the presence of CgTBK1 in the cytoplasm of HeLa cells, suggesting its potential function as the signal transducer between the receptor and transcription factor. We further demonstrated that CgTBK1 interacted with CgSTING in HEK293T cells, providing evidence that CgTBK1 could be activated by direct binding to CgSTING. In summary, we characterized the TBK1 gene in C. gigas and demonstrated its role in the innate immune response to pathogen infections.

Ostreid herpesvirus in wild oysters from the Huelva coast (SW Spain)

This is the first report of ostreid herpesvirus 1 microvariant (OsHV-1 µVar) infecting natural oyster beds located in Huelva (SW Spain). The virus was detected in 3 oyster species present in the intertidal zone: Crassostrea gigas (Thunberg, 1793), C. angulata (Lamarck, 1819) and, for the first time, in Ostrea stentina Payraudeau, 1826. Oysters were identified by a specific polymerase chain reaction (PCR) and posterior restriction fragment length polymorphism (RFLP) analysis based on cytochrome oxidase I (COI) mitochondrial DNA. Results confirmed that C. angulata still remains the dominant oyster population in SW Spain despite the introduction of C. gigas for cultivation in the late 1970s, and its subsequent naturalization. C. angulata shows a higher haplotype diversity than C. gigas. OsHV-1 virus was detected by PCR with C2/C6 pair primers. Posterior RFLP analyses with the restriction enzyme MfeI were done in order to reveal the OsHV-1 µVar. Detections were confirmed by DNA sequencing, and infections were evidenced by in situ hybridization in C. gigas, C. angulata and O. stentina samples. The prevalence was similar among the 3 oyster species but varied between sampling locations, being higher in areas with greater harvesting activities. OsHV-1 µVar accounted for 93% of all OsHV-1 detected.

Ostreid herpesvirus type 1 genomic diversity in wild populations of Pacific oyster Crassostrea gigas from Italian coasts

Ostreid herpesvirus 1 (OsHV-1) is a significant pathogen affecting the young Pacific oyster Crassostrea gigas, worldwide. A new variant, OsHV-1 μVar, has been associated with recurrent mortality events in Europe since 2008. Epidemiological data collection is key for global risk assessment; however little is known about health status and genotypes present in European wild oyster beds. Most studies to date have involved only cultivated individuals during mortality events, and reported low genotype diversity. With this study, conducted along the Italian coasts, we investigated for the first time the presence of OsHV-1 in European natural oyster beds. Analysis of three genomic regions revealed the presence of at least nine different genotypes, including two variants close to the OsHV-1 reference, known since the early 1990s but with no European record reported since 2010, and highlights relevant genotype diversity in natural environment. Phylogenetic analysis distinguished two distinct clusters and geographical distribution of genotypes, with the exception of a variant very closely related to the μVar, which appeared the single genotype present in all the Adriatic sites. Interestingly, these wild symptom free populations could represent, in Europe, an accessible alternative to the import of OsHV-1-resistant oyster strains from the East Pacific, the native area of C. gigas, avoiding the high-risk of non-native marine species and new pathogen introductions.

Two Oyster Species That Show Differential Susceptibility to Virus Infection Also Show Differential Proteomic Responses to Generic dsRNA

Viral diseases are a significant cause of mortality and morbidity in oysters, resulting in significant economic losses. We investigated the proteomic responses of these two species of oysters to generic double-stranded RNAs (poly I:C and poly A:U). Analysis of proteomic data using isobaric tags for relative and absolute quantitaion (iTRAQ) indicated that there were significant differences in the proteomic responses of the two oyster species resulting from this treatment. Gene ontology analysis showed that several biological processes, cellular components, and molecular function were unique to the different data sets. For example, a number of proteins implicated in the TLR signaling pathway were associated with the Saccostrea glomerata data set but were absent in the Crassostra gigas data set. These results suggest that the differences in the proteomic responses to dsRNA may underpin the biological differences in viral susceptibility. Molecular targets previously shown to be expressed in C. gigas in response to OsHV1 infections were not present in our proteomic data sets, although they were present in the RNA extracted from the very same tissues. Taken together, our data indicate that there are substantial disparities between transcriptomic and proteomic responses to dsRNA challenge, and a comprehensive account of the oysters' biological responses to these treatments must take into account that disparity.

In situ localization and tissue distribution of ostreid herpesvirus 1 proteins in infected Pacific oyster, Crassostrea gigas

Immunohistochemistry (IHC) assays were conducted on paraffin sections from experimentally infected spat and unchallenged spat produced in hatchery to determine the tissue distribution of three viral proteins within the Pacific oyster, Crassostrea gigas. Polyclonal antibodies were produced from recombinant proteins corresponding to two putative membrane proteins and one putative apoptosis inhibitor encoded by ORF 25, 72, and 87, respectively. Results were then compared to those obtained by in situ hybridization performed on the same individuals, and showed a substantial agreement according to Landis and Koch numeric scale. Positive signals were mainly observed in connective tissue of gills, mantle, adductor muscle, heart, digestive gland, labial palps, and gonads of infected spat. Positive signals were also reported in digestive epithelia. However, few positive signals were also observed in healthy appearing oysters (unchallenged spat) and could be due to virus persistence after a primary infection. Cellular localization of staining seemed to be linked to the function of the viral protein targeted. A nucleus staining was preferentially observed with antibodies targeting the putative apoptosis inhibitor protein whereas a cytoplasmic localization was obtained using antibodies recognizing putative membrane proteins. The detection of viral proteins was often associated with histopathological changes previously reported during OsHV-1 infection by histology and transmission electron microscopy. Within the 6h after viral suspension injection, positive signals were almost at the maximal level with the three antibodies and all studied organs appeared infected at 28h post viral injection. Connective tissue appeared to be a privileged site for OsHV-1 replication even if positive signals were observed in the epithelium cells of different organs which may be interpreted as a hypothetical portal of entry or release for the virus. IHC constitutes a suited method for analyzing the early infection stages of OsHV-1 infection and a useful tool to investigate interactions between OsHV-1 and its host at a protein level.