Shotgun analysis to identify differences in protein expression between granulocytes and hyalinocytes of the European flat oyster Ostrea edulis

Recovery of wild populations of the European flat oyster Ostrea edulis is important for ecosystem health and conservation of this species, because native oyster populations have dramatically declined or disappeared in most European waters. Diseases have contributed to oyster decline and are important constrains for oyster recovery. Understanding oyster immune system should contribute to design effective strategies to fight oyster diseases. Haemocytes play a pivotal role in mollusc immune responses protecting from infection. Two main types of haemocytes, granulocytes and hyalinocytes, are distinguished in O. edulis. A study aiming to explore differential functions between both haemocyte types and, thus, to enrich the knowledge of Ostrea edulis immune system, was performed by comparing the proteome of the two haemolymph cell types, using a shotgun approach through liquid chromatography (LC) coupled to mass spectrometry (MS). Cells from oyster haemolymph were differentially separated by Percoll density gradient centrifugation. Shotgun LC-MS/MS performance allowed the identification of 145 proteins in hyalinocytes and 138 in the proteome of granulocytes. After a comparative analysis, 55 proteins with main roles in defence were identified, from which 28 were representative of granulocytes and 27 of hyalinocytes, plus 11 proteins shared by both cell types. Different proteins involved in signal transduction, apoptosis, oxidative response, processes related with the cytoskeleton and structure, recognition and wound healing were identified as representatives of each haemocyte type. Important signalling pathways in the immune response such as MAPK, Ras and NF-κβ seemed to be more relevant for granulocytes, while the Wnt signalling pathway, particularly relevant for wound healing, more relevant in hyalinocytes. The differences in proteins involved in recognition and in cytoskeleton and structure suggest differential specialisation in processes of phagocytosis and internalisation of pathogens between haemocyte types. Apoptosis seemed more active in granulocytes. The differences in proteins involved in oxidative response also suggest different redox processes in each cell type.

Effects of seawater salinity and pH on cellular metabolism and enzyme activities in biomineralizing tissues of marine bivalves

Molluscan shell formation is a complex energy demanding process sensitive to the shifts in seawater CaCO₃ saturation due to changes in salinity and pH. We studied the effects of salinity and pH on energy demand and enzyme activities of biomineralizing cells of the Pacific oyster (Crassostrea gigas) and the hard-shell clam (Mercenaria mercenaria). Adult animals were exposed for 14 days to high (30), intermediate (18), or low (10) salinity at either high (8.0-8.2) or low (7.8) pH. Basal metabolic cost as well as the energy cost of the biomineralization-related cellular processes were determined in isolated mantle edge cells and hemocytes. The total metabolic rates were similar in the hemocytes of the two studied species, but considerably higher in the mantle cells of C. gigas compared with those of M. mercenaria. Cellular respiration was unaffected by salinity in the clams' cells, while in oysters' cells the highest respiration rate was observed at intermediate salinity (18). In both studied species, low pH suppressed cellular respiration. Low pH led to an upregulation of Na⁺/K⁺ ATPase activity in biomineralizing cells of oysters and clams. Activities of Ca²⁺ ATPase and H⁺ ATPase, as well as the cellular energy costs of Ca²⁺ and H⁺ transport in the biomineralizing cells were insensitive to the variation in salinity and pH in the two studied species. Variability in cellular response to low salinity and pH indicates that the disturbance of shell formation under these conditions has different underlying mechanisms in the two studied species.

Histopathological and Molecular Study of Pacific Oyster Tissues Provides Insights into V. aestuarianus Infection Related to Oyster Mortality

Consumer preference for healthy and sustainable food products has been steadily increasing in recent years. Bivalve mollusks satisfy these characteristics and have captured ever-increasing market shares. However, the expansion of molluscan culture in worldwide and global trade have favored the spread of pathogens around the world. Combined with environmental changes and intensive production systems this has contributed to the occurrence of mass mortality episodes, thus posing a threat to the production of different species, including the Pacific oyster Crassotrea gigas. In the San Teodoro lagoon, one of the most devoted lagoons to extensive Pacific oyster aquaculture in Sardinia, a mortality outbreak was observed with an estimated 80% final loss of animal production. A study combining cultural, biomolecular and histopathological methods was conducted: (1) to investigate the presence of different Vibrio species and OsHV-1 in selected oyster tissues (digestive gland, gills, and mantle); (2) to quantify Vibrio aestuarianus and to evaluate the severity of hemocyte infiltration in infected tissues; (3) to produce post-amplification data and evaluating ToxR gene as a target for phylogenetic analyses. Results provide new insights into V. aestuarianus infection related to oyster mortality outbreaks and pave the way to the development of tools for oyster management.

Molecular characterization and functional analysis of a Rel gene in the Pacific oyster

