Cooperation and cheating orchestrate Vibrio assemblages and polymicrobial synergy in oysters infected with OsHV-1 virus.

Polymicrobial infections threaten the health of humans and animals but remain understudied in natural systems. We recently described the Pacific Oyster Mortality Syndrome (POMS), a polymicrobial disease affecting oyster production worldwide. In the French Atlantic coast, the disease involves coinfection with ostreid herpesvirus 1 (OsHV-1) and virulent Vibrio. However, it is unknown whether consistent Vibrio populations are associated with POMS in different regions, how Vibrio contribute to POMS, and how they interact with OsHV-1 during pathogenesis. By connecting field-based approaches in a Mediterranean ecosystem, laboratory infection assays and functional genomics, we uncovered a web of interdependencies that shape the structure and function of the POMS pathobiota. We show that Vibrio harveyi and Vibrio rotiferianus are predominant in OsHV-1-diseased oysters and that OsHV-1 drives the partition of the Vibrio community observed in the field. However only V. harveyi synergizes with OsHV-1 by promoting mutual growth and accelerating oyster death. V. harveyi shows high-virulence potential and dampens oyster cellular defenses through a type 3 secretion system, making oysters a more favorable niche for microbe colonization. In addition, V. harveyi produces a key siderophore called vibrioferrin. This important resource promotes the growth of V. rotiferianus, which cooccurs with V. harveyi in diseased oysters, and behaves as a cheater by benefiting from V. harveyi metabolite sharing. Our data show that cooperative behaviors contribute to synergy between bacterial and viral coinfecting partners. Additional cheating behaviors further shape the polymicrobial consortium. Controlling cooperative behaviors or countering their effects opens avenues for mitigating polymicrobial diseases.

Emergence and clonal expansion of Vibrio aestuarianus lineages pathogenic for oysters in Europe.

Crassostrea gigas oysters represent a significant global food source with 4.7 million tons harvested per year. In 2001, the bacterium Vibrio aestuarianus subsp. francensis emerged as a pathogen that causes adult oyster mortality in France and Ireland. Its impact on oyster aquaculture has increased in Europe since its re-emergence in 2012. To better understand the evolutionary mechanisms leading to the emergence and persistence over time of this pathogen, we conducted a survey of mollusc diseases through national reference laboratories across Europe. We analysed 54 new genomes of Vibrio aestuarianus (Va) isolated from multiple environmental compartments since 2001, in areas with and without bivalve mortalities. We used a combination of comparative genomics and population genetics approaches and show that Va has a classical epidemic population structure from which the pathogenic Va francensis subspecies emerged and clonally expanded. Furthermore, we identified a specific cus-cop-containing island conferring copper resistance to Va francensis whose acquisition may have favoured the emergence of pathogenic lineages adapted and specialized to oysters.

First characterization of the parasite Haplosporidium costale in France and development of a real-time PCR assay for its rapid detection in the Pacific oyster, Crassostrea gigas.

The Pacific cupped oyster Crassostrea gigas is one of the most 'globalized' marine invertebrates and its production is predominant in many parts of the world including Europe. However, it is threatened by mortality events associated with pathogenic microorganisms such as the virus OsHV-1 and the bacteria Vibrio aestuarianus. C. gigas is also a host for protozoan parasites including haplosporidians. In contrast with Haplosporidium nelsoni previously detected in Europe, H. costale was considered exotic although its presence in French oysters was suggested in the 1980s based on ultrastructural examination. Here, a combination of light and transmission electron microscopy, PCR and sequencing allowed characterizing the presence of the parasite in the context of low mortality events which occurred in 2019 in France. Histological observation revealed the presence of uninucleated, plasmodial and spore stages within the connective tissues of some oysters. Ultrastructural features were similar to H. costale ones in particular the presence of axe-shaped haplosporosomes in spore cytoplasms. Three fragments of the genome including partial small subunit rRNA gene, the ITS-1, 5.8S and ITS-2 array and part of the actin gene were successfully sequenced and grouped with H. costale homologous sequences. This is the first time that the presence of H. costale was confirmed in C. gigas in France. Furthermore, a TaqMan real-time PCR assay was developed and validated [DSe = 92.6% (78.2-99.8) and DSp = 95.5% (92.3-98.6)] to enable the rapid and specific detection of the parasite. The application of the PCR assay on archived samples revealed that the parasite has been present in French oyster populations at least since 2008. Considering the little information available on this parasite, the newly developed TaqMan assay will be very helpful to investigate the temporal and geographic distribution and the life cycle of the parasite in France and more generally in C. gigas geographic range.

