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.

Microbiota of the Digestive Glands and Extrapallial Fluids of Clams Evolve Differently Over Time Depending on the Intertidal Position.

The Manila clam (Ruditapes philippinarum) is the second most exploited bivalve in the world but remains threatened by diseases and global changes. Their associated microbiota play a key role in their fitness and acclimation capacities. This study aimed at better understanding the behavior of clam digestive glands and extrapallial fluids microbiota at small, but contrasting spatial and temporal scales. Results showed that environmental variations impacted clam microbiota differently according to the considered tissue. Each clam tissue presented its own microbiota and showed different dynamics according to the intertidal position and sampling period. Extrapallial fluids microbiota was modified more rapidly than digestive glands microbiota, for clams placed on the upper and lower intertidal position, respectively. Clam tissues could be considered as different microhabitats for bacteria as they presented different responses to small-scale temporal and spatial variabilities in natural conditions. These differences underlined a more stringent environmental filter capacity of the digestive glands.

Life history of oysters influences Vibrio parahaemolyticus accumulation in Pacific oysters (Crassostrea gigas).

Vibrio parahaemolyticus infection in humans is associated with raw oyster consumption. Evaluation of V. parahaemolyticus presence in oysters is of most interest because of the economic and public health issues that it represents. To explore V. parahaemolyticus accumulation and depuration in adult Crassostrea gigas, we developed a GFP-tagged V. parahaemolyticus strain (IFVp201-gfp+ ), as well as a rapid and efficient quantification method in C. gigas oysters haemolymph by flow cytometry. Impact of the life history of C. gigas on accumulation and depuration of V. parahaemolyticus IFVp201 was subsequently investigated. We found that naive oysters, i.e. grown in controlled facilities with UV treated seawater, accumulated significantly more IFVp201 than environmental oysters, i.e. grown in intertidal environment. We hypothesized that environmental oysters could have been immune primed, thus could limit V. parahaemolyticus accumulation. Meanwhile, both naive and environmental oysters had similar depuration rates.

Seaweeds influence oyster microbiota and disease susceptibility.

A growing awareness of role that microbiota can play in mediating the effects of pathogens on hosts has given rise to the concept of the pathobiome. Recently, we demonstrated that the Pacific oyster mortality syndrome affecting Crassostrea gigas oysters is caused by infection with the Ostreid herpesvirus type 1 (OsHV-1) followed by infection with multiple bacterial taxa. Here we extend the concept of this pathobiome beyond the host species and its bacterial microbiota by investigating how seaweed living in association with oysters influences their response to the disease. We hypothesized that by their mere presence in the environment, different species of seaweeds can positively or negatively influence the risk of disease in oysters by shaping their bacterial microbiota and their immune response. Although seaweed and oysters do not have direct ecological interactions, they are connected by seawater and likely share microbes. To test our hypothesis, oysters were acclimated with green, brown or red algae for 2 weeks and then challenged with OsHV-1. We monitored host survival and pathogen proliferation and performed bacterial microbiota and transcriptome analyses. We found that seaweeds can alter the bacterial microbiota of the host and its response to the disease. More particularly, green algae belonging to the genus Ulva spp. induced bacterial microbiota dysbiosis in oyster and modification of its transcriptional immune response leading to increased susceptibility to the disease. This work provides a better understanding of a marine disease and highlights the importance of considering both macrobiotic and microbiotic interactions for conservation, management and exploitation of marine ecosystems and resources.

Functional Diversification of Oyster Big Defensins Generates Antimicrobial Specificity and Synergy against Members of the Microbiota.