The nuclear factor-κB (NF-κB) signaling pathway plays a crucial role in regulating many physiological processes such as development, inflammation, apoptosis, cell proliferation, differentiation and immune responses. And the NF-κB/Rel family members were considered as the most important transcription factors in the NF-κB signaling pathway. In this study, we cloned a Rel homolog gene (named as CgRel2) from the Pacific oyster, Crassostrea gigas. The 2115-bp open reading frame (ORF) encodes 704 amino acids and CgRel2 possesses a conserved Rel Homology Domain (RHD) at the N-terminus. Phylogenetic analysis revealed that CgRel2 is most closely related to Pinctada fucata dorsal protein. CgRel2 transcripts are widely expressed in all tested tissues, with the highest expression observed in the labial palp and the gill. Moreover, the expression of CgRel2 is significantly upregulated after lipopolysaccharide (LPS), peptidoglycan (PGN), and polyinosinic-polycytidylic acid [poly(I:C)] challenge. CgRel2 transfection into human cell lines activated NF-κB, TNFα and oyster IL-17 (CgIL-17) reporter genes in a dose-dependent manner, while CgRel2 overexpression cannot induce ISRE (Interferon stimulation response element) reporter gene's transcriptional activity. Additionally, the results of co-immunoprecipitation showed that CgRel2 or CgRel1 could interact with oyster IκB1, IκB2 and IκB3 proteins strongly, which may be critical for the immune signaling transduction and the regulation of its immune functions. Together, these results suggest that CgRel2 could respond to pathogenic infection, participate in the immune signal transduction and activate NF-κB, TNFα and CgIL-17 reporter genes. Thus, CgRel2 could play an important role in the oyster immune system.

Genome-wide identification, characterization of RLR genes in Yesso scallop (Patinopecten yessoensis) and functional regulations in responses to ocean acidification

Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs), are crucial sensors with a conserved structure in cytoplasm, inducing the production of cytokines, chemokines and host restriction factors which mediate a variety of intracellular activities to interfere with distinct PAMPs (pathogen-associated molecular patterns) for eliminating pathogens in innate immune system. Although RLR genes have been investigated in most vertebrates and some invertebrates, the systematic identification and characterization of RLR genes have not been reported in scallops. In this study, four RLR genes (PY-10413.4, PY-10413.5, PY-443.7 and PY-443.8, designated PyRLRs) were identified in Yesso scallop (Patinopecten yessoensis) through whole-genome scanning through in silico analysis, including two pairs of tandem duplicate genes located on the same scaffold (PY-10413.4 and PY-10413.5, PY-443.7 and PY-443.8, respectively). Phylogenetic and protein structural analyses were performed to determine the identities and evolutionary relationships of these genes. The expression profiles of PyRLRs were determined in all developmental stages, in healthy adult tissues, and in mantles that simulated ocean acidification (OA) exposure (pH = 6.5 and 7.5) at different time points (3, 6, 12 and 24 h). Spatiotemporal expression patterns suggested the functional roles of PyRLRs in all stages of development and growth of the scallop. Regulation expressions revealed PY-10413.4 and PY-10413.5 with one or two CARD(s) (caspase activation and recruitment domain) were up-regulated expressed at most time points, whereas PY-443.8 and PY-10413.4 without CARD were significantly down-regulated at each time points, suggesting functional differentiations in the two pairs of PyRLRs based on the structural differences in response to OA. Collectively, this study demonstrated gene duplication of RLR family genes and provide primary analysis for versatile roles in the response of the bivalve innate immune system to OA challenge.

IKKε-like plays an important role in the innate immune signaling of the Pacific oyster (Crassostrea gigas)

IκB-related kinase ε (IKKε) plays a crucial role in the activation of nuclear factor κB (NF-κB) by phosphorylating inhibitor of NF-κB (IκB) and in the regulation of interferon (IFN) gene expression by phosphorylating IFN regulatory factors (IRFs). In this study, we cloned an IKKε homologue cDNA (designated as CgIKKε-like) from the Pacific oyster, Crassostrea gigas. The full 2896-bp cDNA sequence comprised a 2163-bp open reading frame (ORF) encoding 720 amino acids. CgIKKε-like is ubiquitously expressed, and its mRNA levels in hemocytes after poly I:C, V. alginolyticus, or OsHV-1 μVar challenge were analyzed by real-time PCR. Compared to that in the control, CgIKKε-like mRNA expression levels were significantly increased at 3 h and peaked at 6 h after OsHV-1 μVar challenge; no obvious changes were observed in expression levels until 24 h after either V. alginolyticus or poly I:C challenge, reaching a maximum at 24 h (p < 0.01) and then rapidly decreasing. CgIKKε-like transfection into human cell lines induced NF-κB and ISRE activation, while transfection with CgIKKε-like deletion mutants abolished NF-κB and ISRE reporter gene activation. Additionally, CgIKKε-like could interact with CgTBK1 and could form homodimers strongly, which may be critical for the immune signaling transduction. Last but not least, we found that CgIKKε-like may increase CgIκBs phosphorylation and could interact with CgIRF8. Together, these results suggest that CgIKKε-like could respond to pathogenic infection, participate in the immune signal transduction and activate NF-κB and ISRE reporter genes. Thus, CgIKKε-like could play an important role in the oyster immune system.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

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.

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

Inefficient immune response is associated with microbial permissiveness in juvenile oysters affected by mass mortalities on field

