Long-range PCR and high-throughput sequencing of Ostreid herpesvirus 1 indicate high genetic diversity and complex evolution process

Heritability of chronic thermal tolerance and genetic correlations with growth traits in the Pacific oyster (Crassostrea gigas)

The increasing seawater temperature during summer months frequently results in severe mortalities in the Pacific oyster Crassostrea gigas around the world, becoming one of the most significant problems challenging the oyster farming industry. In northern China, significant recurrent summer losses of C. gigas have occurred, and its impact on oyster aquaculture has increased in recent years. Selective breeding for improved oyster resistance to high temperature could help to reduce this massive mortality, but the extent of genetic variation underlying this trait is currently unknown. In this study, we constructed 38 full-sib families using the wild C. gigas and estimated the genetic parameters by performing two month-long high-temperature challenge experiments (30 ℃). Experiment 1 was performed in March 2022 followed by experiment 2 in June 2022 (spawning season). In both challenge experiments, there were significant differences in survival among families, suggesting that C. gigas has a different ability to survive under heat stress. Notably, significantly greater mortality was observed for experiment 2, which related to reproductive status and may contribute to additional stress. Thermal tolerance was defined using both binary test survival and time of death traits. Heritability estimates for thermal tolerance were low to moderate (0.16-0.36 for experiment 1 and 0.16-0.33 for experiment 2) using both a Bayesian (MCMCglmm) and a likelihood-based (ASReml-R) approach and estimated heritability of the threshold animal model using ASReml-R (0.16) appeared to be lower compared to MCMCglmm (0.31-0.32). Notably, the genetic and phenotypic correlations for thermal tolerance between two experiments were 0.463 (BS) to 0.491 (TD) and 0.510 (family survival), respectively, which suggested a significant re-ranking of the family breeding values in different time periods. Finally, the genetic and phenotypic correlations were low between growth traits (shell height, shell length, and shell width) and thermal tolerance, suggesting that selection for these traits should be conducted separately. This study reports the first estimation of genetic parameters for chronic thermal tolerance in C. gigas and indicates that this trait is heritable and selective breeding for thermal tolerance is a feasible and promising approach to reduce summer mortality.

Identification and Characterization of Infectious Pathogens Associated with Mass Mortalities of Pacific Oyster (Crassostrea gigas) Cultured in Northern China

The Pacific oyster (Crassostrea gigas) aquaculture industry increased rapidly in China with the introduction and promotion of triploid oysters in recent years. Mass mortalities affecting different life stages of Pacific oysters emerged periodically in several important production areas of Northern China. During 2020 and 2021, we conducted a passive two-year investigation of infectious pathogens linked to mass mortality. Ostreid herpesvirus-1 (OsHV-1) was detected to be associated with mass mortalities of hatchery larvae, but not juveniles and adults in the open sea. Protozoan parasites, such as Marteilia spp., Perkinsus spp. and Bonamia spp. were not detected. Bacterial isolation and identification revealed that Vibrio natriegens and Vibrio alginolyticus were the most frequently (9 out of 13) identified two dominant bacteria associated with mass mortalities. Pseudoalteromonas spp. was identified as the dominant bacteria in three mortality events that occurred during the cold season. Further bacteriological analysis was conducted on two representative isolates of V. natriegens and V. alginolyticus, designated as CgA1-1 and CgA1-2. Multisequence analysis (MLSA) showed that CgA1-1 and CgA1-2 were closely related to each other and nested within the Harveyi clade. Bacteriological investigation revealed faster growth, and more remarkable haemolytic activity and siderophore production capacity at 25 °C than at 15 °C for both CgA1-1 and CgA1-2. The accumulative mortalities of experimental immersion infections were also higher at 25 °C (90% and 63.33%) than at 15 °C (43.33% and 33.33%) using both CgA1-1 and CgA1-2, respectively. Similar clinical and pathological features were identified in samples collected during both naturally and experimentally occurring mortalities, such as thin visceral mass, discolouration, and connective tissue and digestive tube lesions. The results presented here highlight the potential risk of OsHV-1 to hatchery production of larvae, and the pathogenic role of V. natriegens and V. alginolyticus during mass mortalities of all life stages of Pacific oysters in Northern China.

Genotypic Characterization of Infectious Spleen and Kidney Necrosis Virus (ISKNV) in Southeast Asian Aquaculture

Infectious spleen and kidney necrosis virus (ISKNV) is a species within the genus Megalocytivirus (family Iridoviridae), which causes high mortality disease in many freshwater and marine fish species. ISKNV was first reported in Asia and is an emerging threat to aquaculture with increasing global distribution, in part due to its presence in ornamental fish with clinical and subclinical infections. The species ISKNV includes three genotypes: red seabream iridovirus (RSIV), turbot reddish body iridovirus (TRBIV), and ISKNV. There is an increasing overlap in the recognized range of susceptible fish hosts and the geographic distribution of these distinct genotypes. To better understand the disease caused by ISKNV, a nucleic acid hybridization capture enrichment was used prior to sequencing to characterize whole genomes from archived clinical specimens of aquaculture and ornamental fish from Southeast Asia (n = 16). The method was suitable for tissue samples containing 2.50 × 104–4.58 × 109 ISKNV genome copies mg−1. Genome sequences determined using the hybridization capture method were identical to those obtained directly from tissues when there was sufficient viral DNA to sequence without enrichment (n = 2). ISKNV genomes from diverse locations, environments, and hosts had very high similarity and matched established genotype classifications (14 ISKNV genotype Clade 1 genomes with >98.81% nucleotide similarity). Conversely, two different genotypes were obtained at the same time and location (RSIV and ISKNV from grouper, Indonesia with 92.44% nucleotide similarity). Gene-by-gene analysis with representative ISKNV genomes identified 59 core genes within the species (>95% amino acid identity). The 14 Clade 1 ISKNV genomes in this study had 100% aa identity for 92–105 of 122 predicted genes. Despite high overall sequence similarity, phylogenetic analyses using single nucleotide polymorphisms differentiated isolates from different host species, country of origin, and time of collection. Whole genome studies of ISKNV and other megalocytiviruses enable genomic epidemiology and will provide information to enhance disease control in aquaculture.

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.

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.

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.

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.

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.

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

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

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

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