Big defensins are two-domain antimicrobial peptides (AMPs) that have highly diversified in mollusks. Cg-BigDefs are expressed by immune cells in the oyster Crassostrea gigas, and their expression is dampened during the Pacific Oyster Mortality Syndrome (POMS), which evolves toward fatal bacteremia. We evaluated whether Cg-BigDefs contribute to the control of oyster-associated microbial communities. Two Cg-BigDefs that are representative of molecular diversity within the peptide family, namely Cg-BigDef1 and Cg-BigDef5, were characterized by gene cloning and synthesized by solid-phase peptide synthesis and native chemical ligation. Synthetic peptides were tested for antibacterial activity against a collection of culturable bacteria belonging to the oyster microbiota, characterized by 16S sequencing and MALDI Biotyping. We first tested the potential of Cg-BigDefs to control the oyster microbiota by injecting synthetic Cg-BigDef1 into oyster tissues and analyzing microbiota dynamics over 24 h by 16S metabarcoding. Cg-BigDef1 induced a significant shift in oyster microbiota ß-diversity after 6 h and 24 h, prompting us to investigate antimicrobial activities in vitro against members of the oyster microbiota. Both Cg-BigDef1 and Cg-BigDef5 were active at a high salt concentration (400 mM NaCl) and showed broad spectra of activity against bacteria associated with C. gigas pathologies. Antimicrobial specificity was observed for both molecules at an intra- and inter-genera level. Remarkably, antimicrobial spectra of Cg-BigDef1 and Cg-BigDef5 were complementary, and peptides acted synergistically. Overall, we found that primary sequence diversification of Cg-BigDefs has generated specificity and synergy and extended the spectrum of activity of this peptide family.

Species-specific mechanisms of cytotoxicity toward immune cells determine the successful outcome of Vibrio infections.

Vibrio species cause infectious diseases in humans and animals, but they can also live as commensals within their host tissues. How Vibrio subverts the host defenses to mount a successful infection remains poorly understood, and this knowledge is critical for predicting and managing disease. Here, we have investigated the cellular and molecular mechanisms underpinning infection and colonization of 2 virulent Vibrio species in an ecologically relevant host model, oyster, to study interactions with marine Vibrio species. All Vibrio strains were recognized by the immune system, but only nonvirulent strains were controlled. We showed that virulent strains were cytotoxic to hemocytes, oyster immune cells. By analyzing host and bacterial transcriptional responses to infection, together with Vibrio gene knock-outs, we discovered that Vibrio crassostreae and Vibrio tasmaniensis use distinct mechanisms to cause hemocyte lysis. Whereas V. crassostreae cytotoxicity is dependent on a direct contact with hemocytes and requires an ancestral gene encoding a protein of unknown function, r5.7, V. tasmaniensis cytotoxicity is dependent on phagocytosis and requires intracellular secretion of T6SS effectors. We conclude that proliferation of commensal vibrios is controlled by the host immune system, preventing systemic infections in oysters, whereas the successful infection of virulent strains relies on Vibrio species-specific molecular determinants that converge to compromise host immune cell function, allowing evasion of the host immune system.

Differential basal expression of immune genes confers Crassostrea gigas resistance to Pacific oyster mortality syndrome.

BACKGROUND: As a major threat to the oyster industry, Pacific Oyster Mortality Syndrome (POMS) is a polymicrobial disease affecting the main oyster species farmed across the world. POMS affects oyster juveniles and became panzootic this last decade, but POMS resistance in some oyster genotypes has emerged. While we know some genetic loci associated with resistance, the underlying mechanisms remained uncharacterized. So, we developed a comparative transcriptomic approach using basal gene expression profiles between different oyster biparental families with contrasted phenotypes when confronted to POMS (resistant or susceptible). RESULTS: We showed that POMS resistant oysters show differential expression of genes involved in stress responses, protein modifications, maintenance of DNA integrity and repair, and immune and antiviral pathways. We found similarities and clear differences among different molecular pathways in the different resistant families. These results suggest that the resistance process is polygenic and partially varies according to the oyster genotype. CONCLUSIONS: We found differences in basal expression levels of genes related to TLR-NFκB, JAK-STAT and STING-RLR pathways. These differences could explain the best antiviral response, as well as the robustness of resistant oysters when confronted to POMS. As some of these genes represent valuable candidates for selective breeding, we propose future studies should further examine their function.