Since 2008, juvenile Crassostrea gigas oysters have suffered from massive mortalities in European farming areas. This disease of complex etiology is still incompletely understood. Triggered by an elevated seawater temperature, it has been associated to infections by a herpes virus named OsHV-1 as well as pathogenic vibrios of the Splendidus clade. Ruling out the complexity of the disease, most of our current knowledge has been acquired in controlled experiments. Among the many unsolved questions, it is still ignored what role immunity plays in the capacity oysters have to survive an infectious episode. Here we show that juvenile oysters susceptible to the disease mount an inefficient immune response associated with microbial permissiveness and death. We found that, in contrast to resistant adult oysters having survived an earlier episode of mortality, susceptible juvenile oysters never exposed to infectious episodes died by more than 90% in a field experiment. Susceptible oysters were heavily colonized by OsHV-1 herpes virus as well as bacteria including vibrios potentially pathogenic for oysters, which proliferated in oyster flesh and body fluids during the mortality event. Nonetheless, susceptible oysters were found to sense microbes as indicated by an overexpression of immune receptors and immune signaling pathways. However, they did not express important immune effectors involved in antimicrobial immunity and apoptosis and showed repressed expression of genes involved in ROS and metal homeostasis. This contrasted with resistant oysters, which expressed those important effectors, controlled bacterial and viral colonization and showed 100% survival to the mortality event. Altogether, our results demonstrate that the immune response mounted by susceptible oysters lacks some important immune functions and fails in controlling microbial proliferation. This study opens the way to more holistic studies on the "mass mortality syndrome", which are now required to decipher the sequence of events leading to oyster mortalities and determine the relative weight of pathogens, oyster genetics and oyster-associated microbiota in the disease.

Transepithelial migration of mucosal hemocytes in Crassostrea virginica and potential role in Perkinsus marinus pathogenesis

We have recently described the presence of hemocytes associated with mucus covering the pallial organs (mantle, gills, and body wall) 3 of the eastern oyster Crassostrea virginica. These hemocytes, hereby designated "pallial hemocytes" share common general characteristics with circulating hemocytes but also display significant differences particularly in their cell surface epitopes. The specific location of pallial hemocytes as peripheral cells exposed directly to the marine environment confers them a putative sentinel role. The purpose of this study was to gain a better understanding of the source of these pallial hemocytes by evaluating possible exchanges between circulatory and pallial hemocyte populations and whether these exchanges are regulated by pathogen exposure. Bi-directional transepithelial migrations of hemocytes between pallial surfaces and the circulatory system were monitored using standard cell tracking approaches after staining with the vital fluorescent dye carboxyfluorescein diacetate succinimidyl ester (CFSE) in conjunction with fluorescent microscopy and flow cytometry. Results showed bi-directional migration of hemocytes between both compartments and suggest that hemocyte migration from the pallial mucus layer to the circulatory system may occur at a greater rate compared to migration from the circulatory system to the pallial mucus layer, further supporting the role of pallial hemocytes as sentinel cells. Subsequently, the effect of the obligate parasite Perkinsus marinus and the opportunistic pathogen Vibrio alginolyticus on transepithelial migration of oyster hemocytes was investigated. Results showed an increase in hemocyte migration in response to P. marinus exposure. Furthermore, P. marinus cells were acquired by pallial hemocytes before being visible in underlying tissues and the circulatory system suggesting that this parasite could use pallial hemocytes as a vehicle facilitating its access to oyster tissues. These results are discussed in light of new evidence highlighting the role of oyster pallial organs as a portal for the initiation of P. marinus infections in oysters.

Integrating Discovery-driven Proteomics and Selected Reaction Monitoring To Develop a Noninvasive Assay for Geoduck Reproductive Maturation

Geoduck clams (Panopea generosa) are an increasingly important fishery and aquaculture product along the eastern Pacific coast from Baja California, Mexico, to Alaska. These long-lived clams are highly fecund, although sustainable hatchery production of genetically diverse larvae is hindered by the lack of sexual dimorphism, resulting in asynchronous spawning of broodstock, unequal sex ratios, and low numbers of breeders. The development of assays of gonad physiology could indicate sex and maturation stage as well as be used to assess the status of natural populations. Proteomic profiles were determined for three reproductive maturation stages in both male and female clams using data-dependent acquisition (DDA) of gonad proteins. Gonad proteomes became increasingly divergent between males and females as maturation progressed. The DDA data were used to develop targets analyzed with selected reaction monitoring (SRM) in gonad tissue as well as hemolymph. The SRM assay yielded a suite of indicator peptides that can be used as an efficient assay to determine geoduck gonad maturation status. Application of SRM in hemolymph samples demonstrates that this procedure could effectively be used to assess reproductive status in marine mollusks in a nonlethal manner.

High throughput sequencing of RNA transcriptomes in Ruditapes philippinarum identifies genes involved in osmotic stress response

Ruditapes philippinarum, is an economically important marine bivalve species. The ability to cope with low salinity stress is quite important for the survival of aquatic species under natural conditions. In this study, the transcriptional response of the Manila clam to low salinity stress was characterized using RNA sequencing. The transcriptomes of a low salinity-treatment group (FRp1, FRp2), which survived under low salinity stress, and control group (SRp1, SRp2), which was not subjected to low salinity stress, were sequenced with the Illumina HiSeq platform. A total of 196,578 unigenes were generated. GO and KEGG analyses revealed that signal transduction, immune response, cellular component organization or biogenesis, and energy production processes were the most highly enriched pathways among the genes that were differentially expressed under low salinity stress. All these pathways could be assigned to the following biological functions in the low salinity tolerant Manila clam: signal response to low salinity stress, antioxidant response, intracellular free amino acid transport and metabolism, energy production and conversion, cell signaling pathways, and regulation of ionic content and cell volume. In summary, this is the first study using high-throughput sequencing to identify gene targets that could explain osmotic regulation mechanisms under salinity stress in R. philippinarum.

Epipodial Tentacle Gene Expression and Predetermined Resilience to Summer Mortality in the Commercially Important Greenlip Abalone, Haliotis laevigata

