Mass mortalities of unknown aetiology in Pacific oysters Crassostrea gigas in Port Stephens, New South Wales, Australia

Comprehensive Multi-omics Approaches Provide Insights to Summer Mortality in the Clam Meretrix petechialis

Bivalve mass mortalities have been reported worldwide, which not only can be explained as a result of pathogen infection, but may reflect changes in environments. Although these episodes were often reported, there was limited information concerning the molecular responses to various stressors leading to summer mortality. In the present work, RNA sequencing (RNA-seq), tandem mass tagging (TMT)-based quantitative proteomics, and 16S rRNA sequencing were used to explore the natural outbreak of summer mortality in the clam Meretrix petechialis. We identified a total of 172 differentially expressed genes (DEGs) and 222 differentially expressed proteins (DEPs) in the diseased group compared to the normal group. The inconsistent expression profiles of immune DEGs/DEPs may be due to the immune dysregulation of the diseased clams. Notably, 11 solute carrier family genes were found among the top 20 down-regulated genes in the diseased group, indicating that weakened transmembrane transport ability might occur in the diseased clams. Integration analysis of transcriptomic and proteomic results showed that many metabolic processes such as "arginine and proline metabolism" and "tyrosine metabolism" were inhibited in the diseased group, suggesting metabolic inhibition. Moreover, 16S rRNA sequencing revealed that the microbial composition of clam hepatopancreas was disordered in the diseased group. The comparison of DEGs expression between the natural summer mortality event and an artificial challenge experiment involving both Vibrio infection and heat stress revealed 9/15 genes showing similar expression trends between the two conditions, suggesting that the summer mortality might be caused by a combination of high temperature and Vibrio infection. These results would deepen our understanding of summer mortality and provide candidate resistance markers for clam resistance breeding.

Energetic budget of diploid and triploid eastern oysters during a summer die-off

Triploid oysters are widely used in off-bottom aquaculture of eastern oysters, Crassostrea virginica. However, farmers of the Gulf of Mexico (GoM) and Atlantic coast estuaries have observed unresolved, late-spring die-offs of triploid oysters, threatening the sustainability of triploid aquaculture. To investigate this, the physiological processes underlying oyster growth (e.g., feeding, respiration) and mortality of one-year-old diploid and triploid oysters were compared in early summer following an uptick in mortality. It was predicted that higher triploid mortality was the result of energetic imbalances (increased metabolic demands and decreased feeding behavior). Oyster clearance rates, percentage of time valves were open, absorption efficiency, oxygen consumption rates (basal and routine), ammonia excretion rate were measured in the laboratory and scope for growth was calculated. In addition, their condition index, gametogenic stage, Perkinsus marinus infection level, and mortality were measured. Mortality of triploids in the laboratory was greater than for diploids, mirroring mortality observed in a related field study. The physiological parameters measured, however, could not explain triploid mortality. Scope for growth, condition index, and clearance rates of triploids were greater than for diploids, suggesting sufficient energy reserves, while all other measurements where similar between the ploidies. It remains to be determined whether mortality could be caused from disruption of energy homeostasis at the cellular level.

Host-microbiota interactions play a crucial role in oyster adaptation to rising seawater temperature in summer

