Determination of risk factors for herpesvirus outbreak in oysters using a broad-scale spatial epidemiology framework

Diseases of marine fish and shellfish in an age of rapid climate change

A recurring trend in evidence scrutinized over the past few decades is that disease outbreaks will become more frequent, intense, and widespread on land and in water, due to climate change. Pathogens and the diseases they inflict represent a major constraint on seafood production and yield, and by extension, food security. The risk(s) for fish and shellfish from disease is a function of pathogen characteristics, biological species identity, and the ambient environmental conditions. A changing climate can adversely influence the host and environment, while augmenting pathogen characteristics simultaneously, thereby favoring disease outbreaks. Herein, we use a series of case studies covering some of the world's most cultured aquatic species (e.g., salmonids, penaeid shrimp, and oysters), and the pathogens (viral, fungal, bacterial, and parasitic) that afflict them, to illustrate the magnitude of disease-related problems linked to climate change.

Lipidomic insights into the immune response and pearl formation in transplanted pearl oyster Pinctada fucata martensii

During pearl culture, the excess immune responses may induce nucleus rejection and death of pearl oysters after transplantation. To better understand the immune response and pearl formation, lipidomic analysis was applied to investigate changes in the serum lipid profile of pearl oyster Pinctada fucata martensii following transplantation. In total, 296 lipid species were identified by absolute quantitation. During wound healing, the content of TG and DG initially increased and then decreased after 3 days of transplantation with no significant differences, while the level of C22:6 decreased significantly on days 1 and 3. In the early stages of transplantation, sphingosine was upregulated, whereas PC and PUFAs were downregulated in transplanted pearl oyster. PI was upregulated during pearl sac development stages. GP and LC-PUFA levels were upregulated during pearl formation stage. In order to identify enriched metabolic pathways, pathway enrichment analysis was conducted. Five metabolic pathways were found significantly enriched, namely glycosylphosphatidylinositol-anchor biosynthesis, glycerophospholipid metabolism, alpha-linolenic acid metabolism, linoleic acid metabolism and arachidonic acid metabolism. Herein, results suggested that the lipids involved in immune response, pearl sac maturation, and pearl formation in the host pearl oyster after transplantation, which might lead to an improvement in the survival rate and pearl quality of transplanted pearl oyster.

Harsh intertidal environment enhances metabolism and immunity in oyster (Crassostrea gigas) spat

The Pacific oyster Crassostrea gigas is established in the marine intertidal zone, experiencing rapid and highly dynamic environmental changes throughout the tidal cycle. Depending on the bathymetry, oysters face oxygen deprivation, lack of nutrients, and high changes in temperature during alternation of the cycles of emersion/immersion. Here we showed that intertidal oysters at a bathymetry level of 3 and 5 m delayed by ten days the onset of mortality associated with Pacific Oyster Mortality Syndrome (POMS) as compared to subtidal oysters. Intertidal oysters presented a lower growth but similar energetic reserves to subtidal oysters but induced proteomic changes indicative of a boost in metabolism, inflammation, and innate immunity that may have improved their resistance during infection with the Ostreid herpes virus. Our work highlights that intertidal harsh environmental conditions modify host-pathogen interaction and improve oyster health. This study opens new perspectives on oyster farming for mitigation strategies based on tidal height.

Spatial effects in parasite-induced marine diseases of immobile hosts

Emerging marine infectious diseases pose a substantial threat to marine ecosystems and the conservation of their biodiversity. Compartmental models of epidemic transmission in marine sessile organisms, available only recently, are based on non-spatial descriptions in which space is homogenized and parasite mobility is not explicitly accounted for. However, in realistic scenarios epidemic transmission is conditioned by the spatial distribution of hosts and the parasites' mobility patterns, calling for an explicit description of space. In this work, we develop a spatially explicit individual-based model to study disease transmission by waterborne parasites in sessile marine populations. We investigate the impact of spatial disease transmission through extensive numerical simulations and theoretical analysis. Specifically, the effects of parasite mobility into the epidemic threshold and the temporal progression of the epidemic are assessed. We show that larger values of pathogen mobility imply more severe epidemics, as the number of infections increases, and shorter timescales to extinction. An analytical expression for the basic reproduction number of the spatial model, R ~ 0 , is derived as a function of the non-spatial counterpart, R ₀, which characterizes a transition between a disease-free and a propagation phase, in which the disease propagates over a large fraction of the system.