Vibrio aestuarianus subsp. cardii subsp. nov., pathogenic to the edible cockles Cerastoderma edule in France, and establishment of Vibrio aestuarianus subsp. aestuarianus subsp. nov. and Vibrio aestuarianus subsp. francensis subsp. nov.

Cockle mortality events have been reported in northern France since 2012. In the present study, we describe and investigate the implication of a potential bacterial causative agent in cockle mortality. Bacteria isolated from five different cockle mortality events were characterized and studied. Using phenotypic analysis combined with DNA-DNA hybridization (DDH) and whole genome sequencing, the isolates were shown to belong to Vibrio aestuarianus, a species regularly detected in France during oyster mortality events. Comparison of the strains from cockles with strains from French oysters and the type strain showed that the strains from cockles were genetically different to those from oysters and also different to the V. aestuarianus type strain. Moreover, the cockle and oyster strains were classified into two different, but close, groups both separated from the type strain by: (1) analyses of the ldh gene sequences; (2) DDH assays between 12/122 3T3T (LMG 31436T=DSM 109723T), a representative cockle strain, 02/041T (CIP 109791T=LMG 24517T) representative oyster strain and V. aestuarianus type strain LMG 7909T; (3) average nucleotide identity values calculated on the genomes; and (4) phenotypic traits. Finally, results of MALDI-TOF analyses also revealed specific peaks discriminating the three representative strains. The toxicity of representative strains of these cockle isolates was demonstrated by experimental infection of hatchery-produced cockles. The data therefore allow us to propose two novel subspecies of Vibrio aestuarianus: Vibrio aestuarianus subsp. cardii subsp. nov. for the cockle strains and Vibrio aestuarianus subsp. francensis subsp. nov. for the Pacific oyster strains, in addition to an emended description of the species Vibrio aestuarianus.

Modeling the Transmission of Vibrio aestuarianus in Pacific Oysters Using Experimental Infection Data.

Vibrio aestuarianus is a bacterium related to mortality outbreaks in Pacific oysters, Crassostrea gigas, in France, Ireland, and Scotland since 2011. Knowledge about its transmission dynamics is still lacking, impairing guidance to prevent and control the related disease spread. Mathematical modeling is a relevant approach to better understand the determinants of a disease and predict its dynamics in imperfectly observed pathosystems. We developed here the first marine epidemiological model to estimate the key parameters of V. aestuarianus infection at a local scale in a small and closed oyster population under controlled laboratory conditions. Using a compartmental model accounting for free-living bacteria in seawater, we predicted the infection dynamics using dedicated and model-driven collected laboratory experimental transmission data. We estimated parameters and showed that waterborne transmission of V. aestuarianus is possible under experimental conditions, with a basic reproduction number R0 of 2.88 (95% CI: 1.86; 3.35), and a generation time of 5.5 days. Our results highlighted a bacterial dose-dependent transmission of vibriosis at local scale. Global sensitivity analyses indicated that the bacteria shedding rate, the concentration of bacteria in seawater that yields a 50% chance of catching the infection, and the initial bacterial exposure dose W0 were three critical parameters explaining most of the variation in the selected model outputs related to disease spread, i.e., R0, the maximum prevalence, oyster survival curve, and bacteria concentration in seawater. Prevention and control should target the exposure of oysters to bacterial concentration in seawater. This combined laboratory-modeling approach enabled us to maximize the use of information obtained through experiments. The identified key epidemiological parameters should be better refined by further dedicated laboratory experiments. These results revealed the importance of multidisciplinary approaches to gain consistent insights into the marine epidemiology of oyster diseases.

