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

Ostreid herpesvirus 1 latent infection and reactivation in adult Pacific oysters, Crassostrea gigas

Ostreid herpesvirus 1 (OsHV-1) is one of the most economically important pathogens of Pacific oysters. Understanding the pathogenesis of this virus is critical to developing tools to control outbreaks on shellfish farms. OsHV-1 is genetically related to vertebrate herpesviruses, which have a lytic and a latent stage, with the latent stage capable of being reactivated to the lytic stage. Here, OsHV-1 latency in Pacific oysters was investigated in experimentally and naturally infected oysters. Lytic infection in one-year-old oysters injected with the Tomales Bay strain of OsHV-1 was detectable between 1 and 4 days post-injection (dpi) but was not detectable after 5 dpi. The injected oysters shed 1 × 102 to 1 × 104 DNA copies/ml into the water during the 4-day acute phase. Lytic shedding was not detectable in two-year-old oysters injected similarly with the same strain of OsHV-1; however, the OsHV-1 genome was detectable by qPCR in the adductor muscle, gill, mantle, and hemocytes within the first 3 dpi, after which it became undetectable. No OsHV-1 was detectable in the adductor muscle, gill, or mantle from experimentally infected oysters on days 15 and 21 post-injection or from oysters sampled 9 months after surviving an OsHV-1 mortality event; however, OsHV-1 DNA could be detected in hemocytes of both experimentally infected oysters at 21 dpi and naturally infected oysters using nested PCR. In addition, lytic viral gene transcription was detectable in hemocytes of experimentally infected oysters between 1 and 21 dpi and in hemocytes of naturally infected oysters. Furthermore, OsHV-1 reactivation from latency was induced in experimentally infected oysters at 21 dpi and in naturally infected oysters 12 months after an OsHV-1 outbreak.

A common garden experiment supports a genetic component underlying the increased resilience of common cockle (Cerastoderma edule) to the parasite Marteilia cochillia

The common cockle is a valuable bivalve species inhabiting the Atlantic European coasts. The parasite Marteilia cochillia has devastated cockle beds in the southern Galician (NW Spain) rias since 2012. Previous data suggested that cockles from Ría de Arousa acquired some resilience to this parasite through natural selection after consecutive annual marteiliosis outbreaks and candidate markers associated with marteiliosis resilience were identified using population genomics and transcriptomics approaches. Here, a common garden experiment was performed using a naïve stock (from Ría de Muros-Noia) and an affected stock (from Ría de Arousa) to test this hypothesis. Breeders from both stocks were used to produce seed cohorts at hatchery, which were pre-grown in a raft (outdoor nursery stage) and deployed in two shellfish beds affected by marteiliosis in Ría de Arousa (growing-out stage). In both beds, the naïve stock showed high marteiliosis prevalence and was fully depleted in a short period, while the affected stock barely showed evidence of marteiliosis. A set of 45 SNPs putatively associated with marteiliosis resilience were fitted for MassARRAY genotyping to check their role in the differential resilience detected between both stocks. Though no significant differentiation was found between the naïve and the affected stocks with neutral markers, 28 SNPs showed significant divergence between them, suggesting that these SNPs were involved in directional selection during eight generations (to the most) of marteiliosis pressure (long-term selection). Furthermore, signals of selection were also detected in the naïve stock along the marteiliosis outbreak in the growing-out stage (short-term selection) and six SNPs, all shared with the long-term evaluation, showed consistent signals of differentiation according to the infection severity. Some of these SNPs were located within immune genes pertaining to families such as proteasome, ubiquitin, tumor necrosis factor, and glutathione S-transferase. These resilience-associated markers will be useful to recover cockle production in Galicia.

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.

A Deterministic Model for Understanding Nonlinear Viral Dynamics in Oysters

Contamination of oysters with a variety of viruses is one key pathway to trigger outbreaks of massive oyster mortality as well as human illnesses, including gastroenteritis and hepatitis. Much effort has gone into examining the fate of viruses in contaminated oysters, yet the current state of knowledge of nonlinear virus-oyster interactions is not comprehensive because most studies have focused on a limited number of processes under a narrow range of experimental conditions. A framework is needed for describing the complex nonlinear virus-oyster interactions. Here, we introduce a mathematical model that includes key processes for viral dynamics in oysters, such as oyster filtration, viral replication, the antiviral immune response, apoptosis, autophagy, and selective accumulation. We evaluate the model performance for two groups of viruses, those that replicate in oysters (e.g., ostreid herpesvirus) and those that do not (e.g., norovirus), and show that this model simulates well the viral dynamics in oysters for both groups. The model analytically explains experimental findings and predicts how changes in different physiological processes and environmental conditions nonlinearly affect in-host viral dynamics, for example, that oysters at higher temperatures may be more resistant to infection by ostreid herpesvirus. It also provides new insight into food treatment for controlling outbreaks, for example, that depuration for reducing norovirus levels is more effective in environments where oyster filtration rates are higher. This study provides the foundation of a modeling framework to guide future experiments and numerical modeling for better prediction and management of outbreaks.

IMPORTANCE The fate of viruses in contaminated oysters has received a significant amount of attention in the fields of oyster aquaculture, food quality control, and public health. However, intensive studies through laboratory experiments and in situ observations are often conducted under a narrow range of experimental conditions and for a specific purpose in their respective fields. Given the complex interactions of various processes and nonlinear viral responses to changes in physiological and environmental conditions, a theoretical framework fully describing the viral dynamics in oysters is warranted to guide future studies from a top-down design. Here, we developed a process-based, in-host modeling framework that builds a bridge for better communications between different disciplines studying virus-oyster interactions.

Little evidence for genetic variation associated with susceptibility to sea star wasting syndrome in the keystone species Pisaster ochraceus

The keystone species Pisaster ochraceus suffered mass mortalities along the northeast Pacific Ocean from Sea Star Wasting Syndrome (SSWS) outbreaks in 2013-2016. SSWS causation remains of debate, leading to concerns as to whether outbreaks will continue to impact this species. Considering the apparent link between ocean temperature and SSWS, the future of this species and intertidal communities remains uncertain. Surveys of co-occurring apparently normal and wasting P. ochraceus along the central Oregon coast in 2016 allowed us to address whether variation in disease status showed genetic variation that may be associated with differences in susceptibility to SSWS. We performed restriction site-associated DNA sequencing (2bRAD-seq) to genotype ~72,000 single nucleotide polymorphism (SNP) loci across apparently normal and wasting sea stars. Locus-specific analyses of differentiation (F_ST ) between disease-status groups revealed no signal of genetic differences separating the two groups. Using a multivariate approach, we observed weak separation between the groups, but identified 18 SNP loci showing highest discriminatory power between the groups and scanned the genome annotation for linked genes. A total of 34 protein-coding genes were found to be located within 15 kb (measured by linkage disequilibrium decay) of at least one of the 18 SNPs, and 30 of these genes had homologies to annotated protein databases. Our results suggest that the likelihood of developing SSWS symptoms does not have a strong genetic basis. The few genomic regions highlighted had only modest levels of differentiation, but the genes associated with these regions may form the basis for functional studies aiming to understand disease progression.

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.