Selection of Vibrio crassostreae relies on a plasmid expressing a type 6 secretion system cytotoxic for host immune cells.

Pacific oyster mortality syndrome affects juveniles of Crassostrea gigas oysters and threatens the sustainability of commercial and natural stocks of this species. Vibrio crassostreae (V. crassostreae) has been repeatedly isolated from diseased animals, and the majority of the strains have been demonstrated to be virulent for oysters. In this study, we showed that oyster farms exhibited a high prevalence of a virulence plasmid carried by V. crassostreae, while oysters, at an adult stage, were reservoirs of this virulent population. The pathogenicity of V. crassostreae depends on a novel transcriptional regulator, which activates the bidirectional promoter of a type 6 secretion system (T6SS) genes cluster. Both the T6SS and a second chromosomal virulence factor, r5.7, are necessary for virulence but act independently to cause haemocyte (oyster immune cell) cytotoxicity. A phylogenetically closely related T6SS was identified in V. aestuarianus and V. tapetis, which infect adult oysters and clams respectively. We propose that haemocyte cytotoxicity is a lethality trait shared by a broad range of mollusc pathogens, and we speculate that T6SS was involved in parallel evolution of pathogen for molluscs.

High temperature induces transcriptomic changes in Crassostrea gigas that hinder progress of ostreid herpesvirus (OsHV-1) and promote survival.

Of all environmental factors, seawater temperature plays a decisive role in triggering marine diseases. Like fever in vertebrates, high seawater temperature could modulate the host response to pathogens in ectothermic animals. In France, massive mortality of Pacific oysters, Crassostrea gigas, caused by the ostreid herpesvirus 1 (OsHV-1) is markedly reduced when temperatures exceed 24°C in the field. In the present study we assess how high temperature influences the host response to the pathogen by comparing transcriptomes (RNA sequencing) during the course of experimental infection at 21°C (reference) and 29°C. We show that high temperature induced host physiological processes that are unfavorable to the viral infection. Temperature influenced the expression of transcripts related to the immune process and increased the transcription of genes related to the apoptotic process, synaptic signaling and protein processes at 29°C. Concomitantly, the expression of genes associated with catabolism, metabolite transport, macromolecule synthesis and cell growth remained low from the first stage of infection at 29°C. Moreover, viral entry into the host might have been limited at 29°C by changes in extracellular matrix composition and protein abundance. Overall, these results provide new insights into how environmental factors modulate host-pathogen interactions.

High temperature induces transcriptomic changes in Crassostrea gigas that hinders progress of Ostreid herpesvirus (OsHV-1) and promotes survival.

Among all the environmental factors, seawater temperature plays a decisive role in triggering marine diseases. Like fever in vertebrates, high seawater temperature could modulate the host response to the pathogens in ectothermic animals. In France, massive mortality of Pacific oysters Crassostrea gigas caused by the ostreid herpesvirus 1 (OsHV-1) is markedly reduced when temperatures exceed 24°C in the field. In the present study we assess how high temperature influences the host response to the pathogen by comparing transcriptomes (RNA-sequencing) during the course of experimental infection at 21°C (reference) and 29°C. We show that high temperature induced host physiological processes that are unfavorable to the viral infection. Temperature influenced the expression of transcripts related to the immune process and increased the transcription of genes related to apoptotic process, synaptic signaling, and protein processes at 29°C. Concomitantly, the expression of genes associated to catabolism, metabolites transport, macromolecules synthesis and cell growth remained low since the first stage of infection at 29°C. Moreover, viral entry into the host might have been limited at 29°C by changes in extracellular matrix composition and protein abundance. Overall, these results provide new insights into how environmental factors modulate the host-pathogen interactions.

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.

Temperature modulate disease susceptibility of the Pacific oyster Crassostrea gigas and virulence of the Ostreid herpesvirus type 1.