"Summer mortality" is a phenomenon that occurs during warm water temperature spikes that results in the mass mortality of many ecologically and economically important mollusks such as abalone. This study aimed to determine whether the baseline gene expression of abalone before a laboratory-induced summer mortality event was associated with resilience to summer mortality. Tentacle transcriptomes of 35 greenlip abalone (Haliotis laevigata) were sequenced prior to the animals being exposed to an increase in water temperature-simulating conditions which have previously resulted in summer mortality. Abalone derived from three source locations with different environmental conditions were categorized as susceptible or resistant to summer mortality depending on whether they died or survived after the water temperature was increased. We detected two genes showing significantly higher expression in resilient abalone relative to susceptible abalone prior to the laboratory-induced summer mortality event. One of these genes was annotated through the NCBI non-redundant protein database using BLASTX to an anemone (Exaiptasia pallida) Transposon Ty3-G Gag Pol polyprotein. Distinct gene expression signatures were also found between resilient and susceptible abalone depending on the population origin, which may suggest divergence in local adaptation mechanisms for resilience. Many of these genes have been suggested to be involved in antioxidant and immune-related functions. The identification of these genes and their functional roles have enhanced our understanding of processes that may contribute to summer mortality in abalone. Our study supports the hypothesis that prestress gene expression signatures are indicative of the likelihood of summer mortality.

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

Double stranded RNA-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, i.e. 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 massive mortality oyster 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 Fluorescence 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 to 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 Cg-IκB2 specific role. Immune-related genes including Cg-IκB1, Cg-Rel1, Cg-IFI44, Cg-PKR, and Cg-IAP appeared activated in 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 of new reverse genetics tools for assessing immune gene function and avoiding off-target effects in bivalves.

Construction of an Ostrea edulis database from genomic and expressed sequence tags (ESTs) obtained from Bonamia ostreae infected haemocytes: Development of an immune-enriched oligo-microarray

The flat oyster, Ostrea edulis, is one of the main farmed oysters, not only in Europe but also in the United States and Canada. Bonamiosis due to the parasite Bonamia ostreae has been associated with high mortality episodes in this species. This parasite is an intracellular protozoan that infects haemocytes, the main cells involved in oyster defence. Due to the economical and ecological importance of flat oyster, genomic data are badly needed for genetic improvement of the species, but they are still very scarce. The objective of this study is to develop a sequence database, OedulisDB, with new genomic and transcriptomic resources, providing new data and convenient tools to improve our knowledge of the oyster's immune mechanisms. Transcriptomic and genomic sequences were obtained using 454 pyrosequencing and compiled into an O. edulis database, OedulisDB, consisting of two sets of 10,318 and 7159 unique sequences that represent the oyster's genome (WG) and de novo haemocyte transcriptome (HT), respectively. The flat oyster transcriptome was obtained from two strains (naïve and tolerant) challenged with B. ostreae, and from their corresponding non-challenged controls. Approximately 78.5% of 5619 HT unique sequences were successfully annotated by Blast search using public databases. A total of 984 sequences were identified as being related to immune response and several key immune genes were identified for the first time in flat oyster. Additionally, transcriptome information was used to design and validate the first oligo-microarray in flat oyster enriched with immune sequences from haemocytes. Our transcriptomic and genomic sequencing and subsequent annotation have largely increased the scarce resources available for this economically important species and have enabled us to develop an OedulisDB database and accompanying tools for gene expression analysis. This study represents the first attempt to characterize in depth the O. edulis haemocyte transcriptome in response to B. ostreae through massively sequencing and has aided to improve our knowledge of the immune mechanisms of flat oyster. The validated oligo-microarray and the establishment of a reference transcriptome will be useful for large-scale gene expression studies in this species.

The systematic regulation of oyster CgIL17-1 and CgIL17-5 in response to air exposure

As a proinflammatory cytokine, vertebrate interleukin 17 (IL17) plays a vital role in the balance of inflammation and homeostasis, and is involved in a systemic regulation of glucose homeostasis. In the present study, a remarkable increase of glucose concentration was observed in oyster serum after 2 d air exposure, which was followed by a rapid up-regulation of CgIL17-1 and CgIL17-5. After oysters was received an injection of extra glucose, the mRNA expressions of CgIL17-1 and CgIL17-5 were also significantly up-regulated. The histopathological changes of hepatopancreas were observed after the oysters were treated by the recombinant proteins of CgIL17-1 and CgIL17-5 in vivo or subjected to air exposure. A significant decrease of GSK3β (Glycogen synthase kinase-3β) protein was also observed after the injection of CgIL17-1 and CgIL17-5 recombinant proteins in vivo. When the oysters with CgIL17-1 and CgIL17-5 genes knocked down were subjected to air exposure, the decline of GSK3β concentration was slowed down and it could still be obviously detected after 7 d compared with that in the control. Meanwhile, the expression of CgDefensin and CgDFFA was inhibited, while CgIAP was up-regulated when CgIL17-1 and CgIL17-5 genes were knocked down, and the oysters exhibited higher mortality (p < 0.05) at 3 d, whereas lower at the late stage of air exposure compared with that in the controls. The results collectively suggested that once oysters were exposed to air, the synthesis of proinflammatory cytokines CgIL17-1 and CgIL17-5 was induced by the up-regulated glucose concentration in oyster serum, which would be not only a negative feedback to the high glucose concentration through mediating the regulation of GSK3β, but also an inducer on tissue damage and immunocompetence as well as the adaptability to stresses.

Distinct immune responses of juvenile and adult oysters (Crassostrea gigas) to viral and bacterial infections

