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

Bacterial and viral co-infections in aquaculture under climate warming: co-evolutionary implications, diagnosis, and treatment

Climate change and the associated environmental temperature fluctuations are contributing to increases in the frequency and severity of disease outbreaks in both wild and farmed aquatic species. This has a significant impact on biodiversity and also puts global food production systems, such as aquaculture, at risk. Most infections are the result of complex interactions between multiple pathogens, and understanding these interactions and their co-evolutionary mechanisms is crucial for developing effective diagnosis and control strategies. In this review, we discuss current knowledge on bacteria-bacteria, virus-virus, and bacterial and viral co-infections in aquaculture as well as their co-evolution in the context of global warming. We also propose a framework and different novel methods (e.g. advanced molecular tools such as digital PCR and next-generation sequencing) to (1) precisely identify overlooked co-infections, (2) gain an understanding of the co-infection dynamics and mechanisms by knowing species interactions, and (3) facilitate the development multi-pathogen preventive measures such as polyvalent vaccines. As aquaculture disease outbreaks are forecasted to increase both due to the intensification of practices to meet the protein demand of the increasing global population and as a result of global warming, understanding and treating co-infections in aquatic species has important implications for global food security and the economy.

Unveiling CRESS DNA Virus Diversity in Oysters by Virome

Oysters that filter feed can accumulate numerous pathogens, including viruses, which can serve as a valuable viral repository. As oyster farming becomes more prevalent, concerns are mounting about diseases that can harm both cultivated and wild oysters. Unfortunately, there is a lack of research on the viruses and other factors that can cause illness in shellfish. This means that it is harder to find ways to prevent these diseases and protect the oysters. This is part of a previously started project, the Dataset of Oyster Virome, in which we further study 30 almost complete genomes of oyster-associated CRESS DNA viruses. The replication-associated proteins and capsid proteins found in CRESS DNA viruses display varying evolutionary rates and frequently undergo recombination. Additionally, some CRESS DNA viruses have the capability for cross-species transmission. A plethora of unclassified CRESS DNA viruses are detectable in transcriptome libraries, exhibiting higher levels of transcriptional activity than those found in metagenome libraries. The study significantly enhances our understanding of the diversity of oyster-associated CRESS DNA viruses, emphasizing the widespread presence of CRESS DNA viruses in the natural environment and the substantial portion of CRESS DNA viruses that remain unidentified. This study's findings provide a basis for further research on the biological and ecological roles of viruses in oysters and their environment.

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.

Understanding the Dynamic of POMS Infection and the Role of Microbiota Composition in the Survival of Pacific Oysters, Crassostrea gigas

For over a decade, Pacific oyster mortality syndrome (POMS), a polymicrobial disease, induced recurring episodes of massive mortality affecting Crassostrea gigas oysters worldwide. Recent studies evidenced a combined infection of the ostreid herpesvirus (OsHV-1 μVar) and opportunistic bacteria in affected oysters. However, the role of the oyster microbiota in POMS is not fully understood. While some bacteria can protect hosts from infection, even minor changes to the microbial communities may also facilitate infection and worsen disease severity. Using a laboratory-based experimental infection model, we challenged juveniles from 10 biparental oyster families with previously established contrasted genetically based ability to survive POMS in the field. Combining molecular analyses and 16S rRNA gene sequencing with histopathological observations, we described the temporal kinetics of POMS and characterized the changes in microbiota during infection. By associating the microbiota composition with oyster mortality rate, viral load, and viral gene expression, we were able to identify both potentially harmful and beneficial bacterial amplicon sequence variants (ASVs). We also observed a delay in viral infection resulting in a later onset of mortality in oysters compared to previous observations and a lack of evidence of fatal dysbiosis in infected oysters. Overall, these results provide new insights into how the oyster microbiome may influence POMS disease outcomes and open new perspectives on the use of microbiome composition as a complementary screening tool to determine shellfish health and potentially predict oyster vulnerability to POMS. IMPORTANCE For more than a decade, Pacific oyster mortality syndrome (POMS) has severely impacted the Crassostrea gigas aquaculture industry, at times killing up to 100% of young farmed Pacific oysters, a key commercial species that is cultivated globally. These disease outbreaks have caused major financial losses for the oyster aquaculture industry. Selective breeding has improved disease resistance in oysters, but some levels of mortality persist, and additional knowledge of the disease progression and pathogenicity is needed to develop complementary mitigation strategies. In this holistic study, we identified some potentially harmful and beneficial bacteria that can influence the outcome of the disease. These results will contribute to advance disease management and aquaculture practices by improving our understanding of the mechanisms behind genetic resistance to POMS and assisting in predicting oyster vulnerability to POMS.

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.

Characterization of Host Cell Potential Proteins Interacting with OsHV-1 Membrane Proteins

The interaction between viral membrane associate proteins and host cellular surface molecules should facilitate the attachment and entry of OsHV-1 into host cells. Thus, blocking the putative membrane proteins ORF25 and ORF72 of OsHV-1 with antibodies that have previously been reported to subdue OsHV-1 replication in host cells, especially ORF25. In this study, prey proteins in host hemocytes were screened by pull-down assay with recombinant baits ORF25 and ORF72, respectively. Gene Ontology (GO) analysis of these prey proteins revealed that most of them were mainly associated with binding, structural molecule activity and transport activity in the molecular function category. The protein–protein interaction (PPI) network of the prey proteins was constructed by STRING and clustered via K-means. For both ORF25 and ORF72, three clusters of these prey proteins were distinguished that were mainly associated with cytoskeleton assembly, energy metabolism and nucleic acid processing. ORF25 tended to function in synergy with actins, while ORF72 functioned mainly with tubulins. The above results suggest that these two putative membrane proteins, ORF25 and ORF72, might serve a role in the transport of viral particles with the aid of a cytoskeleton inside cells.

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.