A 26-year time series of mortality and growth of the Pacific oyster C. gigas recorded along French coasts

We used a compiled data set from a monitoring network of oyster production coordinated by IFREMER (the French Research Institute for the Exploitation of the Sea). This network monitors the growth and mortality of the Pacific oyster Crassostrea gigas along French coasts since 1993. The archive, although publicly available, has been challenging to use due to changes in protocols and little information on metadata. Here, we describe data collection for almost 30 years, cleaning and processing. For 13 locations, we modeled growth and mortality of spat (less than one-year-old individuals) and half-grown oysters (between one and two-year-old individuals) as a function of time to cope with changes in data acquisition frequency, and produced standardized annual growth and cumulative mortality indicators to improve data usability. This improved database is expected to be used by ecologists interested in the evolution of life-cycle indicators of a marine species under the influence of climate change. It can also be valuable for epidemiologists because mortality data traces the emergence and spread of a massive epizootic.

Measurement(s)	body mass • Cumulative mortality
Technology Type(s)	balance • Count
Factor Type(s)	Location • Sampling year
Sample Characteristic - Organism	Crassostrea gigas
Sample Characteristic - Environment	sea coast
Sample Characteristic - Location	Metropolitan France

The first detection of a novel OsHV-1 microvariant in San Diego, California, USA

The spread, emergence, and adaptation of pathogens causing marine disease has been problematic to fisheries and aquaculture industries for the last several decades creating the need for strategic management and biosecurity practices. The Pacific oyster (Crassostrea gigas), a highly productive species globally, has been a target of disease and mortality caused by a viral pathogen, the Ostreid herpesvirus 1 (OsHV-1) and its microvariants (OsHV-1 µvars). During routine surveillance to establish health history at a shellfish aquaculture nursery system in San Diego, California, the presence of OsHV-1 in Pacific oyster juveniles was detected. Quantification of OsHV-1 in tissues of oysters revealed OsHV-1 viral loads > 10⁶ copies/mg. We characterized and identified the OsHV-1 variant by sequencing of ORFs 4 (C2/C6) and 43 (IA1/IA2), which demonstrated that this variant is a novel OsHV-1 microvariant: OsHV-1 µvar SD. A pilot transmission study indicates that OsHV-1 µvar SD is infectious with high viral loads ~ 7.57 × 10⁶ copies/mg detected in dead individuals. The detection of OsHV-1 µvar SD in a large port mirrors previous studies conducted in Australia where aquaculture farms and feral populations near port locations may be at a higher risk of OsHV-1 emergence. Further research is needed to understand the impacts of OsHV-1 µvar SD, such as transmission studies focusing on potential vectors and characterization of virulence as compared to other OsHV-1 µvars. To increase biosecurity of the global aquaculture industry, active and passive surveillance may be necessary to reduce spread of pathogens and make appropriate management decisions.

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.

Unraveling concordant and varying responses of oyster species to Ostreid Herpesvirus 1 variants