Since 2008, massive mortality events of Pacific oysters (Crassostrea gigas) have been reported worldwide and these disease events are often associated with Ostreid herpesvirus type 1 (OsHV-1). Epidemiological field studies have also reported oyster age and other pathogens of the Vibrio genus are contributing factors to this syndrome. We undertook a controlled laboratory experiment to simultaneously investigate survival and immunological response of juvenile and adult C. gigas at different time-points post-infection with OsHV-1, Vibrio tasmaniensis LGP32 and V. aestuarianus. Our data corroborates epidemiological studies that juveniles are more susceptible to OsHV-1, whereas adults are more susceptible to Vibrio. We measured the expression of 102 immune-genes by high-throughput RT-qPCR, which revealed oysters have different transcriptional responses to OsHV-1 and Vibrio. The transcriptional response in the early stages of OsHV-1 infection involved genes related to apoptosis and the interferon-pathway. Transcriptional response to Vibrio infection involved antimicrobial peptides, heat shock proteins and galectins. Interestingly, oysters in the later stages of OsHV-1 infection had a transcriptional response that resembled an antibacterial response, which is suggestive of the oyster’s microbiome causing secondary infections (dysbiosis-driven pathology). This study provides molecular evidence that oysters can mount distinct immune response to viral and bacterial pathogens and these responses differ depending on the age of the host.

Oysters and Vibrios as a Model for Disease Dynamics in Wild Animals

Disease dynamics in the wild are influenced by a number of ecological and evolutionary factors not addressed by traditional laboratory-based characterization of pathogens. Here we propose the oyster, Crassostrea gigas, as a model for studying the interaction of the environment, bacterial pathogens, and the host in disease dynamics. We show that an important first step is to ask whether the functional unit of pathogenesis is a bacterial clone, a population, or a consortium in order to assess triggers of disease outbreaks and devise appropriate monitoring tools. Moreover, the development of specific-pathogen-free (SPF) oysters has enabled assessment of the infection process under natural conditions. Finally, recent results show the importance of microbial interactions and host genetics in determining oyster health and disease.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

GigaTON: an extensive publicly searchable database providing a new reference transcriptome in the pacific oyster Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is one of the most important aquaculture shellfish resources worldwide. Important efforts have been undertaken towards a better knowledge of its genome and transcriptome, which makes now C. gigas becoming a model organism among lophotrochozoans, the under-described sister clade of ecdysozoans within protostomes. These massive sequencing efforts offer the opportunity to assemble gene expression data and make such resource accessible and exploitable for the scientific community. Therefore, we undertook this assembly into an up-to-date publicly available transcriptome database: the GigaTON (Gigas TranscriptOme pipeliNe) database.

DESCRIPTION

We assembled 2204 million sequences obtained from 114 publicly available RNA-seq libraries that were realized using all embryo-larval development stages, adult organs, different environmental stressors including heavy metals, temperature, salinity and exposure to air, which were mostly performed as part of the Crassostrea gigas genome project. This data was analyzed in silico and resulted into 56621 newly assembled contigs that were deposited into a publicly available database, the GigaTON database. This database also provides powerful and user-friendly request tools to browse and retrieve information about annotation, expression level, UTRs, splice and polymorphism, and gene ontology associated to all the contigs into each, and between all libraries.

CONCLUSIONS

The GigaTON database provides a convenient, potent and versatile interface to browse, retrieve, confront and compare massive transcriptomic information in an extensive range of conditions, tissues and developmental stages in Crassostrea gigas. To our knowledge, the GigaTON database constitutes the most extensive transcriptomic database to date in marine invertebrates, thereby a new reference transcriptome in the oyster, a highly valuable resource to physiologists and evolutionary biologists.

Genome-wide identification and characterization of TRAF genes in the Yesso scallop (Patinopecten yessoensis) and their distinct expression patterns in response to bacterial challenge

The tumor necrosis factor (TNF) receptor associated factors (TRAFs) are the major signal transducers for the TNF receptor superfamily and the interleukin-1 receptor/Toll-like receptor (IL-1R/TLR) superfamily, which regulate a variety of cellular activities and innate immune responses. TRAF genes have been extensively studied in various species, including vertebrates and invertebrates. However, as one of the key component of NF-κB pathway, TRAF genes have not been systematically characterized in marine invertebrates. In this study, we identified and characterized five TRAF genes, PyTRAF2, PyTRAF3, PyTRAF4, PyTRAF6 and PyTRAF7, in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine their identities and evolutionary relationships. In comparison with the TRAF genes from vertebrate species, the structural features were all relatively conserved in the PyTRAF genes. To gain insights into the roles of TRAF genes during scallop innate immune responses, quantitative real-time PCR was used to investigate the expression profiles in the different stages of scallop development, in the healthy adult tissues, and in the hemocytes after bacterial infection with Micrococcus luteus and Vibrio anguillarum. Based on the qRT-PCR analysis, the expression of most of the PyTRAFs was significantly induced in the acute phases (3-6 h) after infection with Gram-positive (M. luteus) and Gram-negative (V. anguillarum) bacteria, and many more dramatic changes in PyTRAFs expression were observed after V. anguillarum challenge. Notably, the strong response in the up-regulation of PyTRAF6 post-bacterial challenge was distinct from that previously reported in scallops and crabs but was similar to that of other shellfish, Echinodermata and even teleost fish. The high level expressions of PyTRAFs in the hemocytes and the gill, and their specific expression patterns after challenges provide insights into the versatile roles and responses of TRAFs in the innate immune system against Gram-negative bacterial pathogens in bivalves.

Molecular Basis for Adaptation of Oysters to Stressful Marine Intertidal Environments