Temperature triggers marine diseases by changing host susceptibility and pathogen virulence. Oyster mortalities associated with the Ostreid herpesvirus type 1 (OsHV-1) have occurred seasonally in Europe when the seawater temperature range reaches 16-24 °C. Here we assess how temperature modulates oyster susceptibility to OsHV-1 and pathogen virulence. Oysters were injected with OsHV-1 suspension incubated at 21 °C, 26 °C and 29 °C and were placed in cohabitation with healthy oysters (recipients) at these three temperatures according to a fractional factorial design. Survival was followed for 14 d and recipients were sampled for OsHV-1 DNA quantification and viral gene expression. The oysters were all subsequently placed at 21 °C to evaluate the potential for virus reactivation, before being transferred to oyster farms to evaluate their long-term susceptibility to the disease. Survival of recipients at 29 °C (86%) was higher than at 21 °C (52%) and 26 °C (43%). High temperature (29 °C) decreased the susceptibility of oysters to OsHV-1 without altering virus infectivity and virulence. At 26 °C, the virulence of OsHV-1 was enhanced. Differences in survival persisted when the recipients were all placed at 21 °C, suggesting that OsHV-1 did not reactivate. Additional oyster mortality followed the field transfer, but the overall survival of oysters infected at 29 °C remained higher.

Long dsRNAs promote an anti-viral response in Pacific oyster hampering ostreid herpesvirus 1 replication.

Double-stranded RNA (dsRNA)-mediated genetic interference (RNAi) is a widely used reverse genetic tool for determining the loss-of-function phenotype of a gene. Here, the possible induction of an immune response by long dsRNA was tested in a marine bivalve (Crassostrea gigas), as well as the specific role of the subunit 2 of the nuclear factor κB inhibitor (IκB2). This gene is a candidate of particular interest for functional investigations in the context of oyster mass mortality events, as Cg-IκB2 mRNA levels exhibited significant variation depending on the amount of ostreid herpesvirus 1 (OsHV-1) DNA detected. In the present study, dsRNAs targeting Cg-IκB2 and green fluorescent protein genes were injected in vivo into oysters before being challenged by OsHV-1. Survival appeared close to 100% in both dsRNA-injected conditions associated with a low detection of viral DNA and a low expression of a panel of 39 OsHV-1 genes as compared with infected control. Long dsRNA molecules, both Cg-IκB2- and GFP-dsRNA, may have induced an anti-viral state controlling the OsHV-1 replication and precluding the understanding of the specific role of Cg-IκB2 Immune-related genes including Cg-IκB1, Cg-Rel1, Cg-IFI44, Cg-PKR and Cg-IAP appeared activated in the dsRNA-injected condition, potentially hampering viral replication and thus conferring a better resistance to OsHV-1 infection. We revealed that long dsRNA-mediated genetic interference triggered an anti-viral state in the oyster, emphasizing the need for new reverse genetics tools for assessing immune gene function and avoiding off-target effects in bivalves.

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.

A core of functional complementary bacteria infects oysters in Pacific Oyster Mortality Syndrome.

BACKGROUND: The Pacific oyster Crassostrea gigas is one of the main cultivated invertebrate species worldwide. Since 2008, oyster juveniles have been confronted with a lethal syndrome known as the Pacific Oyster Mortality Syndrome (POMS). POMS is a polymicrobial disease initiated by a primary infection with the herpesvirus OsHV-1 µVar that creates an oyster immunocompromised state and evolves towards a secondary fatal bacteremia. RESULTS: In the present article, we describe the implementation of an unprecedented combination of metabarcoding and metatranscriptomic approaches to show that the sequence of events in POMS pathogenesis is conserved across infectious environments. We also identified a core bacterial consortium which, together with OsHV-1 µVar, forms the POMS pathobiota. This bacterial consortium is characterized by high transcriptional activities and complementary metabolic functions to exploit host's resources. A significant metabolic specificity was highlighted at the bacterial genus level, suggesting low competition for nutrients between members of the core bacteria. CONCLUSIONS: Lack of metabolic competition between the core bacteria might favor complementary colonization of host tissues and contribute to the conservation of the POMS pathobiota across distinct infectious environments.