Climate change, represented by rising and fluctuating temperature, induces systematic changes in marine organisms and in their bacterial symbionts. However, the role of host-microbiota interactions in the host's response to rising temperature and the underlying mechanisms are incompletely understood in marine organisms. Here, the symbiotic intestinal microbiota and transcriptional responses between diploid and triploid oysters that displayed susceptible and resistant performance under the stress of rising temperature during a summer mortality event were compared to investigate the host-microbiota interactions. The rising and fluctuating temperatures triggered an earlier onset and higher mortality in susceptible oysters (46.7%) than in resistant oysters (17.3%). Correlation analysis between microbial properties and environmental factors showed temperature was strongly correlated with indices of α-diversity and the abundance of top 10 phyla, indicating that temperature significantly shaped the intestinal microbiota of oysters. The microbiota structure of resistant oysters exhibited more rapid changes in composition and diversity compared to susceptible oysters before peak mortality, indicating that resistant oysters possessed a stronger ability to regulate their symbiotic microbiota. Meanwhile, linear discriminant analysis effect size (LefSe) analysis found that the probiotics Verrucomicrobiales and Clostridiales were highly enriched in resistant oysters, and that potential pathogens Betaproteobacteriales and Acidobacteriales were enriched in susceptible oysters. These results implied that the symbiotic microbiota played a significant role in the oysters' adaptation to rising temperature. Accompanying the decrease in unfavorable bacteria before peak mortality, genes related to phagocytosis and lysozymes were upregulated and the xenobiotics elimination pathway was exclusively expressed in resistant oysters, demonstrating the validity of these immunological functions in controlling proliferation of pathogens driven by rising temperature. Compromised immunological functions might lead to proliferation of pathogens in susceptible oysters. This study might uncover a conserved mechanism of adaptation to rising temperature in marine invertebrates from the perspective of interactions between host and symbiotic microbiota.

Phylogenomic diversity of Vibrio species and other Gammaproteobacteria isolated from Pacific oysters (Crassostrea gigas) during a summer mortality outbreak

The Pacific oyster (PO), Crassostrea gigas, is an important commercial marine species but periodically experiences large stock losses due to disease events known as summer mortality. Summer mortality has been linked to environmental perturbations and numerous viral and bacterial agents, indicating this disease is multifactorial in nature. In 2013 and 2014, several summer mortality events occurred within the Port Stephens estuary (NSW, Australia). Extensive culture and molecular-based investigations were undertaken and several potentially pathogenic Vibrio species were identified. To improve species identification and genomically characterise isolates obtained from this outbreak, whole-genome sequencing (WGS) and subsequent genomic analyses were performed on 48 bacterial isolates, as well as a further nine isolates from other summer mortality studies using the same batch of juveniles. Average nucleotide identity (ANI) identified most isolates to the species level and included members of the Photobacterium, Pseudoalteromonas, Shewanella and Vibrio genera, with Vibrio species making up more than two-thirds of all species identified. Construction of a phylogenomic tree, ANI analysis, and pan-genome analysis of the 57 isolates represents the most comprehensive culture-based phylogenomic survey of Vibrios during a PO summer mortality event in Australian waters and revealed large genomic diversity in many of the identified species. Our analysis revealed limited and inconsistent associations between isolate species and their geographical origins, or host health status. Together with ANI and pan-genome results, these inconsistencies suggest that to determine the role that microbes may have in Pacific oyster summer mortality events, isolate identification must be at the taxonomic level of strain. Our WGS data (specifically, the accessory genomes) differentiated bacterial strains, and coupled with associated metadata, highlight the possibility of predicting a strain's environmental niche and level of pathogenicity.

Core Community Persistence Despite Dynamic Spatiotemporal Responses in the Associated Bacterial Communities of Farmed Pacific Oysters

A breakdown in host-bacteria relationships has been associated with the progression of a number of marine diseases and subsequent mortality events. For the Pacific oyster, Crassostrea gigas, summer mortality syndrome (SMS) is one of the biggest constraints to the growth of the sector and is set to expand into temperate systems as ocean temperatures rise. Currently, a lack of understanding of natural spatiotemporal dynamics of the host-bacteria relationship limits our ability to develop microbially based monitoring approaches. Here, we characterised the associated bacterial community of C. gigas, at two Irish oyster farms, unaffected by SMS, over the course of a year. We found C. gigas harboured spatiotemporally variable bacterial communities that were distinct from bacterioplankton in surrounding seawater. Whilst the majority of bacteria-oyster associations were transient and highly variable, we observed clear patterns of stability in the form of a small core consisting of six persistent amplicon sequence variants (ASVs). This core made up a disproportionately large contribution to sample abundance (34 ± 0.14%), despite representing only 0.034% of species richness across the study, and has been associated with healthy oysters in other systems. Overall, our study demonstrates the consistent features of oyster bacterial communities across spatial and temporal scales and provides an ecologically meaningful baseline to track environmental change.

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.