Oysters that occupy estuarine and intertidal habitats have well-developed stress tolerance mechanisms to tolerate harsh and dynamically changing environments. In this review, we summarize common pathways and genomic features in oyster that are responsive to environmental stressors such as temperature, salinity, hypoxia, air exposure, pathogens, and anthropogenic pollutions. We first introduce the key genes involved in several pathways, which constitute the molecular basis for adaptation to stress. We use genome analysis to highlight the strong cellular homeostasis system, a unique adaptive characteristic of oysters. Next, we provide a global view of features of the oyster genome that contribute to stress adaptation, including oyster-specific gene expansion, highly inducible expression, and functional divergence. Finally, we review the consequences of interactions between oysters and the environment from ecological and evolutionary perspectives by discussing mass mortality and adaptive divergence among populations and related species of the genus Crassostrea. We conclude with prospects for future study.

The environmentally tuned transcriptomes of Mytilus mussels

Transcriptomics is a powerful tool for elucidating the molecular mechanisms that underlie the ability of organisms to survive and thrive in dynamic and changing environments. Here, we review the major contributions in this field, and we focus on studies of mussels in the genus Mytilus, which are well-established models for the study of ecological physiology in fluctuating environments. Our review is organized into four main sections. First, we illustrate how the abiotic forces of the intertidal environment drive the rhythmic coupling of gene expression to diel and tidal cycles in Mytilus californianus. Second, we discuss the challenges and pitfalls of conducting transcriptomic studies in field-acclimatized animals. Third, we examine the link between transcriptomic responses to environmental stress and biogeographic distributions in blue mussels, Mytilus trossulus and Mytilus galloprovincialis. Fourth, we present a comparison of transcriptomic datasets and identify 175 genes that share common responses to heat stress across Mytilus species. Taken together, these studies demonstrate that transcriptomics can provide an informative snapshot of the physiological state of an organism within an environmental context. In a comparative framework, transcriptomics can reveal how natural selection has shaped patterns of transcriptional regulation that may ultimately influence biogeography.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

Identification of potential markers and sensitive tissues for low or high salinity stress in an intertidal mud crab (Macrophthalmus japonicus)

Macrophthalmus japonicus is an intertidal mud crab is an ecologically important species in Korea, can tolerate a wide range of natural and anthropogenic stressors. Environmental changes especially salinity cause physiological stress to the marine habitats. Differential gene transcription of M. japonicus tissues provided information about tissue specific responses against salinity. Five potential genes were identified and their transcription levels were determined quantitatively comparison to seawater (SW: 31 ± 1 psu) in M. japonicus gills and hepatopancreas after exposed them to different salinities. Ecdysteroid receptor (Mj-EcR), trypsin (Mj-Tryp), arginine kinase (Mj-AK), lipopolysaccharide and β-1,3-glucan binding protein (Mj-LGBP) and peroxinectin (Mj-Prx) in hepatopancreas up-regulated against different salinities. In contrast, the gills, Mj-EcR, Mj-Tryp and Mj-AK showed late up-regulated responses to 40 psu compared to SW. All genes except Mj-LGBP showed up regulation in the gills as time dependent manner. These genes can be considered as potential markers to assess responses in salinity changes. This study suggests hepatopancreas is a suitable tissue for transcriptional, biochemical and physiological responses analysis on M. japonicus in low and high salinity stress.

Mortality occurrence and pathogen detection in Crassostrea gigas and Mytilus galloprovincialis close-growing in shallow waters (Goro lagoon, Italy)

The complex interactions occurring between farmed bivalves and their potential pathogens in the circumstances of global climate changes are current matter of study, owing to the recurrent production breakdowns reported in Europe and other regions of the world. In the frame of Project FP7-KBBE-2010-4 BIVALIFE, we investigated the occurrence of mortality and potential pathogens during the Spring-Summer transition in Crassostrea gigas and Mytilus galloprovincialis cohabiting in the shallow waters of one northern Italian lagoon (Sacca di Goro, Adriatic Sea) and regarded as susceptible and resistant species, respectively. In 2011, limited bivalve mortality was detected in the open-field trial performed with 6-12 month old spat whereas subsequent trials with 2-3 month old spat produced almost complete (2012) and considerable (2013) oyster mortality. Macroscopical examination and histology excluded the presence of notifiable pathogens but, in the sampling preceding the massive oyster spat mortality of 2012, a μdeleted variant of OsHV-1 DNA was found in wide-ranging amounts in all analyzed oysters in conjunction with substantial levels of Vibrio splendidus and Vibrio aestuarianus. The large oyster spat mortality with borderline OsHV-1 positivity recorded in 2013 supports the multi-factorial etiology of the syndrome. This is the first report of a OsHV-1 (under a form interpreted as the variant μVar) in the Goro lagoon. Transcriptional host footprints are under investigation to better understand the bivalve response to environmental factors, included viral and bacterial pathogens, in relation to the observed mortalities.

Identification of potential general markers of disease resistance in American oysters, Crassostrea virginica through gene expression studies