The Ostreid herpesvirus 1 (OsHV-1) and variants, particularly the microvariants (μVars), are virulent and economically devastating viruses impacting oysters. Since 2008 OsHV-1 μVars have emerged rapidly having particularly damaging effects on aquaculture industries in Europe, Australia and New Zealand. We conducted field trials in Tomales Bay (TB), California where a non-μVar strain of OsHV-1 is established and demonstrated differential mortality of naturally exposed seed of three stocks of Pacific oyster, Crassostrea gigas, and one stock of Kumamoto oyster, C. sikamea. Oysters exposed in the field experienced differential mortality that ranged from 64 to 99% in Pacific oysters (Tasmania>Midori = Willapa stocks), which was much higher than that of Kumamoto oysters (25%). Injection trials were done using French (FRA) and Australian (AUS) μVars with the same oyster stocks as planted in the field and, in addition, two stocks of the Eastern oyster, C. virginica. No mortality was observed in control oysters. One C. virginica stock suffered ~10% mortality when challenged with both μVars tested. Two Pacific oyster stocks suffered 75 to 90% mortality, while one C. gigas stock had relatively low mortality when challenged with the AUS μVar (~22%) and higher mortality when challenged with the French μVar (~72%). Conversely, C. sikamea suffered lower mortality when challenged with the French μVar (~22%) and higher mortality with the AUS μVar (~44%). All dead oysters had higher viral loads (~1000×) as measured by quantitative PCR relative to those that survived. However, some survivors had high levels of virus, including those from species with lower mortality. Field mortality in TB correlated with laboratory mortality of the FRA μVar (69% correlation) but not with that of the AUS μVar, which also lacked correlation with the FRA μVar. The variation in response to OsHV-1 variant challenges by oyster species and stocks demonstrates the need for empirical assessment of multiple OsHV-1 variants.

Fifty years of ecological changes: Regime shifts and drivers in a coastal Mediterranean lagoon during oligotrophication

Thau lagoon is a large Mediterranean coastal lagoons and it supports traditional shellfish farming activities. It has been subject to eutrophication leading to major anoxic events associated with massive mortalities of shellfish stocks. Since the 1970s, improvements have been made to wastewater treatment systems, which have gradually led to oligotrophication of the lagoon. The aim of our study was to determine how the decrease in nutrient inputs resulted in major ecological changes in Thau lagoon, by analysing five decades of time-series (1970-2018) of observations on pelagic and benthic autotrophic communities. We were able to identify two periods during the oligotrophication process. Period 1 (1970-1992) was considered a eutrophic period, characterised by the shift from seagrass dominance to dominance of red macroalgae. Period 2 (1993-2018), characterised by improved eutrophication status, was further divided into three: a transition phase (1993-2003) during which the water column continued to recover but the benthic community lagged behind in recovery and in partial resilience; a regime shift (2003-2006), after which the water column became oligotrophic and seagrass began to recover (2007-2018). Considering anoxia crises as indicators of ecosystem resilience and resistance, we used a generalised linear model to analyse meteorological and environmental data with the aim of identifying the triggers of summer anoxia over the study period. Among the meteorological variables studied, air temperature had the strongest positive effect, followed by the period and wind intensity (both negative effects) and by rainfall in July (positive effect). The risk of triggering anoxia was lower in period 2, evidence for the increasing resistance of the ecosystem to climatic stress throughout the oligotrophication process. At the ecosystem scale and in the long term perspective, the ecological gains related to oligotrophication are especially important in the context of climate change, with more frequent and severe heat waves predicted.

Latitudinal drivers of oyster mortality: deciphering host, pathogen and environmental risk factors

Diseases pose an ongoing threat to aquaculture, fisheries and conservation of marine species, and determination of risk factors of disease is crucial for management. Our objective was to decipher the effects of host, pathogen and environmental factors on disease-induced mortality of Pacific oysters (Crassostrea gigas) across a latitudinal gradient. We deployed young and adult oysters at 13 sites in France and we monitored survival, pathogens and environmental parameters. The young oysters came from either the wild collection or the hatchery while the adults were from the wild only. We then used Cox regression models to investigate the effect of latitude, site, environmental factors and origin on mortality risk and to extrapolate this mortality risk to the distribution limits of the species in Europe. We found that seawater temperature, food level, sea level atmospheric pressure, rainfall and wind speed were associated with mortality risk. Their effect on hatchery oysters was generally higher than on wild animals, probably reflecting that hatchery oysters were free of Ostreid herpesvirus 1 (OsHV-1) whereas those from the wild were asymptomatic carriers. The risk factors involved in young and adult oyster mortalities were different, reflecting distinct diseases. Mortality risk increases from 0 to 90% with decreasing latitude for young hatchery oysters, but not for young wild oysters or adults. Mortality risk was higher in wild oysters than in hatchery ones at latitude > 47.6°N while this was the opposite at lower latitude. Therefore, latitudinal gradient alters disease-induced mortality risk but interacts with the initial health status of the host and the pathogen involved. Practically, we suggest that mortality can be mitigated by using hatchery oysters in north and wild collected oysters in the south.