Regional and oyster microenvironmental scale heterogeneity in the Pacific oyster bacterial community

Different organs of a host represent distinct microenvironments resulting in the establishment of multiple discrete bacterial communities within a host. These discrete bacterial communities can also vary according to geographical location. For the Pacific oyster, Crassostrea gigas, the factors governing bacterial diversity and abundance of different oyster microenvironments are poorly understood. In this study, the factors shaping bacterial abundance, diversity and composition associated with the C. gigas mantle, gill, adductor muscle and digestive gland were characterised using 16S (V3-V4) rRNA amplicon sequencing across six discrete estuaries. Both location and tissue-type, with tissue-type being the stronger determinant, were factors driving bacterial community composition. Bacterial communities from wave-dominated estuaries had similar compositions and higher bacterial abundance despite being geographically distant from one another, possibly indicating that functional estuarine morphology characteristics are a factor shaping the oyster bacterial community. Despite the bacterial community heterogeneity, examinations of the core bacterial community identified Spirochaetaceae bacteria as conserved across all sites and samples. Whereas members of the Vulcaniibacterium, Spirochaetaceae and Margulisbacteria, and Polynucleobacter were regionally conserved members of the digestive gland, gill and mantle bacterial communities, respectively. This indicates that baseline bacterial community profiles for specific locations are necessary when investigating bacterial communities in oyster health.

A New High Throughput Sequencing Assay for Characterizing the Diversity of Natural Vibrio Communities and Its Application to a Pacific Oyster Mortality Event

The Vibrio genus is notable for including several pathogens of marine animals and humans, yet characterization of Vibrio diversity using routine 16S rRNA sequencing methods is often constrained by poor resolution beyond the genus level. Here, a new high throughput sequencing approach targeting the heat shock protein (hsp60) as a phylogenetic marker was developed to more precisely discriminate members of the Vibrio genus in environmental samples. The utility of this new assay was tested using mock communities constructed from known dilutions of Vibrio isolates. Relative to standard and Vibrio-specific 16S rRNA sequencing assays, the hsp60 assay delivered high levels of fidelity with the mock community composition at the species level, including discrimination of species within the Vibrio harveyi clade. This assay was subsequently applied to characterize Vibrio community composition in seawater and delivered substantially improved taxonomic resolution of Vibrio species compared to 16S rRNA analysis. Finally, this assay was applied to examine patterns in the Vibrio community within oysters during a Pacific oyster mortality event. In these oysters, the hsp60 assay identified species-level Vibrio community shifts prior to disease onset, pinpointing V. harveyi as a putative pathogen. Given that shifts in the Vibrio community can precede, cause, and follow disease onset in numerous marine organisms, there is a need for an accurate high throughput assay for defining Vibrio community composition in natural samples. This Vibrio-centric hsp60 sequencing assay offers the potential for precise high throughput characterization of Vibrio diversity, providing an enhanced platform for dissecting Vibrio dynamics in the environment.

Influence of environment on the pathogenesis of Ostreid herpesvirus-1 (OsHV-1) infections in Pacific oysters (Crassostrea gigas) through differential microbiome responses