Epigenetic variations are more substantial than genetic variations in rapid adaptation of oyster to Pacific oyster mortality syndrome.

Disease emergence is accelerating with global changes. Understanding by which mechanisms host populations can rapidly adapt will be crucial for management practices. Pacific oyster mortality syndrome (POMS) imposes a substantial and recurrent selective pressure on oyster populations, and rapid adaptation may arise through genetics and epigenetics. In this study, we used (epi)genome-wide association mapping to show that oysters differentially exposed to POMS displayed genetic and epigenetic signatures of selection. Consistent with higher resistance to POMS, the genes targeted included many genes in several pathways related to immunity. By combining correlation, DNA methylation quantitative trait loci, and variance partitioning, we revealed that a third of phenotypic variation was explained by interactions between the genetic and epigenetic information, ~14% by the genome, and up to 25% by the epigenome alone. Similar to genetically based adaptation, epigenetic mechanisms notably governing immune responses can contribute substantially to the rapid adaptation of hosts to emerging infectious diseases.

m6A Profile Dynamics Indicates Regulation of Oyster Development by m6A-RNA Epitranscriptomes.

The N6-methylation of RNA adenosines (N6-methyladenosine, m6A) is an important regulator of gene expression with critical implications in vertebrate and insect development. However, the developmental significance of epitranscriptomes in lophotrochozoan organisms remains unknown. Using methylated RNA immunoprecipitation sequencing (MeRIP-seq), we generated transcriptome-wide m6A-RNA methylomes covering the whole development of the oyster from oocytes to juveniles. Oyster RNA classes display specific m6A signatures, with messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) exhibiting distinct profiles and being highly methylated compared to transposable element (TE) transcripts. Epitranscriptomes are dynamic and correspond to the chronological steps of development (cleavage, gastrulation, organogenesis, and metamorphosis), with minimal mRNA and lncRNA methylation at the morula stage followed by a global increase. mRNA m6A levels are correlated to transcript levels, and shifts in methylation profiles correspond to expression kinetics. Differentially methylated transcripts cluster according to embryo-larval stages and bear the corresponding developmental functions (cell division, signal transduction, morphogenesis, and cell differentiation). The m6A level of TE transcripts is also regulated and peaks during the gastrulation. We demonstrate that m6A-RNA methylomes are dynamic and associated with gene expression regulation during oyster development. The putative epitranscriptome implication in the cleavage, maternal-to-zygotic transition, and cell differentiation in a lophotrochozoan model brings new insights into the control and evolution of developmental processes.

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.

Early life microbial exposures shape the Crassostrea gigas immune system for lifelong and intergenerational disease protection.

BACKGROUND: The interaction of organisms with their surrounding microbial communities influences many biological processes, a notable example of which is the shaping of the immune system in early life. In the Pacific oyster, Crassostrea gigas, the role of the environmental microbial community on immune system maturation - and, importantly, protection from infectious disease - is still an open question. RESULTS: Here, we demonstrate that early life microbial exposure durably improves oyster survival when challenged with the pathogen causing Pacific oyster mortality syndrome (POMS), both in the exposed generation and in the subsequent one. Combining microbiota, transcriptomic, genetic, and epigenetic analyses, we show that the microbial exposure induced changes in epigenetic marks and a reprogramming of immune gene expression leading to long-term and intergenerational immune protection against POMS. CONCLUSIONS: We anticipate that this protection likely extends to additional pathogens and may prove to be an important new strategy for safeguarding oyster aquaculture efforts from infectious disease. tag the videobyte/videoabstract in this section Video Abstract.