Prevalence and genotypic diversity of ostreid herpesvirus type 1 in Crassostrea gigas cultured in the Gulf of California, Mexico

In bivalve mollusk aquaculture, massive disease outbreaks with high mortality and large economic losses can occur, as in northwest Mexico in the 1990s. A range of pathogens can affect bivalves; one of great concern is ostreid herpesvirus 1 (OsHV-1), of which there are several strains. This virus has been detected in the Gulf of California in occasional or sporadic samplings, but to date, there have been few systematic studies. Monthly samples of Crassostrea gigas, water, and sediment were taken in the La Cruz coastal lagoon and analyzed by PCR. The native mollusk, Dosinia ponderosa, which lives outside the lagoon, was sampled as a control. The virus was found throughout the year only in C. gigas, with prevalence up to 60%. In total, 9 genotype variants were detected, and genetic analysis suggests that linear genotypic evolution has occurred from strain JF894308, present in La Cruz in 2011. There has been no evidence of the entry of new viral genotypes in the recent past, thus confinement of the virus within the lagoons of the Gulf of California could promote a native genotypic diversity in the short term.

Deciphering the effect of food availability, growth and host condition on disease susceptibility in a marine invertebrate

Food provisioning influences disease risk and outcome in animal populations in two ways. On the one hand, unrestricted food supply improves the physiological condition of the host and lowers its susceptibility to infectious disease, reflecting a trade-off between immunity and other fitness-related functions. On the other hand, food scarcity limits the resources available to the pathogen and slows the growth and metabolism of the host on which the pathogen depends to proliferate. Here, we investigated how food availability, growth rate and energetic reserves drive the outcome of a viral disease affecting an ecologically relevant model host, the Pacific oyster, Crassostrea gigas We selected fast- and slow-growing animals, and we exposed them to high and low food rations. We evaluated their energetic reserves, challenged them with a pathogenic virus, monitored daily survival and developed a mortality risk model. Although high food levels and oyster growth were associated with a higher risk of mortality, energy reserves were associated with a lower risk. Food availability acts both as an enabling factor for mortality by increasing oyster growth and as a limiting factor by increasing their energy reserves. This study clarifies how food resources have an impact on susceptibility to disease and indicates how the host's physiological condition could mitigate epidemics. Practically, we suggest that growth should be optimized rather than maximized, considering that trade-offs occur with disease resistance or tolerance.

Fine-scale temporal dynamics of herpes virus and vibrios in seawater during a polymicrobial infection in the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is currently being impacted by a polymicrobial disease that involves early viral infection by ostreid herpesvirus-1 (OsHV-1) followed by a secondary bacterial infection leading to death. A widely used method of inducing infection consists of placing specific pathogen-free oysters ('recipients') in cohabitation in the laboratory with diseased oysters that were naturally infected in the field ('donors'). With this method, we evaluated the temporal dynamics of pathogen release in seawater and the cohabitation time necessary for disease transmission and expression. We showed that OsHV-1 and Vibrio spp. in the seawater peaked concomitantly during the first 48 h and decreased thereafter. We found that 1.5 h of cohabitation with donors was enough time to transmit pathogens to recipients and to induce mortality later, reflecting the highly contagious nature of the disease. Finally, mortality of recipients was associated with increasing cohabitation time with donors until reaching a plateau at 20%. This reflects the cumulative effect of exposure to pathogens. The optimal cohabitation time was 5-6 d, the mortality of recipients occurring 1-2 d earlier.

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.