The oyster microbiome is thought to contribute to the pathogenesis of mass mortality disease in Pacific oysters, associated with OsHV-1. As filter-feeders, oysters host a microbiota that can be influenced by the estuarine environment. This may alter susceptibility to OsHV-1 infections, causing variable mortality. This study aimed at: (1) differences in the microbiome of Pacific oysters with a common origin but grown in geographically distinct estuaries; (2) evaluating changes occurring in the microbiota, especially in Vibrio, and (3) differential responses of the oyster microbiome, in response to an OsHV-1 infection. Pacific oysters sourced from a single hatchery but raised separately in Patonga Creek, Shoalhaven River and Clyde River of NSW, Australia, were used and challenged with OsHV-1. The initial microbiome composition was different in the three batches and changed further, post-injection (p < 0.05). The Patonga oysters with the highest mortality also had higher OsHV-1 and Vibrio quantities compared to the other two batches (p < 0.05). The higher initial bacterial diversity in Patonga oysters decreased in moribund oysters which was not observed in the other two batches (p < 0.05). The microbiome of survivors of OsHV-1 infection and negative control oysters of two batches, did not show any changes with the relevant pre-challenged microbiome. A strong correlation was observed between the OsHV-1 and Vibrio quantities in OsHV-1 infected oysters (r = 0.6; p < 0.001). In conclusion, the Pacific oyster microbiome differed in different batches despite a common hatchery origin. Different microbiomes responded differently with a differential outcome of OsHV-1 challenge. The higher Vibrio load in oysters with higher OsHV-1 content and higher mortality, suggests a role in Vibrio in the pathogenesis of this mortality disease. This study provided insights of the potential of different estuarine environments to shape the Pacific oyster microbiome and how different microbiomes are associated with different outcomes of OsHV-1 infection.

Simulated Marine Heat Wave Alters Abundance and Structure of Vibrio Populations Associated with the Pacific Oyster Resulting in a Mass Mortality Event

Marine heat waves are predicted to become more frequent and intense due to anthropogenically induced climate change, which will impact global production of seafood. Links between rising seawater temperature and disease have been documented for many aquaculture species, including the Pacific oyster Crassostrea gigas. The oyster harbours a diverse microbial community that may act as a source of opportunistic pathogens during temperature stress. We rapidly raised the seawater temperature from 20 °C to 25 °C resulting in an oyster mortality rate of 77.4%. Under the same temperature conditions and with the addition of antibiotics, the mortality rate was only 4.3%, strongly indicating a role for bacteria in temperature-induced mortality. 16S rRNA amplicon sequencing revealed a change in the oyster microbiome when the temperature was increased to 25 °C, with a notable increase in the proportion of Vibrio sequences. This pattern was confirmed by qPCR, which revealed heat stress increased the abundance of Vibrio harveyi and Vibrio fortis by 324-fold and 10-fold, respectively. Our findings indicate that heat stress-induced mortality of C. gigas coincides with an increase in the abundance of putative bacterial pathogens in the oyster microbiome and highlights the negative consequences of marine heat waves on food production from aquaculture.

Variability in the Composition of Pacific Oyster Microbiomes Across Oyster Families Exhibiting Different Levels of Susceptibility to OsHV-1 μvar Disease

Oyster diseases are a major impediment to the profitability and growth of the oyster aquaculture industry. In recent years, geographically widespread outbreaks of disease caused by ostreid herpesvirus-1 microvariant (OsHV-1 μvar) have led to mass mortalities among Crassostrea gigas, the Pacific Oyster. Attempts to minimize the impact of this disease have been largely focused on breeding programs, and although these have shown some success in producing oyster families with reduced mortality, the mechanism(s) behind this protection is poorly understood. One possible factor is modification of the C. gigas microbiome. To explore how breeding for resistance to OsHV-1 μvar affects the oyster microbiome, we used 16S rRNA amplicon sequencing to characterize the bacterial communities associated with 35 C. gigas families, incorporating oysters with different levels of susceptibility to OsHV-1 μvar disease. The microbiomes of disease-susceptible families were significantly different to the microbiomes of disease-resistant families. OTUs assigned to the Photobacterium, Vibrio, Aliivibrio, Streptococcus, and Roseovarius genera were associated with low disease resistance. In partial support of this finding, qPCR identified a statistically significant increase of Vibrio-specific 16S rRNA gene copies in the low disease resistance families, possibly indicative of a reduced host immune response to these pathogens. In addition to these results, examination of the core microbiome revealed that each family possessed a small core community, with OTUs assigned to the Winogradskyella genus and the Bradyrhizobiaceae family consistent members across most disease-resistant families. This study examines patterns in the microbiome of oyster families exhibiting differing levels of OsHV-1 μvar disease resistance and reveals some key bacterial taxa that may provide a protective or detrimental role in OsHV-1 μvar disease outbreaks.