Several diseases have a significant impact on American oyster populations in the Atlantic coasts of North America. Knowledge about the responses of oysters to pathogenic challenge could help in identifying potential markers of disease resistance and biomarkers of the health status of an oyster population. A previous analysis of the transcriptome of resistant and susceptible American oysters in response to challenge with the bacterial pathogen Roseovarius crassostreae, as well as sequencing of suppression subtractive hybridization libraries from oysters challenged with the protozoan parasite Perkinsus marinus, provided a list of genes potentially involved in disease resistance or susceptibility. We investigated the patterns of inducible gene expression of several of these genes in response to experimental challenge with the oyster pathogens R. crassostreae, Vibrio tubiashii, and P. marinus. Oysters showing differential susceptibility to R. crassostreae demonstrated differential patterns of expression of genes coding for immune (serine protease inhibitor-1, SPI1) and stress-related (heat shock protein 70, HSP70; arginine kinase) proteins 30 days after challenge with this bacterial pathogen. Differential patterns of expression of immune (spi1, galectin and a matrix metalloproteinase) and stress-related (hsp70, histone H4, and arginine kinase) genes was observed in hemocytes from adult oysters challenged with P. marinus, but not with V. tubiashii. While levels of spi1 expression in hemocytes collected 8 and 21 days after P. marinus challenge were negatively correlated with parasite load in oysters tissues at the end of the challenge (62 days), levels of expression of hsp70 in hemocytes collected 1-day after challenge were positively correlated with oyster parasite load at 62 days. Our results confirm previous research on the role of serine protease inhibitor-1 in immunity and disease resistance in oysters. They also suggest that HSP70 and histone H4 could be used as a markers of health status or disease susceptibility in oysters.

Molecular signatures at imminent death: hemocyte gene expression profiling of shrimp succumbing to viral and fungal infections

Infectious diseases represent the most serious threat to shrimp farming worldwide. Understanding the molecular mechanisms driving shrimp-pathogen interactions is necessary for developing strategies to control disease outbreaks in shrimp production systems. In the current study, we experimentally reproduced mortality events using standardized infections to characterize the hemocyte transcriptome response of the shrimp Litopenaeus vannamei succumbing to infectious diseases. By using a high-throughput microfluidic RT-qPCR approach, we identified molecular signatures in shrimp during lethal infections caused by the White Spot Syndrome Virus (WSSV) or the filamentous fungus Fusarium solani. We successfully identified gene expression signatures shared by both infections but also pathogen-specific gene responses. Interestingly, whereas lethal WSSV infection induced the expression of antiviral-related genes, the transcript abundance of many antimicrobial effectors was reduced by lethal F. solani infection. To our knowledge, this is the first report of the immune-gene repertoire of infected shrimp at the brink of death.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

Identification of molecular and physiological responses to chronic environmental challenge in an invasive species: the Pacific oyster, Crassostrea gigas

Understanding the environmental responses of an invasive species is critical in predicting how ecosystem composition may be transformed in the future, especially under climate change. In this study, Crassostrea gigas, a species well adapted to the highly variable intertidal environment, was exposed to the chronic environmental challenges of temperature (19 and 24°C) and pH (ambient seawater and a reduction of 0.4 pH units) in an extended 3-month laboratory-based study. Physiological parameters were measured (condition index, shell growth, respiration, excretion rates, O:N ratios, and ability to repair shell damage) alongside molecular analyses. Temperature was by far the most important stressor, as demonstrated by reduced condition indexes and shell growth at 24°C, with relatively little effect detected for pH. Transcriptional profiling using candidate genes and SOLiD sequencing of mantle tissue revealed that classical “stress” genes, previously reported to be upregulated under acute temperature challenges, were not significantly expressed in any of the treatments, emphasizing the different response between acute and longer term chronic stress. The transcriptional profiling also elaborated on the cellular responses underpinning the physiological results, including the identification of the PI3K/AKT/mTOR pathway as a potentially novel marker for chronic environmental challenge. This study represents a first attempt to understand the energetic consequences of cumulative thermal stress on the intertidal C. gigas which could significantly impact on coastal ecosystem biodiversity and function in the future.

Integrative study of physiological changes associated with bacterial infection in Pacific oyster larvae

Background

Bacterial infections are common in bivalve larvae and can lead to significant mortality, notably in hatcheries. Numerous studies have identified the pathogenic bacteria involved in such mortalities, but physiological changes associated with pathogen exposure at larval stage are still poorly understood. In the present study, we used an integrative approach including physiological, enzymatic, biochemical, and molecular analyses to investigate changes in energy metabolism, lipid remodelling, cellular stress, and immune status of Crassostrea gigas larvae subjected to experimental infection with the pathogenic bacteria Vibrio coralliilyticus.

Findings