Ecologically realistic model of infection for exploring the host damage caused by Vibrio aestuarianus.

Although vibrios are frequently associated with marine organisms mortality outbreaks, knowledge on their ecology and pathogenicity is sparse, thus limiting disease management and prophylactic strategies. Here, we investigated V. aestuarianus infection onset and progression in the wild, taking advantage of a 'claire' pond: a semi-closed system with limited seawater renewal, theoretically more adapted to disease transmission. We showed a positive association of the bacteria with oysters, which can constitute a reservoir for the bacteria in the winter. Moreover, passage through oysters was found to be necessary for experimental disease reproduction as vibrios shedding from diseased oysters have higher infectivity than from in vitro grown. We next developed an experimental 'ecologically realistic' infection model in a mesocosm, allowing infection by natural route. By means of this non-invasive protocol, we analysed the pathogenesis of the bacteria and demonstrated the importance of haemolymph for initial colonization and the septicaemic nature of this disease.

Development of a duplex Taqman real-time PCR assay for rapid identification of Vibrio splendidus-related and V. aestuarianus strains from bacterial cultures.

To enable the rapid and accurate identification of Vibrio splendidus-related and V. aestuarianus strains associated with Pacific cupped oyster Crassostrea gigas mortality, we developed a duplex Taqman real-time PCR assay and evaluated its efficacy. This technique proved to be rapid, sensitive, and specific and will be particularly valuable for epidemiologic studies.

Several strains, one disease: experimental investigation of Vibrio aestuarianus infection parameters in the Pacific oyster, Crassostrea gigas.

This study investigated oyster infection dynamics by different strains of Vibrio aestuarianus isolated before and after the apparent re-emergence of this pathogen observed in France in 2011. We conducted experiments to compare minimal infective dose, lethal dose 50 and bacterial shedding for six V. aestuarianus strains. Whatever the strain used, mortality was induced in juvenile oysters by intramuscular injection and reached 90-100% of mortality within 5 days. Moreover, bacterial shedding was comparable among strains and reached its maximum after 20 h (≈10 EXP5 bacteria/mL/animal). Similarly, our first estimations of lethal dose 50 were comparable among strains (minimal infective dose around 0.4 × 10EXP5 bacteria/mL and LD50 around 10EXP5 bacteria/mL) by using seawater containing freshly shed bacteria. These results indicate that, at least with these criteria, despite V. aestuarianus strains genetic diversity, the disease process is similar. The strains isolated after the apparent re-emergence of the bacteria in 2011, do not present a more acute virulence phenotype than the reference strains isolated between 2002 and 2007. Finally, our study provides original and noteworthy data indicating that infected oysters shed bacteria at a level above the threshold of LD50 a few days before they die, meaning that infection is expected to spread in a susceptible population.

Detection and distribution of ostreid herpesvirus 1 in experimentally infected Pacific oyster spat.

High mortality rates are reported in spat and larvae of Pacific oyster Crassostrea gigas and associated with ostreid herpesvirus 1 (OsHV-1) detection in France. Although the viral infection has been experimentally reproduced in oyster larvae and spat, little knowledge is currently available concerning the viral entry and its distribution in organs and tissues. This study compares OsHV-1 DNA and RNA detection and localization in experimentally infected oysters using two virus doses: a low dose that did not induce any mortality and a high dose inducing high mortality. Real time PCR demonstrated significant differences in terms of viral DNA amounts between the two virus doses. RNA transcripts were detected in oysters receiving the highest dose of viral suspension whereas no transcript was observed in oysters injected with the low dose. This study also allowed observing kinetics of viral DNA and RNA detection in different tissues of oyster spat. Finally, viral detection was significantly different in function of tissues (p<0.005), time (p<0.005) with an interaction between tissues and time (p<0.005) for each probe.

Can selection for resistance to OsHV-1 infection modify susceptibility to Vibrio aestuarianus infection in Crassostrea gigas? First insights from experimental challenges using primary and successive exposures.