Characterisation of the Pacific Oyster Microbiome During a Summer Mortality Event

The Pacific oyster, Crassostrea gigas, is a key commercial species that is cultivated globally. In recent years, disease outbreaks have heavily impacted C. gigas stocks worldwide, with many losses incurred during summer. A number of infectious agents have been associated with these summer mortality events, including viruses (particularly Ostreid herpesvirus 1, OsHV-1) and bacteria; however, cases where no known aetiological agent can be identified are common. In this study, we examined the microbiome of disease-affected and disease-unaffected C. gigas during a 2013-2014 summer mortality event in Port Stephens (Australia) where known oyster pathogens including OsHV-1 were not detected. The adductor muscle microbiomes of 70 C. gigas samples across 12 study sites in the Port Stephens estuary were characterised using 16S rRNA (V1-V3 region) amplicon sequencing, with the aim of comparing the influence of spatial location and disease state on the oyster microbiome. Spatial location was found to be a significant determinant of the disease-affected oyster microbiome. Furthermore, microbiome comparisons between disease states identified a significant increase in rare operational taxonomic units (OTUs) belonging to Vibrio harveyi and an unidentified member of the Vibrio genus in the disease-affected microbiome. This is indicative of a potential role of Vibrio species in oyster disease and supportive of previous culture-based examination of this mortality event.

Transgenerational plasticity and antiviral immunity in the Pacific oyster (Crassostrea gigas) against Ostreid herpesvirus 1 (OsHV-1)

The oyster's immune system is capable of adapting upon exposure to a pathogen-associated molecular pattern (PAMP) to have an enhanced secondary response against the same type of pathogen. This has been demonstrated using poly(I:C) to elicit an antiviral response in the Pacific oyster (Crassostrea gigas) against Ostreid herpesvirus (OsHV-1). Improved survival following exposure to poly(I:C) has been found in later life stages (within-generational immune priming) and in the next generation (transgenerational immune priming). The mechanism that the oyster uses to transfer immunity to the next generation is unknown. Here we show that oyster larvae have higher survival to OsHV-1 when their mothers, but not their fathers, are exposed to poly(I:C) prior to spawning. RNA-seq provided no evidence to suggest that parental exposure to poly(I:C) reconfigures antiviral gene expression in unchallenged larvae. We conclude that the improved survival of larvae might occur via maternal provisioning of antiviral compounds in the eggs.

Oyster disease in a changing environment: Decrypting the link between pathogen, microbiome and environment

Shifting environmental conditions are known to be important triggers of oyster diseases. The mechanism(s) behind these synergistic effects (interplay between host, environment and pathogen/s) are often not clear, although there is evidence that shifts in environmental conditions can affect oyster immunity, and pathogen growth and virulence. However, the impact of shifting environmental parameters on the oyster microbiome and how this affects oyster health and susceptibility to infectious pathogens remains understudied. In this review, we summarise the major diseases afflicting oysters with a focus on the role of environmental factors that can catalyse or amplify disease outbreaks. We also consider the potential role of the oyster microbiome in buffering or augmenting oyster disease outbreaks and suggest that a deeper understanding of the oyster microbiome, its links to the environment and its effect on oyster health and disease susceptibility, is required to develop new frameworks for the prevention and management of oyster diseases.