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

Evaluation of long-read sequencing for Ostreid herpesvirus type 1 genome characterization from Magallana gigas infected tissues

Since the 1990s, the Pacific oyster Magallana gigas has faced significant mortality, which has been associated with the detection of the Ostreid Herpesvirus type 1 (OsHV-1). Due to the complex genomic architecture and the presence of multiple genomic isomers, short-read sequencing using Illumina method struggles to accurately assemble tandem and repeat regions and to identify and characterize large structural variations in the OsHV-1 genome. Third-generation sequencing technologies, as long-read real-time nanopore sequencing from Oxford Nanopore Technologies (ONT), offer new possibilities for OsHV-1 whole-genome analysis. Identification of the best method for extraction of high molecular weight (HMW) DNA and development of accurate bioinformatic pipelines for its characterization are now required. To this end, we evaluated and compared six HMW methods and one conventional DNA extraction kit for their ability to extract OsHV-1 DNA from M. gigas-infected tissues. We then evaluated the ability of ONT sequencing to produce an accurate OsHV-1 genome from both whole-genome and "adaptive sampling" (AS) sequencing approaches. Finally, we evaluated the efficiency of bioinformatics tools for de novo assembly and consensus calling to generate accurate OsHV-1 genomes. The HMW DNA extraction kit coupled with ONT sequencing and dedicated bioinformatics tools allowed us to produce accurate OsHV-1 genomes compared to those assembled using Illumina technology. The AS approach allowed up to 60% enrichment for viral data, and the long reads generated by ONT allowed the characterization of OsHV-1 isomers. Together with its portability, this sequencing shows great promise as a diagnostic tool for the characterization of unculturable aquatic viruses directly from host tissues.IMPORTANCEMany aquatic viruses threaten commercially valuable species and cause significant economic losses during outbreaks. To improve our understanding of the origin, transmission patterns and spread of these viruses, additional genomic data are essential. However, genomic characterization of unculturable large DNA viruses is a major challenge. In the present study, we have successfully evaluated the ability of ONT sequencing and adaptive sequencing (AS) to sequence and assemble the complete OsHV-1 genome. Our results show that it is now possible to sequence the whole genome of large DNA viruses directly from infected host tissue, without the need for prior in vitro propagation or prior laboratory steps for virus enrichment.

Laboratory Replication of Ostreid Herpes Virus (OsHV-1) Using Pacific Oyster Tissue Explants

Pacific oysters (Crassostrea or Magallana gigas) are one of the most economically important aquaculture species globally. Over the past two decades, ostreid herpesvirus (OsHV-1) has become a major pathogen of cultured Pacific oysters, resulting in widespread mortality with a global distribution. Experimental use of OsHV-1 is challenging for many reasons, including both complexity of host-pathogen dynamics and a lack of functioning model systems. The goal of this study was to improve the tools available for working with OsHV-1 in both whole animals and in tissue explants established from oysters maintained in controlled laboratory conditions. Tissue explants were taken from oysters originating from two different sources that have different levels of mortality in experimental OsHV-1 infections and were exposed to OsHV-1. A whole-animal infection experiment was run concurrently as a comparison. Quantitative PCR and electron microscopy were used to confirm that the explants were capable of replicating OsHV-1. Furthermore, the quantitative PCR results suggest that the source of the oysters was significant in determining the outcome of infection in the explants, supporting the validity of the explant model for OsHV-1 infection. This tissue explant approach for studying OsHV-1 allows for the control of confounding factors in the disease outcome that is not possible in whole-animal experiments, providing a new tool for the study of OsHV-1 in Pacific oysters.

Comparative genomics of Ascetosporea gives new insight into the evolutionary basis for animal parasitism in Rhizaria

BACKGROUND

Ascetosporea (Endomyxa, Rhizaria) is a group of unicellular parasites infecting aquatic invertebrates. They are increasingly being recognized as widespread and important in marine environments, causing large annual losses in invertebrate aquaculture. Despite their importance, little molecular data of Ascetosporea exist, with only two genome assemblies published to date. Accordingly, the evolutionary origin of these parasites is unclear, including their phylogenetic position and the genomic adaptations that accompanied the transition from a free-living lifestyle to parasitism. Here, we sequenced and assembled three new ascetosporean genomes, as well as the genome of a closely related amphizoic species, to investigate the phylogeny, origin, and genomic adaptations to parasitism in Ascetosporea.

RESULTS

Using a phylogenomic approach, we confirm the monophyly of Ascetosporea and show that Paramyxida group with Mikrocytida, with Haplosporida being sister to both groups. We report that the genomes of these parasites are relatively small (12-36 Mb) and gene-sparse (~ 2300-5200 genes), while containing surprisingly high amounts of non-coding sequence (~ 70-90% of the genomes). Performing gene-tree aware ancestral reconstruction of gene families, we demonstrate extensive gene losses at the origin of parasitism in Ascetosporea, primarily of metabolic functions, and little gene gain except on terminal branches. Finally, we highlight some functional gene classes that have undergone expansions during evolution of the group.

CONCLUSIONS

We present important new genomic information from a lineage of enigmatic but important parasites of invertebrates and illuminate some of the genomic innovations accompanying the evolutionary transition to parasitism in this lineage. Our results and data provide a genetic basis for the development of control measures against these parasites.

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.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

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.

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.

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.

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.