Since 2008, the emergent virus OsHV-1µvar has provoked massive mortality events in Crassostrea gigas spat and juveniles in France. Since 2012, mortality driven by the pathogenic bacteria Vibrio aestuarianus has stricken market-sized adults. A hypothesis to explain the sudden increase in mortality observed in France since 2012 is that selective pressure due to recurrent viral infections could have led to a higher susceptibility of adults to Vibrio infection. In our study, two OsHV-1-resistant lines (AS and BS) and their respective controls (AC and BC) were experimentally challenged in the laboratory to determine their level of susceptibility to V. aestuarianus infection. At the juvenile stage, the selected lines exhibited lower mortality (14 and 33%) than the control lines (71 and 80%), suggesting dual-resistance to OsHV-1 and V. aestuarianus in C. gigas. Interestingly, this pattern was not observed at the adult stage, where higher mortality was detected for AS (68%) and BC (62%) than AC (39%) and BS (49%). These results were confirmed by the analysis of the expression of 31 immune-related genes in unchallenged oysters. Differential gene expression discriminated oysters according to their susceptibility to infection at both the juvenile and adult stages, suggesting that resistance to V. aestuarianus infection resulted in complex interactions between the genotype, stage of development and immunity status. Finally, survivors of the V. aestuarianus challenge at the juvenile stage still exhibited significant mortality at the adult stage during a second and third V. aestuarianus challenge, indicating that these survivors were not genetically resistant.

Whole-genome amplification: a useful approach to characterize new genes in unculturable protozoan parasites such as Bonamia exitiosa.

Bonamia exitiosa is an intracellular parasite (Haplosporidia) that has been associated with mass mortalities in oyster populations in the Southern hemisphere. This parasite was recently detected in the Northern hemisphere including Europe. Some representatives of the Bonamia genus have not been well categorized yet due to the lack of genomic information. In the present work, we have applied Whole-Genome Amplification (WGA) technique in order to characterize the actin gene in the unculturable protozoan B. exitiosa. This is the first protein coding gene described in this species. Molecular analysis revealed that B. exitiosa actin is more similar to Bonamia ostreae actin gene-1. Actin phylogeny placed the Bonamia sp. infected oysters in the same clade where the herein described B. exitiosa actin resolved, offering novel information about the classification of the genus. Our results showed that WGA methodology is a promising and valuable technique to be applied to unculturable protozoans whose genomic material is limited.

Autophagy plays an important role in protecting Pacific oysters from OsHV-1 and Vibrio aestuarianus infections.

Recent mass mortality outbreaks around the world in Pacific oysters, Crassostrea gigas, have seriously affected the aquaculture economy. Although the causes for these mortality outbreaks appear complex, infectious agents are involved. Two pathogens are associated with mass mortality outbreaks, the virus ostreid herpesvirus 1 (OsHV-1) and the bacterium Vibrio aestuarianus. Here we describe the interactions between these 2 pathogens and autophagy, a conserved intracellular pathway playing a key role in innate immunity. We show for the first time that autophagy pathway is present and functional in Pacific oysters and plays an important role to protect animals from infections. This study contributes to better understand the innate immune system of Pacific oysters.

A single regulatory gene is sufficient to alter Vibrio aestuarianus pathogenicity in oysters.

Oyster diseases caused by pathogenic vibrios pose a major challenge to the sustainability of oyster farming. In France, since 2012 a disease affecting specifically adult oysters has been associated with the presence of Vibrio aestuarianus. Here, by combining genome comparison, phylogenetic analyses and high-throughput infections of strains isolated before or during the recent outbreaks, we show that virulent strains cluster into two V. aestuarianus lineages independently of the sampling dates. The bacterial lethal dose was not different between strains isolated before or after 2012. Hence, the emergence of a new highly virulent clonal strain is unlikely. Each lineage comprises nearly identical strains, the majority of them being virulent, suggesting that within these phylogenetically coherent virulent lineages a few strains have lost their pathogenicity. Comparative genomics allowed the identification of a single frameshift in a non-virulent strain. This mutation affects the varS gene that codes for a signal transduction histidine-protein kinase. Genetic analyses confirmed that varS is necessary for infection of oysters and for a secreted metalloprotease expression. For the first time in a Vibrio species, we show here that VarS is a key factor of pathogenicity.