Our results showed that V. coralliilyticus exposure induced (1) limited but significant increase of larvae mortality compared with controls, (2) declined feeding activity, which resulted in energy status changes (i.e. reserve consumption, β-oxidation, decline of metabolic rate), (3) fatty acid remodeling of polar lipids (changes in phosphatidylinositol and lysophosphatidylcholine composition`, non-methylene–interrupted fatty acids accumulation, lower content of major C₂₀ polyunsaturated fatty acids as well as activation of desaturases, phospholipase and lipoxygenase), (4) activation of antioxidant defenses (catalase, superoxide dismutase, peroxiredoxin) and cytoprotective processes (heat shock protein 70, pernin), and (5) activation of the immune response (non-self recognition, NF-κκ signaling pathway, haematopoiesis, eiconosoids and lysophosphatidyl acid synthesis, inhibitor of metalloproteinase and antimicrobial peptides).

Conclusion

Overall, our results allowed us to propose an integrative view of changes induced by a bacterial infection in Pacific oyster larvae, opening new perspectives on the response of marine bivalve larvae to infections.

Transcriptomic responses to salinity stress in the Pacific oyster Crassostrea gigas

BACKGROUND

Low salinity is one of the main factors limiting the distribution and survival of marine species. As a euryhaline species, the Pacific oyster Crassostrea gigas is considered to be tolerant to relative low salinity. The genes that regulate C. gigas responses to osmotic stress were monitored using the next-generation sequencing of whole transcriptome with samples taken from gills. By RNAseq technology, transcript catalogs of up- and down-regulated genes were generated from the oysters exposed to low and optimal salinity seawater.

METHODOLOGY/PRINCIPAL FINDINGS

Through Illumina sequencing, we reported 1665 up-regulated transcripts and 1815 down-regulated transcripts. A total of 45771 protein-coding contigs were identified from two groups based on sequence similarities with known proteins. As determined by GO annotation and KEGG pathway mapping, functional annotation of the genes recovered diverse biological functions and processes. The genes that changed expression significantly were highly represented in cellular process and regulation of biological process, intracellular and cell, binding and protein binding according to GO annotation. The results highlighted genes related to osmoregulation, signaling and interactions of osmotic stress response, anti-apoptotic reactions as well as immune response, cell adhesion and communication, cytoskeleton and cell cycle.

CONCLUSIONS/SIGNIFICANCE

Through more than 1.5 million sequence reads and the expression data of the two libraries, the study provided some useful insights into signal transduction pathways in oysters and offered a number of candidate genes as potential markers of tolerance to hypoosmotic stress for oysters. In addition, the characterization of C. gigas transcriptome will not only provide a better understanding of the molecular mechanisms about the response to osmotic stress of the oysters, but also facilitate research into biological processes to find underlying physiological adaptations to hypoosmotic shock for marine invertebrates.

Defining the limits of physiological plasticity: how gene expression can assess and predict the consequences of ocean change

Anthropogenic stressors, such as climate change, are driving fundamental shifts in the abiotic characteristics of marine ecosystems. As the environmental aspects of our world's oceans deviate from evolved norms, of major concern is whether extant marine species possess the capacity to cope with such rapid change. In what many scientists consider the post-genomic era, tools that exploit the availability of DNA sequence information are being increasingly recognized as relevant to questions surrounding ocean change and marine conservation. In this review, we highlight the application of high-throughput gene-expression profiling, primarily transcriptomics, to the field of marine conservation physiology. Through the use of case studies, we illustrate how gene expression can be used to standardize metrics of sub-lethal stress, track organism condition in natural environments and bypass phylogenetic barriers that hinder the application of other physiological techniques to conservation. When coupled with fine-scale monitoring of environmental variables, gene-expression profiling provides a powerful approach to conservation capable of informing diverse issues related to ocean change, from coral bleaching to the spread of invasive species. Integrating novel approaches capable of improving existing conservation strategies, including gene-expression profiling, will be critical to ensuring the ecological and economic health of the global ocean.

A hemocyte gene expression signature correlated with predictive capacity of oysters to survive Vibrio infections

Background

The complex balance between environmental and host factors is an important determinant of susceptibility to infection. Disturbances of this equilibrium may result in multifactorial diseases as illustrated by the summer mortality syndrome, a worldwide and complex phenomenon that affects the oysters, Crassostrea gigas. The summer mortality syndrome reveals a physiological intolerance making this oyster species susceptible to diseases. Exploration of genetic basis governing the oyster resistance or susceptibility to infections is thus a major goal for understanding field mortality events. In this context, we used high-throughput genomic approaches to identify genetic traits that may characterize inherent survival capacities in C. gigas.

Results

Using digital gene expression (DGE), we analyzed the transcriptomes of hemocytes (immunocompetent cells) of oysters able or not able to survive infections by Vibrio species shown to be involved in summer mortalities. Hemocytes were nonlethally collected from oysters before Vibrio experimental infection, and two DGE libraries were generated from individuals that survived or did not survive. Exploration of DGE data and microfluidic qPCR analyses at individual level showed an extraordinary polymorphism in gene expressions, but also a set of hemocyte-expressed genes whose basal mRNA levels discriminate oyster capacity to survive infections by the pathogenic V. splendidus LGP32. Finally, we identified a signature of 14 genes that predicted oyster survival capacity. Their expressions are likely driven by distinct transcriptional regulation processes associated or not associated to gene copy number variation (CNV).

Conclusions

We provide here for the first time in oyster a gene expression survival signature that represents a useful tool for understanding mortality events and for assessing genetic traits of interest for disease resistance selection programs.

Expression of candidate genes related to metabolism, immunity and cellular stress during massive mortality in the American oyster Crassostrea virginica larvae in relation to biochemical and physiological parameters

Quantification of mRNA of genes related to metabolism, immunity and cellular stress was examined in relation to a massive mortality event during the culture of American oyster larvae, Crassostrea virginica which was probably, in regard to previous microbiological analysis, induced by Vibrio infection. To document molecular changes associated with the mortality event, mRNA levels were compared to biochemical and physiological data, previously described in a companion paper. Among the 18 genes studied, comparatively to the antibiotic control, 10 showed a lower relative gene expression when the massive mortality occurred. Six of them are presumed to be related to metabolism, corroborating the metabolic depression associated with the mortality event suggested by biochemical and physiological analyses. Relationships between the regulation of antioxidant enzyme activities, lipid peroxidation, and the mRNA abundance of genes linked to oxidative stress, cytoprotection, and immune response are also discussed. Finally, we observed an increase in the transcript abundance of two genes involved in apoptosis and cell regulation simultaneously with mortality, suggesting that these processes might be linked.