Ostreid herpesvirus type 1 replication and host response in adult Pacific oysters, Crassostrea gigas.

Since 2008, massive mortality outbreaks associated with OsHV-1 detection have been reported in Crassostrea gigas spat and juveniles in several countries. Nevertheless, adult oysters do not demonstrate mortality in the field related to OsHV-1 detection and were thus assumed to be more resistant to viral infection. Determining how virus and adult oyster interact is a major goal in understanding why mortality events are not reported among adult Pacific oysters. Dual transcriptomics of virus-host interactions were explored by real-time PCR in adult oysters after a virus injection. Thirty-nine viral genes and five host genes including MyD88, IFI44, IkB2, IAP and Gly were measured at 0.5, 10, 26, 72 and 144 hours post infection (hpi). No viral RNA among the 39 genes was detected at 144 hpi suggesting the adult oysters are able to inhibit viral replication. Moreover, the IAP gene (oyster gene) shows significant up-regulation in infected adults compared to control adults. This result suggests that over-expression of IAP could be a reaction to OsHV-1 infection, which may induce the apoptotic process. Apoptosis could be a main mechanism involved in disease resistance in adults. Antiviral activity of haemolymph against herpes simplex virus (HSV-1) was not significantly different between infected adults versus control.

Dual transcriptomics of virus-host interactions: comparing two Pacific oyster families presenting contrasted susceptibility to ostreid herpesvirus 1.

BACKGROUND: Massive mortality outbreaks affecting Pacific oyster (Crassostrea gigas) spat in various countries have been associated with the detection of a herpesvirus called ostreid herpesvirus type 1 (OsHV-1). However, few studies have been performed to understand and follow viral gene expression, as it has been done in vertebrate herpesviruses. In this work, experimental infection trials of C. gigas spat with OsHV-1 were conducted in order to test the susceptibility of several bi-parental oyster families to this virus and to analyze host-pathogen interactions using in vivo transcriptomic approaches. RESULTS: The divergent response of these oyster families in terms of mortality confirmed that susceptibility to OsHV-1 infection has a significant genetic component. Two families with contrasted survival rates were selected. A total of 39 viral genes and five host genes were monitored by real-time PCR. Initial results provided information on (i) the virus cycle of OsHV-1 based on the kinetics of viral DNA replication and transcription and (ii) host defense mechanisms against the virus. CONCLUSIONS: In the two selected families, the detected amounts of viral DNA and RNA were significantly different. This result suggests that Pacific oysters are genetically diverse in terms of their susceptibility to OsHV-1 infection. This contrasted susceptibility was associated with dissimilar host gene expression profiles. Moreover, the present study showed a positive correlation between viral DNA amounts and the level of expression of selected oyster genes.

First description of French V. tubiashii strains pathogenic to mollusk: I. Characterization of isolates and detection during mortality events.

Nine dominant bacterial isolates were obtained from different batches of Crassostrea gigas spat experiencing high mortality rates in a French experimental hatchery/nursery in 2007. Using phenotypic analysis combined with multilocus sequence analysis, the isolates were shown to be genetically close to the Vibrio tubiashii type strain. Based on (1) analyses of the recA gene sequences; (2) the results of DNA-DNA hybridization assays between 07/118 T2 (LMG 27884=CECT 8426), which is a representative strain, and the V. tubiashii type strain (69%); and (3) phenotypic traits, the bacteria were classified in a group close to American V. tubiashii strain. Its virulence (70% of mortalities) and the toxicity of the extracellular products of 07/118 T2 was demonstrated (41% of mortalities). Moreover, a QPCR diagnostic tool targeting the gyrB gene was developed to investigate the epidemiological significance of V. tubiashii in French oyster mortality outbreaks recorded by the national surveillance network. Of the 21 batches originating from hatcheries, only two were positive, whereas V. tubiashii DNA could not be detected in any of the batches of moribund animals collected in field/outdoor facilities. These results demonstrate the existence of a group of virulent V. tubiashii in France that episodically infect C. gigas.